David T. Frayer

The Herschel ATLAS

The Herschel ATLAS is the largest open-time key project that will be carried out on the Herschel Space Observatory. It will survey 570 deg2 of the extragalactic sky, 4 times larger than all the other Herschel extragalactic surveys combined, in five far-infrared and submillimeter bands. We describe the survey, the complementary multiwavelength data sets that will be combined with the Herschel data, and the six major science programs we are undertaking. Using new models based on a previous submillimeter survey of galaxies, we present predictions of the properties of the ATLAS sources in other wave bands.

Mid-Infrared Spectral Diagnosis of Submillimeter Galaxies

We present deep mid-IR spectroscopy with Spitzer of 13 SMGs in the GOODS-N field.We find strong PAH emission in all of our targets, which allows us to measure mid-IR spectroscopic redshifts and place constraints on the contribution from star formation and AGN activity to the mid-IR emission. In the high-S/N composite spectrum, we find that the hot dust continuum from an AGN contributes at most 30% of the mid-IR luminosity. Individually, only 2/13 SMGs have continuum emission dominating the mid-IR luminosity; one of these SMGs, C1, remains undetected in the deep X-ray images but shows a steeply rising continuum in the mid-IR indicative of a Compton-thick AGN. We find that the mid-IR properties of SMGs are distinct from those of 24 μm–selected ULIRGs at z~2; the former are predominantly dominated by star formation, while the latter are a more heterogeneous sample with many showing significant AGN activity.We fit the IRS spectrum and the mid-IR to radio photometry of SMGs with template SEDs to determine the best estimate of the total IR luminosity from star formation. While many SMGs contain an AGN as evinced by their X-ray properties, our multiwavelength analysis shows that the total IR luminosity, L_(IR), in SMGs is dominated by star formation.We find that high-redshift SMGs lie on the relation between L_(IR) and L_(PAH,6.2) (or L_(PAH,7.7) or L_(PAH,11.3))that has been established for local starburst galaxies. This suggests that PAH luminosity can be used as a proxy for the SFR in SMGs. SMGs are consistent with being a short-lived cool phase in a massive merger where the AGN does not appear to have become strong enough to heat the dust and dominate the mid- or far-IR emission.

The Hubble deep field-north SCUBA super-map – IV. Characterizing submillimetre galaxies using deep Spitzer imaging

We present spectral energy distributions (SEDs), Spitzer colours, and infrared (IR) luminosities for 850-μm selected galaxies in the Great Observatories Origins Deep Survey Northern (GOODS-N) field. Using the deep Spitzer Legacy images and new data and reductions of the Very Large Array-Hubble Deep Field (VLA-HDF) radio data, we find statistically secure counterparts for 60 per cent (21/35) of our submillimetre (submm) sample, and identify tentative counterparts for another 12 objects. This is the largest sample of submm galaxies with statistically secure counterparts detected in the radio and with Spitzer. Half of the secure counterparts have spectroscopic redshifts, while the other half have photometric redshifts. We find that in most cases the 850-μm emission is dominated by a single 24-μm source, with a median flux density of 241 μJy, leading to a median 24-to-850-μm flux density ratio of 0.040. A composite rest-frame SED shows that the submm sources peak at longer wavelengths than those of local ultraluminous infrared galaxies (ULIRGs). Using a basic grey-body model, 850-μm selected galaxies appear to be cooler than local ULIRGs of the same luminosity. This demonstrates the strong selection effects, both locally and at high redshift, which may lead to an incomplete census of the ULIRG population. The SEDs of submm galaxies are also different from those of their high-redshift neighbours, the near-IR selected BzK galaxies, whose mid-IR-to-radio SEDs are more like those of local ULIRGs. Using 24-μm, 850-μm and 1.4-GHz observations, we fit templates that span the mid-IR through radio to derive the integrated IR luminosity (LIR) of the submm galaxies and find a median value of LIR(8–1000 μm) = 6.0 × 1012 L. By themselves, 24-μm and radio fluxes are able to predict LIR reasonably well because they are relatively insensitive to temperature. However, the submm flux by itself consistently overpredicts LIR when using spectral templates which obey the local ULIRG temperature–luminosity relation. The shorter Spitzer wavelengths sample the stellar bump at the redshifts of the submm sources, and we find that the Spitzer photometry alone provides a model-independent estimate of the redshift, σ[Δz/(1 + z)] = 0.07. The median redshift for our secure submm counterparts is 2.0. Using X-ray and mid-IR data, only 5 per cent of our secure counterparts (1/21) show strong evidence for an active galactic nucleus dominating the LIR.

The detection of a population of submillimeter-bright, strongly lensed galaxies

Through a Lens Brightly Astronomical sources detected in the submillimeter range are generally thought to be distant, dusty galaxies undergoing a vigorous burst of star formation. They can be detected because the dust absorbs the light from stars and reemits it at longer wavelengths. Their properties are still difficult to ascertain, however, because the combination of interference from dust and the low spatial resolution of submillimeter telescopes prevents further study at other wavelengths. Using data from the Herschel Space Telescope, Negrello et al. (p. 800) showed that by searching for the brightest sources in a wide enough area in the sky it was possible to detect gravitationally lensed submillimeter galaxies with nearly full efficiency. Gravitational lensing occurs when the light of an astronomical object is deflected by a foreground mass. This phenomenon increases the apparent brightness and angular size of the lensed objects, making it easier to study sources that would be otherwise too faint to probe. Data from the Herschel Space Observatory unveils distant, dusty galaxies invisible to optical telescopes. Gravitational lensing is a powerful astrophysical and cosmological probe and is particularly valuable at submillimeter wavelengths for the study of the statistical and individual properties of dusty star-forming galaxies. However, the identification of gravitational lenses is often time-intensive, involving the sifting of large volumes of imaging or spectroscopic data to find few candidates. We used early data from the Herschel Astrophysical Terahertz Large Area Survey to demonstrate that wide-area submillimeter surveys can simply and easily detect strong gravitational lensing events, with close to 100% efficiency.

DEEP SPITZER 24 μm COSMOS IMAGING. I. THE EVOLUTION OF LUMINOUS DUSTY GALAXIES—CONFRONTING THE MODELS

We present the first results obtained from the identification of ~30,000 sources in the Spitzer/24 μm observations of the COSMOS field at S_(24 μm) ≳ 80 μJy. Using accurate photometric redshifts (σ_ z ~ 0.12 at z ~ 2 for 24 μm sources with i ^+ ≳ 25 mag AB) and simple extrapolations of the number counts at faint fluxes, we resolve with unprecedented detail the buildup of the mid-infrared background across cosmic ages. We find that ~50% and ~80% of the 24 μm background intensity originate from galaxies at z ≳ 1 and z ≳ 2, respectively, supporting the scenario where highly obscured sources at very high redshifts (z ≳ 2) contribute only marginally to the cosmic infrared background. Assuming flux-limited selections at optical wavelengths, we also find that the fraction of i ^+-band sources with 24 μm detection strongly increases up to z ~ 2 as a consequence of the rapid evolution that star-forming galaxies have undergone with look-back time. Nonetheless, this rising trend shows a clear break at z ~ 1.3, probably due to k-correction effects implied by the complexity of spectral energy distributions in the mid-infrared. Finally, we compare our results with the predictions from different models of galaxy formation. We note that semianalytical formalisms currently fail to reproduce the redshift distributions observed at 24 μm. Furthermore, the simulated galaxies at S _(24 μm) > 80 μJy exhibit R–K colors much bluer than observed and the predicted K-band fluxes are systematically underestimated at z ≳ 0.5. Unless these discrepancies mainly result from an incorrect treatment of extinction in the models they may reflect an underestimate of the predicted density of high-redshift massive sources with strong ongoing star formation, which would point to more fundamental processes and/or parameters (e.g., initial mass function, critical density to form stars, feedback,...) that are still not fully controlled in the simulations. The most recent backward evolution scenarios reproduce reasonably well the flux/redshift distribution of 24 μm sources up to z ~ 3, although none of them is able to exactly match our results at all redshifts.

SPITZER MID-INFRARED SPECTROSCOPY OF INFRARED LUMINOUS GALAXIES AT z 2. I. THE SPECTRA

We present the mid-infrared spectra obtained with the Spitzer IRS for 52 sources, selected as infrared luminous, z ≳ 1 candidates in the Extragalactic First Look Survey. The sample selection criteria are f_(24 μm) ≳ 0.9 mJy, νf_ν(24 μm)/νf_ν(8 μm) ≳ 3.16, and νf_ν(24 μm)/νf_ν(0.7 μm) ≳ 10. Of the 52 spectra, 47 (90%) produced measurable redshifts based solely on the mid-IR spectral features, with 35/47 (74%) at 1.5 ≲ z ≲ 3.2. Keck spectroscopy of a subsample (17/47) agrees with the mid-IR redshift measurements. The observed spectra fall into three categories: (1) 33% (17/52) have strong PAH emission and are probably powered by star formation with total IR luminosity roughly a factor of 5 higher than the local starburst ULIRGs; (2) 33% (17/52) have only deep silicate absorption at 9.8 μm, indicatiing deeply embedded dusty systems—these data alone cannot determine the energetic nature of the heating sources in these systems; and (3) the remaining 34% are mid-IR continuum-dominated systems with weak PAH emission and/or silicate absorption. This subsample is probably AGNs. We derived monochromatic, rest-frame 5.8 μm, continuum luminosities (νL_ν), ranging from 10^(10.3) to 10^(12.6) L_☉. Our spectra have mid-IR slope α_(5-15 μm) ≳ 2.1, much redder than the median value of 1.3 for the optically selected PG quasars. From the silicate absorption feature, we estimate that roughly two-thirds of the sample have optical depth τ_(9.8 μm) > 1. Their L_(1600 A) and L_(IR) suggest that our sample is among the most luminous and most dust-enshrouded systems of its epoch. Our study has revealed a significant population of dust-enshrouded galaxies at z ~ 2, whose enormous energy output, comparable to that of quasars, is generated by AGNs, as well as starbursts. This IR-luminous population has very little overlap with submillimeter and UV-selected populations.

Molecular Gas in the z = 2.565 Submillimeter Galaxy SMM J14011+0252

We report the detection of CO (3?2) emission from the submillimeter-selected luminous galaxy SMM J14011+0252. The optical counterpart of the submillimeter source has been identified as a merger system with spectral characteristics consistent with a starburst at z=2.565. The CO emission confirms the optical identification of the submillimeter source and implies a molecular gas mass of 5?1010 h?275 M?, after correcting for a lensing amplification factor of 2.75. The large molecular gas mass and the radio emission are consistent with the starburst interpretation of the source. These results are similar to those found for SMM J02399-0136, which was the first submillimeter-selected CO source found at high redshift. The CO detections of these two high-redshift submillimeter galaxies suggest the presence of massive reservoirs of molecular gas, which is consistent with the inferred high rates of star formation (103 M? yr-1). These two systems appear to be associated with merger events which may evolve into present-day luminous elliptical galaxies.

BLAST: the far‐infrared/radio correlation in distant galaxies

We investigate the correlation between far-infrared (FIR) and radio luminosities in distant galaxies, a lynchpin of modern astronomy. We use data from the Balloon-borne Large Aperture Submillimetre Telescope (BLAST), Spitzer, the Large Apex BOlometer CamerA (LABOCA), the Very Large Array and the Giant Metre-wave Radio Telescope (GMRT) in the Extended Chandra Deep Field South (ECDFS). For a catalogue of BLAST 250-μm-selected galaxies, we remeasure the 70–870-μm flux densities at the positions of their most likely 24-μm counterparts, which have a median [interquartile] redshift of 0.74 [0.25, 1.57]. From these, we determine the monochromatic flux density ratio, q_(250)(= log_(10) [ S_(250 μm)/S_(1400 MHz)]), and the bolometric equivalent, q_(IR). At z ≈ 0.6 , where our 250-μm filter probes rest-frame 160-μm emission, we find no evolution relative to q_(160) for local galaxies. We also stack the FIR and submm images at the positions of 24-μm- and radio-selected galaxies. The difference between q_(IR) seen for 250-μm- and radio-selected galaxies suggests that star formation provides most of the IR luminosity in ≲100-μJy radio galaxies, but rather less for those in the mJy regime. For the 24-μm sample, the radio spectral index is constant across 0 < z < 3 , but q_(IR) exhibits tentative evidence of a steady decline such that q_(IR) ∝ (1 +z)^(−0.15±0.03) – significant evolution, spanning the epoch of galaxy formation, with major implications for techniques that rely on the FIR/radio correlation. We compare with model predictions and speculate that we may be seeing the increase in radio activity that gives rise to the radio background.

Herschel-ATLAS: first data release of the Science Demonstration Phase source catalogues

The Herschel Astrophysical Terahertz Large Area Survey (Herschel-ATLAS) is a survey of 550 deg2 with the Herschel Space Observatory in five far-infrared and submillimetre bands. The first data for the survey, observations of a field 4 × 4 deg2 in size, were taken during the Science Demonstration Phase (SDP), and reach a 5σ noise level of 33.5 mJy beam−1 at 250 μm. This paper describes the source extraction methods used to create the corresponding SDP catalogue, which contains 6876 sources, selected at 250 μm, within ∼14 deg2. Spectral and Photometric Imaging REciever (SPIRE) sources are extracted using a new method specifically developed for Herschel data and Photodetector Array Camera and Spectrometer (PACS) counterparts of these sources are identified using circular apertures placed at the SPIRE positions. Aperture flux densities are measured for sources identified as extended after matching to optical wavelengths. The reliability of this catalogue is also discussed, using full simulated maps at the three SPIRE bands. These show that a significant number of sources at 350 and 500 μm have undergone flux density enhancements of up to a factor of ∼2, due mainly to source confusion. Correction factors are determined for these effects. The SDP data set and corresponding catalogue will be available from http://www.h-atlas.org.

THE NATURE OF FAINT SPITZER-SELECTED DUST-OBSCURED GALAXIES

We use deep far-IR, submillimeter, radio, and X-ray imaging and mid-IR spectroscopy to explore the nature of a sample of Spitzer-selected dust-obscured galaxies (DOGs) in GOODS-N. A sample of 79 galaxies satisfy the criteria R − [ 24] > 14 (Vega) down to S_(24) > 100 μJy (median flux density S_(24) = 180 μJy). Twelve of these galaxies have IRS spectra available, which we use to measure redshifts and classify these objects as being dominated by star formation or active galactic nucleus (AGN) activity in the mid-IR. The IRS spectra and Spitzer photometric redshifts confirm that the DOGs lie in a tight redshift distribution around z ~ 2. Based on mid-IR colors, 80% of DOGs are likely dominated by star formation; the stacked X-ray emission from this subsample of DOGs is also consistent with star formation. Since only a small number of DOGs are individually detected at far-IR and submillimeter wavelengths, we use a stacking analysis to determine the average flux from these objects and plot a composite IR (8-1000 μm) spectral energy distribution (SED). The average luminosity of these star-forming DOGs is L_(IR) ~ 1 × 10^(12)L☉. We compare the average star-forming DOG to the average bright (S_(850) > 5 mJy) submillimeter galaxy (SMG); the S_(24) > 100 μJy DOGs are 3 times more numerous but 8 times less luminous in the IR. The far-IR SED shape of DOGs is similar to that of SMGs (average dust temperature of around 30 K), but DOGs have a higher mid-IR-to-far-IR flux ratio. The average star formation-dominated DOG has a star formation rate of 200 M☉ yr^(−1), which, given their space density, amounts to a contribution of 0.01 M☉yr^(−1) Mpc^(−3) (or 5%-10%) to the star formation rate density at z ~ 2.