Source author record

Roberto Maiolino

Roberto Maiolino 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

84works
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

84 published item(s)

preprint2026arXiv

Abundant Population of Broad H$α$ Emitters in the GOODS-N Field Revealed by CONGRESS, FRESCO, and JADES

We present a spectroscopic search for broad H$α$ emitters at z$\approx$3.7-6.5 in the GOODS-N field, utilizing JWST/NIRCam slitless spectroscopy from FRESCO and CONGRESS, complemented by JADES imaging. We identify 19 broad H$α$ emitters with FWHM$>$1000 km/s at z$\approx$4-5.5, including 9 new sources. The broad H$α$ luminosity function (LF) derived from our sample is consistent with those of other JWST-selected broad-line AGN reported in the literature. The black hole masses and AGN bolometric luminosities, inferred from the broad H$α$ components, indicate that most sources are accreting at ~10% of the Eddington limit. We derive their host stellar masses via SED fitting and find higher $M_{BH}/M_{*}$ ratios relative to the local $M_{BH}-M_{*}$ relations, consistent with previous studies. We find that 42% of the sample do not satisfy the widely-used color selection criteria for Little Red Dots (LRDs), with the majority of these sources lacking the characteristic steep rest-optical red slope, indicating that the LRD selection is highly incomplete when selecting AGN galaxies. A comparison of the average SEDs between our sample and LRDs selected in the same field reveals that the steep red slopes observed in some LRDs are likely due to line-boosting effects as previously suggested. Furthermore, we find that 68% of color-selected LRDs with H$α$ detections in the NIRCam/Grism spectra do not exhibit broad-line features. While the limited sensitivity of the grism spectra may hinder the detection of broad-line components in faint sources, our findings still highlight the enigmatic nature of the LRD population.

preprint2026arXiv

GA-NIFS: AGN activity in a Lyα emitter within a triple-AGN system anchored by a passive galaxy at z=3

Massive quenched galaxies at z>3 challenge models of early galaxy evolution, as their rapid formation and abrupt quenching require efficient feedback, often linked to active galactic nuclei (AGN). The quiescent galaxy GS10578 at z=3.1 is a key example of this population. Previous JWST/NIRSpec IFU data revealed an AGN outflow and uncovered a compact pair of AGN separated by ~5 kpc. In addition, VLT/MUSE spectroscopy has identified a third AGN candidate at a projected distance of ~30 kpc, associated with a luminous Lya emitter (LAE2) characterised by high-ionisation UV lines, although rest-frame optical diagnostics were not previously available. We aim to confirm the nature of LAE2 using rest-frame optical diagnostics enabled by new JWST data, and to characterise the physical and ionisation properties of both LAE2 and a distinct nearby Lya emitter (LAE1) that lacks any detectable continuum counterpart. We analyse new NIRSpec IFU observations targeting the optical nebular lines of LAE1 and LAE2, complemented with MUSE data, as part of the GA-NIFS project. We extract integrated and spatially resolved spectra, construct emission-line maps, and use standard diagnostic diagrams to determine ionisation sources and kinematics. LAE2 exhibits line ratios fully consistent with an embedded AGN. Its optical lines display a clumpy morphology and irregular kinematics on sub-kpc scales. Except for Lya, LAE1 remains undetected in all nebular lines and in JWST imaging. The similarity of the LAE1 and LAE2 Lya profiles in both velocity and flux suggests that LAE1 traces resonantly scattered emission rather than in-situ star formation. Our analysis reveals that the environment of GS10578 contains both multi-black-hole activity and gas structures on tens-of-kpc scales, offering new insights into how feedback and satellite interactions influence the late evolutionary stages of quenched massive galaxies.

preprint2026arXiv

JADES: the chemical enrichment pattern of distant galaxies -- $α$ enhancement, silicon depletion, and iron enhancement

We present gas-phase abundances of carbon (C), $α$-elements (O, Ne, Si, and Ar) and iron (Fe) obtained from stacked spectra of high-$z$ star-forming galaxies with the deep Near Infrared Spectrograph medium-resolution data from the James Webb Space Telescope Advanced Deep Extragalactic Survey. Our 564 sources at $z=4$--7 have a median stellar mass of $\log(M_{*}/M_{\odot})=8.46$ and a median star-formation rate of $\log(\mathrm{SFR}/M_{\odot}\,\mathrm{yr^{-1}})=0.30$, placing them close to the star-formation main sequence. We find that the stacked spectrum of all our 564 sources has relatively low [C/O]$=-0.70$, moderate [Ne/O]$=-0.09$, and low [Ar/O]$=-0.28$ values at a low gas-phase metallicity of $12+\log(\mathrm{O/H})=7.71$ ($Z\sim 0.1~Z_\odot$), suggesting dominant yields of core-collapse supernovae evolved from massive stars. The detection of a weak SiIII] emission line in our stacked spectrum provides a silicon-to-oxygen abundance ratio of [Si/O]$=-0.63$, which is lower than that of stars in the Milky Way disc and lower than expected by chemical evolution models, suggesting silicon depletion onto dust grains. Likewise, this Si/O value is lower than that we newly derive for two individual $z>6$ galaxies (GN-z11 and RXCJ2248) with negligible dust attenuation. By performing spectral stacking in bins of $M_{*}$, SFR, specific SFR (sSFR), and ultra-violet (UV) continuum slope $β_{\mathrm{UV}}$, we identify [FeIII] line detections in the high-sSFR bin and the blue-$β_{\mathrm{UV}}$ bin, both of which exhibit supersolar Fe/O ratios, while their C/O, Ar/O, and Si/O ratios are comparable to those of the all-sources stack. Our findings support a chemically young gas composition with rapid dust depletion in the general population of high-$z$ star-forming galaxies, while raising the possibility of anomalous, selective Fe/O enhancement at the very early epoch of star formation.

preprint2026arXiv

On the origins of oxygen: ALMA and JWST characterise the multi-phase, metal-enriched, star-bursting medium within a 'normal' $z > 11$ galaxy

The unexpectedly high abundance of galaxies at $z > 11$ revealed by JWST has sparked a debate on the nature of early galaxies and the physical mechanisms regulating their formation. The Atacama Large Millimeter/submillimeter Array (ALMA) has begun to provide vital insights on their gas and dust content, but so far only for extreme 'blue monsters'. Here we present new, deep ALMA observations of JADES-GS-z11-0, a more typical (sub-$L^*$) $z > 11$ galaxy that bridges the discovery space of JWST and the Hubble Space Telescope. These data confirm the presence of the [O III] 88 $μ$m line at $4.5σ$ significance, precisely at the redshift of several faint emission lines previously seen with JWST/NIRSpec, while the underlying dust continuum remains undetected ($F_ν< 9.0 \, \mathrm{μJy}$), implying an obscured star formation rate (SFR) of $\text{SFR}_\text{IR} \lesssim 6 \, \mathrm{M_\odot \, yr^{-1}}$ and dust mass of $M_\text{dust} \lesssim 1.0 \times 10^{6} \, \mathrm{M_\odot}$ (all $3σ$). The accurate ALMA redshift of $z_\text{[O III]} = 11.1221 \pm 0.0006$ ($\gtrsim \! 5\times$ refined over NIRSpec) helps confirm that redshifts measured purely from the Lyman-$α$ break, even spectroscopically, should properly take into account the effects of potential damped Lyman-$α$ absorption (DLA) systems to avoid systematic overestimates of up to $Δz \approx 0.5$. The [O III] 88 $μ$m luminosity of $L_\text{[O III]} = (1.1 \pm 0.3) \times 10^{8} \, \mathrm{L_\odot}$, meanwhile, agrees well with the scaling relation for local metal-poor dwarfs given the SFR measured by NIRCam, NIRSpec, and MIRI. The spatially resolved MIRI and ALMA emission also underscores that JADES-GS-z11-0 is likely to consist of two low-mass components that are undergoing strong bursts of star formation yet are already pre-enriched in oxygen (~20-30% solar), only 400 Myr after the Big Bang.

preprint2026arXiv

The Advanced X-ray Imaging Satellite (AXIS) Community Science Book

The AXIS Community Science Book represents the collective effort of 592 scientists worldwide to define the transformative science enabled by the Advanced X-ray Imaging Satellite (AXIS), a next-generation X-ray mission selected by NASA's Astrophysics Probe Program for Phase A study. AXIS will advance the legacy of high-angular-resolution X-ray astronomy with ~1.5'' imaging over a wide 24' field of view and an order of magnitude greater collecting area than Chandra in the 0.3-12 keV band. Combining sharp imaging, high throughput, and rapid response capabilities, AXIS will open new windows on virtually every aspect of modern astrophysics, exploring the birth and growth of supermassive black holes, the feedback processes that shape galaxies, the life cycles of stars and exoplanet environments, and the nature of compact stellar remnants, supernova remnants, and explosive transients. This book compiles 138 community-contributed science cases developed by five Science Working Groups focused on AGN and supermassive black holes, galaxy evolution and feedback, compact objects and supernova remnants, stellar physics and exoplanets, and time-domain and multi-messenger astrophysics. Together, these studies establish the scientific foundation for next-generation X-ray exploration in the 2030s and highlight strong synergies with facilities of the 2030s, such as JWST, Roman, Rubin/LSST, SKA, ALMA, ngVLA, and next-generation gravitational-wave and neutrino networks.

preprint2024arXiv

Inside the bubble: exploring the environments of reionisation-era Lyman-$α$ emitting galaxies with JADES and FRESCO

We present a study of the environments of 17 Lyman-$α$ (Ly$α$) emitting galaxies (LAEs) in the reionisation era ($5.8 < z < 8$) identified by JWST/NIRSpec as part of the JWST Advanced Deep Extragalactic Survey (JADES). Unless situated in sufficiently (re)ionised regions, Ly$α$ emission from these galaxies would be strongly absorbed by neutral gas in the intergalactic medium (IGM). We conservatively estimate sizes of the ionised regions required to reconcile the relatively low Ly$α$ velocity offsets ($Δv_\text{Ly$α$}<300\,\mathrm{km\,s^{-1}}$) with moderately high Ly$α$ escape fractions ($f_\mathrm{esc,\,Lyα}>5\%$) observed in our sample of LAEs, suggesting the presence of ionised hydrogen along the line of sight towards at least eight out of 17 LAEs. We find minimum physical `bubble' sizes of the order of $R_\text{ion}\sim0.1$-$1\,\mathrm{pMpc}$ are required in a patchy reionisation scenario where ionised bubbles containing the LAEs are embedded in a fully neutral IGM. Around half of the LAEs in our sample are found to coincide with large-scale galaxy overdensities seen in FRESCO at $z \sim 5.8$-$5.9$ and $z\sim7.3$, suggesting Ly$α$ transmission is strongly enhanced in such overdense regions, and underlining the importance of LAEs as tracers of the first large-scale ionised bubbles. Considering only spectroscopically confirmed galaxies, we find our sample of UV-faint LAEs ($M_\text{UV}\gtrsim-20\,\mathrm{mag}$) and their direct neighbours are generally not able to produce the required ionised regions based on the Ly$α$ transmission properties, suggesting lower-luminosity sources likely play an important role in carving out these bubbles. These observations demonstrate the combined power of JWST multi-object and slitless spectroscopy in acquiring a unique view of the early Universe during cosmic reionisation via the most distant LAEs.

preprint2022arXiv

A Fresh Look at AGN Spectral Energy Distribution Fitting with the XMM-SERVS AGN Sample

We perform spectral energy distribution (SED) fitting to 711 luminous X-ray AGN at 0.7 < z < 4.5 using 10-bands of optical and infra-red photometric data for objects within XMM-SERVS. This fitting provided 510 reliable (reduced $χ^2 < 3$) inferences on AGN and host galaxy properties. The AGN optical (3000Å) luminosity inferred from SED-fitting is found to correlate with the measured X-ray (2-10 keV) luminosity, in good agreement with previous work. Using X-ray hardness as a proxy for AGN obscuration, we also study the differences in the host galaxy properties of obscured and unobscured AGN. Both populations have consistent stellar masses (log$_{10}(M_*/M_{\odot})$ = 10.88 $\pm0.09M_\odot$ and log$_{10}(M_*/M_{\odot})$ = 10.8 $\pm0.1M_\odot$ for unobscured and obscured AGN respectively). We also find evidence for varying AGN emission line properties from a standard AGN template in 18.8% of the sample with a reduced $χ^2 < 3$ where the inclusion of an additional emission line strength free parameter was found to improve the quality of the fit. Comparison of these fits to SDSS spectra showed that emission line properties inferred from broadband photometry were consistent with the results from spectroscopy for 91% of objects. We find that the presence of weaker, more blueshifted emission lines as inferred from the SED fits are associated with more negative values of $α_{ox}$. While the correlation between the hardness of the ionising SED and the emission line properties has been known for some time, we are able to derive this correlation purely from broadband photometry.

preprint2022arXiv

ALMaQUEST -- VII: Star Formation Scaling Relations of Green Valley Galaxies

We utilize the ALMA-MaNGA QUEnch and STar formation (ALMaQUEST) survey to investigate the kpc-scale scaling relations, presented as the resolved star forming main sequence (rSFMS: $Σ_{\rm SFR}$ vs. $Σ_{*}$), the resolved Schmidt-Kennicutt relation (rSK: $Σ_{\rm SFR}$ vs. $Σ_{\rm H_{2}}$), and the resolved molecular gas main sequence (rMGMS: $Σ_{\rm H_{2}}$ vs. $Σ_{*}$), for 11478 star-forming and 1414 retired spaxels (oversampled by a factor of $\sim20$) located in 22 green valley (GV) and 12 main sequence (MS) galaxies. For a given galaxy type (MS or GV), the retired spaxels are found to be offset from the sequences formed by the star-forming spaxels on the rSFMS, rSK, and rMGMS planes, toward lower absolute values of sSFR, SFE, and $f_{\rm H_{2}}$ by $\sim$ 1.1, 0.6, and 0.5 dex. The scaling relations for GV galaxies are found to be distinct from that of the MS galaxies, even if the analyses are restricted to the star-forming spaxels only. It is found that for star-forming spaxels, sSFR, SFE, and $f_{\rm H_{2}}$ in GV galaxies are reduced by $\sim$0.36, 0.14, and 0.21 dex, respectively, compared to those in MS galaxies. Therefore, the suppressed sSFR/SFE/$f_{\rm H_{2}}$ in GV galaxies are associated with not only an increased proportion of retired regions in GV galaxies but also a depletion of these quantities in star-forming regions. Finally, the reduction of SFE and $f_{\rm H_{2}}$ in GV galaxies relative to MS galaxies is seen in both bulge and disk regions (albeit with larger uncertainties), suggesting that statistically, quenching in the GV population may persist from the inner to the outer regions.

preprint2022arXiv

Being KLEVER at cosmic noon: ionised gas outflows are inconspicuous in low-mass star-forming galaxies but prominent in massive AGN hosts

We investigate the presence of ionised gas outflows in a sample of 141 main-sequence star-forming galaxies at $1.2<z<2.6$ from the KLEVER (KMOS Lensed Emission Lines and VElocity Review) survey. Our sample covers an exceptionally wide range of stellar masses, $8.1<\log(M_\star/M_{\odot})<11.3$, pushing outflow studies into the dwarf regime thanks to gravitationally lensed objects. We stack optical rest-frame emission lines (H$β$, [OIII], H$α$ and [NII]) in different mass bins and seek for tracers of gas outflows by using a novel, physically motivated method that improves over the widely used, simplistic double Gaussian fitting. We compare the observed emission lines with the expectations from a rotating disc (disc+bulge for the most massive galaxies) model, whereby significant deviations are interpreted as a signature of outflows. We find clear evidence for outflows in the most massive, $\log(M_\star/M_{\odot}) > 10.8$, AGN-dominated galaxies, suggesting that AGNs may be the primary drivers of these gas flows. Surprisingly, at $\log(M_\star/M_{\odot})\leq 9.6$, the observed line profiles are fully consistent with a rotating disc model, indicating that ionised gas outflows in dwarf galaxies might play a negligible role even during the peak of cosmic star-formation activity. Finally, we find that the observed mass loading factor scales with stellar mass as expected from the TNG50 cosmological simulation, but the ionised gas mass accounts for only 2$\%$ of the predicted value. This suggests that either the bulk of the outflowing mass is in other gaseous phases or the current feedback models implemented in cosmological simulations need to be revised.

preprint2022arXiv

Cold Gas in Massive Galaxies as A Critical Test of Black Hole Feedback Models

Black hole feedback has been widely implemented as the key recipe to quench star formation in massive galaxies in modern semi-analytic models and hydrodynamical simulations. As the theoretical details surrounding the accretion and feedback of black holes continue to be refined, various feedback models have been implemented across simulations, with notable differences in their outcomes. Yet, most of these simulations have successfully reproduced some observations, such as stellar mass function and star formation rate density in the local Universe. We use the recent observation on the change of neutral hydrogen gas mass (including both ${\rm H_2}$ and ${\rm HI}$) with star formation rate of massive central disc galaxies as a critical constraint of black hole feedback models across several simulations. We find that the predictions of IllustrisTNG agree with the observations much better than the other models tested in this work. This favors IllustrisTNG's treatment of active galactic nuclei - where kinetic winds are driven by black holes at low accretion rates - as more plausible amongst those we test. In turn, this also indirectly supports the idea that the massive central disc galaxy population in the local Universe was likely quenched by AGN feedback.

preprint2022arXiv

Dual constraints with ALMA: new [O III] 88 ${\rm μ}$m and dust-continuum observations reveal the ISM conditions of luminous LBGs at $z \sim 7$

We present new [O III] 88 ${\rm μ}$m observations of five bright $z \sim 7$ Lyman-break galaxies spectroscopically confirmed by ALMA through the [C II] 158 ${\rm μ}$m line, unlike recent [O III] detections where Lyman-${\rm α}$ was used. This nearly doubles the sample of Epoch of Reionisation galaxies with robust ($5 σ$) detections of [C II] and [O III]. We perform a multi-wavelength comparison with new deep HST images of the rest-frame UV, whose compact morphology aligns well with [O III] tracing ionised gas. By contrast, we find more spatially extended [C II] emission likely produced in neutral gas, as indicated by a [N II] 205 ${\rm μ}$m non-detection in one source. We find a positive correlation between the equivalent width of the optical [O III] and H${\rm β}$ lines and the [O III]/[C II] ratio, as seen in local metal-poor dwarf galaxies. Cloudy models of a nebula of typical density harbouring a young stellar population with a high ionisation parameter appear to adequately reproduce the far-infrared lines. Surprisingly, however, our models fail to reproduce the strength of [O III] 88 ${\rm μ}$m, unless we assume an ${\rm α}$/Fe enhancement and a near-solar nebular oxygen abundance. On spatially resolved scales, we find [O III]/[C II] shows a tentative anti-correlation with infrared excess, $L_{\rm IR}/L_{\rm UV}$, also seen on global scales in the local Universe. Finally, we introduce the far-infrared spectral energy distribution fitting code MERCURIUS to show that dust-continuum measurements of one source appear to favour a low dust temperature coupled with a high dust mass. This implies a high stellar metallicity yield and may point towards the need of dust production or grain-growth mechanisms beyond supernovae.

preprint2022arXiv

Radio jet-ISM interaction and positive radio-mechanical feedback in Abell 1795

We present XSHOOTER observations with previous ALMA, MUSE and $HST$ observations to study the nature of radio-jet triggered star formation and the interaction of radio jets with the interstellar medium in the brightest cluster galaxy (BCG) in the Abell 1795 cluster. Using $HST$ UV data we determined an ongoing star formation rate of 9.3 M$_\odot$ yr$^{-1}$. The star formation follows the global Kennicutt-Schmidt law, however, it has a low efficiency compared to circumnuclear starbursts in nearby galaxies with an average depletion time of $\sim$1 Gyr. The star formation and molecular gas are offset by $\sim1$ kpc indicating that stars have decoupled from the gas. We detected an arc of high linewidth in ionized gas where electron densities are elevated by a factor of $\sim$4 suggesting a shock front driven by radio jets or peculiar motion of the BCG. An analysis of nebular emission line flux ratios suggests that the gas is predominantly ionized by star formation with a small contribution from shocks. We also calculated the velocity structure function (VSF) of the ionized and molecular gases using velocity maps to characterize turbulent motion in the gas. The ionized gas VSF suggests that the radio jets are driving supersonic turbulence in the gas. Thus radio jets can not only heat the atmosphere on large scales and may quench star formation on longer timescales while triggering star formation in positive feedback on short timescales of a few million years.

preprint2022arXiv

Seeking the growth of the first black hole seeds with JWST

In this paper we provide predictions for the BH population that would be observable with planned JWST surveys at $5 \le z \le 15$. We base our study on the recently developed Cosmic Archaeology Tool (CAT), which allows us to model BH seeds formation and growth, while being consistent with the general population of AGNs and galaxies observed at $4 \le z \le 7$. We find that JWST planned surveys will provide a complementary view on active BHs at $z > 5$, with JADES-Medium/-Deep being capable of detecting the numerous BHs that populate the faint-end of the distribution, COSMOS-Web sampling a large enough area to detect the rarest brightest systems, and CEERS/PRIMER bridging the gap between these two regimes. The relatively small field of view of the above surveys preferentially selects BHs with masses $6 \leq {\rm Log} (M_{\rm BH}/M_\odot) < 8$ at $7 \le z < 10$, residing in relatively metal poor (${\rm Log} (Z/Z_\odot) \ge -2$) and massive ($8\leq {\rm Log} (M_*/M_\odot) < 10$) galaxies. At $z \ge 10$, only JADES-Deep will have the sensitivity to detect growing BHs with masses $4 \leq {\rm Log} (M_{\rm BH}/M_\odot) < 6$, hosted by more metal poor ($-3 \leq {\rm Log} (Z/Z_\odot) < -2$) and less massive ($6 \leq {\rm Log} (M_*/M_\odot) < 8$) galaxies. In our model, the latter population corresponds to heavy BH seeds formed by the direct collapse of super-massive stars in their earliest phases of mass growth. Detecting these systems would provide invaluable insights on the nature and early growth of the first BH seeds.

preprint2022arXiv

The ALMaQUEST Survey IX: The nature of the resolved star forming main sequence

We investigate the nature of the scaling relations between the surface density of star formation rate ($Σ_{\rm SFR}$), stellar mass ($Σ_*$), and molecular gas mass ($Σ_{\rm H_2}$), aiming at distinguishing between the relations that are primary, i.e. more fundamental, and those which are instead an indirect by-product of the other relations. We use the ALMaQUEST survey and analyse the data by using both partial correlations and Random Forest regression techniques. We unambiguously find that the strongest intrinsic correlation is between $Σ_{\rm SFR}$ and $Σ_{\rm H_2}$ (i.e. the resolved Schmidt-Kennicutt relation), followed by the correlation between $Σ_{\rm H_2}$ and $Σ_*$ (resolved Molecular Gas Main Sequence, rMGMS). Once these two correlations are taken into account, we find that there is no evidence for any intrinsic correlation between $Σ_{\rm SFR}$ and $Σ_*$, implying that SFR is entirely driven by the amount of molecular gas, while its dependence on stellar mass (i.e. the resolved Star Forming Main Sequence, rSFMS) simply emerges as a consequence of the relationship between molecular gas and stellar mass.

preprint2022arXiv

The interstellar medium of high-redshift galaxies: Gathering clues from C III] and [C II] lines

A tight relation between [C II] line luminosity and the star formation rate (SFR) has been observed for local galaxies. At high redshift (z > 5), galaxies instead deviate downwards from the local $Σ$_[C II] - $Σ$_SFR relation. This deviation might be caused by different interstellar medium (ISM) properties in galaxies at early epochs. To test this hypothesis, we combined the [C II] and SFR data with C III] line observations and our physical models. We additionally investigated how ISM properties, such as burstiness, $κ_s$, total gas density, $n$, and metallicity, $Z$, affect the deviation from the $Σ$_[C II] - $Σ$_SFR relation in these sources. We present the VLT/X-SHOOTER observations targeting the C III] $λ1909$ line emission in three galaxies at 5.5 < z < 7. We include X-SHOOTER and VLT/MUSE archival data of eight galaxies at 2 < z < 7, and eleven star-forming systems at 6 < z < 7.5, with either C III] or [C II] detection reported in the literature. We detected C III] $λλ1907, 1909$ line emission in HZ10 and we derived the intrinsic, integrated flux of the C III] $λ1909$ line. We constrained the ISM properties for our sample of galaxies, $κ_s$ , $n$, and $Z$, by applying our physically motivated model based on the MCMC algorithm. For the most part, high-z star-forming galaxies show subsolar metallicities. The majority of the sources have $log(κ_s) > 1$, that is, they overshoot the Kennicutt-Schmidt (KS) relation by about one order of magnitude, implying that the whole KS relation might be shifted upwards at early times. Furthermore, all the high-z galaxies of our sample lie below the $Σ$_[C II] - $Σ$_SFR local relation. The total gas density, $n$, shows the strongest correlation with the deviation from the local $Σ$_[C II] - $Σ$_SFR relation.

preprint2022arXiv

The KLEVER survey: Nitrogen abundances at $z\sim$2 and probing the existence of a fundamental nitrogen relation

We present a comparison of the nitrogen-to-oxygen ratio (N/O) in 37 high-redshift galaxies at $z\sim$2 taken from the KMOS Lensed Emission Lines and VElocity Review (KLEVER) Survey with a comparison sample of local galaxies, taken from the Sloan Digital Sky Survey (SDSS). The KLEVER sample shows only a mild enrichment in N/O of $+$0.1 dex when compared to local galaxies at a given gas-phase metallicity (O/H), but shows a depletion in N/O of $-$0.36 dex when compared at a fixed stellar mass (M$_*$). We find a strong anti-correlation in local galaxies between N/O and SFR in the M$_*$-N/O plane, similar to the anti-correlation between O/H and SFR found in the mass-metallicity relation (MZR). We use this anti-correlation to construct a fundamental nitrogen relation (FNR), analogous to the fundamental metallicity relation (FMR). We find that KLEVER galaxies are consistent with both the FMR and the FNR. This suggests that the depletion of N/O in high-$z$ galaxies when considered at a fixed M$_*$ is driven by the redshift-evolution of the mass-metallicity relation in combination with a near redshift-invariant N/O-O/H relation. Furthermore, the existence of an fundamental nitrogen relation suggests that the mechanisms governing the fundamental metallicity relation must be probed by not only O/H, but also N/O, suggesting pure-pristine gas inflows are not the primary driver of the FMR, and other properties such as variations in galaxy age and star formation efficiency must be important.

preprint2022arXiv

The quenching of galaxies, bulges, and disks since cosmic noon: A machine learning approach for identifying causality in astronomical data

We present an analysis of the quenching of star formation in galaxies, bulges, and disks throughout the bulk of cosmic history, from $z=2-0$. We utilise observations from the SDSS and MaNGA at low redshifts. We complement these data with observations from CANDELS at high redshifts. Additionally, we compare the observations to detailed predictions from the LGalaxies semi-analytic model. To analyse the data, we developed a machine learning approach utilising a Random Forest classifier. We first demonstrate that this technique is extremely effective at extracting causal insight from highly complex and inter-correlated model data, before applying it to various observational surveys. Our primary observational results are as follows: At all redshifts studied in this work, we find bulge mass to be the most predictive parameter of quenching, out of the photometric parameter set (incorporating bulge mass, disk mass, total stellar mass, and $B/T$ structure). Moreover, we also find bulge mass to be the most predictive parameter of quenching in both bulge and disk structures, treated separately. Hence, intrinsic galaxy quenching must be due to a stable mechanism operating over cosmic time, and the same quenching mechanism must be effective in both bulge and disk regions. Despite the success of bulge mass in predicting quenching, we find that central velocity dispersion is even more predictive (when available in spectroscopic data sets). In comparison to the LGalaxies model, we find that all of these observational results may be consistently explained through quenching via preventative `radio-mode' active galactic nucleus (AGN) feedback. Furthermore, many alternative quenching mechanisms (including virial shocks, supernova feedback, and morphological stabilisation) are found to be inconsistent with our observational results and those from the literature.

preprint2022arXiv

What drives galaxy quenching? A deep connection between galaxy kinematics and quenching in the local Universe

We develop a 2D inclined rotating disc model, which we apply to the stellar velocity maps of 1862 galaxies taken from the MaNGA survey (SDSS public Data Release 15). We use a random forest classifier to identify the kinematic parameters that are most connected to galaxy quenching. We find that kinematic parameters that relate predominantly to the disc (such as the mean rotational velocity) and parameters that characterise whether a galaxy is rotation- or dispersion-dominated (such as the ratio of rotational velocity to velocity dispersion) are not fundamentally linked to the quenching of star formation. Instead, we find overwhelmingly that it is the absolute level of velocity dispersion (a property that relates primarily to a galaxy's bulge/spheroidal component) that is most important for separating star forming and quenched galaxies. Furthermore, a partial correlation analysis shows that many commonly discussed correlations between galaxy properties and quenching are spurious, and that the fundamental correlation is between quenching and velocity dispersion. In particular, we find that at fixed velocity dispersion, there is only a very weak dependence of quenching on the disc properties, whereby more discy galaxies are slightly more likely to be forming stars. By invoking the tight relationship between black hole mass and velocity dispersion, and noting that black hole mass traces the total energy released by AGN, we argue that these data support a scenario in which quenching occurs by preventive feedback from AGN. The kinematic measurements from this work are publicly available.

preprint2022arXiv

What drives the scatter of local star-forming galaxies in the BPT diagrams? A Machine Learning based analysis

We investigate which physical properties are most predictive of the position of local star forming galaxies on the BPT diagrams, by means of different Machine Learning (ML) algorithms. Exploiting the large statistics from the Sloan Digital Sky Survey (SDSS), we define a framework in which the deviation of star-forming galaxies from their median sequence can be described in terms of the relative variations in a variety of observational parameters. We train artificial neural networks (ANN) and random forest (RF) trees to predict whether galaxies are offset above or below the sequence (via classification), and to estimate the exact magnitude of the offset itself (via regression). We find, with high significance, that parameters primarily associated to variations in the nitrogen-over-oxygen abundance ratio (N/O) are the most predictive for the [N II]-BPT diagram, whereas properties related to star formation (like variations in SFR or EW[H$α$]) perform better in the [S II]-BPT diagram. We interpret the former as a reflection of the N/O-O/H relationship for local galaxies, while the latter as primarily tracing the variation in the effective size of the S$^{+}$ emitting region, which directly impacts the [S II]emission lines. This analysis paves the way to assess to what extent the physics shaping local BPT diagrams is also responsible for the offsets seen in high redshift galaxies or, instead, whether a different framework or even different mechanisms need to be invoked.

preprint2021arXiv

A massive stellar bulge in a regularly rotating galaxy 1.2 billion years after the Big Bang

Cosmological models predict that galaxies forming in the early Universe experience a chaotic phase of gas accretion and star formation, followed by gas ejection due to feedback processes. Galaxy bulges may assemble later via mergers or internal evolution. Here we present submillimeter observations (with spatial resolution of 700 parsecs) of ALESS 073.1, a starburst galaxy at redshift z~5, when the Universe was 1.2 billion years old. This galaxy's cold gas forms a regularly rotating disk with negligible noncircular motions. The galaxy rotation curve requires the presence of a central bulge in addition to a star-forming disk. We conclude that massive bulges and regularly rotating disks can form more rapidly in the early Universe than predicted by models of galaxy formation.

preprint2021arXiv

Galactic Winds across the Gas-Rich Merger Sequence: I. Highly Ionized N V and O VI Outflows in the QUEST Quasars

This program is part of QUEST (Quasar/ULIRG Evolutionary Study) and seeks to examine the gaseous environments of z < 0.3 quasars and ULIRGs as a function of host galaxy properties and age across the merger sequence from ULIRGs to quasars. This first paper in the series focuses on 33 quasars from the QUEST sample and on the kinematics of the highly ionized gas phase traced by the N V 1238, 1243 A and O VI 1032, 1038 A absorption lines in high-quality Hubble Space Telescope (HST) Cosmic Origins Spectrograph (COS) data. N V and O VI outflows are present in about 60% of the QUEST quasars and span a broad range of properties, both in terms of equivalent widths (from 20 mA to 25 A) and kinematics (outflow velocities from a few x 100 km/s up to ~10,000 km/s). The rate of incidence and equivalent widths of the highly ionized outflows are higher among X-ray weak or absorbed sources. The weighted outflow velocity dispersions are highest among the X-ray weakest sources. No significant trends are found between the weighted outflow velocities and the properties of the quasars and host galaxies although this may be due to the limited dynamic range of properties of the current sample. These results will be re-examined in an upcoming paper where the sample is expanded to include the QUEST ULIRGs. Finally, a lower limit of ~0.1% on the ratio of time-averaged kinetic power to bolometric luminosity is estimated in the 2-4 objects with blueshifted P V 1117, 1128 absorption features.

preprint2021arXiv

On the quenching of star formation in observed and simulated central galaxies: Evidence for the role of integrated AGN feedback

In this paper we investigate how massive central galaxies cease their star formation by comparing theoretical predictions from cosmological simulations: EAGLE, Illustris and IllustrisTNG with observations of the local Universe from the Sloan Digital Sky Survey (SDSS). Our machine learning (ML) classification reveals supermassive black hole mass ($M_{\rm BH}$) as the most predictive parameter in determining whether a galaxy is star forming or quenched at redshift $z=0$ in all three simulations. This predicted consequence of active galactic nucleus (AGN) quenching is reflected in the observations, where it is true for a range of indirect estimates of $M_{\rm BH}$ via proxies as well as its dynamical measurements. Our partial correlation analysis shows that other galactic parameters lose their strong association with quiescence, once their correlations with $M_{\rm BH}$ are accounted for. In simulations we demonstrate that it is the integrated power output of the AGN, rather than its instantaneous activity, which causes galaxies to quench. Finally, we analyse the change in molecular gas content of galaxies from star forming to passive populations. We find that both gas fractions ($f_{\rm gas}$) and star formation efficiencies (SFEs) decrease upon transition to quiescence in the observations but SFE is more predictive than $f_{\rm gas}$ in the ML passive/star-forming classification. These trends in the SDSS are most closely recovered in IllustrisTNG and are in direct contrast with the predictions made by Illustris. We conclude that a viable AGN feedback prescription can be achieved by a combination of preventative feedback and turbulence injection which together quench star formation in central galaxies.

preprint2021arXiv

Significant Dust-Obscured Star Formation in Luminous Lyman-break Galaxies at $z$$\sim$$7$$-$$8$

We make use of ALMA continuum observations of $15$ luminous Lyman-break galaxies at $z$$\sim$$7$$-$$8$ to probe their dust-obscured star-formation. These observations are sensitive enough to probe to obscured SFRs of $20$ $M_{\odot}$$/$$yr$ ($3σ$). Six of the targeted galaxies show significant ($\geq$$3$$σ$) dust continuum detections, more than doubling the number of known dust-detected galaxies at $z$$>$$6.5$. Their IR luminosities range from $2.7$$\times$$10^{11}$ $L_{\odot}$ to $1.1$$\times$$10^{12}$ $L_{\odot}$, equivalent to obscured SFRs of $20$ to $105$ $M_{\odot}$$/$$yr$. We use our results to quantify the correlation of the infrared excess IRX on the UV-continuum slope $β_{UV}$ and stellar mass. Our results are most consistent with an SMC attenuation curve for intrinsic $UV$-slopes $β_{UV,intr}$ of $-2.63$ and most consistent with an attenuation curve in-between SMC and Calzetti for $β_{UV,intr}$ slopes of $-2.23$, assuming a dust temperature $T_d$ of $50$ K. Our fiducial IRX-stellar mass results at $z$$\sim$$7$$-$$8$ are consistent with marginal evolution from $z$$\sim$$0$. We then show how both results depend on $T_d$. For our six dust-detected sources, we estimate their dust masses and find that they are consistent with dust production from SNe if the dust destruction is low ($<$$90$%). Finally we determine the contribution of dust-obscured star formation to the star formation rate density for $UV$ luminous ($<$$-$$21.5$ mag: $\gtrsim$$1.7$$L_{UV} ^*$) $z$$\sim$$7$$-$$8$ galaxies, finding that the total SFR density at $z$$\sim$$7$ and $z$$\sim$$8$ from bright galaxies is $0.18_{-0.10}^{+0.08}$ dex and $0.20_{-0.09}^{+0.05}$ dex higher, respectively, i.e. $\sim$$\frac{1}{3}$ of the star formation in $\gtrsim$$1.7$$L_{UV} ^*$ galaxies at $z$$\sim$$7$$-$$8$ is obscured by dust.

preprint2020arXiv

ALMA Observations of Quasar Host Galaxies at $z\simeq4.8$

We present ALMA band-7 data of the [CII] $\lambda157.74\,μ{\rm m}$ emission line and underlying far-infrared (FIR) continuum for twelve luminous quasars at $z \simeq 4.8$, powered by fast-growing supermassive black holes (SMBHs). Our total sample consists of eighteen quasars, twelve of which are presented here for the first time. The new sources consists of six Herschel/SPIRE detected systems, which we define as "FIR-bright" sources, and six Herschel/SPIRE undetected systems, which we define as "FIR-faint" sources. We determine dust masses for the quasars hosts of $M_{dust} \le 0.2-25.0\times 10^8 M_{\odot}$, implying ISM gas masses comparable to the dynamical masses derived from the [CII] kinematics. It is found that on average the MgII line is blueshifted by $\sim 500\,{\rm km\,s}^{-1}$ with respect to the [CII] emission line, which is also observed when complementing our observations with data from the literature. We find that all of our "FIR-bright" subsample and most of the "FIR-faint" objects lie above the main sequence of star forming galaxies at $z \sim 5$. We detect companion sub-millimeter galaxies (SMGs) for two sources, both FIR-faint, with a range of projected distances of $\sim20-60$ kpc and with typical velocity shifts of $\left|Δv\right| \lesssim200\,{\rm km\,s}^{-1}$ from the quasar hosts. Of our total sample of eighteen quasars, 5/18 are found to have dust obscured starforming companions.

preprint2020arXiv

Cool Outflows in Galaxies and their Implications

Neutral-atomic and molecular outflows are a common occurrence in galaxies, near and far. They operate over the full extent of their galaxy hosts, from the innermost regions of galactic nuclei to the outermost reaches of galaxy halos. They carry a substantial amount of material that would otherwise have been used to form new stars. These cool outflows may have a profound impact on the evolution of their host galaxies and environments. This article provides an overview of the basic physics of cool outflows, a comprehensive assessment of the observational techniques and diagnostic tools used to characterize them, a detailed description of the best-studied cases, and a more general discussion of the statistical properties of these outflows in the local and distant universe. The remaining outstanding issues that have not yet been resolved are summarized at the end of the review to inspire new research directions.

preprint2020arXiv

Evidence for the Accretion of Gas in Star-Forming Galaxies: High N/O Abundances in Regions of Anomalously-Low Metallicity

While all models for the evolution of galaxies require the accretion of gas to sustain their growth via on-going star formation, it has proven difficult to directly detect this inflowing material. In this paper we use data of nearby star-forming galaxies in the SDSS IV Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey to search for evidence of accretion imprinted in the chemical composition of the interstellar medium. We measure both the O/H and N/O abundance ratios in regions previously identified as having anomalously low values of O/H. We show that the unusual locations of these regions in the N/O vs. O/H plane indicate that they have been created through the mixing of disk gas having higher metallicity with accreted gas having lower metallicity. Taken together with previous analysis on these anomalously low-metallicity regions, these results imply that accretion of metal-poor gas can probably sustain star formation in present-day late-type galaxies.

preprint2020arXiv

Exploring Chemical Homogeneity in Dwarf Galaxies: A VLT-MUSE study of JKB 18

Deciphering the distribution of metals throughout galaxies is fundamental in our understanding of galaxy evolution. Nearby, low-metallicity, star-forming dwarf galaxies in particular can offer detailed insight into the metal-dependent processes that may have occurred within galaxies in the early Universe. Here we present VLT/MUSE observations of one such system, JKB 18, a blue diffuse dwarf galaxy with a metallicity of only 12+log(O/H)=7.6$\pm$0.2 (~0.08 Zsol). Using high spatial-resolution integral field spectroscopy of the entire system, we calculate chemical abundances for individual HII regions using the direct method and derive oxygen abundance maps using strong-line metallicity diagnostics. With large-scale dispersions in O/H, N/H and N/O of ~0.5-0.6 dex and regions harbouring chemical abundances outside this 1$σ$ distribution, we deem JKB 18 to be chemically inhomogeneous. We explore this finding in the context of other chemically inhomogeneous dwarf galaxies and conclude that neither the accretion of metal-poor gas, short mixing timescales, or self-enrichment from Wolf-Rayet stars are accountable. Using a galaxy-scale, multi-phase, hydrodynamical simulation of a low-mass dwarf galaxy, we find that chemical inhomogeneities of this level may be attributable to the removal of gas via supernovae and the specific timing of the observations with respect to star-formation activity. This study not only draws attention to the fact that dwarf galaxies can be chemically inhomogeneous, but also that the methods used in the assessment of this characteristic can be subject to bias.

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

HIRES, the high-resolution spectrograph for the ELT

HIRES will be the high-resolution spectrograph of the European Extremely Large Telescope at optical and near-infrared wavelengths. It consists of three fibre-fed spectrographs providing a wavelength coverage of 0.4-1.8 mic (goal 0.35-1.8 mic) at a spectral resolution of ~100,000. The fibre-feeding allows HIRES to have several, interchangeable observing modes including a SCAO module and a small diffraction-limited IFU in the NIR. Therefore, it will be able to operate both in seeing and diffraction-limited modes. ELT-HIRES has a wide range of science cases spanning nearly all areas of research in astrophysics and even fundamental physics. Some of the top science cases will be the detection of bio signatures from exoplanet atmospheres, finding the fingerprints of the first generation of stars (PopIII), tests on the stability of Nature's fundamental couplings, and the direct detection of the cosmic acceleration. The HIRES consortium is composed of more than 30 institutes from 14 countries, forming a team of more than 200 scientists and engineers.

preprint2020arXiv

Metal enrichment in the circumgalactic medium and Lyα haloes around quasars at z $\sim$3

Deep observations have detected extended Ly$α$ emission nebulae surrounding tens of quasars at redshift 2 to 6. However, the metallicity of such extended haloes is still poorly understood. We perform a detailed analysis on a large sample of 80 quasars at $z\sim3$ based on MUSE-VLT data. We find clear evidence of extended emission of the UV nebular lines such as CIV or HeII for about 20$\%$ of the sample, while CIII] is only marginally detected in a few objects. By stacking the cubes we detect emission of CIV, HeII and CIII] out to a radius of about 45 kpc. CIV and HeII show a radial decline much steeper than Ly$α$, while CIII] shows a shallower profile similar to Ly$α$ in the inner 45 kpc. We infer that the average metallicity of the circumgalactic gas within the central 30-50~kpc is $\sim$0.5 solar, or even higher. However, we also find evidence of a component of the Ly$α$ haloes, which has much weaker metal emission lines relative to Ly$α$. We suggest that the high metallicity of the circumgalactic medium within the central 30-50 kpc is associated with chemical pre-enrichment by past quasar-driven outflows and that there is a more extended component of the CGM that has much lower metallicity and likely associated with near-pristine gas accreted from the intergalactic medium. We show that our observational results are in good agreement with the expectations of the FABLE zoom-in cosmological simulations.

preprint2020arXiv

The KLEVER Survey: Spatially resolved metallicity maps and gradients in a sample of 1.2 < z < 2.5 lensed galaxies

We present near-infrared observations of 42 gravitationally lensed galaxies obtained in the framework of the KMOS Lensed Emission Lines and VElocity Review (KLEVER) Survey, a program aimed at investigating the spatially resolved properties of the ionised gas in 1.2<z<2.5 galaxies by means of a full coverage of the YJ, H and K near-infrared bands. Detailed metallicity maps and gradients are derived for a sub-sample of 28 galaxies from reconstructed source plane emission line maps, exploiting the variety of different emission line diagnostics provided by the broad wavelength coverage of the survey. About 85% of these galaxies are characterised by metallicity gradients shallower than 0.05 dex/kpc and 89% are consistent with flat slope within 3$σ$ (67% within 1$σ$), suggesting a mild evolution with cosmic time. In the context of cosmological simulations and chemical evolution models, the presence of efficient feedback mechanisms and/or extended star formation profiles on top of the classical "inside-out" scenario of mass assembly is generally required to reproduce the observed flatness of the metallicity gradients beyond z$\sim$1 . Three galaxies with significantly (> 3$σ$) "inverted" gradients are also found, showing an anti-correlation between metallicity and star formation rate density on local scales, possibly suggesting recent episodes of pristine gas accretion or strong radial flows in place. Nevertheless, the individual metallicity maps are characterised by a variety of different morphologies, with flat radial gradients sometimes hiding non-axisymmetric variations on kpc scales which are washed out by azimuthal averages, especially in interacting systems or in those undergoing local episodes of recent star formation.

preprint2020arXiv

What Drives Galaxy Quenching? Resolving Molecular Gas and Star Formation in the Green Valley

We study quenching in seven green valley galaxies on kpc scales by resolving their molecular gas content using \textsuperscript{12}CO(1-0) observations obtained with NOEMA and ALMA, and their star-formation rate using spatially resolved optical spectroscopy from the MaNGA survey. We perform radial stacking of both datasets to increase the sensitivity to molecular gas and star formation, thereby avoiding biases against strongly quenched regions. We find that both spatially resolved gas fraction ($\rm {f_{gas}}$) and star formation efficiency ($\rm {SFE}$) are responsible for quenching green valley galaxies at all radii: both quantities are suppressed with respect to typical star-forming regions. $\rm {f_{gas}}$ and $\rm {SFE}$ have roughly equal influence in quenching the outer disc. We are, however, unable to identify the dominant mechanism in the strongly quenched central regions. We find that $\rm {f_{gas}}$ is reduced by $\rm \sim 1~dex$ in the central regions, but the star formation rate is too low to be measured, leading to upper limits for the $\rm {SFE}$. Moving from the outer disc to central regions, the reduction in $\rm {f_{gas}}$ is driven by an increasing $\rm Σ_{\star}$ profile rather than a decreasing $\rm Σ_{H_{2}}$ profile. The reduced $\rm {f_{gas}}$ may therefore be caused by a decrease in the gas supply rather than molecular gas ejection mechanisms, such as winds driven by active galactic nuclei. We warn more generally that studies investigating $\rm {f_{gas}}$ may be deceiving in inferring the cause of quenching, particularly in the central (bulge-dominated) regions of galaxies.

preprint2019arXiv

Are galactic star formation and quenching governed by local, global or environmental phenomena?

We present an analysis of star formation and quenching in the SDSS-IV MaNGA-DR15, utilising over 5 million spaxels from $\sim$3500 local galaxies. We estimate star formation rate surface densities ($Σ_{\rm SFR}$) via dust corrected $Hα$ flux where possible, and via an empirical relationship between specific star formation rate (sSFR) and the strength of the 4000 Angstrom break (D4000) in all other cases. We train a multi-layered artificial neural network (ANN) and a random forest (RF) to classify spaxels into `star forming' and `quenched' categories given various individual (and groups of) parameters. We find that global parameters (pertaining to the galaxy as a whole) perform collectively the best at predicting when spaxels will be quenched, and are substantially superior to local/ spatially resolved and environmental parameters. Central velocity dispersion is the best single parameter for predicting quenching in central galaxies. We interpret this observational fact as a probable consequence of the total integrated energy from AGN feedback being traced by the mass of the black hole, which is well known to correlate strongly with central velocity dispersion. Additionally, we train both an ANN and RF to estimate $Σ_{\rm SFR}$ values directly via regression in star forming regions. Local/ spatially resolved parameters are collectively the most predictive at estimating $Σ_{\rm SFR}$ in these analyses, with stellar mass surface density at the spaxel location ($Σ_*$) being by far the best single parameter. Thus, quenching is fundamentally a global process but star formation is governed locally by processes within each spaxel.

preprint2019arXiv

Both starvation and outflows drive galaxy quenching

Star-forming galaxies can in principle be transformed into passive systems by a multitude of processes that quench star formation, such as the halting of gas accretion (starvation) or the rapid removal of gas in AGN-driven outflows. However, it remains unclear which processes are the most significant, primary drivers of the SF-passive bimodality. We address this key issue in galaxy evolution by studying the chemical properties of 80,000 local galaxies in SDSS DR7. In order to distinguish between different quenching mechanisms, we analyse the stellar metallicities of star-forming, green valley and passive galaxies. We find that the significant difference in stellar metallicity between passive galaxies and their star-forming progenitors implies that for galaxies at all masses, quenching must have involved an extended phase of starvation. However, some form of gas ejection also has to be introduced into our models to best match the observed properties of local passive galaxies, indicating that, while starvation is likely to be the prerequisite for quenching, it is the combination of starvation and outflows that is responsible for quenching the majority of galaxies. Closed-box models indicate that the duration of the quenching phase is 2-3 Gyr, with an $e$-folding time of 2-4 Gyr, after which further star formation is prevented by an ejective/heating mode. Alternatively, leaky-box models find a longer duration for the quenching phase of 5-7 Gyr and an $e$-folding time of $\sim$1 Gyr, with outflows becoming increasingly important with decreasing stellar mass. Finally, our analysis of local green valley galaxies indicates that quenching is slower in the local Universe than at high-redshift.

preprint2019arXiv

Enhanced UV radiation and dense clumps in Mrk231's molecular outflow

We present interferometric observations of the CN(1-0) line emission in Mrk231 and combine them with previous observations of CO and other H$_2$ gas tracers to study the physical properties of the massive molecular outflow. We find a strong boost of the CN/CO(1-0) line luminosity ratio in the outflow, which is unprecedented compared to any other known Galactic or extragalactic source. For the dense gas phase in the outflow traced by the HCN and CN emissions, we infer $\rm X_{\rm CN}\equiv [CN]/[H_2] > X_{\rm HCN}$ by at least a factor of three, with H$_2$ gas densities of $n_{\rm H_2}\sim10^{5-6}$ cm$^{-3}$. In addition, for the first time, we resolve narrow spectral features in the HCN(1-0) and HCO$^+$(1-0) high-velocity line wings tracing the dense phase of the outflow. The velocity dispersions of these spectral features, $σ_v\sim7-20$ km s$^{-1}$, are consistent with those of massive extragalactic giant molecular clouds detected in nearby starburst nuclei. The H$_2$ gas masses inferred from the HCN data are quite high, $M_{mol}\sim0.3-5\times10^8$ $M_{\odot}$. Our results suggest that massive, denser molecular gas complexes survive embedded into the more diffuse H$_2$ phase of the outflow, and that the chemistry of such outflowing dense clouds is affected by enhanced UV radiation.

preprint2016arXiv

High-Velocity Bipolar Molecular Emission from an AGN Torus

We have detected in ALMA observations CO J = 6 - 5 emission from the nucleus of the Seyfert galaxy NGC 1068. The low-velocity (up to +/- 70 km/s relative to systemic) CO emission resolves into a 12x7 pc structure, roughly aligned with the nuclear radio source. Higher-velocity emission (up to +/- 400 km/s) is consistent with a bipolar outflow in a direction nearly perpendicular (roughly 80 degrees) to the nuclear disk. The position-velocity diagram shows that in addition to the outflow, the velocity field may also contain rotation about the disk axis. These observations provide compelling evidence in support of the disk-wind scenario for the AGN obscuring torus.

preprint2016arXiv

Nitrogen and oxygen abundances in the Local Universe

We present chemical evolution models aimed at reproducing the observed (N/O) vs. (O/H) abundance pattern of star forming galaxies in the Local Universe. We derive gas-phase abundances from SDSS spectroscopy and a complementary sample of low-metallicity dwarf galaxies, making use of a consistent set of abundance calibrations. This collection of data clearly confirms the existence of a plateau in the (N/O) ratio at very low metallicity, followed by an increase of this ratio up to high values as the metallicity increases. This trend can be interpreted as due to two main sources of nitrogen in galaxies: i) massive stars, which produce small amounts of pure primary nitrogen and are responsible for the (N/O) ratio in the low metallicity plateau; ii) low- and intermediate-mass stars, which produce both secondary and primary nitrogen and enrich the interstellar medium with a time delay relative to massive stars, and cause the increase of the (N/O) ratio. We find that the length of the low-metallicity plateau is almost solely determined by the star formation efficiency, which regulates the rate of oxygen production by massive stars. We show that, to reproduce the high observed (N/O) ratios at high (O/H), as well as the right slope of the (N/O) vs. (O/H) curve, a differential galactic wind - where oxygen is assumed to be lost more easily than nitrogen - is necessary. No existing set of stellar yields can reproduce the observed trend without assuming differential galactic winds. Finally, considering the current best set of stellar yields, a bottom-heavy initial mass function is favoured to reproduce the data.

preprint2016arXiv

Outflows and complex stellar kinematics in SDSS star forming galaxies

We investigate the properties of star formation-driven outflows by using a large spectroscopic sample of ~160,000 local "normal" star forming galaxies, drawn from the SDSS, spanning a wide range of star formation rates and stellar masses. The galaxy sample is divided into a fine grid of bins in the M_*-SFR parameter space, for each of which we produce a composite spectrum by stacking together the SDSS spectra of the galaxies contained in that bin. We exploit the high signal-to-noise of the stacked spectra to study the emergence of faint features of optical emission lines that may trace galactic outflows and would otherwise be too faint to detect in individual galaxy spectra. We adopt a novel approach that relies on the comparison between the line-of-sight velocity distribution (LoSVD) of the ionised gas (as traced by the [OIII]5007 and Halpha+[NII]6548,6583 emission lines) and the LoSVD of the stars, which are used as a reference tracing virial motions. Significant deviations of the gas kinematics from the stellar kinematics in the high velocity tail of the LoSVDs are interpreted as a signature of outflows. Our results suggest that the incidence of ionised outflows increases with SFR and sSFR. The outflow velocity (v_out) correlates tightly with the SFR for SFR>1 M_Sun/yr, whereas at lower SFRs the dependence of v_out on SFR is nearly flat. The outflow velocity, although with a much larger scatter, increases also with the stellar velocity dispersion, and we infer velocities as high as v_out~(6-8)*sigma_stars. Strikingly, we detect the signature of ionised outflows only in galaxies located above the main sequence (MS) of star forming galaxies in the M_*-SFR diagram, and the incidence of such outflows increases sharply with the offset from the MS. Our complementary analysis of the stellar kinematics reveals the presence of blue asymmetries of a few 10 km/s in the stellar LoSVDs. [abridged]

preprint2016arXiv

SDSS IV MaNGA - Spatially resolved diagnostic diagrams: A proof that many galaxies are LIERs

We study the spatially resolved excitation properties of the ionised gas in a sample of 646 galaxies using integral field spectroscopy data from SDSS-IV MaNGA. Making use of Baldwin-Philips-Terlevich diagnostic diagrams we demonstrate the ubiquitous presence of extended (kpc scale) low ionisation emission-line regions (LIERs) in both star forming and quiescent galaxies. In star forming galaxies LIER emission can be associated with diffuse ionised gas, most evident as extra-planar emission in edge-on systems. In addition, we identify two main classes of galaxies displaying LIER emission: `central LIER' (cLIER) galaxies, where central LIER emission is spatially extended, but accompanied by star formation at larger galactocentric distances, and `extended LIER' (eLIER) galaxies, where LIER emission is extended throughout the whole galaxy. In eLIER and cLIER galaxies, LIER emission is associated with radially flat, low H$α$ equivalent width of line emission ($<$ 3 Å) and stellar population indices demonstrating the lack of young stellar populations, implying that line emission follows tightly the continuum due to the underlying old stellar population. The H$α$ surface brightness radial profiles are always shallower than $\rm 1/r^{2}$ and the line ratio [OIII]$λ$5007/[OII]$λ$3727,29 (a tracer of the ionisation parameter of the gas) shows a flat gradient. This combined evidence strongly supports the scenario in which LIER emission is not due to a central point source but to diffuse stellar sources, the most likely candidates being hot, evolved (post-asymptotic giant branch) stars. Shocks are observed to play a significant role in the ionisation of the gas only in rare merging and interacting systems.

preprint2016arXiv

SDSS-IV MaNGA IFS Galaxy Survey --- Survey Design, Execution, and Initial Data Quality

The MaNGA Survey (Mapping Nearby Galaxies at Apache Point Observatory) is one of three core programs in the Sloan Digital Sky Survey IV. It is obtaining integral field spectroscopy (IFS) for 10K nearby galaxies at a spectral resolution of R~2000 from 3,622-10,354A. The design of the survey is driven by a set of science requirements on the precision of estimates of the following properties: star formation rate surface density, gas metallicity, stellar population age, metallicity, and abundance ratio, and their gradients; stellar and gas kinematics; and enclosed gravitational mass as a function of radius. We describe how these science requirements set the depth of the observations and dictate sample selection. The majority of targeted galaxies are selected to ensure uniform spatial coverage in units of effective radius (Re) while maximizing spatial resolution. About 2/3 of the sample is covered out to 1.5Re (Primary sample), and 1/3 of the sample is covered to 2.5Re (Secondary sample). We describe the survey execution with details that would be useful in the design of similar future surveys. We also present statistics on the achieved data quality, specifically, the point spread function, sampling uniformity, spectral resolution, sky subtraction, and flux calibration. For our Primary sample, the median r-band signal-to-noise ratio is ~73 per 1.4A pixel for spectra stacked between 1-1.5 Re. Measurements of various galaxy properties from the first year data show that we are meeting or exceeding the defined requirements for the majority of our science goals.

preprint2016arXiv

SDSS-IV MaNGA: A serendipitous observation of a potential gas accretion event

The nature of warm, ionized gas outside of galaxies may illuminate several key galaxy evolutionary processes. A serendipitous observation by the MaNGA survey has revealed a large, asymmetric H$α$ complex with no optical counterpart that extends $\approx8"$ ($\approx6.3$ kpc) beyond the effective radius of a dusty, starbursting galaxy. This H$α$ extension is approximately three times the effective radius of the host galaxy and displays a tail-like morphology. We analyze its gas-phase metallicities, gaseous kinematics, and emission-line ratios, and discuss whether this H$α$ extension could be diffuse ionized gas, a gas accretion event, or something else. We find that this warm, ionized gas structure is most consistent with gas accretion through recycled wind material, which could be an important process that regulates the low-mass end of the galaxy stellar mass function.

preprint2016arXiv

SDSS-IV MaNGA: stellar population gradients as a function of galaxy environment

We study the internal radial gradients of stellar population properties within $1.5\;R_{\rm e}$ and analyse the impact of galaxy environment. We use a representative sample of 721 galaxies with masses ranging between $10^{9}\;M_{\odot}$ to $10^{11.5}\;M_{\odot}$ from the SDSS-IV survey MaNGA. We split this sample by morphology into early-type and late-type galaxies. Using the full spectral fitting code FIREFLY, we derive the light and mass-weighted stellar population properties age and metallicity, and calculate the gradients of these properties. We use three independent methods to quantify galaxy environment, namely the $N^{th}$ nearest neighbour, the tidal strength parameter $Q$ and distinguish between central and satellite galaxies. In our analysis, we find that early-type galaxies generally exhibit shallow light-weighted age gradients in agreement with the literature and mass-weighted median age gradients tend to be slightly positive. Late-type galaxies, instead, have negative light-weighted age gradients. We detect negative metallicity gradients in both early and late-type galaxies that correlate with galaxy mass, with the gradients being steeper and the correlation with mass being stronger in late-types. We find, however, that stellar population gradients, for both morphological classifications, have no significant correlation with galaxy environment for all three characterisations of environment. Our results suggest that galaxy mass is the main driver of stellar population gradients in both early and late-type galaxies, and any environmental dependence, if present at all, must be very subtle.

preprint2015arXiv

Extending Virial Black Hole Mass Estimates to Low-Luminosity or Obscured AGN: the cases of NGC 4395 and MCG -01-24-012

In the last decade, using single epoch (SE) virial based spectroscopic optical observations, it has been possible to measure the black hole (BH) mass on large type 1 Active Galactic Nuclei (AGN) samples. However this kind of measurements can not be applied on those obscured type 2 and/or low luminosity AGN where the nuclear component does not dominate in the optical. We have derived new SE relationships, based on the FWHM and luminosity of the broad line region component of the Pabeta emission line and/or the hard X-ray luminosity in the 14-195 keV band, which have the prospect of better working with low luminosity or obscured AGN. The SE relationships have been calibrated in the 10^5-10^9 M_sol mass range, using a sample of AGN whose BH masses have been previously measured using reverberation mapping techniques. Our tightest relationship between the reverberation-based BH mass and the SE virial product has an intrinsic spread of 0.20 dex. Thanks to these SE relations, in agreement with previous estimates, we have measured a BH mass of M_BH =1.7^+1.3_-0.7 X 10^5 M_sol for the low luminosity, type 1, AGN NGC 4395 (one of the smallest active galactic BH known). We also measured, for the first time, a BH mass of M_BH = 1.5^+1.1_-0.6 X 10^7 M_sol for the Seyfert 2 galaxy MCG -01-24-012.

preprint2015arXiv

From haloes to Galaxies - I. The dynamics of the gas regulator model and the implied cosmic sSFR history

We explore the basic parameters that drive the evolution of the fundamental properties of star forming galaxies within the gas regulator model, or bathtub-model. We derive the general analytic form of the evolution of the key galaxy properties, i.e. gas mass, star formation rate (SFR), stellar mass, specific SFR, gas fraction, gas phase metallicity and stellar metallicity, without assuming that galaxies live in the equilibrium state. We find that the timescale required to reach equilibrium, tau_eq, which is determined by the product of star-formation efficiency and mass-loading factor, is the central parameter that is essentially in control of the evolution of all key galaxy properties. The scatters in most of the key scaling relations, such as the mass-SFR relation and mass-metallicity relation, are primarily governed by tau_eq. Most strikingly, the predicted sSFR evolution is controlled solely by tau_eq (apart from the cosmic time). Although the precise evolution of the sSFR depends on tau_eq, the sSFR history is largely insensitive to different values of tau_eq. The difference between the minimum and maximum sSFR at any epoch is less than a factor of four. The shape of the predicted sSFR history simply mimics that of the specific mass increase rate of the dark matter halos (sMIR_DM) with the typical value of the sSFR around 2*sMIR_DM, in good agreement with the predictions from typical Semi-Analytic Models (SAMs), but both are fundamentally different from the observed sSFR history. This clearly implies that some key process is missing in both typical SAMs and gas regulator model, and we hint at some possible culprit. We emphasize the critical role of tau_eq in controlling the evolution of the galaxy population, especially for gas rich low mass galaxies that are very unlikely to live around the equilibrium state at any epoch and this has been largely ignored in many similar studies.

preprint2015arXiv

Modern yields per stellar generation: the effect of the IMF

Gaseous and stellar metallicities in galaxies are nowadays routinely used to constrain the evolutionary processes in galaxies. This requires the knowledge of the average yield per stellar generation, $y_{\text{Z}}$, i.e. the quantity of metals that a stellar population releases into the interstellar medium (ISM), which is generally assumed to be a fixed fiducial value. Deviations of the observed metallicity from the expected value of $y_{\text{Z}}$ are used to quantify the effect of outflows or inflows of gas, or even as evidence for biased metallicity calibrations or inaccurate metallicity diagnostics. Here we show that $\rm y_{\text{Z}}$ depends significantly on the Initial Mass Function (IMF), varying by up to a factor larger than three, for the range of IMFs typically adopted in various studies. Varying the upper mass cutoff of the IMF implies a further variation of $y_{\text{Z}}$ by an additional factor that can be larger than two. These effects, along with the variation of the gas mass fraction restored into the ISM by supernovae ($R$, which also depends on the IMF), may yield to deceiving results, if not properly taken into account. In particular, metallicities that are often considered unusually high can actually be explained in terms of yield associated with commonly adopted IMFs such as the Kroupa (2001) or Chabrier (2003). We provide our results for two different sets of stellar yields (both affected by specific limitations) finding that the uncertainty introduced by this assumption can be as large as $\sim0.2$ dex. Finally, we show that $y_{\text{Z}}$ is not substantially affected by the initial stellar metallicity as long as $\text{Z}> 10^{-3}~\text{Z}_{\odot}$.

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

P-MaNGA: Full spectral fitting and stellar population maps from prototype observations

MaNGA (Mapping Nearby Galaxies at Apache Point Observatory) is a 6-year SDSS-IV survey that will obtain resolved spectroscopy from 3600 $Å$ to 10300 $Å$ for a representative sample of over 10,000 nearby galaxies. In this paper, we derive spatially resolved stellar population properties and radial gradients by performing full spectral fitting of observed galaxy spectra from P-MaNGA, a prototype of the MaNGA instrument. These data include spectra for eighteen galaxies, covering a large range of morphological type. We derive age, metallicity, dust and stellar mass maps, and their radial gradients, using high spectral-resolution stellar population models, and assess the impact of varying the stellar library input to the models. We introduce a method to determine dust extinction which is able to give smooth stellar mass maps even in cases of high and spatially non-uniform dust attenuation. With the spectral fitting we produce detailed maps of stellar population properties which allow us to identify galactic features among this diverse sample such as spiral structure, smooth radial profiles with little azimuthal structure in spheroidal galaxies, and spatially distinct galaxy sub-components. In agreement with the literature, we find the gradients for galaxies identified as early-type to be on average flat in age, and negative (- 0.15 dex / R$_e$ ) in metallicity, whereas the gradients for late-type galaxies are on average negative in age (- 0.39 dex / R$_e$ ) and flat in metallicity. We demonstrate how different levels of data quality change the precision with which radial gradients can be measured. We show how this analysis, extended to the large numbers of MaNGA galaxies, will have the potential to shed light on galaxy structure and evolution.

preprint2015arXiv

Strangulation as the primary mechanism for shutting down star formation in galaxies

Local galaxies are broadly divided into two main classes, star-forming (gas-rich) and quiescent (passive and gas-poor). The primary mechanism responsible for quenching star formation in galaxies and transforming them into quiescent and passive systems is still unclear. Sudden removal of gas through outflows or stripping is one of the mechanisms often proposed. An alternative mechanism is so-called "strangulation", in which the supply of cold gas to the galaxy is halted. Here we report that the difference between quiescent and star forming galaxies in terms of stellar metallicity (i.e. the fraction of metals heavier than helium in stellar atmospheres) can be used to discriminate efficiently between the two mechanisms. The analysis of the stellar metallicity in local galaxies, from 26,000 spectra, clearly reveals that strangulation is the primary mechanism responsible for quenching star formation, with a typical timescale of 4 billion years, at least for local galaxies with a stellar mass less than 10^11 solar masses. This result is further supported independently by the stellar age difference between quiescent and star-forming galaxies, which indicates that quiescent galaxies of less than 10^11 solar masses are on average observed four billion years after quenching due to strangulation.

preprint2015arXiv

XMM-Newton Observations of Three Interacting Luminous Infrared Galaxies

We investigate the X-ray properties of three interacting luminous infrared galaxy systems. In one of these systems, IRAS 18329+5950, we resolve two separate sources. A second, IRAS 20550+1656, and third, IRAS 19354+4559, have only a single X-ray source detected. We compare the observed emission to PSF profiles and determine that three are extended in emission. One is compact, which is suggestive of an AGN, although all of our profiles have large uncertainties. We then model the spectra to determine soft (0.5--2 keV) and hard (2--10 keV) luminosities for the resolved sources and then compare these to relationships found in the literature between infrared and X-ray luminosities for starburst galaxies. We obtain luminosities of $\log(L_{\textrm{soft}}/\textrm{L}_{\odot}) = 7.32,\:7.06,\:7.68$ and $\log(L_{\textrm{hard}}/\textrm{L}_{\odot}) = 7.33,\: 7.07,\: 7.88$ for IRAS 18329+5950, IRAS 19354+4559, and IRAS 20550+1656, respectively. These are intermediate to two separate predictions in the literature for star-formation-dominated sources. Our highest quality spectrum of IRAS 20550+1656 suggests super-solar abundance of alpha elements at $2σ$ significance, with $\log(\fracα{α_{\odot}}) = [α] = 0.4\pm0.2$. This is suggestive of recent enrichment with Type II supernovae, consistent with a starburst environment. The X-ray properties of the target galaxies are most likely due to starbursts, but we cannot conclusively rule out AGN.

preprint2014arXiv

Black hole accretion preferentially occurs in gas rich galaxies

We have investigated the gas content of a sample of several hundred AGN host galaxies at z$<$1 and compared it with a sample of inactive galaxies, matched in bins of stellar mass and redshift. Gas masses have been inferred from the dust masses, obtained by stacked Herschel far-IR and sub-mm data in the GOODS and COSMOS fields, under reasonable assumptions and metallicity scaling relations for the dust-to-gas ratio. We find that AGNs are on average hosted in galaxies much more gas rich than inactive galaxies. In the vast majority of stellar mass bins, the average gas content of AGN hosts is higher than in inactive galaxies. The difference is up to a factor of ten higher in low stellar mass galaxies, with a significance of 6.5$σ$. In almost half of the AGN sample the gas content is three times higher than in the control sample of inactive galaxies. Our result strongly suggests that the probability of having an AGN activated is simply driven by the amount of gas in the host galaxy; this can be explained in simple terms of statistical probability of having a gas cloud falling into the gravitational potential of the black hole. The increased probability of an AGN being hosted by a star-forming galaxy, identified by previous works, may be a consequence of the relationship between gas content and AGN activity, found in this paper, combined with the Schmidt-Kennicutt law for star formation.

preprint2014arXiv

First CO(17-16) emission line detected in a z > 6 quasar

We report the serendipitous detection of the CO(17-16) emission line toward the quasar SDSSJ114816.64+525150.3 (J1148) at redshift z = 6.4 obtained with the Plateau de Bure Interferometer. The CO(17-16) line is possibly contaminated by OH+ emission, that may account for ~ 35 - 60% of the total flux observed. Photo-Dissociation and X-ray Dominated Regions (PDRs and XDRs) models show that PDRs alone cannot reproduce the high luminosity of the CO(17-16) line relative to low-J CO transitions and that XDRs are required. By adopting a composite PDR+XDR model we derive molecular cloud and radiation field properties in the nuclear region of J1148. Our results show that highly excited CO lines represent a sensitive and possibly unique tool to infer the presence of X-ray faint or obscured supermassive black hole progenitors in high-z galaxies.

preprint2014arXiv

Ionized gas outflows and global kinematics of low-z luminous star forming galaxies

We study the kinematic properties of the ambient ionized ISM and ionized gas outflows in a large and representative sample of local luminous and ultraluminous infrared galaxies (U/LIRGs) (58 systems, 75 galaxies), on the basis of integral field spectroscopy (IFS)-based high S/N integrated spectra at galactic and sub-galactic, i.e. star forming (SF) clumps, scales. Ambient ionized gas. The velocity dispersion of the ionized ISM in U/LIRGs (<sigma> ~ 70 kms-1) is larger than in lower luminosity local star forming galaxies (<sigma> ~ 25 kms-1). While for isolated disc LIRGs star formation appears to sustain turbulence, gravitational energy release associated to interactions and mergers plays an important role driving sigma in the U/LIRG range. We also find that the impact of an AGN in ULIRGs is strong, increasing sigma by a factor 1.5 on average. The observed weak dependency of sigma with SFR surface density for local U/LIRGs is in very good agreement with that measured in some high-z samples. Ionized outflows. The presence of ionized gas outflows in U/LIRGs seems universal based on the detection of a broad, usually blueshifted, Halpha line. AGNs in U/LIRGs are able to generate faster (x2) and more massive (x1.4) ionized gas outflows than pure starbursts. The derived ionized mass loading factors are in general below one, with only a few AGNs above this limit. Only a small fraction of the ionized material from low mass LIRGs (log(Mdyn/Msun) < 10.4) could reach the intergalactic medium, with more massive galaxies retaining the gas. The observed average outflow properties in U/LIRGs are similar to high-z galaxies of comparable SFR. In the bright SF clumps found in LIRGs, ionized gas outflows appear to be very common. For a given SFR surface density, outflows in LIRG clumps would be about one to two orders of magnitude less energetic than those launched by clumps in high-z SF galaxies.

preprint2014arXiv

Overview of the SDSS-IV MaNGA Survey: Mapping Nearby Galaxies at Apache Point Observatory

We present an overview of a new integral field spectroscopic survey called MaNGA (Mapping Nearby Galaxies at Apache Point Observatory), one of three core programs in the fourth-generation Sloan Digital Sky Survey (SDSS-IV) that began on 2014 July 1. MaNGA will investigate the internal kinematic structure and composition of gas and stars in an unprecedented sample of 10,000 nearby galaxies. We summarize essential characteristics of the instrument and survey design in the context of MaNGA's key science goals and present prototype observations to demonstrate MaNGA's scientific potential. MaNGA employs dithered observations with 17 fiber-bundle integral field units that vary in diameter from 12" (19 fibers) to 32" (127 fibers). Two dual-channel spectrographs provide simultaneous wavelength coverage over 3600-10300 A at R~2000. With a typical integration time of 3 hr, MaNGA reaches a target r-band signal-to-noise ratio of 4-8 (per A, per 2" fiber) at 23 AB mag per sq. arcsec, which is typical for the outskirts of MaNGA galaxies. Targets are selected with stellar mass greater than 1e9 Msun using SDSS-I redshifts and i-band luminosity to achieve uniform radial coverage in terms of the effective radius, an approximately flat distribution in stellar mass, and a sample spanning a wide range of environments. Analysis of our prototype observations demonstrates MaNGA's ability to probe gas ionization, shed light on recent star formation and quenching, enable dynamical modeling, decompose constituent components, and map the composition of stellar populations. MaNGA's spatially resolved spectra will enable an unprecedented study of the astrophysics of nearby galaxies in the coming 6 yr.

preprint2013arXiv

Accurate Sky Continuum Subtraction with Fibre-fed Spectrographs

Fibre-fed spectrographs now have throughputs equivalent to slit spectrographs. However, the sky subtraction accuracy that can be reached has often been pinpointed as one of the major issues associated with the use of fibres. Using technical time observations with FLAMES-GIRAFFE, two observing techniques, namely dual staring and cross beam-switching, were tested and the resulting sky subtraction accuracy reached in both cases was quantified. Results indicate that an accuracy of 0.6% on sky subtraction can be reached, provided that the cross beam-switching mode is used. This is very encouraging with regard to the detection of very faint sources with future fibre-fed spectrographs, such as VLT/MOONS or E-ELT/MOSAIC.

preprint2013arXiv

Dust formation in the winds of AGBs: the contribution at low metallicities

We present new models for the evolution of stars with mass in the range 1Msun < M < 7.5Msun, followed from the pre-main-sequence through the asymptotic giant branch phase. The metallicity adopted is $Z=3*10^{-4} (which, with an alpha-enhancement of +0.4, corresponds to [Fe/H]=-2). Dust formation is described by following the growth of dust grains of various types as the wind expands from the stellar surface. Models with mass M>3Msun experience Hot Bottom Burning, thus maintaining the surface C/O below unity. Unlike higher Z models, the scarcity of silicon available in the envelope prevents the formation of silicates in meaningful quantities, sufficient to trigger the acceleration of the wind via radiation pressure on the dust grains formed. No silicate formation occurs below a threshold metallicity of Z=10^{-3}. Low--mass stars, with M< 2.5Msun become carbon stars, forming solid carbon dust in their surroundings. The total dust mass formed depends on the uncertain extent of the inwards penetration of the convective envelope during the Third Dredge--Up episodes following the Thermal Pulses. Carbon grains have sizes 0.08 micron < a_C < 0.12 micron and the total amount of dust formed (increasing with the mass of the star) is M_C=(2-6)*10^{-4}Msun. Our results imply that AGB stars with Z=3*10^{-4} can only contribute to carbon dust enrichment of the interstellar medium on relatively long timescales, > 300 Myr, comparable to the evolutionary time of a 3Msun star. At lower metallicities the scarcity of silicon available and the presence of Hot Bottom Burning even in M< 2Msun, prevents the formation of silicate and carbon grains. We extrapolate our conclusion to more metal--poor environments, and deduce that at Z < 10^{-4} dust enrichment is mostly due to metal condensation in supernova ejecta.

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

SPICA spectroscopic cosmological surveys to unravel galaxy evolution

The main energy-generating mechanisms in galaxies are black hole (BH) accretion and star formation (SF) and the interplay of these processes is driving the evolution of galaxies. MIR/FIR spectroscopy are able to distinguish between BH accretion and SF, as it was shown in the past by infrared spectroscopy from the space by the Infrared Space Observatory and Spitzer. Spitzer and Herschel spectroscopy together can trace the AGN and the SF components in galaxies, with extinction free lines, almost only in the local Universe, except for a few distant objects. One of the major goals of the study of galaxy evolution is to understand the history of the luminosity source of galaxies along cosmic time. This goal can be achieved with far-IR spectroscopic cosmological surveys. SPICA in combination with ground based large single dish submillimeter telescopes, such as CCAT, will offer a unique opportunity to do this. We use galaxy evolution models linked to the observed MIR-FIR counts (including Herschel) to predict the number of sources and their IR lines fluxes, as derived from observations of local galaxies. A shallow survey in an area of 0.5 square degrees, with a typical integration time of 1 hour per pointing, will be able to detect thousands of galaxies in at least three emission lines, using SAFARI, the far-IR spectrometer onboard of SPICA.

preprint2013arXiv

The dependence of the galaxy mass-metallicity relation on environment and the implied metallicity of the IGM

We explore the dependence of the galaxy mass-metallicity relation on environment in SDSS, in terms of both over-density and central/satellite dichotomy. We find that at a given stellar mass, there is a strong dependence of metallicity on over-density for star-forming satellites. Instead, for star-forming centrals no correlation is found. Star-forming satellites at different stellar masses form a tight sequence in the average overdensity - metallicity plane, which covers the entire observed range of metallicities and stellar masses. This remarkable result appears to imply that there exists a universal evolutionary path for all star-forming satellites, regardless of their stellar masses. The strong correlation between over-density and metallicity for star-forming satellites indicates that the gas inflow of satellite galaxies is progressively metal-enriched in denser regions. We interpret our results by employing the gas regulator model and find that the metallicity of the enriched inflow of star-forming satellite galaxies strongly increases with increasing over-density, largely independently of stellar mass. If the metallicity of the inflow of star-forming satellites can represent the metallicity of the IGM, then the implied metallicity of the IGM rises from ~ 0.01 solar metallicity in void-like environments to ~ 0.3 solar metallicity in cluster-like environments, in broad agreement with observations. We show that the observed metallicity difference between star-forming centrals and star-forming satellites becoming smaller towards high stellar masses can be simply explained by the mass-independent enriched inflow, without the need to involve any mass-dependent environmental effect on metallicity. Since satellite galaxies account for at least half of the galaxy population, our findings prompt for a revision of many galaxy evolutionary models, which generally assume pristine gas inflows.

preprint2013arXiv

The Hot and Energetic Universe: A White Paper presenting the science theme motivating the Athena+ mission

This White Paper, submitted to the recent ESA call for science themes to define its future large missions, advocates the need for a transformational leap in our understanding of two key questions in astrophysics: 1) How does ordinary matter assemble into the large scale structures that we see today? 2) How do black holes grow and shape the Universe? Hot gas in clusters, groups and the intergalactic medium dominates the baryonic content of the local Universe. To understand the astrophysical processes responsible for the formation and assembly of these large structures, it is necessary to measure their physical properties and evolution. This requires spatially resolved X-ray spectroscopy with a factor 10 increase in both telescope throughput and spatial resolving power compared to currently planned facilities. Feedback from supermassive black holes is an essential ingredient in this process and in most galaxy evolution models, but it is not well understood. X-ray observations can uniquely reveal the mechanisms launching winds close to black holes and determine the coupling of the energy and matter flows on larger scales. Due to the effects of feedback, a complete understanding of galaxy evolution requires knowledge of the obscured growth of supermassive black holes through cosmic time, out to the redshifts where the first galaxies form. X-ray emission is the most reliable way to reveal accreting black holes, but deep survey speed must improve by a factor ~100 over current facilities to perform a full census into the early Universe. The Advanced Telescope for High Energy Astrophysics (Athena+) mission provides the necessary performance (e.g. angular resolution, spectral resolution, survey grasp) to address these questions and revolutionize our understanding of the Hot and Energetic Universe. These capabilities will also provide a powerful observatory to be used in all areas of astrophysics.

preprint2013arXiv

The LBT Bootes Field Survey: I. The Rest-frame UV and Luminosity Functions and Clustering of Bright Lyman Break Galaxies at z~3

We present a deep LBT/LBC U-band imaging survey (9 deg2) covering the NOAO Bootes field. A total of 14,485 Lyman Break Galaxies (LBGs) at z~3 are selected, which are used to measure the rest-frame UV luminosity function (LF). The large sample size and survey area reduce the LF uncertainties due to Poisson statistics and cosmic variance by >3 compared to previous studies. At the bright end, the LF shows excess power compared to the best-fit Schechter function, which can be attributed to the contribution of $z\sim3$ quasars. We compute the rest-frame near-infrared LF and stellar mass function (SMF) of z~3 LBGs based on the R-band and IRAC [4.5 micro m]-band flux relation. We investigate the evolution of the UV LFs and SMFs between z~7 and z~3, which supports a rising star formation history in the LBGs. We study the spatial correlation function of two bright LBG samples and estimate their average host halo mass. We find a tight relation between the host halo mass and the galaxy star formation rate (SFR),which follows the trend predicted by the baryonic accretion rate onto the halo, suggesting that the star formation in LBGs is fueled by baryonic accretion through the cosmic web. By comparing the SFRs with the total baryonic accretion rates, we find that cosmic star formation efficiency is about 5%-20% and it does not evolve significantly with redshift, halo mass, or galaxy luminosity.

preprint2012arXiv

A UV Ultra-luminous Lyman Break Galaxy at Z=2.78 in NDWFS Bootes Field

We present one of the most ultraviolet (UV) luminous Lyman Break Galaxies (LBGs) (J1432+3358) at z=2.78, discovered in the NOAO Deep Wide-Field Survey (NDWFS) Bootes field. The R-band magnitude of J1432+3358 is 22.29 AB, more than two magnitudes brighter than typical L* LBGs at this redshift. The deep z-band image reveals two components of J1432+3358 separated by 1.0" with flux ratio of 3:1. The high signal-to-noise ratio (S/N) rest-frame UV spectrum shows Lya emission line and interstellar medium absorption lines. The absence of NV and CIV emission lines, the non-detection in X-ray and radio wavelengths and mid-infrared (MIR) colors indicate no or weak active galactic nuclei (AGN) (<10%) in this galaxy. The galaxy shows broader line profile with the full width half maximum (FWHM) of about 1000 km/s and larger outflow velocity (~500 km/s) than those of typical z~3 LBGs. The physical properties are derived by fitting the spectral energy distribution (SED) with stellar synthesis models. The dust extinction, E(B-V)=0.12, is similar to that in normal LBGs. The star formation rates (SFRs) derived from the SED fitting and the dust-corrected UV flux are consistent with each other, ~300 Msun/yr, and the stellar mass is 1.3e11 Msun. The SFR and stellar mass in J1432+3358 are about an order of magnitude higher than those in normal LBGs. The SED-fitting results support that J1432+3358 has a continuous star formation history with the star formation episode of 630 Myr. The morphology of J1432+3358 and its physical properties suggest that J1432+3358 is in an early phase of 3:1 merger process. The unique properties and the low space number density (~1e-7 Mpc^{-3})are consistent with the interpretation that such galaxies are either found in a short unobscured phase of the star formation or that small fraction of intensive star-forming galaxies are unobscured.

preprint2012arXiv

AGN Obscuration and the Unified Model

Unification Models of Active Galactic Nuclei postulate that all the observed differences between Type 1 and Type 2 objects are due to orientation effects with respect to the line-of-sight to the observer. The key ingredient of these models is the obscuring medium, historically envisaged as a toroidal structure on a parsec scale. However, many results obtained in the last few years are clearly showing the need for a more complex geometrical distribution of the absorbing media. In this paper we review the various pieces of evidence for obscuring media on different scales, from the vicinity of the black hole to the host galaxy, in order to picture an updated unification scenario explaining the complex observed phenomenology. We conclude by mentioning some of the open issues.

preprint2012arXiv

ALMA reveals a chemically evolved submillimeter galaxy at z=4.76

The chemical properties of high-z galaxies provide important information to constrain galaxy evolutionary scenarios. However, widely-used metallicity diagnostics based on rest-frame optical emission lines are not usable for heavily dust-enshrouded galaxies (such as Sub-Millimeter Galaxies; SMGs), especially at z>3. Here we focus on the flux ratio of the far-infrared fine-structure emission lines [NII]205um and [CII]158um to assess the metallicity of high-z SMGs. Through ALMA cycle 0 observations, we have detected the [NII]205um emission in a strongly [CII]-emitting SMG, LESS J033229.4-275619 at z=4.76. The velocity-integrated [NII]/[CII] flux ratio is 0.043 +/- 0.008. This is the first measurement of the [NII]/[CII] flux ratio in high-z galaxies, and the inferred flux ratio is similar to the ratio observed in the nearby universe (~0.02-0.07). The velocity-integrated flux ratio and photoionization models suggest that the metallicity in this SMG is consistent with solar, implying the chemical evolution has progressed very rapidly in this system at z=4.76. We also obtain a tight upper limit on the CO(12-11) transition, which translates into CO(12-11)/CO(2-1) <3.8 (3 sigma). This suggests that the molecular gas clouds in LESS J033229.4-275619 are not affected significantly by the radiation field emitted by the AGN in this system.

preprint2012arXiv

An X-ray WHIM metal absorber from a Mpc-scale empty region of space

We report a detection of an absorption line at ~44.8 Å in a > 500 ks Chandra HRC-S/LETG X-ray grating spectrum of the blazar H 2356-309. This line can be identified as intervening CV-Kα absorption, at z\approx0.112, produced by a warm (log T = 5.1 K) intergalactic absorber. The feature is significant at a 2.9σ level (accounting for the number of independent redshift trials). We estimate an equivalent hydrogen column density of log N_H=19.05 (Z/Zsun)^-1 cm^-2. Unlike other previously reported FUV/X-ray metal detections of warm-hot intergalactic medium (WHIM), this CV absorber lies in a region with locally low galaxy density, at ~2.2 Mpc from the closest galaxy at that redshift, and therefore is unlikely to be associated with an extended galactic halo. We instead tentatively identify this absorber with an intervening Warm-Hot Intergalactic Medium filament possibly permeating a large-scale, 30 Mpc extended, structure of galaxies whose redshift centroid, within a cylinder of 7.5 Mpc radius centered on the line of sight to H 2356-309, is marginally consistent (at a 1.8σ level) with the redshift of the absorber.

preprint2012arXiv

Quasar feedback in the early Universe: the case of SDSS J1148+5251

Galaxy-scale gas outflows triggered by active galactic nuclei have been proposed as a key physical process to regulate the co-evolution of nuclear black holes and their host galaxies. The recent detection of a massive gas outflow in one of the most distant quasar, SDSS J1148+5251 at z = 6.4, presented by Maiolino et al. (2012) strongly supports this idea and suggests that strong quasar feedback is already at work at very early times. In a previous work, Valiante et al. (2011), we have presented a hierarchical semi-analytical model, GAMETE/ QSOdust, for the formation and evolution of high-redshift quasars, and we have applied it to the quasar SDSS J1148+5251, with the aim of investigating the star formation history, the nature of the dominant stellar populations and the origin and properties of the large dust mass observed in the host galaxy. A robust prediction of the model is that the evolution of the nuclear black hole and of the host galaxy are tightly coupled by quasar feedback in the form of strong galaxy-scale winds. In the present letter, we show that the gas outflow rate predicted by GAMETE/QSOdust is in good agreement with the lower limit of 3500 Msun/yr inferred by the observations. According to the model, the observed outflow at z = 6.4 is dominated by quasar feedback, as the outflow rate has already considerably depleted the gas content of the host galaxy, leading to a down-turn in the star formation rate at z < 7 - 8. Hence, supernova explosions give a negligible contribution to the observed winds at z = 6.4.

preprint2012arXiv

Scaling relations of metallicity, stellar mass, and star formation rate in metal-poor starbursts: I. A fundamental plane

Most galaxies follow well-defined scaling relations of metallicity (O/H), star formation rate (SFR), and stellar mass. However, low-metallicity starbursts, rare in the Local Universe but more common at high redshift, deviate significantly from these scaling relations. On the "main sequence" of star formation, these galaxies have high SFR for a given M*; and on the mass-metallicity relation, they have excess M* for their low metallicity. In this paper, we characterize O/H, M*, and SFR for these deviant "low-metallicity starbursts", selected from a sample of ~1100 galaxies, spanning almost two orders of magnitude in metal abundance, a factor of ~10^6 in SFR, and of ~10^5 in stellar mass. Our sample includes quiescent star-forming galaxies and blue compact dwarfs at redshift 0, luminous compact galaxies at redshift 0.3, and Lyman Break galaxies at redshifts 1-3.4. Applying a Principal Component Analysis (PCA) to the galaxies in our sample with M*<10^{10} Msun gives a Fundamental Plane (FP) of scaling relations; SFR and stellar mass define the plane itself, and O/H its thickness. The dispersion for our sample in the edge-on view of the plane is 0.17 dex, independently of redshift and including the metal-poor starbursts. The same FP is followed by 55100 galaxies selected from the Sloan Digital Sky Survey, with a dispersion of 0.06dex. In a companion paper, we develop multi-phase chemical evolution models that successfully predict the observed scaling relations and the FP; the deviations from the main scaling relations are caused by a different (starburst or "active") mode of star formation. These scaling relations do not truly evolve, but rather are defined by the different galaxy populations dominant at different cosmological epochs.

preprint2012arXiv

Scaling relations of metallicity, stellar mass, and star formation rate in metal-poor starbursts: II. Theoretical models

Scaling relations of metallicity (O/H), star formation rate (SFR), and stellar mass give important insight on galaxy evolution. They are obeyed by most galaxies in the Local Universe and also at high redshift. In a companion paper, we compiled a sample of ~1100 galaxies from redshift 0 to ~3, spanning almost two orders of magnitude in metal abundance, a factor of $\sim10^6$ in SFR, and of ~10^5 in stellar mass. We have characterized empirically the star-formation "main sequence" (SFMS) and the mass-metallicity relation (MZR) for this sample, and also identified a class of low-metallicity starbursts, rare locally but more common in the distant universe. These galaxies deviate significantly from the main scaling relations, with high SFR and low metal content for a given M*. In this paper, we model the scaling relations and explain these deviations from them with a set of multi-phase chemical evolution models based on the idea that, independently of redshift, initial physical conditions in a galaxy's evolutionary history can dictate its location in the scaling relations. Our models are able to successfully reproduce the O/H, M*, and SFR scaling relations up to z~3, and also successfully predict the molecular cloud fraction as a function of stellar mass. These results suggest that the scaling relations are defined by different modes of star formation: an "active" starburst mode, more common at high redshift, and a quiescent "passive" mode that is predominant locally and governs the main trends.

preprint2011arXiv

Chemical properties in the most distant radio galaxy

We present a deep optical spectrum of TN J0924-2201, the most distant radio galaxy at z = 5.19, obtained with FOCAS on the Subaru Telescope. We successfully detect, for the first time, the CIV1549 emission line from the narrow-line region (NLR). In addition to the emission-line fluxes of Ly alpha and CIV, we set upper limits on the NV and HeII emissions. We use these line detections and upper limits to constrain the chemical properties of TN J0924-2201. By comparing the observed emission-line flux ratios with photoionization models, we infer that the carbon-to-oxygen relative abundance is already [C/O] > -0.5 at a cosmic age of ~ 1.1 Gyr. This lower limit on [C/O] is higher than the ratio expected at the earliest phases of the galaxy chemical evolution, indicating that TN J0924-2201 has already experienced significant chemical evolution at z = 5.19.

preprint2011arXiv

Enhanced [CII] emission in a z=4.76 submillimetre galaxy

We present the detection of bright [CII] emission in the z=4.76 submillimetre galaxy LESS J033229.4-275619 using the Atacama Pathfinder EXperiment. This represents the highest redshift [CII] detection in a submm selected, star-formation dominated system. The AGN contributions to the [CII] and far-infrared (FIR) luminosities are small. We find an atomic mass derived from [CII] comparable to the molecular mass derived from CO. The ratio of the [CII], CO and FIR luminosities imply a radiation field strength G_0~10^3 and a density ~10^4 cm^-3 in a kpc-scale starburst, as seen in local and high redshift starbursts. The high L_[CII]/L_FIR=2.4x10^-3 and the very high L_[CII]/L_CO(1-0) ~ 10^4 are reminiscent of low metallicity dwarf galaxies, suggesting that the highest redshift star-forming galaxies may also be characterised by lower metallicities. We discuss the implications of a reduced metallicity on studies of the gas reservoirs, and conclude that especially at very high redshift, [CII] may be a more powerful and reliable tracer of the interstellar matter than CO.

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

Stellar metallicity of star-forming galaxies at z ~ 3

The stellar metallicity is a direct measure of the amount of metals present in a galaxy, as a large part of the metals lie in its stars. In this paper we investigate new stellar metallicity indicators suitable for high-z galaxies studying the stellar photospheric absorption lines in the rest frame ultraviolet, hence sampling predominantly young hot stars. We defined these new indicators based on the equivalent widths (EW) of selected features using theoretical spectra created with the evolutionary population synthesis code Starburts99. We used them to compute the stellar metallicity for a sample of UV-selected galaxies at z > 3 from the AMAZE survey using very deep (37h per object) VLT/FORS spectra. Moreover, we applied the new metallicity indicators to eight additional high redshift galaxies found in literature. We then compared stellar and gas-phase metallicities measured from the emission lines for all these galaxies, finding that within the errors the two estimates are in good agreement, with possible tendency to have stellar metallicities lower than the gas phase ones. For the first time, we are able to study the stellar mass-metallicity relation at z > 3. We find that the metallicity of young, hot stars in galaxies at z \sim 3 have similar values of the aged stars in local SDSS galaxies, contrary to what observed for the gas phase metallicity.

preprint2011arXiv

Sub-millimeter Detected z~2 Radio-quiet QSOs: Accurate Redshifts, Black Hole Masses, and Inflow/Outflow Velocities

We present near-IR spectroscopy of a sample of luminous (Lbol >10^14 Lsun), sub-millimeter-detected, dusty (M_dust~10^9 Msun), radio-quiet QSOs at z~2. We measure the redshift, width, and luminosity of the Halpha line, and the continuum luminosity near Ha. The Halpha emission line was strongly detected in all sources. Two components -- a very broad (>~5000 km/s) Gaussian and an intermediate-width (>~1500 km/s) Gaussian, were required to fit the Halpha profile of all observed QSOs. Narrow (<~1000 km/s) lines were not detected in the sample QSOs. Relative to the Halpha-derived QSO redshift, the rest-frame UV emission lines in these sub-mm bright QSOs show larger than average blue-shifted velocities, potentially tracing strong -- up to 3000 km/s -- outflows in the Broad Line Region. With the exception of the one QSO which shows exceptionally broad Halpha lines, the black hole masses of the QSO sample are in the range log Mbh = 9.0--9.7 and the Eddington fractions are between 0.5 and ~1. In black hole mass and accretion rate, this sub-mm bright QSO sample is indistinguishable from the Shemmer et al. (2004) optically-bright QSO sample at z~2; the latter is likely dominated by sub-mm dim QSOs. Previous authors have demonstrated a correlation, over six orders of magnitude, between SFR and accretion rate in active galaxies: the sub-mm bright QSOs lie at the upper extremes of both quantities and their SFR is an order of magnitude higher than that predicted from the correlation

preprint2011arXiv

The metallicity properties of zCOSMOS galaxies at 0.2<z<0.8

We study the metallicity properties of galaxies in the zCOSMOS sample between 0.2<z<0.8. At z<0.46, where Ha and [NII] are detected, we find the same dependence of metallicity on stellar mass and Star Formation Rate (SFR), the Fundamental Metallicity Relation, found by Mannucci et al. (2010) in SDSS galaxies on a similar redshift range. We extend this relation to higher redshift, 0.49<z<0.8 where the R23 metallicity index can be measured in our data, finding no evidence for evolution, and a metallicity scatter around the relation of about 0.16 dex. This result confirms, with a much higher level of significance with respect to previous works, the absence of evolution of the FMR during the last half of cosmic history.

preprint2010arXiv

Confirmation of X-Ray Absorption by WHIM in the Sculptor Wall

In a previous paper we reported a 3-sigma detection of an absorption line from the Warm-Hot Intergalactic Medium (WHIM) using the Chandra and XMM X-ray grating spectra of the blazar H2356-309, the sight-line of which intercepts the Sculptor Wall, a large-scale superstructure of galaxies at z ~ 0.03. To verify our initial detection, we obtained a deep (500 ks), follow-up exposure of H2356-309 as part of the Cycle-10 Chandra Large Project Program. From a joint analysis of the Cycle-10 and previous (Cycle-8) Chandra grating data we detect the redshifted OVII WHIM line at a significance level of 3.4-sigma, a substantial improvement over the 1.7-sigma level reported previously when using only the Cycle-8 data. The significance increases to 4.0-sigma when the existing XMM grating data are included in the analysis, thus confirming at higher significance the existence of the line at the redshift of the Sculptor Wall with an equivalent width of 28.5+/-10.5 mA (90% confidence). We obtain a 90% lower limit on the OVII column density of 0.8 10^16 cm^-2 and a 90% upper limit on the Doppler-b parameter of 460 km/s. Assuming the absorber is uniformly distributed throughout the ~ 15 Mpc portion of the blazar's sight-line that intercepts the Sculptor Wall, that the OVII column density is ~ 2 10^16 cm^-2 (corresponding to b > 150 km/s where the inferred column density is only weakly dependent on b), and that the oxygen abundance is 0.1 solar, we estimate a baryon over-density of ~ 30 for the WHIM, which is consistent with the peak of the WHIM mass fraction predicted by cosmological simulations. The clear detection of OVII absorption in the Sculptor Wall demonstrates the viability of using current observatories to study WHIM in the X-ray absorption spectra of blazars behind known large-scale structures.

preprint2010arXiv

Metallicity Diagnostics with Infrared Fine-Structure Lines

Although measuring the gas metallicity in galaxies at various redshifts is crucial to constrain galaxy evolutionary scenarios, only rest-frame optical emission lines have been generally used to measure the metallicity. This has prevented us to accurately measure the metallicity of dust-obscured galaxies, and accordingly to understand the chemical evolution of dusty populations, such as ultraluminous infrared galaxies. Here we propose diagnostics of the gas metallicity based on infrared fine structure emission lines, which are nearly unaffected by dust extinction even the most obscured systems. Specifically, we focus on fine-structure lines arising mostly from HII regions, not in photo-dissociation regions, to minimize the dependence and uncertainties of the metallicity diagnostics from various physical parameters. Based on photoionization models, we show that the emission-line flux ratio of ([OIII]51.80+[OIII]88.33)/[NIII]57.21 is an excellent tracer of the gas metallicity. The individual line ratios [OIII]51.80/[NIII]57.21 or [OIII]88.33/[NIII]57.21 can also be used as diagnostics of the metallicity, but they suffer a stronger dependence on the gas density. The line ratios [OIII]88.33/[OIII]51.80 and [NII]121.7/[NIII]57.21 can be used to measure and, therefore, account for the dependences on the of the gas density and ionization parameter, respectively. We show that these diagnostic fine-structure lines are detectable with Herschel in luminous infrared galaxies out z=0.4. Metallicity measurements with these fine-structure lines will be feasible at relatively high redshift (z=1 or more) with SPICA, the future infrared space observatory.

preprint2010arXiv

Quasar feedback revealed by giant molecular outflows

In the standard scenario for galaxy evolution young star-forming galaxies transform into red bulge-dominated spheroids, where star formation has been quenched. To explain such a transformation, a strong negative feedback generated by accretion onto a super-massive black hole is often invoked. The depletion of gas resulting from quasar-driven outflows should eventually stop star-formation in the host galaxy and lead the black hole to "suicide" for starvation. Direct observational evidence for a major quasar feedback onto the host galaxy is still missing, since outflows previously observed in quasars are generally associated with the ionized component of the gas, which only accounts for a minor fraction of the total gas content, and typically occurring in the central regions. We used the IRAM PdB Interferometer to observe CO(1-0) in Mrk 231, the closest QSO known. Thanks to the wide band we detect broad wings of the CO line, with velocities up to 750 km/s and spatially resolved on the kpc scale. Such broad CO wings trace a giant molecular outflow of about 700 MSun/yr, far larger than the ongoing star-formation rate (~200 MSun/yr) observed in the host galaxy. This wind will totally expel the cold gas reservoir in Mrk 231 in about 1e7 yrs, therefore halting the star-formation activity on the same timescale. The inferred kinetic energy in the molecular outflow is ~1.2e44 erg/s, corresponding to a few percent of the AGN bolometric luminosity, very close to the fraction expected by models ascribing quasar feedback to highly supersonic shocks generated by radiatively accelerated nuclear winds. The direct observational evidence for quasar feedback reported here provides solid support to the scenarios ascribing the observed properties of local massive galaxies to quasar-induced large scale winds.

preprint2010arXiv

The Building Up of the Black Hole Mass - Stellar Mass Relation

We derive the growth of SMBHs relative to the stellar content of their host galaxy predicted under the assumption of BH accretion triggered by galaxy encounters occurring during their merging histories. We show that, within this framework, the ratio Gamma=(M_BH/M_*)(z)/(M_BH/M_*)(z=0) between the Black Hole mass and the galactic stellar mass (normalized to the local value) depends on both BH mass and redshift. While the average value and the spread of Gamma(z) increase with z, such an effect is larger for massive BHs, reaching values Gamma=5 for massive Black Holes (M>10^9 M_{\odot}) at z>4, in agreement with recent observations of high-redshift QSOs; this is due to the effectiveness of interactions in triggering BH accretion in high-density environments at high redshifts. To test such a model against observations, we worked out specific predictions for sub-samples of the simulated galaxies corresponding to the different observational samples for which measurements of Gamma have been obtained. We found that for Broad Line AGNs at 1<z<2 values of Gamma=2 are expected, with a mild trend toward larger value for increasing BH mass. Instead, when we select from our Monte Carlo simulations only extremely gas rich, rapidly star forming galaxies at 2<z<3, we find low values 0.3<Gamma<1.5, consistent with recent observational findings on samples of sub-mm galaxies; in the framework of our model, these objects end up at z=0 in low-to-intermediate mass BHs (M<10^9 M_{\odot}), and they do not represent typical paths leading to local massive galaxies. The latter have formed preferentially through paths passing above the local M_*-M_BH relation. We discuss how the global picture emerging from the model is consistent with a downsizing scenario, where massive BHs accrete a larger fraction of their final mass at high redshifts z>4.

preprint2010arXiv

The mass-metallicity relation of SDSS quasars

Active galactic nuclei (AGNs) are characterized by a clear correlation between luminosity and metallicity (L_AGN-Z_AGN relation). The origin of this correlation is not clear. It may result from a relation between the black hole mass (M_BH) and metallicity, or from a relation between the accretion rate (L/L_Edd) and metallicity. To investigate the origin of the L_AGN-Z_AGN relation, we use optical spectra of 2383 quasars at 2.3 < z < 3.0 from the Sloan Digital Sky Survey. By using this data set, we have constructed composite spectra of 33 subsamples in intervals of both M_BH and L/L_Edd. From these composite spectra we measure emission-line flux ratios that are sensitive to the metallicity of the broad line region (BLR); specifically, NV/CIV, NV/HeII, (SiIV+OIV])/CIV, and AlIII/CIV. We find that there is a significant correlation between M_BH and Z_BLR as inferred from all four metallicity-sensitive emission-line flux ratios. This result strongly suggests that the observed L_AGN-Z_AGN relation is mostly a consequence of the M_BH-Z_AGN relation. The relation between M_BH and Z_BLR is likely a consequence of both the M_BH-M_bul relation and of the mass-metallicity relation in the host galaxy. We also find that L/L_Edd correlates with the emission line flux ratios involving NV (more specifically, NV/CIV and NV/HeII), while it does not correlate with the other two metallicity sensitive emission line flux ratios, i.e., (SiIV+OIV])/CIV and AlIII/CIV. These correlations indicate that the emission-line flux ratios involving NV depend on both metallicity and relative abundance of nitrogen. We suggest that the relation between L/L_Edd and those line ratios involving nitrogen, is caused by a delay of the black hole accretion rate relative to the onset of nuclear star formation of about 10^8 years, which is the timescale required for the nitrogen enrichment.

preprint2010arXiv

The Origin and Evolution of the Mass-Metallicity Relation using GalICS

The GalICS (Galaxies in Cosmological Simulations) semi-analytical model of hierar- chical galaxy formation is used to investigate the effects of different galactic properties, including star formation rate (SFR) and outflows, on the shape of the mass metallic- ity relation and to predict the relation for galaxies at redshift z=2.27 and z=3.54. Our version of GalICS has the chemical evolution implemented in great detail and is less heavily reliant on approximations such as instantaneous recycling. We vary the model parameters controlling both the efficiency and redshift dependence of the SFR as well as the efficiency of supernova feedback. We find that the factors controlling the SFR influence the relation significantly at all redshifts and require a strong redshift dependence, proportional to 1+z, in order to reproduce the observed relation at the low mass end. Indeed, at any redshift, the predicted relation flattens out at the high mass end resulting in a poorer agreement with observations in this regime. We also find that variation of the parameters associated with outflows has a minimal effect on the relation at high redshift but does serve to alter its shape in the more recent past. We thus conclude that the relation is one between SFR and mass and that outflows are only important in shaping the relation at late times. When the relation is stratified by SFR it is apparent that the predicted galaxies with increasing stellar masses have higher SFRs, supporting the view that galaxy downsizing is the origin of the relation. Attempting to reproduce the observed relation, we vary the parameters controlling the efficiency of star formation and its redshift dependence and compare the predicted relations with Erb et al. (2006) at z=2.27 and Maiolino et al. (2008) at z=3.54 in order to find the best-fitting parameters. (Abridged)

preprint2009arXiv

Chemical Evolution of High-Redshift Radio Galaxies

We present new deep optical spectra of 9 high-z radio galaxies (HzRGs) at z > 2.7 obtained with FORS2 on VLT. These rest-frame ultraviolet spectra are used to infer the metallicity of the narrow-line regions (NLRs) in order to investigate the chemical evolution of galaxies in high-z universe. We focus mainly on the CIV/HeII and CIII]/CIV flux ratios that are sensitive to gas metallicity and ionization parameter. Although the NV emission has been widely used to infer the gas metallicity, it is often too weak to be measured accurately for NLRs. By combining our new spectra with data from the literature, we examine the possible redshift evolution of the NLR metallicity for 57 HzRGs at 1 < z < 4. Based on the comparison between the observed emission-line flux ratios and the results of our photoionization model calculations, we find no significant metallicity evolution in NLRs of HzRGs, up to z ~ 4. Our results imply that massive galaxies had almost completed their chemical evolution at much higher redshift (z > 5). Finally, although we detect strong NV emission lines in 5 HzRGs at z > 2.7, we point out that high NV/HeII ratios are not indicative of high metallicities but correspond to high ionization parameters of gas clouds in NLRs.

preprint2009arXiv

Spitzer Infrared Low-Resolution Spectroscopic Study of Buried AGNs in a Complete Sample of Nearby Ultraluminous Infrared Galaxies

We present the results of Spitzer IRS low-resolution infrared 5-35 micron spectroscopy of 17 nearby ULIRGs at z < 0.2, optically classified as non-Seyferts. The presence of optically elusive, but intrinsically luminous, buried AGNs is investigated, based on the strengths of polycyclic aromatic hydrocarbon emission and silicate dust absorption features detected in the spectra. The signatures of luminous buried AGNs, whose intrinsic luminosities range up to ~10^12 Lsun, are found in eight sources. We combine these results with those of our previous research to investigate the energy function of buried AGNs in a complete sample of optically non-Seyfert ULIRGs in the local universe at z < 0.3 (85 sources). We confirm a trend that we previously discovered: that buried AGNs are more common in galaxies with higher infrared luminosities. Because optical Seyferts also show a similar trend, we argue more generally that the energetic importance of AGNs is intrinsically higher in more luminous galaxies, suggesting that the AGN-starburst connections are luminosity-dependent. This may be related to the stronger AGN feedback scenario in currently more massive galaxy systems, as a possible origin of the galaxy downsizing phenomenon.

preprint2006arXiv

Gas metallicity diagnostics in star-forming galaxies

Generally the gas metallicity in distant galaxies can only be inferred by using a few prominent emission lines. Various theoretical models have been used to predict the relationship between emission line fluxes and metallicity, suggesting that some line ratios can be used as diagnostics of the gas metallicity in galaxies. However, accurate empirical calibrations of these emission line flux ratios from real galaxy spectra spanning a wide metallicity range are still lacking. In this paper we provide such empirical calibrations by using the combination of two sets of spectroscopic data: one consisting of low-metallicity galaxies with a measurement of [OIII]4363, and the other one consisting of galaxies in the SDSS DR4 whose metallicity has been determined from various strong emission lines in their spectra. This combined data set constitutes the largest sample of galaxies with information on the gas metallicity available so far and spanning the widest metallicity range. Our empirical diagrams show that the line ratio [OIII]/[OII] is a useful tool to break the degeneracy in the R_23 parameter when no information on the [NII] line is available. The line ratio [NeIII]/[OII] also results to be a useful metallicity indicator for high-z galaxies. Finally, we compare these empirical relations with photoionization models. We find that the empirical R_23-metallicity sequence is strongly discrepant with respect to the trend expected by models with constant ionization parameter, which is interpreted as a consequence of a strong metallicity dependence of the average ionization parameter in galaxies. This result should warn about the use of theoretical models with constant ionization parameter to infer metallicities from observed line ratios. (abbreviated)