M. Sanchez-Portal

The deepest Herschel-PACS far-infrared survey: number counts and infrared luminosity functions from combined PEP/GOODS-H observations

We present results from the deepest Herschel-Photodetector Array Camera and Spectrometer (PACS) far-infrared blank field extragalactic survey, obtained by combining observations of the Great Observatories Origins Deep Survey (GOODS) fields from the PACS Evolutionary Probe (PEP) and GOODS-Herschel key programmes. We describe data reduction and theconstruction of images and catalogues. In the deepest parts of the GOODS-S field, the catalogues reach 3σ depths of 0.9, 0.6 and 1.3 mJy at 70, 100 and 160 μm, respectively, and resolve ~75% of the cosmic infrared background at 100 μm and 160 μm into individually detected sources. We use these data to estimate the PACS confusion noise, to derive the PACS number counts down to unprecedented depths, and to determine the infrared luminosity function of galaxies down to L_(IR) = 10^(11) L⊙ at z ~ 1 and L_(IR) = 10^(12) L⊙ at z ~ 2, respectively. For the infrared luminosity function of galaxies, our deep Herschel far-infrared observations are fundamental because they provide more accurate infrared luminosity estimates than those previously obtained from mid-infrared observations. Maps and source catalogues (>3σ) are now publicly released. Combined with the large wealth of multi-wavelength data available for the GOODS fields, these data provide a powerful new tool for studying galaxy evolution over a broad range of redshifts.

The first Herschel view of the mass-SFR link in high-z galaxies

Aims. We exploit deep observations of the GOODS-N field taken with PACS, the Photodetector Array Camera and Spectrometer, onboard of Herschel, as part of the PACS evolutionary probe guaranteed time (PEP), to study the link between star formation and stellar mass in galaxies to z ∼ 2. Methods. Starting from a stellar mass – selected sample of ∼4500 galaxies with mag4.5 μm < 23.0 (AB), we identify ∼350 objects with a PACS detection at 100 or 160 μ ma nd∼ 1500 with only Spitzer 24 μm counterpart. Stellar masses and total IR luminosities (LIR) are estimated by fitting the spectral energy distributions (SEDs). Results. Consistently with other Herschel results, we find that LIR based only on 24 μm data is overestimated by a median factor ∼ 1. 8a tz ∼ 2, whereas it is underestimated (with our approach) up to a factor ∼ 1. 6a t 0.5 < z < 1.0. We then exploit this calibration to correct LIR based on the MIPS/Spitzer fluxes. These results clearly show how Herschel is fundamental to constrain LIR, and hence the star formation rate (SFR), of high redshift galaxies. Using the galaxies detected with PACS (and/or MIPS), we investigate the existence and evolution of the relations between the SFR, the specific star formation rate (SSFR=SFR/mass) and the stellar mass. Moreover, in order to avoid selection effects, we also repeat this study through a stacking analysis on the PACS images to fully exploit the far-IR information also for the Herschel and Spitzer undetected subsamples. We find that the SSFR-mass relation steepens with redshift, being almost flat at z < 1.0 and reaching a slope of α = −0.50 +0.13 −0.16 at z ∼ 2, at odds with recent works based on radio-stacking analysis at the same redshift. The mean SSFR of galaxies increases with redshift, by a factor ∼15 for

The mean star formation rate of X-ray selected active galaxies and its evolution from z ~ 2.5: results from PEP-Herschel

We study relationships between star-formation rate (SFR) and the accretion luminosity and nuclear obscuration of X-ray selected active galactic nuclei (AGNs) using a combination of deep far-infrared (FIR) and X-ray data in three key extragalactic survey fields (GOODS-South, GOODS-North and COSMOS), as part of the PACS Evolutionary Probe (PEP) program. The use of three fields with differing areas and depths enables us to explore trends between the global FIR luminosity of the AGN hosts and the luminosity of the active nucleus across 4.5 orders of magnitude in AGN luminosity (LAGN) and spanning redshifts from the Local Universe to z = 2.5. Using imaging from the Herschel/PACS instrument in 2−3 bands, we combine FIR detections and stacks of undetected objects to arrive at mean fluxes for subsamples in bins of redshift and X-ray luminosity. We constrain the importance of AGN-heated dust emission in the FIR and confirm that the majority of the FIR emission of AGNs is produced by cold dust heated by star-formation in their host galaxies. We uncover characteristic trends between the mean FIR luminosity (L60) and accretion luminosity of AGNs, which depend both on LAGN and redshift. At low AGN luminosities, accretion and SFR are uncorrelated at all redshifts, consistent with a scenario where most low-luminosity AGNs are primarily fueled by secular processes in their host galaxies. At high AGN luminosities, a significant correlation is observed between L60 and LAGN, but only among AGNs at low and moderate redshifts (z 1) suggesting that the role of mergers is less important at these epochs. At all redshifts, we find essentially no relationship between L60 and nuclear obscuration across five orders of magnitude in obscuring Hydrogen column density (NH), suggesting that various mechanisms are likely to be responsible for obscuring X-rays in active galaxies. We discuss a broad scenario which can account for these trends: one in which two different modes of AGN fueling operate in the low- and high-luminosity regimes of SMBH accretion. We postulate that the dominant mode of accretion among high-luminosity AGNs evolves with redshift. Our study, as well as a body of evidence from the literature and emerging knowledge about the properties of high redshift galaxies, supports this scenario.

A Herschel view of the far-infrared properties of submillimetre galaxies

We study a sample of 61submillimetre galaxies (SMGs) selected from ground-based surveys, with known spectroscopic redshifts and observed with the Herschel Space Observatory as part of the PACS Evolutionary Probe (PEP) and the Herschel Multi-tiered Extragalactic Survey (HerMES) guaranteed time key programmes. Our study makes use of the broad far-infrared and submillimetre wavelength coverage (100−600  μm) only made possible by the combination of observations from the PACS and SPIRE instruments aboard the Herschel Space Observatory. Using a power-law temperature distribution model to derive infrared luminosities and dust temperatures, we measure a dust emissivity spectral index for SMGs of β = 2.0 ± 0.2. Our results unambiguously unveil the diversity of the SMG population. Some SMGs exhibit extreme infrared luminosities of ~10^(13) L_⊙ and relatively warm dust components, while others are fainter (a few times 10^(12) L_⊙) and are biased towards cold dust temperatures. Although at z~2 classical SMGs (>5 mJy at 850 μm) have large infrared luminosities (~10^(13) L_⊙ ), objects only selected on their submm flux densities (without any redshift informations) probe a large range in dust temperatures and infrared luminosities. The extreme infrared luminosities of some SMGs (L_IR ≳ 10^(12.7) L_⊙, 26/61 systems) imply star formation rates (SFRs) of >500 M_⊙ yr^(-1) (assuming a Chabrier IMF and no dominant AGN contribution to the FIR luminosity). Such high SFRs are difficult to reconcile with a secular mode of star formation, and may instead correspond to a merger-driven stage in the evolution of these galaxies. Another observational argument in favour of this scenario is the presence of dust temperatures warmer than that of SMGs of lower luminosities (~40 K as opposed to ~25 K), consistent with observations of local ultra-luminous infrared galaxies triggered by major mergers and with results from hydrodynamic simulations of major mergers combined with radiative transfer calculations. Moreover, we find that luminous SMGs are systematically offset from normal star-forming galaxies in the stellar mass-SFR plane, suggesting that they are undergoing starburst events with short duty cycles, compatible with the major merger scenario. On the other hand, a significant fraction of the low infrared luminosity SMGs have cold dust temperatures, are located close to the main sequence of star formation, and therefore might be evolving through a secular mode of star formation. However, the properties of this latter population, especially their dust temperature, should be treated with caution because at these luminosities SMGs are not a representative sample of the entire star-forming galaxy population.

Enhanced star formation rates in AGN hosts with respect to inactive galaxies from PEP-Herschel observations

We compare the average star formation (SF) activity in X-ray selected AGN hosts with a mass-matched control sample of inactive galaxies, including both star forming and quiescent sources, in the 0.5 10σ confidence level) in the hosts of luminous AGNs. However, when comparing to star forming galaxies only, AGN hosts are found broadly consistent with the locus of their “main sequence”. We investigate the relative far-IR luminosity distributions of active and inactive galaxies, and find a higher fraction of PACS detected, hence normal and highly star forming systems among AGN hosts. Although different interpretations are possible, we explain our findings as a consequence of a twofold AGN growth path: faint AGNs evolve through secular processes, with instantaneous AGN accretion not tightly linked to the current total SF in the host galaxy, while the luminous AGNs co-evolve with their hosts through periods of enhanced AGN activity and star formation, possibly through major mergers. While an increased SF activity with respect to inactive galaxies of similar mass is expected in the latter, we interpret the modest SF offsets measured in low-LX AGN hosts as either a) generated by non-synchronous accretion and SF histories in a merger scenario or b) due to possible connections between instantaneous SF and accretion that can be induced by smaller scale (non-major merger) mechanisms. Far-IR luminosity distributions favour the latter scenario.

Building the cosmic infrared background brick by brick with Herschel/PEP

The cosmic infrared background (CIB) includes roughly half of the energy radiated by all galaxies at all wavelengths across cosmic time, as observed at the present epoch. The PACS Evolutionary Probe (PEP) survey is exploited here to study the CIB and its redshift differential, at 70, 100 and 160 μm, where the background peaks. Combining PACS observations of the GOODS-S, GOODS-N, Lockman Hole and COSMOS areas, we define number counts spanning over more than two orders of magnitude in flux: from ∼ 1m Jy to few hundreds mJy. Stacking of 24 μm sources and P(D) statistics extend the analysis down to ∼0.2 mJy. Taking advantage of the wealth of ancillary data in PEP fields, differential number counts d 2 N/dS /dz and CIB are studied up to z = 5. Based on these counts, we discuss the effects of confusion on PACS blank field observations and provide confusion limits for the three bands considered. While most of the available backward evolution models predict the total PACS number counts with reasonable success, the consistency to redshift distributions and CIB derivatives can still be significantly improved. The new high-quality PEP data highlight the need to include redshift-dependent constraints in future modeling. The total CIB surface brightness emitted above PEP 3σ flux limits is νIν = 4.52 ± 1.18, 8.35 ± 0.95 and 9.49 ± 0.59 [nW m −2 sr −1 ] at 70, 100, and 160 μm, respectively. These values correspond to 58 ± 7% and 74± 5% of the COBE/DIRBE CIB direct measurements at 100 and 160 μm. Employing the P(D) analysis, these fractions increase to ∼65% and ∼89%. More than half of the resolved CIB was emitted at redshift z ≤ 1. The 50%-light redshifts lie at z = 0.58, 0.67 and 0.73 at the three PACS wavelengths. The distribution moves towards earlier epochs at longer wavelengths: while the 70 μ mC IB is mainly produced by z ≤ 1.0 objects, the contribution of z > 1.0 sources reaches 50% at 160 μm. Most of the CIB resolved in the three PACS bands was emitted by galaxies with infrared luminosities in the range 10 11 −10 12 L� .

Dissecting the cosmic infra-red background with Herschel/PEP

The constituents of the cosmic IR background (CIB) are studied at its peak wavelengths (100 and 160 μm) by exploiting Herschel/PACS observations of the GOODS-N, Lockman Hole, and COSMOS fields in the PACS evolutionary probe (PEP) guaranteed-time survey. The GOODS-N data reach 3σ depths of ∼3.0 mJy at 100 μ ma nd∼5.7 mJy at 160 μm. At these levels, source densities are 40 and 18 beams/source, respectively, thus hitting the confusion limit at 160 μm. Differential number counts extend from a few mJy up to 100-200 mJy, and are approximated as a double power law, with the break lying between 5 and 10 mJy. The available ancillary information allows us to split number counts into redshift bins. At z ≤ 0.5 we isolate a class of luminous sources (LIR ∼ 10 11 L� ), whose SEDs resemble late-spiral galaxies, peaking at ∼130 μm restframe and significantly colder than what is expected on the basis of pre-Herschel models. By integrating number counts over the whole covered flux range, we obtain a surface brightness of 6.36± 1.67 and 6.58± 1.62 [nW m −2 sr −1 ] at 100 and 160 μm, resolving ∼45% and ∼52% of the CIB, respectively. When stacking 24 μm sources, the inferred CIB lies within 1.1σ and 0.5σ from direct measurements in the two bands, and fractions increase to 50% and 75%. Most of this resolved CIB fraction was radiated at z ≤ 1.0, with 160 μm sources found at higher redshift than 100 μm ones.

The evolution of the dust temperatures of galaxies in the SFR–M∗ plane up to z ~ 2

We study the evolution of the dust temperature of galaxies in the SFR− M∗ plane up to z ~ 2 using far-infrared and submillimetre observations from the Herschel Space Observatory taken as part of the PACS Evolutionary Probe (PEP) and Herschel Multi-tiered Extragalactic Survey (HerMES) guaranteed time key programmes. Starting from a sample of galaxies with reliable star-formation rates (SFRs), stellar masses (M_∗) and redshift estimates, we grid the SFR− M_∗parameter space in several redshift ranges and estimate the mean dust temperature (T_(dust)) of each SFR–M_∗ − z bin. Dust temperatures are inferred using the stacked far-infrared flux densities (100–500  μm) of our SFR–M_∗ − z bins. At all redshifts, the dust temperature of galaxies smoothly increases with rest-frame infrared luminosities (L_(IR)), specific SFRs (SSFR; i.e., SFR/M_∗), and distances with respect to the main sequence (MS) of the SFR− M_∗ plane (i.e., Δlog (SSFR)_(MS) = log [SSFR(galaxy)/SSFR_(MS)(M_∗,z)]). The T_(dust) − SSFR and T_(dust) − Δlog (SSFR)_(MS) correlations are statistically much more significant than the T_(dust) − LIR one. While the slopes of these three correlations are redshift-independent, their normalisations evolve smoothly from z = 0 and z ~ 2. We convert these results into a recipe to derive T_(dust) from SFR, M_∗ and z, valid out to z ~ 2 and for the stellar mass and SFR range covered by our stacking analysis. The existence of a strong T_(dust) − Δlog (SSFR)_(MS) correlation provides us with several pieces of information on the dust and gas content of galaxies. Firstly, the slope of the T_(dust) − Δlog (SSFR)_(MS) correlation can be explained by the increase in the star-formation efficiency (SFE; SFR/M_(gas)) with Δlog (SSFR)_(MS) as found locally by molecular gas studies. Secondly, at fixed Δlog (SSFR)_(MS), the constant dust temperature observed in galaxies probing wide ranges in SFR and M_∗ can be explained by an increase or decrease in the number of star-forming regions with comparable SFE enclosed in them. And thirdly, at high redshift, the normalisation towards hotter dust temperature of the T_(dust) − Δlog (SSFR)_(MS) correlation can be explained by the decrease in the metallicities of galaxies or by the increase in the SFE of MS galaxies. All these results support the hypothesis that the conditions prevailing in the star-forming regions of MS and far-above-MS galaxies are different. MS galaxies have star-forming regions with low SFEs and thus cold dust, while galaxies situated far above the MS seem to be in a starbursting phase characterised by star-forming regions with high SFEs and thus hot dust.

The star-formation rates of 1.5 < z < 2.5 massive galaxies

The star formation rate (SFR) is a key parameter in the study of galaxy evolution. The accuracy of SFR measurements at z ∼ 2 has been questioned following a disagreement between observations and theoretical models. The latter predict SFRs at this redshift that are typically a factor 4 or more lower than the measurements. We present star-formation rates based on calorimetric measurements of the far-infrared (FIR) luminosities for massive 1.5 12.2 L� . The SFGs and AGNs tend to exhibit the same 24 μm excess. The UV SFRs are in closer agreement with the FIR-based SFRs. Using a Calzetti UV extinction correction results in a mean excess of up to 0.3 dex and a scatter of 0.35 dex from the FIR SFRs. The previous UV SFRs are thus confirmed and the mean excess, while narrowing the gap, is insufficient to explain the discrepancy between the observed SFRs and simulation predictions.

THE IMPACT OF EVOLVING INFRARED SPECTRAL ENERGY DISTRIBUTIONS OF GALAXIES ON STAR FORMATION RATE ESTIMATES

We combine Herschel-Photodetector Array Camera and Spectrometer (PACS) data from the PACS Evolutionary Probe (PEP) program with Spitzer 24 μm and 16 μm photometry and ultra deep Infrared Spectrograph (IRS) mid-infrared spectra to measure the mid- to far-infrared spectral energy distribution (SED) of 0.7 < z < 2.5 normal star-forming galaxies (SFGs) around the main sequence (the redshift-dependent relation of star formation rate (SFR) and stellar mass). Our very deep data confirm from individual far-infrared detections that z ~ 2 SFRs are overestimated if based on 24 μm fluxes and SED templates that are calibrated via local trends with luminosity. Galaxies with similar ratios of rest-frame νL_ν(8) to 8-1000 μm infrared luminosity (LIR) tend to lie along lines of constant offset from the main sequence. We explore the relation between SED shape and offset in specific star formation rate (SSFR) from the redshift-dependent main sequence. Main-sequence galaxies tend to have a similar νL_ν(8)/LIR regardless of LIR and redshift, up to z ~ 2.5, and νL_ν(8)/LIR decreases with increasing offset above the main sequence in a consistent way at the studied redshifts. We provide a redshift-independent calibration of SED templates in the range of 8-60 μm as a function of Δlog(SSFR) offset from the main sequence. Redshift dependency enters only through the evolution of the main sequence with time. Ultra deep IRS spectra match these SED trends well and verify that they are mostly due to a change in ratio of polycyclic aromatic hydrocarbon (PAH) to LIR rather than continua of hidden active galactic nuclei (AGNs). Alternatively, we discuss the dependence of νL_ν(8)/LIR on LIR. The same νL_ν(8)/LIR is reached at increasingly higher LIR at higher redshift, with shifts relative to local by 0.5 and 0.8 dex in log(LIR) at redshifts z ~ 1 and z ~ 2. Corresponding SED template calibrations are provided for use if no stellar masses are on hand. For most of those z ~ 2 SFGs that also host an AGN, the mid-infrared is dominated by the star-forming component.