Source author record

F. Pozzi

F. Pozzi appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

93works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

93 published item(s)

preprint2025arXiv

Active galactic nuclei-heated dust revealed in "little red dots"

Little red dots (LRDs) are a puzzling population of extragalactic sources whose origin is highly debated. In this {work}, we performed a comprehensive stacking analysis of NIRCam, MIRI, and ALMA images of a large and homogeneously selected sample of LRDs from multiple JWST Legacy fields. We report clear evidence of hot-dust emission in the median stacked spectral energy distribution (SED) that features a rising near-infrared continuum up to rest-frame $λ_{\rm rest}$$\sim$ 3$μ$m, which is best explained by a standard dusty active galactic nucleus (AGN) structure. Although LRDs are likely to be a heterogeneous population, our findings suggest that most ($\gtrsim$50 %) LRDs show AGN-heated dust emission, regardless of whether the optical and ultraviolet (UV) continua are stellar or AGN-dominated. In either case, the best-fit dusty-AGN SED, combined with the lack of X-ray detection in the deep Chandra stacks, suggests that Compton-thick ($N_{\rm H}$$>$3$\times$10$^{24}$ cm$^{-2}$) gas obscuration is common, and likely confined within the dust sublimation radius ($R$$_{\rm sub}$$\sim$0.1 pc). Therefore, we argue that AGN-heated dust does not directly obscure either the optical-UV continuum or the broad-line region emission, in order to explain the observed blue UV slopes and prominent Balmer features. While a gas-dust displacement is in line with several models, the formation scenario (in-situ or ex-situ) of this pre-enriched hot dust remains unclear.

preprint2022arXiv

A new estimate of the cosmic star formation density from a radio-selected sample, and the contribution of $H$-dark galaxies at $z \geq 3$

The Star Formation Rate Density (SFRD) history of the Universe is well constrained up to redshift $z \sim 2$. At earlier cosmic epochs, the picture has been largely inferred from UV-selected galaxies (e.g. Lyman-break galaxies, LBGs). However, LBGs' inferred SFRs strongly depend on the assumed dust extinction correction, which is not well-constrained at high-$z$, while observations in the radio domain are not affected by this issue. In this work we measure the SFRD from a 1.4 GHz-selected sample of $\sim$600 galaxies in the GOODS-N field up to redshift $\sim 3.5$. We take into account the contribution of Active Galactic Nuclei from the Infrared-Radio correlation. We measure the radio luminosity function, fitted with a modified Schechter function, and derive the SFRD. The cosmic SFRD shows a rise up to $z \sim 2$ and then an almost flat plateau up to $z \sim 3.5$. Our SFRD is in agreement with the ones from other FIR/radio surveys and a factor 2 higher than those from LBG samples. We also estimate that galaxies lacking a counterpart in the HST/WFC3 H-band ($H$-dark) make up $\sim 25\%$ of the $ϕ$-integrated SFRD relative to the full sample at z $\sim 3.2$, and up to $58\%$ relative to LBG samples.

preprint2022arXiv

ALPINE: A Large Survey to Understand Teenage Galaxies

A multiwavelength study of galaxies is important to understand their formation and evolution. Only in the recent past, thanks to the Atacama Large (Sub) Millimeter Array (ALMA), were we able to study the far-infrared (IR) properties of galaxies at high redshifts. In this article, we summarize recent research highlights and their significance to our understanding of early galaxy evolution from the ALPINE survey, a large program with ALMA to observe the dust continuum and 158um C+ emission of normal star-forming galaxies at z = 4-6. Combined with ancillary data at UV through near-IR wavelengths, ALPINE provides the currently largest multiwavelength sample of post-reionization galaxies and has advanced our understanding of (i) the demographics of C+ emission; (ii) the relation of star formation and C+ emission; (iii) the gas content; (iv) outflows and enrichment of the intergalactic medium; and (v) the kinematics, emergence of disks, and merger rates in galaxies at z > 4. ALPINE builds the basis for more detailed measurements with the next generation of telescopes, and places itself as an important post-reionization baseline sample to allow a continuous study of galaxies over 13 billion years of cosmic time.

preprint2022arXiv

IR characteristic emission and dust properties of star-forming galaxies at 4.5 $<$ z $<$ 6.2

The luminosity functions at z < 4 - 5 suggest that most galaxies have a relatively low stellar mass (logM_star = 10) and a low dust attenuation (A_FUV = 1.0). The physical properties of these objects are quite homogeneous. We used an approach where we combined their rest-frame far-infrared and submillimeter emissions and utilized the universe and the redshift as a spectrograph to increase the amount of information in a collective way. From a subsample of 27 ALMA-detected galaxies at z > 4.5, we built an infrared spectral energy distribution composite template. It was used to fit, with CIGALE, the 105 galaxies (detections and upper limits) in the sample from the FUV to the FIR. The derived physical parameters provide information to decipher the nature of the dust cycle and of the stellar populations in these galaxies. The derived IR composite template is consistent with the galaxies in the studied sample. A delayed star formation history with tau_main = 500 Myrs is slightly favored by the statistical analysis as compared to a delayed with a final burst or a continuous star formation history. The position of the sample in the star formation rate (SFR)- M_star diagram is consistent with previous papers. The redshift evolution of the log M_star versus A_FUV relation is in agreement with evolution in the redshift of this relation. This evolution is necessary to explain the cosmic evolution of the average dust attenuation of galaxies. Evolution is also observed in the L_dust/ L_FUV (IRX) versus UV slope beta_FUV diagram: younger galaxies have bluer beta_FUV. We modeled the shift of galaxies in the IRX versus the beta_FUV diagram with the mass-weighted age as a free parameter, and we provide an equation to make predictions.

preprint2022arXiv

The ALMA REBELS Survey: Average [CII] $158\,\rm{μm}$ sizes of Star-Forming Galaxies from $z\sim 7$ to $z\sim 4$

We present the average [CII] $158\,\rm{μm}$ emission line sizes of UV-bright star-forming galaxies at $z\sim7$. Our results are derived from a stacking analysis of [CII] $158\,\rm{μm}$ emission lines and dust continua observed by ALMA, taking advantage of the large program Reionization Era Bright Emission Line Survey (REBELS). We find that the average [CII] emission at $z\sim7$ has an effective radius $r_e$ of $2.2\pm0.2\,\rm{kpc}$. It is $\gtrsim2\times$ larger than the dust continuum and the rest-frame UV emission, in agreement with recently reported measurements for $z\lesssim6$ galaxies. Additionally, we compared the average [CII] size with $4<z<6$ galaxies observed by the ALMA Large Program to INvestigate [CII] at Early times (ALPINE). By analysing [CII] sizes of $4<z<6$ galaxies in two redshift bins, we find an average [CII] size of $r_{\rm e}=2.2\pm0.2\,\rm{kpc}$ and $r_{\rm e}=2.5\pm0.2\,\rm{kpc}$ for $z\sim5.5$ and $z\sim4.5$ galaxies, respectively. These measurements show that star-forming galaxies, on average, show no evolution in the size of the [CII] $158\,{\rm μm}$ emitting regions at redshift between $z\sim7$ and $z\sim4$. This finding suggest that the star-forming galaxies could be morphologically dominated by gas over a wide redshift range.

preprint2022arXiv

The molecular gas properties in local Seyfert 2 galaxies

We present a multi-wavelength study of the molecular gas properties of a sample of local Seyfert 2 galaxies to assess if, and to what extent, the presence of an active galactic nucleus (AGN) can affect the Interstellar Medium (ISM) properties in a sample of 33 local Seyfert 2 galaxies. We compare the molecular gas content (MH2), derived from new and archival low-J CO line measurements of a sample of AGN and a control sample of star-forming galaxies (SFGs). Both the AGN and the control sample are characterised in terms of host-galaxy properties (e.g., stellar and dust masses, Mstar and Mdust, respectively; and star formation rate, SFR). We also investigate the effect of AGN activity onto the emission of polycyclic aromatic hydrocarbon (PAH) molecules in the mid-infrared (MIR), a waveband where the dust-reprocessed emission from the obscured AGN contributes the most. The AGN hosted in less massive galaxies (i.e., log(Mstar/Msun)<10.5; log(Mdust/Msun)<7.5) show larger molecular gas content with respect to SFGs matched in stellar and dust mass. When comparing their depletion time(tdep~MH2/SFR), AGN show tdep~0.3-1.0 Gyr, similar to those observed in the control sample of SFGs. Seyfert 2 galaxies show fainter PAH luminosity at increasingly larger dominance of the nuclear activity in the MIR. We find no clear evidence for a systematic reduction of the molecular gas reservoir at galactic scale in Seyfert galaxies with respect to SFGs. This is in agreement with recent studies showing that molecular gas content only is reduced in regions of sub-kpc size, where the emission from the accreting supermassive black hole dominates. Nonetheless, we show that the impact of AGN activity on the ISM is clearly visible as suppression of the PAH luminosity.

preprint2020arXiv

Dust and gas content of high-redshift galaxies hosting obscured AGN in the CDF-S

Obscured AGN represent a significant fraction of the entire AGN population, especially at high redshift (~70% at z=3--5). They are often characterized by the presence of large gas and dust reservoirs that are thought to sustain and possibly obscure vigorous star formation processes that make these objects shine at far-IR and sub-mm wavelengths. We exploit ALMA Cycle 4 observations of the continuum (~2.1mm) and high-J CO emission of a sample of six X-ray selected SMGs hosting an obscured AGN at z_spec>2.5 in the 7 Ms CDF-S. We measured the masses and sizes of the dust and molecular gas and we derived the gas density and column density on the basis of a uniform sphere geometry. Finally, we compared the measured column densities with those derived from the Chandra X-ray spectra. We detected both the continuum and line emission for three sources for which we measured both the flux density and size. For the undetected sources, we derived an upper limit on the flux density. We found that the detected galaxies are rich in gas and dust (molecular gas mass in the range <0.5 - 2.7 x 10^10 M_sun for α_CO=0.8 and up to ~2 x 10^11~M_sun for α_CO=6.5, and dust mass <0.9 - 4.9 x 10^8 M_sun) and compact (gas major axis 2.1-3.0 kpc, dust major axis 1.4-2.7 kpc). The column densities associated with the ISM are on the order of 10^(23-24) cm-2, which is comparable with those derived from the X-ray spectra. For the detected sources we also derived dynamical masses in the range 0.8 - 3.7 x 10^10 M_sun. We conclude that the ISM of high redshift galaxies can substantially contribute to nuclear obscuration up to the Compton-thick (>10^24 cm-2) regime. In addition, we found that all the detected sources show a velocity gradient reminding one rotating system, even though two of them show peculiar features in their morphology that can be associated with a chaotic, possibly merging, structure.

preprint2020arXiv

The ALPINE-ALMA [C II] Survey: [C II]158micron Emission Line Luminosity Functions at $z \sim 4-6$

We present the [CII]158$μ$m line luminosity functions (LFs) at $z\sim4-6$ using the ALMA observations of 118 sources, which are selected to have UV luminosity $M_{1500A}<-20.2$ and optical spectroscopic redshifts in COSMOS and ECDF-S. Of the 118 targets, 75 have significant [CII] detections and 43 are upper limits. This is by far the largest sample of [CII] detections which allows us to set constraints to the volume density of [CII] emitters at $z\sim4-6$. But because this is a UV-selected sample, we are missing [CII]-bright but UV-faint sources making our constraints strict lower limits. Our derived LFs are statistically consistent with the $z\sim0$ [CII] LF at $10^{8.25} - 10^{9.75}L_\odot$. We compare our results with the upper limits of the [CII] LF derived from serendipitous sources in the ALPINE maps (Loiacono et al. 2020). We also infer the [CII] LFs based on published far-IR and CO LFs at $z\sim4-6$. Combining our robust lower limits with these additional estimates, we set further constraints to the true number density of [CII] emitters at $z\sim 4 - 6$. These additional LF estimates are largely above our LF at $L_{[CII]}>10^9L_{\odot}$, suggesting that UV-faint but [CII]-bright sources likely make a significant contributions to the [CII] emitter volume density. When we include all the LF estimates, we find that available model predictions underestimate the number densities of [CII] emitters at $z\sim4-6$. Finally, we set a constraint on the molecular gas mass density at $z\sim4-6$, with $ρ_{mol} \sim (2-7)\times10^7M_\odot$\,Mpc$^{-3}$. This is broadly consistent with previous studies.

preprint2020arXiv

The ALPINE-ALMA [CII] Survey: Multi-Wavelength Ancillary Data and Basic Physical Measurements

We present the ancillary data and basic physical measurements for the galaxies in the ALMA Large Program to Investigate C+ at Early Times (ALPINE) survey - the first large multi-wavelength survey which aims at characterizing the gas and dust properties of 118 main-sequence galaxies at redshifts 4.4 < z < 5.9 via the measurement of [CII]-emission at 158 micro-meter (64% at >3.5$σ$) and the surrounding far-infrared (FIR) continuum in conjunction with a wealth of optical and near-infrared data. We outline in detail the spectroscopic data and selection of the galaxies as well as the ground- and space-based imaging products. In addition, we provide several basic measurements including stellar masses, star formation rates (SFR), rest-frame ultra-violet (UV) luminosities, UV continuum slopes ($β$), and absorption line redshifts, as well as H-alpha emission derived from Spitzer colors. We find that the ALPINE sample is representative of the 4 < z < 6 galaxy population selected by photometric methods and only slightly biased towards bluer colors ($Δβ$ ~ 0.2). Using [CII] as tracer of the systemic redshift (confirmed for one galaxy at z=4.5 out of 118 for which we obtained optical [OII]$λ$3727A emission), we confirm red shifted Ly-alpha emission and blue shifted absorption lines similar to findings at lower redshifts. By stacking the rest-frame UV spectra in the [CII] rest-frame we find that the absorption lines in galaxies with high specific SFR are more blue shifted, which could be indicative of stronger winds and outflows.

preprint2019arXiv

The ALPINE-ALMA [CII] survey: Star-formation-driven outflows and circumgalactic enrichment in the early Universe

We study the efficiency of galactic feedback in the early Universe by stacking the [C II] 158 um emission in a large sample of normal star-forming galaxies at 4 < z < 6 from the ALMA Large Program to INvestigate [C II] at Early times (ALPINE) survey. Searching for typical signatures of outflows in the high-velocity tails of the stacked [C II] profile, we observe (i) deviations from a single-component Gaussian model in the combined residuals and (ii) broad emission in the stacked [C II] spectrum, with velocities of |v|<~ 500 km/s. The significance of these features increases when stacking the subset of galaxies with star formation rates (SFRs) higher than the median (SFRmed = 25 Msun/yr), thus confirming their star-formation-driven nature. The estimated mass outflow rates are comparable to the SFRs, yielding mass-loading factors of the order of unity (similarly to local star-forming galaxies), suggesting that star-formation-driven feedback may play a lesser role in quenching galaxies at z > 4. From the stacking analysis of the datacubes, we find that the combined [C II] core emission (|v|< 200 km/s) of the higher-SFR galaxies is extended on physical sizes of ~ 30 kpc (diameter scale), well beyond the analogous [C II] core emission of lower-SFR galaxies and the stacked far-infrared continuum. The detection of such extended metal-enriched gas, likely tracing circumgalactic gas enriched by past outflows, corroborates previous similar studies, confirming that baryon cycle and gas exchanges with the circumgalactic medium are at work in normal star-forming galaxies already at early epochs.

preprint2016arXiv

Characterizing elusive, faint dusty star-forming galaxies: a lensed, optically undetected ALMA galaxy at z~3.3

We present the serendipitous ALMA detection of a faint submillimeter galaxy (SMG) lensed by a foreground z~1 galaxy. By optimizing the source detection to deblend the system, we accurately build the full spectral energy distribution of the distant galaxy from the I814 band to radio wavelengths. It is extremely red, with a I-K colour larger than 2.5. We estimate a photometric redshift of 3.28 and determine the physical parameters. The distant galaxy turns out to be magnified by the foreground lens by a factor of ~1.5, which implies an intrinsic Ks-band magnitude of ~24.5, a submillimeter flux at 870um of ~2.5 mJy and a SFR of ~150-300Msun/yr, depending on the adopted tracer. These values place our source towards the faint end of the distribution of observed SMGs, and in particular among the still few faint SMGs with a fully characterized spectral energy distribution, which allows us not only to accurately estimate its redshift but also to measure its stellar mass and other physical properties. The galaxy studied in this work is a representative of the population of faint SMGs, of which only few objects are known to date, that are undetected in optical and therefore are not typically accounted for when measuring the cosmic star formation history (SFH). This faint galaxy population thus likely represents an important and missing piece in our understanding of the cosmic SFH. Its observation and characterization is of major importance to achieve a solid picture of galaxy evolution.

preprint2016arXiv

Tracing black hole accretion with SED decomposition and IR lines: from local galaxies to the high-z Universe

We present new estimates of AGN accretion and star-formation luminosity in galaxies obtained for the local 12-$μ$m sample of Seyfert galaxies (12MGS), by performing a detailed broad-band spectral energy distribution (SED) decomposition including the emission of stars, dust heated by star formation and a possible AGN dusty torus. Thanks to the availability of data from the X-rays to the sub-millimetre, we constrain and test the contribution of the stellar, AGN and star-formation components to the SEDs. The availability of Spitzer-IRS low resolution mid-infrared (mid-IR) spectra is crucial to constrain the dusty torus component at its peak wavelengths. The results of SED-fitting are also tested against the available information in other bands: the reconstructed AGN bolometric luminosity is compared to those derived from X-rays and from the high excitation IR lines tracing AGN activity like [Ne V] and [O IV]. The IR luminosity due to star-formation (SF) and the intrinsic AGN bolometric luminosity are shown to be strongly related to the IR line luminosity. Variations of these relations with different AGN fractions are investigated, showing that the relation dispersions are mainly due to different AGN relative contribution within the galaxy. Extrapolating these local relations between line and SF or AGN luminosities to higher redshifts, by means of recent Herschel galaxy evolution results, we then obtain mid- and far-IR line luminosity functions useful to estimate how many star-forming galaxies and AGN we expect to detect in the different lines at different redshifts and luminosities with future IR facilities (e.g., JWST, SPICA).

preprint2015arXiv

Are z>2 Herschel galaxies proto-spheroids?

We present a backward approach for the interpretation of the evolution of the near-infrared and the far-infrared luminosity functions across the redshift range 0<z<3. In our method, late-type galaxies are treated by means of a parametric phenomenological method based on PEP/HerMES data up to z~4, whereas spheroids are described by means of a physically motivated backward model. The spectral evolution of spheroids is modelled by means of a single-mass model, associated to a present-day elliptical with K-band luminosity comparable to the break of the local early-type luminosity function. The formation of proto-spheroids is assumed to occurr across the redshift range 1< z < 5. The key parameter is represented by the redshift z_0.5 at which half proto-spheroids are already formed. A statistical study indicates for this parameter values between z_0.5=1.5 and z_0.5=3. We assume as fiducial value z_0.5~2, and show that this assumption allows us to describe accourately the redshift distributions and the source counts. By assuming z_0.5 ~ 2 at the far-IR flux limit of the PEP-COSMOS survey, the PEP-selected sources observed at z>2 can be explained as progenitors of local spheroids caught during their formation. We also test the effects of mass downsizing by dividing the spheroids into three populations of different present-day stellar masses. The results obtained in this case confirm the validity of our approach, i.e. that the bulk of proto-spheroids can be modelled by means of a single model which describes the evolution of galaxies at the break of the present-day early type K-band LF.

preprint2015arXiv

CO luminosity function from Herschel-selected galaxies and the contribution of AGN

We derive the CO luminosity function (LF) for different rotational transitions (i.e. (1-0), (3-2), (5-4)) starting from the Herschel LF by Gruppioni et al. and using appropriate $L'_{\rm CO} - L_{\rm IR}$ conversions for different galaxy classes. Our predicted LFs fit the data so far available at $z\approx0$ and $2$. We compare our results with those obtained by semi-analytical models (SAMs): while we find a good agreement over the whole range of luminosities at $z\approx0$, at $z\approx1$ and $z\approx2$ the tension between our LFs and SAMs in the faint and bright ends increases. We finally discuss the contribution of luminous AGN ($L_{X}>10^{44}\,\rm{erg\,s^{-1}}$) to the bright end of the CO LF concluding that they are too rare to reproduce the actual CO luminosity function at $z\approx2$.

preprint2015arXiv

Mapping the average AGN accretion rate in the SFR-M* plane for Herschel selected galaxies at 0<z<2.5

We study the relation of AGN accretion, star formation rate (SFR), and stellar mass (M$_*$) using a sample of $\approx$ 8600 star-forming galaxies up to z=2.5 selected with \textit{Herschel} imaging in the GOODS and COSMOS fields. For each of them we derive SFR and M$_*$, both corrected, when necessary, for emission from an active galactic nucleus (AGN), through the decomposition of their spectral energy distributions (SEDs). About 10 per cent of the sample are detected individually in \textit{Chandra} observations of the fields. For the rest of the sample we stack the X-ray maps to get average X-ray properties. After subtracting the X-ray luminosity expected from star formation and correcting for nuclear obscuration, we derive the average AGN accretion rate for both detected sources and stacks, as a function of M$_{*}$, SFR and redshift. The average accretion rate correlates with SFR and with M$_*$. The dependence on SFR becomes progressively more significant at z$>$0.8. This may suggest that SFR is the original driver of these correlations. We find that average AGN accretion and star formation increase in a similar fashion with offset from the star-forming "main-sequence". Our interpretation is that accretion onto the central black hole and star formation broadly trace each other, irrespective of whether the galaxy is evolving steadily on the main-sequence or bursting.

preprint2015arXiv

Non-linearity and environmental dependence of the star forming galaxies Main Sequence

Using data from four deep fields (COSMOS, AEGIS, ECDFS, and CDFN), we study the correlation between the position of galaxies in the star formation rate (SFR) versus stellar mass plane and local environment at $z<1.1$. To accurately estimate the galaxy SFR, we use the deepest available Spitzer/MIPS 24 and Herschel/PACS datasets. We distinguish group environments ( $M_{halo}\sim$10$^{12.5-14.2}$$M_{\odot}$) based on the available deep X-ray data and lower halo mass environments based on the local galaxy density. We confirm that the Main Sequence (MS) of star forming galaxies is not a linear relation and there is a flattening towards higher stellar masses ( $M_*>10^{10.4-10.6}$ $M_{\odot}$), across all environments. At high redshift ( $0.5<z<1.1$ ), the MS varies little with environment. At low redshift ( $0.15<z<0.5$ ), group galaxies tend to deviate from the mean MS towards the region of quiescence with respect to isolated galaxies and less-dense environments. We find that the flattening of the MS toward low SFR is due to an increased fraction of bulge dominated galaxies at high masses. Instead, the deviation of group galaxies from the MS at low redshift is caused by a large fraction of red disk dominated galaxies which are not present in the lower density environments. Our results suggest that above a mass threshold ( $\sim10^{10.4}-10^{10.6}$$M_{\odot}$ ) stellar mass, morphology and environment act together in driving the evolution of the SF activity towards lower level. The presence of a dominating bulge and the associated quenching processes are already in place beyond $z\sim$1. The environmental effects appear, instead, at lower redshifts and have a long time-scale.

preprint2015arXiv

Predictions for surveys with the SPICA Mid-infrared Instrument

We present predictions for number counts and redshift distributions of galaxies detectable in continuum and in emission lines with the Mid-infrared (MIR) Instrument (SMI) proposed for the Space Infrared Telescope for Cosmology and Astrophysics (SPICA). We have considered 24 MIR fine-structure lines, four Polycyclic Aromatic Hydrocarbon (PAH) bands (at 6.2, 7.7, 8.6 and 11.3$μ$m) and two silicate bands (in emission and in absorption) at 9.7$μ$m and 18.0$μ$m. Six of these lines are primarily associated with Active Galactic Nuclei (AGNs), the others with star formation. A survey with the SMI spectrometers of 1 hour integration per field-of-view (FoV) over an area of $1\,\hbox{deg}^2$ will yield $5\,σ$ detections of $\simeq 140$ AGN lines and of $\simeq 5.2\times10^{4}$ star-forming galaxies, $\simeq 1.6\times10^{4}$ of which will be detected in at least two lines. The combination of a shallow ($20.0\,\hbox{deg}^{2}$, $1.4\times10^{-1}$ h integration per FoV) and a deep survey ($6.9\times10^{-3}\,\hbox{deg}^{2}$, $635$ h integration time), with the SMI camera, for a total of $\sim$1000 h, will accurately determine the MIR number counts of galaxies and of AGNs over five orders of magnitude in flux density, reaching values more than one order of magnitude fainter than the deepest Spitzer $24\,μ$m surveys. This will allow us to determine the cosmic star formation rate (SFR) function down to SFRs more than 100 times fainter than reached by the Herschel Observatory.

preprint2015arXiv

Relationship between star formation rate and black hole accretion at z=2: The different contributions in quiescent, normal and starburst galaxies

We investigate the co-evolution of black-hole-accretion-rate (BHAR) and star-formation-rate (SFR) in $1.5<z<2.5$ galaxies displaying a greater diversity of star-forming properties compared to previous studies. We combine X-ray stacking and far-IR photometry of stellar mass-limited samples of normal star-forming, starburst and quiescent/quenched galaxies in the COSMOS field. We corroborate the existence of a strong correlation between BHAR (i.e. the X-ray luminosity, L_X), and stellar mass (M*) for normal star-forming galaxies, although find a steeper relation than previously reported. We find that starbursts show a factor of 3 enhancement in BHAR compared to normal SF galaxies (against a factor of 6 excess in SFR), while quiescents show a deficit of a factor 5.5 at a given mass. One possible interpretation of this is that the starburst phase does not coincide with cosmologically relevant BH growth, or that starburst-inducing mergers are more efficient at boosting SFR than BHAR. Contrary to studies based on smaller samples, we find the BHAR/SFR ratio of main sequence (MS) galaxies is not mass invariant, but scales weakly as M*^(0.43\pm0.09}, implying faster BH growth in more massive galaxies at $z\sim2$. Furthermore, BHAR/SFR during the starburst is a factor of 2 lower than in MS galaxies, at odds with the predictions of hydrodynamical simulations of merger galaxies that foresee a sudden enhancement of L_X/SFR during the merger. Finally, we estimate that the bulk of the accretion density of the Universe at $z\sim2$ is associated with normal star-forming systems, with only 6(+/-1)% and 11(+/-1)% associated with starburst and quiescent galaxies, respectively.

preprint2015arXiv

Star Formation in Herschel's Monsters versus Semi-Analytic Models

We present a direct comparison between the observed star formation rate functions (SFRF) and the state-of-the-art predictions of semi-analytic models (SAM) of galaxy formation and evolution. We use the PACS Evolutionary Probe Survey (PEP) and Herschel Multi-tiered Extragalactic Survey (HerMES) data-sets in the COSMOS and GOODS-South fields, combined with broad-band photometry from UV to sub-mm, to obtain total (IR+UV) instantaneous star formation rates (SFRs) for individual Herschel galaxies up to z~4, subtracted of possible active galactic nucleus (AGN) contamination. The comparison with model predictions shows that SAMs broadly reproduce the observed SFRFs up to z~2, when the observational errors on the SFR are taken into account. However, all the models seem to under-predict the bright-end of the SFRF at z>2. The cause of this underprediction could lie in an improper modelling of several model ingredients, like too strong (AGN or stellar) feedback in the brighter objects or too low fall-back of gas, caused by weak feedback and outflows at earlier epochs.

preprint2015arXiv

The most obscured AGN in the COSMOS field

Highly obscured active galactic nuclei (AGN) are common in nearby galaxies, but are difficult to observe beyond the local Universe, where they are expected to significantly contribute to the black hole accretion rate density. Furthermore, Compton-thick (CT) absorbers (NH>10^24 cm^-2) suppress even the hard X-ray (2-10 keV) AGN nuclear emission, and therefore the column density distribution above 10^24 cm^-2 is largely unknown. We present the identification and multi-wavelength properties of a heavily obscured (NH>~10^25 cm^-2), intrinsically luminous (L(2-10keV)>10^44 erg s^-1) AGN at z=0.353 in the COSMOS field. Several independent indicators, such as the shape of the X-ray spectrum, the decomposition of the spectral energy distribution and X-ray/[NeV] and X-ray/6μm luminosity ratios, agree on the fact that the nuclear emission must be suppressed by a 10^25 cm^-2 column density. The host galaxy properties show that this highly obscured AGN is hosted in a massive star-forming galaxy, showing a barred morphology, which is known to correlate with the presence of CT absorbers. Finally, asymmetric and blueshifted components in several optical high-ionization emission lines indicate the presence of a galactic outflow, possibly driven by the intense AGN activity (L(Bol)/L(Edd) = 0.3-0.5). Such highly obscured, highly accreting AGN are intrinsically very rare at low redshift, whereas they are expected to be much more common at the peak of the star formation and BH accretion history, at z~2-3. We demonstrate that a fully multi-wavelength approach can recover a sizable sample of such peculiar sources in large and deep surveys such as COSMOS.

preprint2015arXiv

The star formation rate cookbook at 1 < z < 3: Extinction-corrected relations for UV & [OII]λ3727 luminosities

We use a spectroscopic sample of 286 star-forming galaxies (SFGs) at 1<z<3 from the GMASS survey to study different star formation rate (SFR) estimators. Infrared (IR) data are used to derive empirical calibrations to correct ultraviolet (UV) and [OII]λ3727 luminosities for dust extinction and dust-corrected estimates of SFR. In the selection procedure we fully exploit the available spectroscopic information. On the basis of three continuum indices, we are able to identify and exclude from the sample galaxies in which old stellar populations might bring a non-negligible contribution to IR luminosity (LIR) and continuum reddening. Using Spitzer-MIPS and Herschel-PACS data we derive LIR for two-thirds of our sample. The LIR/LUV ratio is used as a probe of effective attenuation (AIRX) to search for correlations with continuum and spectroscopic features. The relation between AIRX and UV continuum slope (β) was tested for our sample and found to be broadly consistent with the literature results at the same redshift, though with a larger dispersion with respect to UV-selected samples. We find a correlation between the rest-frame equivalent width (EW) of the [OII]λ3727 line and β, which is the main result of this work. We therefore propose the [OII]λ3727 line EW as a dust attenuation probe and calibrate it through AIRX, though the assumption of a reddening curve is still needed to derive the actual attenuation towards the [OII]λ3727 line. We tested the issue of differential attenuation towards stellar continuum and nebular emission: our results are in line with the traditional prescription of extra attenuation towards nebular lines. A set of relations is provided that allows the recovery of the total unattenuated SFR from UV and [OII]λ3727 luminosities. (Abridged)

preprint2014arXiv

A Multiwavelength Consensus on the Main Sequence of Star-Forming Galaxies at z~2

We compare various star formation rate (SFR) indicators for star-forming galaxies at $1.4<z<2.5$ in the COSMOS field. The main focus is on the SFRs from the far-IR (PACS-Herschel data) with those from the ultraviolet, for galaxies selected according to the BzK criterion. FIR-selected samples lead to a vastly different slope of the SFR-stellar mass ($M_*$) relation, compared to that of the dominant main sequence population as measured from the UV, since the FIR selection picks predominantly only a minority of outliers. However, there is overall agreement between the main sequences derived with the two SFR indicators, when stacking on the PACS maps the BzK-selected galaxies. The resulting logarithmic slope of the SFR-{$M_*$} relation is $\sim0.8-0.9$, in agreement with that derived from the dust-corrected UV-luminosity. Exploiting deeper 24$μ$m-Spitzer data we have characterized a sub-sample of galaxies with reddening and SFRs poorly constrained, as they are very faint in the $B$ band. The combination of Herschel with Spitzer data have allowed us to largely break the age/reddening degeneracy for these intriguing sources, by distinguishing whether a galaxy is very red in B-z because of being heavily dust reddened, or whether because star formation has been (or is being) quenched. Finally, we have compared our SFR(UV) to the SFRs derived by stacking the radio data and to those derived from the H$α$ luminosity of a sample of star-forming galaxies at $1.4<z<1.7$. The two sets of SFRs are broadly consistent as they are with the SFRs derived from the UV and by stacking the corresponding PACS data in various mass bins.

preprint2014arXiv

ALMA reveals a warm and compact starburst around a heavily obscured supermassive black hole at z=4.75

We report ALMA Cycle 0 observations at 1.3mm of LESS J033229.4-275619 (XID403), an Ultraluminous Infrared Galaxy at $z=4.75$ in the Chandra Deep Field South hosting a Compton-thick QSO. The source is not resolved in our data at a resolution of $\sim$0.75 arcsec, placing an upper-limit of 2.5 kpc to the half-light radius of the continuum emission from heated-dust. After deconvolving for the beam size, however, we found a $\sim3σ$ indication of an intrinsic source size of $0.27\pm0.08$ arcsec (Gaussian FWHM), which would correspond to $r_{half}\sim0.9\pm0.3$ kpc. We build the far-IR SED of XID403 by combining datapoints from both ALMA and Herschel and fit it with a modified blackbody spectrum. For the first time, we measure the dust temperature $T_d=58.5\pm5.3$ K in this system, which is comparable to what has been observed in other high-z submillimeter galaxies. The measured star formation rate is SFR=$1020\pm150$ $M_{\odot}$ yr$^{-1}$, in agreement with previous estimates at lower S/N. Based on the measured SFR and source size, we constrain the SFR surface density to be $Σ_{SFR}>26\;M_{\odot}$yr$^{-1}$kpc$^{-2}$ ($\sim200\;M_{\odot}$yr$^{-1}$kpc$^{-2}$ for $r_{half}\sim0.9$ kpc). The compactness of this starburst is comparable to what has been observed in other local and high-z starburst galaxies. If the gas mass measured from previous [CII] and CO(2-1) observations at low resolution is confined within the same dust region, assuming $r_{half}\sim0.9\pm0.3$ kpc, this would produce a column density of $N_H\sim0.3-1.1\times10^{24}$cm$^{-2}$ towards the central SMBH, similar to the column density of $\approx1.4\times10^{24}$cm$^{-2}$ measured from the X-rays. Then, in principle, if both gas and dust were confined on sub-kpc scales, this would be sufficient to produce the observed X-ray column density without any need of a pc-scale absorber [abridged].

preprint2014arXiv

Compton Thick AGN in the XMM-COSMOS survey

Heavily obscured, Compton Thick (CT, NH>10^24 cm^-2) AGN may represent an important phase in AGN/galaxy co-evolution and are expected to provide a significant contribution to the cosmic X-ray background (CXB). Through direct X-ray spectra analysis, we selected 39 heavily obscured AGN (NH>3x10^23 cm^-2) in the 2 deg^2 XMM-COSMOS survey. After selecting CT AGN based on the fit of a simple absorbed two power law model to the XMM data, the presence of CT AGN was confirmed in 80% of the sources using deeper Chandra data and more complex models. The final sample of CT AGN comprises 10 sources spanning a large range of redshift and luminosity. We collected the multi-wavelength information available for all these sources, in order to study the distribution of SMBH and host properties, such as BH mass (M_BH), Eddington ratio (λ_Edd), stellar mass (M*), specific star formation rate (sSFR) in comparison with a sample of unobscured AGN. We find that highly obscured sources tend to have significantly smaller M_BH and higher λ_edd with respect to unobscured ones, while a weaker evolution in M* is observed. The sSFR of highly obscured sources is consistent with the one observed in the main sequence of star forming galaxies, at all redshift. We also present optical spectra, spectral energy distribution (SED) and morphology for the sample of 10 CT AGN: all the available optical spectra are dominated by the stellar component of the host galaxy, and a highly obscured torus component is needed in the SED of the CT sources. Exploiting the high resolution Hubble-ACS images available, we conclude that these highly obscured sources have a significantly larger merger fraction with respect to other X-ray selected samples of AGN. Finally we discuss implications in the context of AGN/galaxy co-evolutionary models, and compare our results with the predictions of CXB synthesis models.

preprint2014arXiv

Exploring the relationship between black hole accretion and star formation with blind mid-/far-infrared spectroscopic surveys

We present new estimates of redshift-dependent luminosity functions of IR lines detectable by SPICA/SAFARI and excited both by star formation and by AGN activity. The new estimates improve over previous work by using updated evolutionary models and dealing in a self consistent way with emission of galaxies as a whole, including both the starburst and the AGN component. New relationships between line and AGN bolometric luminosity have been derived and those between line and IR luminosities of the starburst component have been updated. These ingredients were used to work out predictions for the source counts in 11 mid/far-IR emission lines partially or entirely excited by AGN activity. We find that the statistics of the emission line detection of galaxies as a whole is mainly determined by the star formation rate, because of the rarity of bright AGNs. We also find that the slope of the line integral number counts is flatter than 2 implying that the number of detections at fixed observing time increases more by extending the survey area than by going deeper. We thus propose a wide spectroscopic survey of 1 hour integration per field-of-view over an area of 5 deg$^{2}$ to detect (at 5$σ$) $\sim$760 AGNs in [OIV]25.89$μ$m $-$ the brightest AGN mid-infrared line $-$ out to $z\sim2$. Pointed observations of strongly lensed or hyper-luminous galaxies previously detected by large area surveys such as those by Herschel and by the SPT can provide key information on the galaxy-AGN co-evolution out to higher redshifts.

preprint2014arXiv

Herschel Far-IR counterparts of SDSS galaxies: Analysis of commonly used Star Formation Rate estimates

We study a hundred of galaxies from the spectroscopic Sloan Digital Sky Survey with individual detections in the Far-Infrared Herschel PACS bands (100 or 160 $μ$m) and in the GALEX Far-UltraViolet band up to z$\sim$0.4 in the COSMOS and Lockman Hole fields. The galaxies are divided into 4 spectral and 4 morphological types. For the star forming and unclassifiable galaxies we calculate dust extinctions from the UV slope, the H$α$/H$β$ ratio and the $L_{\rm IR}/L_{\rm UV}$ ratio. There is a tight correlation between the dust extinction and both $L_{\rm IR}$ and metallicity. We calculate SFR$_{total}$ and compare it with other SFR estimates (H$α$, UV, SDSS) finding a very good agreement between them with smaller dispersions than typical SFR uncertainties. We study the effect of mass and metallicity, finding that it is only significant at high masses for SFR$_{Hα}$. For the AGN and composite galaxies we find a tight correlation between SFR and L$_{IR}$ ($σ\sim$0.29), while the dispersion in the SFR - L$_{UV}$ relation is larger ($σ\sim$0.57). The galaxies follow the prescriptions of the Fundamental Plane in the M-Z-SFR space.

preprint2014arXiv

Large scale clustering measurements with photometric redshifts: comparing the dark matter halos of X-ray AGN, star-forming and passive galaxies at z~1

We combine multiwavelength data in the AEGIS-XD and C-COSMOS surveys to measure the typical dark matter halo mass of X-ray selected AGN [Lx(2-10keV)>1e42 erg/s] in comparison with far-infrared selected star-forming galaxies detected in the Herschel/PEP survey (PACS Evolutionary Probe; Lir>1e11 solar) and quiescent systems at z~1. We develop a novel method to measure the clustering of extragalactic populations that uses photometric redshift Probability Distribution Functions in addition to any spectroscopy. This is advantageous in that all sources in the sample are used in the clustering analysis, not just the subset with secure spectroscopy. The method works best for large samples. The loss of accuracy because of the lack of spectroscopy is balanced by increasing the number of sources used to measure the clustering. We find that X-ray AGN, far-infrared selected star-forming galaxies and passive systems in the redshift interval 0.6<z<1.4 are found in halos of similar mass, $\log M_{DMH}/(M_{\odot}\,h^{-1})\approx13.0$. We argue that this is because the galaxies in all three samples (AGN, star-forming, passive) have similar stellar mass distributions, approximated by the J-band luminosity. Therefore all galaxies that can potentially host X-ray AGN, because they have stellar masses in the appropriate range, live in dark matter haloes of $\log M_{DMH}/(M_{\odot}\,h^{-1})\approx13.0$ independent of their star-formation rates. This suggests that the stellar mass of X-ray AGN hosts is driving the observed clustering properties of this population. We also speculate that trends between AGN properties (e.g. luminosity, level of obscuration) and large scale environment may be related to differences in the stellar mass of the host galaxies.

preprint2014arXiv

The evolution of galaxy star formation activity in massive halos

There is now a large consensus that the current epoch of the Cosmic Star Formation History (CSFH) is dominated by low mass galaxies while the most active phase at 1<z<2 is dominated by more massive galaxies, which undergo a faster evolution. Massive galaxies tend to inhabit very massive halos such as galaxy groups and clusters. We aim to understand whether the observed "galaxy downsizing" could be interpreted as a "halo downsizing", whereas the most massive halos, and their galaxy populations, evolve more rapidly than the halos of lower mass. Thus, we study the contribution to the CSFH of galaxies inhabiting group-sized halos. This is done through the study of the evolution of the Infra-Red (IR) luminosity function of group galaxies from redshift 0 to ~1.6. We use a sample of 39 X-ray selected groups in the Extended Chandra Deep Field South (ECDFS), the Chandra Deep Field North (CDFN), and the COSMOS field, where the deepest available mid- and far-IR surveys have been conducted with Spitzer MIPS and Hersche PACS. Groups at low redshift lack the brightest, rarest, and most star forming IR-emitting galaxies observed in the field. Their IR-emitting galaxies contribute <10% of the comoving volume density of the whole IR galaxy population in the local Universe. At redshift >~1, the most IR-luminous galaxies (LIRGs and ULIRGs) are preferentially located in groups, and this is consistent with a reversal of the star-formation rate vs .density anti-correlation observed in the nearby Universe. At these redshifts, group galaxies contribute 60-80% of the CSFH, i.e. much more than at lower redshifts. Below z~1, the comoving number and SFR densities of IR-emitting galaxies in groups decline significantly faster than those of all IR-emitting galaxies. Our results are consistent with a "halo downsizing" scenario and highlight the significant role of "environment" quenching in shaping the CSFH.

preprint2014arXiv

The evolution of star formation activity in galaxy groups

We study the evolution of the total star formation (SF) activity, total stellar mass and halo occupation distribution in massive halos by using one of the largest X-ray selected sample of galaxy groups with secure spectroscopic identification in the major blank field surveys (ECDFS, CDFN, COSMOS, AEGIS). We provide an accurate measurement of SFR for the bulk of the star-forming galaxies using very deep mid-infrared Spitzer MIPS and far-infrared Herschel PACS observations. For undetected IR sources, we provide a well-calibrated SFR from SED fitting. We observe a clear evolution in the level of SF activity in galaxy groups. The total SF activity in the high redshift groups (0.5<z<1.1) is higher with respect to the low redshift (0.15<z<0.5) sample at any mass by 0.8+/-0.12 dex. A milder difference (0.35+/-0.1 dex) is observed between the low redshift bin and the groups at z~0. We show that the level of SF activity is declining more rapidly in the more massive halos than in the more common lower mass halos. We do not observe any evolution in the halo occupation distribution and total stellar mass- halo mass relations in groups. The picture emerging from our findings suggests that the galaxy population in the most massive systems is evolving faster than galaxies in lower mass halos, consistently with a "halo downsizing" scenario.

preprint2014arXiv

The far-infrared/radio correlation and radio spectral index of galaxies in the SFR-M* plane up to z 2

[Abridged] We study the evolution of the radio spectral index and far-infrared/radio correlation (FRC) across the star-formation rate-stellar masse (i.e. SFR-M*) plane up to z 2. We start from a M*-selected sample of galaxies with reliable SFR and redshift estimates. We then grid the SFR-M* plane in several redshift ranges and measure the infrared luminosity, radio luminosity, radio spectral index, and ultimately the FRC index (i.e. qFIR) of each SFR-M*-z bin. The infrared luminosities of our SFR-M*-z bins are estimated using their stacked far-infrared flux densities inferred from observations obtained with Herschel. Their radio luminosities and radio spectral indices (i.e. alpha, where Snu nu^-alpha) are estimated using their stacked 1.4GHz and 610MHz flux densities from the VLA and GMRT, respectively. Our far-infrared and radio observations include the most widely studied blank extragalactic fields -GOODS-N/S, ECDFS, and COSMOS- covering a sky area of 2deg^2. Using this methodology, we constrain the radio spectral index and FRC index of star-forming galaxies with M*>10^10Msun and 0<z<2.3. We find that alpha^1.4GHz_610MHz does not evolve significantly with redshift or with the distance of a galaxy with respect to the main sequence (MS) of the SFR-M* plane (i.e. Delta_log(SSFR)_MS=log[SSFR(galaxy)/SSFR_MS(M*,z)]). Instead, star-forming galaxies have a radio spectral index consistent with a canonical value of 0.8, which suggests that their radio spectra are dominated by non-thermal optically thin synchrotron emission. We find that qFIR displays a moderate but statistically significant redshift evolution as qFIR(z)=(2.35+/-0.08)*(1+z)^(-0.12+/-0.04), consistent with some previous literature. Finally, we find no significant correlation between qFIR and Delta_log(SSFR)_MS, though a weak positive trend, as observed in one of our redshift bins, cannot be firmly ruled out using our dataset.

preprint2014arXiv

The PEP Survey: Infrared Properties of Radio-Selected AGN

By exploiting the VLA-COSMOS and the Herschel-PEP surveys, we investigate the Far Infrared (FIR) properties of radio-selected AGN. To this purpose, from VLA-COSMOS we considered the 1537, F[1.4 GHz]>0.06 mJy sources with a reliable redshift estimate, and sub-divided them into star-forming galaxies and AGN solely on the basis of their radio luminosity. The AGN sample is complete with respect to radio selection at all z<~3.5. 832 radio sources have a counterpart in the PEP catalogue. 175 are AGN. Their redshift distribution closely resembles that of the total radio-selected AGN population, and exhibits two marked peaks at z~0.9 and z~2.5. We find that the probability for a radio-selected AGN to be detected at FIR wavelengths is both a function of radio power and redshift, whereby powerful sources are more likely to be FIR emitters at earlier epochs. This is due to two distinct effects: 1) at all radio luminosities, FIR activity monotonically increases with look-back time and 2) radio activity of AGN origin is increasingly less effective at inhibiting FIR emission. Radio-selected AGN with FIR emission are preferentially located in galaxies which are smaller than those hosting FIR-inactive sources. Furthermore, at all z<~2, there seems to be a preferential (stellar) mass scale M ~[10^{10}-10^{11}] Msun which maximizes the chances for FIR emission. We find such FIR (and MIR) emission to be due to processes indistinguishable from those which power star-forming galaxies. It follows that radio emission in at least 35% of the entire AGN population is the sum of two contributions: AGN accretion and star-forming processes within the host galaxy.

preprint2014arXiv

The role of massive halos in the Star Formation History of the Universe

The most striking feature of the Cosmic Star Formation History (CSFH) of the Universe is a dramatic drop of the star formation (SF) activity, since z~1. In this work we investigate if the very same process of assembly and growth of structures is one of the major drivers of the observed decline. We study the contribution to the CSFH of galaxies in halos of different masses. This is done by studying the total SFR-halo mass-redshift plane from redshift 0 to redshift z~1.6 in a sample of 57 groups and clusters by using the deepest available mid- and far-infrared surveys conducted with Spitzer MIPS and Herschel PACS and SPIRE. Our results show that low mass groups provide a 60-80% contribution to the CSFH at z~1. Such contribution declines faster than the CSFH in the last 8 billion years to less than 10% at z<0.3, where the overall SF activity is sustained by lower mass halos. More massive systems provide only a marginal contribution (<10%) at any epoch. A simplified abundance matching method shows that the large contribution of low mass groups at z~1 is due to a large fraction (>50%) of very massive, highly star forming Main Sequence galaxies. Below z~1 a quenching process must take place in massive halos to cause the observed faster suppression of their SF activity. Such process must be a slow one though, as most of the models implementing a rapid quenching of the SF activity in accreting satellites significantly underpredicts the observed SF level in massive halos at any redshift. Starvation or the transition from cold to hot accretion would provide a quenching timescale of 1 Gyrs more consistent with the observations. Our results suggest a scenario in which, due to the structure formation process, more and more galaxies experience the group environment and, thus, the associated quenching process. This leads to the progressive suppression of their SF activity shaping the CSFH below z~1.

preprint2014arXiv

Tracing the cosmic growth of super massive black holes to z ~ 3 with Herschel

We study a sample of Herschel-PACS selected galaxies within the GOODS-South and the COSMOS fields in the framework of the PACS Evolutionary Probe (PEP) project. Starting from the rich multi-wavelength photometric data-sets available in both fields, we perform a broad-band Spectral Energy Distribution (SED) decomposition to disentangle the possible active galactic nucleus (AGN) contribution from that related to the host galaxy. We find that 37 per cent of the Herschel-selected sample shows signatures of nuclear activity at the 99 per cent confidence level. The probability to reveal AGN activity increases for bright ($L_{\rm 1-1000} > 10^{11} \rm L_{\odot}$) star-forming galaxies at $z>0.3$, becoming about 80 per cent for the brightest ($L_{\rm 1-1000} > 10^{12} \rm L_{\odot}$) infrared (IR) galaxies at $z \geq 1$. Finally, we reconstruct the AGN bolometric luminosity function and the super-massive black hole growth rate across cosmic time up to $z \sim 3$ from a Far-Infrared (FIR) perspective. This work shows general agreement with most of the panchromatic estimates from the literature, with the global black hole growth peaking at $z \sim 2$ and reproducing the observed local black hole mass density with consistent values of the radiative efficiency $ε_{\rm rad}$ ($\sim$0.07).

preprint2013arXiv

Cosmological model dependence of the galaxy luminosity function: far-infrared results in the Lemaitre-Tolman-Bondi model

This is the first paper of a series aiming at investigating galaxy formation and evolution in the giant-void class of the Lemaitre-Tolman-Bondi (LTB) models that best fits current cosmological observations. Here we investigate the Luminosity Function (LF) methodology, and how its estimates would be affected by a change on the cosmological model assumed in its computation. Are the current observational constraints on the allowed Cosmology enough to yield robust LF results? We use the far-infrared source catalogues built on the observations performed with the Herschel/PACS instrument, and selected as part of the PACS evolutionary probe (PEP) survey. Schechter profiles are obtained in redshift bins up to z approximately 4, assuming comoving volumes in both the standard model, that is, Friedmann-Lemaitre-Robertson-Walker metric with a perfect fluid energy-momentum tensor, and non-homogeneous LTB dust models, parametrized to fit the current combination of results stemming from the observations of supernovae Ia, the cosmic microwave background, and baryonic acoustic oscillations. We find that the luminosity functions computed assuming both the standard model and LTB void models show in general good agreement. However, the faint-end slope in the void models shows a significant departure from the standard model up to redshift 0.4. We demonstrate that this result is not artificially caused by the used LF estimator which turns out to be robust under the differences in matter-energy density profiles of the models. The differences found in the LF slopes at the faint end are due to variation in the luminosities of the sources, which depend on the geometrical part of the model. It follows that either the standard model is over-estimating the number density of faint sources or the void models are under-estimating it.

preprint2013arXiv

Herschel reveals the obscured star formation in HiZELS Hα emitters at z=1.47

We describe the far-infrared (FIR; rest-frame 8--1000μm) properties of a sample of 443 Hα-selected star-forming galaxies in the COSMOS and UDS fields detected by the HiZELS imaging survey. Sources are identified using narrow-band filters in combination with broad-band photometry to uniformly select Hα (and [OII] if available) emitters in a narrow redshift slice at z = 1.47+/-0.02. We use a stacking approach in Spitzer, Herschel (from PEP and HerMES surveys) and AzTEC images to describe their typical FIR properties. We find that HiZELS galaxies with observed Hα luminosities of ~ 10^{8.1-9.1} Lo have bolometric FIR luminosities of typical LIRGs, L_FIR ~ 10^{11.48+/-0.05} Lo. Combining the Hα and FIR luminosities, we derive median SFR = 32+/-5 Mo/yr and Hα extinctions of A(Hα) = 1.0+/-0.2 mag. Perhaps surprisingly, little difference is seen in typical HiZELS extinction levels compared to local star-forming galaxies. We confirm previous empirical stellar mass (M*) to A(Hα) relations and the little or no evolution up to z = 1.47. For HiZELS galaxies, we provide an empirical parametrisation of the SFR as a function of (u-z)_rest colours and 3.6μm photometry. We find that the observed Hα luminosity is a dominant SFR tracer when (u-z)_rest ~< 0.9 mag or when 3.6μm photometry > 22 mag (Vega) or when M* < 10^9.7 Mo. We do not find any correlation between the [OII]/Hα and FIR luminosity, suggesting that this emission line ratio does not trace the extinction of the most obscured star-forming regions. The luminosity-limited HiZELS sample tends to lie above of the so-called `main sequence' for star-forming galaxies, especially at low M*. This work suggests that obscured star formation is linked to the assembly of M*, with deeper potential wells in massive galaxies providing dense, heavily obscured environments in which stars can form rapidly.

preprint2013arXiv

Lyman break and UV-selected galaxies at z ~ 1: II. PACS-100um/160um FIR detections

We report the PACS-100um/160um detections of a sample of 42 GALEX-selected and FIR-detected Lyman break galaxies (LBGs) at z ~ 1 located in the COSMOS field and analyze their ultra-violet (UV) to far-infrared (FIR) properties. The detection of these LBGs in the FIR indicates that they have a dust content high enough so that its emission can be directly detected. According to a spectral energy distribution (SED) fitting with stellar population templates to their UV-to-near-IR observed photometry, PACS-detected LBGs tend to be bigger, more massive, dustier, redder in the UV continuum, and UV-brighter than PACS-undetected LBGs. PACS-detected LBGs at z ~ 1 are mostly disk-like galaxies and are located over the green-valley and red sequence of the color-magnitude diagram of galaxies at their redshift. By using their UV and IR emission, we find that PACS-detected LBGs tend to be less dusty and have slightly higher total star-formation rates (SFRs) than other PACS-detected UV-selected galaxies within their same redshift range. As a consequence of the selection effect due to the depth of the FIR observations employed, all our PACS-detected LBGs are LIRGs. However, none of them are in the ULIRG regime, where the FIR observations are complete. The finding of ULIRGs-LBGs at higher redshifts suggests an evolution of the FIR emission of LBGs with cosmic time. In an IRX-$β$ diagram, PACS-detected LBGs at z ~ 1 tend to be located around the relation for local starburst similarly to other UV-selected PACS-detected galaxies at their same redshift. Consequently, the dust-correction factors obtained with their UV continuum slope allow to determine their total SFR, unlike at higher redshifts. However, the dust attenuation derived from UV to NIR SED fitting overestimates the total SFR for most of our PACS-detected LBGs in age-dependent way: the overestimation factor is higher in younger galaxies.

preprint2013arXiv

On the redshift evolution (0 < z < 4) of dust attenuation and of the total (UV+IR) star formation rate density

Using new homogeneous LFs in the FUV and in the FIR Herschel/PEP and Herschel/HerMES, we study the evolution of the dust attenuation with redshift. With this information in hand, we are able to estimate the redshift evolution of the total (FUV + FIR) star formation rate density SFRD_TOT. By integrating SFRD_TOT, we follow the mass building and analyze the redshift evolution of the stellar mass density (SMD). This letter aims at providing a complete view of star formation from the local universe to z = 4 and, using assumptions on earlier star formation history, compares this evolution to what was known before in an attempt to draw a homogeneous picture of the global evolution of star formation in galaxies. The main conclusions of this letter are: 1) the dust attenuation A_FUV is found to increase from z = 0 to z \sim 1.2 and then starts to decrease up to our last data point at z = 3.6; 2) the estimated SFRD confirms published results up to z = 2. At z > 2, we observe either a plateau or a small increase up to z = 3 and then a likely decrease up to z = 3.6; 3) the peak of A_FUV is delayed with respect to the plateau of SFRD_TOT and a likely origin might be found in the evolution of the bright ends of the FUV and FIR LFs; 4) using assumptions (namely exponential rise and linear rise with time) for the evolution of the star formation density from z = 3.6 to z_form = 10, we integrate SFRD_TOT and find a good agreement with the published SMDs.

preprint2013arXiv

Panchromatic spectral energy distributions of Herschel sources

(abridged) Far-infrared Herschel photometry from the PEP and HerMES programs is combined with ancillary datasets in the GOODS-N, GOODS-S, and COSMOS fields. Based on this rich dataset, we reproduce the restframe UV to FIR ten-colors distribution of galaxies using a superposition of multi-variate Gaussian modes. The median SED of each mode is 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 AGN. The defined Gaussian grouping is also used to identify rare 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 eight other 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 micron detected sources in PEP GOODS fields. 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.

preprint2013arXiv

Reversal or no reversal: the evolution of the star formation rate-density relation up to z~1.6

We investigate the evolution of the star formation rate (SFR)-density relation in the Extended Chandra Deep Field South (ECDFS) and the Great Observatories Origin Deep Survey (GOODS) fields up to z~1.6. In addition to the "traditional method", in which the environment is defined according to a statistical measurement of the local galaxy density, we use a "dynamical" approach, where galaxies are classified according to three different environment regimes: group, "filament-like", and field. Both methods show no evidence of a SFR-density reversal. Moreover, group galaxies show a mean SFR lower than other environments up to z~1, while at earlier epochs group and field galaxies exhibit consistent levels of star formation (SF) activity. We find that processes related to a massive dark matter halo must be dominant in the suppression of the SF below z~1, with respect to purely density-related processes. We confirm this finding by studying the distribution of galaxies in different environments with respect to the so-called Main Sequence (MS) of star-forming galaxies. Galaxies in both group and "filament-like" environments preferentially lie below the MS up to z~1, with group galaxies exhibiting lower levels of star-forming activity at a given mass. At z>1, the star-forming galaxies in groups reside on the MS. Groups exhibit the highest fraction of quiescent galaxies up to z~1, after which group, "filament-like", and field environments have a similar mix of galaxy types. We conclude that groups are the most efficient locus for star-formation quenching. Thus, a fundamental difference exists between bound and unbound objects, or between dark matter haloes of different masses.

preprint2013arXiv

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-PACS (Photodetector Array Camera and Spectrometer) far-infrared blank field extragalactic survey, obtained by combining observations of the GOODS (Great Observatories Origins Deep Survey) fields from the PACS Evolutionary Probe (PEP) and GOODS-Herschel key programmes. We describe data reduction and the construction of images and catalogues. In the deepest parts of the GOODS-S field, the catalogues reach 3-sigma depths of 0.9, 0.6 and 1.3 mJy at 70, 100 and 160 um, respectively, and resolve ~75% of the cosmic infrared background at 100um and 160um 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 LIR=10^11 Lsun at z~1 and LIR=10^12 Lsun 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-sigma) 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.

preprint2013arXiv

The dust content of QSO hosts at high redshift

Infrared observations of high-z quasar (QSO) hosts indicate the presence of large masses of dust in the early universe. When combined with other observables, such as neutral gas masses and star formation rates, the dust content of z~6 QSO hosts may help constraining their star formation history. We have collected a database of 58 sources from the literature discovered by various surveys and observed in the FIR. We have interpreted the available data by means of chemical evolution models for forming proto-spheroids, investigating the role of the major parameters regulating star formation and dust production. For a few systems, given the derived small dynamical masses, the observed dust content can be explained only assuming a top-heavy initial mass function, an enhanced star formation efficiency and an increased rate of dust accretion. However, the possibility that, for some systems, the dynamical mass has been underestimated cannot be excluded. If this were the case, the dust mass can be accounted for by standard model assumptions. We provide predictions regarding the abundance of the descendants of QSO hosts; albeit rare, such systems should be present and detectable by future deep surveys such as Euclid already at z>4.

preprint2013arXiv

The evolution of the dust and gas content in galaxies

We use deep Herschel PACS and SPIRE observations in GOODSS, GOODSN and COSMOS to estimate the average dust mass (Mdust) of galaxies on a redshift-stellar mass (Mstar)-SFR grid. We study the scaling relations between Mdust, Mstar and SFR at z<=2.5. No clear evolution of Mdust is observed at fixed SFR and Mstar. We find a tight correlation between SFR and Mdust, likely a consequence of the Schmidt-Kennicutt (S-K) law. The Mstar-Mdust correlation observed by previous works flattens or sometimes disappears when fixing the SFR. Most of it likely derives from the combination of the Mdust-SFR and Mstar-SFR correlations. We then investigate the gas content as inferred by converting Mdust by assuming that the dust/gas ratio scales linearly with the gas metallicity. All galaxies in the sample follow, within uncertainties, the same SFR-Mgas relation (integrated S-K law), which broadly agrees with CO-based results for the bulk of the population, despite the completely different approaches. The majority of galaxies at z~2 form stars with an efficiency (SFE=SFR/Mgas) ~5 times higher than at z~0. It is not clear what fraction of such variation is an intrinsic redshift evolution and what fraction arises from selection effects. The gas fraction (fgas) decreases with Mstar and increases with SFR, and does not evolve with z at fixed Mstar and SFR. We explain these trends by introducing a universal relation between fgas, Mstar and SFR, non-evolving out to z~2.5. Galaxies move across this relation as their gas content evolves in time. We use the 3D fundamental fgas-Mstar-SFR relation and the redshift evolution of the Main Sequence to estimate the evolution of fgas in the average population of galaxies as a function of z and Mstar, and we find evidence a downsizing scenario.

preprint2013arXiv

The evolution of the dust temperatures of galaxies in the SFR$-M_{\ast}$ plane up to $z$$\,\thicksim\,$$2$

[Abridged] We study the evolution of the dust temperatures of galaxies in the SFR-M* plane up to z~2 using observations from the Herschel Space Observatory. 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 Tdust of each SFR-M*-z bin. Dust temperatures are inferred using the stacked far-infrared flux densities of our SFR-M*-z bins. At all redshifts, Tdust increases with infrared luminosities (LIR), specific SFRs (SSFR; i.e., SFR/M*) and distances with respect to the main sequence (MS) of the SFR-M* plane (i.e., D_SSFR_MS=log[SSFR(galaxy)/SSFR_MS(M*,z)]). The Tdust-SSFR and Tdust-D_SSFR_MS correlations are statistically more significant than the Tdust-LIR one. While the slopes of these three correlations are redshift-independent, their normalizations evolve from z=0 and z~2. We convert these results into a recipe to derive Tdust from SFR, M* and z. The existence of a strong Tdust-D_SSFR_MS correlation provides us with information on the dust and gas content of galaxies. (i) The slope of the Tdust-D__SSFR_MS correlation can be explained by the increase of the star-formation efficiency (SFE; SFR/Mgas) with D_SSFR_MS as found locally by molecular gas studies. (ii) At fixed D_SSFR_MS, the constant Tdust observed in galaxies probing large ranges in SFR and M* can be explained by an increase or decrease of the number of star-forming regions with comparable SFE enclosed in them. (iii) At high redshift, the normalization towards hotter temperature of the Tdust-D_SSFR_MS correlation can be explained by the decrease of the metallicities of galaxies or by the increase of 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.

preprint2013arXiv

The Herschel census of infrared SEDs through cosmic time

Using Herschel data from the deepest SPIRE and PACS surveys (HerMES and PEP) in COSMOS and GOODS (N+S), we examine the dust properties of IR-luminous (L_IR>10^10 L_sun) galaxies at 0.1<z<2 and determine how these evolve with cosmic time. The unique angle of this work is the rigorous analysis of survey selection effects, making this the first study of the star-formation-dominated, IR-luminous population within a framework almost entirely free of selection biases. We find that IR-luminous galaxies have SEDs with broad far-IR peaks characterised by cool/extended dust emission and average dust temperatures in the 25-45K range. Hot (T>45K) SEDs and cold (T<25K), cirrus-dominated SEDs are rare, with most sources being within the range occupied by warm starbursts such as M82 and cool spirals such as M51. We observe a luminosity-temperature (L-T) relation, where the average dust temperature of log[L_IR/L_sun]=12.5 galaxies is about 10K higher than that of their log[L_IR/L_sun]=10.5 counterparts. However, although the increased dust heating in more luminous systems is the driving factor behind the L-T relation, the increase in dust mass and/or starburst size with luminosity plays a dominant role in shaping it. Our results show that the dust conditions in IR-luminous sources evolve with cosmic time: at high redshift, dust temperatures are on average up to 10K lower than what is measured locally. This is manifested as a flattening of the L-T relation, suggesting that (U)LIRGs in the early Universe are typically characterised by a more extended dust distribution and/or higher dust masses than local equivalent sources. Interestingly, the evolution in dust temperature is luminosity dependent, with the fraction of LIRGs with T<35K showing a 2-fold increase from z~0 to z~2, whereas that of ULIRGs with T<35K shows a 6-fold increase.

preprint2013arXiv

The Herschel PEP/HerMES Luminosity Function. I: Probing the Evolution of PACS selected Galaxies to z~4

We exploit the deep and extended far infrared data sets (at 70, 100 and 160 um) of the Herschel GTO PACS Evolutionary Probe (PEP) Survey, in combination with the HERschel Multi tiered Extragalactic Survey (HerMES) data at 250, 350 and 500 um, to derive the evolution of the restframe 35 um, 60 um, 90 um, and total infrared (IR) luminosity functions (LFs) up to z~4. We detect very strong luminosity evolution for the total IR LF combined with a density evolution. In agreement with previous findings, the IR luminosity density increases steeply to z~1, then flattens between z~1 and z~3 to decrease at z greater than 3. Galaxies with different SEDs, masses and sSFRs evolve in very different ways and this large and deep statistical sample is the first one allowing us to separately study the different evolutionary behaviours of the individual IR populations contributing to the IR luminosity density. Galaxies occupying the well established SFR/stellar mass main sequence (MS) are found to dominate both the total IR LF and luminosity density at all redshifts, with the contribution from off MS sources (0.6 dex above MS) being nearly constant (~20% of the total IR luminosity density) and showing no significant signs of increase with increasing z over the whole 0.8<z<2.2 range. Sources with mass in the 10< log(M/Msun) <11 range are found to dominate the total IR LF, with more massive galaxies prevailing at the bright end of the high-z LF. A two-fold evolutionary scheme for IR galaxies is envisaged: on the one hand, a starburst-dominated phase in which the SMBH grows and is obscured by dust, is followed by an AGN dominated phase, then evolving toward a local elliptical. On the other hand, moderately starforming galaxies containing a low-luminosity AGN have various properties suggesting they are good candidates for systems in a transition phase preceding the formation of steady spiral galaxies.

preprint2013arXiv

The Herschel-PEP survey: evidence for downsizing in the hosts of dusty star-forming systems

By making use of Herschel-PEP observations of the COSMOS and Extended Groth Strip fields, we have estimated the dependence of the clustering properties of FIR-selected sources on their 100um fluxes. Our analysis shows a tendency for the clustering strength to decrease with limiting fluxes: r0(S100um >8 mJy)~4.3 Mpc and r0(S100um >5 mJy)~5.8 Mpc. These values convert into minimum halo masses Mmin~10^{11.6} Msun for sources brighter than 8 mJy and Mmin~10^{12.4} Msun for S100um > 5 mJy galaxies. We show such an increase of the clustering strength to be due to an intervening population of z~2 sources, which are very strongly clustered and whose relative contribution, equal to about 10% of the total counts at S100um > 2 mJy, rapidly decreases for brighter flux cuts. By removing such a contribution, we find that z <~ 1 FIR galaxies have approximately the same clustering properties, irrespective of their flux level. The above results were then used to investigate the intrinsic dependence on cosmic epoch of the clustering strength of dusty star-forming galaxies between z~0 and z~2.5. In order to remove any bias in the selection process, the adopted sample only includes galaxies observed at the same rest-frame wavelength, lambda~60 um, which have comparable luminosities and therefore star-formation rates (SFR>~100 Msun/yr). Our analysis shows that the same amount of (intense) star forming activity takes place in extremely different environments at the different cosmological epochs. For z<~1 the hosts of such star forming systems are small, Mmin~10^{11} Msun, isolated galaxies. High (z~2) redshift star formation instead seems to uniquely take place in extremely massive/cluster-like halos, Mmin~10^{13.5} Msun, which are associated with the highest peaks of the density fluctuation field at those epochs. (abridged)

preprint2013arXiv

The incidence of obscuration in active galactic nuclei

We study the incidence of nuclear obscuration on a complete sample of 1310 AGN selected on the basis of their rest-frame 2-10 keV X-ray flux from the XMM-COSMOS survey, in the redshift range 0.3<z<3.5. We classify the AGN as obscured or un-obscured on the basis of either the optical spectral properties and the overall SED or the shape of the X-ray spectrum. The two classifications agree in about 70% of the objects, and the remaining 30% can be further subdivided into two distinct classes: at low luminosities X-ray un-obscured AGN do not always show signs of broad lines or blue/UV continuum emission in their optical spectra, most likely due to galaxy dilution effects; at high luminosities broad line AGN may have absorbed X-ray spectra, which hints at an increased incidence of small-scale (sub-parsec) dust-free obscuration. We confirm that the fraction of obscured AGN is a decreasing function of the intrinsic X-ray luminosity, while the incidence of absorption shows significant evolution only for the most luminous AGN, which appear to be more commonly obscured at higher redshift. We find no significant difference between the mean stellar masses and star formation rates of obscured and un-obscured AGN hosts. We conclude that the physical state of the medium responsible for obscuration in AGN is complex, and mainly determined by the radiation environment (nuclear luminosity) in a small region enclosed within the gravitational sphere of influence of the central black hole, but is largely insensitive to the wider scale galactic conditions.

preprint2013arXiv

The Mean Star Formation Rate of X-ray selected Active Galaxies and its Evolution from z=2.5: Results from PEP-Herschel

(Abridged) We study relationships between the SFR and the nuclear properties of X-ray selected AGNs out to z=2.5, using far-IR data in three extragalactic deep fields as part of the PACS Evolutionary Probe (PEP) program. Guided by studies of intrinsic infra-red AGN SEDs, we show that the majority of the FIR emission in AGNs is produced by cold dust heated by star-formation. We uncover characteristic redshift-dependent trends between the mean FIR luminosity (L_fir) and accretion luminosity (L_agn) of AGNs. 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 L_fir and L_agn, but only among AGNs at low and moderate redshifts (z<1). We interpret this as a signature of the increasing importance of major-mergers in driving both the growth of super-massive black holes (SMBHs) and global star-formation in their hosts at high AGN luminosities. However, we also find that the enhancement of SFR in luminous AGNs weakens or disappears at high redshifts (z>1). This suggests that the role of mergers in SMBH-galaxy co-evolution is less important at these epochs. At all redshifts, we find essentially no relationship between L_fir and nuclear obscuration across five orders of magnitude in obscuring column density, suggesting that various different mechanisms are likely to be responsible for obscuring X-rays in active galaxies. We explain our results within a scenario in which two different modes of SMBH fueling operate among low- and high-luminosity AGNs. We postulate, guided by emerging knowledge about the properties of high redshift galaxies, that the dominant mode of accretion among high-luminosity AGNs evolves with redshift.

preprint2012arXiv

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

We study a sample of 61 submillimetre galaxies (SMGs) selected from ground-based surveys, with known spectroscopic redshifts and observed with Herschel as part of the PACS Evolutionary Probe (PEP) and the Herschel Multi-tiered Extragalactic Survey (HerMES) key programmes. We use the broad far-infrared wavelength coverage (100-600um) provided by the combination of PACS and SPIRE observations. Using a power-law temperature distribution model to derive infrared luminosities and dust temperatures, we measure a dust emissivity spectral index for SMGs of beta=2.0+/-0.2. Our results unveil the diversity of the SMG population. Some SMGs exhibit extreme infrared luminosities of ~10^13 Lsun and relatively warm dust components, while others are fainter (~10^12 Lsun) and are biased towards cold dust temperatures. The extreme infrared luminosities of some SMGs (LIR>10^12.7 Lsun, 26/61 systems) imply SFRs of >500Msun yr^-1. 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 (~40K as opposed to ~25K), consistent with observations of local ULIRGs triggered by major mergers and with results from hydrodynamic simulations of major mergers combined with radiative transfer calculations. Luminous SMGs are also 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 thus might be evolving through a secular mode of star formation. [abridged]

preprint2012arXiv

A Population of z>2 Far-Infrared Herschel-SPIRE selected Starbursts

We present spectroscopic observations for a sample of 36 Herschel-SPIRE 250-500um selected galaxies (HSGs) at 2<z<5 from the Herschel Multi-tiered Extragalactic Survey (HerMES). Redshifts are confirmed as part of a large redshift survey of Herschel-SPIRE-selected sources covering ~0.93deg^2 in six extragalactic legacy fields. Observations were taken with the Keck I Low Resolution Imaging Spectrometer (LRIS) and the Keck II DEep Imaging Multi-Object Spectrograph (DEIMOS). Precise astrometry, needed for spectroscopic follow-up, is determined by identification of counterparts at 24um or 1.4GHz using a cross-identification likelihood matching method. Individual source luminosities range from log(L_IR/Lsun)=12.5-13.6 (corresponding to star formation rates 500-9000Msun/yr, assuming a Salpeter IMF), constituting some of the most intrinsically luminous, distant infrared galaxies yet discovered. We present both individual and composite rest-frame ultraviolet spectra and infrared spectral energy distributions (SEDs). The selection of these HSGs is reproducible and well characterized across large areas of sky in contrast to most z>2 HyLIRGs in the literature which are detected serendipitously or via tailored surveys searching only for high-z HyLIRGs; therefore, we can place lower limits on the contribution of HSGs to the cosmic star formation rate density at (7+-2)x10^(-3)Msun/yr h^3Mpc^(-3) at z~2.5, which is >10% of the estimated total star formation rate density (SFRD) of the Universe from optical surveys. The contribution at z~4 has a lower limit of 3x10^(-3)Msun/yr h^3 Mpc^(-3), ~>20% of the estimated total SFRD. This highlights the importance of extremely infrared-luminous galaxies with high star formation rates to the build-up of stellar mass, even at the earliest epochs.

preprint2012arXiv

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 500um, taken with the Keck I Low Resolution Imaging Spectrometer (LRIS) and the Keck II DEep Imaging Multi-Object Spectrograph (DEIMOS). The redshift distribution of these SPIRE sources from the Herschel Multitiered Extragalactic Survey (HerMES) 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 (<delta_z>/(1+z)>=0.29) than is seen in non-infrared selected samples, likely due to enhanced star formation rates and dust obscuration in infrared-selected galaxies. We estimate that the vast majority (72-83%) of z<2 Herschel-selected galaxies would drop out of traditional submillimeter surveys at 0.85-1mm. 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 24um 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. [abridged]

preprint2012arXiv

Accreting SMBHs in the COSMOS field and the connection to their host galaxies

Using the wide multi-band photometry available in the COSMOS field we explore the host galaxy properties of a large sample of Active Galactic Nuclei (AGN) obtained by combining X-ray and optical spectroscopic selections. Based on a careful study of their Spectral Energy Distribution (SED), which has been parametrized using a 2-component (AGN+galaxy) model fit, we derived dust-corrected rest-frame magnitudes, colors, stellar masses and star formation rates (SFRs). We find that AGN hosts span a large range of stellar masses and SFRs. No color-bimodality is seen at any redshift in the AGN hosts, which are found to be mainly massive, red galaxies. Once accounting for the color-mass degeneracy in well defined mass-matched samples, we find a residual marginal enhancement of AGN incidence in redder galaxies with lower specific star formation rates, and we argue that this result might emerge because of our ability to properly account for AGN light contamination and dust extinction. Interestingly, we find that the probability for a galaxy to host a black hole growing at any given "specific accretion rate" (i.e. the ratio of X-ray luminosity to the host stellar mass) is almost independent of the host galaxy mass, while it decreases as a power-law with Lx/M. By analyzing the normalization of such probability distribution, we show how the incidence of AGN increases with redshift as rapidly as (1+z)^4, in close resemblance with the overall evolution of the specific star formation rate of the entire galaxy population. Although AGN activity and star formation in galaxies do appear to have a common triggering mechanism, at least in a statistical sense, within the COSMOS sample we do not find strong evidence of any 'smoking gun' signaling powerful AGN influence on the star-forming properties of their hosts galaxies.

preprint2012arXiv

Bolometric luminosities and Eddington ratios of X-ray selected Active Galactic Nuclei in the XMM-COSMOS survey

Bolometric luminosities and Eddington ratios of both X-ray selected broad-line (Type-1) and narrow-line (Type-2) AGN from the XMM-Newton survey in the COSMOS field are presented. The sample is composed by 929 AGN (382 Type-1 AGN and 547 Type-2 AGN) and it covers a wide range of redshifts, X-ray luminosities and absorbing column densities. About 65% of the sources are spectroscopically identified as either Type-1 or Type-2 AGN (83% and 52% respectively), while accurate photometric redshifts are available for the rest of the sample. The study of such a large sample of X-ray selected AGN with a high quality multi-wavelength coverage from the far-infrared (now with the inclusion of Herschel data at 100 micron and 160 micron) to the optical-UV allows us to obtain accurate estimates of bolometric luminosities, bolometric corrections and Eddington ratios. The kbol-Lbol relations derived in the present work are calibrated for the first time against a sizable AGN sample, and rely on observed redshifts, X-ray luminosities and column density distributions. We find that kbol is significantly lower at high Lbol with respect to previous estimates by Marconi et al. (2004) and Hopkins et al. (2007). Black hole masses and Eddington ratios are available for 170 Type-1 AGN, while black hole masses for Type-2 AGN are computed for 481 objects using the black hole mass-stellar mass relation and the morphological information. We confirm a trend between kbol and lambda_Edd, with lower hard X-ray bolometric corrections at lower Eddington ratios for both Type-1 and Type-2 AGN. We find that, on average, Eddington ratio increases with redshift for all Types of AGN at any given Mbh, while no clear evolution with redshift is seen at any given Lbol.

preprint2012arXiv

Comparison of star formation rates from Halpha and infrared luminosities as seen by Herschel

We empirically test the relation between the SFR(LIR) derived from the infrared luminosity, LIR, and the SFR(Ha) derived from the Ha emission line luminosity using simple conversion relations. We use a sample of 474 galaxies at z = 0.06 - 0.46 with both Ha detection (from 20k zCOSMOS survey) and new far-IR Herschel data (100 and 160 μm). We derive SFR(Ha) from the Ha extinction corrected emission line luminosity. We find a very clear trend between E(B - V) and LIR that allows to estimate extinction values for each galaxy even if the Ha emission line measurement is not reliable. We calculate the LIR by integrating from 8 up to 1000 μm the SED that is best fitting our data. We compare SFR(Ha) with the SFR(LIR). We find a very good agreement between the two SFR estimates, with a slope of m = 1.01 \pm 0.03 in the SFR(LIR) vs SFR(Ha) diagram, a normalization constant of a = -0.08 \pm 0.03 and a dispersion of sigma = 0.28 dex.We study the effect of some intrinsic properties of the galaxies in the SFR(LIR)-SFR(Ha) relation, such as the redshift, the mass, the SSFR or the metallicity. The metallicity is the parameter that affects most the SFR comparison. The mean ratio of the two SFR estimators log[SFR(LIR)/SFR(Ha)] varies by approx. 0.6 dex from metal-poor to metal-rich galaxies (8.1 < log(O/H) + 12 < 9.2). This effect is consistent with the prediction of a theoretical model for the dust evolution in spiral galaxies. Considering different morphological types, we find a very good agreement between the two SFR indicators for the Sa, Sb and Sc morphologically classified galaxies, both in slope and normalization. For the Sd, irregular sample (Sd/Irr), the formal best-fit slope becomes much steeper (m = 1.62 \pm 0.43), but it is still consistent with 1 at the 1.5 sigma level, because of the reduced statistics of this sub-sample.

preprint2012arXiv

Cosmic Evolution of Star-Formation Enhancement in Close Major-Merger Galaxy Pairs Since z = 1

The infrared (IR) emission of M_* galaxies (10^{10.4} < M_{star} < 10^{11.0} M_\sun) in galaxy pairs, derived using data obtained in Herschel (PEP/HerMES) and Spitzer (S-COSMOS) surveys, is compared to that of single disk galaxies in well matched control samples to study the cosmic evolution of the star-formation enhancement induced by galaxy-galaxy interaction. Both the mean IR SED and mean IR luminosity of star-forming galaxies (SFGs) in SFG+SFG (S+S) pairs in the redshift bin of 0.6 < z < 1 are consistent with no star-formation enhancement. SFGs in S+S pairs in a lower redshift bin of 0.2 < z < 0.6 show marginal evidence for a weak star-formation enhancement. Together with the significant and strong sSFR enhancement shown by SFGs in a local sample of S+S pairs (obtained using previously published Spitzer observations), our results reveal a trend for the star-formation enhancement in S+S pairs to decrease with increasing redshift. Between z=0 and z=1, this decline of interaction-induced star-formation enhancement occurs in parallel with the dramatic increase (by a factor of ~10) of the sSFR of single SFGs, both can be explained by the higher gas fraction in higher z disks. SFGs in mixed pairs (S+E pairs) do not show any significant star-formation enhancement at any redshift. The difference between SFGs in S+S pairs and in S+E pairs suggests a modulation of the sSFR by the inter-galactic medium IGM in the dark matter halos (DMH) hosting these pairs.

preprint2012arXiv

Deep observations of CO line emission from star-forming galaxies in a cluster candidate at z=1.5

We report results from a deep Jansky Very Large Array (JVLA) search for CO 1-0 line emission from galaxies in a candidate galaxy cluster at z~1.55 in the COSMOS field. We target 4 galaxies with optical spectroscopic redshifts in the range z=1.47-1.59. Two of these 4 galaxies, ID51613 and ID51813, are nominally detected in CO line emission at the 3-4 sigma level. We find CO luminosities of 2.4x10^10 K km/s pc^2 and 1.3x10^10 K km/s pc^2, respectively. Taking advantage from the clustering and 2-GHz bandwidth of the JVLA, we perform a search for emission lines in the proximity of optical sources within the field of view of our observations. We limit our search to galaxies with K<23.5 (AB) and z_phot=1.2-1.8. We find 2 bright optical galaxies to be associated with significant emission line peaks (>4 sigma) in the data cube, which we identify with the CO line emission. To test the reliability of the line peaks found, we performed a parallel search for line peaks using a Bayesian inference method. Monte Carlo simulations show that such associations are statistically significant, with probabilities of chance association of 3.5% and 10.7% for ID 51207 and ID 51380, respectively. Modeling of their optical/IR SEDs indicates that the CO detected galaxies and candidates have stellar masses and SFRs in the range (0.3-1.1)x10^11 M_sun and 60-160 M_sun/yr, with SFEs comparable to that found in other star-forming galaxies at similar redshifts. By comparing the space density of CO emitters derived from our observations with the space density derived from previous CO detections at z~1.5, and with semi-analytic predictions for the CO luminosity function, we suggest that the latter tend to underestimate the number of CO galaxies detected at high-redshift. Finally, we argue about the benefits of future blind CO searches in clustered fields with upcoming submm/radio facilities.

preprint2012arXiv

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 mass-matched control inactive galaxies,including star forming and quiescent sources, at 0.5<z<2.5. Recent observations carried out by PACS, the 60-210um Herschel photometric camera, in GOODS-S, GOODS-N and COSMOS allow us to unbiasedly estimate the far-IR luminosity, and hence the SF properties, of the two samples. Accurate AGN host stellar masses are measured by decomposing their total emission into the stellar and nuclear components. We find a higher average SF activity in AGN hosts with respect to non-AGNs. The level of SF enhancement is modest (~0.26dex at ~3sigma) at low X-ray luminosities (Lx<~10^43.5erg/s) and more pronounced (0.56dex at >10sigma) for bright AGNs. However, when comparing to star forming galaxies only, AGN hosts are 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, while luminous AGNs co-evolve with their hosts through periods of enhanced AGN activity and SF, possibly through major mergers. While an increased SF with respect to non-AGNs 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.

preprint2012arXiv

Herschel observations of a z~2 stellar mass selected galaxy sample drawn from the GOODS NICMOS Survey

We present a study of the far-IR properties of a stellar mass selected sample of 1.5 < z < 3 galaxies with log(M_*/M_sun) > 9.5 drawn from the GOODS NICMOS Survey (GNS), the deepest H-band Hubble Space Telescope survey of its type prior to the installation of WFC3. We use far-IR and sub-mm data from the PACS and SPIRE instruments on-board Herschel, taken from the PACS Evolutionary Probe (PEP) and Herschel Multi-Tiered Extragalactic Survey (HerMES) key projects respectively. We find a total of 22 GNS galaxies, with median log(M_*/M_sun) = 10.8 and z = 2.0, associated with 250 um sources detected with SNR > 3. We derive mean total IR luminosity log L_IR (L_sun) = 12.36 +/- 0.05 and corresponding star formation rate SFR_(IR+UV) = (280 +/- 40) M_sun/yr for these objects, and find them to have mean dust temperature T_dust ~ 35 K. We find that the SFR derived from the far-IR photometry combined with UV-based estimates of unobscured SFR for these galaxies is on average more than a factor of 2 higher than the SFR derived from extinction corrected UV emission alone, although we note that the IR-based estimate is subject to substantial Malmquist bias. To mitigate the effect of this bias and extend our study to fainter fluxes, we perform a stacking analysis to measure the mean SFR in bins of stellar mass. We obtain detections at the 2-4 sigma level at SPIRE wavelengths for samples with log(M_*/M_sun) > 10. In contrast to the Herschel detected GNS galaxies, we find that estimates of SFR_(IR+UV) for the stacked samples are comparable to those derived from extinction corrected UV emission, although the uncertainties are large. We find evidence for an increasing fraction of dust obscured star formation with stellar mass, finding SFR_IR/SFR_UV \propto M_*^{0.7 +/- 0.2}, which is likely a consequence of the mass--metallicity relation.

preprint2012arXiv

PACS-Herschel FIR detections of Lyman-alpha emitters at 2.0<z<3.5

In this work we analyze the physical properties of a sample of 56 spectroscopically selected star-forming (SF) Ly$α$ emitting galaxies at 2.0$\lesssim$z$\lesssim$3.5 using both a spectral energy distribution (SED) fitting procedure from rest-frame UV to mid-IR and direct 160$μ$m observations taken with the Photodetector Array Camera & Spectrometer (PACS) instrument onboard \emph{Herschel Space Observatory}. We define LAEs as those Ly$α$ emitting galaxies whose rest-frame Ly$α$ equivalent widths (Ly$α$ EW$_{rest-frame}$) are above 20Å, the typical threshold in narrow-band searches. Ly$α$ emitting galaxies with Ly$α$ EW$_{rest-frame}$ are called non-LAEs. As a result of an individual SED fitting for each object, we find that the studied sample of LAEs contains galaxies with ages mostly below 100Myr and a wide variety of dust attenuations, SFRs, and stellar masses. The heterogeneity in the physical properties is also seen in the morphology, ranging from bulge-like galaxies to highly clumpy systems. In this way, we find that LAEs at 2.0$\lesssim$z$\lesssim$3.5 are very diverse, and do not have a bimodal nature, as suggested in previous works. Furthermore, the main difference between LAEs and non-LAEs is their dust attenuation, because LAEs are not as dusty as non-LAEs. On the FIR side, four galaxies of the sample (two LAEs and two non-LAEs) have PACS-FIR counterparts. Their total IR luminosity place all of them in the ULIRG regime and are all dusty objects, with A$_{1200}$$\gtrsim$4mag. This is an indication from direct FIR measurements that dust and Ly$α$ emission are not mutually exclusive. This population of red and dusty LAEs is not seen at z$\sim$0.3, suggesting an evolution with redshift of the IR nature of galaxies selected via their Ly$α$ emission.

preprint2012arXiv

Physical properties of Lyman-alpha emitters at $z\sim 0.3$ from UV-to-FIR measurements

The analysis of the physical properties of low-redshift Ly$α$ emitters (LAEs) can provide clues in the study of their high-redshift analogues. At $z \sim 0.3$, LAEs are bright enough to be detected over almost the entire electromagnetic spectrum and it is possible to carry out a more precise and complete study than at higher redshifts. In this study, we examine the UV and IR emission, dust attenuation, SFR and morphology of a sample of 23 GALEX-discovered star-forming (SF) LAEs at $z \sim 0.3$ with direct UV (GALEX), optical (ACS) and FIR (PACS and MIPS) data. Using the same UV and IR limiting luminosities, we find that LAEs at $z\sim 0.3$ tend to be less dusty, have slightly higher total SFRs, have bluer UV continuum slopes, and are much smaller than other galaxies that do not exhibit Ly$α$ emission in their spectrum (non-LAEs). These results suggest that at $z \sim 0.3$ Ly$α$ photons tend to escape from small galaxies with low dust attenuation. Regarding their morphology, LAEs belong to Irr/merger classes, unlike non-LAEs. Size and morphology represent the most noticeable difference between LAEs and non-LAEs at $z \sim 0.3$. Furthermore, the comparison of our results with those obtained at higher redshifts indicates that either the Ly$α$ technique picks up different kind of galaxies at different redshifts or that the physical properties of LAEs are evolving with redshift.

preprint2012arXiv

The AGN content in luminous IR galaxies at z\sim2 from a global SED analysis including Herschel data

We use Herschel-PACS far-infrared data, combined with previous multi-band information and mid-IR spectra, to properly account for the presence of an active nucleus and constrain its energetic contribution in luminous infrared (IR) sources at z\sim2. The sample is composed of 24 sources in the GOODS-South field, with typical IR luminosity of 10^{12} Lo. Data from the 4 Ms Chandra X-ray imaging in this field are also used to identify and characterize AGN emission. We reproduce the observed spectral energy distribution (SED), decomposed into a host-galaxy and an AGN component. A smooth-torus model for circum-nuclear dust is used to account for the direct and re-processed contribution from the AGN. We confirm that galaxies with typical L_{8-1000um}\sim10^{12}Lo at z\sim2 are powered predominantly by star-formation. An AGN component is present in nine objects (\sim35% of the sample) at the 3sigma confidence level, but its contribution to the 8-1000 um emission accounts for only \sim5% of the energy budget. The AGN contribution rises to \sim23% over the 5-30 um range (in agreement with Spitzer IRS results) and to \sim60% over the narrow 2-6 um range. The presence of an AGN is confirmed by X-ray data for 3 (out of nine) sources, with X-ray spectral analysis indicating the presence of significant absorption, i.e. NH\sim10^{23}-10^{24} cm^{-2}. An additional source shows indications of obscured AGN emission from X-ray data. The comparison between the mid-IR--derived X-ray luminosities and those obtained from X-ray data suggests that obscuration is likely present also in the remaining six sources that harbour an AGN according to the SED-fitting analysis.

preprint2012arXiv

The complex physics of dusty star-forming galaxies at high redshifts as revealed by Herschel and Spitzer

We combine far-infrared photometry from Herschel (PEP/HERMES) with deep mid-infrared spectroscopy from Spitzer to investigate the nature and the mass assembly history of a sample of 31 Luminous and Ultraluminous Infrared Galaxies at z~1 and 2 selected in GOODS-S with 24 $μ$m fluxes between 0.2 and 0.5 mJy. We model the data with a self-consistent physical model (GRASIL) which includes a state-of-the-art treatment of dust extinction and reprocessing. We find that all of our galaxies appear to require massive populations of old (>1 Gyr) stars and, at the same time, to host a moderate ongoing activity of SF (SFR < 100 M$_{\odot}$/yr). The bulk of the stars appear to have been formed a few Gyr before the observation in essentially all cases. Only five galaxies of the sample require a recent starburst superimposed on a quiescent star formation history (SFH). We also find discrepancies between our results and those based on optical-only SED fitting for the same objects; by fitting their observed Spectral Energy Distributions with our physical model we find higher extinctions (by $Δ$A_{V} ~ 0.81 and 1.14) and higher stellar masses (by $Δ$Log(M*) ~ 0.16 and 0.36 dex) for z~1 and z~2 (U)LIRGs, respectively. The stellar mass difference is larger for the most dust obscured objects. We also find lower SFRs than those computed from L_{IR} using the Kennicutt relation due to the significant contribution to the dust heating by intermediate-age stellar populations through 'cirrus' emission (~73% and ~66% of total L_{IR} for z~1 and z~2 (U)LIRGs, respectively).

preprint2011arXiv

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 micron, 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 ~1 mJy to few hundreds mJy. Stacking of 24 micron 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 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. The total CIB surface brightness emitted above PEP 3 sigma flux limits is 4.52 +/- 1.18, 8.35 +/- 0.95 and 9.49 +/- 0.59 [nW/m2/sr] at 70, 100, and 160 micron, respectively. These values correspond to 58 +/- 7% and 74 +/- 5% of the COBE/DIRBE CIB direct measurements at 100 and 160 micron. 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 micron CIB is mainly produced by z<=1.0 objects, the contribution of z>1.0 sources reaches 50% at 160 micron. Most of the CIB resolved in the three PACS bands was emitted by galaxies with infrared luminosities in the range 1e11-1e12 L(sun).

preprint2011arXiv

FIR measurements of Ly-$α$ emitters at z$\lesssim$1.0: dust attenuation from PACS-\emph{Herschel}

One remaining open question regarding the physical properties of Ly$α$ emitters (LAEs) is their dust content and its evolution with redshift. The variety of results is large and with those reported by now is difficult to establish clear relations between dust, other fundamental parameters of galaxies (star-formation rate, metallicity or age) and redshift. In this Letter, we report \emph{Herschel} PACS-100$μ$m, PACS-160$μ$m and \emph{Spitzer} MIPS-24$μ$m detections of a sample of spectroscopically GALEX-selected LAEs at z$\sim$0.3 and $\sim$1.0. Five out of ten and one out of two LAEs are detected in, at least, one PACS band at z$\sim$0.3 and $\sim$1.0, respectively. These measurements have a great importance given that they allow us to quantify, for the first time, the dust content of LAEs from direct FIR observations. MIPS-24$μ$m detections allow us to determine IR properties of the PACS-undetected LAEs. We obtain that mid-IR/FIR detected star-forming (SF) LAEs at z$\sim$0.3 have dust content within 0.75$\lesssim$ $A_{1200Å}$ $\lesssim$2.0, with a median value of A$_{1200\textrmÅ}$$\sim$1.1. This range broadens out to 0.75$\lesssim$ $A_{1200Å}$ $\lesssim$2.5 when considering those LAEs at z$\sim$1.0. Only one SF LAE is undetected both in MIPS-24$μ$m and PACS, with $A_{1200Å}$ $\lesssim$0.75. These results seem to be larger than those reported for high-redshift LAEs and, therefore, although an evolutionary trend is not clearly seen, it could point out that low-redshift LAEs are dustier than high-redshift ones. However, the diverse methods used could introduce a systematic offset in the results.

preprint2011arXiv

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.

preprint2011arXiv

On the nature of the absorber in IRAS 09104+4109: the X-ray and mid-infrared view

We present a long (~76 ks) Chandra observation of IRAS 09104+4109, a hyper-luminous galaxy, optically classified as a Type 2 AGN hosted in a cD galaxy in a cluster at z=0.442. We also report on the results obtained by fitting its broad-band spectral energy distribution. The Compton-thick nature of this source (which has been often referred to as an "archetype" of Compton-thick Type 2 quasars) was formerly claimed on the basis of its marginal detection in the PDS instrument onboard BeppoSAX, being then disputed using XMM-Newton data. Both Chandra analysis and optical/mid-IR spectral fitting are consistent with the presence of heavy (~1-5 10^{23} cm^{-2}) but not extreme (Compton-thick) obscuration. However, using the mid-IR and the [OIII] emission as proxies of the nuclear hard X-ray luminosity suggests the presence of heavier obscuration. The 54-month Swift BAT map shows excess hard X-ray emission likely related to a nearby (z=0.009) Type 2 AGN, close enough to IRAS 09104+4109 to significantly enhance and contaminate its emission in the early BeppoSAX PDS data.

preprint2011arXiv

PACS Evolutionary Probe (PEP) - A Herschel Key Program

Deep far-infrared photometric surveys studying galaxy evolution and the nature of the cosmic infrared background are a key strength of the Herschel mission. We describe the scientific motivation for the PACS Evolutionary Probe (PEP) guaranteed time key program and its role in the complement of Herschel surveys, and the field selection which includes popular multiwavelength fields such as GOODS, COSMOS, Lockman Hole, ECDFS, EGS. We provide an account of the observing strategies and data reduction methods used. An overview of first science results illustrates the potential of PEP in providing calorimetric star formation rates for high redshift galaxy populations, thus testing and superseeding previous extrapolations from other wavelengths, and enabling a wide range of galaxy evolution studies.

preprint2011arXiv

The evolution of quiescent galaxies at high redshift (z > 1.4)

We have studied the evolution of high redshift quiescent galaxies over an effective area of ~1.7 deg^2 in the COSMOS field. Galaxies have been divided according to their star-formation activity and the evolution of the different populations has been investigated in detail. We have studied an IRAC (mag_3.6 < 22.0) selected sample of ~18000 galaxies at z > 1.4 with multi-wavelength coverage. We have derived accurate photometric redshifts (sigma=0.06) and other important physical parameters through a SED-fitting procedure. We have divided our sample into actively star-forming, intermediate and quiescent galaxies depending on their specific star formation rate. We have computed the galaxy stellar mass function of the total sample and the different populations at z=1.4-3.0. We have studied the properties of high redshift quiescent galaxies finding that they are old (1-4 Gyr), massive (log(M/M_sun)~10.65), weakly star forming stellar populations with low dust extinction (E(B-V) < 0.15) and small e-folding time scales (tau ~ 0.1-0.3 Gyr). We observe a significant evolution of the quiescent stellar mass function from 2.5 < z < 3.0 to 1.4 < z < 1.6, increasing by ~ 1 dex in this redshift interval. We find that z ~ 1.5 is an epoch of transition of the GSMF. The fraction of star-forming galaxies decreases from 60% to 20% from z ~ 2.5-3.0 to z ~ 1.4-1.6 for log(M/M_sun) > 11, while the quiescent population increases from 10% to 50% at the same redshift and mass intervals. We compare the fraction of quiescent galaxies derived with that predicted by theoretical models and find that the Kitzbichler & White (2007) model is the one that better reproduces the data. Finally, we calculate the stellar mass density of the star-forming and quiescent populations finding that there is already a significant number of quiescent galaxies at z > 2.5 (rho~6.0 MsunMpc^-3).

preprint2011arXiv

The evolution of the star formation activity per halo mass up to redshift ~ 1.6 as seen by Herschel

Star formation in massive galaxies is quenched at some point during hierarchical mass assembly. To understand where and when the quenching processes takes place, we study the evolution of the total star formation rate per unit total halo mass (Σ(SFR/M)) in three different mass scales: low mass halos (field galaxies), groups, and clusters, up to a redshift ~1.6. We use deep far-infrared PACS data at 100 and 160 um to accurately estimate the total star formation rate of the Luminous Infrared Galaxy population of 9 clusters with mass ~10^{15} M_{\odot}, and 9 groups/poor clusters with mass ~ 5 x 10^{13} M_{\odot}. Estimates of the field Σ(SFR/M) are derived from the literature, by dividing the star formation rate density by the mean comoving matter density of the universe. The field Σ(SFR/M) increases with redshift up to z~1 and it is constant thereafter. The evolution of the Σ(SFR/M)-z relation in galaxy systems is much faster than in the field. Up to redshift z~0.2, the field has a higher Σ(SFR/M) than galaxy groups and galaxy clusters. At higher redshifts, galaxy groups and the field have similar Σ(SFR/M), while massive clusters have significantly lower Σ(SFR/M) than both groups and the field. There is a hint of a reversal of the SFR activity vs. environment at z~1.6, where the group Σ(SFR/M) lies above the field Σ(SFR/M)-z relation. We discuss possible interpretations of our results in terms of the processes of downsizing, and star-formation quenching.

preprint2011arXiv

The impact of evolving infrared spectral energy distributions of galaxies on star formation rate estimates

We combine Herschel-PACS data from the PEP program with Spitzer 24 um and 16 um photometry and ultra deep 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 around the main sequence (the redshift-dependent relation of star formation rate and stellar mass). Our deep data confirm from individual far-infrared detections that z~2 star formation rates are overestimated if based on 24 um fluxes and SED templates that are calibrated via local trends with luminosity. Galaxies with similar ratios of rest-frame nuLnu(8) to 8-1000 um 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 nuLnu(8)/LIR regardless of LIR and redshift, up to z~2.5, and nuLnu(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 um 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 PAH to LIR rather than continua of hidden AGN. Alternatively, we discuss the dependence of nuLnu(8)/LIR on LIR. Same nuLnu(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 in hand. For most of those z~2 star forming galaxies that also host an AGN, the mid-infrared is dominated by the star forming component.

preprint2011arXiv

The lesser role of starbursts for star formation at z=2

Two main modes of star formation are know to control the growth of galaxies: a relatively steady one in disk-like galaxies, defining a tight star formation rate (SFR)-stellar mass sequence, and a starburst mode in outliers to such a sequence which is generally interpreted as driven by merging. Such starburst galaxies are rare but have much higher SFRs, and it is of interest to establish the relative importance of these two modes. PACS/Herschel observations over the whole COSMOS and GOODS-South fields, in conjunction with previous optical/near-IR data, have allowed us to accurately quantify for the first time the relative contribution of the two modes to the global SFR density in the redshift interval 1.5<z<2.5, i.e., at the cosmic peak of the star formation activity. The logarithmic distributions of galaxy SFRs at fixed stellar mass are well described by Gaussians, with starburst galaxies representing only a relatively minor deviation that becomes apparent for SFRs more than 4 times higher than on the main sequence. Such starburst galaxies represent only 2% of mass-selected star forming galaxies and account for only 10% of the cosmic SFR density at z~2. Only when limited to SFR>1000M(sun)/yr, off-sequence sources significantly contribute to the SFR density (46+/-20%). We conclude that merger-driven starbursts play a relatively minor role for the formation of stars in galaxies, whereas they may represent a critical phase towards the quenching of star formation and morphological transformation in galaxies.

preprint2011arXiv

The PEP survey: clustering of infrared-selected galaxies and structure formation at z~2 in the GOODS South

ABRIDGED-This paper presents the first direct estimate of the 3D clustering properties of far-infrared sources up to z~3. This has been possible thanks to the Pacs Evolutionary Probe (PEP) survey of the GOODS South field performed with the PACS instrument onboard the Herschel Satellite. An analysis of the two-point correlation function over the whole redshift range spanned by the data reports for the correlation length, r_0~6.3 Mpc and r_0~6.7 Mpc, respectively at 100um and 160um, corresponding to dark matter halo masses M>~10^{12.4} M_sun. Objects at z~2 instead seem to be more strongly clustered, with r_0~19 Mpc and r_0~17 Mpc in the two considered PACS channels. This dramatic increase of the correlation length between z~1 and z~2 is connected with the presence of a wide, M>~10^{14} M_sun, filamentary structure which includes more than 50% of the sources detected at z~2. An investigation of the properties of such sources indicates the possibility for boosted star-forming activity in those which reside within the overdense environment with respect of more isolated galaxies found in the same redshift range. Lastly, we also present our results on the evolution of the relationship between luminous and dark matter in star-forming galaxies between z~1 and z~2. We find that the increase of (average) stellar mass in galaxies <M*> between z~1 and z~2 is about a factor 10 lower than that of the dark matter haloes hosting such objects (<M*>[z~1]/<M*>[z~2] ~ 0.4 vs M_{halo}[z~1]/M_{halo}[z~2] ~ 0.04). Our findings agree with the evolutionary picture of downsizing whereby massive galaxies at z~2 were more actively forming stars than their z~1 counterparts, while at the same time contained a lower fraction of their mass in the form of luminous matter.

preprint2011arXiv

The redshift and nature of AzTEC/COSMOS 1: A starburst galaxy at z=4.6

Based on broad/narrow-band photometry and Keck DEIMOS spectroscopy we report a redshift of z=4.64-0.08+0.06 for AzTEC/COSMOS 1, the brightest sub-mm galaxy in the AzTEC/COSMOS field. In addition to the COSMOS-survey X-ray to radio data, we report observations of the source with Herschel/PACS (100, 160 micron), CSO/SHARC II (350 micron), CARMA and PdBI (3 mm). We do not detect CO(5-4) line emission in the covered redshift ranges, 4.56-4.76 (PdBI/CARMA) and 4.94-5.02 (CARMA). If the line is within this bandwidth, this sets 3sigma upper limits on the gas mass to <~8x10^9 M_Sol and <~5x10^10 M_Sol, respectively (assuming similar conditions as observed in z~2 SMGs). This could be explained by a low CO-excitation in the source. Our analysis of the UV-IR spectral energy distribution of AzTEC 1 shows that it is an extremely young (<~50 Myr), massive (M*~10^11 M_Sol), but compact (<~2 kpc) galaxy forming stars at a rate of ~1300 M_Sol/yr. Our results imply that AzTEC 1 is forming stars in a 'gravitationally bound' regime in which gravity prohibits the formation of a superwind, leading to matter accumulation within the galaxy and further generations of star formation.

preprint2010arXiv

Discovery of Compton-thick quasars in the Sloan Digital Sky Survey

We present new and archival Chandra snapshot (10 ks each) observations of 15 optically identified (from the Sloan Digital Sky Survey, SDSS) Type 2 quasars at z=0.40-0.73. When combined with existing X-ray data, this work provides complete X-ray coverage for all 25 radio-quiet Type 2 quasars with logL_[OIII]>9.28 L_sun from Zakamska et al. (2003). Two targets out of 15 were not detected by Chandra and most of the remaining sources are X-ray weak, with nine having less than 10 counts in the 0.5-8keV band. Low-to-moderate quality spectral analysis was limited to three sources, whose properties are consistent with the presence of column densities in the range NH~10^22-10^23 cm^-2 in the source rest frame. If the [OIII] luminosity is a reliable proxy for the intrinsic X-ray luminosity, the current X-ray data indicate that Compton-thick quasars may hide among ~65 per cent of the SDSS Type 2 quasar population (L_{X, meas}/L_{X, [OIII]}<0.01); however, since the Type 2 quasar sample is selected on [OIII] luminosity, the estimated Compton-thick quasar fraction may be overestimated. Using archival Spitzer observations, we find that ~50 per cent of SDSS Type 2 quasars appear to be obscured by Compton-thick material based on both the L_{X, meas}/L_{X, mid-IR} (where mid-IR corresponds to rest-frame 12.3 micron) and L_{X, meas}/L_{X, [OIII]} ratios. We use this information to provide an estimate of the Compton-thick quasar number density at z=0.3-0.8, which we find is in broad agreement with the expectations from X-ray background models.

preprint2010arXiv

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 um) 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 sigma depths of ~3.0 mJy at 100 um and ~5.7 mJy at 160 um. At these levels, source densities are 40 and 18 beams/source, respectively, thus hitting the confusion limit at 160 um. 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 (L(IR)~1e11 Lsun), whose SEDs resemble late-spiral galaxies, peaking at ~130 um 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 um, resolving ~45% and ~52% of the CIB, respectively. When stacking 24 um sources, the inferred CIB lies within 1.1 sigma and 0.5 sigma 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 um sources found at higher redshift than 100 um ones.

preprint2010arXiv

Evolution of Dust Temperature of Galaxies through Cosmic Time as seen by Herschel

We study the dust properties of galaxies in the redshift range 0.1<z<2.8 observed by the Herschel Space Observatory in the field of the Great Observatories Origins Deep Survey-North as part of PEP and HerMES key programmes. Infrared (IR) luminosity (L_IR) and dust temperature (T_dust) of galaxies are derived from the spectral energy distribution (SED) fit of the far-infrared (FIR) flux densities obtained with PACS and SPIRE instruments onboard Herschel. As a reference sample, we also obtain IR luminosities and dust temperatures of local galaxies at z<0.1 using AKARI and IRAS data in the field of the Sloan Digital Sky Survey. We compare the L_IR-T_dust relation between the two samples and find that: the median T_dust of Herschel-selected galaxies at z>0.5 with L_IR>5x10^{10} L_\odot, appears to be 2-5 K colder than that of AKARI-selected local galaxies with similar luminosities; and the dispersion in T_dust for high-z galaxies increases with L_IR due to the existence of cold galaxies that are not seen among local galaxies. We show that this large dispersion of the L_IR-T_dust relation can bridge the gap between local star-forming galaxies and high-z submillimeter galaxies (SMGs). We also find that three SMGs with very low T_dust (<20 K) covered in this study have close neighbouring sources with similar 24-\mum brightness, which could lead to an overestimation of FIR/(sub)millimeter fluxes of the SMGs.

preprint2010arXiv

Far-infrared properties of submillimeter and optically faint radio galaxies

We use deep observations obtained with the Photodetector Array Camera and Spectrometer (PACS) onboard the Herschel space observatory to study the far-infrared (FIR) properties of submm and optically faint radio galaxies (SMGs and OFRGs). From literature we compiled a sample of 35 securely identified SMGs and nine OFRGs located in the GOODS-N and the A2218 fields. This sample is cross-matched with our PACS 100 um and 160 um multi-wavelength catalogs. About half of the galaxies in our sample are detected with PACS. The dust temperatures and the infrared luminosities of our galaxies are derived by fitting their PACS and SCUBA 850 um (only the upper limits for the OFRGs) flux densities with a single modified (beta=1.5) black body function. The median dust temperature of our SMG sample is T=36+/-8K while for our OFRG sample it is T=47+/-3K. For both samples, median dust temperatures derived from Herschel data agree well with previous estimates. In particular, Chapman et al. (2005) found a dust temperature of T=36+/-7K for a large sample of SMGs assuming the validity of the FIR/radio correlation. The agreement between our studies confirms that the local FIR/radio correlation effectively holds at high redshift even though we find <q>=2.17+/-0.19, a slightly lower value than that observed in local systems. The median IR luminosities of SMGs and OFRGs are 4.6*10^12 Lsun and 2.6*10^12 Lsun, respectively. We note that for both samples the IR luminosity estimates from the radio part of the spectral energy distribution are accurate, while estimates from the mid-IR are considerably (x3) more uncertain. Our observations confirm the remarkably high luminosities of SMGs and thus imply median star-formation rates of 960Msun yr^-1 for SMGs with S(850um)>5mJy and 460Msun yr^-1 for SMGs with S(850um)>2mJy, assuming a Chabrier IMF and no dominant AGN contribution to the far-infrared luminosity.

preprint2010arXiv

Herschel deep far-infrared counts through the Abell 2218 cluster-lens

Gravitational lensing by massive galaxy clusters allows study of the population of intrinsically faint infrared galaxies that lie below the sensitivity and confusion limits of current infrared and submillimeter telescopes. We present ultra-deep PACS 100 and 160 microns observations toward the cluster lens Abell 2218, to penetrate the Herschel confusion limit. We derive source counts down to a flux density of 1 mJy at 100 microns and 2 mJy at 160 microns, aided by strong gravitational lensing. At these levels, source densities are 20 and 10 beams/source in the two bands, approaching source density confusion at 160 microns. The slope of the counts below the turnover of the Euclidean-normalized differential curve is constrained in both bands and is consistent with most of the recent backwards evolutionary models. By integrating number counts over the flux range accessed by Abell 2218 lensing (0.94-35 mJy at 100 microns and 1.47-35 mJy at 160 microns, we retrieve a cosmic infrared background (CIB) surface brightness of ~8.0 and ~9.9 nW m^-2 sr^-1, in the respective bands. These values correspond to 55% (+/- 24%) and 77% (+/- 31%) of DIRBE direct measurements. Combining Abell 2218 results with wider/shallower fields, these figures increase to 62% (+/- 25%) and 88% (+/- 32%) CIB total fractions, resolved at 100 and 160 microns, disregarding the high uncertainties of DIRBE absolute values.

preprint2010arXiv

Herschel FIR counterparts of selected Ly-alpha emitters at z~2.2. Fast evolution since z~3 or missed obscured AGNs?

Ly-alpha emitters (LAEs) are seen everywhere in the redshift domain from local to z~7. Far-infrared (FIR) counterparts of LAEs at different epochs could provide direct clues on dust content, extinction, and spectral energy distribution (SED) for these galaxies. We search for FIR counterparts of LAEs that are optically detected in the GOODS-North field at redshift z~2.2 using data from the Herschel Space Telescope with the Photodetector Array Camera and Spectrometer (PACS). The LAE candidates were isolated via color-magnitude diagram using the medium-band photometry from the ALHAMBRA Survey, ancillary data on GOODS-North, and stellar population models. According to the fitting of these spectral synthesis models and FIR/optical diagnostics, most of them seem to be obscured galaxies whose spectra are AGN-dominated. From the analysis of the optical data, we have observed a fraction of AGN or composite over source total number of ~0.75 in the LAE population at z~2.2, which is marginally consistent with the fraction previously observed at z=2.25 and even at low redshift (0.2<z<0.45), but significantly different from the one observed at redshift ~3, which could be compatible either with a scenario of rapid change in the AGN fraction between the epochs involved or with a non detection of obscured AGN in other z=2-3 LAE samples due to lack of deep FIR observations. We found three robust FIR (PACS) counterparts at z~2.2 in GOODS-North. This demonstrates the possibility of finding dust emission in LAEs even at higher redshifts.

preprint2010arXiv

Herschel PEP: 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<z<2.5, normal star-forming galaxies (SFGs), which do not depend on extinction corrections and/or extrapolations of spectral energy distributions. The measurements are based on observations in GOODS-N with the Photodetector Array Camera & Spectrometer (PACS) onboard Herschel, as part of the PACS Evolutionary Probe (PEP) project, that resolve for the first time individual SFGs at these redshifts at FIR wavelengths. We compare FIR-based SFRs to the more commonly used 24 micron and UV SFRs. We find that SFRs from 24 micron alone are higher by a factor of ~4-7.5 than the true SFRs. This overestimation depends on luminosity: gradually increasing for log L(24um)>12.2 L_sun. The SFGs and AGNs tend to exhibit the same 24 micron 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.

preprint2010arXiv

Herschel unveils a puzzling uniformity of distant dusty galaxies

The Herschel Space Observatory enables us to accurately measure the bolometric output of starburst galaxies and active galactic nuclei (AGN) by directly sampling the peak of their far-infrared (IR) emission. Here we examine whether the spectral energy distribution (SED) and dust temperature of galaxies have strongly evolved since z~2.5. We use Herschel deep extragalactic surveys from 100 to 500um to compute total IR luminosities in galaxies down to the faintest levels, using PACS and SPIRE in the GOODS-North field (PEP and HerMES key programs). We show that measurements in the SPIRE bands can be used below the statistical confusion limit if information at higher spatial resolution is used to identify isolated galaxies whose flux is not boosted by bright neighbors. Below z~1.5, mid-IR extrapolations are correct for star-forming galaxies with a dispersion of only 40% (0.15dex), therefore similar to z~0 galaxies. This narrow distribution is puzzling when considering the range of physical processes that could have affected the SED of these galaxies. Extrapolations from only one of the 160um, 250um or 350um bands alone tend to overestimate the total IR luminosity. This may be explained by the lack of far-IR constraints around and above ~150um (rest-frame) on local templates. We also note that the dust temperature of luminous IR galaxies around z~1 is mildly colder by 10-15% than their local analogs and up to 20% for ULIRGs at z~1.6. Above z=1.5, distant galaxies are found to exhibit a substantially larger mid- over far-IR ratio, which could either result from stronger broad emission lines or warm dust continuum heated by a hidden AGN. Two thirds of the AGNs identified in the field with a measured redshift exhibit the same behavior as purely star-forming galaxies. Hence a large fraction of AGNs harbor star formation at very high SFR and in conditions similar to purely star-forming galaxies.

preprint2010arXiv

PEP: first Herschel probe of dusty galaxy evolution up to z~3

We exploit the deepest existing far-infrared (FIR) data obtained so far by Herschel at 100 and 160 um in the GOODS-N, as part of the PACS Evolutionary Probe (PEP) survey, to derive for the first time the evolution of the rest-frame 60-um, 90-um, and total IR luminosity functions (LFs) of galaxies and AGNs from z=0 to unprecedented high redshifts (z~2-3). The PEP LFs were computed using the 1/Vmax method. The FIR sources were classified by means of a detailed broad- band SED-fitting analysis and spectral characterisation. Based on the best-fit model results, k-correction and total IR (8-1000 um) luminosity were obtained for each source. LFs (monochromatic and total) were then derived for various IR populations separately in different redshift bins and compared to backward evolution model predictions. We detect strong evolution in the LF to at least z~2. Objects with SEDs similar to local spiral galaxies are the major contributors to the star formation density (SFD) at z< 0.3, then, as redshift increases, moderate SF galaxies - most likely containing a low-luminosity AGN - start dominating up to z ~= 1.5. At >1.5 the SFD is dominated by the contributions of starburst galaxies. In agreement with previous findings, the comoving IR LD derived from our data evolves approximately as (1 + z)^(3.8+/-0.3) up to z~1, there being some evidence of flattening up to z~2.

preprint2010arXiv

Star formation in AGN hosts in GOODS-N

Sensitive Herschel far-infrared observations can break degeneracies that were inherent to previous studies of star formation in high-z AGN hosts. Combining PACS 100 and 160um observations of the GOODS-N field with 2Msec Chandra data, we detect ~20% of X-ray AGN individually at >3sig. The host far-infrared luminosity of AGN with L2-10~10^43erg/s increases with redshift by an order of magnitude from z=0 to z~1. In contrast, there is little dependence of far-infrared luminosity on AGN luminosity, for L2-10<~10^44erg/s AGN at z>~1. We do not find a dependence of far-infrared luminosity on X-ray obscuring column, for our sample which is dominated by L2-10<10^44erg/s AGN. In conjunction with properties of local and luminous high-z AGN, we interpret these results as reflecting the interplay between two paths of AGN/host coevolution. A correlation of AGN luminosity and host star formation is traced locally over a wide range of luminosities and also extends to luminous high z AGN. This correlation reflects an evolutionary connection, likely via merging. For lower AGN luminosities, star formation is similar to that in non-active massive galaxies and shows little dependence on AGN luminosity. The level of this secular, non-merger driven star formation increasingly dominates over the correlation at increasing redshift.

preprint2010arXiv

The dust content of high-z submillimeter galaxies revealed by Herschel

We use deep observations taken with the Photodetector Array Camera and Spectrometer (PACS), on board the Herschel satellite as part of the PACS evolutionary probe (PEP) guaranteed project along with submm ground-based observations to measure the dust mass of a sample of high-z submillimeter galaxies (SMGs). We investigate their dust content relative to their stellar and gas masses, and compare them with local star-forming galaxies. High-z SMGs are dust rich, i.e. they have higher dust-to-stellar mass ratios compared to local spiral galaxies (by a factor of 30) and also compared to local ultraluminous infrared galaxies (ULIRGs, by a factor of 6). This indicates that the large masses of gas typically hosted in SMGs have already been highly enriched with metals and dust. Indeed, for those SMGs whose gas mass is measured, we infer dust-to-gas ratios similar or higher than local spirals and ULIRGs. However, similarly to other strongly star-forming galaxies in the local Universe and at high-z, SMGs are characterized by gas metalicities lower (by a factor of a few) than local spirals, as inferred from their optical nebular lines, which are generally ascribed to infall of metal-poor gas. This is in contrast with the large dust content inferred from the far-IR and submm data. In short, the metalicity inferred from the dust mass is much higher (by more than an order of magnitude) than that inferred from the optical nebular lines. We discuss the possible explanations of this discrepancy and the possible implications for the investigation of the metalicity evolution at high-z.

preprint2010arXiv

The far-infrared/radio correlation as probed by Herschel

We set out to determine the ratio, q(IR), of rest-frame 8-1000um flux, S(IR), to monochromatic radio flux, S(1.4GHz), for galaxies selected at far-IR and radio wavelengths, to search for signs that the ratio evolves with redshift, luminosity or dust temperature, and to identify any far-IR-bright outliers - useful laboratories for exploring why the far-IR/radio correlation is generally so tight when the prevailing theory suggests variations are almost inevitable. We use flux-limited 250-um and 1.4-GHz samples, obtained in GOODS-N using Herschel (HerMES; PEP) and the VLA. We determine bolometric IR output using ten bands spanning 24-1250um, exploiting data from PACS and SPIRE, as well as Spitzer, SCUBA, AzTEC and MAMBO. We also explore the properties of an L(IR)-matched sample, designed to reveal evolution of q(IR) with z, spanning log L(IR) = 11-12 L(sun) and z=0-2, by stacking into the radio and far-IR images. For 1.4-GHz-selected galaxies, we see tentative evidence of a break in the flux ratio, q(IR), at L(1.4GHz) ~ 10^22.7 W/Hz, where AGN are starting to dominate the radio power density, and of weaker correlations with z and T(d). From our 250-um-selected sample we identify a small number of far-IR-bright outliers, and see trends of q(IR) with L(1.4GHz), L(IR), T(d) and z, noting that some of these are inter-related. For our L(IR)-matched sample, there is no evidence that q(IR) changes significantly as we move back into the epoch of galaxy formation: we find q(IR) goes as (1+z)^gamma, where gamma = -0.04 +/- 0.03 at z=0-2; however, discounting the least reliable data at z<0.5 we find gamma = -0.26 +/- 0.07, modest evolution which may be related to the radio background seen by ARCADE2, perhaps driven by <10uJy radio activity amongst ordinary star-forming galaxies at z>1.

preprint2010arXiv

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

We exploit deep observations of the GOODS-N field taken with PACS, on board of Herschel, as part of the PEP guaranteed time, to study the link between star formation and stellar mass in galaxies to z~2. Starting from a stellar mass-selected sample of ~4500 galaxies with mag[4.5mu]<23 (AB), we identify ~350 objects with a PACS detection at 100 or 160mu and ~1500 with only Spitzer 24 mu counterpart. Stellar masses and total IR luminosities (LIR) are estimated by fitting the SEDs. Consistently with other Herschel results, we find that LIR based only on 24 mu data is overestimated by a median factor ~1.8 at z~2, whereas it is underestimated (with our approach) up to a factor ~1.6 at 0.5<z<1.0. We then exploit this calibration to correct LIR based on the MIPS fluxes. These results clearly show how Herschel is fundamental to constrain LIR, and hence the 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 star formation rate (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 alpha=-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 massive M>10^(11) (M_sun) galaxies from z=0 to z=2, and seems to flatten at z>1.5 in this mass range. Moreover, the most massive galaxies have the lowest SSFR at any z, implying that they have formed their stars earlier and more rapidly than their low mass counterparts.

preprint2010arXiv

The HELLAS2XMM survey. XIII. Multi-component analysis of the spectral energy distribution of obscured AGN

We combine near-to-mid-IR Spitzer data with shorter wavelength observations (optical to X-rays) to get insights on the properties of a sample of luminous, obscured Active Galactic Nuclei (AGN). We aim at modeling their broad-band Spectral Energy Distributions (SEDs) in order to estimate the main parameters related to the dusty torus. The sample comprises 16 obscured high-redshift (0.9<z<2.1) xray luminous quasars (L_2-10 ~ 10^44 erg s-1) selected from the HELLAS2XMM survey. The SEDs are described by a multi-component model including a stellar component, an AGN component and a starburst. The majority (~80%) of the sources show moderate optical depth (tau_9.7um<3) and the derived column densities N_H are consistent with the xray inferred values (10^22 <N_H< 3x10^23 cm-2) for most of the objects, confirming that the sources are moderately obscured Compton-thin AGN. Accretion luminosities in the range 5x10^44 < Lbol < 4x10^46 erg s-1 are inferred. We compare model luminosities with those obtained by integrating the observed SED, finding that the latter are lower by a factor of ~2 in the median. The discrepancy can be as high as an order of magnitude for models with high optical depth (tau_9.7um=10). The ratio between the luminosities obtained by the fitting procedure and from the observed SED suggest that, at least for Type~2 AGN, observed bolometric luminosities are likely to underestimate intrinsic ones and the effect is more severe for highly obscured sources. Bolometric corrections from the hard X-ray band are computed and have a median value of k_2-10kev ~ 20. The obscured AGN in our sample are characterized by relatively low Eddington ratios (median lambda_Edd~0.08). On average, they are consistent with the Eddington ratio increasing at increasing bolometric correction (e.g. Vasudevan & Fabiam 2009).

preprint2010arXiv

The M_BH-M_star relation of obscured AGNs at high redshift

We report the detection of broad Halpha emission in three X-ray selected obscured AGNs at z=1-2. By exploiting the Halpha width and the intrinsic X-ray luminosity, we estimate their black hole masses, which are in the range 0.1-3x10^9 Msun. By means of multi-band photometric data, we measure the stellar mass of their host galaxy and, therefore, infer their M_BH/M_star ratio. These are the first obscured AGNs at high-z, selected based on their black hole accretion (i.e. on the basis of their X-ray luminosity), that can be located on the M_BH-M_star relation at high-z. All of these obscured high-z AGNs are fully consistent with the local M_BH-M_star relation. This result conflicts with those for other samples of AGNs in the same redshift range, whose M_BH/M_star ratio departs significantly from the value observed in local galaxies. We suggest that the obscured AGNs in our sample are in an advanced evolutionary stage, have already settled onto the local M_BH-M_star relation, and whose nuclear activity has been temporarily revived by recent galaxy interactions.

preprint2009arXiv

Spectroscopic Cosmological Surveys in the Far-IR

We show the feasibility of spectroscopic cosmological surveys with the SAFARI instrument onboard of SPICA. The work is done through simulations that make use of both empirical methods, i.e. the use of observed luminosity functions and theoretical models for galaxy formation and evolution. The relations assumed between the line emission to trace AGN and star formation activity have been derived from the observations of local samples of galaxies. The results converge to indicate the use of blind spectroscopy with the SAFARI FTS at various resolutions to study galaxy evolution from the local to the distant (z~3) Universe. Specifically, two different and independent galaxy evolution models predict about 7-10 sources to be spectroscopically detected in more than one line in a 2'x 2'SAFARI field of view, down to the expected flux limits of SAFARI, with about 20% of sources to be detected at z>2. SPICA-SAFARI will be therefore excellent at detecting high-z sources and at assessing in a direct way their nature (e.g whether mainly AGN or Star Formation powered) thanks to blind spectroscopy.

preprint2009arXiv

Spectroscopic identifications of Spitzer sources in the SWIRE/XMM-NEWTON/ELAIS-S1 field: a large fraction of Active Galactic Nuclei with high F(24micron)/F(R) ratio

We present a catalog of optical spectroscopic identifications of sources detected by Spitzer at 3.6 or 24 micron down to 10 and 280 microJy, respectively, in the SWIRE/XMM-Newton/ELAIS-S1 field and classified via line width analysis and diagnostic diagrams. A total of 1376 sources down to R~24.2 mag have been identified (1362 detected at 3.6 micron, 419 at 24 micron, and 405 at both) by low-resolution optical spectroscopy carried out with FORS2, VIMOS, and EFOSC2 at the Very Large Telescope and 3.6m ESO telescopes. The spectroscopic campaigns have been carried out over the central 0.6 square degrees area of ELAIS-S1 which, in particular, has also been observed by XMM-Newton and Chandra. We find the first direct optical spectroscopic evidence that the fraction of active galactic nuclei (AGN; mostly AGN2) increases with increasing F(24 micron)/F(R) ratio, reaching values of 70(+/-20)% in the range 316<F(24 micron)F(R)<1000. We present an IRAC-MIPS color-color diagram able to separate AGN1 from obscured AGN2 candidates. After having corrected for the spectroscopic incompleteness of our sample, it results that the AGN fraction at F(24 micron)=0.8 mJy is ~22(+/-7)% and decreases slowly to ~19(+/-5)% down to F(24 micron)=0.3 mJy.

preprint2008arXiv

Direct observation of a localization transition in quasi-periodic photonic lattices

The localization of waves in non-periodic media is a universal phenomenon, occurring in a variety of different quantum and classical systems, including condensed-matter, Bose-Einstein condensates in optical lattices, quantum chaotic systems, sound waves and light. A localization phase transition is expected to occur in three dimensional disordered systems as the strength of disorder crosses a critical value. Recently, a crossover from an extended to a localized phase has been observed in low-dimensional photonic lattices and Bose-Einstein condensates. Other experiments studied the critical behaviour near the transition in three dimensions via transmission measurements. However, no direct observation of a localization transition for light has been reported. In 1979 Aubry and Andre predicted that for a certain class of quasi-periodic potentials, a localization phase transition can occur already in one-dimension. This strongly correlated potential is markedly different from the disordered case were the potential is uncorrelated. Here we report an experiment that realizes the Aubry-Andre model in quasi-periodic photonic lattices. We observe the signature of a localization phase transition by directly measuring the expansion rates of initially narrow wave packets propagating in the lattice. Below the transition point, all the modes of the system are extended and therefore an initially narrow wave-packet eventually spreads across the entire lattice. Above the critical point, all modes are localized and expansion is suppressed. In addition, we study the effect of weak nonlinear interactions on light propagation below and above the transition.

preprint2008arXiv

The Contribution of AGN and Star-Forming Galaxies to the Mid-Infrared as Revealed by their Spectral Energy Distributions

We present the broad-band Spectral Energy Distributions (SEDs) of the largest available highly (72%) complete spectroscopic sample of mid-infrared (MIR) selected galaxies and AGN at intermediate redshift. The sample contains 203 extragalactic sources from the 15-micron survey in the ELAIS-SWIRE field S1, all with measured spectroscopic redshift. Most of these sources have full multi-wavelength coverage from the far-UV to the far-infrared and lie in the redshift range 0.1<z<1.3. Due to its size, this sample allows us for the first time to characterise the spectral properties of the sources responsible for the strong evolution observed in the MIR. Based on SED-fitting technique we have classified the MIR sources, identifying AGN signatures in about 50% of them. This fraction is significantly higher than that derived from optical spectroscopy (~29%) and is due in particular to the identification of AGN activity in objects spectroscopically classified as galaxies. It is likely that in most of our objects, the AGN is either obscured or of low-luminosity, and thus it does not dominate the energetic output at any wavelength, except in the MIR, showing up just in the range where the host galaxy SED has a minimum. The fraction of AGN strongly depends on the flux density, with that derived through the SED-fitting being about 20% at S(15)~0.5-1 mJy and gradually increasing up to 100% at S(15)>10 mJy, while that obtained from optical spectroscopy never being >30%, even at the higher flux densities. The results of this work will be very useful for updating all the models aimed at interpreting the deep infrared survey data and, in particular, for constraining the nature and the role of dust-obscured systems in the intermediate/high-redshift Universe.

preprint2004arXiv

The European Large Area ISO Survey (ELAIS): The Final Band-merged Catalogue

We present the final band-merged ELAIS catalogue at 6.7, 15, 90, and 175 $μ$m, and the associated data at u,g,r,i,z,J,H,K, and 20cm. The origin of the survey, infrared and radio observations, data-reduction and optical identifications are briefly reviewed, and a summary of the area covered, and completeness limit for each infrared band is given. A detailed discussion of the band-merging and optical association strategy is given. The total catalogues consists of 2860 sources. For extragalactic sources observed in 3 or more infrared bands, colour-colour diagrams are presented and discussed in terms of the contributing infrared populations. Spectral energy distributions are shown for selected sources and compared with cirrus, M82 and Arp220 starburst, and AGN dust torus models. Spectroscopic redshifts are tabulated, where available. For the N1 and N2 areas, the INT ugriz Wide Field Survey, permits photometric redshifts to be estimated for galaxies and quasars. These agree well with the spectroscopic redshifts, within the uncertainty of the photometric method. The redshift distribution is given for selected ELAIS bands and colour-redshift diagrams are discussed. There is a high proportion of ultraluminous infrared galaxies in the ELAIS Catalogue ($> 10 %$ of 15 $μ$m sources), many with Arp220-like colours. 10 hyperluminous infrared galaxies and 10 EROs are found in the survey. The large numbers of ultraluminous galaxies imply very strong evolution in the star-formation rate between z = 0 and 1.