Source author record

Emanuele Daddi

Emanuele Daddi 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

71works
6topics
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

71 published item(s)

preprint2026arXiv

Galactic bars are already mature at Cosmic Noon: bar strength and flatness at z ~ 1.5

In this work, we explore the nature of $z>1$ galactic bars. Once thought to be highly transient, our results demonstrate otherwise. Our sample consists of nine massive ($>10^{10.5}\,\rm M_{\odot}$) star-forming barred-spiral galaxies at $z_{\rm spec} \sim 1.5$. Using rest-frame near-IR (F444W) JWST/NIRCam imaging, we apply ellipse fitting along with 1D and 2D morphological modeling to directly measure bar properties. We find that five galaxies host flat surface brightness profiles (bar Sérsic index $<0.4$), indicative of highly evolved, "mature" bars. By contrast, only two galaxies show exponential profiles, characteristic of young bars, and these are also shorter in absolute length than the flat bars. We therefore conclude that a large fraction of bars at this epoch have already matured, thereby indicating the presence of well-settled disks required to facilitate bar formation and sustained evolution well before $z\sim1.5$. To assess the gravitational impact of the bars, we calculate the maximum transverse-to-radial force ratio ($Q_{b}$). We find that $Q_{b}$ values are comparable to, or weaker than, those of bars in the local Universe, Seven of the nine bars show only a marginal increase in strength with maturity (from exponential to flat bars). Contrarily however, the remaining two bars are flat, but have the lowest $Q_{b}$ values in our sample. We hence propose that the mature bars at $z\sim 1.5$ may experience phases of weakening due to rapid gas inflows and/or minor mergers. In conclusion, our work sheds light on the rapidly evolving nature of high-z bars and paves the way for larger statistical studies.

preprint2026arXiv

Galaxies in the Epoch of Reionization Are All Bark and No Bite -- Plenty of Ionizing Photons, Low Escape Fractions

Early results from JWST suggest that epoch-of-reionization (EoR) galaxies produce copious ionizing photons, which, if they escape efficiently, could cause reionization to occur too early. We study this problem using \jwst\ imaging and prism spectroscopy for 412 galaxies at 4.5 < z < 9.0. We fit these data simultaneously with stellar-population and nebular-emission models that include a parameter for the fraction of ionizing photons that escape the galaxy, $f_\mathrm{esc}$. We find that the ionization production efficiency, $ξ_\mathrm{ion}$ = Q(H) / L(UV), increases with redshift and decreasing UV luminosity, but shows significant scatter, $σ( \log ξ_\mathrm{ion})$ = 0.3 dex. The inferred escape fractions averaged over the population are low, ranging from $\langle f_\mathrm{esc} \rangle$ = $2.6\pm 1.4$\% at 6 < z < 9 to $6.5\pm 2.2$\% at 4.5 < z < 6 with weak or no indication of evolution with redshift. This implies that in our models most of the ionizing photons need to be absorbed to account for the nebular emission. We compute the impact of our results on reionization, including the distributions for $ξ_\mathrm{ion}$ and $f_\mathrm{esc}$, and the evolution and uncertainty of the UV luminosity function. Considering galaxies brighter than M(UV) < -16 mag, we would produce an IGM hydrogen-ionized fraction of $x_e = 0.5$ at 5.3 < z < 5.8, possibly too late compared to constraints from from QSO sightlines. Including fainter galaxies, M(UV) < -14 mag, we obtain $x_e = 0.5$ at 6.0 < z < 8.1, fully consistent with QSO and CMB data. This implies that EoR galaxies produce plenty of ionizing photons, but these do not efficiently escape. This may be a result of high gas column densities combined with burstier star-formation histories, which limit the time massive stars are able to clear channels through the gas for ionizing photons to escape.

preprint2022arXiv

AGN Selection Methods Have Profound Impacts on the Distributions of Host Galaxy Properties

We present a comparative study of X-ray and IR AGNs at $z\approx2$ to highlight the important AGN selection effects on the distributions of host galaxy properties. Compared with non-AGN star-forming galaxies (SFGs) on the main sequence, X-ray AGNs have similar median star formation (SF) properties, but their incidence (q$_{\rm{AGN}}$) is higher among galaxies with either enhanced or suppressed SF, and among galaxies with larger stellar mass surface density, regardless if it is measured within half-light radius ($Σ_e$) or central 1kpc ($Σ_{\rm{1kpc}}$). Unlike X-ray AGNs, IR AGNs are less massive, and have enhanced SF and similar distributions of colors, $Σ_e$ and $Σ_{\rm{1kpc}}$ relative to non-AGN SFGs. Given that $Σ_e$ and $Σ_{\rm{1kpc}}$ strongly correlate with M$_*$, we introduce the fractional mass within central 1kpc ($\rm{\frac{M_{1kpc}}{M_*}}$), which only weakly depends on M$_*$, to quantify galaxy compactness. Both AGN populations have similar $\rm{\frac{M_{1kpc}}{M_*}}$ distributions compared to non-AGN SFGs'. While q$_{\rm{AGN}}$ increases with $Σ_e$ and $Σ_{\rm{1kpc}}$, it remains constant with $\rm{\frac{M_{1kpc}}{M_*}}$, indicating that the trend of increasing q$_{\rm{AGN}}$ with $\rmΣ$ is driven by M$_*$ more than morphology. While our findings are not in conflict with the scenario of AGN quenching, they do not imply it either, because the incidence of AGNs hosted in transitional galaxies depends crucially on AGN selections. Additionally, despite the relatively large uncertainty of AGN bolometric luminosities, their very weak correlation, if any, with SF activities, regardless of AGN selections, also argues against a direct causal link between the presences of AGNs and the quenching of massive galaxies at $z\sim2$.

preprint2022arXiv

Coincidence between morphology and star-formation activity through cosmic time: the impact of the bulge growth

The origin of the quenching in galaxies is still highly debated. Different scenarios and processes are proposed. We use multi-band (400-1600 nm) bulge-disc decompositions of massive galaxies in the redshift range 0<z<2 to explore the distribution and the evolution of galaxies in the log SFR-logM* plane as a function of the stellar mass-weighted bulge-to-total ratio (BTM) and also for internal galaxy components (bulge/disc) separately. We find evidence of a clear link between the presence of a bulge and the flattening of the Main Sequence in the high-mass end. All bulgeless galaxies (BTM<0.2) lie on the main-sequence, and there is little evidence of a quenching channel without bulge growth. Galaxies with a significant bulge component (BTM>0.2) are equally distributed in number between star forming and passive regions. The vast majority of bulges in the Main Sequence galaxies are quiescent, while star formation is localized in the disc component. Our current findings underline a strong correlation between the presence of the bulge and the star formation state of the galaxy. A bulge, if present, is often quiescent, independently of the morphology or the star formation activity of the host galaxy. Additionally, if a galaxy is quiescent, with a large probability, is hosting a bulge. Conversely, if the galaxy has a disky shape is highly probable to be star forming.

preprint2022arXiv

Starbursts with suppressed velocity dispersion revealed in a forming cluster at z=2.51

One of the most prominent features of galaxy clusters is the presence of a dominant population of massive ellipticals in their cores. Stellar archaeology suggests that these gigantic beasts assembled most of their stars in the early Universe via starbursts. However, the role of dense environments and their detailed physical mechanisms in triggering starburst activities remain unknown. Here we report spatially resolved Atacama Large Millimeter/submillimeter Array (ALMA) observations of the CO $J= 3-2$ emission line, with a resolution of about 2.5 kiloparsecs, toward a forming galaxy cluster core with starburst galaxies at $z=2.51$. In contrast to starburst galaxies in the field often associated with galaxy mergers or highly turbulent gaseous disks, our observations show that the two starbursts in the cluster exhibit dynamically cold (rotation-dominated) gas-rich disks. Their gas disks have extremely low velocity dispersion ($σ_{\mathrm{0}} \sim 20-30$ km s$^{-1}$), which is three times lower than their field counterparts at similar redshifts. The high gas fraction and suppressed velocity dispersion yield gravitationally unstable gas disks, which enables highly efficient star formation. The suppressed velocity dispersion, likely induced by the accretion of corotating and coplanar cold gas, might serve as an essential avenue to trigger starbursts in massive halos at high redshifts.

preprint2021arXiv

A titanic interstellar medium ejection from a massive starburst galaxy at redshift 1.4

Feedback-driven winds from star formation or active galactic nuclei might be a relevant channel for the abrupt quenching star formation in massive galaxies. However, both observations and simulations support the idea that these processes are non-conflictingly co-evolving and self-regulating. Furthermore, evidence of disruptive events that are capable of fast quenching is rare, and constraints on their statistical prevalence are lacking. Here we present a massive starburst galaxy at z=1.4 which is ejecting $46 \pm 13$\% of its molecular gas mass at a startling rate of $\gtrsim 10,000$ M$_{\odot}{\rm yr}^{-1}$. A broad component that is red-shifted from the galaxy emission is detected in four (low- and high-J) CO and [CI] transitions and in the ionized phase, which ensures a robust estimate of the expelled gas mass. The implied statistics suggest that similar events are potentially a major star-formation quenching channel. However, our observations provide compelling evidence that this is not a feedback-driven wind, but rather material from a merger that has been probably tidally ejected. This finding challenges some literature studies in which the role of feedback-driven winds might be overstated.

preprint2021arXiv

CO excitation, molecular gas density and interstellar radiation field in local and high-redshift galaxies

We study the Carbon Monoxide (CO) excitation, mean molecular gas density and interstellar radiation field (ISRF) intensity in a comprehensive sample of 76 galaxies from local to high redshift (z~0-6), selected based on detections of their CO transitions J=2-1 and 5-4 and their optical/infrared/(sub-)millimeter spectral energy distributions (SEDs). We confirm the existence of a tight correlation between CO excitation as traced by the CO(5-4)/(2-1) line ratio (R52), and the mean ISRF intensity U as derived from infrared SED fitting using dust SED templates. By modeling the molecular gas density probability distribution function (PDF) in galaxies and predicting CO line ratios with large velocity gradient radiative transfer calculations, we present a framework linking global CO line ratios to the mean molecular hydrogen gas density nH2 and kinetic temperature Tkin. Mapping in this way observed R52 ratios to nH2 and Tkin probability distributions, we obtain positive U-nH2 and U-Tkin correlations, which imply a scenario in which the ISRF in galaxies is mainly regulated by Tkin and (non-linearly) by nH2. A small fraction of starburst galaxies showing enhanced nH2 could be due to merger-driven compaction. Our work demonstrates that ISRF and CO excitation are tightly coupled, and that density-PDF modeling is a promising tool for probing detailed ISM properties inside galaxies.

preprint2021arXiv

COLDz: Probing Cosmic Star Formation With Radio Free-free Emission

Radio free-free emission is considered to be one of the most reliable tracers of star formation in galaxies. However, as it constitutes the faintest part of the radio spectrum -- being roughly an order of magnitude less luminous than radio synchrotron emission at the GHz frequencies typically targeted in radio surveys -- the usage of free-free emission as a star formation rate tracer has mostly remained limited to the local Universe. Here we perform a multi-frequency radio stacking analysis using deep Karl G. Jansky Very Large Array observations at 1.4, 3, 5, 10 and 34 GHz in the COSMOS and GOODS-North fields to probe free-free emission in typical galaxies at the peak of cosmic star formation. We find that $z \sim 0.5 - 3$ star-forming galaxies exhibit radio emission at rest-frame frequencies of $\sim 65 - 90$ GHz that is $\sim 1.5 - 2\times$ fainter than would be expected from a simple combination of free-free and synchrotron emission, as in the prototypical starburst galaxy M82. We interpret this as a deficit in high-frequency synchrotron emission, while the level of free-free emission is as expected from M82. We additionally provide the first constraints on the cosmic star formation history using free-free emission at $0.5 \lesssim z \lesssim 3$, which are in good agreement with more established tracers at high redshift. In the future, deep multi-frequency radio surveys will be crucial in order to accurately determine the shape of the radio spectrum of faint star-forming galaxies, and to further establish radio free-free emission as a tracer of high-redshift star formation.

preprint2021arXiv

Feedback Factory: Multiple faint radio-jets detected in a cluster at z=2

We report the detection of multiple faint radio sources, that we identify as AGN-jets, within CLJ1449+0856 at z=2 using 3 GHz VLA observations. We study the effects of radio-jet based kinetic feedback at high redshifts, which has been found to be crucial in low redshift clusters to explain the observed thermodynamic properties of their ICM. We investigate this interaction at an epoch featuring high levels of AGN activity and a transitional phase of ICM in regards to the likelihood of residual cold-gas accretion. We measure a total flux of $\rm 30.6 \pm 3.3~μJy$ from the 6 detected jets. Their power contribution is estimated to be $1.2 ~(\pm 0.6)~ \times 10^{44} ~\rm ergs~ s^{-1}$, although this value could be up to $4.7 ~ \times 10^{44} ~\rm ergs~ s^{-1}$. This is a factor $\sim 0.25 - 1.0$ of the previously estimated instantaneous energy injection into the ICM of CLJ1449+0856 from AGN outflows and star formation, that have already been found to be sufficient in globally offsetting the cooling flows in the cluster core. In line with the already detected abundance of star formation, this mode of feedback being distributed over multiple sites, contrary to a single central source observed at low redshifts, points to accretion of gas into the cluster centre. This also suggests a 'steady state' of the cluster featuring non cool-core like behaviour. Finally, we also examine the TIR-radio luminosity ratio for the known sample of galaxies within the cluster core and find that dense environments do not have any serious consequence on the compliance of galaxies to the IR-radio correlation.

preprint2021arXiv

Hierarchical fragmentation in high redshift galaxies revealed by hydrodynamical simulations

High-redshift star-forming galaxies have very different morphologies compared to nearby ones. Indeed, they are often dominated by bright star-forming structures of masses up to $10^{8-9}$ $\mathrm{M}_\odot$ dubbed «giant clumps». However, recent observations questioned this result by showing only low-mass structures or no structure at all. We use Adaptative Mesh Refinement hydrodynamical simulations of galaxies with parsec-scale resolution to study the formation of structures inside clumpy high-redshift galaxies. We show that in very gas-rich galaxies star formation occurs in small gas clusters with masses below $10^{7-8}$ $\mathrm{M}_\odot$ that are themselves located inside giant complexes with masses up to $10^8$ and sometimes $10^9$ $\mathrm{M}_\odot$. Those massive structures are similar in mass and size to the giant clumps observed in imaging surveys, in particular with the Hubble Space Telescope. Using mock observations of simulated galaxies, we show that at very high resolution with instruments like the Atacama Large Millimeter Array or through gravitational lensing, only low-mass structures are likely to be detected, and their gathering into giant complexes might be missed. This leads to the non-detection of the giant clumps and therefore introduces a bias in the detection of these structures. We show that the simulated giant clumps can be gravitationally bound even when undetected in mocks representative for ALMA observations and HST observations of lensed galaxies. We then compare the top-down fragmentation of an initially warm disc and the bottom-up fragmentation of an initially cold disc to show that the process of formation of the clumps does not impact their physical properties.

preprint2021arXiv

The stellar mass versus stellar metallicity relation of star-forming galaxies at $1.6\le z\le3.0$ and implications for the evolution of the $α$-enhancement

We measure the relationship between stellar mass and stellar metallicity, the stellar mass--metallicity relation (MZR), for 1336 star-forming galaxies at $1.6\le z\le3.0$ (<z>=2.2) using rest-frame far-ultraviolet spectra from the zCOSMOS-deep survey. High signal-to-noise composite spectra containing stellar absorption features are fit with population synthesis model spectra of a range of metallicity. We find stellar metallicities, which mostly reflect iron abundances, scaling as $(Z_{Fe,\ast}/Z_{Fe,\odot})=-(0.81\pm0.01)+(0.32+0.03)\log(M_\ast/10^{10}M_\odot)$ across the mass range of $10^9\lesssim M_\ast/M_\odot\lesssim10^{11}$, being $\approx6\times$ lower than seen locally at the same masses. The instantaneous oxygen-to-iron ratio ($α$-enhancement) inferred using the gas-phase oxygen MZRs, is on average found to be [O/Fe]$\approx0.47$, being higher than the local [O/Fe]$\approx0$. The observed changes in [O/Fe] and [Fe/H] are reproduced in simple flow-through gas-regulator models with steady star-formation histories (SFHs) that follow the evolving main sequence. Our models show that the [O/Fe] is determined almost entirely by the instantaneous specific star formation rate alone while being independent of the SFHs, mass, and the gas-regulation characteristics of the systems. We find that the locations of $\sim10^{10}M_\odot$ galaxies at z~2 in the [O/Fe]--metallicity planes are in remarkable agreement with the sequence of low-metallicity thick-disk stars in our Galaxy. This manifests a beautiful concordance between the results of Galactic archaeology and observations of high-redshift Milky Way progenitors. However, there remains a question of how and when the old metal-rich, low-$α$/Fe stars seen in the bulge had formed by z~2 because such a stellar population is not seen in our data and difficult to explain in the context of our models.

preprint2020arXiv

A Census of Sub-kiloparsec Resolution Metallicity Gradients in Star-forming Galaxies at Cosmic Noon from HST Slitless Spectroscopy

We present hitherto the largest sample of gas-phase metallicity radial gradients measured at sub-kiloparsec resolution in star-forming galaxies in the redshift range of $z\in[1.2, 2.3]$. These measurements are enabled by the synergy of slitless spectroscopy from the Hubble Space Telescope near-infrared channels and the lensing magnification from foreground galaxy clusters. Our sample consists of 76 galaxies with stellar mass ranging from 10$^7$ to 10$^{10}$ $M_\odot$, instantaneous star-formation rate in the range of [1, 100] $M_\odot$/yr, and global metallicity [$\frac{1}{12}$, 2] solar. At 2-$σ$ confidence level, 15/76 galaxies in our sample show negative radial gradients, whereas 7/76 show inverted gradients. Combining ours and all other metallicity gradients obtained at similar resolution currently available in the literature, we measure a negative mass dependence of $Δ\log({\rm O/H})/Δr~ [\mathrm{dex~kpc^{-1}}] = \left(-0.020\pm0.007\right) + \left(-0.016\pm0.008\right) \log(M_\ast/10^{9.4} M_\odot)$ with the intrinsic scatter being $σ=0.060\pm0.006$ over four orders of magnitude in stellar mass. Our result is consistent with strong feedback, not secular processes, being the primary governor of the chemo-structural evolution of star-forming galaxies during the disk mass assembly at cosmic noon. We also find that the intrinsic scatter of metallicity gradients increases with decreasing stellar mass and increasing specific star-formation rate. This increase in the intrinsic scatter is likely caused by the combined effect of cold-mode gas accretion and merger-induced starbursts, with the latter more predominant in the dwarf mass regime of $M_\ast\lesssim10^9 M_\odot$.

preprint2020arXiv

COLDz: A High Space Density of Massive Dusty Starburst Galaxies ~1 Billion Years after the Big Bang

We report the detection of CO($J$=2$\to$1) emission from three massive dusty starburst galaxies at $z$$>$5 through molecular line scans in the NSF's Karl G. Jansky Very Large Array (VLA) CO Luminosity Density at High Redshift (COLDz) survey. Redshifts for two of the sources, HDF 850.1 ($z$=5.183) and AzTEC-3 ($z$=5.298), were previously known. We revise a previous redshift estimate for the third source GN10 ($z$=5.303), which we have independently confirmed through detections of CO $J$=1$\to$0, 5$\to$4, 6$\to$5, and [CII] 158 $μ$m emission with the VLA and the NOrthern Extended Milllimeter Array (NOEMA). We find that two currently independently confirmed CO sources in COLDz are "optically dark", and that three of them are dust-obscured galaxies at $z$$>$5. Given our survey area of $\sim$60 arcmin$^2$, our results appear to imply a $\sim$6-55 times higher space density of such distant dusty systems within the first billion years after the Big Bang than previously thought. At least two of these $z$$>$5 galaxies show star-formation rate surface densities consistent with so-called "maximum" starbursts, but we find significant differences in CO excitation between them. This result may suggest that different fractions of the massive gas reservoirs are located in the dense, star-forming nuclear regions - consistent with the more extended sizes of the [CII] emission compared to the dust continuum and higher [CII]-to-far-infrared luminosity ratios in those galaxies with lower gas excitation. We thus find substantial variations in the conditions for star formation between $z$$>$5 dusty starbursts, which typically have dust temperatures $\sim$57%$\pm$25% warmer than starbursts at $z$=2-3 due to their enhanced star formation activity.

preprint2020arXiv

Deceptively cold dust in the massive starburst galaxy GN20 at $z\sim4$

We present new observations, carried out with IRAM NOEMA, of the atomic neutral carbon transitions [CI](1-0) at 492 GHz and [CI](2-1) at 809 GHz of GN20, a well-studied star-bursting galaxy at $z=4.05$. The high luminosity line ratio [CI](2-1)/[CI](1-0) implies an excitation temperature of $48^{+14}_{-9}$ K, which is significantly higher than the apparent dust temperature of $T_{\rm d}=33\pm2$ K ($β=1.9$) derived under the common assumption of an optically thin far-infrared dust emission, but fully consistent with $T_{\rm d}=52\pm5$ K of a general opacity model where the optical depth ($τ$) reaches unity at a wavelength of $λ_0=170\pm23$ $μ$m. Moreover, the general opacity solution returns a factor of $\sim 2\times$ lower dust mass and, hence, a lower molecular gas mass for a fixed gas-to-dust ratio, than with the optically thin dust model. The derived properties of GN20 thus provide an appealing solution to the puzzling discovery of starbursts appearing colder than main-sequence galaxies above $z>2.5$, in addition to a lower dust-to-stellar mass ratio that approaches the physical value predicted for starburst galaxies.

preprint2020arXiv

Dust, gas, and metal content in star-forming galaxies at $z\sim3.3$ revealed with ALMA and Near-IR spectroscopy

We conducted sub-millimeter observations with the Atacama Large Millimeter/sub-millimeter Array (ALMA) of star-forming galaxies at $z\sim3.3$, whose gas-phase metallicities have been previously measured. We investigate the dust and gas contents of the galaxies at $z\sim3.3$ and study how galaxies are interacting with their circumgalactic/intergalactic medium at this epoch by probing their gas mass fractions and gas-phase metallicities. Single-band dust continuum emission tracing dust mass and the relation between the gas-phase metallicity and gas-to-dust mass ratio are used to estimate the gas masses. The estimated gas mass fractions and depletion timescales are $f_{\rm gas}=$ 0.20-0.75 and $t_{\rm dep}=$ 0.09-1.55 Gyr, respectively. Although the galaxies appear to tightly distribute around the star-forming main sequence at $z\sim3.3$, both quantities show a wider spread at a fixed stellar mass than expected from the scaling relation, suggesting a large diversity of fundamental gas properties among star-forming galaxies apparently on the main sequence. Comparing gas mass fraction and gas-phase metallicity between the star-forming galaxies at $z\sim3.3$ and at lower redshifts, star-forming galaxies at $z\sim3.3$ appear to be more metal-poor than local galaxies with similar gas mass fractions. Using the gas regulator model to interpret this offset, we find that it can be explained by a higher mass-loading factor, suggesting that the mass-loading factor in outflows increases at earlier cosmic times.

preprint2020arXiv

Early forecasts of the evolution of the COVID-19 outbreaks and quantitative assessment of the effectiveness of countering measures

We discovered that the time evolution of the inverse fractional daily growth of new infections, N/dN, in the current outbreak of COVID-19 is accurately described by a universal function, namely the two-parameter Gumbel cumulative function, in all countries that we have investigated. While the two Gumbel parameters, as determined bit fits to the data, vary from country to country (and even within different regions of the same country), reflecting the diversity and efficacy of the adopted containment measures, the functional form of the evolution of N/dN appears to be universal. The result of the fit in a given region or country appears to be stable against variations of the selected time interval. This makes it possible to robustly estimate the two parameters from the data data even over relatively small time periods. In turn, this allows one to predict with large advance and well-controlled confidence levels, the time of the peak in the daily new infections, its magnitude and duration (hence the total infections), as well as the time when the daily new infections decrease to a pre-set value (e.g. less than about 2 new infections per day per million people), which can be very useful for planning the reopening of economic and social activities. We use this formalism to predict and compare these key features of the evolution of the COVID-19 disease in a number of countries and provide a quantitative assessment of the degree of success in in their efforts to countain the outbreak.

preprint2020arXiv

From haloes to galaxies -- II. The fundamental relations in star formation and quenching

Star formation and quenching are two of the most important processes in galaxy formation and evolution. We explore in the local Universe the interrelationships among key integrated galaxy properties, including stellar mass $M_*$, star formation rate (SFR), specific SFR (sSFR), molecular gas mass $M_{\rm H_2}$, star formation efficiency (SFE) of the molecular gas and molecular gas to stellar mass ratio $μ$. We aim to identify the most fundamental scaling relations among these key galaxy properties and their interrelationships. We show the integrated $M_{\rm H_2}$-SFR, SFR-$M_*$ and $M_{\rm H_2}$-$M_*$ relation can be simply transformed from the $μ$-sSFR, SFE-$μ$ and SFE-sSFR relation, respectively. The transformation, in principle, can increase or decrease the scatter of each relation. Interestingly, we find the latter three relations all have significantly smaller scatter than the former three corresponding relations. We show the probability to achieve the observed small scatter by accident is extremely close to zero. This suggests that the smaller scatters of the latter three relations are driven by a more fundamental physical connection among these quantities. We then show the large scatters in the former relations are due to their systematic dependence on other galaxy properties, and on star formation and quenching process. We propose the sSFR-$μ$-SFE relation as the Fundamental Formation Relation (FFR), which governs the star formation and quenching process, and provides a simple framework to study galaxy evolution. Other scaling relations, including integrated Kennicutt-Schmidt law, star-forming main sequence and molecular gas main sequence, can all be derived from the FFR.

preprint2020arXiv

The ALMA Spectroscopic Survey in the HUDF: Deep 1.2 mm continuum number counts

We present the results from the 1.2 mm continuum image obtained as part of the ALMA Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS). The 1.2 mm continuum image has a size of 2.9 (4.2) arcmin$^2$ within a primary beam response of 50% (10%) and a rms value of $9.3\thinspace{\rmμJy\thinspace beam^{-1}}$. We detect 35 sources at high significance (Fidelity $\geq0.5$), 32 of these have well characterized near-infrared HST counterparts. We estimate the 1.2 mm number counts to flux levels of $<30\thinspace{\rmμJy}$ in two different ways: we first use the detected sources to constrain the number counts and find a significant flattening of the counts below $S_ν\sim 0.1$ mJy. In a second approach, we constrain the number counts by using a probability of deflection statistics (P(D)) analysis. For this latter approach, we describe new methods to accurately measure the noise in interferometric imaging (employing jack-knifing in the cube and in the visibility plane). This independent measurement confirms the flattening of the number counts. Our analysis of the differential number counts shows that we are detecting $\sim$93% ($\sim$100% if we include the lower fidelity detections) of the total continuum dust emission associated to galaxies in the HUDF. The ancillary data allows us to study the dependence of the 1.2 mm number counts on redshift ($z=0-4$), galaxy dust mass (${\rm M}_{\rm dust}=10^{7}-10^{9}{\rm M}_{\odot}$), stellar mass (${\rm M}_{*}=10^{9}-10^{12}{\rm M}_{\odot}$), and star-formation rate (${\rm SFR}=1-1000\thinspace{\rm M}_{\odot}\thinspace{\rm yr^{-1}}$). In an accompanying paper we show that the number counts are crucial to constrain galaxy evolution models and the understanding of star-forming galaxies at high redshift.

preprint2020arXiv

The ALMA Spectroscopic Survey in the HUDF: The Cosmic Dust and Gas Mass Densities in Galaxies up to $z\sim3$

Using the deepest 1.2 mm continuum map to date in the Hubble Ultra Deep Field obtained as part of the ALMA Spectroscopic Survey (ASPECS) large program, we measure the cosmic density of dust and implied gas (H$_{2}+$H I) mass in galaxies as a function of look-back time. We do so by stacking the contribution from all $H$-band selected galaxies above a given stellar mass in distinct redshift bins, $ρ_{\rm dust}(M_\ast>M,z)$ and $ρ_{\rm gas}(M_\ast>M,z)$. At all redshifts, $ρ_{\rm dust}(M_\ast>M,z)$ and $ρ_{\rm gas}(M_\ast>M,z)$ grow rapidly as $M$ decreases down to $10^{10}\,M_\odot$, but this growth slows down towards lower stellar masses. This flattening implies that at our stellar mass-completeness limits ($10^8\,M_\odot$ and $10^{8.9}\,M_\odot$ at $z\sim0.4$ and $z\sim3$), both quantities converge towards the total cosmic dust and gas mass densities in galaxies. The cosmic dust and gas mass densities increase at early cosmic time, peak around $z\sim2$, and decrease by a factor $\sim4$ and 7, compared to the density of dust and molecular gas in the local universe, respectively. The contribution of quiescent galaxies -- i.e., with little on-going star-formation -- to the cosmic dust and gas mass densities is minor ($\lesssim10\%$). The redshift evolution of the cosmic gas mass density resembles that of the star-formation rate density, as previously found by CO-based measurements. This confirms that galaxies have relatively constant star-formation efficiencies (within a factor $\sim2$) across cosmic time. Our results also imply that by $z\sim0$, a large fraction ($\sim90\%$) of dust formed in galaxies across cosmic time has been destroyed or ejected to the intergalactic medium.

preprint2020arXiv

The properties of the interstellar medium of galaxies across time as traced by the neutral atomic carbon [CI]

We report ALMA observations of the neutral atomic carbon transitions [CI] and multiple CO lines in a sample of $\sim30$ main sequence galaxies at $z\sim1$, including novel information on [CI](2-1) and CO(7-6) for 7 of such normal objects. We complement our observations with a collection of $>200$ galaxies with coverage of similar transitions, spanning the $z=0-4$ redshift interval and a variety of ambient conditions from local to high-redshift starbursts. We find systematic variations in the [CI]/IR and [CI]/high-$J$ ($J=7$) CO luminosity ratios among the various samples. We interpret these differences as increased dense molecular gas fractions and star formation efficiencies in the strongest high-redshift starbursts with respect to normal main sequence galaxies. We further report constant $L'_{\rm [CI]2-1}$/$L'_{\rm [CI]1-0}$ ratios across the galaxy populations and redshifts, suggesting that gas temperatures $T_{\rm exc}$ traced by [CI] do not strongly vary. We find only a mild correlation with $T_{\rm dust}$ and that, generally, $T_{\rm exc} \lesssim T_{\rm dust}$. We fit the line ratios with classical PDR models, retrieving consistently larger densities and intensities of the UV radiation fields in submm galaxies than in main sequence and local objects. However, these simple models fall short in representing the complexity of a multiphase interstellar medium and should be treated with caution. Finally, we compare our observations with the Santa Cruz semi-analytical model of galaxy evolution, recently extended to simulate submm emission. While we confirm the success in reproducing the CO lines, we find systematically larger [CI] luminosities at fixed IR luminosity than predicted theoretically. This highlights the necessity of improving our understanding of the mechanisms regulating the [CI] emission on galactic scales. We release our data compilation to the community.

preprint2019arXiv

Probing black hole accretion tracks, scaling relations and radiative efficiencies from stacked X-ray active galactic nuclei

The masses of supermassive black holes at the centres of local galaxies appear to be tightly correlated with the mass and velocity dispersions of their galactic hosts. However, the local Mbh-Mstar relation inferred from dynamically measured inactive black holes is up to an order-of-magnitude higher than some estimates from active black holes, and recent work suggests that this discrepancy arises from selection bias on the sample of dynamical black hole mass measurements. In this work we combine X-ray measurements of the mean black hole accretion luminosity as a function of stellar mass and redshift with empirical models of galaxy stellar mass growth, integrating over time to predict the evolving Mbh-Mstar relation. The implied relation is nearly independent of redshift, indicating that stellar and black hole masses grow, on average, at similar rates. Matching the de-biased local Mbh-Mstar relation requires a mean radiative efficiency ~0.15, in line with theoretical expectations for accretion onto spinning black holes. However, matching the "raw" observed relation for inactive black holes requires a mean radiative efficiency around 0.02, far below theoretical expectations. This result provides independent evidence for selection bias in dynamically estimated black hole masses, a conclusion that is robust to uncertainties in bolometric corrections, obscured active black hole fractions, and kinetic accretion efficiency. For our fiducial assumptions, they favour moderate-to-rapid spins of typical supermassive black holes, to achieve a mean radiative efficiency ~0.12-0.20. Our approach has similarities to the classic Soltan analysis, but by using galaxy-based data instead of integrated quantities we are able to focus on regimes where observational uncertainties are minimized.

preprint2019arXiv

Resolving the Interstellar Medium in Ultraluminous Infrared QSO Hosts with ALMA

We present ALMA observations of the CO(1-0) line and 3-mm continuum emission in eight ultraluminous infrared (IR) quasi-stellar objects (QSOs) at z = 0.06-0.19. All eight IR QSO hosts are clearly resolved in their CO molecular gas emission with a median source size of 3.2 kpc, and seven out of eight sources are detected in 3-mm continuum, which is found to be more centrally concentrated with respect to molecular gas with sizes of 0.4-1.0 kpc. Our observations reveal a diversity of CO morphology and kinematics for the IR QSO systems which can be roughly classified into three categories, rotating gas disk with ordered velocity gradient, compact CO peak with disturbed velocity, and multiple CO distinct sources undergoing a merger between luminous QSO and a companion galaxy separated by a few kpc. The molecular gas in three of IR QSO hosts are found to be rotation-dominated with the ratio of the maximum rotation velocity to the local velocity dispersion of $V_{\rm rot}/σ=4-6$. Basic estimates of the dynamical masses within the CO-emitting regions give masses between $7.4\times10^9$ and $6.9\times10^{10}$ $M_\odot$. We find an increasing trend between BH mass accretion rate and star formation rate (SFR) over three orders of magnitude in far-IR luminosity/SFR, in line with the correlation between QSO bolometric luminosity and SF activity, indicative of a likely direct connection between AGN and SF activity over galaxy evolution timescales.

preprint2019arXiv

The ALMA Spectroscopic Survey in the HUDF: CO emission lines and 3 mm continuum sources

The ALMA SPECtroscopic Survey in the {\it Hubble} Ultra Deep Field is an ALMA large program that obtained a frequency scan in the 3\,mm band to detect emission lines from the molecular gas in distant galaxies. We here present our search strategy for emission lines and continuum sources in the HUDF. We compare several line search algorithms used in the literature, and critically account for the line-widths of the emission line candidates when assessing significance. We identify sixteen emission lines at high fidelity in our search. Comparing these sources to multi-wavelength data we find that all sources have optical/infrared counterparts. Our search also recovers candidates that have lower significance that can be used statistically to derive, e.g. the CO luminosity function. We apply the same detection algorithm to obtain a sample of six 3 mm continuum sources. All of these are also detected in the 1.2 mm continuum with optical/near-infrared counterparts. We use the continuum sources to compute 3 mm number counts in the sub-mJy regime, and find them to be higher by an order of magnitude than expected for synchrotron-dominated sources. However, the number counts are consistent with those derived at shorter wavelengths (0.85--1.3\,mm) once extrapolating to 3\,mm with a dust emissivity index of $β=1.5$, dust temperature of 35\,K and an average redshift of $z=2.5$. These results represent the best constraints to date on the faint end of the 3 mm number counts.

preprint2016arXiv

ALMA spectroscopic survey in the Hubble Ultra Deep Field: CO luminosity functions and the evolution of the cosmic density of molecular gas

In this paper we use ASPECS, the ALMA Spectroscopic Survey in the {\em Hubble} Ultra Deep Field (UDF) in band 3 and band 6, to place blind constraints on the CO luminosity function and the evolution of the cosmic molecular gas density as a function of redshift up to $z\sim 4.5$. This study is based on galaxies that have been solely selected through their CO emission and not through any other property. In all of the redshift bins the ASPECS measurements reach the predicted `knee' of the CO luminosity function (around $5\times10^{9}$ K km/s pc$^2$). We find clear evidence of an evolution in the CO luminosity function with respect to $z\sim 0$, with more CO luminous galaxies present at $z\sim 2$. The observed galaxies at $z\sim 2$ also appear more gas-rich than predicted by recent semi-analytical models. The comoving cosmic molecular gas density within galaxies as a function of redshift shows a factor 3-10 drop from $z \sim 2$ to $z \sim 0$ (with significant error bars), and possibly a decline at $z>3$. This trend is similar to the observed evolution of the cosmic star formation rate density. The latter therefore appears to be at least partly driven by the increased availability of molecular gas reservoirs at the peak of cosmic star formation ($z\sim2$).

preprint2016arXiv

ALMA Spectroscopic Survey in the Hubble Ultra Deep Field: Molecular gas reservoirs in high-redshift galaxies

We study the molecular gas properties of high-$z$ galaxies observed in the ALMA Spectroscopic Survey (ASPECS) that targets a $\sim1$ arcmin$^2$ region in the Hubble Ultra Deep Field (UDF), a blind survey of CO emission (tracing molecular gas) in the 3mm and 1mm bands. Of a total of 1302 galaxies in the field, 56 have spectroscopic redshifts and correspondingly well-defined physical properties. Among these, 11 have infrared luminosities $L_{\rm{}IR}>10^{11}$ L$_\odot$, i.e. a detection in CO emission was expected. Out these, 7 are detected at various significance in CO, and 4 are undetected in CO emission. In the CO-detected sources, we find CO excitation conditions that are lower than typically found in starburst/SMG/QSO environments. We use the CO luminosities (including limits for non-detections) to derive molecular gas masses. We discuss our findings in context of previous molecular gas observations at high redshift (star-formation law, gas depletion times, gas fractions): The CO-detected galaxies in the UDF tend to reside on the low-$L_{\rm{}IR}$ envelope of the scatter in the $L_{\rm{}IR}-L'_{\rm{}CO}$ relation, but exceptions exist. For the CO-detected sources, we find an average depletion time of $\sim$ 1 Gyr, with significant scatter. The average molecular-to-stellar mass ratio ($M_{\rm{}H2}$/$M_*$) is consistent with earlier measurements of main sequence galaxies at these redshifts, and again shows large variations among sources. In some cases, we also measure dust continuum emission. On average, the dust-based estimates of the molecular gas are a factor $\sim$2-5$\times$ smaller than those based on CO. Accounting for detections as well as non-detections, we find large diversity in the molecular gas properties of the high-redshift galaxies covered by ASPECS.

preprint2016arXiv

ALMA Spectroscopic Survey in the Hubble Ultra Deep Field: Survey Description

We present the rationale for and the observational description of ASPECS: The ALMA SPECtroscopic Survey in the Hubble Ultra-Deep Field (UDF), the cosmological deep field that has the deepest multi-wavelength data available. Our overarching goal is to obtain an unbiased census of molecular gas and dust continuum emission in high-redshift (z$>$0.5) galaxies. The $\sim$1$'$ region covered within the UDF was chosen to overlap with the deepest available imaging from HST. Our ALMA observations consist of full frequency scans in band 3 (84-115 GHz) and band 6 (212-272 GHz) at approximately uniform line sensitivity ($L'_{\rm CO}\sim$2$\times$10$^{9}$ K km/s pc$^2$), and continuum noise levels of 3.8 $μ$Jy beam$^{-1}$ and 12.7 $μ$Jy beam$^{-1}$, respectively. The molecular surveys cover the different rotational transitions of the CO molecule, leading to essentially full redshift coverage. The [CII] emission line is also covered at redshifts $6.0<z<8.0$. We present a customized algorithm to identify line candidates in the molecular line scans, and quantify our ability to recover artificial sources from our data. Based on whether multiple CO lines are detected, and whether optical spectroscopic redshifts as well as optical counterparts exist, we constrain the most likely line identification. We report 10 (11) CO line candidates in the 3mm (1mm) band, and our statistical analysis shows that $<$4 of these (in each band) are likely spurious. Less than 1/3 of the total CO flux in the low-J CO line candidates are from sources that are not associated with an optical/NIR counterpart. We also present continuum maps of both the band 3 and band 6 observations. The data presented here form the basis of a number of dedicated studies that are presented in subsequent papers.

preprint2016arXiv

ALMA Spectroscopic Survey in the Hubble Ultra Deep Field: The Infrared Excess of UV-selected z=2-10 galaxies as a function of UV-continuum Slope and Stellar Mass

We make use of deep 1.2mm-continuum observations (12.7microJy/beam RMS) of a 1 arcmin^2 region in the Hubble Ultra Deep Field to probe dust-enshrouded star formation from 330 Lyman-break galaxies spanning the redshift range z=2-10 (to ~2-3 Msol/yr at 1sigma over the entire range). Given the depth and area of ASPECS, we would expect to tentatively detect 35 galaxies extrapolating the Meurer z~0 IRX-beta relation to z>~2 (assuming T_d~35 K). However, only 6 tentative detections are found at z>~2 in ASPECS, with just three at >3sigma. Subdividing z=2-10 galaxies according to stellar mass, UV luminosity, and UV-continuum slope and stacking the results, we only find a significant detection in the most massive (>10^9.75 Msol) subsample, with an infrared excess (IRX=L_{IR}/L_{UV}) consistent with previous z~2 results. However, the infrared excess we measure from our large selection of sub-L* (<10^9.75 Msol) galaxies is 0.11(-0.42)(+0.32) and 0.14(-0.14)(+0.15) at z=2-3 and z=4-10, respectively, lying below even an SMC IRX-beta relation (95% confidence). These results demonstrate the relevance of stellar mass for predicting the IR luminosity of z>~2 galaxies. We furthermore find that the evolution of the IRX-stellar mass relationship depends on the evolution of the dust temperature. If the dust temperature increases monotonically with redshift (as (1+z)^0.32) such that T_d~44-50 K at z>=4, current results are suggestive of little evolution in this relationship to z~6. We use these results to revisit recent estimates of the z>~3 SFR density. One less obvious implication is in interpreting the high Halpha EWs seen in z~5 galaxies: our results imply that star-forming galaxies produce Lyman-continuum photons at twice the efficiency (per unit UV luminosity) as implied in conventional models. Star-forming galaxies can then reionize the Universe, even if the escape fraction is <10%.

preprint2016arXiv

Discovery of a galaxy cluster with a violently starbursting core at z=2.506

We report the discovery of a remarkable concentration of massive galaxies with extended X-ray emission at $z_{spec} = 2.506$, which contains 11 massive ($M_{*} \gtrsim 10^{11} M_{\odot}$) galaxies in the central 80kpc region (11.6$σ$ overdensity). We have spectroscopically confirmed 17 member galaxies with 11 from CO and the remaining ones from $Hα$. The X-ray luminosity, stellar mass content and velocity dispersion all point to a collapsed, cluster-sized dark matter halo with mass $M_{200c} = 10^{13.9\pm0.2} M_{\odot}$, making it the most distant X-ray-detected cluster known to date. Unlike other clusters discovered so far, this structure is dominated by star-forming galaxies (SFGs) in the core with only 2 out of the 11 massive galaxies classified as quiescent. The star formation rate (SFR) in the 80kpc core reaches $\sim$3400 $M_{\odot}$ yr$^{-1}$ with a gas depletion time of $\sim 200$ Myr, suggesting that we caught this cluster in rapid build-up of a dense core. The high SFR is driven by both a high abundance of SFGs and a higher starburst fraction ($\sim25\%$, compared to 3\%-5\% in the field). The presence of both a collapsed, cluster-sized halo and a predominant population of massive SFGs suggests that this structure could represent an important transition phase between protoclusters and mature clusters. It provides evidence that the main phase of massive galaxy passivization will take place after galaxies accrete onto the cluster, providing new insights into massive cluster formation at early epochs. The large integrated stellar mass at such high redshift challenges our understanding of massive cluster formation.

preprint2016arXiv

Observational evidence of a slow downfall of star formation efficiency in massive galaxies during the last 10 Gyr

In this paper we study the causes of the reported mass-dependence of the slope of SFR-M* relation, the so-called "Main Sequence" of star-forming galaxies, and discuss its implication on the physical processes that shaped the star formation history of massive galaxies over cosmic time. We use the CANDELS near-IR imaging from the Hubble Space Telescope to perform the bulge-to-disk decomposition of distant galaxies and measure for the first time the slope of the SFR-Mdisk relation at z=1. We find that this relation follows very closely the shape of the SFR-M* correlation, still with a pronounced flattening at the high-mass end. This is clearly excluding, at least at z=1, the secular growth of quiescent bulges in star-forming galaxies as the main driver for the change of slope of the Main Sequence. Then, by stacking the Herschel data available in the CANDELS field, we estimate the total gas mass and the star formation efficiency at different positions on the SFR-M* relation. We find that the relatively low SFRs observed in massive galaxies (M* > 5e10 Msun) are caused by a decreased star formation efficiency, by up to a factor of 3 as compared to lower stellar mass galaxies, and not by a reduced gas content. The trend at the lowest masses is likely linked to the dominance of atomic over molecular gas. We argue that this stellar-mass-dependent SFE can explain the varying slope of the Main Sequence since z=1.5, hence over 70% of the Hubble time. The drop of SFE occurs at lower masses in the local Universe (M* > 2e10 Msun) and is not present at z=2. Altogether this provides evidence for a slow downfall of the star formation efficiency in massive Main Sequence galaxies. The resulting loss of star formation is found to be rising starting from z=2 to reach a level comparable to the mass growth of the quiescent population by z=1. We finally discuss the possible physical origin of this phenomenon.

preprint2015arXiv

A Physical Approach to the Identification of High-z Mergers: Morphological Classification in the Stellar Mass Domain

At z>1, the distinction between merging and 'normal' star-forming galaxies based on single band morphology is often hampered by the presence of large clumps which result in a disturbed, merger-like appearance even in rotationally supported disks. In this paper we discuss how a classification based on canonical, non-parametric structural indices measured on resolved stellar mass maps, rather than on single-band images, reduces the misclassification of clumpy but not merging galaxies. We calibrate the mass-based selection of mergers using the MIRAGE hydrodynamical numerical simulations of isolated and merging galaxies which span a stellar mass range of $10^{9.8}$-$10^{10.6}M_{sun}$ and merger ratios between 1:1-1:6.3. These simulations are processed to reproduce the typical depth and spatial resolution of observed HUDF data. We test our approach on a sample of real z~2 galaxies with kinematic classification into disks or mergers and on ~100 galaxies in the HUDF field with photometric/spectroscopic redshift between 1.5$\leqslant z \leqslant$3 and $M>10^{9.4}M_{sun}$. We find that a combination of the asymmetry $A_{MASS}$ and $M_{20, MASS}$ indices measured on the stellar mass maps can efficiently identify real (major) mergers with $\lesssim$20% contamination from clumpy disks in the merger sample. This mass-based classification cannot be reproduced in star-forming galaxies by $H-$band measurements alone, which instead result in a contamination from clumpy galaxies that can be as high as 50%. Moreover, we find that the mass-based classification always results in a lower contamination from clumpy galaxies than an $H-$band classification, regardless of the depth of the imaging used (e.g., CANDELS vs. HUDF).

preprint2015arXiv

GOODS-HERSCHEL: star formation, dust attenuation and the FIR-radio correlation on the Main Sequence of star-forming galaxies up to z~4

We use deep panchromatic datasets in the GOODS-N field, from GALEX to the deepest Herschel far-infrared and VLA radio continuum imaging, to explore, using mass-complete samples, the evolution of the star formation activity and dust attenuation of star-forming galaxies to z~4. Our main results can be summarized as follows: i) the slope of the SFR-M correlation is consistent with being constant, and equal to ~0.8 at least up to z~1.5, while its normalization keeps increasing with redshift; ii) for the first time here we are able to explore the FIR-radio correlation for a mass-selected sample of star-forming galaxies: the correlation does not evolve up to z~4; iii) we confirm that galaxy stellar mass is a robust proxy for UV dust attenuation in star-forming galaxies, with more massive galaxies being more dust attenuated, strikingly we find that this attenuation relation evolves very weakly with redshift, the amount of dust attenuation increasing by less than 0.3 magnitudes over the redshift range [0.5-4] for a fixed stellar mass, as opposed to a tenfold increase of star formation rate; iv) the correlation between dust attenuation and the UV spectral slope evolves in redshift, with the median UV spectral slope of star-forming galaxies becoming bluer with redshift. By z~3, typical UV slopes are inconsistent, given the measured dust attenuation, with the predictions of commonly used empirical laws. Finally, building on existing results, we show that gas reddening is marginally larger (by a factor of around 1.3) than stellar reddening at all redshifts probed, and also that the amount of dust attenuation at a fixed ISM metallicity increases with redshift. We speculate that our results support evolving ISM conditions of typical star-forming galaxies such that at z~1.5 Main Sequence galaxies have ISM conditions getting closer to those of local starbursts.

preprint2015arXiv

High-J CO Versus Far-Infrared Relations in Normal and Starburst Galaxies

We present correlations between 9 CO transition ($J=4-3$ to $12-11$) and beam-matched far-infrared (Far-IR) luminosities ($L_{\mathrm{FIR},\,b}$) among 167 local galaxies, using {\it{Herschel}} Spectral and Photometric Imaging Receiver Fourier Transform Spectrometer (SPIRE; FTS) spectroscopic data and Photoconductor Array Camera and Spectrometer (PACS) photometry data. We adopt entire-galaxy FIR luminosities ($L_{\mathrm{FIR},\,e}$) from the {\it{IRAS}} Revised Bright Galaxy Sample and correct to $L_{\mathrm{FIR},\,b}$ using PACS images to match the varying FTS beam sizes. All 9 correlations between $L'_{\mathrm{CO}}$ and $L_{\mathrm{FIR},\,b}$ are essentially linear and tight ($σ$=0.2-0.3 dex dispersion), even for the highest transition, $J=12-11$. This supports the notion that the star formation rate (SFR) is linearly correlated with the dense molecular gas ($n_{\mathrm{H}_2}\gtrsim10^{4-6}\,cm^{-3}$). We divide the entire sample into three subsamples and find that smaller sample sizes can induce large differences in the correlation slopes. We also derive an average CO spectral line energy distribution (SLED) for the entire sample and discuss the implied average molecular gas properties for these local galaxies. We further extend our sample to high-{\it{z}} galaxies with CO $J=5-4$ data from the literature as an example, including submillimeter galaxies (SMGs) and "normal" star-forming BzKs. BzKs have similar FIR/CO(5-4) ratios as that of local galaxies, an follow well the locally-determined correlation, whereas SMG ratios fall around or slightly above the local correlation with large uncertainties. Finally, by including Galactic CO($J=10-9$) data as well as very limited high-{\it{z}} CO $J=10-9$ data, we verify that the CO(10-9) -- FIR correlation successfully extends to Galactic young stellar objects, suggesting that linear correlations are valid over 15 orders of magnitude.

preprint2015arXiv

Old age and super-solar metallicity in a massive z~1.4 early-type galaxy from VLT/X-Shooter spectroscopy

We present the first estimate of age, stellar metallicity and chemical abundance ratios, for an individual early-type galaxy at high-redshift (z = 1.426) in the COSMOS field. Our analysis is based on observations obtained with the X-Shooter instrument at the VLT, which cover the visual and near infrared spectrum at high (R >5000) spectral resolution. We measure the values of several spectral absorptions tracing chemical species, in particular Magnesium and Iron, besides determining the age-sensitive D4000 break. We compare the measured indices to stellar population models, finding good agreement. We find that our target is an old (t > 3 Gyr), high-metallicity ([Z/H] > 0.5) galaxy which formed its stars at z_{form} > 5 within a short time scale ~0.1 Gyr, as testified by the strong [α/Fe] ratio ( > 0.4), and has passively evolved in the first > 3-4 Gyr of its life. We have verified that this result is robust against the choice and number of fitted spectral features, and stellar population model. The result of an old age and high-metallicity has important implications for galaxy formation and evolution confirming an early and rapid formation of the most massive galaxies in the Universe.

preprint2015arXiv

SHARDS: A global view of the star formation activity at z~0.84 and z~1.23

In this paper, we present a comprehensive analysis of star-forming galaxies (SFGs) at intermediate redshifts (z~1). We combine the ultra-deep optical spectro-photometric data from the Survey for High-z Absorption Red and Dead Sources (SHARDS) with deep UV-to-FIR observations in the GOODS-N field. Exploiting two of the 25 SHARDS medium-band filters, F687W17 and F823W17, we select [OII] emission line galaxies at z~0.84 and z~1.23 and characterize their physical properties. Their rest-frame equivalent widths (EW$_{\mathrm{rf}}$([OII])), line fluxes, luminosities, star formation rates (SFRs) and dust attenuation properties are investigated. The evolution of the EW$_{\mathrm{rf}}$([OII]) closely follows the SFR density evolution of the universe, with a trend of EW$_{\mathrm{rf}}$([OII])$\propto$(1+z)$^3$ up to redshift z~1, followed by a possible flattening. The SF properties of the galaxies selected on the basis of their [OII] emission are compared with complementary samples of SFGs selected by their MIR and FIR emission, and also with a general mass-selected sample of galaxies at the same redshifts. We demonstrate observationally that the UVJ diagram (or, similarly, a cut in the specific SFR) is only partially able to distinguish the quiescent galaxies from the SFGs. The SFR-M$_*$ relation is investigated for the different samples, yelding a logarithmic slope ~1, in good agreement with previous results. The dust attenuations derived from different SFR indicators (UV(1600), UV(2800), [OII], IR) are compared and show clear trends with respect to both the stellar mass and total SFR, with more massive and highly star-forming galaxies being affected by stronger dust attenuation.

preprint2015arXiv

Star formation and quenching among the most massive galaxies at z~1.7

We have conducted a detailed object-by-object study of a mass-complete (M*>10^11 M_sun) sample of 56 galaxies at 1.4 < z < 2 in the GOODS-South field, showing that an accurate de-blending in MIPS/24um images is essential to properly assign to each galaxy its own star formation rate (SFR), whereas an automatic procedure often fails. This applies especially to galaxies with SFRs below the Main Sequence (MS) value, which may be in their quenching phase. After that, the sample splits evenly between galaxies forming stars within a factor of 4 of the MS rate (~45%), and sub-MS galaxies with SFRs ~10-1000 times smaller (~55%). We did not find a well defined class of intermediate, transient objects below the MS, suggesting that the conversion of a massive MS galaxy into a quenched remnant may take a relatively short time (<1 Gyr), though a larger sample should be analyzed in the same way to set precise limits on the quenching timescale. X-ray detected AGNs represent a ~30% fraction of the sample, and are found among both star-forming and quenched galaxies. The morphological analysis revealed that ~50% of our massive objects are bulge-dominated, and almost all MS galaxies with a relevant bulge component host an AGN. We also found sub-MS SFRs in many bulge-dominated systems, providing support to the notion that bulge growth, AGN activity and quenching of star formation are closely related to each other.

preprint2015arXiv

The Herschel view of the dominant mode of galaxy growth from z=4 to the present day

We present an analysis of the deepest Herschel images in four major extragalactic fields GOODS-North, GOODS-South, UDS and COSMOS obtained within the GOODS-Herschel and CANDELS-Herschel key programs. The picture provided by 10497 individual far-infrared detections is supplemented by the stacking analysis of a mass-complete sample of 62361 star-forming galaxies from the CANDELS-HST H band-selected catalogs and from two deep ground-based Ks band-selected catalogs in the GOODS-North and the COSMOS-wide fields, in order to obtain one of the most accurate and unbiased understanding to date of the stellar mass growth over the cosmic history. We show, for the first time, that stacking also provides a powerful tool to determine the dispersion of a physical correlation and describe our method called "scatter stacking" that may be easily generalized to other experiments. We demonstrate that galaxies of all masses from z=4 to 0 follow a universal scaling law, the so-called main sequence of star-forming galaxies. We find a universal close-to-linear slope of the logSFR-logM* relation with evidence for a flattening of the main sequence at high masses (log(M*/Msun) > 10.5) that becomes less prominent with increasing redshift and almost vanishes by z~2. This flattening may be due to the parallel stellar growth of quiescent bulges in star-forming galaxies. Within the main sequence, we measure a non varying SFR dispersion of 0.3 dex. The specific SFR (sSFR=SFR/M*) of star-forming galaxies is found to continuously increase from z=0 to 4. Finally we discuss the implications of our findings on the cosmic SFR history and show that more than 2/3 of present-day stars must have formed in a regime dominated by the main sequence mode. As a consequence we conclude that, although omnipresent in the distant Universe, galaxy mergers had little impact in shaping the global star formation history over the last 12.5 Gyr.

preprint2015arXiv

The influence of wavelength, flux, and lensing selection effects on the redshift distribution of dusty, star-forming galaxies

We interpret the large variety of redshift distributions of galaxies found by far-infrared and (sub-)millimeter deep surveys depending on their depth and wavelength using the Béthermin et al. (2012) phenomenological model of galaxy evolution. This model reproduces without any new parameter tuning the observed redshift distributions from 100 $μ$m to 1.4 mm, and especially the increase of the median redshift with survey wavelength. This median redshift varies also significantly with the depth of the surveys, and deeper surveys do necessarily not probe higher redshifts. Paradoxically, at fixed wavelength and flux limit, the lensed sources are not always at higher redshift. We found that the higher redshift of 1.4 mm-selected south pole telescope (SPT) sources compared to other SMG surveys is not only caused by the lensing selection, but also by the longer wavelength. This SPT sample is expected to be dominated by a population of lensed main-sequence galaxies and a minor contribution ($\sim$10\%) of unlensed extreme starbursts.

preprint2015arXiv

The morphologies of massive galaxies from z~3 - Witnessing the 2 channels of bulge growth

[abridged] We quantify the morphological evolution of z~0 massive galaxies ($M*/M_\odot\sim10^{11}$) from z~3 in the 5 CANDELS fields. The progenitors are selected using abundance matching techniques to account for the mass growth. The morphologies strongly evolve from z~3. At z<1, the population matches the massive end of the Hubble sequence, with 30% of spheroids, 50% of galaxies with equally dominant disk and bulge components and 20% of disks. At z~2-3 there is a majority of irregular systems (~60-70%) with still 30% of spheroids. We then analyze the SFRs, gas fractions and structural properties for the different morphologies independently. Our results suggest two distinct channels for the growth of bulges in massive galaxies. Around 30-40% were already bulges at z~2.5, with low average SFRs and gas-fractions (10-15%), high Sersic indices (n>3-4) and small effective radii ($R_e$~1 kpc) pointing towards an early formation through gas-rich mergers or VDI. Between z~ 2.5 and z~0, they rapidly increase their size by a factor of ~4-5, become all passive but their global morphology remains unaltered. The structural evolution is independent of the gas fractions, suggesting that it is driven by ex-situ events. The remaining 60% experience a gradual morphological transformation, from clumpy disks to more regular bulge+disks systems, essentially happening at z>1. It results in the growth of a significant bulge component (n~3) for 2/3 of the systems possibly through the migration of clumps while the remaining 1/3 keeps a rather small bulge (n~1.5-2). The transition phase between disturbed and relaxed systems and the emergence of the bulge is correlated with a decrease of the star formation activity and the gas fractions. The growth of the effective radii scales roughly with $H(z)^{-1}$ and it is therefore consistent with the expected growth of disks in galaxy haloes.

preprint2014arXiv

AGN Emission Line Diagnostics and the Mass-Metallicity Relation up to Redshift z~2: the Impact of Selection Effects and Evolution

Emission line diagnostic diagrams probing the ionization sources in galaxies, such as the Baldwin-Phillips-Terlevich (BPT) diagram, have been used extensively to distinguish AGN from purely star-forming galaxies. Yet, they remain poorly understood at higher redshifts. We shed light on this issue with an empirical approach based on a z~0 reference sample built from ~300,000 SDSS galaxies, from which we mimic selection effects due to typical emission line detection limits at higher redshift. We combine this low-redshift reference sample with a simple prescription for luminosity evolution of the global galaxy population to predict the loci of high-redshift galaxies on the BPT and Mass-Excitation (MEx) diagnostic diagrams. The predicted bivariate distributions agree remarkably well with direct observations of galaxies out to z~1.5, including the observed stellar mass-metallicity (MZ) relation evolution. As a result, we infer that high-redshift star-forming galaxies are consistent with having "normal" ISM properties out to z~1.5, after accounting for selection effects and line luminosity evolution. Namely, their optical line ratios and gas-phase metallicities are comparable to that of low-redshift galaxies with equivalent emission-line luminosities. In contrast, AGN narrow-line regions may show a shift toward lower metallicities at higher redshift. While a physical evolution of the ISM conditions is not ruled out for purely star-forming galaxies, and may be more important starting at z>2, we find that reliably quantifying this evolution is hindered by selections effects. The recipes provided here may serve as a basis for future studies toward this goal. Code to predict the loci of galaxies on the BPT and MEx diagnostic diagrams, and the MZ relation as a function of emission line luminosity limits, is made publicly available.

preprint2014arXiv

Clustering, host halos and environment of z$\sim$2 galaxies as a function of their physical properties

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

preprint2014arXiv

Dust and gas in luminous proto-cluster galaxies at z=4.05: the case for different cosmic dust evolution in normal and starburst galaxies

We measure the dust and gas content of the three sub-millimeter galaxies (SMGs) in the GN20 proto-cluster at $z=4.05$ using new IRAM Plateau de Bure interferometer (PdBI) CO(4-3) and 1.2-3.3 mm continuum observations. All these three SMGs are heavily dust obscured, with UV-based star formation rate (SFR) estimates significantly smaller than the ones derived from the bolometric infrared (IR), consistent with the spatial offsets revealed by HST and CO imaging. Based also on evaluations of the specific SFR, CO-to-H$_2$ conversion factor and gas depletion timescale, we classify all the three galaxies as starbursts (SBs), although with a lower confidence for GN20.2b that might be a later stage merging event. We place our measurements in the context of the evolutionary properties of main sequence (MS) and SB galaxies. ULIRGs have 3-5 times larger $L'_{\rm CO}/M_{\rm dust}$ and $M_{\rm dust}/M_\star$ ratios than $z=0$ MS galaxies, but by $z\sim2$ the difference appears to be blurred, probably due to differential metallicity evolution. SB galaxies appear to slowly evolve in their $L'_{\rm CO}/M_{\rm dust}$ and $M_{\rm dust}/M_\star$ ratios all the way to $z>6$ (consistent with rapid enrichment of SB events), while MS galaxies rapidly increase in $M_{\rm dust}/M_\star$ from $z=0$ to 2 (due to gas fraction increase, compensated by a decrease of metallicities). While no IR/submm continuum detection is available for indisputably normal massive galaxies at $z>2.5$, we show that if metallicity indeed decrease rapidly for these systems at $z>3$ as claimed in the literature, we should expect a strong decrease of their $M_{\rm dust}/M_\star$, consistent with recent PdBI and ALMA upper limits. We conclude that the $M_{\rm dust}/M_\star$ ratio could be a powerful tool for distinguishing starbursts from normal galaxies at $z>4$.

preprint2014arXiv

Evolution of the dust emission of massive galaxies up to z=4 and constraints on their dominant mode of star formation

We aim to measure the average dust and molecular gas content of massive star-forming galaxies ($\rm > 3 \times 10^{10}\,M_\odot$) up to z=4 in the COSMOS field to determine if the intense star formation observed at high redshift is induced by major mergers or caused by large gas reservoirs. Firstly, we measured the evolution of the average spectral energy distributions as a function of redshift using a stacking analysis of Spitzer, Herschel, LABOCA, and AzTEC data for two samples of galaxies: normal star-forming objects and strong starbursts, as defined by their distance to the main sequence. We found that the mean intensity of the radiation field $< U >$ heating the dust (strongly correlated with dust temperature) increases with increasing redshift up to z$\sim$4 in main-sequence galaxies. We can reproduce this evolution with simple models that account for the decrease of the gas metallicity with redshift. No evolution of $< U >$ with redshift is found in strong starbursts. We then deduced the evolution of the molecular gas fraction (defined here as $\rm M_{\rm mol}/(M_{\rm mol}+M_\star)$) with redshift and found a similar, steeply increasing trend for both samples. At z$\sim$4, this fraction reaches $\sim$60%. The average position of the main-sequence galaxies is on the locus of the local, normal star-forming disks in the integrated Schmidt-Kennicutt diagram (star formation rate versus mass of molecular gas), suggesting that the bulk of the star formation up to z=4 is dominated by secular processes.

preprint2014arXiv

Polycyclic Aromatic Hydrocarbon and Mid-Infrared Continuum Emission in a z>4 Submillimeter Galaxy

We report the detection of 6.2um polycyclic aromatic hydrocarbon (PAH) and rest-frame 4-7um continuum emission in the z=4.055 submillimeter galaxy GN20, using the Infrared Spectrograph (IRS) on board the Spitzer Space Telescope. This represents the first detection of PAH emission at z>4. The strength of the PAH emission feature is consistent with a very high star formation rate of ~1600Msun/yr. We find that this intense starburst powers at least ~1/3 of the faint underlying 6um continuum emission, with an additional, significant (and perhaps dominant) contribution due to a power-law-like hot dust source, which we interpret to likely be a faint, dust-obscured active galactic nucleus (AGN). The inferred 6um AGN continuum luminosity is consistent with a sensitive upper limit on the hard X-ray emission as measured by the Chandra X-Ray Observatory if the previously undetected AGN is Compton-thick. This is in agreement with the finding at optical/infrared wavelengths that the galaxy and its nucleus are heavily dust-obscured. Despite the strong power-law component enhancing the mid-infrared continuum emission, the intense starburst associated with the photon-dominated regions that give rise to the PAH emission appears to dominate the total energy output in the infrared. GN20 is one of the most luminous starburst galaxies known at any redshift, embedded in a rich protocluster of star-forming galaxies. This investigation provides an improved understanding of the energy sources that power such exceptional systems, which represent the extreme end of massive galaxy formation at early cosmic times.

preprint2013arXiv

A deep search for molecular gas in two massive Lyman break galaxies at z=3 and 4: vanishing CO-emission due to low metallicity?

We present deep IRAM Plateau de Bure Interferometer (PdBI) observations, searching for CO-emission toward two massive, non-lensed Lyman break galaxies (LBGs) at z=3.216 and 4.058. With one low significance CO detection (3.5 sigma) and one sensitive upper limit, we find that the CO lines are >~ 3-4 times weaker than expected based on the relation between IR and CO luminosities followed by similarly, massive galaxies at z=0-2.5. This is consistent with a scenario in which these galaxies have low metallicity, causing an increased CO-to-H_2 conversion factor, i.e., weaker CO-emission for a given molecular (H_2) mass. The required metallicities at z>3 are lower than predicted by the fundamental metallicity relation (FMR) at these redshifts, consistent with independent evidence. Unless our galaxies are atypical in this respect, detecting molecular gas in normal galaxies at z>3 may thus remain challenging even with ALMA.

preprint2013arXiv

Empirical predictions for (sub-)millimeter line and continuum deep fields

[abridged] Modern (sub-)millimeter/radio interferometers will enable us to measure the dust and molecular gas emission from galaxies that have luminosities lower than the Milky Way, out to high redshifts and with unprecedented spatial resolution and sensitivity. This will provide new constraints on the star formation properties and gas reservoir in galaxies throughout cosmic times through dedicated deep field campaigns targeting the CO/[CII] lines and dust continuum emission. In this paper, we present empirical predictions for such (sub-)millimeter line and continuum deep fields. We base these predictions on the deepest available optical/near-infrared ACS and NICMOS data on the Hubble Ultra Deep Field. Using a physically-motivated spectral energy distribution model, we fit the observed optical/near-infrared emission of 13,099 galaxies with redshifts up to z=5, and obtain median likelihood estimates of their stellar mass, star formation rate, dust attenuation and dust luminosity. We derive statistical constraints on the dust emission in the infrared and (sub-)millimeter which are consistent with the observed optical/near-infrared emission in terms of energy balance. This allows us to estimate, for each galaxy, the (sub-)millimeter continuum flux densities in several ALMA, PdBI/NOEMA and JVLA bands. Using empirical relations between the observed CO/[CII] line luminosities and the infrared luminosity, we infer the flux of the CO(1-0) and [CII] lines from the estimated infrared luminosity of each galaxy in our sample. We then predict the fluxes of higher CO transition lines CO(2-1) to CO(7-6) bracketing two extreme gas excitation scenarios. We use our predictions to discuss possible deep field strategies with ALMA. The predictions presented in this study will serve as a direct benchmark for future deep field campaigns in the (sub-)millimeter regime.

preprint2013arXiv

On the effect of the Cosmic Microwave Background in high-redshift (sub-)millimeter observations

Modern (sub-)millimeter interferometers enable the measurement of the cool gas and dust emission of high-redshift galaxies (z>5). However, at these redshifts the cosmic microwave background (CMB) temperature is higher, approaching, and even exceeding, the temperature of cold dust and molecular gas observed in the local Universe. In this paper, we discuss the impact of the warmer CMB on (sub-)millimeter observations of high-redshift galaxies. The CMB affects the observed (sub-)millimeter dust continuum and the line emission (e.g. carbon monoxide, CO) in two ways: (i) it provides an additional source of (both dust and gas) heating; and (ii) it is a non-negligible background against which the line and continuum emission are measured. We show that these two competing processes affect the way we interpret the dust and gas properties of high-redshift galaxies using spectral energy distribution models. We quantify these effects and provide correction factors to compute what fraction of the intrinsic dust (and line) emission can be detected against the CMB as a function of frequency, redshift and temperature. We discuss implications on the derived properties of high-redshift galaxies from (sub-)millimeter data. Specifically, the inferred dust and molecular gas masses can be severely underestimated for cold systems if the impact of the CMB is not properly taken into account.

preprint2013arXiv

Probing the Interstellar Medium of z~1 Ultra-luminous Infrared Galaxies through Interferometric Observations of CO and Spitzer Mid-infrared Spectroscopy

We explore the relationship between gas, dust and star formation in a sample of 12 ultra-luminous infrared galaxies (ULIRGs) at high redshift compared to a similar sample of local galaxies. We present new CO observations and/or Spitzer mid-IR spectroscopy for 6 70 micron selected galaxies at z~1 in order to quantify the properties of the molecular gas reservoir, the contribution of an active galactic nuclei (AGN) to the mid-IR luminosity and the star formation efficiency (SFE=LIR/L'CO). The mid-IR spectra show strong polycyclic aromatic hydrocarbon (PAH) emission and our spectral decomposition suggests that the AGN makes a minimal contribution (<25%) to the mid-IR luminosity. The 70 micron selected ULIRGs which we find to be spectroscopic close pairs, are observed to have high SFE, similar to local ULIRGs and high redshift submillimeter galaxies, consistent with enhanced IR luminosity due to an ongoing major merger. Combined with existing observations of local and high redshift ULIRGs, we further compare the PAH, IR and CO luminosities. We show that the ratio LPAH6.2/LIR decreases with increasing IR luminosity for both local and high redshift galaxies but the trend for high redshift galaxies is shifted to higher IR luminosities; the average LPAH6.2/LIR ratio at a given LIR is ~3 times higher at high redshift. When we normalize by the molecular gas, we find this trend to be uniform for galaxies at all redshifts and that the molecular gas is correlated with the PAH dust emission.The similar trends seen in the [CII] to molecular gas ratios in other studies suggests that PAH emission, like [CII], continues to be a good tracer of photodissociation regions even at high redshift. Together the CO, PAH and far-IR fine structure lines should be useful for constraining the interstellar medium conditions in high redshift galaxies.

preprint2013arXiv

SHARDS: stellar populations and star formation histories of a mass-selected sample of 0.65<z<1.1 galaxies

We report on results from the analysis of a stellar mass-selected (log M*>9.0) sample of 1644 galaxies at 0.65<z<1.1 with ultra-deep (m<26.5) optical medium-band (R~50) photometry from the Survey for High-z Absorption Red and Dead Sources (SHARDS). The spectral resolution of SHARDS allows us to consistently measure the strength of the 4000 Angstrom spectral break [Dn(4000), an excellent age indicator for the stellar populations of quiescent galaxies] for all galaxies at z~0.9 down to log M*9. The Dn(4000) index cannot be resolved from broad-band photometry, and measurements from optical spectroscopic surveys are typically limited to galaxies at least x10 more massive. When combined with the rest-frame U-V colour, Dn(4000) provides a powerful diagnostic of the extinction affecting the stellar population that is relatively insensitive to degeneracies with age, metallicity or star formation history. We use this novel approach to estimate de-reddened colours and light-weighted stellar ages for individual sources. We explore the relationships linking stellar mass, (U-V), and Dn(4000) for the sources in the sample, and compare them to those found in local galaxies. The main results are: a) both Dn(4000) and (U-V) correlate with M*. The dispersion in Dn(4000) values at a given M* increases with M*, while the dispersion for (U-V) decreases due to the higher average extinction prevalent in massive star-forming galaxies. b) for massive galaxies, we find a smooth transition between the blue cloud and red sequence in the intrinsic U-V colour, in contrast with other recent results. c) at a fixed stellar age, we find a positive correlation between extinction and stellar mass. d) the fraction of sources with declining or halted star formation increases steeply with the stellar mass, from ~5% at log M*~9.0-9.5 to ~80% at log M*>11, in agreement with downsizing scenarios.

preprint2013arXiv

The long lives of giant clumps and the birth of outflows in gas-rich galaxies at high redshift

Star-forming disk galaxies at high redshift are often subject to violent disk instability, characterized by giant clumps whose fate is yet to be understood. The main question is whether the clumps disrupt within their dynamical timescale (<50Myr), like molecular clouds in today's galaxies, or whether they survive stellar feedback for more than a disk orbital time (~300Myr) in which case they can migrate inward and help building the central bulge. We present 3.5-7pc resolution AMR simulations of high-redshift disks including photo-ionization, radiation pressure, and supernovae feedback (Renaud et al. 2013, and Perret et al., this astro-ph issue). Our modeling of radiation pressure determines the mass loading and initial velocity of winds from basic physical principles. We find that the giant clumps produce steady outflow rates comparable to and sometimes somewhat larger than their star formation rate, with velocities largely sufficient to escape galaxy. The clumps also lose mass, especially old stars, by tidal stripping, and the stellar populations contained in the clumps hence remain relatively young (<=200Myr), as observed. The clumps survive gaseous outflows and stellar loss, because they are wandering in gas-rich turbulent disks from which they can re-accrete gas at high rates compensating for outflows and tidal stripping, overall keeping realistic and self-regulated gaseous and stellar masses. Our simulations produce gaseous outflows with velocities, densities and mass loading consistent with observations, and at the same time suggest that the giant clumps survive for hundreds of Myr and complete their migration to the center of high-redshift galaxies, without rapid dispersion and reformation of clumps.

preprint2013arXiv

The redshift-evolution of the distribution of star formation among dark matter halos as seen in the infrared

[Abridged] Recent studies revealed a strong correlation between the star formation rate (SFR) and stellar mass of star-forming galaxies, the so-called star-forming main sequence. An empirical modeling approach (2-SFM) which distinguishes between the main sequence and rarer starburst galaxies is capable of reproducing most statistical properties of infrared galaxies. In this paper, we extend this approach by establishing a connection between stellar mass and halo mass with the technique of abundance matching. Based on a few, simple assumptions and a physically motivated formalism, our model successfully predicts the (cross-)power spectra of the cosmic infrared background (CIB), the cross-correlation between CIB and cosmic microwave background (CMB) lensing, and the correlation functions of bright, resolved infrared galaxies measured by Herschel, Planck, ACT and SPT. We use this model to infer the redshift distribution these observables, as well as the level of correlation between CIB-anisotropies at different wavelengths. We also predict that more than 90% of cosmic star formation activity occurs in halos with masses between 10^11.5 and 10^13.5 Msun. Taking into account subsequent mass growth of halos, this implies that the majority of stars were initially (at z>3) formed in the progenitors of clusters, then in groups at 0.5<z<3 and finally in Milky-Way-like halos at z<0.5. At all redshifts, the dominant contribution to the star formation rate density stems from halos of mass ~10^12 Msun, in which the instantaneous star formation efficiency is maximal (~70%). The strong redshift-evolution of SFR in the galaxies that dominate the CIB is thus plausibly driven by increased accretion from the cosmic web onto halos of this characteristic mass scale. Material available at http://irfu.cea.fr/Sap/Phocea/Page/index.php?id=537.

preprint2013arXiv

WFC3 grism confirmation of the distant cluster Cl J1449+0856 at z=2.00: Quiescent and star-forming galaxy populations

We present deep Hubble Space Telescope Wide Field Camera 3 slitless spectroscopic observations of the distant cluster Cl J1449+0856. These cover a single pointing with 18 orbits of G141 spectroscopy and F140W imaging, allowing us to derive secure redshifts down to m_140~25.5 AB and 3sigma line fluxes of 5*10^(-18) erg/s/cm^2. In particular, we were able to spectroscopically confirm 12 early-type galaxies in the field up to z~3, 6 of which in the cluster core, which represents the first direct spectroscopic confirmation of passive galaxies in a z=2 cluster environment. With 140 redshifts in a ~6 arcmin^2 field, we can trace the spatial and redshift galaxy distribution in the cluster core and background field. We find two strong peaks at z=2.00 and z=2.07, where only one was seen in our previously published ground-based data. Thanks to the spectroscopic confirmation of the cluster ETGs, we can now re-evaluate the redshift of Cl J1449+0856 at z=2.00, rather than z=2.07, with the background overdensity being revealed to be sparse and "sheet"-like. This presents an interesting case of chance alignment of two close yet unrelated structures, each one preferentially selected by different observing strategies. With 6 quiescent or early-type spectroscopic members and 20 star-forming ones, Cl J1449+0856 is now reliably confirmed to be at z=2.00. The identified members can now allow for a detailed study of galaxy properties in the densest environment at z=2.

preprint2013arXiv

Widespread and Hidden Active Galactic Nuclei in Star-Forming Galaxies at Redshift > 0.3

We characterize the incidence of active galactic nuclei (AGNs) is 0.3 < z < 1 star-forming galaxies by applying multi-wavelength AGN diagnostics (X-ray, optical, mid-infrared, radio) to a sample of galaxies selected at 70-micron from the Far-Infrared Deep Extragalactic Legacy survey (FIDEL). Given the depth of FIDEL, we detect "normal" galaxies on the specific star formation rate (sSFR) sequence as well as starbursting systems with elevated sSFR. We find an overall high occurrence of AGN of 37+/-3%, more than twice as high as in previous studies of galaxies with comparable infrared luminosities and redshifts but in good agreement with the AGN fraction of nearby (0.05 < z < 0.1) galaxies of similar infrared luminosities. The more complete census of AGNs comes from using the recently developed Mass-Excitation (MEx) diagnostic diagram. This optical diagnostic is also sensitive to X-ray weak AGNs and X-ray absorbed AGNs, and reveals that absorbed active nuclei reside almost exclusively in infrared-luminous hosts. The fraction of galaxies hosting an AGN appears to be independent of sSFR and remains elevated both on the sSFR sequence and above. In contrast, the fraction of AGNs that are X-ray absorbed increases substantially with increasing sSFR, possibly due to an increased gas fraction and/or gas density in the host galaxies.

preprint2012arXiv

A unified empirical model for infrared galaxy counts based on observed physical evolution of distant galaxies

We reproduce the mid-infrared to radio galaxy counts with a new empirical model based on our current understanding of the evolution of main-sequence (MS) and starburst (SB) galaxies. We rely on a simple Spectral Energy Distribution (SED) library based on Herschel observations: a single SED for the MS and another one for SB, getting warmer with redshift. Our model is able to reproduce recent measurements of galaxy counts performed with Herschel, including counts per redshift slice. This agreement demonstrates the power of our 2 Star-Formation Modes (2SFM) decomposition for describing the statistical properties of infrared sources and their evolution with cosmic time. We discuss the relative contribution of MS and SB galaxies to the number counts at various wavelengths and flux densities. We also show that MS galaxies are responsible for a bump in the 1.4 GHz radio counts around 50 μJy. Material of the model (predictions, SED library, mock catalogs...) is available online at http://irfu.cea.fr/Sap/Phocea/Page/index.php?id=537.

preprint2012arXiv

An Observed Link between Active Galactic Nuclei and Violent Disk Instabilities in High-Redshift Galaxies

We provide evidence for a correlation between the presence of giant clumps and the occurrence of active galactic nuclei (AGN) in disk galaxies. Giant clumps of 10^8-9 Msun arise from violent gravitational instability in gas-rich galaxies, and it has been proposed that this instability could feed supermassive black holes (BH). We use emission line diagnostics to compare a sample of 14 clumpy (unstable) disks and a sample of 13 smoother (stable) disks at redshift z~0.7. The majority of clumpy disks in our sample have a high probability of containing AGN. Their [OIII] emission line is strongly excited, inconsistent with low-metallicity star formation alone. [NeIII] excitation is also higher. Stable disks rarely have such properties. Stacking ultra sensitive Chandra observations (4 Ms) reveals an X-ray excess in clumpy galaxies, which cannot be solely due to star formation and confirms the presence of AGN. The clumpy galaxies in our intermediate-redshift sample have properties typical of gas-rich disk galaxies rather than mergers, being in particular on the Main Sequence of star formation. This suggests that our findings apply to the physically-similar and numerous gas-rich unstable disks at z>1. Using the observed [OIII] and X-ray luminosities, we conservatively estimate that AGN hosted by clumpy disks have typical bolometric luminosities of the order of a few 10^43 erg/s, BH growth rates ~10^-2 Msun/yr, and that these AGN are substantially obscured in X-rays. This moderate-luminosity mode could be sufficient to provide a large fraction of today's BH mass over a couple of Gyr given that our observations suggest a high duty cycle (>10%), accretion bursts with higher luminosities being possible over shorter phases. The observed evolution of disk instabilities with mass and redshift could explain the simultaneous downsizing of star formation and of BH growth.

preprint2012arXiv

Evidence for a wide range of UV obscuration in z ~ 2 dusty galaxies from the GOODS-Herschel survey

Dusty galaxies at z ~ 2 span a wide range of relative brightness between rest-frame mid-infrared (8um) and ultraviolet wavelengths. We attempt to determine the physical mechanism responsible for this diversity. Dust-obscured galaxies (DOGs), which have rest-frame mid-IR to UV flux density ratios > 1000, might be abnormally bright in the mid-IR, perhaps due to prominent AGN and/or PAH emission, or abnormally faint in the UV. We use far-infrared data from the GOODS-Herschel survey to show that most DOGs with 10^12 L_Sun < L_IR < 10^13 L_Sun are not abnormally bright in the mid-IR when compared to other dusty galaxies with similar IR (8--1000um) luminosities. We observe a relation between the median IR to UV luminosity ratios and the median UV continuum power-law indices for these galaxies, and we find that only 24% have specific star formation rates which indicate the dominance of compact star-forming regions. This circumstantial evidence supports the idea that the UV- and IR-emitting regions in these galaxies are spatially coincident, which implies a connection between the abnormal UV faintness of DOGs and dust obscuration. We conclude that the range in rest-frame mid-IR to UV flux density ratios spanned by dusty galaxies at z ~ 2 is due to differing amounts of UV obscuration. Of galaxies with these IR luminosities, DOGs are the most obscured. We attribute differences in UV obscuration to either: 1) differences in the degree of alignment between the spatial distributions of dust and massive stars, or 2) differences in the total dust content.

preprint2012arXiv

GMASS ultradeep spectroscopy of galaxies at z~2 - VII. Sample selection and spectroscopy

The populations of both quiescent and actively star-forming galaxies at 1<z<2 are still under-represented in our spectroscopic census of galaxies throughout the history of the Universe. In the light of galaxy formation models, however, the evolution of galaxies at these redshifts is of pivotal importance and merits further investigation. We therefore designed a spectroscopic observing campaign of a sample of both massive, quiescent and star-forming galaxies at z>1.4, called Galaxy Mass Assembly ultra-deep Spectroscopic Survey (GMASS). To determine redshifts and physical properties, such as metallicity, dust content, dynamical masses, and star formation history, we performed ultra-deep spectroscopy with the red-sensitive optical spectrograph FORS2 at the VLT. Our sample consists of objects, within the CDFS/GOODS area, detected at 4.5 micron, to be sensitive to stellar mass rather than star formation intensity. The spectroscopic targets were selected with a photometric redshift constraint (z>1.4) and magnitude constraints (B(AB)<26, I(AB)<26.5), which should ensure that these are faint, distant, and fairly massive galaxies. We present the sample selection, survey design, observations, data reduction, and spectroscopic redshifts. Up to 30 hours of spectroscopy of 174 spectroscopic targets and 70 additional objects enabled us to determine 210 redshifts, of which 145 are at z>1.4. From the redshifts and photometry, we deduce that the BzK selection criteria are efficient (82%) and suffer low contamination (11%). Several papers based on the GMASS survey show its value for studies of galaxy formation and evolution. We publicly release the redshifts and reduced spectra. In combination with existing and on-going additional observations in CDFS/GOODS, this data set provides a legacy for future studies of distant galaxies.

preprint2012arXiv

The early early type: discovery of a passive galaxy at z=3

We present the discovery of a massive, quiescent galaxy at z=2.99. We have obtained a HST/WFC3 spectrum of this object and measured its redshift from the detection of a deep 4000A break consistent with an old population and a high metallicity. By stellar population modeling of both its grism spectrum and broad-band photometry, we derive an age of ~0.7 Gyr, implying a formation redshift of z>4, and a mass >10^11 Msun. Although this passive galaxy is the most distant confirmed so far, we find that it is slightly less compact than other z>2 early-types of similar mass, being overall more analogous to those z~1.6 field early-type galaxies. The discovery of this object shows that early-type galaxies are detectable to at least z=3 and suggests that the diversity of structural properties found in z=1.4-2 ellipticals to earlier epochs could have its origin in a variety of formation histories among their progenitors.

preprint2011arXiv

Black Hole growth and AGN obscuration by instability-driven inflows in high-redshift disk galaxies fed by cold streams

Disk galaxies at high redshift have been predicted to maintain high gas surface densities due to continuous feeding by intense cold streams leading to violent gravitational instability, transient features and giant clumps. Gravitational torques between the perturbations drive angular momentum out and mass in, and the inflow provides the energy for keeping strong turbulence. We use analytic estimates of the inflow for a self-regulated unstable disk at a Toomre stability parameter Q~1, and isolated galaxy simulations capable of resolving the nuclear inflow down to the central parsec. We predict an average inflow rate ~10 Msun/yr through the disk of a 10^11 Msun galaxy, with conditions representative of z~2 stream-fed disks. The inflow rate scales with disk mass and (1+z)^{3/2}. It includes clump migration and inflow of the smoother component, valid even if clumps disrupt. This inflow grows the bulge, while only a fraction ~ 10^-3 of it needs to accrete onto a central black hole (BH), in order to obey the observed BH-bulge relation. A galaxy of 10^11 Msun at z~2 is expected to host a BH of ~10^8 Msun, accreting on average with moderate sub-Eddington luminosity L_X ~ 10^42-43 erg/s, accompanied by brighter episodes when dense clumps coalesce. We note that in rare massive galaxies at z~6, the same process may feed 10^9 Msun BH at the Eddington rate. High central gas column densities can severely obscure AGN in high-redshift disks, possibly hindering their detection in deep X-ray surveys.

preprint2011arXiv

CANDELS: The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey

The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) is designed to document the first third of galactic evolution, over the approximate redshift (z) range 8--1.5. It will image >250,000 distant galaxies using three separate cameras on the Hubble Space Telescope, from the mid-ultraviolet to the near-infrared, and will find and measure Type Ia supernovae at z>1.5 to test their accuracy as standardizable candles for cosmology. Five premier multi-wavelength sky regions are selected, each with extensive ancillary data. The use of five widely separated fields mitigates cosmic variance and yields statistically robust and complete samples of galaxies down to a stellar mass of 10^9 M_\odot to z \approx 2, reaching the knee of the ultraviolet luminosity function (UVLF) of galaxies to z \approx 8. The survey covers approximately 800 arcmin^2 and is divided into two parts. The CANDELS/Deep survey (5σ point-source limit H=27.7 mag) covers \sim 125 arcmin^2 within GOODS-N and GOODS-S. The CANDELS/Wide survey includes GOODS and three additional fields (EGS, COSMOS, and UDS) and covers the full area to a 5σ point-source limit of H \gtrsim 27.0 mag. Together with the Hubble Ultra Deep Fields, the strategy creates a three-tiered "wedding cake" approach that has proven efficient for extragalactic surveys. Data from the survey are nonproprietary and are useful for a wide variety of science investigations. In this paper, we describe the basic motivations for the survey, the CANDELS team science goals and the resulting observational requirements, the field selection and geometry, and the observing design. The Hubble data processing and products are described in a companion paper.

preprint2011arXiv

CANDELS: The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey - The Hubble Space Telescope Observations, Imaging Data Products and Mosaics

This paper describes the Hubble Space Telescope imaging data products and data reduction procedures for the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). This survey is designed to document the evolution of galaxies and black holes at $z\sim1.5-8$, and to study Type Ia SNe beyond $z>1.5$. Five premier multi-wavelength sky regions are selected, each with extensive multiwavelength observations. The primary CANDELS data consist of imaging obtained in the Wide Field Camera 3 / infrared channel (WFC3/IR) and UVIS channel, along with the Advanced Camera for Surveys (ACS). The CANDELS/Deep survey covers \sim125 square arcminutes within GOODS-N and GOODS-S, while the remainder consists of the CANDELS/Wide survey, achieving a total of \sim800 square arcminutes across GOODS and three additional fields (EGS, COSMOS, and UDS). We summarize the observational aspects of the survey as motivated by the scientific goals and present a detailed description of the data reduction procedures and products from the survey. Our data reduction methods utilize the most up to date calibration files and image combination procedures. We have paid special attention to correcting a range of instrumental effects, including CTE degradation for ACS, removal of electronic bias-striping present in ACS data after SM4, and persistence effects and other artifacts in WFC3/IR. For each field, we release mosaics for individual epochs and eventual mosaics containing data from all epochs combined, to facilitate photometric variability studies and the deepest possible photometry. A more detailed overview of the science goals and observational design of the survey are presented in a companion paper.

preprint2011arXiv

Molecular Gas at High Redshift

In order to understand galaxy evolution through cosmic times it is critical to derive the properties of the molecular gas content of galaxies, i.e. the material out of which stars ultimately form. The last decade has seen rapid progress in this area, with the detection of massive molecular gas reservoirs at high redshifts in submillimeter-selected galaxies and quasars. In the latter case, molecular gas reservoirs have been quantified out to redshifts z>6, i.e. towards the end of cosmic reionization when the universe was less than one Gyr old. The recent discovery of molecular gas in more normal galaxies have extended these studies from the most extreme objects in the universe (SFR~1000M_sun/yr; quasars and submillimeter galaxies) to more 'normal' starforming systems at redshifts 1.5-2.5 (with SFR~100M_sun/yr). However, detecting the molecular gas reservoirs of high-redshift galaxies that only have moderate star formation rates (~<10 M_sun/yr, similar to the faint galaxies seen in the Hubble Ultra Deep Field) will likely have to await the completion of ALMA.

preprint2011arXiv

On Star Formation Rates and Star Formation Histories of Galaxies out to z ~ 3

We compare multi-wavelength SFR indicators out to z~3 in GOODS-South. Our analysis uniquely combines U-to-8um photometry from FIREWORKS, MIPS 24um and PACS 70, 100, and 160um photometry from the PEP survey, and Ha spectroscopy from the SINS survey. We describe a set of conversions that lead to a continuity across SFR indicators. A luminosity-independent conversion from 24um to total infrared luminosity yields estimates of LIR that are in the median consistent with the LIR derived from PACS photometry, albeit with significant scatter. Dust correction methods perform well at low to intermediate levels of star formation. They fail to recover the total amount of star formation in systems with large SFR_IR/SFR_UV ratios, typically occuring at the highest SFRs (SFR_UV+IR \gtrsim 100 Msun/yr) and redshifts (z \gtrsim 2.5) probed. Finally, we confirm that Ha-based SFRs at 1.5<z<2.6 are consistent with SFR_SED and SFR_UV+IR provided extra attenuation towards HII regions is taken into account (Av,neb = Av,continuum / 0.44). With the cross-calibrated SFR indicators in hand, we perform a consistency check on the star formation histories inferred from SED modeling. We compare the observed SFR-M relations and mass functions at a range of redshifts to equivalents that are computed by evolving lower redshift galaxies backwards in time. We find evidence for underestimated stellar ages when no stringent constraints on formation epoch are applied. We demonstrate how resolved SED modeling, or alternatively deep UV data, may help to overcome this bias. The age bias is most severe for galaxies with young stellar populations, and reduces towards older systems. Finally, our analysis suggests that SFHs typically vary on timescales that are long (at least several 100 Myr) compared to the galaxies' dynamical time.

preprint2011arXiv

The zCOSMOS-SINFONI Project I: Sample Selection and Natural-Seeing Observations

The zCOSMOS SINFONI project is aimed at studying the physical and kinematical properties of a sample of massive z~1.4-2.5 star-forming galaxies, through SINFONI near-IR integral field spectroscopy (IFS), combined with the multi-wavelength information from the zCOSMOS (COSMOS) survey. The project is based on 1 hour of natural-seeing observations per target, and Adaptive Optics (AO) follow-up for a major part of the sample, which includes 30 galaxies selected from the zCOSMOS/VIMOS spectroscopic survey. This first paper presents the sample selection, and the global physical characterization of the target galaxies from multicolor photometry, i.e., star formation rate (SFR), stellar mass, age, etc. The Halpha integrated properties such as, flux, velocity dispersion, and size, are derived from the natural-seeing observations, while the follow up AO observations will be presented in the next paper of this series. Our sample appears to be well representative of star-forming galaxies at z~2, covering a wide range in mass and SFR. The Halpha integrated properties of the 25 Halpha detected galaxies are similar to those of other IFS samples at the same redshifts. Good agreement is found among the SFRs derived from Halpha luminosity and other diagnostic methods, provided the extinction affecting the Halpha luminosity is about twice that affecting the continuum. A preliminary kinematic analysis, based on the maximum observed velocity difference across the source, and on the integrated velocity dispersion, indicates that the sample splits nearly 50-50 into rotation-dominated and velocity dispersion-dominated galaxies, in good agreement with previous surveys.

preprint2010arXiv

A Multi-wavelength View of the Star Formation Activity at z~3

We present a multi-wavelength, UV-to-radio analysis for a sample of massive (M$_{\ast}$ $\sim$ 10$^{10}$ M$_\odot$) IRAC- and MIPS 24$μ$m-detected Lyman Break Galaxies (LBGs) with spectroscopic redshifts z$\sim$3 in the GOODS-North field (L$_{\rm UV}$$>1.8\times$L$^{\ast}_{z=3}$). For LBGs without individual 24$μ$m detections, we employ stacking techniques at 24$μ$m, 1.1mm and 1.4GHz, to construct the average UV-to-radio spectral energy distribution and find it to be consistent with that of a Luminous Infrared Galaxy (LIRG) with L$\rm_{IR}$=4.5$^{+1.1}_{-2.3}$$\times 10^{11}$ L$_{\odot}$ and a specific star formation rate (SSFR) of 4.3 Gyr$^{-1}$ that corresponds to a mass doubling time $\sim$230 Myrs. On the other hand, when considering the 24$μ$m-detected LBGs we find among them galaxies with L$\rm_{IR}> 10^{12}$ L$_{\odot}$, indicating that the space density of $z\sim$3 UV-selected Ultra-luminous Infrared Galaxies (ULIRGs) is $\sim$(1.5$\pm$0.5)$\times 10^{-5}$ Mpc$^{-3}$. We compare measurements of star formation rates (SFRs) from data at different wavelengths and find that there is tight correlation (Kendall's $τ>$ 99.7%) and excellent agreement between the values derived from dust-corrected UV, mid-IR, mm and radio data for the whole range of L$\rm_{IR}$ up to L$\rm_{IR}$ $\sim$ 10$^{13}$ L$_{\odot}$. This range is greater than that for which the correlation is known to hold at z$\sim$2, possibly due to the lack of significant contribution from PAHs to the 24$μ$m flux at $z\sim$3. The fact that this agreement is observed for galaxies with L$\rm_{IR}$ $>$ 10$^{12}$ L$_{\odot}$ suggests that star-formation in UV-selected ULIRGs, as well as the bulk of star-formation activity at this redshift, is not embedded in optically thick regions as seen in local ULIRGs and submillimeter-selected galaxies at $z=2$.

preprint2010arXiv

A Wide Area Survey for High-redshift Massive Galaxies. II. Near-infrared Spectroscopy of BzK-selected Massive Star-forming Galaxies

Results are presented from NIR spectroscopy of a sample of BzK-selected, massive star-forming galaxies (sBzKs) at 1.5<z<2.3 that were obtained with OHS/CISCO at Subaru and with SINFONI at VLT. Among the 28 sBzKs observed, Ha emission was detected in 14 objects, and for 11 of them the [NII]6583 was also measured. Multiwavelength photometry was also used to derive stellar masses and extinction parameters, whereas Ha and [NII] have allowed us to estimate SFR, metallicities, ionization mechanisms, and dynamical masses. In order to enforce agreement between SFRs from Ha with those derived from rest-frame UV and MIR, additional obscuration for the emission lines (that originate in HII regions) was required compared to the extinction derived from the slope of the UV continuum. We have also derived the stellar mass-metallicity relation, as well as the relation between stellar mass and specific SFR, and compared them to the results in other studies. At a given stellar mass, the sBzKs appear to have been already enriched to metallicities close to those of local star-forming galaxies of similar mass. The sBzKs presented here tend to have higher metallicities compared to those of UV-selected galaxies, indicating that NIR selected galaxies tend to be a chemically more evolved population. The sBzKs show specific SFRs that are systematically higher, by up to ~2 orders of magnitude, compared to those of local galaxies of the same mass. The empirical correlations between stellar mass and metallicity, and stellar mass and specific SFR are then compared with those of population synthesis models constructed either with the simple closed-box assumption, or within an infall scenario. Within the assumptions that are built-in such models, it appears that a short timescale for the star-formation (~100 Myr) and large initial gas mass appear to be required if one wants to reproduce both relations simultaneously.

preprint2010arXiv

Cold molecular gas in massive disk galaxies at z=1.5

We report the detection of the CO J=1-0 emission line in three near-infrared selected star-forming galaxies at z~1.5 with the Very Large Array (VLA) and the Green Bank telescope (GBT). These observations directly trace the bulk of molecular gas in these galaxies. We find H_2 gas masses of 8.3 \pm 1.9 x 10^{10} M_sun, 5.6 \pm 1.4 x 10^{10} M_sun and 1.23 \pm 0.34 x 10^{11} M_sun for BzK-4171, BzK-21000 and BzK-16000, respectively, assuming a conversion alpha_CO=3.6 M_sun (K km s^{-1} pc^{2})^{-1}. We combined our observations with previous CO 2-1 detections of these galaxies to study the properties of their molecular gas. We find brightness temperature ratios between the CO 2-1 and CO 1-0 emission lines of 0.80_{-0.22}^{+0.35}, 1.22_{-0.36}^{+0.61} and 0.41_{-0.13}^{+0.23} for BzK-4171, BzK-21000 and BzK-16000, respectively. At the depth of our observations it is not possible to discern between thermodynamic equilibrium or sub-thermal excitation of the molecular gas at J=2. However, the low temperature ratio found for BzK-16000 suggests sub-thermal excitation of CO already at J=2. For BzK-21000, a Large Velocity Gradient model of its CO emission confirms previous results of the low-excitation of the molecular gas at J=3. From a stacked map of the CO 1-0 images, we measure a CO 2-1 to CO 1-0 brightness temperature ratio of 0.92_{-0.19}^{+0.28}. This suggests that, on average, the gas in these galaxies is thermalized up to J=2, has star-formation efficiencies of ~100 L_sun (K km s^{-1} pc^2)^{-1} and gas consumption timescales of ~0.4 Gyr, unlike SMGs and QSOs at high redshifts.

preprint2009arXiv

Destruction of Molecular Gas Reservoirs in Early-Type Galaxies by Active Galactic Nucleus Feedback

Residual star formation at late times in early-type galaxies and their progenitors must be suppressed in order to explain the population of red, passively evolving systems we see today. Likewise, residual or newly accreted reservoirs of molecular gas that are fuelling star formation must be destroyed. This suppression of star formation in early-type galaxies is now commonly attributed to AGN feedback wherein the reservoir of gas is heated and expelled during a phase of accretion onto the central supermassive black hole. However, direct observational evidence for a link between the destruction of this molecular gas and an AGN phase has been missing so far. We present new mm-wavelength observations from the IRAM 30m telescope of a sample of low redshift SDSS early-type galaxies currently undergoing this process of quenching of late-time star formation. Our observations show that the disappearance of the molecular gas coincides within less than 100 Myr with the onset of accretion onto the black hole and is too rapid to be due to star formation alone. Since our sample galaxies are not associated to powerful quasar activity or radio jets, we conclude that low-luminosity AGN episodes are sufficient to suppress residual star formation in early-type galaxies. This `suppression mode' of AGN feedback is very different from the `truncation mode' linked to powerful quasar activity during early phases of galaxy formation.

preprint2009arXiv

Early-type galaxies in the PEARS survey: Probing the stellar populations at moderate redshift

Using HST/ACS slitless grism spectra from the PEARS program, we study the stellar populations of morphologically selected early-type galaxies in the GOODS North and South fields. The sample - extracted from a visual classification of the (v2.0) HST/ACS images and restricted to redshifts z>0.4 - comprises 228 galaxies (F775W<24 ABmag) out to z~1.3 over 320 arcmin2, with a median redshift zM=0.75. This work significantly increases our previous sample from the GRAPES survey in the HUDF (18 galaxies over ~11 arcmin2; Pasquali et al. 2006b). The grism data allow us to separate the sample into `red' and `blue' spectra, with the latter comprising 15% of the total. Three different grids of models parameterising the star formation history are used to fit the low-resolution spectra. Over the redshift range of the sample - corresponding to a cosmic age between 5 and 10 Gyr - we find a strong correlation between stellar mass and average age, whereas the **spread** of ages (defined by the RMS of the distribution) is roughly ~1 Gyr and independent of stellar mass. The best-fit parameters suggest it is formation epoch and not formation timescale, that best correlates with mass in early-type galaxies. This result - along with the recently observed lack of evolution of the number density of massive galaxies - motivates the need for a channel of (massive) galaxy formation bypassing any phase in the blue cloud, as suggested by the simulations of Dekel et al. (2009).

preprint2007arXiv

Redshifts of Emission Line Objects in the Hubble Ultra Deep Field

We present redshifts for 115 emission line objects in the Hubble Ultra Deep Field (HUDF) identified through the GRism ACS Program for Extragalactic Science (GRAPES) project using the slitless grism spectroscopy mode of the ACS Camera on the Hubble Space Telescope (HST). The sample was selected by an emission line search on all extracted 1-dimensional GRAPES spectra. We identify the emission lines using line wavelength ratios where multiple lines are detected in the grism wavelength range (5800A < lambda < 9600A), and using photometric redshift information where multiple lines are unavailable. We then derive redshifts using the identified lines. Our redshifts are accurate to delta(z) = 0.009, based on both statistical uncertainty estimates and comparison with published ground-based spectra. Over 40% of our sample is fainter than typical magnitude limits for ground-based spectroscopy (with i_{AB}>25 mag). Such emission lines would likely remain undiscovered without our deep survey. The emission line objects fall into 3 categories: 1) Most are low to moderate redshift galaxies (0 < z < 2), including many actively star forming galaxies with strong HII regions; 2) 9 are high redshift (4 < z < 7) Lyman-alpha emitters; and 3) at least 3 are candidate AGNs.

preprint2006arXiv

The Hubble Deep Field North SCUBA Super-map IV - Characterizing submillimetre galaxies using deep Spitzer imaging

We present SEDs, Spitzer colours, and IR luminosities for 850 micron selected galaxies in the GOODS-N field. Using the deep Spitzer Legacy images and new data and reductions of the VLA-HDF radio data, we find statistically secure counterparts for 60 per cent (21/35) of our submm sample, and identify tentative counterparts for another 12 objects. This is the largest sample of submm galaxies with statistically secure counterparts detected in the radio and with Spitzer. We find that in most cases the 850 micron emission is dominated by a single 24 micron source. A composite rest-frame SED shows that the submm sources peak at longer wavelengths than those of local ULIRGs of the same luminosity and therefore appear to be cooler. The SEDs of submm galaxies are also different from those of their high redshift neighbours, the near-IR selected BzK galaxies, whose mid-IR to radio SEDs are more like those of local ULIRGs. Using 24 micron. 850 micron and 1.4 GHz observations, we fit templates that span the mid-IR through radio to derived the integrated IR luminosity of the submm galaxies. By themselves, 24 micron and radio fluxes are able to predict LIR reasonably well because they are relatively insensitive to temperature. However, the submm flux by itself consistently overpredicts LIR when using spectral templates which obey the local ULIRG temperature-luminosity relation. The shorter Spitzer wavelengths sample the stellar bump at the redshifts of the submm sources, and we find that the Spitzer photometry alone provides a model independent estimate of the redshift. Using X-ray and mid-IR data, only 5 per cent of our secure counterparts show strong evidence for an active galactic nucleus dominating the IR luminosity.

preprint2004arXiv

Stellar Populations and Kinematics of Red Galaxies at z>2: Implications for the Formation of Massive Galaxies

We recently identified a substantial population of galaxies at z>2 with red rest-frame optical colors. These distant red galaxies (DRGs) are efficiently selected by the simple observed color criterion J-K>2.3. In this paper we present NIR spectroscopy with Keck/NIRSPEC of six DRGs at 2.4<z<3.2. We detect continuum emission and emission lines of all observed galaxies. Equivalent widths of H alpha are 20-30 Ang, smaller than measured for LBGs and nearby LIRGs, and comparable to normal nearby galaxies. The modest equivalent widths imply that the galaxies either have a decreasing star formation rate, or that they are very dusty. Fitting both the photometry and the H alpha lines, we find continuum extinction A_V=1-2 mag, ages 1-2.5 Gyr, star formation rates 200-400 solar masses/yr, and stellar masses 1-5x10^11 solar masses for models with constant star formation rates. From [NII]/H alpha ratios we infer that the metallicities are high, 1-1.5 x Solar. For four galaxies we can determine line widths from the optical emission lines. The widths are high, ranging from 130-240 km/s, and by combining data for LBGs and DRGs we find significant correlations between linewidth and restframe U-V color, and between linewidth and stellar mass. The latter correlation has a similar slope and offset as the ``baryonic Tully-Fisher relation'' for nearby galaxies. The median dynamical mass is ~2x10^11 solar masses, supporting the high stellar masses inferred from the photometry. We find that the median M/L_V ~ 0.8, a factor of ~5 higher than measured for LBGs. We infer from our small sample that DRGs are dustier, more metal rich, more massive, and have higher ages than z=3 LBGs of the same rest-frame V-band luminosity. Their high M/L ratios imply that they contribute significantly to the stellar mass density at z~2.5. [ABRIDGED]