The contribution of starbursts and normal galaxies to infrared luminosity functions at z < 2

We present a parameter-less approach to predict the shape of the infrared (IR) luminosity function (LF) at redshifts z < 2. It requires no tuning and relies on only three observables: (1) the redshift evolution of the stellar mass function for star-forming galaxies, (2) the evolution of the specific star formation rate (sSFR) of main-sequence galaxies, and (3) the double-Gaussian decomposition of the sSFR-distribution at fixed stellar mass into a contribution (assumed redshift- and mass-invariant) from main-sequence and starburst activity. This self-consistent and simple framework provides a powerful tool for predicting cosmological observables: observed IR LFs are successfully matched at all z < 2, suggesting a constant or only weakly redshift-dependent contribution (8-14%) of starbursts to the star formation rate density. We separate the contributions of main-sequence and starburst activity to the global IR LF at all redshifts. The luminosity threshold above which the starburst component dominates the IR LF rises from log(LIR/Lsun) = 11.4 to 12.8 over 0 < z < 2, reflecting our assumed (1+z)^2.8-evolution of sSFR in main-sequence galaxies.Comment: 7 pages, 4 figures & 1 table. Accepted for publication in ApJL. Minor typos corrected in v2 following receipt of proof

Modeling the connection between ultraviolet and infrared galaxy populations across cosmic times

Using a phenomenological approach, we self-consistently model the redshift evolution of the ultraviolet (UV) and infrared (IR) luminosity functions across cosmic time, as well as a range of observed IR properties of UV-selected galaxy population. This model is an extension of the 2SFM (2 star-formation modes) formalism, which is based on the observed "main-sequence" of star-forming galaxies, i.e. a strong correlation between their stellar mass and their star formation rate (SFR), and a secondary population of starbursts with an excess of star formation. The balance between the UV light from young, massive stars and the dust-reprocessed IR emission is modeled following the empirical relation between the attenuation (IRX for IR excess hereafter) and the stellar mass, assuming a scatter of 0.4\,dex around this relation. We obtain a good overall agreement with the measurements of the IR luminosity function up to z~3 and the UV luminosity functions up to z~6, and show that a scatter on the IRX-M relation is mandatory to reproduce these observables. We also naturally reproduce the observed, flat relation between the mean IRX and the UV luminosity at LUV>109.5 L⊙. Finally, we perform predictions of the UV properties and detectability of IR-selected samples and the vice versa, and discuss the results in the context of the UV-rest-frame and sub-millimeter surveys of the next decade

The contribution of starbursts and normal galaxies to IR luminosity functions and the molecular gas content of the Universe at z<2

We present a parameter-less approach capable of predicting the shape of the infrared luminosity function at redshifts z ≤2. It relies on three observables: (1) the redshift evolution of the stellar mass function for star-forming galaxies, (2) the evolution of the specific star formation rate of main-sequence galaxies, and (3) the double-Gaussian decomposition of the specific star formation rate distribution at fixed stellar mass into the contributions (assumed to be redshift- and mass-invariant) from main-sequence and starburst activity. Using this self-consistent and simple framework, we identify the contributions of main-sequence and starburst activity to the global infrared luminosity function and find a constant or only weakly redshift-dependent contribution (8–14%) of starbursts to the star formation rate density at z ≤2. Over the same redshift range, we also infer the evolution of the cosmic abundance of molecular gas in star-forming galaxies, based on the relations between star formation rate and molecular gas mass followed by normal and starburst galaxies

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)

GOODS-HerschelHerschel: identification of the individual galaxies responsible for the 80-290μ\mum cosmic infrared background

We propose a new method of pushing $Herschel$ to its faintest detection limits using universal trends in the redshift evolution of the far infrared over 24$\mu$m colours in the well-sampled GOODS-North field. An extension to other fields with less multi-wavelength information is presented. This method is applied here to raise the contribution of individually detected $Herschel$ sources to the cosmic infrared background (CIRB) by a factor 5 close to its peak at 250$\mu$m and more than 3 in the 350$\mu$m and 500$\mu$m bands. We produce realistic mock $Herschel$ images of the deep PACS and SPIRE images of the GOODS-North field from the GOODS-$Herschel$ Key Program and use them to quantify the confusion noise at the position of individual sources, i.e., estimate a "local confusion noise". Two methods are used to identify sources with reliable photometric accuracy extracted using 24$\mu$m prior positions. The clean index (CI), previously defined but validated here with simulations, which measures the presence of bright 24$\mu$m neighbours and the photometric accuracy index (PAI) directly extracted from the mock $Herschel$ images. After correction for completeness, thanks to our mock $Herschel$ images, individually detected sources make up as much as 54% and 60% of the CIRB in the PACS bands down to 1.1 mJy at 100$\mu$m and 2.2 mJy at 160$\mu$m and 55, 33, and 13% of the CIRB in the SPIRE bands down to 2.5, 5, and 9 mJy at 250$\mu$m, 350$\mu$m, and 500$\mu$m, respectively. The latter depths improve the detection limits of $Herschel$ by factors of 5 at 250$\mu$m, and 3 at 350$\mu$m and 500$\mu$m as compared to the standard confusion limit. Interestingly, the dominant contributors to the CIRB in all $Herschel$ bands appear to be distant siblings of the Milky Way ($z$$\sim$0.96 for $\lambda$$<$300$\mu$m) with a stellar mass of $M_{\star}$$\sim$9$\times$10$^{10}$M$_{\odot}$.Comment: 22 pages, 16 figures. Accepted for publication by Astronomy and Astrophysic

HerMES: point source catalogues from Herschel-SPIRE observations II

Key Programme on the Herschel Space Observatory. With a wedding cake survey strategy, it consists of nested fields with varying depth and area totalling ∼380 deg2. In this paper, we present deep point source catalogues extracted from Herschel-Spectral and Photometric Imaging Receiver (SPIRE) observations of all HerMES fields, except for the later addition of the 270 deg2 HerMES Large-Mode Survey (HeLMS) field. These catalogues constitute the second Data Release (DR2) made in 2013 October. A sub-set of these catalogues, which consists of bright sources extracted from Herschel-SPIRE observations completed by 2010 May 1 (covering ∼74 deg2) were released earlier in the first extensive data release in 2012 March. Two different methods are used to generate the point source catalogues, the SUSSEXTRACTOR point source extractor used in two earlier data releases (EDR and EDR2) and a new source detection and photometry method. The latter combines an iterative source detection algorithm, STARFINDER, and a De-blended SPIRE Photometry algorithm. We use end-to-end Herschel-SPIRE simulations with realistic number counts and clustering properties to characterize basic properties of the point source catalogues, such as the completeness, reliability, photometric and positional accuracy. Over 500 000 catalogue entries in HerMES fields (except HeLMS) are released to the public through the HeDAM (Herschel Database in Marseille) website (http://hedam.lam.fr/HerMES)

Satellite content and quenching of star formation in galaxy groups at z ~ 1.8

We study the properties of satellites in the environment of massive star-forming galaxies at z ~ 1.8 in the COSMOS field, using a sample of 215 galaxies on the main sequence of star formation with an average mass of ~1011M⊙. At z> 1.5, these galaxies typically trace halos of mass ≳1013M⊙. We use optical-near-infrared photometry to estimate stellar masses and star formation rates (SFR) of centrals and satellites down to ~ 6 × 109M⊙. We stack data around 215 central galaxies to statistically detect their satellite halos, finding an average of ~3 galaxies in excess of the background density. We fit the radial profiles of satellites with simple β-models, and compare their integrated properties to model predictions. We find that the total stellar mass of satellites amounts to ~68% of the central galaxy, while spectral energy distribution modeling and far-infrared photometry consistently show their total SFR to be 25-35% of the central's rate. We also see significant variation in the specific SFR of satellites within the halo with, in particular, a sharp decrease at <100 kpc. After considering different potential explanations, we conclude that this is likely an environmental signature of the hot inner halo. This effect can be explained in the first order by a simple free-fall scenario, suggesting that these low-mass environments can shut down star formation in satellites on relatively short timescales of ~0.3 Gyr

HerMES: A Statistical Measurement of the Redshift Distribution of Herschel-SPIRE Sources Using the Cross-correlation Technique

The wide-area imaging surveys with the Herschel Space Observatory at submillimeter (sub-mm) wavelengths have now resulted in catalogs of the order of one-hundred-thousand dusty, starburst galaxies. These galaxies capture an important phase of galaxy formation and evolution, but, unfortunately, the redshift distribution of these galaxies, N(z), is still mostly uncertain due to limitations associated with counterpart identification at optical wavelengths and spectroscopic follow-up. We make a statistical estimate of N(z) using a clustering analysis of sub-mm galaxies detected at each of 250, 350 and 500 μm from the Herschel Multi-tiered Extragalactic Survey centered on the Boötes field. We cross-correlate Herschel galaxies against galaxy samples at optical and near-IR wavelengths from the Sloan Digital Sky Survey, the NOAO Deep Wide Field Survey, and the Spitzer Deep Wide Field Survey. We create optical and near-IR galaxy samples based on their photometric or spectroscopic redshift distributions and test the accuracy of those redshift distributions with similar galaxy samples defined with catalogs from the Cosmological Evolution Survey (COSMOS), which has superior spectroscopic coverage. We model the clustering auto- and cross-correlations of Herschel and optical/IR galaxy samples to estimate N(z) and clustering bias factors. The S_(350) > 20 mJy galaxies have a bias factor varying with redshift as b(z) = 1.0^(+1.0)_(–0.5)(1 + z)^1.2^(+0.3)_(–0.7). This bias and the redshift dependence is broadly in agreement with galaxies that occupy dark matter halos of mass in the range of 1012 to 10^(13) M_☉. We find that galaxy selections in all three Spectral and Photometric Imaging Receiver (SPIRE) bands share a similar average redshift, with = 1.8 ± 0.2 for 250 μm selected samples, and = 1.9 ± 0.2 for both 350 and 500 μm samples, while their distributions behave differently. For 250 μm selected galaxies we find the a larger number of sources with z ≤ 1 when compared with the subsequent two SPIRE bands, with 350 and 500 μm selected SPIRE samples having peaks in N(z) at progressively higher redshifts. We compare our clustering-based N(z) results to sub-mm galaxy model predictions in the literature, and with an estimate of N(z) using a stacking analysis of COSMOS 24 μm detections

Suppressed CO emission and high G/D ratios in z=2 galaxies with sub-solar gas-phase metallicity

We study a population of significantly sub-solar enrichment galaxies at z=1.99, to investigate how molecular gas, dust and star-formation relate in low-metallicity galaxies at the peak epoch of star-formation. We target our sample with several deep ALMA and VLA datasets, and find no individual detections of CO[4-3], CO[1-0] or dust, in stark contrast to the >60% detection rate expected for solar-enrichment galaxies with these MS Halpha SFRs. We find that both low and high density molecular gas (traced by CO[1-0] and CO[4-3] respectively) are affected by the low enrichment, showing sample average (stacked) luminosity deficits >0.5-0.7 dex below expectations. This is particularly pertinent for the use of high-J CO emission as a proxy of instantaneous star-formation rate. Our individual galaxy data and stacked constraints point to a strong inverse dependence ∝Zγ of gas-to-dust ratios (G/D) and CO-to-H2 conversion factors (aco) on metallicity at z∼2, with γG/D<-2.2 and γαCO<-0.8, respectively. We quantify the importance of comparing G/D and aco vs. metallicity trends from the literature on a common, suitably normalised metallicity scale. When accounting for systematic offsets between different metallicity scales, our z∼2 constraints on these scaling relations are consistent with the corresponding relations for local galaxies. However, among those local relations, we favour those with a steep/double power-law dependence of G/D on metallicity. Finally, we discuss the implications of these findings for (a) gas mass measurements for sub-M∗ galaxies, and (b) efforts to identify the characteristic galaxy mass scale contributing most to the comoving molecular gas density at z=2