Millimetre Survey of z\sim 2 ULIRGs

A Millimetre Survey of Starburst Dominated Ultraluminous Infrared Galaxies at


We present millimetre observations of a sample of 12 high redshift ultraluminous infrared galaxies (ULIRGs) in the Extended Growth Strip (EGS). These objects were initially selected on the basis of their observed mid–IR colours ( and ) to lie at high redshift , and subsequent 20–38 mid–IR spectroscopy confirms that they lie in a narrow redshift window centered on . We detect 9/12 of the objects in our sample at high significance () with a mean 1200 flux of mJy. Our millimetre photometry, combined with existing far–IR photometry from the Far–IR Deep Extragalactic Legacy (FIDEL) Survey and accurate spectroscopic redshifts, places constraints both sides of the thermal dust peak. This allows us to estimate the dust properties, including the far–IR luminosity, dust temperature, and dust mass. We find that our sample is similar to other high– and intermediate– ULIRGs, and local systems, but has a different dust selection function than submillimeter–selected galaxies. Finally, we use existing 20cm radio continuum imaging to test the far–IR/radio correlation at high redshift. We find that our sample is consistent with the local relation, implying little evolution. Furthermore, this suggests that our sample selection method is efficient at identifying ultraluminous, starburst–dominated systems within a very narrow redshift range centered at .

galaxies: starburst – galaxies: formation – galaxies: high-redshift – submillimetre – radio continuum: galaxies – infrared: galaxies

1 Introduction

First discovered over 30 years ago, infrared (IR) luminous (LIRGs: ) and ultraluminous (ULIRGs: ) galaxies are among the most extreme objects in the universe, with energy outputs rivaling those of bright quasars (see Sanders & Mirabel, 1996, for a review). Locally, theoretical modeling (Mihos & Hernquist, 1994) and both CO spectroscopic (e.g., Solomon et al., 1997) and optical/near–IR imaging (Veilleux et al., 2002) indicate that they are powered in large part by star formation (SF) induced by major mergers. In the context of a merger–driven model of galaxy evolution, they represent the starburst prelude to the rapid, self–regulated growth of a nuclear supermassive black hole (e.g., Silk & Rees, 1998; Murray et al., 2005; Di Matteo et al., 2005; Hopkins et al., 2007; Younger et al., 2008a), a bright quasar, and eventually a “red and dead” early–type massive galaxy (Sanders et al., 1988a, b; Hopkins et al., 2006, 2008b, 2008a).

Though locally they contribute a very small fraction of the infrared luminosity density, at high redshift LIRGs and ULIRGs take on increasing cosmological importance and may dominate cosmic SF at (Blain et al., 1999, 2002; Le Floc’h et al., 2005). There has also been recent observational evidence for a significant populations of hyperluminous sources at still greater redshift (Wang et al., 2007; Younger et al., 2007, 2008b; Dannerbauer et al., 2008; Younger et al., 2008c). Several techniques have been developed for selecting high redshift ULIRGs, including: direct far–IR selection via (sub)millimeter surveys (e.g., Hughes et al., 1998; Barger et al., 1998; Greve et al., 2004; Pope et al., 2006; Scott et al., 2008), dust–obscured galaxies (DOGs) selected on the basis of their 24–R colour (e.g., Houck et al., 2005; Yan et al., 2005), and those selected based on their mid–IR colours (Farrah et al., 2008, Huang et al., in preparation). Understanding the nature of these different populations, and the engine driving their extreme luminosities, is crucial to a thorough understanding of SF and galaxy evolution at high redshift.

Object Name r’
[mag] [Jy] [mJy] [mJy] [mJy] [Jy]
EGS142301.5+523223 egs1
EGS142148.5+521435 egs4
EGS141928.1+524342 egs10
EGS141917.4+524922 egs11
EGS141920.4+525038 egs12
EGS141900.3+524948 egs14
EGS141829.5+523630 egs21
EGS141822.5+523938 egs23
EGS141834.6+524506 egs24
EGS141746.2+523322 egs26
EGS141836.8+524604 egs24a
EGS142219.8+531950 egsb2
Table 1: The observed r’ band fluxes and far–IR SEDs of the ULIRGs in our sample. Upper limits are at the level. The 70,160, and 20 cm photometry is from Huang et al. (in preparation). The 1200 MAMBO observations and data reduction are described in § 3.

At the same time, a remarkably tight correlation between the radio and far–IR luminosity of galaxies, spanning several decades in FIR luminosity, has been known for more than 20 years (Helou et al., 1985; Condon & Broderick, 1991; Condon, 1992; Yun et al., 2001). Though it is thought that this relationship arises from synchrotron loses associated with massive star forming regions (Hummel, 1986; Hummel et al., 1988; Condon, 1992), the detailed physics are relatively poorly understood (see also the discussion in Thompson et al., 2006). Therefore, the persistence of this relation at intermediate to high redshift – where the conditions of galaxy formation and star forming environment were likely quite different than they are locally – is a subject of great interest in current research (Garrett, 2002; Gruppioni et al., 2003; Appleton et al., 2004; Kovács et al., 2006; Boyle et al., 2007; Vlahakis et al., 2007; Ibar et al., 2008; Sajina et al., 2008).

In this work, we present far–IR observations of a sample of 12 mid–IR selected ULIRGs at , first identified by Huang et al. (in preparation). Our far–IR photometry, which spans both sides of the thermal peak, allows us to more robustly estimate the luminosity and dust properties – e.g., dust temperature and mass – of these objects. The relative uniformity of this sample, combined with mid–IR spectroscopic redshifts that indicate a remarkably tight redshift range, makes it ideal for such a study. It furthermore provides a promising platform for probing the far–IR/radio correlation at .

Figure 1: The correlation between the millimetre–to–24 flux density ratio () and the obscuration parameter of Lutz et al. (2005): . Upper limits are at the level are indicated by inverted, filled triangles. We find the all of our sources fall above the mean measured by Lutz et al. (2005, shaded region and dashed line), with a mean of (dotted line) with an intrinsic dispersion of . Therefore, our sources are more far–IR luminous at fixed 24 independent of obscuration, which is consistent with a larger contribution from a starburst starburst over a dusty AGN.

This paper is organized as follows: in § 2 we outline our sample selection, in § 3 we detail our observations, in § 4 we present the results of our observations, in § 5 we describe our far–IR spectral energy distribution fitting method, in § 6 we present estimates of the dust properties for objects in our sample, in § 7 we examine the far–IR/radio correlation at high redshift, and in § 8 we conclude. Throughout this work, we assume a flat concordance cosmology with , and (Spergel et al., 2003, 2007; Tegmark et al., 2004, 2006).

2 Sample Selection

It has been shown (Sawicki, 2002; Huang et al., 2004; Papovich, 2008) that an efficient method for selecting high redshift galaxies is to track the 1.6 stellar emission bump via mid–IR colors. Using the Infrared Array Camera (IRAC: Fazio et al., 2004) on board the Spitzer Space Telescope (Werner et al., 2004), applying the color selection criteria


where [3.6] denotes the AB magnitude in the 3.6 band (and likewise for the 4.5 and 8.0 bands), will identify massive galaxies in the redshift range – the redshift range where the peak of the 1.6 stellar bump is redshifted into the 4.5 and 5.8 bands. Furthermore, since this is effectively a rest frame near–IR selection criterion, it is relatively insensitive to dust reddening and is roughly stellar–mass limited (Huang et al., 2004, 2005; Conselice et al., 2007a).

Huang et al. (in preparation) used this technique, combined with a 24 flux limit of mJy using the Multiband Imaging Photometer for Spitzer (MIPS: Rieke et al., 2004), to select high redshift, massive, IR luminous systems for mid–IR spectroscopic followup using the Infrared Spectrograph (IRS: Houck et al., 2004) – both also on board Spitzer. This yielded a sample of 12 objects in the Extended Groth Strip (EGS), all with in a remarkably tight redshift range centered around . This was also a considerably more uniform sample of high redshift IR luminous galaxies than other Spitzer selected samples; see Huang et al. (in preparartion) for a more detailed comparison to the various other selection techniques. Our sample was drawn from the Huang et al. (in preparation) EGS sample with mid–IR derived spectroscopy redshifts.

3 Observations

The optical and near– to mid–IR observations used for this work are described in detail in Davis et al. (2007), Huang et al. (2008, in preparation), and references therein. Briefly, we make extensive use of deep imaging with a limiting magnitude of 26.5 AB, IRAC imaging at (3.6, 4.5, 5.8, 8.0) to flux limits of (1,1, 1.2, 6.3, 6.9) (see also Barmby et al., 2004, 2008), and MIPS imaging at 24 to a depth of 77 (Sanders et al., 2007). We also make use of extremely deep far–IR MIPS imaging at 70 and 160 that combined existing observations (PID 00008, PI Fazio) and the Far–Infrared Deep Extragalactic Survey (FIDELS; PID 30948, PI Dickinson)2, and deep 1.4 GHz radio imaging to a depth of 50 (Ivison et al., 2007). Photometry for our sources was first presented by Huang et al. (2008, in preparation).

Figure 2: The millimetre detectability as a function of the 70– (: left) and 160–to–millimetre (: right) flux density ratio. Upper and lower limits are indicated by arrows. We find a strong negative correlation at the level (see Table 2) between the millimetre flux and , reflecting the fact that warmer sources are less likely to have a significant cold dust component. We do not, however, find a significant correlation between the millimetre flux and .

The millimetre observations presented here were performed at the Institut de Radioastronomie Millimetrique (IRAM) 30 metre telescope at Pico Veleta during the Winter/Spring of 2008, using the 117 element Max–Plank Bolometer Array 2 (MAMBO2: Kreysa et al., 1998). All our observations were performed in good weather conditions, with and 1.9 Hz (corresponding to a 0.52s wobbler “period”) chopped sky noise mJy/beam (HPBW=11 arcsec). Primary pointing, focus, and flux calibration was performed on either Mars ( Jy) or another bright calibrator source – including K3–50A ( Jy) and CW–LEO ( Jy). Secondary pointing was performed on J1419+544 ( Jy), a bright quasar only 1 degree away. On–off observations were obtained in 20 minute scans with a wobbler throw of 35 arcseconds, which were repeated until either a significant detection (S/N ) or a noise level of mJy. This typically necessitated hours of integration time per source (full range 1–3 hours). The data were reduced using the MOPSIC Package developed by R. Zylka, using the standard reduction pipeline.

4 Results

In all, 9 of the 12 targets had significant detections () at 1200– egs11 was detected at , egs24 was detected at , and egsb2 was consistent with the sky noise. Excluding egsb2, the mean observed flux density mJy and mid– to far–IR colour with an intrinsic dispersion of . The average 1200 flux density and detection rate is similar to the IRAC selected starburst sample Lonsdale et al. (2008). This average flux density ratio noticeably lower than is typical of bright SMGs (: Egami et al., 2004; Ivison et al., 2004). However, while it is consistent with the detected sources in their sample, this is a significantly higher average value than has been observed in other Spitzer selected samples of high redshift ULIRGs: (Lutz et al., 2005) followed up a sample of 40 obscured 24 sources (see also Houck et al., 2005; Yan et al., 2005; Weedman et al., 2006a, b) with MAMBO at 1.2mm and found an overall with , which is consistent with the expectation from an IR bright active galactic nucleus (AGN: Elvis et al., 1994). Sajina et al. (2008) found similar millimetre fluxes for a comparable sample (Yan et al., 2007). This is not unexpected, as obscured 24 sources are thought to contain a large fraction of AGN (Weedman et al., 2006b).

In Figure 1, we show the obscuration parameter from Lutz et al. (2005), , where is taken from the imaging (Davis et al., 2007). These authors found that when they restricted their sample to the most obscured sources with , the average increased to with . Though they could not examine it directly due to the low detection rate of sources in their sample, this suggests a correlation between and . For our sources with , a Spearman rank correlation coefficient analysis finds a mild correlation with at the level. The low significance of this correlation may be due to small sample size, but it is broadly consistent with the findings of Lutz et al. (2005) and suggests that the most heavily obscured sources efficiently reprocess the attenuated stellar radiation into cold ( K) thermal dust emission.

Given the narrow redshift range of our sources, we can also use the 70 and 160–to–millimetre flux density ratios to investigate the correlation between strong millimetre emission and the dust properties of objects in our sample. At , the 70 and 160 MIPS channels probe the rest–frame 24 and 50 emission. These wavelengths probably probe AGN–heated dust (e.g., Elvis et al., 1994; Urry & Padovani, 1995) and cold dust from a starburst (Sanders & Mirabel, 1996) respectively. In Figure 2 we show the millimetre flux as a function of both and respectively, and in Table 2 we present the results of a Pearson rank correlation coefficient analysis. We find a strong negative correlation between the millimetre flux and ; sources with more hot relative to cold dust – and therefore likely a larger AGN contribution – are less likely to be detected at 1200. We do not however, find a strong correlation between and . We believe that this is mostly likely a result of the relatively poor sensitivity of the 160 FIDELS survey: if an object has enough cold dust to be detected at 160 at all, it is also likely to be a strong millimetre source.

Parameter Significance
-0.76 2.0
-0.10 0.3
-0.04 0.1
0.26 0.8
-0.62 1.6
0.31 0.7
Table 2: The results of a Spearman rank correlation coefficient () analysis of the trends of different PAH parameters (see Figure 3 with the observed millimeter flux density .

It been observed that the PAH luminosity is correlated with the far–IR (Rigopoulou et al., 1999; Wu et al., 2005; Brandl et al., 2006; Calzetti et al., 2007; Desai et al., 2007; Shi et al., 2007), and therefore we would expect strong PAH sources to show strong millimetre emission. To investigate this for our sample, which contains all strong PAH emitters, in Figure 3 we show the correlations between the millimetre flux () and PAH luminosities ( and ) and equivalent widths ( and ) at 7.7 and 11.3. The results of a Spearman rank correlation coefficient analysis are presented in Table 2. All of these results are low significance, likely due to the effects of a relatively small and homogenous sample. They are broadly consistent – considering the low significance of the inverse correlation – with observations of local systems that found a correlation between the PAH luminosity and ongoing star formation (e.g., Wu et al., 2005; Brandl et al., 2006; Calzetti et al., 2007). In particular, we find a positive correlation between 11.3 PAH luminosity and equivalent width and the millimetre flux at a and level of significance respectively. Since the 11.3 feature is more sensitive to obscuration (Rigopoulou et al., 1999), sources with strong observed features at this wavelength are likely to be more intense starbursts, and thus more likely to be detected in the far–IR. However, we find a significant negative correlation between the 7.7 equivalent width and at significance. This differs from the results of Sajina et al. (2008) who find a positive correlation () between and at confidence for a sample of high redshift obscured Spitzer selected ULIRGs (Yan et al., 2007). While the correlation seen for our sources is, again, low–significance, it suggests that the relationship between strong PAH emission and far–IR luminosity is somewhat different in their sample – a reasonable result considering the far different mid– to far–IR properties of similar objects (Lutz et al., 2005) are significantly different from those presented here.

Figure 3: The correlation between the observed millimeter flux density and PAH luminosities (left: and ) and equivalent widths (right: and ) at 7.7 (top) and 11.3 (bottom). The results of a Pearson rank correlation coefficient analysis are shown in Table 2.

5 Spectral Energy Distribution Fitting

Figure 4: The observed far–IR SEDs (see Table 1) of the ULIRGs in our sample, along with the fitted model. The model parameters and derived quantities are listed in each panel, as well as tabulated in Table 3, and the fitting procedure summarized in § 5. The solid lines have a fixed dust emissivity , while the dotted lines are the same for .
Figure 5: Same as Figure 4.

Our millimetre observations constrain the shape of the far–IR SED of our sources with photometry on both sides of the thermal dust peak. Given a choice of a simple far–IR SED parameterization, this allows us to characterize the dust properties in a sample of very luminous, high redshift systems from the far–IR directly, without invoking additional assumptions – for example assuming the local far–IR/radio correlation (e.g., Chapman et al., 2005). For this work, we use a single temperature greybody fitting form first introduced by Hildebrand (1983), in which the thermal dust spectrum is approximated as , where is the Planck function, is the dust emissivity with , and is the effective – or emission weighted – dust temperature. This yields a dust spectrum that for takes the form:


This is then typically paired with a power–law component at short wavelengths , to account for a subdominant warmer dust component arising either from the starburst or an AGN, which is matched to the long wavelength SED at to ensure a continuous first derivative. Blain et al. (2003) find that this form provides a good match to the SEDs of observed systems from over a large range in total IR luminosities (see also Blain, 1999; Carilli & Yun, 1999; Dunne et al., 2000; Yun & Carilli, 2002; Sajina et al., 2006; Yang et al., 2007, for additional discussion of far–IR fitting methods)

This model was fit to the 70, 160, and 1200 data using a Levenberg-Marquardt least-squares minimization routine, the results of which are summarized in Figures 4 and 5 and Table 3. Since is degenerate with (Blain et al., 2003; Sajina et al., 2006) for sparsely sampled far–IR SEDs, we fix which is consistent with the results of SED fitting of local and high redshift systems (Dunne et al., 2000; Yun & Carilli, 2002; Yang & Phillips, 2007); using a higher value of will lead to systematically lower dust temperatures (Blain et al., 2003; Sajina et al., 2006). Furthermore, the accurate redshifts supplied by mid–IR spectroscopy were essential in removing the additional degeneracy between and redshift (Carilli & Yun, 1999; Yun & Carilli, 2002; Blain et al., 2003). From these fitted models we can estimate and the far–IR luminosity. We can also estimate the total dust mass in these systems according to:


where is the total dust mass, is the observed flux density, is the luminosity distance, is the rest frame dust mass absorption coefficient at the observed wavelength (for which we use the Milky Way dust model of Weingartner & Draine, 2001), and is the Planck function at the rest wavelength. At long wavelengths – e.g., 1200 observed, in the rest frame – observations are directly probing populations of dust that are preferentially colder than , which is effectively an emission weighted average over dust populations at a large range of physical temperatures. Therefore, we use the fitted models to estimate the peak flux and wavelength, which is dominated by dust at and simultaneously minimizes the effects of uncertainty in (), for input into equation 4.

Finally, we compute the parameter for the far–IR/radio correlation as defined by Helou et al. (1985):


where is the integrated IR luminosity between 40 and 120, and is the rest–frame 1.4 GHz continuum radio luminosity. Since our sources are at , their observed 1.4 GHz flux corresponds to a rest–frame frequency of 4.2 GHz. Therefore, we applied a –correction similar to that used in Yang et al. (2007), assuming a power–law synchrotron , where we assume as in typical local star–forming systems (Condon, 1983, 1992). The results, as with the far–IR luminosity and dust properties, are listed in Table 3, and discussed in the following sections.

[K] [] [] []
egs1 1.88
egs4 1.90
egs10 1.94
egs11 1.80
egs12 2.01
egs14 1.86
egs21 3.01
egs23 1.78
egs24 1.85
egs26 1.80
egs24a 2.12
egsb2 1.59
Table 3: Derived parameters for our high redshift ULIRG sample, including: the mid–IR spectroscopic redshift (from Huang et al. in preparation) and the effective dust temperature (), far–IR luminosity from 40–1000 () and 40–120 (), dust mass (), and far–IR/radio parameter (see § 7 and Condon, 1992). For details on the fitting procedure, see § 5 and Blain et al. (2003).

6 Dust Properties

Our far–IR and millimetre photometry that constrains both sides of the cold dust peak allows a determination of the far–IR luminosity, dust temperature, and dust mass from the far–IR alone. As a result, the dust properties in our high redshift ULIRG sample can be determined with a minimum of additional assumptions. While the properties of dust in these extreme environments and at high redshift is of itself an interesting result, it is particularly so in the context of analogous and related objects at different redshifts and with different selection techniques. Therefore, we assembled results from five samples with extensive far–IR observations from the literature for comparison: another Spitzer–selected high redshift sample (HZS), a low redshift sample (LZS), intermediate redshift sample (IZS), a submillimetre selected sample (SMS), and finally a sample of quasars – which are thought to follow the ULIRG phase in an evolutionary sequence (Sanders et al., 1988a, b, c; Hopkins et al., 2006, 2008b, 2008a) – with extensive far–IR observations (QS). Our far–IR SED parameterization is consistent with those used for the analysis of all samples discussed below – including the HZS, IZS, SMS, and QS – and therefore there are no hidden systematic biases.

The HZS consists of bright, obscured 24 sources sected on the basis of their ratio (Houck et al., 2005; Yan et al., 2005; Weedman et al., 2006a, b), which have mid–IR spectrscopic (Yan et al., 2007; Sajina et al., 2007) and MAMBO 1200 photometric followup (Lutz et al., 2005; Sajina et al., 2008). As we note in § 4, Lutz et al. (2005) find a signficantly lower average ratio than the objects in our sample, even among the most obscured () subsample. The authors find a median dust temperature of K ( K for 14 strong–PAH srouces, and K for 34 weak–PAH sources), as compared to K for our sample of 12 objects, with comparable far–IR luminosities and dust masses. This difference may arise from the nature of their sources, combined with their SED fitting technique; the HZS has shown a significant contribution to the far–IR by an AGN (see also § 7 Weedman et al., 2006b; Sajina et al., 2008), which the authors account for by fitting a combination of models to the mid–IR (see Sajina et al., 2007) and far–IR separately in order to separate out the cold dust population. However, it is most likely driven largely by their choice of , which leads to cooler fitted dust temperatures than the adopted here (Blain et al., 2003; Sajina et al., 2006) – indeed, when we adopt we find a median dust temperature of K while the total IR luminosity decreases by . Furthermore, since their mid–IR models will implicitly include some of the warm dust heated by the starburst, their two–component fitting method may also lead to somewhat lower dust temperatures. Therefore, while theirs is a rather different population of objects, we find that the dust properties in both are similar. This is an interesting result, considering the large fraction of obscured AGN in such samples (Weedman et al., 2006b; Sajina et al., 2008); the cold dust appears to be similar in both AGN and starburst dominated systems.

The LZS is the SCUBA Local Universe Galaxy Survey (SLUGS Dunne et al., 2000), which consists of 104 local () infrared luminous galaxies selected from the 60 flux limited IRAS Bright Source Catalog (BSC: Soifer et al., 1989) with complete 450 and 850 followup with the Submillimetre Common User Bolometer Array (SCUBA: Holland et al., 1999) at the James Clerk Maxwell Telescope (JCMT). The IZS (Yang et al., 2007) was drawn from the FSC–FIRST Survey (Stanford et al., 2000), which consists of a matched sample of 60 IRAS Faint Source Catalog (FSC: Moshir et al., 1992) and 1.4 GHz sources from the Faint Images of the Radio Sky at Twenty Centimeters (FIRST: Becker et al., 1995) Survey, with complete 350 followup from the SHARC–II camera (Dowell et al., 2003) at the Caltech Submillimetre Observatory (CSO). The SMS (Kovács et al., 2006) consists of SMGs identified at 850 in blank field SCUBA surveys with 20cm radio counterparts, optical redshifts (Chapman et al., 2005), and complete 350 followup with SHARCH-II (for a similar sample, see Coppin et al., 2008). Finally, the QS (Haas et al., 2000, 2003) consists of 64 quasars selected from the Palomar Green Survey (PG: Schmidt & Green, 1983) with far–IR coverage from 5–200 from ISOPHOT (Lemke et al., 1996) on board the ISO satellite and some (sub)millimetre coverage with MAMBO and SCUBA.

Figure 6: The temperature–luminosity relation , where is taken from 40-1000 . The fitting method is detailed in § 5, and the values for our sample tabulated in Table 3. Included are results from our sample (red circles) as compared to local (D00: Dunne et al., 2000, black triangles), intermediate redshift (Y07: Yang et al., 2007, open stars), high redshift submillimeter galaxies (K06: Kovács et al., 2006, black squares), and PG quasars (H03 Haas et al., 2003, blue circles). The dashed lines indicated tracks in constant dust mass (see Equation 7 and § 6) with (left to right) , , , and , assuming a fixed emissivity of . The dotted lines are the same for .
Figure 7: Same Figure 6 for the luminosity–dust mass relation. The dashed lines are tracks in constant dust temperature with (left to right) , 30, 40, 50, and 60 K.
Figure 8: Same Figure 6 for the luminosity–dust mass relation. The dashed lines are tracks in constant far–IR luminosity with (left to right) , , and .

For single temperature dust models, the total far–IR luminosity can be related to the dust mass and temperature via the dust mass absorption coefficient (De Breuck et al., 2003; Yang et al., 2007):


where is the thermal dust spectrum at rest–frame wavelength . For greybody thermal emission, this yields:


where is the dust emissivity, is the gamma function, and is the Riemann Zeta function. As a toy model, we use the Milky Way dust model of Weingartner & Draine (2001) at a reference wavelength of 125 ( g cm), and make the assumption that and does not vary significantly between objects – which is approximately the case locally (e.g., Dunne et al., 2000; Klaas et al., 2001; Yun & Carilli, 2002; Yang & Phillips, 2007). Also, for all the following comparisons the total far–IR luminosity is taken from 40–1000 () unless otherwise stated.

In Figure 6 we show the relation for our sources as compared to the LZS, IZS, SMS, and QS. The dashed lines represent tracks (left to right) in constant , , and . We find the objects in our sample all fall within a relatively narrow range of , though they vary by an order of magnitude in . This is similar to both the LZS and IZS, but characteristically different from the SMS, which resembles a constant selection with a broad distribution extending to lower temperatures. The QS, by contrast, does not follow a clear trend in , thought those with high luminosities tend to have higher dust temperatures and lower dust masses than U/HyLIRGs of comparable luminosity.

This is further illustrated in Figures 7 and 8, where we present the and scalings respectively. In Figure 7, the dashed lines represent tracks in constant (left to right) , 40, 50, and 60 K, while in Figure 8 they represent tracks (left to right) in constant , , and . These two projections also show the dust selection function of the different samples samples, whereby the LZS (roughtly), IZS and our sample follow a track in constant temperature, the SMS is constant in dust mass, and the QS is preferentially warmer and high luminosity with lower dust mass.

7 The Far–IR/Radio Correlation at High Redshift

It has long been recognized that there exists a tight correlation between the far–IR and synchrotron radio emission in star–forming galaxies, spanning nearly four decades in luminosity (Helou et al., 1985; Condon & Broderick, 1991; Condon, 1992; Yun et al., 2001). Though the origin of this correlation is an active subject of ongoing research, the popular wisdom states that it arises from a “cosmic conspiracy” that fixes the ratio of synchrotron losses due to magnetic fields associated with massive star formation and inverse–Compton from the ambient radiation field (Hummel, 1986; Hummel et al., 1988; Condon, 1992). This hypothesis is apparently born out in observations of local starbursts, which show a roughly constant far–IR/radio correlation over four decades in magnetic field energy density (Condon et al., 1991).

However, should the far–IR/radio correlation arise from this fixed ratio of synchrotron and inverse–Compton losses, one might expect it to evolve with redshift. The physics that govern this ratio are complex and not entirely understood (see e.g., Voelk, 1989; Thompson et al., 2006). Therefore, it is important to check if it remains valid and/or constant at high redshift and in more luminous systems.

At intermediate redshift , there is some evidence for the constacy of the far–IR/radio correlation (Garrett, 2002; Gruppioni et al., 2003; Appleton et al., 2004; Boyle et al., 2007). However, at high redshift , the picture is somewhat less clear (Kovács et al., 2006; Vlahakis et al., 2007; Ibar et al., 2008; Sajina et al., 2008), and the effects of sample selection – in particular, contamination by radio AGN – are uncertain. Our sample offers an attractive opportunity to probe the far–IR/radio correlation in a relatively narrow redshift range, and for a well–defined sample.

Figure 9: The evolution of the far–IR/radio correlation with redshift parameterized by (see Equation 5 and § 5 for definitions). The labeling is the same as Figure 6, and the grey shaded region indicates the local value and associated dispersion from Yun et al. (2001).
Figure 10: The correlation between the far–IR/radio correlation coefficient (see Equation 5 and § 5 for definitions) with the far–IR luminosity . The labeling is the same as Figure 6, and the grey shaded region indicates the local value and associated dispersion from Yun et al. (2001).

Overall, our sample has a mean (see Equation 5) with intrinsic dispersion – well in line with results for local systems. In fact, only one object (egsb2) is inconsistent with the local , and it is known to harbor an X–ray luminous AGN (Huang et al. in prepration). This supports the conclusion of Huang et al. (2008, in preparation) that AGNs contribute little to the bolometric luminosity of these objects. Furthermore, it is substantially different from the resutls of Sajina et al. (2008) for obscured 24 sources, which cover a wide range in both for strong and weak PAH emitters – the mean with for their “strong–PAH” sources with radio detections, and with for their “weak–PAH” sources. This is, however, not surprising considering this sample selection is known to be biased towards obscured AGN (Weedman et al., 2006b). We take this as evidence to suggest that our selection criteria much more efficiently picks starburst–dominated systems in a narrow redshift range.

The radio/far–IR properties of this sample are also consistent with the AGN contribution inferred from their X–ray/far–IR properties. The rest–frame 25–60 flux density ratio versus – where is the rest–frame full band 2–10 keV X–ray luminosity – has been shown to be sensitive to AGN activity (de Grijp et al., 1985, 1987; Risaliti et al., 2000). Huang et al. (in preparation) find that only three objects in this sample have even marginal X–ray detections: egs1 ( ergs s, , ), egs14 ( ergs s, , ), and egsb2 ( ergs s, , ). Of these, only egsb2 is consistent with a significant contribution from an obscured AGN (see Figure 5 of Risaliti et al., 2000). Also, while egs24 is not detected in the X–rays, its 25–60 flux density ratio suggests the presence of a Compton–thick AGN ( cm; Risaliti et al., 2000). The remaining 8 objects in our sample have a stacked X–ray flux corresponding to ergs s, along with () and , which are all consistent with a pure starburst.

In Figures 9 and 10 we present the scaling with redshift and far–IR luminosity for our sample (red circles), local systems (grey shaded region with ; Yun et al., 2001), the IZS (open stars), and the SMS (black squares). Immediately it appears that the local far–IR/radio correlation applies over a wide redshift range over more than two decades in luminosity. The apparent inverse trend of with both and redshift in the IZS is probably due to their sample selection, which is more likely to pick out IR bright AGN at higher redshifts (and luminosities). By contrast, both our sample and the SMS show a clear trend towards higher at higher : the Spearman rank correlation coefficients are at and at for our sample and the SMS respectively. This trend is not seen in surveys of local systems spanning (Yun et al., 2001). However, it is not clear whether or not this effect is due to sample selection – our sample and the SMS are selected in very different ways, and show different trends in dust properties (see § 6) – or the strong redshift evolution of the bright end of the galaxy IR luminosity function (Sanders & Mirabel, 1996; Le Floc’h et al., 2005; Hopkins et al., 2008b) and that of massive galaxies (e.g., Conselice et al., 2007b) versus the comparatively weaker evolution of flat–spectrum radio sources (Dunlop & Peacock, 1990). Such a selection effect would potentially bias systems with luminosities of to radio–loud AGN as star formation moves preferentially to higher luminosity objects.

8 Conclusion

We present far–IR observations, including 1200 photometry with MAMBO, of high redshift ULIRGs selected by Huang et al. (in preparation) on the basis of their observed mid–IR colour. All the objects in the sample have mid–IR spectroscopic redshift and lie in a remarkably tight redshift window centered at . We successfully detected (S/N ) 9 of 12 total sources. When combined with existing MIPS photometry at 70 and 160, this provides constraints on both sides of the cold dust peak, allowing a more robust determination of the far–IR SED.

We fit a model to the photometry, from which we estimate both the far–IR luminosity and dust properties – including temperature and mass. When compare to other samples, including IR luminous systems at low and intermediate redshift, submillimetre selected galaxies, and QSOs, and find that our sample, along with the low and intermediate redshift sources, is roughly a constant dust temperature selection. Submillimetre sources, by contrast, are more similar to a constant dust mass selection. Finally, we use existing radio observations at 20cm to test the far–IR/radio correlation at high redshift. We find that all but one of the sources in our sample is consistent with the local relation, suggesting that it remains valid at high redshift. This furthermore indicates that the Huang et al. (in preparation) method efficiently selects massive, high redshift, starburst–dominated systems in a remarkably tight redshift range.


Thanks to Lars Hernquist, Philip Hopkins, T. J. Cox, Dusan Keres, Chris Hayward, and Stephanie Bush for helpful discussions, and to the referee, Stephen Serjeant, for his helpful suggestions for improving this manuscript. This work is based in part on observations made with the Spitzer Space Telescope, which is operated by the Jet Propulsion Laboratory, California Institute of Technology under a contract with NASA. Support for this work was provided by NASA through an award issued by JPL/Caltech. This work includes observations made with IRAM, which is supported by INSU/CNRS (France), MPG (Germany) and IGN (Spain). Thanks to the staff of IRAM Granada for their support in taking the MAMBO observations presented here.


  1. pagerange: 000–000
  2. For a description of the FIDELS observing strategy and data reduction, see


  1. Appleton, P. N. et al. 2004, ApJS, 154, 147
  2. Barger, A. J. et al. 1998, Nature, 394, 248
  3. Barmby, P., Huang, J. ., Ashby, M. L. N., Eisenhardt, P. R. M., Fazio, G. G., & Wright, E. L. 2008, ApJ, in press [astro-ph/0803.0748], 803
  4. Barmby, P. et al. 2004, ApJS, 154, 97
  5. Becker, R. H., White, R. L., & Helfand, D. J. 1995, ApJ, 450, 559
  6. Blain, A. W. 1999, MNRAS, 309, 955
  7. Blain, A. W., Barnard, V. E., & Chapman, S. C. 2003, MNRAS, 338, 733
  8. Blain, A. W., Smail, I., Ivison, R. J., & Kneib, J.-P. 1999, MNRAS, 302, 632
  9. Blain, A. W., Smail, I., Ivison, R. J., Kneib, J.-P., & Frayer, D. T. 2002, Phys. Rep., 369, 111
  10. Boyle, B. J. et al. 2007, MNRAS, 376, 1182
  11. Brandl, B. R. et al. 2006, ApJ, 653, 1129
  12. Calzetti, D. et al. 2007, ArXiv e-prints, 705
  13. Carilli, C. L. & Yun, M. S. 1999, ApJ, 513, L13
  14. Chapman, S. C., Blain, A. W., Smail, I., & Ivison, R. J. 2005, ApJ, 622, 772
  15. Condon, J. J. 1983, ApJS, 53, 459
  16. —. 1992, ARA&A, 30, 575
  17. Condon, J. J. & Broderick, J. J. 1991, AJ, 102, 1663
  18. Condon, J. J., Huang, Z.-P., Yin, Q. F., & Thuan, T. X. 1991, ApJ, 378, 65
  19. Conselice, C. J. et al. 2007a, ApJ, 660, L55
  20. —. 2007b, MNRAS, 381, 962
  21. Coppin, K. et al. 2008, MNRAS, 384, 1597
  22. Dannerbauer, H., Walter, F., & Morrison, G. 2008, ApJ, 673, L127
  23. Davis, M. et al. 2007, ApJ, 660, L1
  24. De Breuck, C. et al. 2003, A&A, 401, 911
  25. de Grijp, M. H. K., Lub, J., & Miley, G. K. 1987, A&AS, 70, 95
  26. de Grijp, M. H. K., Miley, G. K., Lub, J., & de Jong, T. 1985, Nature, 314, 240
  27. Desai, V. et al. 2007, ApJ, 669, 810
  28. Di Matteo, T., Springel, V., & Hernquist, L. 2005, Nature, 433, 604
  29. Dowell, C. D. et al. 2003, in Presented at the Society of Photo-Optical Instrumentation Engineers (SPIE) Conference, Vol. 4855, Millimeter and Submillimeter Detectors for Astronomy. Edited by Phillips, Thomas G.; Zmuidzinas, Jonas. Proceedings of the SPIE, Volume 4855, pp. 73-87 (2003)., ed. T. G. Phillips & J. Zmuidzinas, 73–87
  30. Dunlop, J. S. & Peacock, J. A. 1990, MNRAS, 247, 19
  31. Dunne, L. et al. 2000, MNRAS, 315, 115
  32. Egami, E. et al. 2004, ApJS, 154, 130
  33. Elvis, M. et al. 1994, ApJS, 95, 1
  34. Farrah, D. et al. 2008, ApJ, 677, 957
  35. Fazio, G. G. et al. 2004, ApJS, 154, 10
  36. Garrett, M. A. 2002, A&A, 384, L19
  37. Greve, T. R. et al. 2004, MNRAS, 354, 779
  38. Gruppioni, C. et al. 2003, MNRAS, 341, L1
  39. Haas, M. et al. 2000, A&A, 354, 453
  40. —. 2003, A&A, 402, 87
  41. Helou, G., Soifer, B. T., & Rowan-Robinson, M. 1985, ApJ, 298, L7
  42. Hildebrand, R. H. 1983, QJRAS, 24, 267
  43. Holland, W. S. et al. 1999, MNRAS, 303, 659
  44. Hopkins, P. F., Cox, T. J., Kereš, D., & Hernquist, L. 2008a, ApJS, 175, 390
  45. Hopkins, P. F., Hernquist, L., Cox, T. J., & Kereš, D. 2008b, ApJS, 175, 356
  46. Hopkins, P. F., Hernquist, L., Cox, T. J., Robertson, B., & Krause, E. 2007, ApJ, 669, 45
  47. Hopkins, P. F. et al. 2006, ApJS, 163, 1
  48. Houck, J. R. et al. 2004, ApJS, 154, 18
  49. —. 2005, ApJ, 622, L105
  50. Huang, J.-S. et al. 2004, ApJS, 154, 44
  51. —. 2005, ApJ, 634, 137
  52. Hughes, D. H. et al. 1998, Nature, 394, 241
  53. Hummel, E. 1986, A&A, 160, L4
  54. Hummel, E., Davies, R. D., Pedlar, A., Wolstencroft, R. D., & van der Hulst, J. M. 1988, A&A, 199, 91
  55. Ibar, E. et al. 2008, MNRAS, 386, 953
  56. Ivison, R. J. et al. 2004, ApJS, 154, 124
  57. —. 2007, ApJ, 660, L77
  58. Klaas, U. et al. 2001, A&A, 379, 823
  59. Kovács, A. et al. 2006, ApJ, 650, 592
  60. Kreysa, E. et al. 1998, in Proc. SPIE Vol. 3357, p. 319-325, Advanced Technology MMW, Radio, and Terahertz Telescopes, Thomas G. Phillips; Ed., ed. T. G. Phillips, 319–325
  61. Le Floc’h, E. et al. 2005, ApJ, 632, 169
  62. Lemke, D. et al. 1996, A&A, 315, L64
  63. Lonsdale, C. J. et al. 2008, ApJ, in press [astro-ph/0810.1293]
  64. Lutz, D. et al. 2005, ApJ, 625, L83
  65. Mihos, J. C. & Hernquist, L. 1994, ApJ, 431, L9
  66. Moshir, M., Kopman, G., & Conrow, T. A. O. 1992, IRAS Faint Source Survey, Explanatory supplement version 2 (Pasadena: Infrared Processing and Analysis Center, California Institute of Technology, 1992, edited by Moshir, M.; Kopman, G.; Conrow, T. a.o.)
  67. Murray, N., Quataert, E., & Thompson, T. A. 2005, ApJ, 618, 569
  68. Papovich, C. 2008, ApJ, 676, 206
  69. Pope, A. et al. 2006, MNRAS, 370, 1185
  70. Rieke, G. H., , et al. 2004, ApJS, 154, 25
  71. Rigopoulou, D., Spoon, H. W. W., Genzel, R., Lutz, D., Moorwood, A. F. M., & Tran, Q. D. 1999, AJ, 118, 2625
  72. Risaliti, G., Gilli, R., Maiolino, R., & Salvati, M. 2000, A&A, 357, 13
  73. Sajina, A., Scott, D., Dennefeld, M., Dole, H., Lacy, M., & Lagache, G. 2006, MNRAS, 369, 939
  74. Sajina, A., Yan, L., Armus, L., Choi, P., Fadda, D., Helou, G., & Spoon, H. 2007, ApJ, 664, 713
  75. Sajina, A. et al. 2008, ApJ, 683, 659
  76. Sanders, D. B. & Mirabel, I. F. 1996, ARA&A, 34, 749
  77. Sanders, D. B., Soifer, B. T., Elias, J. H., Madore, B. F., Matthews, K., Neugebauer, G., & Scoville, N. Z. 1988a, ApJ, 325, 74
  78. Sanders, D. B., Soifer, B. T., Elias, J. H., Neugebauer, G., & Matthews, K. 1988b, ApJ, 328, L35
  79. Sanders, D. B., Soifer, B. T., & Scoville, N. Z. 1988c, Science, 239, 625
  80. Sanders, D. B. et al. 2007, ApJS, 172, 86
  81. Sawicki, M. 2002, AJ, 124, 3050
  82. Schmidt, M. & Green, R. F. 1983, ApJ, 269, 352
  83. Scott, K. S. et al. 2008, MNRAS, 385, 2225
  84. Shi, Y. et al. 2007, ApJ, 669, 841
  85. Silk, J. & Rees, M. J. 1998, A&A, 331, L1
  86. Soifer, B. T., Boehmer, L., Neugebauer, G., & Sanders, D. B. 1989, AJ, 98, 766
  87. Solomon, P. M., Downes, D., Radford, S. J. E., & Barrett, J. W. 1997, ApJ, 478, 144
  88. Spergel, D. N. et al. 2003, ApJS, 148, 175
  89. —. 2007, ApJS, 170, 377
  90. Stanford, S. A., Stern, D., van Breugel, W., & De Breuck, C. 2000, ApJS, 131, 185
  91. Tegmark, M. et al. 2004, Phys. Rev. D, 69, 103501
  92. —. 2006, Phys. Rev. D, 74, 123507
  93. Thompson, T. A., Quataert, E., Waxman, E., Murray, N., & Martin, C. L. 2006, ApJ, 645, 186
  94. Urry, C. M. & Padovani, P. 1995, PASP, 107, 803
  95. Veilleux, S., Kim, D.-C., & Sanders, D. B. 2002, ApJS, 143, 315
  96. Vlahakis, C., Eales, S., & Dunne, L. 2007, MNRAS, 379, 1042
  97. Voelk, H. J. 1989, A&A, 218, 67
  98. Wang, W.-H. et al. 2007, ApJ, 670, L89
  99. Weedman, D. et al. 2006a, ApJ, 653, 101
  100. Weedman, D. W., Le Floc’h, E., Higdon, S. J. U., Higdon, J. L., & Houck, J. R. 2006b, ApJ, 638, 613
  101. Weingartner, J. C. & Draine, B. T. 2001, ApJ, 548, 296
  102. Werner, M. W. et al. 2004, ApJS, 154, 1
  103. Wu, H. et al. 2005, ApJ, 632, L79
  104. Yan, L. et al. 2005, ApJ, 628, 604
  105. —. 2007, ApJ, 658, 778
  106. Yang, M., Greve, T. R., Dowell, C. D., & Borys, C. 2007, ApJ, 660, 1198
  107. Yang, M. & Phillips, T. 2007, ApJ, 662, 284
  108. Younger, J. D., Hopkins, P. F., Cox, T. J., & Hernquist, L. 2008a, ApJ, 686, 815
  109. Younger, J. D. et al. 2007, ApJ, 671, 1531
  110. —. 2008b, MNRAS, 387, 707
  111. —. 2008c, ApJ, 688, 59
  112. Yun, M. S. & Carilli, C. L. 2002, ApJ, 568, 88
  113. Yun, M. S., Reddy, N. A., & Condon, J. J. 2001, ApJ, 554, 803
This is a comment super asjknd jkasnjk adsnkj
The feedback cannot be empty
Comments 0
The feedback cannot be empty
Add comment

You’re adding your first comment!
How to quickly get a good reply:
  • Offer a constructive comment on the author work.
  • Add helpful links to code implementation or project page.