M. Viero

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.

Over half of the far-infrared background light comes from galaxies at z ≥ 1.2

Submillimetre surveys during the past decade have discovered a population of luminous, high-redshift, dusty starburst galaxies. In the redshift range 1 ≤ z ≤ 4, these massive submillimetre galaxies go through a phase characterized by optically obscured star formation at rates several hundred times that in the local Universe. Half of the starlight from this highly energetic process is absorbed and thermally re-radiated by clouds of dust at temperatures near 30 K with spectral energy distributions peaking at 100 μm in the rest frame. At 1 ≤ z ≤ 4, the peak is redshifted to wavelengths between 200 and 500 μm. The cumulative effect of these galaxies is to yield extragalactic optical and far-infrared backgrounds with approximately equal energy densities. Since the initial detection of the far-infrared background (FIRB), higher-resolution experiments have sought to decompose this integrated radiation into the contributions from individual galaxies. Here we report the results of an extragalactic survey at 250, 350 and 500 μm. Combining our results at 500 μm with those at 24 μm, we determine that all of the FIRB comes from individual galaxies, with galaxies at z ≥ 1.2 accounting for 70% of it. As expected, at the longest wavelengths the signal is dominated by ultraluminous galaxies at z > 1.

BLAST: RESOLVING THE COSMIC SUBMILLIMETER BACKGROUND

The Balloon-borne Large Aperture Submillimeter Telescope (BLAST) has made 1 deg2, deep, confusion-limited maps at three different bands, centered on the Great Observatories Origins Deep Survey South Field. By calculating the covariance of these maps with catalogs of 24 μm sources from the Far-Infrared Deep Extragalactic Legacy Survey, we have determined that the total submillimeter intensities are 8.60 ± 0.59, 4.93 ± 0.34, and 2.27 ± 0.20 nW m–2 sr–1 at 250, 350, and 500 μm, respectively. These numbers are more precise than previous estimates of the cosmic infrared background (CIB) and are consistent with 24 μm-selected galaxies generating the full intensity of the CIB. We find that the fraction of the CIB that originates from sources at z ≥ 1.2 increases with wavelength, with 60% from high-redshift sources at 500 μm. At all BLAST wavelengths, the relative intensity of high-z sources is higher for 24 μm-faint sources than that for 24 μm-bright sources. Galaxies identified as active galactic nuclei (AGNs) by their Infrared Array Camera colors are 1.6-2.6 times brighter than the average population at 250-500 μm, consistent with what is found for X-ray-selected AGNs. BzK-selected galaxies are found to be moderately brighter than typical 24 μm-selected galaxies in the BLAST bands. These data provide high-precision constraints for models of the evolution of the number density and intensity of star-forming galaxies at high redshift.

The suppression of star formation by powerful active galactic nuclei.

The old, red stars that constitute the bulges of galaxies, and the massive black holes at their centres, are the relics of a period in cosmic history when galaxies formed stars at remarkable rates and active galactic nuclei (AGN) shone brightly as a result of accretion onto black holes. It is widely suspected, but unproved, that the tight correlation between the mass of the black hole and the mass of the stellar bulge results from the AGN quenching the surrounding star formation as it approaches its peak luminosity. X-rays trace emission from AGN unambiguously, whereas powerful star-forming galaxies are usually dust-obscured and are brightest at infrared and submillimetre wavelengths. Here we report submillimetre and X-ray observations that show that rapid star formation was common in the host galaxies of AGN when the Universe was 2–6 billion years old, but that the most vigorous star formation is not observed around black holes above an X-ray luminosity of 1044 ergs per second. This suppression of star formation in the host galaxy of a powerful AGN is a key prediction of models in which the AGN drives an outflow, expelling the interstellar medium of its host and transforming the galaxy’s properties in a brief period of cosmic time.

HerMES: deep number counts at 250 μm, 350 μm and 500 μm in the COSMOS and GOODS-N fields and the build-up of the cosmic infrared background

Aims. The Spectral and Photometric Imaging REceiver (SPIRE) onboard the Herschel space telescope has provided confusion limited maps of deep fields at 250 μm, 350 μm, and 500 μm, as part of the Herschel Multi-tiered Extragalactic Survey (HerMES). Unfortunately, due to confusion, only a small fraction of the cosmic infrared background (CIB) can be resolved into individually-detected sources. Our goal is to produce deep galaxy number counts and redshift distributions below the confusion limit at SPIRE wavelengths (~20 mJy), which we then use to place strong constraints on the origins of the cosmic infrared background and on models of galaxy evolution. Methods. We individually extracted the bright SPIRE sources (>20 mJy) in the COSMOS field with a method using the positions, the flux densities, and the redshifts of the 24 μm sources as a prior, and derived the number counts and redshift distributions of the bright SPIRE sources. For fainter SPIRE sources (<20 mJy), we reconstructed the number counts and the redshift distribution below the confusion limit using the deep 24 μm catalogs associated with photometric redshift and information provided by the stacking of these sources into the deep SPIRE maps of the GOODS-N and COSMOS fields. Finally, by integrating all these counts, we studied the contribution of the galaxies to the CIB as a function of their flux density and redshift. Results. Through stacking, we managed to reconstruct the source counts per redshift slice down to ~2 mJy in the three SPIRE bands, which lies about a factor 10 below the 5σ confusion limit. Our measurements place tight constraints on source population models. None of the pre-existing models are able to reproduce our results at better than 3-σ. Finally, we extrapolate our counts to zero flux density in order to derive an estimate of the total contribution of galaxies to the CIB, finding 10.1_(-2.3)^(+2.6) nW m^(-2) sr^(-1), 6.5_(-1.6)^(+1.7) nW m^(-2) sr^(-1), and 2.8_(-0.8)^(+0.9) nW m^(-2) sr^(-1) at 250 μm, 350 μm, and 500 μm, respectively. These values agree well with FIRAS absolute measurements, suggesting our number counts and their extrapolation are sufficient to explain the CIB. We find that half of the CIB is emitted at z = 1.04, 1.20, and 1.25, respectively. Finally, combining our results with other works, we estimate the energy budget contained in the CIB between 8 μm and 1000 μm: 26_(-3)^(+7) nW m^(-2) sr^(-1).

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.

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.

A joint analysis of BLAST 250–500 μm and LABOCA 870 μm observations in the Extended Chandra Deep Field-South

We present a joint analysis of the overlapping Balloon-borne Large Aperture Submillimetre Telescope (BLAST) 250, 350, 500 μm, and LABOCA 870 μm observations [from the LABOCA ECDFS Submm Survey (LESS) survey] of the Extended Chandra Deep Field-South. Out to z∼ 3, the BLAST filters sample near the peak wavelength of thermal far-infrared (FIR) emission from galaxies (rest-frame wavelengths ∼60–200 μm), primarily produced by dust heated through absorption in star-forming clouds. However, identifying counterparts to individual BLAST peaks is very challenging, given the large beams [full-width at half-maximum (FWHM) 36–60 arcsec]. In contrast, the ground-based 870 μm observations have a significantly smaller 19 arcsec FWHM beam, and are sensitive to higher redshifts (z∼ 1–5, and potentially beyond) due to the more favourable negative K-correction. We use the LESS data, as well as deep Spitzer and VLA imaging, to identify 118 individual sources that produce significant emission in the BLAST bands. We characterize the temperatures and FIR luminosities for a subset of 69 sources which have well-measured submillimetre (submm) spectral energy distributions (SEDs) and redshift measurements out to z∼ 3. For flux-limited sub-samples in each BLAST band, and a dust emissivity index β= 2.0, we find a median temperature T= 30 K (all bands) as well as median redshifts: z= 1.1 (interquartile range 0.2–1.9) for S250 > 40 mJy; z= 1.3 (interquartile range 0.6–2.1) for S350 > 30 mJy; and z= 1.6 (interquartile range 1.3–2.3) for S500 > 20 mJy. Taking into account the selection effects for our survey (a bias towards detecting lower-temperature galaxies), we find no evidence for evolution in the local FIR–temperature correlation out to z∼ 2.5. Comparing with star-forming galaxy SED templates, about 8 per cent of our sample appears to exhibit significant excesses in the radio and/or mid-IR, consistent with those sources harbouring active galactic nuclei (AGN). Since our statistical approach differs from most previous studies of submm galaxies, we describe the following techniques in two appendices: our ‘matched filter’ for identifying sources in the presence of point-source confusion; and our approach for identifying counterparts using likelihood ratios. This study is a direct precursor to future joint FIR/submm surveys, for which we outline a potential identification and SED measurement strategy.

HerMES: Cosmic Infrared Background Anisotropies and the Clustering of Dusty Star-Forming Galaxies

We present measurements of the auto- and cross-frequency power spectra of the cosmic infrared background (CIB) at 250, 350, and 500 μm (1200, 860, and 600 GHz) from observations totaling ~70 deg2 made with the SPIRE instrument aboard the Herschel Space Observatory. We measure a fractional anisotropy δI/I = 14% ± 4%, detecting signatures arising from the clustering of dusty star-forming galaxies in both the linear (2-halo) and nonlinear (1-halo) regimes; and that the transition from the 2- to 1-halo terms, below which power originates predominantly from multiple galaxies within dark matter halos, occurs at k θ ~ 0.10-0.12 arcmin–1 (l ~ 2160-2380), from 250 to 500 μm. New to this paper is clear evidence of a dependence of the Poisson and 1-halo power on the flux-cut level of masked sources—suggesting that some fraction of the more luminous sources occupy more massive halos as satellites, or are possibly close pairs. We measure the cross-correlation power spectra between bands, finding that bands which are farthest apart are the least correlated, as well as hints of a reduction in the correlation between bands when resolved sources are more aggressively masked. In the second part of the paper, we attempt to interpret the measurements in the framework of the halo model. With the aim of fitting simultaneously with one model the power spectra, number counts, and absolute CIB level in all bands, we find that this is achievable by invoking a luminosity-mass relationship, such that the luminosity-to-mass ratio peaks at a particular halo mass scale and declines toward lower and higher mass halos. Our best-fit model finds that the halo mass which is most efficient at hosting star formation in the redshift range of peak star-forming activity, z ~ 1-3, is log(M peak/M ☉) ~ 12.1 ± 0.5, and that the minimum halo mass to host infrared galaxies is log(M min/M ☉) ~ 10.1 ± 0.6.

Panchromatic spectral energy distributions of Herschel sources

Combining far-infrared Herschel photometry from the PACS Evolutionary Probe (PEP) and Herschel Multi-tiered Extragalactic Survey (HerMES) guaranteed time programs with ancillary datasets in the GOODS-N, GOODS-S, and COSMOS fields, it is possible to sample the 8–500 μm spectral energy distributions (SEDs) of galaxies with at least 7–10 bands. Extending to the UV, optical, and near-infrared, the number of bands increases up to 43. We reproduce the distribution of galaxies in a carefully selected restframe ten colors space, based on this rich data-set, using a superposition of multivariate Gaussian modes. We use this model to classify galaxies and build median SEDs of each class, which are then fitted with a modified version of the magphys code that combines stellar light, emission from dust heated by stars and a possible warm dust contribution heated by an active galactic nucleus (AGN). The color distribution of galaxies in each of the considered fields can be well described with the combination of 6–9 classes, spanning a large range of far- to near-infrared luminosity ratios, as well as different strength of the AGN contribution to bolometric luminosities. The defined Gaussian grouping is used to identify rare or odd sources. The zoology of outliers includes Herschel-detected ellipticals, very blue z ~ 1 Ly-break galaxies, quiescent spirals, and torus-dominated AGN with star formation. Out of these groups and outliers, a new template library is assembled, consisting of 32 SEDs describing the intrinsic scatter in the restframe UV-to-submm colors of infrared galaxies. This library is tested against L(IR) estimates with and without Herschel data included, and compared to eightother popular methods often adopted in the literature. When implementing Herschel photometry, these approaches produce L(IR) values consistent with each other within a median absolute deviation of 10–20%, the scatter being dominated more by fine tuning of the codes, rather than by the choice of SED templates. Finally, the library is used to classify 24 μm detected sources in PEP GOODS fields on the basis of AGN content, L(60)/L(100) color and L(160)/L(1.6) luminosity ratio. AGN appear to be distributed in the stellar mass (M_∗) vs. star formation rate (SFR) space along with all other galaxies, regardless of the amount of infrared luminosity they are powering, with the tendency to lie on the high SFR side of the “main sequence”. The incidence of warmer star-forming sources grows for objects with higher specific star formation rates (sSFR), and they tend to populate the “off-sequence” region of the M_∗ − SFR − z space.