J. L. Wardlow

A dust-obscured massive maximum-starburst galaxy at a redshift of 6.34

Massive present-day early-type (elliptical and lenticular) galaxies probably gained the bulk of their stellar mass and heavy elements through intense, dust-enshrouded starbursts—that is, increased rates of star formation—in the most massive dark-matter haloes at early epochs. However, it remains unknown how soon after the Big Bang massive starburst progenitors exist. The measured redshift (z) distribution of dusty, massive starbursts has long been suspected to be biased low in z owing to selection effects, as confirmed by recent findings of systems with redshifts as high as ∼5 (refs 2–4). Here we report the identification of a massive starburst galaxy at z = 6.34 through a submillimetre colour-selection technique. We unambiguously determined the redshift from a suite of molecular and atomic fine-structure cooling lines. These measurements reveal a hundred billion solar masses of highly excited, chemically evolved interstellar medium in this galaxy, which constitutes at least 40 per cent of the baryonic mass. A ‘maximum starburst’ converts the gas into stars at a rate more than 2,000 times that of the Milky Way, a rate among the highest observed at any epoch. Despite the overall downturn in cosmic star formation towards the highest redshifts, it seems that environments mature enough to form the most massive, intense starbursts existed at least as early as 880 million years after the Big Bang.

The LABOCA survey of the Extended Chandra Deep Field-South: a photometric redshift survey of submillimetre galaxies

We derive photometric redshifts from 17-band optical to mid-infrared photometry of 78 robust radio, 24-mu m and Spitzer IRAC counterparts to 72 of the 126 submillimetre galaxies (SMGs) selected at 870 mu m by LABOCA observations in the Extended Chandra Deep Field-South (ECDF-S). We test the photometric redshifts of the SMGs against the extensive archival spectroscopy in the ECDF-S. The median photometric redshift of identified SMGs is z = 2.2 +/- 0.1, the standard deviation is sigma(z) = 0.9 and we identify 11 (similar to 15 per cent) high-redshift (z >= 3) SMGs. A statistical analysis of sources in the error circles of unidentified SMGs identifies a population of possible counterparts with a redshift distribution peaking at z = 2.5 +/- 0.2, which likely comprises similar to 60 per cent of the unidentified SMGs. This confirms that the bulk of the undetected SMGs are coeval with those detected in the radio/mid-infrared. We conclude that at most similar to 15 per cent of all the SMGs are below the flux limits of our IRAC observations and thus may lie at z greater than or similar to 3 and hence at most similar to 30 per cent of all SMGs have z greater than or similar to 3. We estimate that the full S(870 mu m) > 4mJy SMG population has a median redshift of 2.5 +/- 0.5. In contrast to previous suggestions, we find no significant correlation between submillimetre flux and redshift. The median stellar mass of the SMGs derived from spectral energy distribution fitting is (9.1 +/- 0.5) x 10(10)M(circle dot) although we caution that the uncertainty in the star formation histories results in a factor of similar to 5 uncertainty in these stellarmasses. Using a single temperature modified blackbody fit with beta = 1.5, the median characteristic dust temperature of SMGs is 37.4 +/- 1.4K. The infrared luminosity function shows that SMGs at z = 2-3 typically have higher far-infrared luminosities and luminosity density than those at z = 1-2. This is mirrored in the evolution of the star formation rate density (SFRD) for SMGs which peaks at z similar to 2. The maximum contribution of bright SMGs to the global SFRD (similar to 5 per cent for SMGs with S(870 mu m) greater than or similar to 4mJy or similar to 50 per cent extrapolated to SMGs with S(870 mu m) > 1mJy) also occurs at z similar to 2.

An ALMA survey of sub-millimetre galaxies in the extended chandra deep field south: The far-infrared properties of SMGs

We exploit Atacama Large Millimeter Array (ALMA) 870 mu m observations of sub-millimetre sources in the Extended Chandra Deep Field South to investigate the far-infrared properties of high-redshift sub-millimetre galaxies (SMGs). Using the precisely located 870 mu m ALMA positions of 99 SMGs, together with 24 mu m and radio imaging, we deblend the Herschel/SPIRE imaging to extract their far-infrared fluxes and colours. The median redshifts for ALMA LESS (ALESS) SMGs which are detected in at least two SPIRE bands increases with wavelength of the peak in their spectral energy distributions (SEDs), with z = 2.3 +/- 0.2, 2.5 +/- 0.3 and 3.5 +/- 0.5 for the 250, 350 and 500 mu m peakers, respectively. 34 ALESS SMGs do not have a >3 sigma counterpart at 250, 350 or 500 mu m. These galaxies have a median photometric redshift derived from the rest-frame UV-mid-infrared SEDs of z = 3.3 +/- 0.5, which is higher than the full ALESS SMG sample; z = 2.5 +/- 0.2. We estimate the far-infrared luminosities and characteristic dust temperature of each SMG, deriving L-IR = (3.0 +/- 0.3) x 10(12) L-circle dot (SFR = 300 +/- 30 M-circle dot yr(-1)) and T-d = 32 +/- 1 K. The characteristic dust temperature of these high-redshift SMGs is Delta T-d = 3-5K lower than comparably luminous galaxies at z = 0, reflecting the more extended star formation in these systems. We show that the contribution of S-870 mu m >= 1 mJy SMGs to the cosmic star formation budget is 20 per cent of the total over the redshift range z similar to 1-4. Adopting an appropriate gas-to-dust ratio, we estimate a typical molecular mass of the ALESS SMGs of M-H2 = (4.2 +/- 0.4) x 10(10) M-circle dot. Finally, we show that SMGs with S-870 mu m > 1 mJy (L-IR greater than or similar to 10(12) L-circle dot) contain similar to 10 per cent of the z similar to 2 volume-averaged H-2 mass density.

An ALMA survey of submillimetre galaxies in the Extended Chandra Deep Field South: high-resolution 870 μm source counts

We report the first counts of faint submillimetre galaxies (SMGs) in the 870-mu m band derived from arcsecond-resolution observations with the Atacama Large Millimeter Array (ALMA). We have used ALMA to map a sample of 122 870-mu m-selected submillimetre sources drawn from the 0 degrees.5x0 degrees.5 the Large Apex BOlometer CAmera (LABOCA) Extended Chandra Deep Field South submillimetre survey (LESS). These ALMA maps have an average depth of sigma 870(mu m) similar to 0.4 mJy, some approximately three times deeper than the original LABOCA survey and critically the angular resolution is more than an order of magnitude higher, FWHM of similar to 1.5 arcsec compared to similar to 19 arcsec for the LABOCA discovery map. This combination of sensitivity and resolution allows us to precisely pinpoint the SMGs contributing to the submillimetre sources from the LABOCA map, free from the effects of confusion. We show that our ALMA-derived SMG counts broadly agree with the submillimetre source counts from previous, lower resolution single-dish surveys, demonstrating that the bulk of the submillimetre sources are not caused by blending of unresolved SMGs. The difficulty which well-constrained theoretical models have in reproducing the high surface densities of SMGs, thus remains. However, our observations do show that all of the very brightest sources in the LESS sample, S-870 (mu m) greater than or similar to 12 mJy, comprise emission from multiple, fainter SMGs, each with 870-mu m fluxes of less than or similar to 9 mJy. This implies a natural limit to the star formation rate in SMGs of less than or similar to 10(3) M-circle dot yr(-1), which in turn suggests that the space densities of z > 1 galaxies with gas masses in excess of similar to 5 x 10(10) M-circle dot is <10(-5) Mpc(-3). We also discuss the influence of this blending on the identification and characterization of the SMG counterparts to these bright submillimetre sources and suggest that it may be responsible for previous claims that they lie at higher redshifts than fainter SMGs.

The LABOCA Survey of the Extended Chandra Deep Field South : Clustering of submillimetre galaxies

We present a measurement of the spatial clustering of submillimetre galaxies (SMGs) at z= 1-3. Using data from the 870μm Large APEX Bolometer Camera (LABOCA) submillimetre survey of the Extended Chandra Deep Field-South, we employ a novel technique to measure the cross-correlation between SMGs and galaxies, accounting for the full probability distributions for photometric redshifts of the galaxies. From the observed projected two-point cross-correlation function we derive the linear bias and characteristic dark matter halo masses for the SMGs. We detect clustering in the cross-correlation between SMGs and galaxies at the >4σ level. Accounting for the clustering of galaxies from their autocorrelation function, we estimate an autocorrelation length for SMGs of r o = 7.7 -2.3 +1.8 h -1 Mpc assuming a power-law slope γ= 1.8, and derive a corresponding dark matter halo mass of log(M halo[h -1M ⊙]) = 12.8 -0.5 +0.3. Based on the evolution of dark matter haloes derived from simulations, we show that that the z= 0 descendants of SMGs are typically massive (~2-3L *) elliptical galaxies residing in moderate- to high-mass groups (log(M halo[h -1M ⊙]) = 13.3 -0.5 +0.3). From the observed clustering we estimate an SMG lifetime of ~100Myr, consistent with lifetimes derived from gas consumption times and star formation time-scales, although with considerable uncertainties. The clustering of SMGs at z~ 2 is consistent with measurements for optically selected quasi-stellar objects (QSOs), supporting evolutionary scenarios in which powerful starbursts and QSOs occur in the same systems. Given that SMGs reside in haloes of characteristic mass ~6 × 10 12h -1M ⊙, we demonstrate that the redshift distribution of SMGs can be described remarkably well by the combination of two effects: the cosmological growth of structure and the evolution of the molecular gas fraction in galaxies. We conclude that the powerful starbursts in SMGs likely represent a short-lived but universal phase in massive galaxy evolution, associated with the transition between cold gas-rich, star-forming galaxies and passively evolving systems.

Herschel-atlas galaxy counts and high-redshift luminosity functions : The formation of massive early-type galaxies

Exploiting the Herschel Astrophysical Terahertz Large Area Survey Science Demonstration Phase survey data, we have determined the luminosity functions (LFs) at rest-frame wavelengths of 100 and 250 μm and at several redshifts z gsim 1, for bright submillimeter galaxies with star formation rates (SFRs) gsim 100 M ☉ yr–1. We find that the evolution of the comoving LF is strong up to z ≈ 2.5, and slows down at higher redshifts. From the LFs and the information on halo masses inferred from clustering analysis, we derived an average relation between SFR and halo mass (and its scatter). We also infer that the timescale of the main episode of dust-enshrouded star formation in massive halos (M H gsim 3 × 1012 M ☉) amounts to ~7 × 108 yr. Given the SFRs, which are in the range of 102-103 M ☉ yr–1, this timescale implies final stellar masses of the order of 1011-1012 M ☉. The corresponding stellar mass function matches the observed mass function of passively evolving galaxies at z gsim 1. The comparison of the statistics for submillimeter and UV-selected galaxies suggests that the dust-free, UV bright phase is gsim 102 times shorter than the submillimeter bright phase, implying that the dust must form soon after the onset of star formation. Using a single reference spectral energy distribution (SED; the one of the z ≈ 2.3 galaxy SMM J2135-0102), our simple physical model is able to reproduce not only the LFs at different redshifts >1 but also the counts at wavelengths ranging from 250 μm to ≈1 mm. Owing to the steepness of the counts and their relatively broad frequency range, this result suggests that the dispersion of submillimeter SEDs of z > 1 galaxies around the reference one is rather small.

An ALMA Survey of Submillimeter Galaxies in the Extended Chandra Deep Field South: The Redshift Distribution and Evolution of Submillimeter Galaxies

We present the first photometric redshift distribution for a large sample of 870 mu m submillimeter galaxies (SMGs) with robust identifications based on observations with ALMA. In our analysis we consider 96 SMGs in the Extended Chandra Deep Field South, 77 of which have 4-19 band photometry. We model the SEDs for these 77 SMGs, deriving a median photometric redshift of z(phot) = 2.3 +/- 0.1. The remaining 19 SMGs have insufficient photometry to derive photometric redshifts, but a stacking analysis of Herschel observations confirms they are not spurious. Assuming that these SMGs have an absolute H-band magnitude distribution comparable to that of a complete sample of z similar to 1-2 SMGs, we demonstrate that they lie at slightly higher redshifts, raising the median redshift for SMGs to zphot = 2.5 +/- 0.2. Critically we show that the proportion of galaxies undergoing an SMG-like phase at z >= 3 is at most 35% +/- 5% of the total population. We derive a median stellar mass of M star = (8 +/- 1) x 10(10) M circle dot, although there are systematic uncertainties of up to 5 x for individual sources. Assuming that the star formation activity in SMGs has a timescale of similar to 100 Myr, we show that their descendants at z similar to 0 would have a space density and MH distribution that are in good agreement with those of local ellipticals. In addition, the inferred mass-weighted ages of the local ellipticals broadly agree with the look-back times of the SMG events. Taken together, these results are consistent with a simple model that identifies SMGs as events that form most of the stars seen in the majority of luminous elliptical galaxies at the present day.

Gravitational Lens Models Based on Submillimeter Array Imaging of Herschel-selected Strongly Lensed Sub-millimeter Galaxies at z > 1.5

Strong gravitational lenses are now being routinely discovered in wide-field surveys at (sub-)millimeter wavelengths. We present Submillimeter Array (SMA) high-spatial resolution imaging and Gemini-South and Multiple Mirror Telescope optical spectroscopy of strong lens candidates discovered in the two widest extragalactic surveys conducted by the Herschel Space Observatory: the Herschel-Astrophysical Terahertz Large Area Survey (H-ATLAS) and the Herschel Multi-tiered Extragalactic Survey (HerMES). From a sample of 30 Herschel sources with S 500 > 100 mJy, 21 are strongly lensed (i.e., multiply imaged), 4 are moderately lensed (i.e., singly imaged), and the remainder require additional data to determine their lensing status. We apply a visibility-plane lens modeling technique to the SMA data to recover information about the masses of the lenses as well as the intrinsic (i.e., unlensed) sizes (r half) and far-infrared luminosities (L FIR) of the lensed submillimeter galaxies (SMGs). The sample of lenses comprises primarily isolated massive galaxies, but includes some groups and clusters as well. Several of the lenses are located at z lens > 0.7, a redshift regime that is inaccessible to lens searches based on Sloan Digital Sky Survey spectroscopy. The lensed SMGs are amplified by factors that are significantly below statistical model predictions given the 500 μm flux densities of our sample. We speculate that this may reflect a deficiency in our understanding of the intrinsic sizes and luminosities of the brightest SMGs. The lensed SMGs span nearly one decade in L FIR (median L FIR = 7.9 × 1012 L ☉) and two decades in FIR luminosity surface density (median ΣFIR = 6.0 × 1011 L ☉ kpc–2). The strong lenses in this sample and others identified via (sub-)mm surveys will provide a wealth of information regarding the astrophysics of galaxy formation and evolution over a wide range in redshift.

A submillimetre galaxy at z=4.76 in the LABOCA survey of the Extended Chandra Deep Field South

We report on the identification of the highest redshift submillimetre-selected source currently known LESS J033229.4−275619. This source was detected in the Large Apex Bolometer Camera (LABOCA) Extended Chandra Deep Field-South (ECDF-S) Submillimetre Survey (LESS), a sensitive 870-μm survey (σ_(870 μm)∼ 1.2 mJy) of the full 30 × 30 arcmin_2 ECDF-S with the LABOCA on the Atacama Pathfinder Experiment telescope. The submillimetre emission is identified with a radio counterpart for which optical spectroscopy provides a redshift of z= 4.76 . We show that the bolometric emission is dominated by a starburst with a star formation rate of ∼1000 M_⊙ yr^(−1), although we also identify a moderate luminosity active galactic nucleus (AGN) in this galaxy. Thus it has characteristics similar to those of z∼ 2 submillimetre galaxies (SMGs), with a mix of starburst and obscured AGN signatures. This demonstrates that ultraluminous starburst activity is not just restricted to the hosts of the most luminous (and hence rare) quasi-stellar objects at z∼ 5 , but was also occurring in less extreme galaxies at a time when the Universe was less than 10 per cent of its current age. Assuming that we are seeing the major phase of star formation in this galaxy, then we demonstrate that it would be identified as a luminous distant red galaxy at z∼ 3 and that the current estimate of the space density of z > 4 SMGs is only sufficient to produce ≳10 per cent of the luminous red galaxy population at these early times. However, this leaves open the possibility that some of these galaxies formed through less intense, but more extended star formation events. If the progenitors of all of the luminous red galaxies at z∼ 3 go through an ultraluminous starburst at z≳ 4 then the required volume density of z > 4 SMGs will exceed that predicted by current galaxy formation models by more than an order of magnitude.

A redshift survey of Herschel far-infrared selected starbursts and implications for obscured star formation

We present Keck spectroscopic observations and redshifts for a sample of 767 Herschel-SPIRE selected galaxies (HSGs) at 250, 350, and 500 μm, taken with the Keck I Low Resolution Imaging Spectrometer and the Keck II DEep Imaging Multi-Object Spectrograph. The redshift distribution of these SPIRE sources from the Herschel Multitiered Extragalactic Survey peaks at z = 0.85, with 731 sources at z < 2 and a tail of sources out to z ~ 5. We measure more significant disagreement between photometric and spectroscopic redshifts (〈Δz/(1 + z_(spec))〉 = 0.29) than is seen in non-infrared selected samples, likely due to enhanced star formation rates and dust obscuration in infrared-selected galaxies. The infrared data are used to directly measure integrated infrared luminosities and dust temperatures independent of radio or 24 μm flux densities. By probing the dust spectral energy distribution (SED) at its peak, we estimate that the vast majority (72%-83%) of z < 2 Herschel-selected galaxies would drop out of traditional submillimeter surveys at 0.85-1 mm. We find that dust temperature traces infrared luminosity, due in part to the SPIRE wavelength selection biases, and partially from physical effects. As a result, we measure no significant trend in SPIRE color with redshift; if dust temperature were independent of luminosity or redshift, a trend in SPIRE color would be expected. Composite infrared SEDs are constructed as a function of infrared luminosity, showing the increase in dust temperature with luminosity, and subtle change in near-infrared and mid-infrared spectral properties. Moderate evolution in the far-infrared (FIR)/radio correlation is measured for this partially radio-selected sample, with q_(IR) ∝(1 + z)^(–0.30±0.02) at z < 2. We estimate the luminosity function and implied star formation rate density contribution of HSGs at z < 1.6 and find overall agreement with work based on 24 μm extrapolations of the LIRG, ULIRG, and total infrared contributions. This work significantly increased the number of spectroscopically confirmed infrared-luminous galaxies at z » 0 and demonstrates the growing importance of dusty starbursts for galaxy evolution studies and the build-up of stellar mass throughout cosmic time.