The far infra-red SEDs of main sequence and starburst galaxies

We compare observed far infra-red/sub-millimetre (FIR/sub-mm) galaxy spectral energy distributions (SEDs) of massive galaxies ($M_{\star}\gtrsim10^{10}$ $h^{-1}$M$_{\odot}$) derived through a stacking analysis with predictions from a new model of galaxy formation. The FIR SEDs of the model galaxies are calculated using a self-consistent model for the absorption and re-emission of radiation by interstellar dust based on radiative transfer calculations and global energy balance arguments. Galaxies are selected based on their position on the specific star formation rate (sSFR) - stellar mass ($M_{\star}$) plane. We identify a main sequence of star-forming galaxies in the model, i.e. a well defined relationship between sSFR and $M_\star$, up to redshift $z\sim6$. The scatter of this relationship evolves such that it is generally larger at higher stellar masses and higher redshifts. There is remarkable agreement between the predicted and observed average SEDs across a broad range of redshifts ($0.5\lesssim z\lesssim4$) for galaxies on the main sequence. However, the agreement is less good for starburst galaxies at $z\gtrsim2$, selected here to have elevated sSFRs$>10\times$ the main sequence value. We find that the predicted average SEDs are robust to changing the parameters of our dust model within physically plausible values. We also show that the dust temperature evolution of main sequence galaxies in the model is driven by star formation on the main sequence being more burst-dominated at higher redshifts.Comment: 20 pages, 13 figures. Accepted to MNRA

Star formation properties of sub-mJy radio sources

We investigate the star formation properties of ~800 sources detected in one of the deepest radio surveys at 1.4 GHz. Our sample spans a wide redshift range (~0.1 - 4) and about four orders of magnitude in star formation rate (SFR). It includes both star forming galaxies (SFGs) and active galactic nuclei (AGNs), further divided into radio-quiet and radio-loud objects. We compare the SFR derived from the far infrared luminosity, as traced by Herschel, with the SFR computed from their radio emission. We find that the radio power is a good SFR tracer not only for pure SFGs but also in the host galaxies of RQ AGNs, with no significant deviation with redshift or specific SFR. Moreover, we quantify the contribution of the starburst activity in the SFGs population and the occurrence of AGNs in sources with different level of star formation. Finally we discuss the possibility of using deep radio survey as a tool to study the cosmic star formation history.Comment: 18 pages, 14 figures, 1 table (available in its entirety as ancillary data

ALMA observations of atomic carbon in z~4 dusty star-forming galaxies

We present ALMA [CI]($1-0$) (rest frequency 492 GHz) observations for a sample of 13 strongly-lensed dusty star-forming galaxies originally discovered at 1.4mm in a blank-field survey by the South Pole Telescope. We compare these new data with available [CI] observations from the literature, allowing a study of the ISM properties of $\sim 30$ extreme dusty star-forming galaxies spanning a redshift range $2 < z < 5$. Using the [CI] line as a tracer of the molecular ISM, we find a mean molecular gas mass for SPT-DSFGs of $6.6 \times 10^{10}$ M$_{\odot}$. This is in tension with gas masses derived via low-$J$ $^{12}$CO and dust masses; bringing the estimates into accordance requires either (a) an elevated CO-to-H$_2$ conversion factor for our sample of $\alpha_{\rm CO} \sim 2.5$ and a gas-to-dust ratio $\sim200$, or (b) an high carbon abundance $X_{\rm CI} \sim 7\times10^{-5}$. Using observations of a range of additional atomic and molecular lines (including [CI], [CII], and multiple transitions of CO), we use a modern Photodissociation Region code (3D-PDR) to assess the physical conditions (including the density, UV radiation field strength, and gas temperature) within the ISM of the DSFGs in our sample. We find that the ISM within our DSFGs is characterised by dense gas permeated by strong UV fields. We note that previous efforts to characterise PDR regions in DSFGs may have significantly underestimated the density of the ISM. Combined, our analysis suggests that the ISM of extreme dusty starbursts at high redshift consists of dense, carbon-rich gas not directly comparable to the ISM of starbursts in the local Universe.Comment: 21 pages, 12 figures. Accepted for publication in MNRA

The Cosmic Far-Infrared Background Buildup Since Redshift 2 at 70 and 160 microns in the COSMOS and GOODS fields

The Cosmic Far-Infrared Background (CIB) at wavelengths around 160 {\mu}m corresponds to the peak intensity of the whole Extragalactic Background Light, which is being measured with increasing accuracy. However, the build up of the CIB emission as a function of redshift, is still not well known. Our goal is to measure the CIB history at 70 {\mu}m and 160 {\mu}m at different redshifts, and provide constraints for infrared galaxy evolution models. We use complete deep Spitzer 24 {\mu}m catalogs down to about 80 {\mu}Jy, with spectroscopic and photometric redshifts identifications, from the GOODS and COSMOS deep infrared surveys covering 2 square degrees total. After cleaning the Spitzer/MIPS 70 {\mu}m and 160 {\mu}m maps from detected sources, we stacked the far-IR images at the positions of the 24 {\mu}m sources in different redshift bins. We measured the contribution of each stacked source to the total 70 and 160 {\mu}m light, and compare with model predictions and recent far-IR measurements made with Herschel/PACS on smaller fields. We have detected components of the 70 and 160 {\mu}m backgrounds in different redshift bins up to z ~ 2. The contribution to the CIB is maximum at 0.3 <= z <= 0.9 at 160{\mu}m (and z <= 0.5 at 70 {\mu}m). A total of 81% (74%) of the 70 (160) {\mu}m background was emitted at z < 1. We estimate that the AGN relative contribution to the far-IR CIB is less than about 10% at z < 1.5. We provide a comprehensive view of the CIB buildup at 24, 70, 100, 160 {\mu}m. IR galaxy models predicting a major contribution to the CIB at z < 1 are in agreement with our measurements, while our results discard other models that predict a peak of the background at higher redshifts. Our results are available online http://www.ias.u-psud.fr/irgalaxies/ .Comment: Accepted in Astronomy & Astrophysic

Clustering, host halos and environment of z∼\sim2 galaxies as a function of their physical properties

Using a sample of 25683 star-forming and 2821 passive galaxies at $z\sim2$, selected in the COSMOS field following the BzK color criterion, we study the hosting halo mass and environment of galaxies as a function of their physical properties. Spitzer and Herschel provide accurate SFR estimates for starburst galaxies. We measure the auto- and cross-correlation functions of various galaxy sub-samples and infer the properties of their hosting halos using both an HOD model and the linear bias at large scale. We find that passive and star-forming galaxies obey a similarly rising relation between the halo and stellar mass. The mean host halo mass of star forming galaxies increases with the star formation rate between 30 and 200 M$_\odot$.yr$^{-1}$, but flattens for higher values, except if we select only main-sequence galaxies. This reflects the expected transition from a regime of secular co-evolution of the halos and the galaxies to a regime of episodic starburst. We find similar large scale biases for main-sequence, passive, and starburst galaxies at equal stellar mass, suggesting that these populations live in halos of the same mass. We detect an excess of clustering on small scales for passive galaxies and showed, by measuring the large-scale bias of close pairs, that this excess is caused by a small fraction ($\sim16$) of passive galaxies being hosted by massive halos ($\sim 3 \times 10^{13}$ M$_\odot$) as satellites. Finally, extrapolating the growth of halos hosting the z$\sim$2 population, we show that M$_\star \sim 10^{10}$ M$_\odot$ galaxies at z$\sim$2 will evolve, on average, into massive (M$_\star \sim 10^{11}$ M$_\odot$), field galaxies in the local Universe and M$_\star \sim 10^{11}$ M$_\odot$ galaxies at z=2 into local, massive, group galaxies. The most massive main-sequence galaxies and close pairs of massive, passive galaxies end up in today's clusters.Comment: 18 pages, 16 figures, Accepted by A&

SPT0346-52: Negligible AGN Activity in a Compact, Hyper-starburst Galaxy at z = 5.7

We present Chandra ACIS-S and ATCA radio continuum observations of the strongly lensed dusty, star-forming galaxy SPT-S J034640-5204.9 (hereafter SPT0346-52) at $z$ = 5.656. This galaxy has also been observed with ALMA, HST, Spitzer, Herschel, APEX, and the VLT. Previous observations indicate that if the infrared (IR) emission is driven by star formation, then the inferred lensing-corrected star formation rate ($\sim$ 4500 M_{\sun} yr$^{-1}$) and star formation rate surface density $\Sigma_{\rm SFR}$ ($\sim$ 2000 M_{\sun} {yr^{-1}} {kpc^{-2}}) are both exceptionally high. It remained unclear from the previous data, however, whether a central active galactic nucleus (AGN) contributes appreciably to the IR luminosity. The {\it Chandra} upper limit shows that SPT0346-52 is consistent with being star-formation dominated in the X-ray, and any AGN contribution to the IR emission is negligible. The ATCA radio continuum upper limits are also consistent with the FIR-to-radio correlation for star-forming galaxies with no indication of an additional AGN contribution. The observed prodigious intrinsic IR luminosity of (3.6 $\pm$ 0.3) $\times$ 10$^{13}$ L_{\sun} originates almost solely from vigorous star formation activity. With an intrinsic source size of 0.61 $\pm$ 0.03 kpc, SPT0346-52 is confirmed to have one of the highest $\Sigma_{SFR}$ of any known galaxy. This high $\Sigma_{SFR}$, which approaches the Eddington limit for a radiation pressure supported starburst, may be explained by a combination of very high star formation efficiency and gas fraction.Comment: 8 pages, 6 figures, accepted for publication in Ap

Genesis of the dusty Universe: modeling submillimetre source counts

We model the evolution of IR galaxies using a phenomenological approach to match the observed source counts at different IR wavelengths. We introduce a new algorithm for reproducing source counts based on direct integration of probability distributions rather than Monte-Carlo sampling. We construct a simple model for the evolution of the luminosity function and the colour distribution of IR galaxies which utilizes a minimum number of free parameters. Moreover we analyze how each of these parameters is constrained by observational data. The model is based on pure luminosity evolution and adopts the Dale & Helou SED templates. We find that the 850um source counts and their redshift distribution depend strongly on the shape of the luminosity evolution function, but only weakly on the details of the SEDs. We derive the best-fit evolutionary model using the 850um counts and redshift distribution as constraints. Moreover our best-fit shows a flattening of the faint end of the luminosity function towards high redshifts and requires a colour evolution which implies the typical dust temperatures of objects with the same luminosities to decrease with redshift. We compare our best-fit model to observed source counts at shorter and longer wavelengths which indicates our model reproduces the 70um and 1100um source counts remarkably well, but under-produces the counts at intermediate wavelengths. Analysis reveals that the discrepancy arises at low redshifts, indicating that revision of the adopted SED library towards lower dust temperatures (at a fixed infrared luminosity) is required. This modification is equivalent to a population of cold galaxies existing at low redshifts, as also indicated by recent Herschel results, which are underrepresented in IRAS sample. We show that the modified model successfully reproduces the source counts in a wide range of IR and submm wavelengths.Comment: 21 pages, 11 figures, 2 tables. Accepted for publication in MNRAS. Supplementary information could be found at http://www.strw.leidenuniv.nl/genesis

Planck's Dusty GEMS: Gravitationally lensed high-redshift galaxies discovered with the Planck survey

We present an analysis of 11 bright far-IR/submm sources discovered through a combination of the Planck survey and follow-up Herschel-SPIRE imaging. Each source has a redshift z=2.2-3.6 obtained through a blind redshift search with EMIR at the IRAM 30-m telescope. Interferometry obtained at IRAM and the SMA, and optical/near-infrared imaging obtained at the CFHT and the VLT reveal morphologies consistent with strongly gravitationally lensed sources. Additional photometry was obtained with JCMT/SCUBA-2 and IRAM/GISMO at 850 um and 2 mm, respectively. All objects are bright, isolated point sources in the 18 arcsec beam of SPIRE at 250 um, with spectral energy distributions peaking either near the 350 um or the 500 um bands of SPIRE, and with apparent far-infrared luminosities of up to 3x10^14 L_sun. Their morphologies and sizes, CO line widths and luminosities, dust temperatures, and far-infrared luminosities provide additional empirical evidence that these are strongly gravitationally lensed high-redshift galaxies. We discuss their dust masses and temperatures, and use additional WISE 22-um photometry and template fitting to rule out a significant contribution of AGN heating to the total infrared luminosity. Six sources are detected in FIRST at 1.4 GHz. Four have flux densities brighter than expected from the local far-infrared-radio correlation, but in the range previously found for high-z submm galaxies, one has a deficit of FIR emission, and 6 are consistent with the local correlation. The global dust-to-gas ratios and star-formation efficiencies of our sources are predominantly in the range expected from massive, metal-rich, intense, high-redshift starbursts. An extensive multi-wavelength follow-up programme is being carried out to further characterize these sources and the intense star-formation within them.Comment: A&A accepte

The Rest-Frame Submillimeter Spectrum of High-Redshift, Dusty, Star-Forming Galaxies

We present the average rest-frame spectrum of high-redshift dusty, star-forming galaxies from 250-770GHz. This spectrum was constructed by stacking ALMA 3mm spectra of 22 such sources discovered by the South Pole Telescope and spanning z=2.0-5.7. In addition to multiple bright spectral features of 12CO, [CI], and H2O, we also detect several faint transitions of 13CO, HCN, HNC, HCO+, and CN, and use the observed line strengths to characterize the typical properties of the interstellar medium of these high-redshift starburst galaxies. We find that the 13CO brightness in these objects is comparable to that of the only other z>2 star-forming galaxy in which 13CO has been observed. We show that the emission from the high-critical density molecules HCN, HNC, HCO+, and CN is consistent with a warm, dense medium with T_kin ~ 55K and n_H2 >~ 10^5.5 cm^-3. High molecular hydrogen densities are required to reproduce the observed line ratios, and we demonstrate that alternatives to purely collisional excitation are unlikely to be significant for the bulk of these systems. We quantify the average emission from several species with no individually detected transitions, and find emission from the hydride CH and the linear molecule CCH for the first time at high redshift, indicating that these molecules may be powerful probes of interstellar chemistry in high-redshift systems. These observations represent the first constraints on many molecular species with rest-frame transitions from 0.4-1.2mm in star-forming systems at high redshift, and will be invaluable in making effective use of ALMA in full science operations.Comment: 19 pages, 10 figures (2 in appendices); accepted for publication in Ap

Modelling Galaxy and AGN Evolution in the IR: Black Hole Accretion versus Star-Formation Activity

We present a new backward evolution model for galaxies and AGNs in the infrared (IR). What is new in this model is the separate study of the evolutionary properties of the different IR populations (i.e. spiral galaxies, starburst galaxies, low-luminosity AGNs, "unobscured" type 1 AGNs and "obscured" type 2 AGNs) defined through a detailed analysis of the spectral energy distributions (SEDs) of large samples of IR selected sources. The evolutionary parameters have been constrained by means of all the available observables from surveys in the mid- and far-IR (source counts, redshift and luminosity distributions, luminosity functions). By decomposing the SEDs representative of the three AGN classes into three distinct components (a stellar component emitting most of its power in the optical/near-IR, an AGN component due to hot dust heated by the central black hole peaking in the mid-IR, and a starburst component dominating the far-IR spectrum) we have disentangled the AGN contribution to the monochromatic and total IR luminosity emitted by the different populations considered in our model from that due to star-formation activity. We have then obtained an estimate of the total IR luminosity density (and star-formation density - SFD - produced by IR galaxies) and the first ever estimate of the black hole mass accretion density (BHAR) from the IR. The derived evolution of the BHAR is in agreement with estimates from X-rays, though the BHAR values we derive from IR are slightly higher than the X-ray ones. Finally, we have simulated source counts, redshift distributions and SFD and BHAR that we expect to obtain with the future cosmological Surveys in the mid-/far-IR that will be performed with JWST-MIRI and SPICA-SAFARI.Comment: 19 pages, 15 figures, 3 tables. Accepted for publication in MNRA