Source author record

R. Maiolino

R. 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

174works
9topics
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

174 published item(s)

preprint2026arXiv

GA-NIFS: high prevalence of dusty and metal-enriched outflows in massive and luminous star-forming galaxies at $z\sim3-9$

We present a search for and characterization of ionized outflows in 15 star-forming systems at $z\sim3-9$ with no evidence of Active Galactic Nuclei (AGN), observed with JWST/NIRSpec IFU as part of the GA-NIFS program. The targets often show satellites and complex substructure, from which we isolate 40 galaxies/regions. The sample probes the high-mass end of the galaxy population, with most sources having log$_{10}$~(M$_\star$/M$_\odot$)=$9.5-11$, extending previous studies on high-z star formation driven outflows that mainly focused on lower-mass galaxies. Using the [OIII]5007 and H$α$ emission lines, we identify broad kinematic components consistent with galactic outflows in 14 galaxies/regions. We find that the outflowing gas is more dust attenuated (by $A_{\rm V}$$\sim0.59$ mag on average) and metal-enriched (0.13 dex) than the interstellar medium (ISM) of the host galaxies, but its velocities are insufficient to escape the galaxies and reach the circumgalactic medium, suggesting that outflows mainly redistribute dust and metals around their hosts. The outflows identified in this study display velocity dispersions within $σ_{\rm out}=130-340$~km~s$^{-1}$ and outflow velocities $v_{\rm out}=170-600$~km~s$^{-1}$, and, when combined with less luminous and less massive star-forming galaxies from previous works, reveal a statistically significant correlation between $v_{\rm out}$ and star formation rate (SFR). The typically low mass-loading factors ($η=\dot{M}_{\rm out}/SFR$$\leq1$, in 9 out of 14 the outflows) indicate that these outflows do not strongly suppress star formation. Overall, our results suggest that ejective feedback via ionized outflows is inefficient in massive, luminous star-forming galaxies within the first 2 Gyr of the Universe.

preprint2024arXiv

HYPERION. Interacting companion and outflow in the most luminous $z>6$ quasar

We present ALMA deep observations of the [CII] 158 $μ$m emission line and the continuum at 253 GHz and 99 GHz towards SDSS J0100+2802 at $z\simeq 6.3$, the most luminous quasi-stellar object (QSO) at z$>$6. J0100+2802 belongs to the HYPERION sample of luminous QSOs at $z\sim 6-7.5$. The observations (at 2.2" resolution in Band 3 and 0.9" resolution in Band 6) are optimized to detect extended emission around the QSO. We detect an interacting, tidally disrupted companion both in [CII], peaking at $z\sim 6.332$, and in continuum, stretching on scales up to 20 kpc from the quasar, with a knotty morphology. The higher velocity dispersion in the direction of the companion emission and the complex morphology of tidally stretched galaxy suggest a possible ongoing or future merger. For the newly-detected companion we derive a range of dust mass, $M_{\rm dust}=(0.3-2.6)\times 10^7\ \rm M_\odot$, and of star formation rate, SFR$=[35-344]\ \rm M_\odot$. This shows that both the QSO and its companion are gas rich and that a major merger may be at the origin of the boosted star formation. We also detect a broad blueshifted component in the [CII] spectrum, that we interpret as a gaseous outflow for which we estimate a mass outflow rate in the range $\dot{M}_{\rm out}=(118-269)\ \rm M_\odot\ yr^{-1}$. J0100+2802 was recently found to reside in a strong overdensity, however this close companion remained undetected by both previous higher resolution ALMA observations and by JWST-NIRCAM imaging. Our results highlight the importance of deep medium-resolution ALMA observations for the study of QSOs and their environment at the Epoch of Reionization.

preprint2023arXiv

Radio antenna design for sky-averaged 21 cm cosmology experiments: the REACH case

Following the reported detection of an absorption profile associated with the 21~cm sky-averaged signal from the Cosmic Dawn by the EDGES experiment in 2018, a number of experiments have been set up to verify this result. This paper discusses the design process used for global 21~cm experiments, focusing specifically on the Radio Experiment for the Analysis of Cosmic Hydrogen (REACH). This experiment will seek to understand and compensate for systematic errors present using detailed modelling and characterization of the instrumentation. There is detailed the quantitative figures of merit and numerical modelling used to assist the design process of the REACH dipole antenna (one of the 2 antenna designs for REACH Phase I). This design process produced a 2.5:1 frequency bandwidth dipole. The aim of this design was to balance spectral smoothness and low impedance reflections with the ability to describe and understand the antenna response to the sky signal to inform the critically important calibration during observation and data analysis.

preprint2022arXiv

Bayesian hierarchical modelling of the $\mathrm{M_{\star}}$-SFR relation from 1<z<6 in ASTRODEEP

The Hubble Frontier Fields represent the opportunity to probe the high-redshift evolution of the main sequence of star-forming galaxies to lower masses than possible in blank fields thanks to foreground lensing of massive galaxy clusters. We use the BEAGLE SED-fitting code to derive stellar masses, $\mathrm{M_{\star}}=\log(M/\mathrm{M_{\odot}})$, SFRs, $Ψ=\log(ψ/\mathrm{M_{\odot}}\,\mathrm{yr}^{-1})$ and redshifts from galaxies within the ASTRODEEP catalogue. We fit a fully Bayesian hierarchical model of the main sequence over $1.25<z<6$ of the form $Ψ= α_\mathrm{9.7}(z) + β(\mathrm{M_{\star}}-9.7) + \mathcal{N}(0,σ^2)$ while explicitly modelling the outlier distribution. The redshift-dependent intercept at $\mathrm{M_{\star}}=9.7$ is parametrized as $α_\mathrm{9.7}(z) = \log[N (1+z)^γ] + 0.7$. Our results agree with an increase in normalization of the main sequence to high redshifts that follows the redshift-dependent rate of accretion of gas onto dark matter halos with $γ=2.40^{+0.18}_{-0.18}$. We measure a slope and intrinsic scatter of $β=0.79^{+0.03}_{-0.04}$ and $σ=0.26^{+0.02}_{-0.02}$. We find that the sampling of the SED provided by the combination of filters (Hubble + ground-based Ks-band + Spitzer 3.6 and 4.5 $\mathrm{μm}$) is insufficient to constrain $\mathrm{M_{\star}}$ and $Ψ$ over the full dynamic range of the observed main sequence, even at the lowest redshifts studied. While this filter set represents the best current sampling of high-redshift galaxy SEDs out to $z>3$, measurements of the main sequence to low masses and high redshifts still strongly depend on priors employed in SED fitting (as well as other fitting assumptions). Future data-sets with JWST should improve this.

preprint2022arXiv

Detection of companion galaxies around hot dust-obscured hyper-luminous galaxy W0410-0913

The phase transition between galaxies and quasars is often identified with the rare population of hyper-luminous, hot dust-obscured galaxies. Galaxy formation models predict these systems to grow via mergers, that can deliver large amounts of gas toward their centers, induce intense bursts of star formation and feed their supermassive black holes. Here we report the detection of 24 galaxies emitting Lyman-alpha emission on projected physical scales of about 400 kpc around the hyper-luminous hot dust-obscured galaxy W0410-0913, at redshift z = 3.631, using Very Large Telescope observations. While this indicates that W0410-0913 evolves in a very dense environment, we do not find clear signs of mergers that could sustain its growth. Data suggest that if mergers occurred, as models expect, these would involve less massive satellites, with only a moderate impact on the internal interstellar medium of W0410-0913, which is sustained by a rotationally-supported fast-rotating molecular disk, as Atacama Large Millimeter Array observations suggest.

preprint2022arXiv

Diagnostics for PopIII galaxies and Direct Collapse Black Holes in the early universe

Forthcoming observational facilities will make the exploration of the early universe routine, likely probing large populations of galaxies at very low metallicities. It will therefore be important to have diagnostics that can solidly identify and distinguish different classes of objects in such low metallicity regimes. We use new photoionisation models to develop diagnostic diagrams involving various nebular lines. We show that combinations of these diagrams allow the identification and discrimination of the following classes of objects in the early universe: PopIII and Direct Collapse Black Holes (DCBH) in pristine environments, PopIII and DCBH embedded in slightly enriched ISM (Z~10^{-5}-10^{-4}), (metal poor) PopII and AGN in enriched ISM. Diagnostics involving rest-frame optical lines (that will be accessible by JWST) have a better discriminatory power, but also rest-frame UV diagnostics can provide very useful information. Interestingly, we find that metal lines such as [OIII] 5007A and CIV 1549A can remain relatively strong (about a factor of 0.1-1 relative H-beta and HeII 1640A, respectively), even in extremely metal poor environments (Z~10^{-5}-10^{-4}), which could be embedding PopIII galaxies and DCBH.

preprint2022arXiv

Evidence of extended cold molecular gas and dust halos around $\mathbf{z\sim2.3}$ Extremely Red Quasars with ALMA

Large-scale outflows are believed to be an important mechanism in the evolution of galaxies. We can determine the impact of these outflows by studying either current galaxy outflows and their effect in the galaxy or by studying the effect of past outflows on the gas surrounding the galaxy. In this work, we examine the CO(7-6), [CI]\,($^{3} \rm P_{1} \rightarrow {\rm ^3 P}_{0}$), H$_2$O 2$_{11}$--2$_{02}$ and dust continuum emission of 15 extremely red quasars (ERQs) at z$\sim$2.3 using ALMA. By investigating the radial surface brightness profiles of both the individual sources and the stacked emission, we detect extended cold gas and dust emission on scales of $\sim$14 kpc in CO(7-6), [CI](2-1), and dust continuum. This is the first time that the presence of a large amount of molecular gas was detected on large, circum-galactic medium scales around quasar host galaxies using [CI] extended emission. We estimate the dust and molecular gas mass of these halos to be 10$^{7.6}$ and 10$^{10.6}$ M$_\odot$, indicating significant dust and molecular gas reservoirs around these extreme quasars. By estimating the timescale at which this gas can reach these distances by molecular gas outflows (7-32 Myr), we conclude that these halos are a relic of past AGN or starburst activity, rather than an effect of the current episode of extreme quasar activity.

preprint2022arXiv

The Near-Infrared Spectrograph (NIRSpec) on the James Webb Space Telescope I. Overview of the instrument and its capabilities

We provide an overview of the design and capabilities of the near-infrared spectrograph (NIRSpec) onboard the James Webb Space Telescope. NIRSpec is designed to be capable of carrying out low-resolution ($R\!=30\!-330$) prism spectroscopy over the wavelength range $0.6-5.3\!~μ$m and higher resolution ($R\!=500\!-1340$ or $R\!=1320\!-3600$) grating spectroscopy over $0.7-5.2\!~μ$m, both in single-object mode employing any one of five fixed slits, or a 3.1$\times$3.2 arcsec$^2$ integral field unit, or in multiobject mode employing a novel programmable micro-shutter device covering a 3.6$\times$3.4~arcmin$^2$ field of view. The all-reflective optical chain of NIRSpec and the performance of its different components are described, and some of the trade-offs made in designing the instrument are touched upon. The faint-end spectrophotometric sensitivity expected of NIRSpec, as well as its dependency on the energetic particle environment that its two detector arrays are likely to be subjected to in orbit are also discussed.

preprint2022arXiv

The Near-Infrared Spectrograph (NIRSpec) on the James Webb Space Telescope II. Multi-object spectroscopy (MOS)

We provide an overview of the capabilities and performance of the Near-Infrared Spectrograph (NIRSpec) on the James Webb Space Telescope (JWST) when used in its multi-object spectroscopy (MOS) mode employing a novel Micro Shutter Array (MSA) slit device. The MSA consists of four separate 98 arcsec $\times$ 91 arcsec quadrants each containing $365\times171$ individually addressable shutters whose open areas on the sky measure 0.20 arcsec $\times$ 0.46 arcsec on a 0.27 arcsec $\times$ 0.53 arcsec pitch. This is the first time that a configurable multi-object spectrograph has been available on a space mission. The levels of multiplexing achievable with NIRSpec MOS mode are quantified and we show that NIRSpec will be able to observe typically fifty to two hundred objects simultaneously with the pattern of close to a quarter of a million shutters provided by the MSA. This pattern is fixed and regular, and we identify the specific constraints that it yields for NIRSpec observation planning. We also present the data processing and calibration steps planned for the NIRSpec MOS data. The significant variation in size of the mostly diffraction-limited instrument point spread function over the large wavelength range of 0.6-5.3 $μ$m covered by the instrument, combined with the fact that most targets observed with the MSA cannot be expected to be perfectly centred within their respective slits, makes the spectrophotometric and wavelength calibration of the obtained spectra particularly complex. These challenges notwithstanding, the sensitivity and multiplexing capabilities anticipated of NIRSpec in MOS mode are unprecedented, and should enable significant progress to be made in addressing a wide range of outstanding astrophysical problems.

preprint2022arXiv

The Near-Infrared Spectrograph (NIRSpec) on the James Webb Space Telescope III. Integral-field spectroscopy

The Near-Infrared Spectrograph (NIRSpec) on the James Webb Space Telescope (JWST) offers the first opportunity to use integral-field spectroscopy from space at near-infrared wavelengths. More specifically, NIRSpec's integral-field unit can obtain spectra covering the wavelength range $0.6 - 5.3~μ$m for a contiguous 3.1 arcsec $\times$ 3.2 arcsec sky area at spectral resolutions of $R \approx 100$, 1000, and 2700. In this paper we describe the optical and mechanical design of the NIRSpec integral-field spectroscopy mode, together with its expected performance. We also discuss a few recommended observing strategies, some of which are driven by the fact that NIRSpec is a multipurpose instrument with a number of different observing modes, which are discussed in companion papers. We briefly discuss the data processing steps required to produce wavelength- and flux-calibrated data cubes that contain the spatial and spectral information. Lastly, we mention a few scientific topics that are bound to benefit from this highly innovative capability offered by JWST/NIRSpec.

preprint2022arXiv

The Near-Infrared Spectrograph (NIRSpec) on the James Webb Space Telescope IV. Capabilities and predicted performance for exoplanet characterization

The Near-Inrared Spectrograph (NIRSpec) on the James Webb Space Telescope (JWST) is a very versatile instrument, offering multiobject and integral field spectroscopy with varying spectral resolution ($\sim$30 to $\sim$3000) over a wide wavelength range from 0.6 to 5.3 micron, enabling scientists to study many science themes ranging from the first galaxies to bodies in our own Solar System. In addition to its integral field unit and support for multiobject spectroscopy, NIRSpec features several fixed slits and a wide aperture specifically designed to enable high precision time-series and transit as well as eclipse observations of exoplanets. In this paper we present its capabilities regarding time-series observations, in general, and transit and eclipse spectroscopy of exoplanets in particular. Due to JWST's large collecting area and NIRSpec's excellent throughput, spectral coverage, and detector performance, this mode will allow scientists to characterize the atmosphere of exoplanets with unprecedented sensitivity.

preprint2022arXiv

The PUMA project. III. Incidence and properties of ionised gas disks in ULIRGs, associated velocity dispersion and its dependence on starburstiness

A classical scenario suggests that ULIRGs transform colliding spiral galaxies into a spheroid dominated early-type galaxy. Recent high-resolution simulations have instead shown that, under some circumstances, rotation disks can be preserved during the merging process or rapidly regrown after coalescence. Our goal is to analyze in detail the ionised gas kinematics in a sample of ULIRGs to infer the incidence of gas rotational dynamics in late-stage interacting galaxies and merger remnants. We analysed MUSE data of a sample of 20 nearby (z<0.165) ULIRGs, as part of the "Physics of ULIRGs with MUSE and ALMA" (PUMA) project. We found that 27% individual nuclei are associated with kpc-scale disk-like gas motions. The rest of the sample displays a plethora of gas kinematics, dominated by winds and merger-induced flows, which make the detection of rotation signatures difficult. On the other hand, the incidence of stellar disk-like motions is ~2 times larger than gaseous disks, as the former are probably less affected by winds and streams. The eight galaxies with a gaseous disk present relatively high intrinsic gas velocity dispersion (sigma = 30-85 km/s), rotationally-supported motions (with gas rotation velocity over velocity dispersion vrot/sigma > 1-8), and dynamical masses in the range (2-7)x1e10 Msun. By combining our results with those of local and high-z disk galaxies from the literature, we found a significant correlation between sigma and the offset from the main sequence (MS), after correcting for their evolutionary trends. Our results confirm the presence of kpc-scale rotating disks in interacting galaxies and merger remnants, with an incidence going from 27% (gas) to ~50% (stars). The ULIRGs gas velocity dispersion is up to a factor of ~4 higher than in local normal MS galaxies, similar to high-z starbursts as presented in the literature.

preprint2022arXiv

Widespread, strong outflows in XQR-30 quasars at the Reionisation epoch

Luminous quasars powered by accretion onto billion solar mass black holes already exist at the epoch of Reionisation, when the Universe was 0.5-1 Gyr old. These objects likely reside in over-dense regions of the Universe, and will grow to form today's giant galaxies. How their huge black holes formed in such short times is debated, particularly as they lie above the local black hole mass-galaxy dynamical mass correlation, thus following the black hole-dominance growth path. It is unknown what slowed down the black hole growth, leading towards the symbiotic growth observed in the local Universe, and when this process started, although black hole feedback is a likely driver. This deadlock is due to the lack of large, homogeneous samples of high-redshift quasars with high-quality, broad-band spectroscopic information. Here we report results from a VLT/X-shooter survey of 30 quasars at redshift 5.8$\le$z$\le$6.6 (XQR-30). About 50% of their spectra reveal broad blue-shifted absorption line (BAL) throughs, tracing powerful ionised winds. The BAL fraction in z$\gtrsim$6 quasars is 2-3 times higher than in quasars at z~2-4.5. XQR-30 BAL quasars exhibit extreme outflow velocities, up to 17% of the light speed, rarely observed at lower redshift. These outflows inject large amounts of energy into the galaxy interstellar medium, which can contrast nuclear gas accretion, slowing down the black-hole growth. The star-formation rate in high-z quasar hosts is generally $>$100 M$_\odot$/yr, so these galaxies are growing at a fast rate. The BAL phase may then mark the beginning of significant feedback, acting first on black hole growth and possibly later on galaxy growth. The red optical colors of BAL quasars at z$\gtrsim$6 indeed suggest that these systems are dusty and may be caught during an initial quenching phase of obscured accretion.

preprint2021arXiv

MAGNUM survey: Compact jets causing large turmoil in galaxies -- Enhanced line widths perpendicular to radio jets as tracers of jet-ISM interaction

Outflows accelerated by AGN are commonly observed in the form of coherent, mildly collimated high-velocity gas directed along the AGN ionisation cones and kinetically powerful (>$10^{44-45}$ erg/s) jets. Recent works found that outflows can also be accelerated by low-power (<$10^{44}$ erg/s) jets, and the most recent cosmological simulations indicate that these are the dominant source of feedback on sub-kpc scales. We study the relation between radio jets and the distribution and kinematics of the ionised gas in IC 5063, NGC 5643, NGC 1068, and NGC 1386 as part of our MAGNUM survey of nearby Seyfert galaxies. All these objects host a small-scale (<1 kpc) low-power (<$10^{44}$ erg/s) radio jet that has small inclinations (<45°) with respect to the galaxy disc. We employed seeing-limited optical integral field spectroscopic observations from MUSE at VLT to obtain flux, kinematic, and excitation maps of the extended ionised gas, that we compared with archival radio images and Chandra X-ray observations. We detect a strong (up to >800-1000 km/s), extended (>1 kpc) and shock-excited emission-line velocity spread perpendicular to the AGN ionisation cones and jets in all four targets. These broad and symmetric line profiles are not associated with a single coherent velocity of the gas, differently from the 'classical' asymmetric-line outflow observed along the ionisation cones and jets. We interpret the observed phenomenon as due to the action of the jets perturbing the gas in the galaxy disc. These intense and extended velocity spreads perpendicular to AGN jets and cones are indeed currently only observed in galaxies hosting a low-power jet whose inclination is sufficiently low with respect to the galaxy disc to impact on and strongly affect its material. In line with cosmological simulations, our results demonstrate that low-power jets are indeed capable of affecting the host galaxy.

preprint2020arXiv

AION: An Atom Interferometer Observatory and Network

We outline the experimental concept and key scientific capabilities of AION (Atom Interferometer Observatory and Network), a proposed UK-based experimental programme using cold strontium atoms to search for ultra-light dark matter, to explore gravitational waves in the mid-frequency range between the peak sensitivities of the LISA and LIGO/Virgo/ KAGRA/INDIGO/Einstein Telescope/Cosmic Explorer experiments, and to probe other frontiers in fundamental physics. AION would complement other planned searches for dark matter, as well as probe mergers involving intermediate mass black holes and explore early universe cosmology. AION would share many technical features with the MAGIS experimental programme in the US, and synergies would flow from operating AION in a network with this experiment, as well as with other atom interferometer experiments such as MIGA, ZAIGA and ELGAR. Operating AION in a network with other gravitational wave detectors such as LIGO, Virgo and LISA would also offer many synergies.

preprint2020arXiv

Chemical abundances of Seyfert 2 AGNs-III. Reducing the oxygen abundance discrepancy

We investigate the discrepancy between oxygen abundance estimations for narrow-line regions (NLRs) of Active Galactic Nuclei (AGNs) type Seyfert 2 derived by using direct estimations of the electron temperature (Te-method) and those derived by using photoionization models. In view of this, observational emission-line ratios in the optical range (3000 < λ(Å) < 7000) of Seyfert 2 nuclei compiled from the literature were reproduced by detailed photoionization models built with the Cloudy code. We find that the derived discrepancies are mainly due to the inappropriate use of the relations between temperatures of the low (t2) and high (t3) ionization gas zones derived for H II regions in AGN chemical abundance studies. Using a photoionization model grid, we derived a new expression for t2 as a function of t3 valid for Seyfert 2 nuclei. The use of this new expression in the AGN estimation of the O/H abundances based on Te-method produces O/H abundances slightly lower (about 0.2 dex) than those derived from detailed photoionization models. We also find that the new formalism for the Te-method reduces by about 0.4 dex the O/H discrepancies between the abundances obtained from strong emission-line calibrations and those derived from direct estimations.

preprint2020arXiv

MOONRISE: The Main MOONS GTO Extragalactic Survey

The MOONS instrument possesses an exceptional combination of large multiplexing, high sensitivity, broad simultaneous spectral coverage (from optical to near-infrared bands), large patrol area and high fibre density. These properties provide the unprecedented potential of enabling, for the very first time, SDSS-like surveys around Cosmic Noon (z~1-2.5), when the star formation rate in the Universe peaked. The high-quality spectra delivered by MOONS will sample the same nebular and stellar diagnostics observed in extensive surveys of local galaxies, providing an accurate and consistent description of the evolution of various physical properties of galaxies, and hence a solid test of different scenarios of galaxy formation and transformation. Most importantly, by spectroscopically identifying hundreds of thousands of galaxies at high redshift, the MOONS surveys will be capable of determining the environments in which primeval galaxies lived and will reveal how such environments affected galaxy evolution. In this article, we specifically focus on the main Guaranteed Time Observation (GTO) MOONS extragalactic survey, MOONRISE, by providing an overview of its scientific goals and observing strategy.

preprint2020arXiv

MOONS Surveys of the Milky Way and its Satellites

The study of resolved stellar populations in the Milky Way and other Local Group galaxies can provide us with a fossil record of their chemo-dynamical and star-formation histories over timescales of many billions of years. In the galactic components and stellar systems of the Milky Way and its satellites, individual stars can be resolved. Therefore, they represent a unique laboratory in which to investigate the details of the processes behind the formation and evolution of the disc and dwarf/irregular galaxies. MOONS at the VLT represents a unique combination of an efficient infrared multi-object spectrograph and a large-aperture 8-m-class telescope which will sample the cool stellar populations of the dense central regions of the Milky Way and its satellites, delivering accurate radial velocities, metallicities, and other chemical abundances for several millions of stars over its lifetime (see Cirasuolo et al., this issue). MOONS will observe up to 1000 targets across a 25-arcminute field of view in the optical and near-infrared (0.6-1.8 micron) simultaneously. A high-resolution (R~19700) setting in the H band has been designed for the accurate determination of stellar abundances such as alpha, light, iron-peak and neutron-capture elements.

preprint2020arXiv

Physics of ULIRGs with MUSE and ALMA: The PUMA project I. Properties of the survey and first MUSE data results

Ultraluminous infrared galaxies (ULIRGs) are characterised by extreme starburst (SB) and AGN activity, and are therefore ideal laboratories for studying the outflow phenomena. We have recently started a project called Physics of ULIRGs with MUSE and ALMA (PUMA), which is a survey of 25 nearby (z < 0.165) ULIRGs observed with the integral field spectrograph MUSE and the interferometer ALMA. This sample includes systems with both AGN and SB nuclear activity in the pre- and post-coalescence phases of major mergers. The main goals of the project are to study the prevalence of multi-phase outflows as a function of the galaxy properties, to constrain the driving mechanisms of the outflows (e.g. distinguish between SB and AGN winds), and to identify feedback effects on the host galaxy. In this first paper, we present details on the sample selection, MUSE observations, and derive first data products. MUSE data were analysed to study the dynamical status of each of the 21 ULIRGs observed so far, taking the stellar kinematics and the morphological properties inferred from MUSE narrow-band images into account. We also located the ULIRG nuclei, using near-IR (HST) and mm (ALMA) data, and studied their optical spectra to infer the ionisation state through BPT diagnostics, and outflows in both ionised and neutral gas. We show that the morphological and stellar kinematic classifications are consistent: post-coalescence systems are more likely associated with ordered motions, while interacting (binary) systems are dominated by non-ordered and streaming motions. We also find broad and asymmetric [OIII] and NaID profiles in almost all nuclear spectra, with line widths in the range 300-2000 km/s, possibly associated with AGN- and SB-driven winds. This result reinforces previous findings that indicated that outflows are ubiquitous during the pre- and post-coalescence phases of major mergers.

preprint2020arXiv

SDSS IV MaNGA -- Metallicity and ionisation parameter in local star-forming galaxies from Bayesian fitting to photoionisation models

We measured gas-phase metallicity, ionisation parameter and dust extinction for 1795 representative local star-forming galaxies using integral field spectroscopy from the SDSS-IV MaNGA survey. We self-consistently derive these quantities by comparing observed line fluxes with photoionisation models using a Bayesian framework. We also present the first comprehensive study of the [SIII]$λλ$9069,9532 nebular lines, which have long been predicted to be ideal tracers of the ionisation parameter. Unfortunately, we find that current photoionisation models substantially over-predict [SIII] lines intensity, while broadly reproducing other observed optical line ratios. We discuss how to nonetheless make use of the information provided by [SIII] lines by setting a prior on the ionisation parameter. Following this approach, we derive spatially-resolved maps and radial profiles of metallicity and ionisation parameter. The metallicity radial profiles are comparable with previous works, with metallicity declining toward the outer parts and a flattening in the central regions, in agreement with infall models of galaxy formation, that predict that spiral discs build up through accretion of material, which leads to an inside-out growth. On the other hand, ionisation parameter radial profiles are flat for low-mass galaxies, while their slope becomes positive as galaxy mass increases. However, the ionisation parameter maps we obtain are clumpy, especially for low-mass galaxies. Ionisation parameter is tightly correlated with the H$α$ equivalent width [EW(H$α$)], following a nearly universal relation, which we attribute to the change of the spectral shape of ionising sources due to ageing of HII regions. We derive a positive correlation between ionisation parameter and metallicity at fixed EW(H$α$), in disagreement with previous theoretical works expecting an anti-correlation.

preprint2020arXiv

The ALMA view of the high-redshift relation between supermassive black holes and their host galaxies

This work aims at studying the $M_{BH}-M_{dyn}$ relation of a sample of $2<z<7$ quasars by constraining their host galaxy masses through full kinematical modeling of the cold gas kinematics, thus avoiding all possible biases and effects introduced by the rough virial estimates usually adopted so far. For this purpose we retrieved public observations of $72$ quasar host galaxies observed in ${\rm [CII]_{158μm}}$ or ${\rm CO}$ transitions with the Atacama Large Millimeter Array (ALMA). We then selected those quasars whose line emission is spatially resolved and performed a kinematic analysis on ALMA observations. We estimated the dynamical mass of the systems by modeling the gas kinematics with a rotating disc taking into account geometrical and instrumental effects. Our dynamical mass estimates, combined with $M_{BH}$ obtained from literature and our own new ${\rm CIV}\lambda1550$ observations, have allowed us to investigate the $ M_{BH}/M_{dyn}$ in the early Universe. Overall we obtained a sample of $10$ quasars at $z\sim2-7$ in which line emission is detected with high S/N ($> 5-10$) and the gas kinematics is spatially resolved and dominated by ordered rotation. The estimated dynamical masses place $6$ out of $10$ quasars above the local relation yielding to a $M_{BH}/M_{dyn}$ ratios $\sim10\times$ higher than those estimated in low-$z$ galaxies. On the other hand, we found that $4$ quasars at $z\sim 4-6$ have dynamical-to-BH mass ratios consistent with what is observed in early-type galaxies in the local Universe.

preprint2020arXiv

The ALPINE-ALMA [CII] Survey: Size of Individual Star-Forming Galaxies at z=4-6 and their Extended Halo Structure

We present the physical extent of [CII] 158um line-emitting gas from 46 star-forming galaxies at z=4-6 from the ALMA Large Program to INvestigate CII at Early Times (ALPINE). Using exponential profile fits, we measure the effective radius of the [CII] line (r_e,[CII]) for individual galaxies and compare them with the rest-frame ultra-violet (UV) continuum (r_e,UV) from Hubble Space Telescope images. The effective radius r_e,[CII] exceeds r_e,UV by factors of ~2-3 and the ratio of r_e,[CII]/r_e,UV increases as a function of M_star. We do not find strong evidence that [CII] line, the rest-frame UV, and FIR continuum are always displaced over ~ 1-kpc scale from each other. We identify 30% of isolated ALPINE sources as having an extended [CII] component over 10-kpc scales detected at 4.1$σ$-10.9$σ$ beyond the size of rest-frame UV and far-infrared (FIR) continuum. One object has tentative rotating features up to ~10-kpc, where the 3D model fit shows the rotating [CII]-gas disk spread over 4 times larger than the rest-frame UV-emitting region. Galaxies with the extended [CII] line structure have high star-formation rate (SFR), stellar mass (M_star), low Lya equivalent-width, and more blue-shifted (red-shifted) rest-frame UV metal absorption (Lya line), as compared to galaxies without such extended [CII] structures. Although we cannot rule out the possibility that a selection bias towards luminous objects may be responsible for such trends, the star-formation driven outflow also explains all these trends. Deeper observations are essential to test whether the extended [CII] line structures are ubiquitous to high-z star-forming galaxies.

preprint2020arXiv

The WISSH QSOs project IX. Cold gas content and environment of luminous QSOs at z~2.4-4.7

Sources at the brightest end of QSO luminosity function during the peak epoch of star formation and black hole accretion (z~2-4, i.e. Cosmic noon) are privileged sites to study the feeding & feedback cycle of massive galaxies. We perform the first systematic study of cold gas properties in the most luminous QSOs, by characterising their host-galaxies and environment. We analyse ALMA, NOEMA and JVLA observations of FIR continuum, CO and [CII] emission lines in eight QSOs ($L_{\rm Bol}>3\times10^{47}$ erg/s) from the WISSH sample at z~2.4-4.7. We report a 100% emission line detection rate and a 80% detection rate in continuum emission, and we find CO emission to be consistent with the steepest CO ladders observed so far. Sub-mm data reveal presence of (one or more) bright companion galaxies around 80% of WISSH QSOs, at projected distances of 6-130 kpc. We observe a variety of sizes for the molecular gas reservoirs (1.7-10 kpc), associated with rotating disks with disturbed kinematics. WISSH QSOs typically show lower CO luminosity and higher star formation efficiency than FIR matched, z~0-3 main-sequence galaxies, implying that, given the observed SFR ~170-1100 $M_\odot$/yr, molecular gas is converted into stars on <50 Myr. Most targets show extreme dynamical to black-hole mass ratios $M_{\rm dyn}/M_{\rm BH}\sim3-10$, two orders of magnitude smaller than local relations. The molecular gas fraction in WISSH hosts is lower by a factor of ~10-100 than in star forming galaxies with similar $M_*$. WISSH QSOs undergo an intense growth phase of both the central SMBH and host-galaxy. They pinpoint high-density sites where giant galaxies assemble and mergers play a major role in the build-up of the final host-galaxy mass. The observed low molecular gas fraction and short depletion timescale are likely due to AGN feedback, as traced by fast AGN-driven ionised outflows in all our targets.

preprint2020arXiv

Universal bolometric corrections for AGN over 7 luminosity decades

The AGN bolometric correction is a key element to understand BH demographics and compute accurate BH accretion histories from AGN luminosities. However, current estimates still differ from each other by up to a factor of two to three, and rely on extrapolations at the lowest and highest luminosities. Here we revisit this fundamental issue presenting general hard X-ray ($K_{X}$) and optical ($K_{O}$) bolometric corrections, computed combining several AGN samples spanning the widest (about 7 dex) luminosity range ever used for this kind of studies. We analysed a total of $\sim 1000$ type 1 and type 2 AGN for which a dedicated SED-fitting has been carried out. We provide a bolometric correction separately for type 1 and type 2 AGN; the two bolometric corrections results to be in agreement in the overlapping luminosity range and therefore, for the first time, a universal bolometric correction for the whole AGN sample (both type 1 and type 2) has been computed. We found that $K_{X}$ is fairly constant at $log(L_{BOL}/L_{\odot}) < 11$, while it increases up to about one order of magnitude at $log(L_{BOL}/L_{\odot}) \sim 14.5$. A similar increasing trend has been observed when its dependence on either the Eddington ratio or the BH mass is considered, while no dependence on redshift up to $z\sim3.5$ has been found. On the contrary, the optical bolometric correction appears to be fairly constant (i.e. $K_{O} \sim 5$) whatever is the independent variable. We also verified that our bolometric corrections correctly predict the AGN bolometric luminosity functions. According to this analysis, our bolometric corrections can be applied to the whole AGN population in a wide range of luminosity and redshift.

preprint2019arXiv

AGN feedback in a galaxy merger: Multi-phase, galaxy-scale outflows including a fast molecular gas blob ~6 kpc away from IRAS F08572+3915

To understand the role that AGN feedback plays in galaxy evolution we need in-depth studies of the multi-phase structure and energetics of galaxy-wide outflows. In this work we present new, deep ($\sim$50 hr) NOEMA CO(1-0) line observations of the molecular gas in the powerful outflow driven by the AGN in the ultra-luminous infrared galaxy IRAS F08572+3915. We spatially resolve the outflow, finding that its most likely configuration is a wide-angle bicone aligned with the kinematic major axis of the rotation disk. The molecular gas in the wind reaches velocities up to approximately $\pm$1200 km s$^{-1}$ and transports nearly 20% of the molecular gas mass in the system. We detect a second outflow component located $\sim$6 kpc north-west from the galaxy moving away at $\sim$900 km s$^{-1}$, which could be the result of a previous episode of AGN activity. The total mass and energetics of the outflow, which includes contributions from the ionized, neutral, warm and cold molecular gas phases is strongly dominated by the cold molecular gas. In fact, the molecular mass outflow rate is higher than the star formation rate, even if we only consider the gas in the outflow that is fast enough to escape the galaxy, which accounts for about $\sim$40% of the total mass of the outflow. This results in an outflow depletion time for the molecular gas in the central $\sim$1.5 kpc region of only $\sim3$ Myr, a factor of $\sim2$ shorter than the depletion time by star formation activity.

preprint2019arXiv

Molecular Gas Inflows and Outflows in Ultraluminous Infrared Galaxies at $z\sim0.2$ and one QSO at $z=6.1$

Aims. We aim to search and characterize inflows and outflows of molecular gas in four ultraluminous infrared galaxies (ULIRGs) at $z\sim0.2-0.3$ and one distant QSO at $z=6.13$. Methods. We use Herschel PACS and ALMA Band 7 observations of the hydroxyl molecule (OH) line at rest-frame wavelength 119 $μ$m which in absorption can provide unambiguous evidence for inflows or outflows of molecular gas in nuclear regions of galaxies. Our study contributes to double the number of OH observations of luminous systems at $z\sim0.2-0.3$, and push the search for molecular outflows based on the OH transition to $z\sim6$. Results. We detect OH high-velocity absorption wings in three of the four ULIRGs. In two cases, IRAS F20036-1547 and IRAS F13352+6402, the blueshifted absorption profiles indicate the presence of powerful and fast molecular gas outflows. Consistent with an inside-out quenching scenario, these outflows are depleting the central reservoir of molecular gas at a similar rate than the intense star formation activity. In the case of the starburst-dominated system IRAS 10091+4704, we detect an inverted P-Cygni profile that is unique among ULIRGs and indicates the presence of a fast ($\sim400$ km s$^{-1}$) inflow of molecular gas at a rate of $\sim100~M_{\odot}~{\rm yr}^{-1}$ towards the central region. Finally, we tentatively detect ($\sim3σ$) the OH doublet in absorption in the $z=6.13$ QSO ULAS J131911+095051. The OH feature is blueshifted with a median velocity that suggests the presence of a molecular outflow, although characterized by a modest molecular mass loss rate of $\sim200~M_{\odot}~{\rm yr}^{-1}$. This value is comparable to the small mass outflow rates found in the stacking of the [CII] spectra of other $z\sim6$ QSOs and suggests that ejective feedback in this phase of the evolution of ULAS J131911+095051 has subsided.

preprint2019arXiv

Molecular outflows in local galaxies: Method comparison and a role of intermittent AGN driving

We report new detections and limits from a NOEMA and ALMA CO(1-0) search for molecular outflows in 13 local galaxies with high FIR surface brightness, and combine with results from the literature. CO line ratios and outflow structure provide some constraints on the conversion from observables to quantities such as molecular mass outflow rates. Ratios between outflow emission in higher J CO transitions and in CO(1-0) typically are consistent with excitation Ri1<~1. For IRAS 13120-5453, however, R31=2.10 indicates optically thin CO in the outflow. Like much of the outflow literature, we use alpha(CO) = 0.8, and we present arguments for using C=1 in deriving molecular mass outflow rates Mdot = C*M*v/R. We compare the two main methods for molecular outflow detection: CO mm interferometry and Herschel OH spectroscopy. For 26 sources studied with both methods, we find 80% agreement in detecting vout>~150km/s outflows, and non-matches can be plausibly ascribed to outflow geometry and SNR. For 12 bright ULIRGs with detailed OH-based outflow modeling, CO outflows are detected in all but one. Outflow masses, velocities, and sizes for these 11 sources agree well between the two methods, and modest remaining differences may relate to the different but overlapping regions sampled by CO emission and OH absorption. Outflow properties correlate better with AGN luminosity and with bolometric luminosity than with FIR surface brightness. The most massive outflows are found for systems with current AGN activity, but significant outflows in non-AGN systems must relate to star formation or to AGN activity in the recent past. We report scaling relations for the increase of outflow mass, rate, momentum rate, and kinetic power with bolometric luminosity. Short ~10^6yr flow times and some sources with resolved multiple outflow episodes support a role of intermittent driving, likely by AGN. (abridged)

preprint2019arXiv

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

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

preprint2016arXiv

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

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

preprint2016arXiv

Detection of faint broad emission lines in type 2 AGN: I. Near infrared observations and spectral fitting

We present medium resolution near infrared spectroscopic observations of 41 obscured and intermediate class AGN (type 2, 1.9 and 1.8; AGN2) with redshift $z \lesssim$0.1, selected from the Swift/BAT 70-month catalogue. The observations have been carried out in the framework of a systematic study of the AGN2 near infrared spectral properties and have been executed using ISAAC/VLT, X-shooter/VLT and LUCI/LBT, reaching an average S/N ratio of $\sim$30 per resolution element. For those objects observed with X-shooter we also obtained simultaneous optical and UV spectroscopy. We have identified a component from the broad line region in 13 out of 41 AGN2, with FWHM ${\rm > 800 }$ km/s. We have verified that the detection of the broad line region components does not significantly depend on selection effects due to the quality of the spectra, the X-ray or near infrared fluxes, the orientation angle of the host galaxy or the hydrogen column density measured in the X-ray band. The average broad line region components found in AGN2 has a significantly (a factor 2) smaller FWHM if compared with a control sample of type 1 AGN.

preprint2016arXiv

EELT-HIRES the high-resolution spectrograph for the E-ELT

The first generation of E-ELT instruments will include an optical-infrared High Resolution Spectrograph, conventionally indicated as EELT-HIRES, which will be capable of providing unique breakthroughs in the fields of exoplanets, star and planet formation, physics and evolution of stars and galaxies, cosmology and fundamental physics. A 2-year long phase A study for EELT-HIRES has just started and will be performed by a consortium composed of institutes and organisations from Brazil, Chile, Denmark, France, Germany, Italy, Poland, Portugal, Spain, Sweden, Switzerland and United Kingdom. In this paper we describe the science goals and the preliminary technical concept for EELT-HIRES which will be developed during the phase A, as well as its planned development and consortium organisation during the study.

preprint2016arXiv

Fast outflows and star formation quenching in quasar host galaxies

Negative feedback from active galactic nuclei (AGN) is considered a key mechanism in shaping galaxy evolution. Fast, extended outflows are frequently detected in the AGN host galaxies at all redshifts and luminosities, both in ionised and molecular gas. However, these outflows are only "potentially" able to quench star formation and we are still missing a decisive evidence of negative feedback in action. Here we present Spectrograph for INtegral Field Observations in the Near Infrared (SINFONI) H- and K-band integral-field spectroscopic observations of two quasars at $z\sim$2.4 characterised by fast, extended outflows detected through the [OIII]$λ$5007 line (Carniani et al. 2015). The high signal-to-noise ratio of our observations allows us to identify faint narrow (FWHM $< 500$ km/s), and spatially extended components in [OIII]$λ$5007 and H$α$ emission associated with star formation in the host galaxy. Such star-formation powered emission is spatially anti-correlated with the fast outflows. The ionised outflows therefore appear to be able to suppress star formation in the region where the outflow is expanding. However the detection of narrow, spatially-extended H$α$ emission indicates star formation rates of at least $\sim 50-90$ M$_{\odot}$/yr, suggesting either that AGN feedback does not affect the whole galaxy or that many feedback episodes are required before star formation is completely quenched. On the other hand, the narrow H$α$ emission extending along the edges of the outflow cone may lead also to a positive feedback interpretation. Our results highlight the possible twofold role of galaxy-wide outflows in host galaxy evolution.

preprint2016arXiv

First observational support for overlapping reionized bubbles generated by a galaxy overdensity

We present the analysis of deep HST multi-band imaging of the BDF field specifically designed to identify faint companions around two of the few Ly-alpha emitting galaxies spectroscopically confirmed at z~7 (Vanzella et al. 2011). Although separated by only 4.4 proper Mpc these galaxies cannot generate HII regions large enough to explain visibility of their Ly-alpha line, thus requiring a population of fainter ionizing sources in their vicinity. We use deep HST and VLT-Hawk-I data to select z~7 Lyman break galaxies around the emitters. We select 6 new robust z~7 LBGs at Y~26.5-27.5 whose average spectral energy distribution is consistent with the objects being at the redshift of the close-by Ly-alpha emitters. The resulting number density of z~7 LBGs in the BDF field is a factor ~3-4 higher than expected in random pointings of the same size. We compare these findings with cosmological hydrodynamic plus radiative transfer simulations of a universe with a half neutral IGM: we find that indeed Ly-alpha emitter pairs are only found in completely ionized regions characterized by significant LBG overdensities. Our findings match the theoretical prediction that the first ionization fronts are generated within significant galaxy overdensities and support a scenario where faint, "normal" star-forming galaxies are responsible for reionization.

preprint2016arXiv

Galaxy metallicities depend primarily on stellar mass and molecular gas mass

In this work we present an analysis of the behaviour of galaxies in a four-dimensional parameter space defined by stellar mass, metallicity, star formation rate, and molecular gas mass. We analyse a combined sample of 227 galaxies, which draws from a number of surveys across the redshift range 0 < z < 2 (> 90% of the sample at z~0), and covers > 3 decades in stellar mass.Using Principle Component Analysis, we demonstrate that galaxies in our sample lie on a 2-dimensional plane within this 4D parameter space, indicative of galaxies that exist in an equilibrium between gas inflow and outflow. Furthermore, we find that the metallicity of galaxies depends only on stellar mass and molecular gas mass. In other words, gas-phase metallicity has a negligible dependence on star formation rate, once the correlated effect of molecular gas content is accounted for. The well-known `fundamental metallicity relation', which describes a close and tight relationship between metallicity and SFR (at fixed stellar mass) is therefore entirely a by-product of the underlying physical relationship with molecular gas mass (via the Schmidt-Kennicutt relation).

preprint2016arXiv

Neutral gas outflows in nearby [U]LIRGs via optical NaD feature

We studied the properties of the neutral gas in a sample of 38 local [U]LIRGs, which mainly covers the less explored LIRG luminosity range. This study is based on the analysis of the spatially integrated and spatially resolved spectra of the NaD feature obtained with the VIMOS/VLT-IFS. Analyzing spatially integrated spectra, we find that the contribution of stars to the observed NaD-EWs is generally small (<35%), and therefore this feature is dominated by ISM absorption. After subtracting the stellar contribution, we find that the pure-ISM integrated spectra generally show blueshifted NaD profiles, indicating neutral gas outflow velocities V in the range 65-260 km/s. Excluding galaxies with AGNs, V shows a dependency with the SFR of the type V \propto SFR^0.15 (in rather good agreement with previous results). The spatially resolved analysis was performed for 40 galaxies, 22 of which have neutral gas velocity fields dominated by noncircular motions with signatures of cone-like winds. However, a modest number of targets (11/40) show disk rotation. We found that the wind masses are in the range (0.4-7.5)x10^8 Msun, reaching up to 3% of the dynamical mass of the host. The mass rates are typically only 0.2-0.4 times the SFR indicating that, the mass loss is too small to slow down the SF significantly. In the majority of cases, the velocity of the outflowing gas is not sufficient to escape the host potential well and, therefore, most of the gas rains back into the disk. On average V/v(esc) is higher in less massive galaxies, confirming that the galaxy mass has a primary role in shaping the recycling of gas and metals. The comparison between the wind power and kinetic power of the starburst associated with SNe indicates that the starburst could drive the outflows in nearly all the [U]LIRGs galaxies, as the wind power is generally lower than 20% of the kinetic power supplied by the starburst.

preprint2016arXiv

New fully empirical calibrations of strong-line metallicity indicators in star forming galaxies

We derive new empirical calibrations for strong-line diagnostics of gas phase metallicity in local star forming galaxies by uniformly applying the Te method over the full metallicity range probed by the Sloan Digital Sky Survey (SDSS). To measure electron temperatures at high metallicity, where the auroral lines needed are not detected in single galaxies, we stacked spectra of more than 110,000 galaxies from the SDSS in bins of log[O II]/H$β$ and log[O III]/H$β$. This stacking scheme does not assume any dependence of metallicity on mass or star formation rate, but only that galaxies with the same line ratios have the same oxygen abundance. We provide calibrations which span more than 1 dex in metallicity and are entirely defined on a consistent absolute Te metallicity scale for galaxies. We apply our calibrations to the SDSS sample and find that they provide consistent metallicity estimates to within 0.05 dex.

preprint2016arXiv

Tracing the reionization epoch with ALMA: [CII] emission in z~7 galaxies

We present new results on [CII]158$μ$ m emission from four galaxies in the reionization epoch. These galaxies were previously confirmed to be at redshifts between 6.6 and 7.15 from the presence of the Ly$α$ emission line in their spectra. The Ly$α$ emission line is redshifted by 100-200 km/s compared to the systemic redshift given by the [CII] line. These velocity offsets are smaller than what is observed in z~3 Lyman break galaxies with similar UV luminosities and emission line properties. Smaller velocity shifts reduce the visibility of Ly$α$ and hence somewhat alleviate the need for a very neutral IGM at z~7 to explain the drop in the fraction of Ly$α$ emitters observed at this epoch. The galaxies show [CII] emission with L[CII]=0.6-1.6 x10$^8 L_\odot$: these luminosities place them consistently below the SFR-L[CII] relation observed for low redshift star forming and metal poor galaxies and also below z =5.5 Lyman break galaxies with similar star formation rates. We argue that previous undetections of [CII] in z~7 galaxies with similar or smaller star formation rates are due to selection effects: previous targets were mostly strong Ly$α$ emitters and therefore probably metal poor systems, while our galaxies are more representative of the general high redshift star forming population .

preprint2015arXiv

ALMA constraints on the faint millimetre source number counts and their contribution to the cosmic infrared background

We have analysed 18 ALMA continuum maps in Bands 6 and 7, with rms down to 7.8$μ$Jy, to derive differential number counts down to 60$μ$Jy and 100$μ$Jy at $λ=$1.3 mm and $λ=$1.1 mm, respectively. The area covered by the combined fields is $\rm 9.5\times10^{-4}deg^2$ at 1.1mm and $\rm 6.6\times10^{-4}deg^{2}$ at 1.3mm. We improved the source extraction method by requiring that the dimension of the detected sources be consistent with the beam size. This method enabled us to remove spurious detections that have plagued the purity of the catalogues in previous studies. We detected 50 faint sources with S/N$>$3.5 down to 60$μ$Jy, hence improving the statistics by a factor of four relative to previous studies. The inferred differential number counts are $\rm dN/d(Log_{10}S)=1\times10^5~deg^2$ at a 1.1 mm flux $S_{λ= 1.1~mm} = 130~μ$Jy, and $\rm dN/d(Log_{10}S)=1.1\times10^5~deg^2$ at a 1.3 mm flux $\rm S_{λ= 1.3~mm} = 60~μ$Jy. At the faintest flux limits, i.e. 30$μ$Jy and 40$μ$Jy, we obtain upper limits on the differential number counts of $\rm dN/d(Log_{10}S) < 7\times10^5~deg^2$ and $\rm dN/d(Log_{10}S)<3\times10^5~deg^2$, respectively. Our results provide a new lower limit to CIB intensity of 17.2${\rm Jy\ deg^{-2}}$ at 1.1mm and of 12.9${\rm Jy\ deg^{-2}}$ at 1.3mm. Moreover, the flattening of the integrated number counts at faint fluxes strongly suggests that we are probably close to the CIB intensity. Our data imply that galaxies with SFR$<40~M_{\odot}/yr$ certainly contribute less than 50% to the CIB while more than 50% of the CIB must be produced by galaxies with $\rm SFR>40~M_{\odot}/yr$. The differential number counts are in nice agreement with recent semi-analytical models of galaxy formation even as low as our faint fluxes. Consequently, this supports the galaxy evolutionary scenarios and assumptions made in these models.

preprint2015arXiv

Evidence for feedback in action from the molecular gas content in the z~1.6 outflowing QSO XID2028

Gas outflows are believed to play a pivotal role in shaping galaxies, as they regulate both star formation and black hole growth. Despite their ubiquitous presence, the origin and the acceleration mechanism of such powerful and extended winds is not yet understood. Direct observations of the cold gas component in objects with detected outflows at other wavelengths are needed to assess the impact of the outflow on the host galaxy interstellar medium (ISM). We observed with the Plateau de Bure Interferometer an obscured quasar at z~1.5, XID2028, for which the presence of an ionised outflow has been unambiguously signalled by NIR spectroscopy. The detection of CO(3-2) emission in this source allows us to infer the molecular gas content and compare it to the ISM mass derived from the dust emission. We then analyze the results in the context of recent insights on scaling relations, which describe the gas content of the overall population of star-forming galaxies at a similar redshifts. The Star formation efficiency (~100) and gas mass (M_gas=2.1-9.5x10^{10} M_sun) inferred from the CO(3-2) line depend on the underlying assumptions on the excitation of the transition and the CO-to-H2 conversion factor. However, the combination of this information and the ISM mass estimated from the dust mass suggests that the ISM/gas content of XID2028 is significantly lower than expected for its observed M$_\star$, sSFR and redshift, based on the most up-to-date calibrations (with gas fraction <20% and depletion time scale <340 Myr). Overall, the constraints we obtain from the far infrared and millimeter data suggest that we are observing QSO feedback able to remove the gas from the host

preprint2015arXiv

Galaxy gas flows inferred from a detailed, spatially resolved metal budget

We use the most extensive integral field spectroscopic map of a local galaxy, NGC 628, combined with gas and stellar mass surface density maps, to study the distribution of metals in this galaxy out to 3 effective radii ($\rm R_e$). At each galactocentric distance, we compute the metal budget and thus constrain the mass of metals lost. We find that in the disc about 50% of the metals have been lost throughout the lifetime of the galaxy. The fraction of metals lost is higher in the bulge ($\sim$70%) and decreases towards the outer disc ($\rm \sim 3 \ R_e$). In contrast to studies based on the gas kinematics, which are only sensitive to ongoing outflow events, our metal budget analysis enables us to infer the average outflow rate during the galaxy lifetime. By using simple physically motivated models of chemical evolution we can fit the observed metal budget at most radii with an average outflow loading factor of order unity, thus clearly demonstrating the importance of outflows in the evolution of disc galaxies of this mass range ($\rm log(M_\star/M_\odot) \sim 10)$. The observed gas phase metallicity is higher than expected from the metal budget and suggests late-time accretion of enriched gas, likely raining onto the disc from the metal-enriched halo.

preprint2015arXiv

Ionised outflows in z $\sim$ 2.4 quasar host galaxies

AGN-driven outflows are invoked by galaxy evolutionary models to quench star formation and to explain the origin of the relations observed locally between super massive black holes and their host galaxies. This work aims to detect the presence of extended ionised outflows in luminous quasars where we expect the maximum activity both in star formation and in black hole accretion. Currently, there are only a few studies based on spatially resolved observations of outflows at high redshift, $z>2$. We analyse a sample of six luminous (${\rm L>10^{47} \ erg/s}$) quasars at $z\sim2.4$, observed in H-band using the near-IR integral field spectrometer SINFONI at VLT. We perform a kinematic analysis of the [OIII] emission line at $λ= 5007Å$. [OIII] has a complex gas kinematic, with blue-shifted velocities of a few hundreds of km/s and line widths up to 1500 km/s. Using the spectroastrometric method we infer size of the ionised outflows of up to $\sim$2 kpc. The properties of the ionised outflows, mass outflow rate, momentum rate and kinetic power, are correlated with the AGN luminosity. The increase in outflow rate with increasing AGN luminosity is consistent with the idea that a luminous AGN pushes away the surrounding gas through fast outflows driven by radiation pressure, which depends on the emitted luminosity. We derive mass outflow rates of about 6-700 M$_{\odot}$/yr for our sample, which are lower than those observed in molecular outflows. Indeed physical properties of ionised outflows show dependences on AGN luminosity which are similar to those of molecular outflows but indicating that the mass of ionised gas is smaller than that of the molecular one. Alternatively, this discrepancy between ionised and molecular outflows could be explained with different acceleration mechanisms.

preprint2015arXiv

Is there any evidence that ionised outflows quench star formation in type 1 quasars at z<1?

The aim of this paper is to test the basic model of negative AGN feedback. According to this model, once the central black hole accretes at the Eddington limit and reaches a certain critical mass, AGN driven outflows blow out gas, suppressing star formation in the host galaxy and self-regulating black hole growth. We consider a sample of 224 quasars selected from the SDSS at z<1 observed in the infrared band by Herschel. We evaluate the star formation rate in relation to several outflow signatures traced by the [OIII]4959,5007 and [OII]3726,3729 emission lines in about half of the sample with high quality spectra. Most of the quasars show asymmetric and broad wings in [OIII], which we interpret as outflow signatures. We separate the quasars in two groups, ``weakly'' and ``strongly'' outflowing, using three different criteria. When we compare the mean star formation rate in five redshift bins in the two groups, we find that the SFRs are comparable or slightly larger in the strongly outflowing quasars. We estimate the stellar mass from SED fitting and the quasars are distributed along the star formation main sequence, although with a large scatter. The scatter from this relation is uncorrelated with respect to the kinematic properties of the outflow. Moreover, for quasars dominated in the infrared by starburst or by AGN emission, we do not find any correlation between the star formation rate and the velocity of the outflow, a trend previously reported in the literature for pure starburst galaxies. We conclude that the basic AGN negative feedback scenario seems not to agree with our results. Although we use a large sample of quasars, we did not find any evidence that the star formation rate is suppressed in the presence of AGN driven outflows on large scale. A possibility is that feedback is effective over much longer timescales than those of single episodes of quasar activity.

preprint2015arXiv

Mapping metals at high redshift with far-infrared lines

Cosmic metal enrichment is one of the key physical processes regulating galaxy formation and the evolution of the intergalactic medium (IGM). However, determining the metal content of the most distant galaxies has proven so far almost impossible; also, absorption line experiments at $z\sim6$ become increasingly difficult because of instrumental limitations and the paucity of background quasars. With the advent of ALMA, far-infrared emission lines provide a novel tool to study early metal enrichment. Among these, the [CII] line at 157.74 $μ$m is the most luminous line emitted by the interstellar medium of galaxies. It can also resonant scatter CMB photons inducing characteristic intensity fluctuations ($ΔI/I_{CMB}$) near the peak of the CMB spectrum, thus allowing to probe the low-density IGM. We compute both [CII] galaxy emission and metal-induced CMB fluctuations at $z\sim 6$ by using Adaptive Mesh Refinement cosmological hydrodynamical simulations and produce mock observations to be directly compared with ALMA BAND6 data ($ν_{obs}\sim 272$ GHz). The [CII] line flux is correlated with $M_{UV}$ as $\log(F_{peak}/μ{\rm Jy})=-27.205-2.253\,M_{UV}-0.038\,M_{UV}^2$. Such relation is in very good agreement with recent ALMA observations (e.g. Maiolino et al. 2015; Capak et al. 2015) of $M_{UV}<-20$ galaxies. We predict that a $M_{UV}=-19$ ($M_{UV}=-18$) galaxy can be detected at $4σ$ in $\simeq40$ (2000) hours, respectively. CMB resonant scattering can produce $\simeq\pm 0.1\,μ$Jy/beam emission/absorptions features that are very challenging to be detected with current facilities. The best strategy to detect these signals consists in the stacking of deep ALMA observations pointing fields with known $M_{UV}\simeq-19$ galaxies. This would allow to simultaneously detect both [CII] emission from galactic reionization sources and CMB fluctuations produced by $z\sim6$ metals.

preprint2015arXiv

Molecular depletion times and the CO-to-H2 conversion factor in metal-poor galaxies

Tracing molecular hydrogen content with carbon monoxide in low-metallicity galaxies has been exceedingly difficult. Here we present a new effort, with IRAM 30-m observations of 12CO(1-0) of a sample of 8 dwarf galaxies having oxygen abundances ranging from 12+logO/H=7.7 to 8.4. CO emission is detected in all galaxies, including the most metal-poor galaxy of our sample (0.1 Zsun); to our knowledge this is the largest number of 12CO(1-0) detections ever reported for galaxies with 12+logO/H<=8 (0.2 Zsun) outside the Local Group. We calculate stellar masses (Mstar) and star-formation rates (SFRs), and analyze our results by combining our observations with galaxy samples from the literature. Extending previous results for a correlation of the molecular gas depletion time, tau(dep), with Mstar and specific SFR (sSFR), we find a variation in tau(dep) of a factor of 200 or more (from <50 Myr to 10 Gyr) over a spread of 1000 in sSFR and Mstar. We exploit the variation of tau(dep) to constrain the CO-to-H2 mass conversion factor alpha(CO) at low metallicity, and assuming a power-law variation find alpha(CO) \propto (Z/Zsun)^1.9, similar to results based on dust continuum measurements compared with gas mass. By including HI measurements, we show that the fraction of total gas mass relative to the baryonic mass is higher in galaxies that are metal poor, of low mass, and of high sSFR. Finally, comparisons of the data with star-formation models of the molecular gas phases suggest that, at metallicities Z/Zsun<=0.2, there are some discrepancies with model predictions.

preprint2015arXiv

Molecular gas as the driver of fundamental galactic relations

There has been much recent work dedicated to exploring secondary correlations in the mass-metallicity relation, with significant dependence on both the star formation rate and HI content being demonstrated. Previously, a paucity of molecular gas data (combined with sample selection bias) hampered the investigation of any such relation with molecular gas content. In this work, we assemble a sample of 221 galaxies from a variety of surveys in the redshift range 0 < z < 2, to explore the connection between molecular gas content and metallicity. We explore the effect of gas mass on the mass-metallicity relation, finding that the offset from the relation is negatively correlated against both molecular and total gas mass. We then employ a principle component analysis technique to explore secondary dependences in the mass-metallicity relation, finding that the secondary dependence with gas mass is significantly stronger than with star formation rate, and as such the underlying `Fundamental Metallicity Relation' is between stellar mass, metallicity, and gas mass. In particular, the metallicity dependence on SFR is simply a byproduct of the dependence on the molecular gas content, via the Schmidt-Kennicutt relation. Finally, we note that our principle component analysis finds essentially no connection between gas-phase metallicity and the efficiency of star formation.

preprint2015arXiv

New observations of z~7 galaxies: evidence for a patchy reionization

We present new results from our search for z~7 galaxies from deep spectroscopic observations of candidate z-dropouts in the CANDELS fields. Despite the extremely low flux limits achieved by our sensitive observations, only 2 galaxies have robust redshift identifications, one from its Lyalpha emission line at z=6.65, the other from its Lyman-break, i.e. the continuum discontinuity at the Lyalpha wavelength consistent with a redshift 6.42, but with no emission line. In addition, for 23 galaxies we present deep limits in the Lyalpha EW derived from the non detections in ultra-deep observations. Using this new data as well as previous samples, we assemble a total of 68 candidate z~7 galaxies with deep spectroscopic observations, of which 12 have a line detection. With this much enlarged sample we can place solid constraints on the declining fraction of Ly$α$ emission in z~7 Lyman break galaxies compared to z~6, both for bright and faint galaxies. Applying a simple analytical model, we show that the present data favor a patchy reionization process rather than a smooth one.

preprint2015arXiv

P-MaNGA: Emission Lines Properties - Gas Ionisation and Chemical Abundances from Prototype Observations

MaNGA (Mapping Nearby Galaxies at Apache Point Observatory) is a SDSS-IV survey that will obtain spatially resolved spectroscopy from 3600 Å to 10300 Å for a representative sample of over 10000 nearby galaxies. In this paper we present the analysis of nebular emission line properties in 14 galaxies obtained with P-MaNGA, a prototype of the MaNGA instrument. Using spatially resolved diagnostic diagrams we find extended star formation in galaxies that are centrally dominated by Seyfert/LINER-like emission, which illustrates that galaxy characterisations based on single fibre spectra are necessarily incomplete. We observe extended LINER-like emission (up to $\rm 1 R_{e}$) in three galaxies. We make use of the $\rm EW(H α)$ to argue that the observed emission is consistent with ionisation from hot evolved stars. We derive stellar population indices and demonstrate a clear correlation between $\rm D_n(4000)$ and $\rm EW(H δ_A)$ and the position in the ionisation diagnostic diagram: resolved galactic regions which are ionised by a Seyfert/LINER-like radiation field are also devoid of recent star formation and host older and/or more metal rich stellar populations. We also detect extraplanar LINER-like emission in two highly inclined galaxies, and identify it with diffuse ionised gas. We investigate spatially resolved metallicities and find a positive correlation between metallicity and star formation rate (SFR) surface density. We study the relation between N/O vs O/H on resolved scales. We find that, at given N/O, regions within individual galaxies are spread towards lower metallicities, deviating from the sequence defined by galactic central regions as traced by Sloan $3''$ fibre spectra. We suggest that the observed dispersion can be a tracer for gas flows in galaxies: infalls of pristine gas and/or the effect of a galactic fountain.

preprint2015arXiv

Predictions for surveys with the SPICA Mid-infrared Instrument

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

preprint2015arXiv

The assembly of "normal" galaxies at z=7 probed by ALMA

We report new deep ALMA observations aimed at investigating the [CII]158um line and continuum emission in three spectroscopically confirmed Lyman Break Galaxies at 6.8<z<7.1, i.e. well within the re-ionization epoch. With Star Formation Rates of SFR ~ 5-15 Msun/yr these systems are much more representative of the high-z galaxy population than other systems targeted in the past by millimeter observations. For the galaxy with the deepest observation we detect [CII] emission at redshift z=7.107, fully consistent with the Lyalpha redshift, but spatially offset by 0.7" (4 kpc) from the optical emission. At the location of the optical emission, tracing both the Lyalpha line and the far-UV continuum, no [CII] emission is detected in any of the three galaxies, with 3sigma upper limits significantly lower than the [CII] emission observed in lower reshift galaxies. These results suggest that molecular clouds in the central parts of primordial galaxies are rapidly disrupted by stellar feedback. As a result, [CII] emission mostly arises from more external accreting/satellite clumps of neutral gas. These findings are in agreement with recent models of galaxy formation. Thermal far-infrared continuum is not detected in any of the three galaxies. However, the upper limits on the infrared-to-UV emission ratio do not exceed those derived in metal- and dust-poor galaxies.

preprint2015arXiv

The MAGNUM survey: Positive feedback in the nuclear region of NGC 5643 suggested by MUSE

We study the ionization and kinematics of the ionized gas in the nuclear region of the barred Seyfert 2 galaxy NGC~5643 using MUSE integral field observations in the framework of the MAGNUM (Measuring Active Galactic Nuclei Under MUSE Microscope) survey. The data were used to identify regions with different ionization conditions and to map the gas density and the dust extinction. We find evidence for a double sided ionization cone, possibly collimated by a dusty structure surrounding the nucleus. At the center of the ionization cone, outflowing ionized gas is revealed as a blueshifted, asymmetric wing of the [OIII] emission line, up to projected velocity v(10)~-450 km/s. The outflow is also seen as a diffuse, low luminosity radio and X-ray jet, with similar extension. The outflowing material points in the direction of two clumps characterized by prominent line emission with spectra typical of HII regions, located at the edge of the dust lane of the bar. We propose that the star formation in the clumps is due to `positive feedback' induced by gas compression by the nuclear outflow, providing the first candidate for outflow induced star formation in a Seyfert-like radio quiet AGN. This suggests that positive feedback may be a relevant mechanism in shaping the black hole-host galaxy coevolution.

preprint2015arXiv

The multi-phase winds of Markarian 231: from the hot, nuclear, ultra-fast wind to the galaxy-scale, molecular outflow

We present the best sensitivity and angular resolution maps of the molecular disk and outflow of Mrk 231, as traced by CO observations obtained with IRAM/PdBI, and we analyze archival Chandra and NuSTAR observations. We constrain the physical properties of both the molecular disk and outflow, the presence of a highly-ionized ultra-fast nuclear wind, and their connection. The molecular outflow has a size of ~1 kpc, and extends in all directions around the nucleus, being more prominent along the south-west to north-east direction, suggesting a wide-angle biconical geometry. The maximum projected velocity of the outflow is nearly constant out to ~1 kpc, thus implying that the density of the outflowing material decreases from the nucleus outwards as $r^{-2}$. This suggests that either a large part of the gas leaves the flow during its expansion or that the bulk of the outflow has not yet reached out to ~1 kpc, thus implying a limit on its age of ~1 Myr. We find $\dot M_{OF}=[ 500-1000]~ M_{\odot}~yr^{-1}$ and $\dot E_{kin,OF}=[7-10]\times 10^{43}$ erg s$^{-1}$. Remarkably, our analysis of the X-ray data reveals a nuclear ultra-fast outflow (UFO) with velocity -20000 km s$^{-1}$, $\dot M_{UFO}=[0.3- 2.1] ~M_\odot yr^{-1}$, and momentum load $\dot P_{UFO}/\dot P_{rad}=[0.2-1.6]$.We find $\dot E_{kin,UFO}\sim \dot E_{kin,OF}$ as predicted for outflows undergoing an energy conserving expansion. This suggests that most of the UFO kinetic energy is transferred to mechanical energy of the kpc-scale outflow, strongly supporting that the energy released during accretion of matter onto super-massive black holes is the ultimate driver of giant massive outflows. We estimate a momentum boost $\dot P_{OF}/\dot P_{UFO}\approx [30-60]$. The ratios $\dot E_{kin, UFO}/L_{bol,AGN} =[ 1-5]\%$ and $\dot E_{kin,OF}/L_{bol,AGN} = [1-3]\%$ agree with the requirements of the most popular models of AGN feedback.

preprint2014arXiv

ALLSMOG: an APEX Low-redshift Legacy Survey for MOlecular Gas. I - molecular gas scaling relations, and the effect of the CO/H2 conversion factor

We present ALLSMOG, the APEX Low-redshift Legacy Survey for MOlecular Gas. ALLSMOG is a survey designed to observe the CO(2-1) emission line with the APEX telescope, in a sample of local galaxies (0.01 < z < 0.03), with stellar masses in the range 8.5 < log(M*/Msun) < 10. This paper is a data release and initial analysis of the first two semesters of observations, consisting of 42 galaxies observed in CO(2-1). By combining these new CO(2-1) emission line data with archival HI data and SDSS optical spectroscopy, we compile a sample of low-mass galaxies with well defined molecular gas masses, atomic gas masses, and gas-phase metallicities. We explore scaling relations of gas fraction and gas consumption timescale, and test the extent to which our findings are dependent on a varying CO/H2 conversion factor. We find an increase in the H2/HI mass ratio with stellar mass which closely matches semi-analytic predictions. We find a mean molecular gas fraction for ALLSMOG galaxies of MH2/M* = (0.09 - 0.13), which decreases with stellar mass. We measure a mean molecular gas consumption timescale for ALLSMOG galaxies of 0.4 - 0.7 Gyr. We also confirm the non-universality of the molecular gas consumption timescale, which varies (with stellar mass) from ~100 Myr to ~2 Gyr. Importantly, we find that the trends in the H2/HI mass ratio, gas fraction, and the non-universal molecular gas consumption timescale are all robust to a range of recent metallicity-dependent CO/H2 conversion factors.

preprint2014arXiv

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

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

preprint2014arXiv

Blowin' in the wind: both `negative' and `positive' feedback in an obscured high-z Quasar

Quasar feedback in the form of powerful outflows is invoked as a key mechanism to quench star formation in galaxies, preventing massive galaxies to over-grow and producing the red colors of ellipticals. On the other hand, some models are also requiring `positive' AGN feedback, inducing star formation in the host galaxy through enhanced gas pressure in the interstellar medium. However, finding observational evidence of the effects of both types of feedback is still one of the main challenges of extragalactic astronomy, as few observations of energetic and extended radiatively-driven winds are available. Here we present SINFONI near infrared integral field spectroscopy of XID2028, an obscured, radio-quiet z=1.59 QSO detected in the XMM-COSMOS survey, in which we clearly resolve a fast (1500 km/s) and extended (up to 13 kpc from the black hole) outflow in the [OIII] lines emitting gas, whose large velocity and outflow rate are not sustainable by star formation only. The narrow component of Ha emission and the rest frame U band flux from HST-ACS imaging enable to map the current star formation in the host galaxy: both tracers independently show that the outflow position lies in the center of an empty cavity surrounded by star forming regions on its edge. The outflow is therefore removing the gas from the host galaxy (`negative feedback'), but also triggering star formation by outflow induced pressure at the edges (`positive feedback'). XID2028 represents the first example of a host galaxy showing both types of feedback simultaneously at work.

preprint2014arXiv

Constraining the Nature of Two Ly-alpha Emitters detected by ALMA at z = 4.7

We report optical VLT FORS2 spectroscopy of the two Ly-alpha emitters (LAEs) companions to the quasi-stellar object (QSO) - sub-millimetre galaxy (SMG) system BRI1202-0725 at z = 4.7, which have recently been detected in the [CII]158um line by the Atacama Large Millimetre/Sub-millimetre Array (ALMA). We detect Ly-alpha emission from both sources and so confirm that these Ly-alpha emitter candidates are physically associated with the BRI1202- 0725 system. We also report the lack of detection of any high ionisation emission lines (N V, Si IV, C IV and He II) and find that these systems are likely not photoionised by the quasar, leaving in situ star formation as the main powering source of these LAEs. We also find that both LAEs have Ly-alpha emission much broader (1300 km/s) than the [CII] emission and broader than most LAEs. In addition, both LAEs have roughly symmetric Ly-alpha profiles implying that both systems are within the HII sphere produced by the quasar. This is the first time that the proximity zone of a quasar is probed by exploiting nearby Ly-alpha emitters. We discuss the observational properties of these galaxies in the context of recent galaxy formation models.

preprint2014arXiv

Constraints on the star-formation rate of z~3 LBGs with measured metallicity in the CANDELS GOODS-South field

We analyse 14 LBGs at z~2.8-3.8 constituting the only sample where both a spectroscopic measurement of their metallicity and deep IR observations (CANDELS+HUGS survey) are available. Fixing the metallicity of population synthesis models to the observed values, we determine best-fit physical parameters under different assumptions about the star-formation history and also consider the effect of nebular emission. For comparison we determine the UV slope of the objects, and use it to estimate their SFR_UV99 by correcting the UV luminosity following Meurer et al. (1999). A comparison between SFR obtained through SED-fitting (SFR_fit) and the SFR_UV99 shows that the latter are underestimated by a factor 2-10, regardless of the assumed SFH. Other SFR indicators (radio, far-IR, X-ray, recombination lines) coherently indicate SFRs a factor of 2-4 larger than SFR_UV99 and in closer agreement with SFR_fit. This discrepancy is due to the solar metallicity implied by the usual beta-A1600 conversion factor. We propose a refined relation, appropriate for sub-solar metallicity LBGs: A1600 = 5.32+1.99beta. This relation reconciles the dust-corrected UV with the SED-fitting and the other SFR indicators. We show that the fact that z~3 galaxies have sub-solar metallicity implies an upward revision by a factor of ~1.5-2 of the global SFRD, depending on the assumptions about the age of the stellar populations. We find very young best-fit ages (10-500 Myrs) for all our objects. From a careful examination of the uncertainties in the fit and the amplitude of the Balmer break we conclude that there is little evidence of the presence of old stellar population in at least half of the LBGs in our sample, suggesting that these objects are probably caught during a huge star-formation burst, rather than being the result of a smooth evolution.

preprint2014arXiv

Dynamics and Metallicity of Far-IR Selected Galaxies

We present near-infrared integral field spectroscopy of ten Herschel selected galaxies at z~1.5. From detailed mapping of the H$α$ and [NII] emission lines we trace the dynamics, star formation rates, metallicities and also investigate gas fractions for these galaxies. For a few galaxies the distribution of star formation as traced by H$α$ only traces a small fraction of the stellar disc, which could be tracing recent minor merging events. The rest-frame ultraviolet (UV) continuum emission often has a distribution completely different from H$α$, which warns about the use of UV-SED based star formation tracers in these systems. Our analysis of galaxy dynamics shows that minor dynamical disruptions (e.g. minor merging) are generally not enough to cause a deviation from the established 'Main Sequence' relation. Most galaxies are found to follow the fundamental metallicity relation (FMR), although with large scatter. One galaxy, (a small satellite galaxy of a massive companion) is found to deviate strongly from the FMR. This deviation is in nice agreement with the correlation recently discovered in local galaxies between gas metallicity and environment, which has been ascribed to enriched inter-galactic medium (IGM) in dense environments, and therefore suggests that here the IGM was already significantly enriched by z~1.5.

preprint2014arXiv

Galaxy-wide outflows in z$\sim$1.5 luminous obscured QSOs revealed through NIR slit-resolved spectroscopy

Aims. The co-evolution of galaxies and super massive black holes (SMBHs) requires that some sort of feedback mechanism is operating during the active galactic nuclei (AGN) phases. AGN driven winds are the most likely candidates for such feedback mechanism, but direct observational evidence of their existence and of their effects on the host galaxies are still scarce and their physical origin is still hotly debated. Methods. X-shooter observations of a sample of X-ray selected, obscured quasars at z$\sim$1.5, selected on the basis of their observed red colors and X-ray-to-optical flux ratio, have shown the presence of outflowing ionized gas identified by broad [OIII] emission lines in 6 out of 8 objects, confirming the efficiency of the selection criteria. Here we present slit-resolved spectroscopy for the two brightest sources, XID2028 and XID5321, to study the complex emission and absorption line kinematics. Results. We detect outflow extended out to $\sim$ 10 kpc from the central black hole (BH), both as blueshifted and redshifted emission. Interestingly, we also detect kpc scale outflows in the [OII] emission lines and in the neutral gas component, traced by the sodium D and magnesium absorption lines, confirming that a substantial amount of the outflowing mass is in the form of neutral gas. Conclusions. The measured gas velocities and the outflow kinetic powers, inferred under reasonable assumptions on the geometry and physical properties of these two systems, favor an AGN origin for the observed winds.

preprint2014arXiv

Gas reservoir of a hyper-luminous QSO at z=2.6

Understanding the relationship between the formation and evolution of galaxies and their central super massive black holes (SMBH) is one of the main topics in extragalactic astrophysics. Links and feedback may reciprocally affect both black hole and galaxy growth. Observations of the CO line at redshifts of 2-4 are crucial to investigate the gas mass, star formation activity and accretion onto SMBHs, as well as the effect of AGN feedback. Potential correlations between AGN and host galaxy properties can be highlighted by observing extreme objects. Despite their luminosity, hyper-luminous QSOs at z=2-4 are still little studied at mm wavelengths. We targeted CO(3-2) in ULAS J1539+0557, an hyper-luminos QSO (Lbol> 10^48 erg/s) at z=2.658, selected through its unusual red colors in the UKIDSS Large Area Survey (ULAS). We find a molecular gas mass of 4.1+-0.8 10^10 Msun, and a gas fraction of 0.4-0.1, depending mostly on the assumed source inclination. We also find a robust lower limit to the star-formation rate (SFR=250-1600 Msun/yr) and star-formation efficiency (SFE=25-350 Lsun/(K km s-1 pc2) by comparing the observed optical-near-infrared spectral energy distribution with AGN and galaxy templates. The black hole gas consumption timescale, M(H_2)/dM(accretion)/dt, is ~160 Myr, similar or higher than the gas consumption timescale. The gas content and the star formation efficiency are similar to those of other high-luminosity, highly obscured QSOs, and at the lower end of the star-formation efficiency of unobscured QSOs, in line with predictions from AGN-galaxy co-evolutionary scenarios. Further measurements of the (sub)-mm continuum in this and similar sources are mandatory to obtain a robust observational picture of the AGN evolutionary sequence.

preprint2014arXiv

Half-Megasecond Chandra Spectral Imaging of the Hot Circumgalactic Nebula around Quasar Mrk 231

A deep 400-ksec ACIS-S observation of the nearest quasar known, Mrk 231, is combined with archival 120-ksec data obtained with the same instrument and setup to carry out the first ever spatially resolved spectral analysis of a hot X-ray emitting circumgalactic nebula around a quasar. The 65 x 50 kpc X-ray nebula shares no resemblance with the tidal debris seen at optical wavelengths. One notable exception is the small tidal arc 3.5 kpc south of the nucleus where excess soft X-ray continuum emission and Si XIII 1.8 keV line emission are detected, consistent with star formation and its associated alpha-element enhancement, respectively. An X-ray shadow is also detected at the location of the 15-kpc northern tidal tail. The hard X-ray continuum emission within 6 kpc of the center is consistent with being due entirely to the bright central AGN. The soft X-ray spectrum of the outer (>6 kpc) portion of the nebula is best described as the sum of two thermal components with T~3 and ~8 million K and spatially uniform super-solar alpha element abundances, relative to iron. This result implies enhanced star formation activity over ~10^8 yrs accompanied with redistribution of the metals on large scale. The low-temperature thermal component is not present within 6 kpc of the nucleus, suggesting extra heating in this region from the circumnuclear starburst, the central quasar, or the wide-angle quasar-driven outflow identified from optical IFU spectroscopy on a scale of >3 kpc. Significant azimuthal variations in the soft X-ray intensity are detected in the inner region where the outflow is present. The soft X-ray emission is weaker in the western quadrant, coincident with a deficit of Halpha and some of the largest columns of neutral gas outflowing from the nucleus. Shocks created by the interaction of the wind with the ambient ISM may heat the gas to high temperatures at this location. (abridged)

preprint2014arXiv

Karl G. Jansky Very Large Array observations of cold dust and molecular gas in starbursting quasar host galaxies at z~4.5

We present Karl G. Jansky Very Large Array (VLA) observations of 44 GHz continuum and CO J=2-1 line emission in BR1202-0725 at z=4.7 (a starburst galaxy and quasar pair) and BRI1335-0417 at z=4.4 (also hosting a quasar). With the full 8 GHz bandwidth capabilities of the upgraded VLA, we study the (rest-frame) 250 GHz thermal dust continuum emission for the first time along with the cold molecular gas traced by the Low-J CO line emission. The measured CO J=2-1 line luminosities of BR1202-0725 are L'(CO) = (8.7+/-0.8)x10^10 K km/s pc^2 and L'(CO) = (6.0+/-0.5)x10^10 K km/s pc^2 for the submm galaxy (SMG) and quasar, which are equal to previous measurements of the CO J=5-4 line luminosities implying thermalized line emission and we estimate a combined cold molecular gas mass of ~9x10^10 Msun. In BRI1335-0417 we measure L'(CO) = (7.3+/-0.6)x10^10 K km/s pc^2. We detect continuum emission in the SMG BR1202-0725 North (S(44GHz) = 51+/-6 microJy), while the quasar is detected with S(44GHz) = 24+/-6 microJy and in BRI1335-0417 we measure S(44GHz) = 40+/-7 microJy. Combining our continuum observations with previous data at (rest-frame) far-infrared and cm-wavelengths, we fit three component models in order to estimate the star-formation rates. This spectral energy distribution fitting suggests that the dominant contribution to the observed 44~GHz continuum is thermal dust emission, while either thermal free-free or synchrotron emission contributes less than 30%.

preprint2014arXiv

Metallicity evolution, metallicity gradients and gas fractions at z~3.4

We used near-infrared integral field spectroscopic observations from the AMAZE and LSD programs to constrain the metallicity in a sample of 40 star forming galaxies at 3<z<5 (most of which at z~3.4). We measure metallicities by exploiting strong emission line diagnostics. We found that a significant fraction of star-forming galaxies at z~3.4 deviate from the Fundamental Metallicity Relation (FMR), with a metallicity up to a factor of ten lower than expected according to the FMR. This deviation does not correlate with the dynamical properties of the galaxy or with the presence of interactions. To investigate the origin of the metallicity deviations in more detail, we also infer information on the gas content, by inverting the Schmidt-Kennicutt relation. In agreement with recent CO observational data, we found that, in contrast with the steeply rising trend at 0<z<2, the gas fraction in massive galaxies remains constant, with indication of a marginal decline, at 2<z<3.5. When combined with the metallicity information, we infer that to explain both the low metallicity and gas content in z~3.4 galaxies, both prominent outflows and massive pristine gas inflows are needed. In ten galaxies we can also spatially resolve the metallicity distribution. We found that the metallicity generally anticorrelates with the distribution of star formation and with the gas surface density. We discuss these findings in terms of pristine gas inflows towards the center, and outflows of metal rich gas from the center toward the external regions. (Abridged)

preprint2014arXiv

Spatially resolved kinematics, galactic wind, and quenching of star formation in the luminous infrared galaxy IRAS F11506-3851

We present a multi-wavelength integral field spectroscopic study of the low-z LIRG IRAS F11506-3851, on the basis of VIMOS and SINFONI (ESO-VLT) observations. The morphology and the 2D kinematics of the gaseous (neutral and ionized) and stellar components have been mapped using the NaD doublet, the H$α$ line, and the near-IR CO(2-0) and CO(3-1) bands. The kinematics of the ionized gas and the stars are dominated by rotation, with large observed velocity amplitudes and centrally peaked velocity dispersion maps. The stars lag behind the warm gas and represent a dynamically hotter system, as indicated by the observed dynamical ratios. Thanks to these IFS data we have disentangled the contribution of the stars and the ISM to the NaD feature, finding that it is dominated by the absorption of neutral gas clouds in the ISM. The neutral gas 2D kinematics shows a complex structure dominated by two components. On the one hand, the thick slowly rotating disk lags significantly compared to the ionized gas and the stars, with an irregular and off-center velocity dispersion map. On the other hand, a kpc-scale neutral gas outflow is observed along the semi-minor axis of the galaxy, as revealed by large blueshifted velocities (30-154 km/s). We derive an outflowing mass rate in neutral gas of about 48 $\dot{M_{\rm w}}$/yr. Although this implies a global mass loading factor of 1.4, the 2D distribution of the ongoing SF suggests a much larger value of mass loading factor associated with the inner regions (R$<$200 pc), where the current SF represents only 3 percent of the total. All together these results strongly suggest that we are witnessing (nuclear) quenching due to SF feedback in IRAS F11506-3851. However, the relatively large mass of molecular gas detected in the nuclear region via the H2 1-0 S(1) line suggests that further episodes of SF may take place again.

preprint2014arXiv

The CAMbridge Emission Line Surveyor (CAMELS)

The CAMbridge Emission Line Surveyor (CAMELS) is a pathfinder program to demonstrate on-chip spectrometry at millimetre wavelengths. CAMELS will observe at frequencies from 103-114.7 GHz, providing 512 channels with a spectral resolution of R = 3000. In this paper we describe the science goals of CAMELS, the current system level design for the instrument and the work we are doing on the detailed designs of the individual components. In addition, we will discuss our efforts to understand the impact that the design and calibration of the filter bank on astronomical performance. The shape of the filter channels, the degree of overlap and the nature of the noise all effect how well the parameters of a spectral line can be recovered. We have developed a new and rigorous method for analysing performance, based on the concept of Fisher information. This can in be turn coupled to a detailed model of the science case, allowing design trade-offs to be properly investigated.

preprint2014arXiv

The hidden quasar nucleus of a WISE-selected, hyperluminous, dust-obscured galaxy at z ~ 2.3

We present the first X-ray spectrum of a Hot dust-obscured galaxy (DOG), namely W1835+4355 at z ~ 2.3. Hot DOGs represent a very rare population of hyperluminous (>= 10^47 erg/s), dust-enshrouded objects at z > 2 recently discovered in the WISE All Sky Survey. The 40 ks XMM-Newton spectrum reveals a continuum as flat (Gamma ~ 0.8) as typically seen in heavily obscured AGN. This, along with the presence of strong Fe Kalpha emission, clearly suggests a reflection-dominated spectrum due to Compton-thick absorption. In this scenario, the observed luminosity of L(2-10 keV) ~ 2 x 10^44 erg/s is a fraction (<10%) of the intrinsic one, which is estimated to be >~ 5 x 10^45 erg/s by using several proxies. The Herschel data allow us to constrain the SED up to the sub-mm band, providing a reliable estimate of the quasar contribution (~ 75%) to the IR luminosity as well as the amount of star formation (~ 2100 Msun/yr). Our results thus provide additional pieces of evidence that associate Hot DOGs with an exceptionally dusty phase during which luminous quasars and massive galaxies co-evolve and a very efficient and powerful AGN-driven feedback mechanism is predicted by models.

preprint2014arXiv

The host galaxies of X-ray selected Active Galactic Nuclei to z=2.5: Structure, star-formation and their relationships from CANDELS and Herschel/PACS

We study the relationship between the structure and star-formation rate (SFR) of X-ray selected low and moderate luminosity active galactic nuclei (AGNs) in the two Chandra Deep Fields, using Hubble Space Telescope imaging from the Cosmic Assembly Near Infrared Extragalactic Legacy Survey (CANDELS) and deep far-infrared maps from the PEP+GOODS-Herschel survey. We derive detailed distributions of structural parameters and FIR luminosities from carefully constructed control samples of galaxies, which we then compare to those of the AGNs. At z~1, AGNs show slightly diskier light profiles than massive inactive (non-AGN) galaxies, as well as modestly higher levels of gross galaxy disturbance (as measured by visual signatures of interactions and clumpy structure). In contrast, at z~2, AGNs show similar levels of galaxy disturbance as inactive galaxies, but display a red central light enhancement, which may arise due to a more pronounced bulge in AGN hosts or due to extinguished nuclear light. We undertake a number of tests of these alternatives, but our results do not strongly favour one interpretation over the other. The mean SFR and its distribution among AGNs and inactive galaxies are similar at z>1.5. At z<1, however, clear and significant enhancements are seen in the SFRs of AGNs with bulge-dominated light profiles. These trends suggest an evolution in the relation between nuclear activity and host properties with redshift, towards a minor role for mergers and interactions at z>1.5.

preprint2014arXiv

The M_BH - M_star relation for X-ray obscured, red QSOs at 1.2< z <2.6

We present near-infrared spectra, obtained with SINFONI and XShooter observations at ESO VLT, of nine dusty, red QSOs at 1.2<z<2.6. The sources are hard X-ray detected, characterized by cold absorption (N_H>10^{21} - 10^{22} cm^{-2}) and show a broad Ha component in the NIR spectra. We complement this sample with twelve additional sources taken from the literature with similar properties resulting in a total sample of 21 X-ray obscured, intermediate type (1.8-1.9), dusty reddened QSOs. From the broad Ha line we have computed the BH masses through the virial formula and derived Eddington ratios. Moreover, from optical/IR multi-component SED fitting we have derived the stellar mass of their host galaxies and their SFRs. We find that most of the sources in our sample are hosted in starburst and main sequence star-forming galaxies with Eddington ratios lambda>0.1. We find a strong trend with the BH mass i.e. less massive objects are scattered below and above the local relation while the most massive ones are mainly located above it. We also studied the evolution of these sources on the M_BH-M_star plane compared to a sample of optically blue type--1 QSOs and we find that obscured red QSOs show a ratio of M_BH to M_star that increases with redshift which is consistent with or slightly lower than what has been found for blue QSOs. These sources may represent the blow-out phase at the end of the rapid BH growth and immediately preceding the classical blue QSOs typically sampled in optical surveys. They in fact show evidence of outflows in the ionized gas component, but their BH has already fully formed.

preprint2014arXiv

Varying [CII]/[NII] line ratios in the interacting system BR1202-0725 at z=4.7

We study the properties of the interstellar medium in the interacting system BR1202-0725 at z=4.7 via its [NII] and [CII] fine-structure line emission. This system consists of a QSO, a sub-mm galaxy (SMG), and two Ly-alpha emitters (LAEs). Such a diversity in galaxy properties makes BR1202-0725 a unique laboratory of star formation and galaxy evolution at high redshift. We present ionized nitrogen ([NII] 205 micron) observations of this system, obtained with the IRAM Plateau de Bure Interferometer. We find no [NII] emission at the quasar location, but tentative [NII] line detections associated with the SMG and one of the LAEs. Together with available ionized carbon ([CII] 158 micron) ALMA observations of this system, we find the following: The [CII]/[NII] luminosity ratio is >5.5 for the QSO and the SMG, but it is as low as ~2 in the LAE, suggesting that, in this source, most of the [CII] emission is associated with the ionized medium (HII regions) rather than the neutral one (PDRs). This study demonstrates the importance of combined studies of multiple fine--structure lines in order to pin down the physical properties of the interstellar medium in distant galaxies.

preprint2014arXiv

Very extended cold gas, star formation and outflows in the halo of a bright QSO at z>6

Past observations of QSO host galaxies at z >6 have found cold gas and star formation on compact scales of a few kiloparsecs. We present new high sensitivity IRAM PdBI follow-up observations of the [CII] 158micron emission line and FIR continuum in the host galaxy of SDSS J1148+5152, a luminous QSO at redshift 6.4189. We find that a large fraction of the gas traced by [CII] is at high velocities, up to ~1400 km/s relative to the systemic velocity, confirming the presence of a major quasar-driven outflow indicated by previous observations. The outflow has a complex morphology and reaches a maximum projected radius of ~30 kpc. The extreme spatial extent of the outflow allows us, for the first time in an external galaxy, to estimate mass-loss rate, kinetic power and momentum rate of the outflow as a function of the projected distance from the nucleus and the dynamical time-scale. These trends reveal multiple outflow events during the past 100 Myr, although the bulk of the mass, energy and momentum appear to have been released more recently, within the past ~20 Myr. Surprisingly, we discover that also the quiescent gas at systemic velocity is extremely extended. More specifically, we find that, while 30% of the [CII] within v\in(-200, 200) km/s traces a compact component that is not resolved by our observations, 70% of the [CII] emission in this velocity range is extended, with a projected FWHM size of 17.4+-1.4 kpc. We detect FIR continuum emission associated with both the compact and the extended [CII] components, although the extended FIR emission has a FWHM of 11+-3 kpc, thus smaller than the extended [CII] source. Overall, our results indicate that the cold gas traced by [CII] is distributed up to r~30 kpc. A large fraction of extended [CII] is likely associated with star formation on large scales, but the [CII] source extends well beyond the FIR continuum.

preprint2014arXiv

X-shooter reveals powerful outflows in z~1.5 X-ray selected obscured quasi stellar objects

We present X-shooter at the Vewry Large Telescope observations of a sample of 10 luminous, X-ray obscured QSOs at z$\sim1.5$ from the XMM-COSMOS survey, expected to be caught in the transitioning phase from starburst to active galactic nucleus (AGN)-dominated systems. The main selection criterion is X-ray detection at bright fluxes (L$>=10^{44}$ erg s$^{-1}$) coupled to red optical-to-NIR-to-MIR colors. Thanks to its large wavelength coverage, X-shooter allowed us to determine accurate redshifts from the presence of multiple emission lines for five out of six targets for which we had only a photometric redshift estimate, with an 80 percent success rate, significantly larger than what is observed in similar programs of spectroscopic follow-up of red QSOs. We report the detection of broad and shifted components in the [OIII]$λλ$5007,4959 complexes for six out of eight sources with these lines observable in regions free from strong atmospheric absorptions. The full width half-maximum (FWHM) associated with the broad components are in the range FWHM$\sim900-1600$ km s$^{-1}$, larger than the average value observed in SDSS Type 2 AGN samples at similar observed [OIII] luminosity, but comparable to those observed for QSO/ultraluminous infrared galaxies systems for which the presence of kpc scale outflows have been revealed through integral field unit spectroscopy. Although the total outflow energetics (inferred under reasonable assumptions) may be consistent with winds accelerated by stellar processes, we favour an AGN origin for the outflows given the high outflow velocities oberved (v$>1000$ km s$^{-1}$) and the presence of strong winds also in objects undetected in the far infrared.

preprint2013arXiv

A Community Science Case for E-ELT HIRES

Building on the experience of the high-resolution community with the suite of VLT high-resolution spectrographs, which has been tremendously successful, we outline here the (science) case for a high-fidelity, high-resolution spectrograph with wide wavelength coverage at the E-ELT. Flagship science drivers include: the study of exo-planetary atmospheres with the prospect of the detection of signatures of life on rocky planets; the chemical composition of planetary debris on the surface of white dwarfs; the spectroscopic study of protoplanetary and proto-stellar disks; the extension of Galactic archaeology to the Local Group and beyond; spectroscopic studies of the evolution of galaxies with samples that, unlike now, are no longer restricted to strongly star forming and/or very massive galaxies; the unraveling of the complex roles of stellar and AGN feedback; the study of the chemical signatures imprinted by population III stars on the IGM during the epoch of reionization; the exciting possibility of paradigm-changing contributions to fundamental physics. The requirements of these science cases can be met by a stable instrument with a spectral resolution of R~100,000 and broad, simultaneous spectral coverage extending from 370nm to 2500nm. Most science cases do not require spatially resolved information, and can be pursued in seeing-limited mode, although some of them would benefit by the E-ELT diffraction limited resolution. Some multiplexing would also be beneficial for some of the science cases. (Abridged)

preprint2013arXiv

A fundamental relation between the metallicity, gas content, and stellar mass of local galaxies

Recent results have suggested that the well known mass-metallicity relation has a strong dependence on the star formation rate, to the extent that a three dimensional `fundamental metallicity relation' exists which links the three parameters with minimal scatter. In this work, we use a sample of 4253 local galaxies observed in atomic hydrogen from the ALFALFA survey to demonstrate, for the first time, that a similar fundamental relation (the HI-FMR) also exists between stellar mass, gas-phase metallicity, and HI mass. This latter relation is likely more fundamental, driving the relation between metallicity, SFR and mass. At intermediate masses, the behaviour of the gas fundamental metallicity relation is very similar to that expressed via the star formation rate. However, we find that the dependence of metallicity on HI content persists to the highest stellar masses, in contrast to the `saturation' of metallicity with SFR. It is interesting to note that the dispersion of the relation is very low at intermediate stellar masses (9< log(M*/Msun) <11), suggesting that in this range galaxies evolve smoothy, in an equilibrium between gas inflow, outflow and star formation. At high and low stellar masses, the scatter of the relation is significantly higher, suggesting that merging events and/or stochastic accretion and star formation may drive galaxies outside the relation. We also assemble a sample of galaxies observed in CO. However, due to a small sample size, strong selection bias, and the influence of a metallicity-dependent CO/H2 conversion factor, the data are insufficient to test any influence of molecular gas on metallicity.

preprint2013arXiv

A GIANO-TNG high resolution IR spectrum of the airglow emission

A flux-calibrated high resolution spectrum of the airglow emission is a practical lambda-calibration reference for astronomical spectral observations. It is also useful for constraining the molecular parameters of the OH molecule and the physical conditions in the upper mesosphere. methods: We use the data collected during the first technical commissioning of the GIANO spectrograph at the Telescopio Nazionale Galileo (TNG). The high resolution (R~50,000) spectrum simultaneously covers the 0.95-2.4 micron wavelength range. Relative flux calibration is achieved by the simultaneous observation of spectrophotometric standard star. results: We derive a list of improved positions and intensities of OH infrared lines. The list includes Lambda-split doublets many of which are spectrally resolved. Compared to previous works, the new results correct errors in the wavelengths of the Q-branch transitions. The relative fluxes of OH lines from different vibrational bands show remarkable deviations from theoretical predictions: the Deltav=3,4 lines are a factor of 2 and 4 brighter than expected. We also find evidence of a significant fraction (1-4%) of OH molecules with ``non-thermal'' population of high-J levels. Finally we list wavelengths and fluxes of 153 lines not attributable to OH. Most of these can be associated to O2, while 37 lines in the H band are not identified. The O2 and unidentified lines in the H band account for ~5% of the total airglow flux in this band.

preprint2013arXiv

A molecular line scan in the Hubble Deep Field North

We present a molecular line scan in the Hubble Deep Field North (HDF-N) that covers the entire 3mm window (79-115 GHz) using the IRAM Plateau de Bure Interferometer. Our CO redshift coverage spans z<0.45, 1<z<1.9 and all z>2. We reach a CO detection limit that is deep enough to detect essentially all z>1 CO lines reported in the literature so far. We have developed and applied different line searching algorithms, resulting in the discovery of 17 line candidates. We estimate that the rate of false positive line detections is ~2/17. We identify optical/NIR counterparts from the deep ancillary database of the HDF-N for seven of these candidates and investigate their available SEDs. Two secure CO detections in our scan are identified with star-forming galaxies at z=1.784 and at z=2.047. These galaxies have colors consistent with the `BzK' color selection and they show relatively bright CO emission compared with galaxies of similar dust continuum luminosity. We also detect two spectral lines in the submillimeter galaxy HDF850.1 at z=5.183. We consider an additional 9 line candidates as high quality. Our observations also provide a deep 3mm continuum map (1-sigma noise level = 8.6 $μ$Jy/beam). Via a stacking approach, we find that optical/MIR bright galaxies contribute only to <50% of the SFR density at 1<z<3, unless high dust temperatures are invoked. The present study represents a first, fundamental step towards an unbiased census of molecular gas in `normal' galaxies at high-z, a crucial goal of extragalactic astronomy in the ALMA era.

preprint2013arXiv

A Molecular Line Scan in the Hubble Deep Field North: Constraints on the CO Luminosity Function and the Cosmic H2 Density

We present direct constraints on the CO luminosity function at high redshift and the resulting cosmic evolution of the molecular gas density, $ρ_{\rm H2}$(z), based on a blind molecular line scan in the Hubble Deep Field North (HDF-N) using the IRAM Plateau de Bure Interferometer. Our line scan of the entire 3mm window (79-115 GHz) covers a cosmic volume of ~7000 Mpc$^3$, and redshift ranges z<0.45, 1.01<z<1.89 and z>2. We use the rich multiwavelength and spectroscopic database of the HDF-N to derive some of the best constraints on CO luminosities in high redshift galaxies to date. We combine the blind CO detections in our molecular line scan (presented in a companion paper) with stacked CO limits from galaxies with available spectroscopic redshifts (slit or mask spectroscopy from Keck and grism spectroscopy from HST) to give first blind constraints on high-z CO luminosity functions and the cosmic evolution of the H2 mass density $ρ_{\rm H2}$(z) out to redshifts z~3. A comparison to empirical predictions of $ρ_{\rm H2}$(z) shows that the securely detected sources in our molecular line scan already provide significant contributions to the predicted $ρ_{\rm H2}$(z) in the redshift bins <z>~1.5 and <z>~2.7. Accounting for galaxies with CO luminosities that are not probed by our observations results in cosmic molecular gas densities $ρ_{\rm H2}$(z) that are higher than current predictions. We note however that the current uncertainties (in particular the luminosity limits, number of detections, as well as cosmic volume probed) are significant, a situation that is about to change with the emerging ALMA observatory.

preprint2013arXiv

ALMA observations of cool dust in a low-metallicity starburst, SBS0335-052

We present Atacama Large Millimeter/submillimeter Array (ALMA) Cycle 0 Band 7 observations of an extremely metal-poor dwarf starburst galaxy in the Local Universe, SBS0335-052 (12+log(O/H)~7.2). With these observations, dust is detected at 870micron (ALMA Band 7), but 87% of the flux in this band is due to free-free emission from the starburst. We have compiled a spectral energy distribution (SED) of SBS0335-052 that spans almost 6 orders of magnitude in wavelength and fit it with a spherical dust shell heated by a single-age stellar population; the best-fit model gives a dust mass of (3.8+/-0.6)x10^4 Msun. We have also constructed a SED including Herschel archival data for IZw18, another low-metallicity dwarf starburst (12+log(O/H)=7.17), and fit it with a similar model to obtain a dust mass of (3.4+/-1.0)x10^2 Msun. Compared with their atomic gas mass, the dust mass of SBS0335-052 far exceeds the prediction of a linear trend of dust-to-gas mass ratio with metallicity, while IZw18 falls far below. We use gas scaling relations to assess a putative missing gas component in both galaxies and find that the missing, possibly molecular, gas in SBS0335-052 is a factor of 6 times higher than the value inferred from the observed HI column density; in IZw18 the missing component is 4 times smaller. Ultimately, despite their similarly low metallicity, the differences in gas and dust column densities in SBS0335-052 and IZw18 suggest that metal abundance does not uniquely define star-formation processes. At some level, self-shielding and the survival of molecules may depend just as much on gas and dust column density as on metallicity. The effects of low metallicity may at least be partially compensated for by large column densities in the interstellar medium.

preprint2013arXiv

Chronos: A NIR spectroscopic galaxy survey. From the formation of galaxies to the peak of activity

Chronos is our response to ESA's call for white papers to define the science for the future L2, L3 missions. Chronos targets the formation and evolution of galaxies, by collecting the deepest NIR spectroscopic data, from the formation of the first galaxies at z~10 to the peak of formation activity at z~1-3. The strong emission from the atmospheric background makes this type of survey impossible from a ground-based observatory. The spectra of galaxies represent the equivalent of a DNA fingerprint, containing information about the past history of star formation and chemical enrichment. The proposed survey will allow us to dissect the formation process of galaxies including the timescales of quenching triggered by star formation or AGN activity, the effect of environment, the role of infall/outflow processes, or the connection between the galaxies and their underlying dark matter haloes. To provide these data, the mission requires a 2.5m space telescope optimised for a campaign of very deep NIR spectroscopy. A combination of a high multiplex and very long integration times will result in the deepest, largest, high-quality spectroscopic dataset of galaxies from z=1 to 12, spanning the history of the Universe, from 400 million to 6 billion years after the big bang, i.e. covering the most active half of cosmic history.

preprint2013arXiv

Fast Molecular Outflows in Luminous Galaxy Mergers: Evidence for Quasar Feedback from Herschel

We report the results from a systematic search for molecular (OH-119 um) outflows with Herschel-PACS in a sample of 43 nearby (z < 0.3) galaxy mergers, mostly ultraluminous infrared galaxies (ULIRGs) and QSOs. We find that the character of the OH feature (strength of the absorption relative to the emission) correlates with that of the 9.7-um silicate feature, a measure of obscuration in ULIRGs. Unambiguous evidence for molecular outflows, based on the detection of OH absorption profiles with median velocities more blueshifted than -50 km/sec, is seen in 26 (70%) of the 37 OH-detected targets, suggesting a wide-angle (~145 degrees) outflow geometry. Conversely, unambiguous evidence for molecular inflows, based on the detection of OH absorption profiles with median velocities more redshifted than +50 km/sec, is seen in only 4 objects, suggesting a planar or filamentary geometry for the inflowing gas. Terminal outflow velocities of ~-1000 km/sec are measured in several objects, but median outflow velocities are typically ~-200 km s^{-1}. While the outflow velocities show no statistically significant dependence on the star formation rate, they are distinctly more blueshifted among systems with large AGN fractions and luminosities [log (L_AGN / L_sun) > 11.8 +/- 0.3]. The quasars in these systems play a dominant role in driving the molecular outflows. In contrast, the most AGN dominated systems, where OH is seen purely in emission, show relatively modest OH line widths, despite their large AGN luminosities, perhaps indicating that molecular outflows subside once the quasar has cleared a path through the obscuring material.

preprint2013arXiv

GIANO-TNG spectroscopy of red supergiants in the young star cluster RSGC2

The inner disk of the Galaxy has a number of young star clusters dominated by red supergiants that are heavily obscured by dust extinction and observable only at infrared wavelengths. These clusters are important tracers of the recent star formation and chemical enrichment history in the inner Galaxy. During the technical commissioning and as a first science verification of the GIANO spectrograph at the Telescopio Nazionale Galileo, we secured high-resolution (R~50,000) near-infrared spectra of three red supergiants in the young Scutum cluster RSGC2. Taking advantage of the full YJHK spectral coverage of GIANO in a single exposure, we were able to identify several tens of atomic and molecular lines suitable for chemical abundance determinations. By means of spectral synthesis and line equivalent width measurements, we obtained abundances of Fe and other iron-peak elements such as V, Cr, Ni, of alpha (O, Mg, Si, Ca and Ti) and other light elements (C, N, Na, Al, K, Sc), and of some s-process elements (Y, Sr). We found iron abundances between half and one third solar and solar-scaled [X/Fe] abundance patterns of iron-peak, alpha and most of the light elements, consistent with a thin-disk chemistry. We found a depletion of [C/Fe] and enhancement of [N/Fe], consistent with CN burning, and low 12C/13C abundance ratios (between 9 and 11), requiring extra-mixing processes in the stellar interiors during the post-main sequence evolution. Finally, we found a slight [Sr/Fe] enhancement and a slight [Y/Fe] depletion (by a factor of <=2), with respect to solar.

preprint2013arXiv

Gravitational Lens Models Based on Submillimeter Array Imaging of Herschel-selected Strongly Lensed Sub-millimeter Galaxies at z>1.5

Strong gravitational lenses are now being routinely discovered in wide-field surveys at (sub)millimeter wavelengths. We present Submillimeter Array (SMA) high-spatial resolution imaging and Gemini-South and Multiple Mirror Telescope optical spectroscopy of strong lens candidates discovered in the two widest extragalactic surveys conducted by the Herschel Space Observatory: the Herschel-Astrophysical Terahertz Large Area Survey (H-ATLAS) and the Herschel Multi-tiered Extragalactic Survey (HerMES). From a sample of 30 Herschel sources with S_500>100 mJy, 21 are strongly lensed (multiply imaged), 4 are moderately lensed (singly imaged), and the remainder require additional data to determine their lensing status. We apply a visibility-plane lens modeling technique to the SMA data to recover information about the masses of the lenses as well as the intrinsic (i.e., unlensed) sizes (r_half) and far-infrared luminosities (L_FIR) of the lensed submillimeter galaxies (SMGs). The sample of lenses comprises primarily isolated massive galaxies, but includes some groups and clusters as well. Several of the lenses are located at z_lens>0.7, a redshift regime that is inaccessible to lens searches based on Sloan Digital Sky Survey spectroscopy. The lensed SMGs are amplified by factors that are significantly below statistical model predictions given the 500um flux densities of our sample. We speculate that this may reflect a deficiency in our understanding of the intrinsic sizes and luminosities of the brightest SMGs. The lensed SMGs span nearly one decade in L_FIR (median L_FIR=7.9x10^12 L_sun) and two decades in FIR luminosity surface density (median Sigma_FIR=6.0x10^11 L_sun kpc^-2). The strong lenses in this sample and others identified via (sub-)mm surveys will provide a wealth of information regarding the astrophysics of galaxy formation and evolution over a wide range in redshift.

preprint2013arXiv

High resolution mapping of CO(1-0) in NGC6240

We present high spatial resolution CO(1-0) mapping of the luminous infrared galaxy NGC6240 obtained with the IRAM - Plateau de Bure Interferometer. This source is a well-known early-stage merging system hosting two AGN. We find a broad CO(1-0) line profile, with maximum velocity 800 km/s and FWZI=1400 km/s, and displaying several kinematic components, revealing the complexity of the gas dynamics in this system. We detect a blueshifted CO emission with velocity in the range -200 and -500 km/s, peaked around the southern AGN, at the same position where the H2 outflow is located, and with a mass loss rate of ~500 Msun/yr. We interpret this blueshifted component as a outflow, originating from the southern nucleus. The spatial and spectral match strongly suggests that the CO outflow is connected to the H2 superwind located around the southern nucleus, and to the large scale CO outflow, with similar velocities, extended on scales of 10 kpc. The large mass loading factor (dM/dt /SFR ~ 10) of the molecular gas suggests that the outflow is likely driven by both SNa winds and the radiation of the southern AGN. We discovered a nuclear, redshifted CO emission peaking in the mid point of the two nuclei, as it is the case for the CO emission at the systemic velocity. The large velocity dispersion, which reaches its maximum (~500 Km/s) in the mid point between the two nuclei, suggests that the gas might be highly turbulent in this region, although the presence of an unresolved rotation component cannot be ruled out.

preprint2013arXiv

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

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

preprint2013arXiv

Massive Molecular Outflows and Evidence for AGN Feedback from CO Observations

We study the properties of massive, galactic-scale outflows of molecular gas and investigate their impact on galaxy evolution. We present new IRAM PdBI CO(1-0) observations of local ULIRGs and QSO hosts: clear signature of massive and energetic molecular outflows, extending on kpc scales, is found in the CO(1-0) kinematics of four out of seven sources, with measured outflow rates of several 100 Msun/yr. We combine these new observations with data from the literature, and explore the nature and origin of massive molecular outflows within an extended sample of 19 local galaxies. We find that starburst-dominated galaxies have an outflow rate comparable to their SFR, or even higher by a factor of ~ 2-4, implying that starbursts can indeed be effective in removing cold gas from galaxies. Nevertheless, our results suggest that the presence of an AGN can boost the outflow rate by a large factor, which is found to increase with the L_AGN/L_bol ratio. The gas depletion time-scales due to molecular outflows are anti-correlated with the presence and luminosity of an AGN in these galaxies, and range from a few hundred million years in starburst galaxies, down to just a few million years in galaxies hosting powerful AGNs. In quasar hosts the depletion time-scales due to the outflow are much shorter than the depletion time-scales due to star formation. We estimate the outflow kinetic power and find that, for galaxies hosting powerful AGNs, it corresponds to about 5% of the AGN luminosity, as expected by models of AGN feedback. Moreover, we find that momentum rates of about 20 L_AGN/c are common among the AGN-dominated sources in our sample. For "pure" starburst galaxies our data tentatively support models in which outflows are mostly momentum-driven by the radiation pressure from young stars onto dusty clouds.

preprint2013arXiv

Metals in the IGM approaching the re-ionization epoch: results from X-shooter at the VLT

We present the results of observations taken with the X-shooter spectrograph devoted to the study of quasars at z ~ 6. This paper focuses on the properties of metals at high redshift traced, in particular, by the C IV doublet absorption systems. Six objects were observed with resolutions ~27 and 34 km/s in the visual, and 37.5 and 53.5 km/s in the near infrared. We detected 102 C IV lines in the range: 4.35 < z < 6.2 of which 27 are above z ~ 5. Thanks to the characteristics of resolution and spectral coverage of X-shooter, we could also detect 25 Si IV doublets associated with the C IV at z>5. The column density distribution function of the C IV line sample is observed to evolve in redshift for z>~ 5.3, with respect to the normalization defined by low redshift (1.5 < z <4) C IV lines. This behaviour is reflected in the redshift evolution of the C IV cosmic mass density, Omega_CIV, of lines with column density in the range 13.4 < log N(C IV) < 15, which is consistent with a drop of a factor of ~2 for z>~ 5.3. Considering only the stronger C IV lines (13.8 < log N(CIV) < 15), Omega_CIV gently rises by a factor of ~10 between z ~ 6.2 and z~ 1.5 with a possible flattening toward z ~ 0. The increase is well fitted by a power law: Omega_CIV = (2+-1)x10^{-8} [(1+z)/4]^{-3.1+-0.1}. An insight into the properties of the C IV absorbers and their evolution with redshift is obtained by comparing the observed column densities of associated C IV, Si IV and C II absorptions with the output of a set of CLOUDY photo-ionization models. As already claimed by cosmological simulations, we find that C IV is a good tracer of the metallicity in the low-density intergalacitc medium gas at z ~ 5-6 while at z ~ 3 it arises in gas with over density delta ~100.

preprint2013arXiv

Nuclear Activity is more prevalent in Star-Forming Galaxies

We explore the question of whether low and moderate luminosity Active Galactic Nuclei (AGNs) are preferentially found in galaxies that are undergoing a transition from active star formation to quiescence. This notion has been suggested by studies of the UV-to-optical colors of AGN hosts, which find them to be common among galaxies in the so-called "Green Valley", a region of galaxy color space believed to be composed mostly of galaxies undergoing star-formation quenching. Combining the deepest current X-ray and Herschel. PACS far-infrared (FIR) observations of the two Chandra Deep Fields (CDFs) with redshifts, stellar masses and rest-frame photometry derived from the extensive and uniform multi-wavelength data in these fields, we compare the rest-frame U-V color distributions and SFR distributions of AGNs and carefully constructed samples of inactive control galaxies. The UV-to-optical colors of AGNs are consistent with equally massive inactive galaxies at redshifts out to z~2, but we show that such colors are poor tracers of star formation. While the FIR distributions of both star-forming AGNs and star-forming inactive galaxies are statistically similar, we show that AGNs are preferentially found in star-forming host galaxies, or, in other words, AGNs are less likely to be found in weakly star-forming or quenched galaxies. We postulate that, among X-ray selected AGNs of low and moderate accretion luminosities, the supply of cold gas primarily determines the accretion rate distribution of the nuclear black holes.

preprint2013arXiv

Strongly star-forming rotating disks in a complex merging system at z = 4,7 as revealed by ALMA

We performed a kinematical analysis of the [CII] line emission of the BR 1202-0725 system at z~4,7 using ALMA observations. The most prominent sources of this system are a quasar and a submillimeter galaxy, separated by a projected distance of about 24 kpc and characterized by very high SFR, higher than 1000 Msun/yr. However, the ALMA observations reveal that these galaxies apparently have undisturbed rotating disks, which is at variance with the commonly accepted scenario in which strong star formation activity is induced by a major merger. We also detected faint components which, after spectral deblending, were spatially resolved from the main QSO and SMG emissions. The relative velocities and positions of these components are compatible with orbital motions within the gravitational potentials generated by the QSO host galaxy and the SMG, suggesting that they are smaller galaxies in interaction or gas clouds in accretion flows of tidal streams. We did not find any clear spectral evidence for outflows caused by AGN or stellar feedback. This suggests that the high star formation rates might be induced by interactions or minor mergers with these companions, which do not affect the large-scale kinematics of the disks, however. Our kinematical analysis also indicates that the QSO and the SMG have similar Mdyn, mostly in the form of molecular gas, and that the QSO host galaxy and the SMG are seen close to face-on with slightly different disk inclinations: the QSO host galaxy is seen almost face-on (i~15), while the SMG is seen at higher inclinations (i~25). Finally, the ratio between the black hole mass of the QSO, obtained from XShooter spectroscopy, and the Mdyn of the host galaxy is similar to value found in very massive local galaxies, suggesting that the evolution of black hole galaxy relations is probably better studied with dynamical than with stellar host galaxy masses.

preprint2013arXiv

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

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

preprint2013arXiv

The evolution of the dust and gas content in galaxies

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

preprint2013arXiv

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

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

preprint2013arXiv

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

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

preprint2013arXiv

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

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

preprint2013arXiv

The Hot and Energetic Universe: Astrophysics of feedback in local AGN

Understanding the astrophysics of feedback in active galactic nuclei (AGN) is key to understanding the growth and co-evolution of supermassive black holes and galaxies. AGN-driven winds/outflows are potentially the most effective way of transporting energy and momentum from the nuclear scales to the host galaxy, quenching star formation by sweeping away the gas reservoir. Key questions in this field are: 1) how do accretion disks around black holes launch winds/outflows, and how much energy do these carry? 2) How are the energy and metals accelerated in winds/outflows transferred and deposited into the circumgalactic medium? X-ray observations are a unique way to address these questions because they probe the phase of the outflows which carries most of the kinetic energy. We show how a high throughput, high spectral resolution instrument like the X-ray Integral Field Unit (X-IFU) on Athena+ will allow us to address these questions by determining the physical parameters (ionization state, density, temperature, abundances, velocities, geometry, etc.) of the outflows on a dynamical time-scale, in a broad sample of nearby bright AGN. The X-IFU will also allow direct spectral imaging of the impact of these winds on the host galaxy for local AGN, forming a template for understanding AGN at higher redshifts where wind shocks cannot be resolved.

preprint2013arXiv

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

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

preprint2013arXiv

The Mean Star-Forming Properties of QSO Host Galaxies

Quasi-stellar objects (QSOs) occur in galaxies in which supermassive black holes (SMBHs) are growing substantially through rapid accretion of gas. Many popular models of the co-evolutionary growth of galaxies and SMBHs predict that QSOs are also sites of substantial recent star formation, mediated by important processes, such as major mergers, which rapidly transform the nature of galaxies. A detailed study of the star-forming properties of QSOs is a critical test of such models. We present a far-infrared Herschel/PACS study of the mean star formation rate (SFR) of a sample of spectroscopically observed QSOs to z~2 from the COSMOS extragalactic survey. This is the largest sample to date of moderately luminous AGNs studied using uniform, deep far-infrared photometry. We study trends of the mean SFR with redshift, black hole mass, nuclear bolometric luminosity and specific accretion rate (Eddington ratio). To minimize systematics, we have undertaken a uniform determination of SMBH properties, as well as an analysis of important selection effects within spectroscopic QSO samples that influence the interpretation of SFR trends. We find that the mean SFRs of these QSOs are consistent with those of normal massive star-forming galaxies with a fixed scaling between SMBH and galaxy mass at all redshifts. No strong enhancement in SFR is found even among the most rapidly accreting systems, at odds with several co-evolutionary models. Finally, we consider the qualitative effects on mean SFR trends from different assumptions about the star-forming properties of QSO hosts and redshift evolution of the SMBH-galaxy relationship. While limited currently by uncertainties, valuable constraints on AGN-galaxy co-evolution can emerge from our approach.

preprint2012arXiv

[CII] line emission in massive star-forming galaxies at z=4.7

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of the [CII] 157.7micron fine structure line and thermal dust continuum emission from a pair of gas-rich galaxies at z=4.7, BR1202-0725. This system consists of a luminous quasar host galaxy and a bright submm galaxy (SMG), while a fainter star-forming galaxy is also spatially coincident within a 4" (25 kpc) region. All three galaxies are detected in the submm continuum, indicating FIR luminosities in excess of 10^13 Lsun for the two most luminous objects. The SMG and the quasar host galaxy are both detected in [CII] line emission with luminosities, L([CII]) = (10.0 +/- 1.5)x10^9 Lsun and L([CII]) = (6.5+/-1.0)x10^9 Lsun, respectively. We estimate a luminosity ratio, L([CII])/L(FIR) = (8.3+/-1.2)x10^-4 for the starburst SMG to the North, and L([CII])/L(FIR) = (2.5+/-0.4)x10^-4 for the quasar host galaxy, in agreement with previous high-redshift studies that suggest lower [CII]-to-FIR luminosity ratios in quasars than in starburst galaxies. The third fainter object with a flux density, S(340GHz) = 1.9+/-0.3 mJy, is coincident with a Ly-Alpha emitter and is detected in HST ACS F775W and F814W images but has no clear counterpart in the H-band. Even if this third companion does not lie at a similar redshift to BR1202-0725, the quasar and the SMG represent an overdensity of massive, infrared luminous star-forming galaxies within 1.3 Gyr of the Big Bang.

preprint2012arXiv

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

We study a sample of 61 submillimetre galaxies (SMGs) selected from ground-based surveys, with known spectroscopic redshifts and observed with Herschel as part of the PACS Evolutionary Probe (PEP) and the Herschel Multi-tiered Extragalactic Survey (HerMES) key programmes. We use the broad far-infrared wavelength coverage (100-600um) provided by the combination of PACS and SPIRE observations. Using a power-law temperature distribution model to derive infrared luminosities and dust temperatures, we measure a dust emissivity spectral index for SMGs of beta=2.0+/-0.2. Our results unveil the diversity of the SMG population. Some SMGs exhibit extreme infrared luminosities of ~10^13 Lsun and relatively warm dust components, while others are fainter (~10^12 Lsun) and are biased towards cold dust temperatures. The extreme infrared luminosities of some SMGs (LIR>10^12.7 Lsun, 26/61 systems) imply SFRs of >500Msun yr^-1. Such high SFRs are difficult to reconcile with a secular mode of star formation, and may instead correspond to a merger-driven stage in the evolution of these galaxies. Another observational argument in favour of this scenario is the presence of dust temperatures warmer than that of SMGs of lower luminosities (~40K as opposed to ~25K), consistent with observations of local ULIRGs triggered by major mergers and with results from hydrodynamic simulations of major mergers combined with radiative transfer calculations. Luminous SMGs are also offset from normal star-forming galaxies in the stellar mass-SFR plane, suggesting that they are undergoing starburst events with short duty cycles, compatible with the major merger scenario. On the other hand, a significant fraction of the low infrared luminosity SMGs have cold dust temperatures, are located close to the main sequence of star formation, and thus might be evolving through a secular mode of star formation. [abridged]

preprint2012arXiv

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

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

preprint2012arXiv

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

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

preprint2012arXiv

Dust formation around AGB and SAGB stars: a trend with metallicity?

We calculate the dust formed around AGB and SAGB stars of metallicity Z=0.008 by following the evolution of models with masses in the range 1M<M<8M throughthe thermal pulses phase, and assuming that dust forms via condensation of molecules within a wind expanding isotropically from the stellar surface. We find that, because of the strong Hot Bottom Burning (HBB) experienced, high mass models produce silicates, whereas lower mass objects are predicted to be surrounded by carbonaceous grains; the transition between the two regimes occurs at a threshold mass of 3.5M. These fndings are consistent with the results presented in a previous investigation, for Z=0.001. However, in the present higher metallicity case, the production of silicates in the more massive stars continues for the whole AGB phase, because the HBB experienced is softer at Z=0.008 than at Z=0.001, thus the oxygen in the envelope, essential for the formation of water molecules, is never consumed completely. The total amount of dust formed for a given mass experiencing HBB increases with metallicity, because of the higher abundance of silicon, and the softer HBB, both factors favouring a higher rate of silicates production. This behaviour is not found in low mass stars,because the carbon enrichment of the stellar surface layers, due to repeated Third Drege Up episodes, is almost independent of the metallicity. Regarding cosmic dust enrichment by intermediate mass stars, we find that the cosmic yield at Z=0.008 is a factor 5 larger than at Z=0.001. In the lower metallicity case carbon dust dominates after about 300 Myr, but at Z=0.008 the dust mass is dominated by silicates at all times,with a prompt enrichment occurring after about 40 Myr, associated with the evolution of stars with masses M =7.5 -8M.

preprint2012arXiv

Dust temperature and CO-to-H2 conversion factor variations in the SFR-M* plane

Deep Herschel imaging and 12CO(2-1) line luminosities from the IRAM PdBI are combined for a sample of 17 galaxies at z>1 from the GOODS-N field. The sample includes galaxies both on and above the main sequence (MS) traced by star-forming galaxies in the SFR-M* plane. The far-infrared data are used to derive dust masses, Mdust. Combined with an empirical prescription for the dependence of the gas-to-dust ratio on metallicity (GDR), the CO luminosities and Mdust values are used to derive for each galaxy the CO-to-H2 conversion factor, alpha_co. Like in the local Universe, the value of alpha_co is a factor of ~5 smaller in starbursts compared to normal star-forming galaxies (SFGs). We also uncover a relation between alpha_co and dust temperature (Tdust; alpha_co decreasing with increasing Tdust) as obtained from modified blackbody fits to the far-infrared data. While the absolute normalization of the alpha_co(Tdust) relation is uncertain, the global trend is robust against possible systematic biases in the determination of Mdust, GDR or metallicity. Although we cannot formally distinguish between a step and a smooth evolution of alpha_co with the dust temperature, we can conclude that in galaxies of near-solar metallicity, a critical value of Tdust=30K can be used to determine whether the appropriate alpha_co is closer to the starburst value (1.0 Msun(K kms pc^2)^-1, if Tdust>30K) or closer to the Galactic value (4.35 Msun (K kms pc^2)^-1, if Tdust<30K). This indicator has the great advantage of being less subjective than visual morphological classifications of mergers/SFGs, which can be difficult at high z because of the clumpy nature of SFGs. In the absence of far-infrared data, the offset of a galaxy from the main sequence (i.e., log[SSFR(galaxy)/SSFR_MS(M*,z)]) can be used to identify galaxies requiring the use of an alpha_co conversion factor lower than the Galactic value.

preprint2012arXiv

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

We compare the average star formation (SF) activity in X-ray selected AGN hosts with mass-matched control inactive galaxies,including star forming and quiescent sources, at 0.5<z<2.5. Recent observations carried out by PACS, the 60-210um Herschel photometric camera, in GOODS-S, GOODS-N and COSMOS allow us to unbiasedly estimate the far-IR luminosity, and hence the SF properties, of the two samples. Accurate AGN host stellar masses are measured by decomposing their total emission into the stellar and nuclear components. We find a higher average SF activity in AGN hosts with respect to non-AGNs. The level of SF enhancement is modest (~0.26dex at ~3sigma) at low X-ray luminosities (Lx<~10^43.5erg/s) and more pronounced (0.56dex at >10sigma) for bright AGNs. However, when comparing to star forming galaxies only, AGN hosts are broadly consistent with the locus of their `main sequence'. We investigate the relative far-IR luminosity distributions of active and inactive galaxies, and find a higher fraction of PACS detected, hence normal and highly star forming systems among AGN hosts. Although different interpretations are possible, we explain our findings as a consequence of a twofold AGN growth path: faint AGNs evolve through secular processes, with instantaneous AGN accretion not tightly linked to the current total SF in the host, while luminous AGNs co-evolve with their hosts through periods of enhanced AGN activity and SF, possibly through major mergers. While an increased SF with respect to non-AGNs of similar mass is expected in the latter, we interpret the modest SF offsets measured in low-Lx AGN hosts as either a) generated by non-synchronous accretion and SF histories in a merger scenario or b) due to possible connections between instantaneous SF and accretion that can be induced by smaller scale (non-major merger) mechanisms. Far-IR luminosity distributions favour the latter scenario.

preprint2012arXiv

Evidence for Low Extinction in Actively Star Forming Galaxies at z$>$6.5

We present a search for the [CII] 158micron fine structure line (a main cooling line of the interstellar medium) and the underlying far-infrared (FIR) continuum in three high-redshift (6.6<z<8.2) star-forming galaxies using the IRAM Plateau de Bure interferometer. We targeted two Lyman-Alpha-selected galaxies (Lyman-Alpha-Emitters, LAEs) with moderate UV-based star formation rates (SFR~20 M_sun/yr; Himiko at z=6.6 and IOK-1 at z=7.0) and a Gamma Ray Burst (GRB) host galaxy (GRB 090423 at z~8.2). Based on our 3 sigma rest-frame FIR continuum limits, previous (rest-frame) UV continuum measurements and spectral energy distribution (SED) fitting, we rule out SED shapes similar to highly obscured galaxies (e.g. Arp220, M82) and less extreme dust-rich nearby spiral galaxies (e.g. M51) for the LAEs. Conservatively assuming a SED shape typical of local spiral galaxies we derive upper limits for the FIR-based star formation rates (SFRs) of ~70 M_sun/yr, ~50 M_sun/yr and ~40 M_sun/yr for Himiko, IOK-1 and GRB 090423, respectively. For the LAEs these limits are only a factor ~3 higher than the published UV-based SFRs (uncorrected for extinction). This indicates that the dust obscuration in the z>6 LAEs studied here is lower by a factor of a few than what has recently been found in some LAEs at lower redshift (2<z<3.5) with similar UV-based SFRs. A low obscuration in our z>6 LAE sample is consistent with recent rest-frame UV studies of z~7 Lyman-Break-Galaxies (LBGs).

preprint2012arXiv

Evidence of strong quasar feedback in the early Universe

Most theoretical models invoke quasar driven outflows to quench star formation in massive galaxies, this feedback mechanism is required to account for the population of old and passive galaxies observed in the local universe. The discovery of massive, old and passive galaxies at z=2, implies that such quasar feedback onto the host galaxy must have been at work very early on, close to the reionization epoch. We have observed the [CII]158um transition in SDSSJ114816.64+525150.3 that, at z=6.4189, is one of the most distant quasars known. We detect broad wings of the line tracing a quasar-driven massive outflow. This is the most distant massive outflow ever detected and is likely tracing the long sought quasar feedback, already at work in the early Universe. The outflow is marginally resolved on scales of about 16 kpc, implying that the outflow can really affect the whole galaxy, as required by quasar feedback models. The inferred outflow rate, dM/dt > 3500 Msun/yr, is the highest ever found. At this rate the outflow can clean the gas in the host galaxy, and therefore quench star formation, in a few million years.

preprint2012arXiv

Ionized nitrogen at high redshift

We present secure [NII[ detections in two mm-bright, strongly lensed objects at high redshift, APM08279+5255 (z=3.911) and MM18423+5938 (z=3.930), using the IRAM Plateau de Bure Interferometer. Due to its ionization energy [NII] is a good tracer of the ionized gas phase in the interstellar medium. The measured fluxes are S([NII])=(4.8+/-0.8) Jy km/s and (7.4+/-0.5) Jy km/s respectively, yielding line luminosities of L([NII]) =(1.8+/-0.3) x 10^9 μ^{-1} Lsun for APM08279+5255 and L([NII]) =(2.8+/-0.2) x 10^9 μ^{-1} Lsun for MM18423+5938. Our high-resolution map of the [NII] and 1 mm continuum emission in MM18423+5938 clearly resolves an Einstein ring in this source, and reveals a velocity gradient in the dynamics of the ionized gas. A comparison of these maps with high-resolution EVLA CO observations enables us to perform the first spatially-resolved study of the dust continuum-to-molecular gas surface brightness (Sigma_{FIR} Sigma_{CO}^N, which can be interpreted as the star formation law) in a high-redshift object. We find a steep relation (N=1.4+/-0.2), consistent with a starbursting environment. We measure a [NII]/FIR luminosity ratio in APM0828+5255 and MM18423+5938 of 9.0 x 10^{-6} and 5.8 x 10^{-6}, respectively. This is in agreement with the decrease of the [NII]/FIR ratio at high FIR luminosities observed in local galaxies.

preprint2012arXiv

MOONS: a Multi-Object Optical and Near-infrared Spectrograph for the VLT

MOONS is a new conceptual design for a Multi-Object Optical and Near-infrared Spectrograph for the Very Large Telescope (VLT), selected by ESO for a Phase A study. The baseline design consists of 1000 fibers deployable over a field of view of 500 square arcmin, the largest patrol field offered by the Nasmyth focus at the VLT. The total wavelength coverage is 0.8um-1.8um and two resolution modes: medium resolution and high resolution. In the medium resolution mode (R=4,000-6,000) the entire wavelength range 0.8um-1.8um is observed simultaneously, while the high resolution mode covers simultaneously three selected spectral regions: one around the CaII triplet (at R=8,000) to measure radial velocities, and two regions at R=20,000 one in the J-band and one in the H-band, for detailed measurements of chemical abundances. The grasp of the 8.2m Very Large Telescope (VLT) combined with the large multiplex and wavelength coverage of MOONS - extending into the near-IR - will provide the observational power necessary to study galaxy formation and evolution over the entire history of the Universe, from our Milky Way, through the redshift desert and up to the epoch of re-ionization at z>8-9. At the same time, the high spectral resolution mode will allow astronomers to study chemical abundances of stars in our Galaxy, in particular in the highly obscured regions of the Bulge, and provide the necessary follow-up of the Gaia mission. Such characteristics and versatility make MOONS the long-awaited workhorse near-IR MOS for the VLT, which will perfectly complement optical spectroscopy performed by FLAMES and VIMOS.

preprint2012arXiv

Near-infrared spectroscopy of a nitrogen-loud quasar SDSS J1707+6443

We present near-infrared spectroscopy of the z=3.2 quasar SDSS J1707+6443, obtained with MOIRCS on the Subaru Telescope. This quasar is classified as a "nitrogen-loud" quasar because of the fairly strong NIII] and NIV] semi-forbidden emission lines from the broad-line region (BLR) observed in its rest-frame UV spectrum. However, our rest-frame optical spectrum from MOIRCS shows strong [OIII] emission from the narrow-line region (NLR) suggesting that, at variance with the BLR, NLR gas is not metal-rich. In order to reconcile these contradictory results, there may be two alternative possibilities; (1) the strong nitrogen lines from the BLR are simply due to a very high relative abundance of nitrogen rather than to a very high BLR metallicity, or (2) the BLR metallicity is not representative of the metallicity of the host galaxy, better traced by the NLR. In either case, the strong broad nitrogen lines in the UV spectrum are not indication of a chemically enriched host galaxy. We estimated the black hole mass and Eddington ratio of this quasar from the velocity width of both CIV and H_beta, that results in log(M_BH/M_sun) = 9.50 and log(L_bol/L_Edd) = -0.34. The relatively high Eddington ratio is consistent with our earlier result that strong nitrogen emission from BLRs is associated with high Eddington ratios. Finally, we detected significant [NeIII] emission from the NLR, implying a quite high gas density of n~10^6 cm^-3 and suggesting a strong coupling between quasar activity and dense interstellar clouds in the host galaxy.

preprint2012arXiv

NGC6240: extended CO structures and their association with shocked gas

We present deep CO observations of NGC6240 performed with the IRAM Plateau de Bure Interferometer (PdBI). NGC6240 is the prototypical example of a major galaxy merger in progress, caught at an early stage, with an extended, strongly-disturbed butterfly-like morphology and the presence of a heavily obscured active nucleus in the core of each progenitor galaxy. The CO line shows a skewed profile with very broad and asymmetric wings detected out to velocities of -600 km/s and +800 km/s with respect to the systemic velocity. The PdBI maps reveal the existence of two prominent structures of blueshifted CO emission. One extends eastward, i.e. approximately perpendicular to the line connecting the galactic nuclei, over scales of ~7 kpc and shows velocities up to -400 km/s. The other extends southwestward out to ~7 kpc from the nuclear region, and has a velocity of -100 km/s with respect to the systemic one. Interestingly, redshifted emission with velocities 400 to 800 km/s is detected around the two nuclei, extending in the east-west direction, and partly overlapping with the eastern blue-shifted structure, although tracing a more compact region of size ~1.7 kpc. The overlap between the southwestern CO blob and the dust lanes seen in HST images, which are interpreted as tidal tails, indicates that the molecular gas is deeply affected by galaxy interactions. The eastern blueshifted CO emission is co-spatial with an Halpha filament that is associated with strong H2 and soft X-ray emission. The analysis of Chandra X-ray data provides strong evidence for shocked gas at the position of the Halpha emission. Its association with outflowing molecular gas supports a scenario where the molecular gas is compressed into a shock wave that propagates eastward from the nuclei. If this is an outflow, the AGN are likely the driving force.

preprint2012arXiv

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

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

preprint2012arXiv

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

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

preprint2012arXiv

Resolved [CII] emission in a lensed quasar at z=4.4

We present one of the first resolved maps of the [CII] 158 micron line, a powerful tracer of the star forming inter-stellar medium, at high redshift. We use the new IRAM PdBI receivers at 350 GHz to map this line in BRI 0952-0115, the host galaxy of a lensed quasar at z=4.4 previously found to be very bright in [CII] emission. The [CII] emission is clearly resolved and our data allow us to resolve two [CII] lensed images associated with the optical quasar images. We find that the star formation, as traced by [CII], is distributed over a region of ~ 1 kpc in size near the quasar nucleus, and we infer a star formation surface density >150 Msun/yr/kpc^2, similar to that observed in local ULIRGs. We also reveal another [CII] component, extended over ~ 12 kpc, and located at ~ 10 kpc from the quasar. We suggest that this component is a companion disk galaxy, in the process of merging with the quasar host, whose rotation field is distorted by the interaction with the quasar host, and where star formation, although intense, is more diffuse. These observations suggest that galaxy merging at high-z can enhance star formation at the same time in the form of more compact regions, in the vicinity of the accreting black hole, and in more extended star forming galaxies.

preprint2012arXiv

Suzaku reveals X-ray continuum piercing the nuclear absorber in Markarian 231

We report the results from a 2011 Suzaku observation of the nearby low-ionization BAL quasar/ULIRG Markarian 231. These data reveal that the X-ray spectrum has undergone a large variation from the 2001 XMM-Newton and BeppoSAX observations. We interpret this finding according to a scenario whereby the X-ray continuum source is obscured by a two-component partial-covering absorber with NH ~10^22 and ~10^24 cm^-2, respectively. The observed spectral change is mostly explained by a progressive appearance of the primary continuum at <10 keV due to the decrease of the covering fraction of the denser absorption component. The properties of the X-ray obscuration in Mrk 231 match well with those of the X-ray shielding gas predicted by the theoretical models for an efficient radiatively-driven acceleration of the BAL wind. In particular, the X-ray absorber might be located at the extreme base of the outflow. We measure a 2-10 keV luminosity of L(2-10) = 3.3 x 10^43 erg s^-1 for the 2011 data set, i.e. an increase of 30% with respect to the 2001 value.

preprint2012arXiv

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

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

preprint2012arXiv

The anatomy of an extreme starburst within 1.3Gyr of the Big Bang revealed by ALMA

We present further analysis of the [CII] 158$μ$m fine structure line and thermal dust continuum emission from the archetype extreme starburst/AGN group of galaxies in the early Universe, BRI 1202-0725 at $z=4.7$, using the Atacama Large Millimeter Array. The group is long noted for having a closely separated (26kpc in projection) FIR-hyperluminous quasar host galaxy and an optically obscured submm galaxy (SMG). A short ALMA test observation reveals a rich laboratory for the study of the myriad processes involved in clustered massive galaxy formation in the early Universe. Strong [CII] emission from the SMG and the quasar have been reported earlier by Wagg et al. (2012) based on these observations. In this letter, we examine in more detail the imaging results from the ALMA observations, including velocity channel images, position-velocity plots, and line moment images. We present detections of [CII] emission from two Ly$α$-selected galaxies in the group, demonstrating the relative ease with which ALMA can detect the [CII] emission from lower star formation rate galaxies at high redshift. Imaging of the [CII] emission shows a clear velocity gradient across the SMG, possibly indicating rotation or a more complex dynamical system on a scale $\sim 10$kpc. There is evidence in the quasar spectrum and images for a possible outflow toward the southwest, as well as more extended emission (a 'bridge'), between the quasar and the SMG, although the latter could simply be emission from Ly$α$-1 blending with that of the quasar at the limited spatial resolution of the current observations. These results provide an unprecedented view of a major merger of gas rich galaxies driving extreme starbursts and AGN accretion during the formation of massive galaxies and supermassive black holes within 1.3 Gyr of the Big Bang.

preprint2012arXiv

The blue UV slopes of z~4 Lyman break galaxies: implications for the corrected star formation rate density

We plan to analyse dust extinction in Lyman Break Galaxies (LBGs) by introducing a new and more reliable approach to their selection and to the characterization of their distribution of UV slopes beta, using deep IR images from HST. We exploit deep WFC3 IR observations of the ERS and HUDF fields over GOODS-South, combined with HST-ACS optical data, to select z~4 LBGs through a new (B-V) vs. (V-H) colour diagram. The UV slope of the selected galaxies is robustly determined by a linear fit over their observed I, Z, Y, J magnitudes, coherently with the original definition of beta. The same fit is used to determine their rest-frame UV magnitudes M1600 through a simple interpolation. We estimate the effect of observational uncertainties with detailed simulations that we also exploit, under a parametric maximum-likelihood approach, to constrain the probability density function of UV slopes PDF(beta) as a function of rest-frame magnitude. We find 142 and 25 robust LBGs in the ERS and HUDF fields respectively, limiting our sample to S/N(H)>10 objects. Our newly defined criteria improve the selection of z~4 LBGs and allow us to exclude red interlopers at lower redshift, especially z~3-3.5 objects. We find that z~4 LBGs are characterized by blue UV slopes, suggesting a low dust extinction: all L<L* galaxies have an average UV slope <beta>\simeq-2.1, while brighter objects only are slightly redder (<beta> -1.9). We find an intrinsic dispersion ~ 0.3 for PDF(beta) at all magnitudes. The SFRD at z~4 corrected according to these estimates turns out to be lower than previously found: log(SFRD)\simeq-1.09 M_sun/yr/Mpc^3. Finally, we discuss how the UV slope of z~4 galaxies changes as a function of the dust-corrected UV magnitude (i.e. SFR) showing that most galaxies with a high SFR (> 80 M_sun/yr) are highly extincted objects. [Abridged]

preprint2012arXiv

The Intense Starburst HDF850.1 in a Galaxy Overdensity at z=5.2 in the Hubble Deep Field

The Hubble Deep Field (HDF) is a region in the sky that provides one of the deepest multi-wavelength views of the distant universe and has led to the detection of thousands of galaxies seen throughout cosmic time. An early map of the HDF at a wavelength of 850 microns that is sensitive to dust emission powered by star formation revealed the brightest source in the field, dubbed HDF850.1. For more than a decade, this source remained elusive and, despite significant efforts, no counterpart at shorter wavelengths, and thus no redshift, size or mass, could be identified. Here we report, using a millimeter wave molecular line scan, an unambiguous redshift determination for HDF850.1 of z=5.183. This places HDF850.1 in a galaxy overdensity at z~5.2 in the HDF, corresponding to a cosmic age of only 1.1 Gyr after the Big Bang. This redshift is significantly higher than earlier estimates and higher than most of the >100 sub-millimeter bright galaxies identified to date. The source has a star formation rate of 850 M_sun/yr and is spatially resolved on scales of 5 kpc, with an implied dynamical mass of ~1.3x10^11 M_sun, a significant fraction of which is present in the form of molecular gas. Despite our accurate redshift and position, a counterpart arising from starlight remains elusive.

preprint2012arXiv

The physics and the structure of the quasar-driven outflow in Mrk 231

Massive AGN-driven outflows are invoked by AGN-galaxy co-evolutionary models to suppress both star formation and black hole accretion. Massive molecular outflows have recently been revealed in some AGN hosts. However, the physical properties and structure of these AGN-driven molecular outflows are still poorly constrained. Here we present new IRAM PdBI observations of Mrk231, the closest quasar known, targeting both the CO(1-0) and CO(2-1) transitions. We detect broad wings in both transitions, tracing a massive molecular outflow with velocities up to 800 km/s. The wings are spatially resolved at high significance level (5-11 sigma), indicating that the molecular outflow extends on the kpc scale. The CO(2-1)/CO(1-0) ratio of the red broad wings is consistent with the ratio observed in the narrow core, while the blue broad wing is less excited than the core. The latter result suggests that quasar driven outflow models invoking shocks (which would predict higher gas excitation) are not appropriate to describe the bulk of the outflow in Mrk231. However, we note that within the central 700 pc the CO(2-1)/CO(1-0) ratio of the red wing is slightly, but significantly, higher than in the line core, suggesting that shocks may play a role in the central region. We also find that the average size of the outflow anticorrelates with the critical density of the transition used as a wind tracer. This indicates that, although diffuse and dense clumps coexist in the outflowing gas, dense outflowing clouds have shorter lifetime and that they evaporate into the diffuse component along the outflow or, more simply, that diffuse clouds are more efficiently accelerated to larger distances by radiation pressure.

preprint2012arXiv

Variations in the Fundamental constants in the QSO Host J1148+5251 at z = 6.4 and the BR1202-0725 System at z = 4.7

We use sensitive observations of three high redshift sources; [CII] fine structure and CO(2-1) rotational transitions for the z=6.4 Quasar host galaxy (QSO) J1148+5251, and [CII] and CO(5-4) transitions from the QSO BR1202-0725 and its sub-millimeter companion (SMG) galaxy at z=4.7. We use these observations to place constraints on the quantity Dz = z(CO) - z(CII) for each source where z(CO) and z(CII) are the observed redshifts of the CO rotational transition and [CII] fine structure transition respectively, using a combination of approaches; 1) By modelling the emission line profiles using `shapelets' to compare both the emission redshifts and the line profiles themselves, in order to make inferences about the intrinsic velocity differences between the molecular and atomic gas, and 2) By performing a marginalisation over all model parameters in order to calculate a non-parametric estimate of Dz. We derive 99% confidence intervals for the marginalised posterior of Dz of (-1.9 pm 1.3) x10^-3, (-3 pm 8) x10^-4 and (-2 pm 4) x10^-3 for J1148+5251, and the BR1202-0725 QSO and SMG respectively. We show the [CII] and CO(2-1) line profiles for J1148+5251 are consistent with each other within the limits of the data, whilst the [CII] and CO(5-4) line profiles from the BR1202-0725 QSO and SMG respectively have 65 and >99.9% probabilities of being inconsistent, with the CO(5-4) lines ~ 30% wider than the [CII] lines. Therefore whilst the observed values of Dz can correspond to variations in the quantity Delta F/F with cosmic time, where F=alpha^2/mu, with alpha the fine structure constant, and mu the proton-to-electron mass ratio, of both (-3.3 pm 2.3) x10^-4 for a look back time of 12.9 Gyr and of (-5 pm 15) x10^-5 for a look back time of 12.4 Gyr we propose that they are the result of the two species of gas being spatially separated as indicated by the inconsistencies in their line profiles.

preprint2011arXiv

A dynamical mass estimator for high z galaxies based on spectroastrometry

Galaxy dynamical masses are important physical quantities to constrain galaxy evolutionary models, especially at high redshifts. However, at z~2 the limited signal to noise ratio and spatial resolution of the data usually do not allow spatially resolved kinematical modeling and very often only virial masses can be estimated from line widths. But even such estimates require a good knowledge of galaxy size, which may be smaller than the spatial resolution. Spectroastrometry is a technique which combines spatial and spectral resolution to probe spatial scales significantly smaller than the spatial resolution of the observations. Here we apply it to the case of high-z galaxies and present a method based on spectroastrometry to estimate dynamical masses of high z galaxies, which overcomes the problem of size determination with poor spatial resolution. We construct and calibrate a "spectroastrometric" virial mass estimator, modifying the "classical" virial mass formula. We apply our method to the [O III] or Hα emission line detected in z~2-3 galaxies from AMAZE, LSD and SINS samples and we compare the spectroastrometric estimator with dynamical mass values resulting from full spatially resolved kinematical modeling. The spectroastrometric estimator is found to be a good approximation of dynamical masses, presenting a linear relation with a residual dispersion of only 0.15 dex. This is a big improvement compared to the "classical" virial mass estimator which has a non linear relation and much larger dispersion (0.47 dex) compared to dynamical masses. By applying our calibrated estimator to 16 galaxies from the AMAZE and LSD samples, we obtain masses in the ~10^7-10^10 M\odot range extending the mass range attainable with dynamical modeling.

preprint2011arXiv

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

The cosmic infrared background (CIB) includes roughly half of the energy radiated by all galaxies at all wavelengths across cosmic time, as observed at the present epoch. The PACS Evolutionary Probe (PEP) survey is exploited here to study the CIB and its redshift differential, at 70, 100 and 160 micron, where the background peaks. Combining PACS observations of the GOODS-S, GOODS-N, Lockman Hole and COSMOS areas, we define number counts spanning over more than two orders of magnitude in flux: from ~1 mJy to few hundreds mJy. Stacking of 24 micron sources and P(D) statistics extend the analysis down to ~0.2 mJy. Taking advantage of the wealth of ancillary data in PEP fields, differential number counts and CIB are studied up to z=5. Based on these counts, we discuss the effects of confusion on PACS blank field observations and provide confusion limits for the three bands considered. The total CIB surface brightness emitted above PEP 3 sigma flux limits is 4.52 +/- 1.18, 8.35 +/- 0.95 and 9.49 +/- 0.59 [nW/m2/sr] at 70, 100, and 160 micron, respectively. These values correspond to 58 +/- 7% and 74 +/- 5% of the COBE/DIRBE CIB direct measurements at 100 and 160 micron. Employing the P(D) analysis, these fractions increase to ~65% and ~89%. More than half of the resolved CIB was emitted at redshift z<=1. The 50%-light redshifts lie at z=0.58, 0.67 and 0.73 at the three PACS wavelengths. The distribution moves towards earlier epochs at longer wavelengths: while the 70 micron CIB is mainly produced by z<=1.0 objects, the contribution of z>1.0 sources reaches 50% at 160 micron. Most of the CIB resolved in the three PACS bands was emitted by galaxies with infrared luminosities in the range 1e11-1e12 L(sun).

preprint2011arXiv

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

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

preprint2011arXiv

Is GRB 050904 at z=6.3 absorbed by dust?

Claim of dust extinction for this GRB has been debated in the past. We suggest that the discrepant results occur primarily because most of previous studies have not simultaneously investigated the X-ray to near-IR spectral energy distribution of this GRB. The difficulty with this burst is that the X-ray afterglow is dominated by strong flares at early times and is poorly monitored at late times. In addition, the Z band photometry, which is the most sensitive to dust extinction, has been found to be affected by strong systematics. In this paper we carefully re-analyze the Swift/XRT afterglow observations of this GRB, using extensive past studies of X-ray flare properties when computing the X-ray afterglow flux level and exploiting the recent reanalysis of the optical (UV rest frame) data of the same GRB. We extract the X-ray to optical/near-IR afterglow SED for the three epochs where the best spectral coverage is available: 0.47, 1.25, and 3.4 days after the trigger. A spectral power-law model has been fitted to the extracted SEDs. We discuss that no spectral breaks or chromatic temporal breaks are expected in the epochs of interest. To fit any UV rest-frame dust absorption, we tested the Small Magellanic Cloud (SMC) extinction curve, the mean extinction curve (MEC) found for a sample of QSO at $z>4$ and its corresponding attenuation curve, as well as a starburst attenuation curve, and the extinction curve consistent with a supernova dust origin (SN-type). The SMC extinction curve and the SN-type one provide good fit to the data at all epochs, with an average amount of dust absorption at $λ_{rest} = 3000 Å$ of $A_{3000} = 0.25\pm 0.07$ mag. These results indicate that the primeval galaxy at $z = 6.3$ hosting this GRB has already enriched its ISM with dust.

preprint2011arXiv

Massive molecular outflows and negative feedback in ULIRGs observed by Herschel-PACS

Mass outflows driven by stars and active galactic nuclei are a key element in many current models of galaxy evolution. They may produce the observed black hole-galaxy mass relation and regulate and quench both star formation in the host galaxy and black hole accretion. However, observational evidence of such feedback processes through outflows of the bulk of the star forming molecular gas is still scarce. Here we report the detection of massive molecular outflows, traced by the hydroxyl molecule (OH), in far-infrared spectra of ULIRGs obtained with Herschel-PACS as part of the SHINING key project. In some of these objects the (terminal) outflow velocities exceed 1000 km/s, and their outflow rates (up to ~1200 M_sol/yr) are several times larger than their star formation rates. We compare the outflow signatures in different types of ULIRGs and in starburst galaxies to address the issue of the energy source (AGN or starburst) of these outflows. We report preliminary evidence that ULIRGs with a higher AGN luminosity (and higher AGN contribution to L_IR) have higher terminal velocities and shorter gas depletion time scales. The outflows in the observed ULIRGs are able to expel the cold gas reservoirs from the centres of these objects within ~1E6-1E8 years.

preprint2011arXiv

Observational evidence of quasar feedback quenching star formation at high redshift

Most galaxy evolutionary models require quasar feedback to regulate star formation in their host galaxies. In particular, at high redshift, models expect that feedback associated with quasar-driven outflows is so efficient that the gas in the host galaxy is largely swept away or heated up, hence suppressing star formation in massive galaxies. We observationally investigate this phenomenon by using VLT-SINFONI integral field spectroscopy of the luminous quasar 2QZJ002830.4-281706 at z=2.4. The spectra sample the optical emission lines redshifted into the near-IR. The [OIII]5007 emission-line kinematics map reveals a massive outflow on scales of several kpc. The detection of narrow Halpha emission reveals star formation in the quasar host galaxy, with SFR=100 Msun/yr. However, the star formation is not distributed uniformly, but is strongly suppressed in the region with the highest outflow velocity and highest velocity dispersion. This result indicates that star formation in this region is strongly quenched by the quasar outflow, which is cleaning the galaxy disk of its molecular gas. This is one of the first direct observational proofs of quasar feedback quenching the star formation at high redshift.

preprint2011arXiv

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

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

preprint2011arXiv

Optical-NIR spectra of quasars close to reionization (z~ 6)

X-shooter, with its characteristics of resolution, spectral coverage and efficiency, provides a unique opportunity to obtain spectra of the highest-redshift quasars (z ~ 6) that will allow us to carry out successful investigations on key cosmological issues, from the details of the re-ionization process, to the evolution of the first galaxies and AGNs. In this paper, we present the spectra of three z ~ 6 quasars: one obtained during the commissioning of X-shooter and two in the context of our ongoing GTO programme. Combining this sample with data in the literature, we update the value of the C IV cosmic mass density in the range 4.5 < z < 5, confirming the constant trend with redshift between 2.5 and 5.

preprint2011arXiv

ORIGIN: Metal Creation and Evolution from the Cosmic Dawn

ORIGIN is a proposal for the M3 mission call of ESA aimed at the study of metal creation from the epoch of cosmic dawn. Using high-spectral resolution in the soft X-ray band, ORIGIN will be able to identify the physical conditions of all abundant elements between C and Ni to red-shifts of z=10, and beyond. The mission will answer questions such as: When were the first metals created? How does the cosmic metal content evolve? Where do most of the metals reside in the Universe? What is the role of metals in structure formation and evolution? To reach out to the early Universe ORIGIN will use Gamma-Ray Bursts (GRBs) to study their local environments in their host galaxies. This requires the capability to slew the satellite in less than a minute to the GRB location. By studying the chemical composition and properties of clusters of galaxies we can extend the range of exploration to lower redshifts (z ~ 0.2). For this task we need a high-resolution spectral imaging instrument with a large field of view. Using the same instrument, we can also study the so far only partially detected baryons in the Warm-Hot Intergalactic Medium (WHIM). The less dense part of the WHIM will be studied using absorption lines at low redshift in the spectra for GRBs.

preprint2011arXiv

PACS Evolutionary Probe (PEP) - A Herschel Key Program

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

preprint2011arXiv

Spectroscopic confirmation of two Lyman break galaxies at redshift beyond 7

We report the spectroscopic confirmation of two Lyman break galaxies at redshift > 7. The galaxies were observed as part of an ultra-deep spectroscopic campaign with FORS2 at the ESO/VLT for the confirmation of z~7 ``z--band dropout'' candidates selected from our VLT/Hawk-I imaging survey. Both galaxies show a prominent emission line at 9735A and 9858A respectively: the lines have fluxes of ~ 1.6-1.2 x 10^(-17) erg/s/cm2 and exhibit a sharp decline on the blue side and a tail on the red side. The asymmetry is quantitatively comparable to the observed asymmetry in z ~ 6 Ly-alpha lines, where absorption by neutral hydrogen in the IGM truncates the blue side of the emission line profile. We carefully evaluate the possibility that the galaxies are instead at lower redshift and we are observing either [OII], [OIII] or H-alpha emission: however from the spectroscopic and the photometric data we conclude that there are no other plausible identifications, except for Ly-alpha at redshift >7, implying that these are two of the most robust redshift determination for galaxies in the reionization epoch. Based on their redshifts and broad--band photometry, we derive limits on the star formation rate and on the ultraviolet spectral slopes of the two galaxies. We argue that these two galaxies alone are unlikely to have ionized the IGM in their surroundings.

preprint2011arXiv

Spectroscopic confirmation of z~7 LBGs: probing the earliest galaxies and the epoch of reionization

We present the final results from our ultra-deep spectroscopic campaign with FORS2 at the ESO/VLT for the confirmation of z~7 "z--band dropout" candidates selected from our VLT/Hawk-I imaging survey over three independent fields. In particular we report on two newly discovered galaxies at redshift ~6.7 in the NTT deep field: both galaxies show a Ly-alpha emission line with rest-frame EWs of the order 15-20 A and luminosities of 2-4 X 10^{42} erg/s. We also present the results of ultra-deep observations of a sample of i-dropout galaxies, from which we set a solid upper limit on the fraction of interlopers. Out of the 20 z-dropouts observed we confirm 5 galaxies at 6.7 < z < 7.1. This is systematically below the expectations drawn on the basis of lower redshift observations: in particular there is a significant lack of objects with intermediate Ly-alpha EWs (between 20 and 55 A). We conclude that the trend for the fraction of Ly-alpha emission in LBGs that is constantly increasing from z~3 to z~6 is most probably reversed from z~6 to z~7. Explaining the observed rapid change in the LAE fraction among the drop-out population with reionization requires a fast evolution of the neutral fraction of hydrogen in the Universe. Assuming that the Universe is completely ionized at z=6 and adopting the semi-analytical models of Dijkstra et al. (2011), we find that our data require a change of the neutral hydrogen fraction of the order Delta chi_{HI} ~ 0.6 in a time Delta z ~ 1, provided that the escape fraction does not increase dramatically over the same redshift interval.

preprint2011arXiv

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

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

preprint2011arXiv

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

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

preprint2011arXiv

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

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

preprint2011arXiv

The transition from carbon dust to silicates production in low-metallicity AGB and SAGB stars

We compute the mass and composition of dust produced by stars with masses in the range 1Msun<M<8 Msun and with a metallicity of Z=0.001 during their AGB and Super AGB phases. Stellar evolution is followed from the pre-main sequence phase using the code ATON which provides, at each timestep, the thermodynamics and the chemical stucture of the wind. We use a simple model to describe the growth of the dust grains under the hypothesis of a time-independent, spherically symmetric stellar wind. We find that the total mass of dust injected by AGB stars in the interstellar medium does not increase monotonically with stellar mass and ranges between a minimum of 10^{-6}Msun for the 1.5Msun stellar model, up to 2x10^{-4} Msun, for the 6Msun case. Dust composition depends on the stellar mass: low-mass stars (M < 3Msun) produce carbon-rich dust, whereas more massive stars, experiencing Hot Bottom Burning, never reach the carbon-star stage, and produce silicates and iron. This is in partial disagreement with previous investigations in the literature, which are based on synthetic AGB models and predict that, when the initial metallicity is Z=0.001, C-rich dust is formed at all stellar masses. The differences are due to the different modelling of turbulent convection in the super-adiabaticity regime. Also in this case the treatment of super-adiabatic convection shows up as the most relevant issue affecting the dust-formation process. We also investigate Super AGB stars with masses 6.5Msun<M<8 Msun that evolve over a ONe core.Due to a favourable combination of mass loss and Hot Bottom Burning, these stars are predicted to be the most efficient silicate-dust producers, releasing [2 - 7]x 10^{-4} Msun masses of dust. We discuss the robustness of these predictions and their relevance for the nature and evolution of dust at early cosmic times.

preprint2010arXiv

"Comets" orbiting a black hole

We use a long (300 ksec), continuous Suzaku X-ray observation of the active nucleus in NGC1365 to investigate the structure of the circumnuclear BLR clouds through their occultation of the X-ray source. The variations of the absorbing column density and of the covering factor indicate that the clouds surrounding the black hole are far from having a spherical geometry (as sometimes assumed), instead they have a strongly elongated and cometary shape, with a dense head (n=10^11 cm^-3) and an expanding, dissolving tail. We infer that the cometary tails must be longer than a few times 10^13 cm and their opening angle must be smaller than a few degrees. We suggest that the cometary shape may be a common feature of BLR clouds in general, but which has been difficult to recognize observationally so far. The cometary shape may originate from shocks and hydrodynamical instabilities generated by the supersonic motion of the BLR clouds into the intracloud medium. As a consequence of the mass loss into their tail, we infer that the BLR clouds probably have a lifetime of only a few months, implying that they must be continuously replenished. We also find a large, puzzling discrepancy (two orders of magnitude) between the mass of the BLR inferred from the properties of the absorbing clouds and the mass of the BLR inferred from photoionization models; we discuss the possible solutions to this discrepancy.

preprint2010arXiv

A fundamental relation between mass, SFR and metallicity in local and high redshift galaxies

We show that the mass-metallicity relation observed in the local universe is due to a more general relation between stellar mass M*, gas-phase metallicity and SFR. Local galaxies define a tight surface in this 3D space, the Fundamental Metallicity Relation (FMR), with a small residual dispersion of ~0.05 dex in metallicity, i.e, ~12%. At low stellar mass, metallicity decreases sharply with increasing SFR, while at high stellar mass, metallicity does not depend on SFR. High redshift galaxies, up to z~2.5 are found to follow the same FMR defined by local SDSS galaxies, with no indication of evolution. The evolution of the mass-metallicity relation observed up to z=2.5 is due to the fact that galaxies with progressively higher SFRs, and therefore lower metallicities, are selected at increasing redshifts, sampling different parts of the same FMR. By introducing the new quantity mu_alpha=log(M*)-alpha log(SFR), with alpha=0.32, we define a projection of the FMR that minimizes the metallicity scatter of local galaxies. The same quantity also cancels out any redshift evolution up to z~2.5, i.e, all galaxies have the same range of values of mu_0.32. At z>2.5, evolution of about 0.6 dex off the FMR is observed, with high-redshift galaxies showing lower metallicities. The existence of the FMR can be explained by the interplay of infall of pristine gas and outflow of enriched material. The former effect is responsible for the dependence of metallicity with SFR and is the dominant effect at high-redshift, while the latter introduces the dependence on stellar mass and dominates at low redshift. The combination of these two effects, together with the Schmidt-Kennicutt law, explains the shape of the FMR and the role of mu_0.32. The small metallicity scatter around the FMR supports the smooth infall scenario of gas accretion in the local universe.

preprint2010arXiv

Dissecting the cosmic infra-red background with Herschel/PEP

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

preprint2010arXiv

Dynamical properties of AMAZE and LSD galaxies from gas kinematics and the Tully-Fisher relation at z~3

We present a SINFONI integral field kinematical study of 33 galaxies at z~3 from the AMAZE and LSD projects which are aimed at studying metallicity and dynamics of high-redshift galaxies. The number of galaxies analyzed in this paper constitutes a significant improvement compared to existing data in the literature and this is the first time that a dynamical analysis is obtained for a relatively large sample of galaxies at z~3. 11 galaxies show ordered rotational motions (~30% of the sample), in these cases we estimate dynamical masses by modeling the gas kinematics with rotating disks and exponential mass distributions. We find dynamical masses in the range 2 \times 10^9 M\odot - 2 \times 10^11 M\odot with a mean value of ~ 2 \times 10^10 M\odot. By comparing observed gas velocity dispersion with that expected from models, we find that most rotating objects are dynamically "hot", with intrinsic velocity dispersions of the order of ~90 km s-1. The median value of the ratio between the maximum disk rotational velocity and the intrinsic velocity dispersion for the rotating objects is 1.6, much lower than observed in local galaxies value (~10) and slightly lower than the z~2 value (2 - 4). Finally we use the maximum rotational velocity from our modeling to build a baryonic Tully-Fisher relation at z~3. Our measurements indicate that z~3 galaxies have lower stellar masses (by a factor of ten on average) compared to local galaxies with the same dynamical mass. However, the large observed scatter suggests that the Tully-Fisher relation is not yet "in place" at these early cosmic ages, possibly due to the young age of galaxies. A smaller dispersion of the Tuly-Fisher relation is obtained by taking into account the velocity dispersion with the use of the S_0.5 indicator, suggesting that turbulent motions might have an important dynamical role.

preprint2010arXiv

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

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

preprint2010arXiv

Far-infrared properties of submillimeter and optically faint radio galaxies

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

preprint2010arXiv

Gas accretion as the origin of chemical abundance gradients in distant galaxies

It has recently been suggested that galaxies in the early Universe can grow through the accretion of cold gas, and that this may have been the main driver of star formation and stellar mass growth. Because the cold gas is essentially primordial, it has a very low abundance of elements heavier than helium (metallicity). As it is funneled to the centre of a galaxy, it will lead the central gas having an overall lower metallicity than gas further from the centre, because the gas further out has been enriched by supernovae and stellar winds, and not diluted by the primordial gas. Here we report chemical abundances across three rotationally-supported star-forming galaxies at z~3, only 2 Gyr after the Big Bang. We find an 'inverse' gradient, with the central, star forming regions having a lower metallicity than less active ones, opposite to what is seen in local galaxies. We conclude that the central gas has been diluted by the accretion of primordial gas, as predicted by 'cold flow' models.

preprint2010arXiv

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

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

preprint2010arXiv

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

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

preprint2010arXiv

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

The star formation rate (SFR) is a key parameter in the study of galaxy evolution. The accuracy of SFR measurements at z~2 has been questioned following a disagreement between observations and theoretical models. The latter predict SFRs at this redshift that are typically a factor 4 or more lower than the measurements. We present star-formation rates based on calorimetric measurements of the far-infrared (FIR) luminosities for massive 1.5<z<2.5, normal star-forming galaxies (SFGs), which do not depend on extinction corrections and/or extrapolations of spectral energy distributions. The measurements are based on observations in GOODS-N with the Photodetector Array Camera & Spectrometer (PACS) onboard Herschel, as part of the PACS Evolutionary Probe (PEP) project, that resolve for the first time individual SFGs at these redshifts at FIR wavelengths. We compare FIR-based SFRs to the more commonly used 24 micron and UV SFRs. We find that SFRs from 24 micron alone are higher by a factor of ~4-7.5 than the true SFRs. This overestimation depends on luminosity: gradually increasing for log L(24um)>12.2 L_sun. The SFGs and AGNs tend to exhibit the same 24 micron excess. The UV SFRs are in closer agreement with the FIR-based SFRs. Using a Calzetti UV extinction correction results in a mean excess of up to 0.3 dex and a scatter of 0.35 dex from the FIR SFRs. The previous UV SFRs are thus confirmed and the mean excess, while narrowing the gap, is insufficient to explain the discrepancy between the observed SFRs and simulation predictions.

preprint2010arXiv

Herschel unveils a puzzling uniformity of distant dusty galaxies

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

preprint2010arXiv

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

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

preprint2010arXiv

SAGACE: the Spectroscopic Active Galaxies And Clusters Explorer

The SAGACE experiment consists of a mm/sub-mm telescope with a 3-m diameter primary mirror, coupled to a cryogenic multi-beam differential spectrometer. SAGACE explores the sky in the 100-760 GHz frequency range, using four diffraction-limited bolometer arrays. The instrument is designed to perform spectroscopic surveys of the Sunyaev-Zeldovich effects of thousands of galaxy clusters, of the spectral energy distribution of active galactic nuclei, and of the [CII] line of a thousand galaxies in the redshift desert. In 2008 a full phase-A study for a national small mission was completed and delivered to the Italian Space Agency (ASI). We have shown that taking advantage of the differential operation of the Fourier Transform Spectrometer, this ambitious instrument can operate from a Molniya orbit, and can be built and operated within the tight budget of a small mission.

preprint2010arXiv

Star formation in AGN hosts in GOODS-N

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

preprint2010arXiv

Studying the WHIM Content of the Galaxy Large-Scale Structures along the Line of Sight to H 2356-309

We make use of a 500ks Chandra HRC-S/LETG spectrum of the blazar H2356-309, combined with a lower S/N spectrum of the same target, to search for the presence of warm-hot absorbing gas associated with two Large-Scale Structures (LSSs) crossed by this sightline at z=0.062 (the Pisces-Cetus Supercluster, PCS) and at z=0.128 ("Farther Sculptor Wall", FSW). No statistically significant (>=3sigma) individual absorption is detected from any of the strong He- or H-like transitions of C, O and Ne at the redshifts of the structures. However we are still able to constrain the physical and geometrical parameters of the associated putative absorbing gas, by performing joint spectral fit of marginal detections and upper limits of the strongest expected lines with our self-consistent hybrid ionization WHIM spectral model. At the redshift of the PCS we identify a warm phase with logT=5.35_-0.13^+0.07 K and log N_H =19.1+/-0.2 cm^-2 possibly coexisting with a hotter and less significant phase with logT=6.9^+0.1_-0.8 K and log N_H=20.1^+0.3_-1.7 cm^-2 (1sigma errors). For the FSW we estimate logT=6.6_-0.2^+0.1 K and log N_H=19.8_-0.8^+0.4 cm^-2. Our constraints allow us to estimate the cumulative number density per unit redshifts of OVII WHIM absorbers. We also estimate the cosmological mass density obtaining Omega_b(WHIM)=(0.021^+0.031_-0.018) (Z/Z_sun)^-1, consistent with the mass density of the intergalactic 'missing baryons' for high metallicities.

preprint2010arXiv

The bright end of the z ~ 7 UV Luminosity Function from a wide and deep HAWK-I survey

(Abridged) We present here the second half of an ESO Large Programme, which exploits the unique combination of area and sensitivity provided in the near-IR by the camera Hawk-I at the VLT. We have obtained - 30 observing hours with Hawk-I in the Y-band of two high galactic latitude fields. We combined the Y-band data with deep J and K Hawk-I observations, and with FORS1/FORS2 U, B, V, R, I, and Z observations to select z-drop galaxies having Z - Y > 1, no optical detection and flat Y - J and Y - K colour terms. We detect 8 high-quality candidates in the magnitude range Y = 25.5 - 26.5 that we add to the z-drop candidates selected in two Hawk-I pointings over the GOODS-South field. We use this full sample of 15 objects found in -161 arcmin^2 of our survey to constrain the average physical properties and the evolution of the number density of z ~ 7 LBGs. A stacking analysis yields a best-fit SED with photometric redshift z= 6.85 +0.20 -0.15 and an E(B-V)=0.05 +0.15 -0.05. We compute a binned estimate of the z ~ 7 LF and explore the effects of photometric scatter and model uncertainties on the statistical constraints. After accounting for the expected incompleteness through MonteCarlo simulations, we strengthen our previous finding that a Schechter luminosity function constant from z=6 to z=7 is ruled out at a >99% confidence level, even including the effects of cosmic variance. For galaxies brighter than M_1500= -19.0, we derive a luminosity density rho_UV = 1.5^{+2.1}{-0.8} x 10^25 erg/s/Hz/Mpc^3, implying a decrease by a factor 3.5 from z=6 to z=6.8. We find that, under standard assumptions, the emission rate of ionizing photons coming from UV bright galaxies is lower by at least a factor of two than the value required for reionization. Finally, we exploit deep Hawk-I J and K band observations to derive an upper limit on the number density of M1500<~ -22.0 LBGs at z-8 (Y-dropouts).

preprint2010arXiv

The chemical evolution of elliptical galaxies with stellar and QSO dust production

The presence of dust strongly affects the way we see galaxies and also the chemical abundances we measure in gas. It is therefore important to study he chemical evolution of galaxies by taking into account dust evolution. We aim at performing a detailed study of abundance ratios of high redshift objects and their dust properties. We focus on Lyman-Break galaxies (LBGs) and Quasar (QSO) hosts and likely progenitors of low- and high-mass present-day elliptical galaxies, respectively. We have adopted a chemical evolution model for elliptical galaxies taking account the dust production from low and intermediate mass stars, supernovae Ia, supernovae II, QSOs and both dust destruction and accretion processes. By means of such a model we have followed the chemical evolution of ellipticals of different baryonic masses. Our model complies with chemical downsizing. We made predictions for the abundance ratios versus metallicity trends for models of differing masses that can be used to constrain the star formation rate, initial mass function and dust mass in observed galaxies. We predict the existence of a high redshift dust mass-stellar mass relationship. We have found a good agreement with the properties of LBGs if we assume that they formed at redshift z=2-4. In particular, a non-negligible amount of dust is needed to explain the observed abundance pattern. We studied the QSO SDSS J114816, one of the most distant QSO ever observed (z=6.4), and we have been able to reproduce the amount of dust measured in this object. The dust is clearly due to the production from supernovae and the most massive AGB stars as well as from the grain growth in the interstellar medium. The QSO dust is likely to dominate only in the very central regions of the galaxies and during the early development of the galactic wind.

preprint2010arXiv

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

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

preprint2010arXiv

The extinction law at high redshift and its implications

We analyze the optical-near infrared spectra of 33 quasars with redshifts 3.9<z<6.4 to investigate the properties of dust extinction at these cosmic epochs. The SMC extinction curve has been shown to reproduce the dust reddening of most quasars at z<2.2; we investigate whether this curve also provides a good description of dust extinction at higher redshifts. We fit the observed spectra with synthetic absorbed quasar templates obtained by varying the intrinsic slope (alpha), the absolute extinction (A3000), and by using a grid of empirical and theoretical extinction curves. We find that seven quasars in our sample are affected by substantial extinction (A3000 >0.8), and characterized by very steep intrinsic slopes (alpha<-2.3). All of the individual quasars require extinction curve deviating from that of the SMC, with a tendency to flatten at lambda<2000 A (in the rest frame of the source). We obtain a mean extinction curve at z>4, both by performing a simultaneous fit of all quasars and by averaging the extinction curves inferred for individual quasars. In the case of broad absorption line quasars the mean extinction curve deviates from the SMC at a confidence level >95%. The difference between extinction curves in quasars at z>4 and in quasars at lower redshift is indicative of either a different dust production mechanism at high redshift, or a different mechanism for processing dust into the ISM. We suggest that the same transitions may also apply to normal, star-forming galaxies at z>4. In particular, the observed change in the average spectral slope of galaxies at z>4 may be partially ascribed to a variation in the extinction curve, rather than a lower dust content at high redshift. In this scenario, the extinction curve inferred at z>4 would imply a cosmic star-formation corrected for dust attenuation a factor of ~2 higher than estimated in the past.

preprint2010arXiv

The far-infrared/radio correlation as probed by Herschel

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

preprint2010arXiv

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

We exploit deep observations of the GOODS-N field taken with PACS, on board of Herschel, as part of the PEP guaranteed time, to study the link between star formation and stellar mass in galaxies to z~2. Starting from a stellar mass-selected sample of ~4500 galaxies with mag[4.5mu]<23 (AB), we identify ~350 objects with a PACS detection at 100 or 160mu and ~1500 with only Spitzer 24 mu counterpart. Stellar masses and total IR luminosities (LIR) are estimated by fitting the SEDs. Consistently with other Herschel results, we find that LIR based only on 24 mu data is overestimated by a median factor ~1.8 at z~2, whereas it is underestimated (with our approach) up to a factor ~1.6 at 0.5<z<1.0. We then exploit this calibration to correct LIR based on the MIPS fluxes. These results clearly show how Herschel is fundamental to constrain LIR, and hence the SFR, of high redshift galaxies. Using the galaxies detected with PACS (and/or MIPS), we investigate the existence and evolution of the relations between the star formation rate (SFR), the specific star formation rate (SSFR=SFR/mass) and the stellar mass. Moreover, in order to avoid selection effects, we also repeat this study through a stacking analysis on the PACS images to fully exploit the far-IR information also for the Herschel and Spitzer undetected subsamples. We find that the SSFR-mass relation steepens with redshift, being almost flat at z<1.0 and reaching a slope of alpha=-0.50^(+0.13)_(-0.16) at z~2, at odds with recent works based on radio-stacking analysis at the same redshift. The mean SSFR of galaxies increases with redshift, by a factor ~15 for massive M>10^(11) (M_sun) galaxies from z=0 to z=2, and seems to flatten at z>1.5 in this mass range. Moreover, the most massive galaxies have the lowest SSFR at any z, implying that they have formed their stars earlier and more rapidly than their low mass counterparts.

preprint2010arXiv

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

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

preprint2010arXiv

The M_BH-M_star relation of obscured AGNs at high redshift

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

preprint2010arXiv

Witnessing the key early phase of quasar evolution: an obscured AGN pair in the interacting galaxy IRAS 20210+1121

We report the discovery of an active galactic nucleus (AGN) pair in the interacting galaxy system IRAS 20210+1121 at z = 0.056. An XMM-Newton observation reveals the presence of an obscured (Nh ~ 5 x 10^{23} cm^-2), Seyfert-like (L_{2-10 keV} = 4.7 x 10^{42} erg/s) nucleus in the northern galaxy, which lacks unambiguous optical AGN signatures. Our spectral analysis also provides strong evidence that the IR-luminous southern galaxy hosts a Type 2 quasar embedded in a bright starburst emission. In particular, the X-ray primary continuum from the nucleus appears totally depressed in the XMM-Newton band as expected in case of a Compton-Thick absorber, and only the emission produced by Compton scattering ('reflection') of the continuum from circumnuclear matter is seen. As such, IRAS 20210+1121 seems to provide an excellent opportunity to witness a key, early phase in the quasar evolution predicted by the theoretical models of quasar activation by galaxy collisions.

preprint2010arXiv

X-ray absorption by Broad Line Region Clouds in Mrk 766

We present a new analysis of a 9-day long XMM-Newton monitoring of the Narrow Line Seyfert 1 galaxy Mrk 766. We show that the strong changes in spectral shape which occurred during this observation can be interpreted as due to Broad Line Region clouds crossing the line of sight to the X-ray source. Within the occultation scenario, the spectral and temporal analysis of the eclipses provides precise estimates of the geometrical structure, location and physical properties of the absorbing clouds. In particular, we show that these clouds have cores with column densities of at least a few 10^23 cm^-2 and velocities in the plane of the sky of the order of thousands km/s. The three different eclipses monitored by XMM-Newton suggest a broad range in cloud velocities (by a factor ~4-5). Moreover, two iron absorption lines clearly associated with each eclipse suggest the presence of highly ionized gas around the obscuring clouds, and an outflow component of the velocity spanning from 3,000 to 15,000 km/s

preprint2009arXiv

A strong excess in the 20-100 keV emission of NGC 1365

We present a new Suzaku observation of the obscured AGN in NGC 1365, revealing an unexpected excess of X-rays above 20 keV of at least a factor ~2 with respect to the extrapolation of the best-fitting 3-10 keV model. Additional Swift-BAT and Integral-IBIS observations show that the 20-100 keV is concentrated within ~1.5 arcmin from the center of the galaxy, and is not significantly variable on time scales from days to years. A comparison of this component with the 3-10 keV emission, which is characterized by a rapidly variable absorption, suggests a complex structure of the circumnuclear medium, consisting of at least two distinct components with rather different physical properties, one of which covering >80% of the source with a column density NH~3-4x10^24 cm^(-2). An alternative explanation is the presence of a double active nucleus in the center of NGC 1365.

preprint2009arXiv

Early metal enrichment in high-redshift quasars

Quasars are powerful systems whose spectrum is rich of metal features that allow us to investigate the chemical evolution of galaxies at very high redshift, even close to the reionization epoch. I review the main observational constraints on the metallicity of quasars host galaxies at high redshift and discuss the implications and issues for models of galaxy evolution in the early universe.

preprint2009arXiv

Evidence for a fast evolution of the UV luminosity function beyond redshift 6 from a deep HAWK-I survey of the GOODS-S field

We perform a deep search for galaxies in the redshift range 6.5<z<7.5, to measure the evolution of the number density of luminous galaxies in this redshift range and derive useful constraints on the evolution of their Luminosity Function. We present here the first results of an ESO Large Program, that exploits the unique combination of area and sensitivity provided in the near-IR by the camera Hawk-I at the VLT. We have obtained two Hawk-I pointings on the GOODS South field for a total of 32 observing hours, covering ~90 arcmin2. The images reach Y=26.7 mags for the two fields. We have used public ACS images in the z band to select z-dropout galaxies with the colour criteria Z-Y>1, Y-J<1.5 and Y-K<2. The other public data in the UBVRIJHK bands are used to reject possible low redshift interlopers. The output has been compared with extensive Monte Carlo simulations to quantify the observational effects of our selection criteria as well as the effects of photometric errors. We detect 7 high quality candidates in the magnitude range Y=25.5-26.7. This interval samples the critical range for M* at z>6 (M_1500 ~- 19.5 to -21.5). After accounting for the expected incompleteness, we rule out at a 99% confidence level a Luminosity Function constant from z=6 to z=7, even including the effects of cosmic variance. For galaxies brighter than M_1500=-19.0 we derive a luminosity density rho_UV=1.5^{+2.0}_{-0.9} 10^25 erg/s/Hz/Mpc3, implying a decrease by a factor 3.5 from z=6 to z~6.8. On the basis of our findings, we make predictions for the surface densities expected in future surveys surveys, based on ULTRA-VISTA, HST-WFC3 or JWST-NIRCam, evaluating the best observational strategy to maximise their impact.

preprint2009arXiv

Exploring the active galactic nucleus and starburst content of local ultraluminous infrared galaxies through 5-8 micron spectroscopy

We present a 5-8 micron analysis of the Spitzer-IRS spectra of 71 ultraluminous infrared galaxies (ULIRGs) with redshift z < 0.15, devoted to the study of the role of active galactic nuclei (AGN) and starbursts (SB) as the power source of the extreme infrared emission. Around 5 micron an AGN is much brighter (by a factor 30) than a starburst of equal bolometric luminosity. This allows us to detect the presence of even faint accretion-driven cores inside ULIRGs: signatures of AGN activity are found in 70 per cent of our sample (50/71 sources). Through a simple analytical model we are also able to obtain a quantitative estimate of the AGN/SB contribution to the overall energy output of each source. Although the main fraction of ULIRG luminosity is confirmed to arise from star formation events, the AGN contribution is non-negligible (23 per cent) and is shown to increase with luminosity. The existence of a rather heterogeneous pattern in the composition and geometrical structure of the dust among ULIRGs is newly supported by the comparison between individual absorption features and continuum extinction.

preprint2009arXiv

LSD: Lyman-break galaxies Stellar populations and Dynamics. I: Mass, metallicity and gas at z~3.1

We present the first results of a project, LSD, aimed at obtaining spatially-resolved, near-infrared spectroscopy of a complete sample of Lyman-Break Galaxies at z~3. Deep observations with adaptive optics resulted in the detection of the main optical lines, such as [OII], Hbeta and [OIII], which are used to study sizes, SFRs, morphologies, gas-phase metallicities, gas fractions and effective yields. Optical, near-IR and Spitzer/IRAC photometry is used to measure stellar mass. We obtain that morphologies are usually complex, with the presence of several peaks of emissions and companions that are not detected in broad-band images. Typical metallicities are 10-50% solar, with a strong evolution of the mass-metallicity relation from lower redshifts. Stellar masses, gas fraction, and evolutionary stages vary significantly among the galaxies, with less massive galaxies showing larger fractions of gas. In contrast with observations in the local universe, effective yields decrease with stellar mass and reach solar values at the low-mass end of the sample. This effect can be reproduced by gas infall with rates of the order of the SFRs. Outflows are present but are not needed to explain the mass-metallicity relation. We conclude that a large fraction of these galaxies are actively creating stars after major episodes of gas infall or merging.

preprint2009arXiv

On the observed distributions of black hole masses and Eddington ratios from radiation pressure corrected virial indicators

The application of the virial theorem to the Broad Line Region of Active Galactic Nuclei allows Black Hole mass estimates for large samples of objects at all redshifts. In a recent paper we showed that ionizing radiation pressure onto BLR clouds affects virial BH mass estimates and we provided empirically calibrated corrections. More recently, a new test of the importance of radiation forces has been proposed: the MBH-sigma relation has been used to estimate MBH for a sample of type-2 AGN and virial relations (with and without radiation pressure) for a sample of type-1 AGN extracted from the same parent population. The observed L/LEdd distribution based on virial BH masses is in good agreement with that based on MBH-sigma only if radiation pressure effects are negligible, otherwise significant discrepancies are observed. In this paper we investigate the effects of intrinsic dispersions associated to the virial relations providing MBH, and we show that they explain the discrepancies between the observed L/LEdd distributions of type-1 and type-2 AGN. These discrepancies in the L/LEdd distributions are present regardless of the general importance of radiation forces, which must be negligible only for a small fraction of sources with large L/LEdd. Average radiation pressure corrections should then be applied in virial MBH estimators until their dependence on observed source physical properties has been fully calibrated. Finally, the comparison between MBH and L/LEdd distributions derived from sigma-based and virial estimators can constrain the variance of BLR physical properties in AGN.

preprint2009arXiv

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

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

preprint2009arXiv

Strong [CII] emission at high redshift

We report the detection of the [CII]157.74um fine-structure line in the lensed galaxy BRI 0952-0115 at z=4.43, using the APEX telescope. This is the first detection of the [CII] line in a source with L_FIR < 10^13 L_sun at high redshift. The line is very strong compared to previous [CII] detections at high-z (a factor of 5-8 higher in flux), partly due to the lensing amplification. The L_[CII]/L_FIR ratio is 10^-2.9, which is higher than observed in local galaxies with similar infrared luminosities. Together with previous observations of [CII] at high redshift, our result suggests that the [CII] emission in high redshift galaxies is enhanced relative to local galaxies of the same infrared luminosity. This finding may result from selection effects of the few current observations of [CII] at high redshift, and in particular the fact that non detections may have not been published (although the few published upper limits are still consistent with the [CII] enhancement scenario). If the trend is confirmed with larger samples, it would indicate that high-z galaxies are characterized by different physical conditions with respect to their local counterparts. Regardless of the physical origin of the trend, this effect would increase the potential of the [CII]158um line to search and characterize high-z sources.

preprint2009arXiv

The evolution of the mass-metallicity relation in galaxies of different morphological types

By means of chemical evolution models for ellipticals, spirals and irregular galaxies, we aim at investigating the physical meaning and the redshift evolution of the mass-metallicity relation as well as how this relation is connected with galaxy morphology. {abridged} We assume that galaxy morphologies do not change with cosmic time. We present a method to account for a spread in the epochs of galaxy formation and to refine the galactic mass grid. (abridged) We compare our predictions to observational results obtained for galaxies between redshifts 0.07 and 3.5. We reproduce the mass-metallicity (MZ) relation mainly by means of an increasing efficiency of star formation with mass in galaxies of all morphological types, without any need to invokegalactic outflows favoring the loss of metals in the less massive galaxies. Our predictions can help constraining the slope and the zero point of the observed local MZ relation, both affected by uncertainties related to the use of different metallicity calibrations. We show how, by considering the MZ, the O/H vs star formation rate (SFR), and the SFR vs galactic mass diagrams at various redshifts, it is possible to constrain the morphology of the galaxies producing these relations. Our results indicate that the galaxies observed at z=3.5 should be mainly proto-ellipticals, whereas at z=2.2 the observed galaxies consist of a morphological mix of proto-spirals and proto-ellipticals. At lower redshifts, the observed MZ relation is well reproduced by considering both spirals and irregulars. (abridged)

preprint2009arXiv

The Mass-Metallicity relation in galaxies of different morphological types

By means of chemical evolution models of different morphological types, we study the mass-metallicity (MZ) relation and its evolution with redshift. Our aim is to understand the role of galaxies of different morphological types in the MZ relation at various redshift. One major result is that at high redshift, the majority of the galaxies falling on the MZ plot are apparently proto-ellipticals. Finally, we show some preliminary results of a study of the MZ relation in a framework of hierarchical galaxy formation.

preprint2008arXiv

AMAZE. I. The evolution of the mass-metallicity relation at z>3

We present initial results of an ESO-VLT large programme (AMAZE) aimed at determining the evolution of the mass-metallicity relation at z>3 by means of deep near-IR spectroscopy. Gas metallicities are measured, for an initial sample of nine star forming galaxies at z~3.5, by means of optical nebular lines redshifted into the near-IR. Stellar masses are accurately determined by using Spitzer-IRAC data, which sample the rest-frame near-IR stellar light in these distant galaxies. When compared with previous surveys, the mass-metallicity relation inferred at z~3.5 shows an evolution much stronger than observed at lower redshifts. The evolution is prominent even in massive galaxies, indicating that z~3 is an epoch of major action in terms of star formation and metal enrichment also for massive systems. There are also indications that the metallicity evolution of low mass galaxies is stronger relative to high mass systems, an effect which can be considered the chemical version of the galaxy downsizing. The mass-metallicity relation observed at z~3.5 is difficult to reconcile with the predictions of some hierarchical evolutionary models. Such discrepancies suggest that at z>3 galaxies are assembled mostly with relatively un-evolved sub-units, i.e. small galaxies with low star formation efficiency. The bulk of the star formation and metallicity evolution probably occurs once small galaxies are already assembled into bigger systems.

preprint2008arXiv

Spectral decomposition of starbursts and AGNs in 5-8 micron Spitzer IRS spectra of local ULIRGs

We present an analysis of the 5-8 micron Spitzer-IRS spectra of a sample of 68 local Ultraluminous Infrared Galaxies (ULIRGs). Our diagnostic technique allows a clear separation of the active galactic nucleus (AGN) and starburst (SB) components in the observed mid-IR emission, and a simple analytic model provides a quantitative estimate of the AGN/starburst contribution to the bolometric luminosity. We show that AGNs are ~30 times brighter at 6 micron than starbursts with the same bolometric luminosity, so that even faint AGNs can be detected. Star formation events are confirmed as the dominant power source for extreme infrared activity, since ~85% of ULIRG luminosity arises from the SB component. Nonetheless an AGN is present in the majority (46/68) of our sources.

preprint2007arXiv

The evolution of the mass-metallicity relation at z~3

We present preliminary results of an ESO-VLT large programme (AMAZE) aimed at determining the evolution of the mass-metallicity relation at z~3 by means of deep near-IR spectroscopy. Gas metallicities and stellar masses are measured for an initial sample of nine star forming galaxies at z~3.3. When compared with previous surveys, the mass-metallicity relation inferred at z~3.3 shows an evolution significantly stronger than observed at lower redshifts. There are also some indications that the metallicity evolution of low mass galaxies is stronger relative to high mass systems, an effect which can be considered as the chemical version of the galaxy downsizing. The mass-metallicity relation observed at z~3.3 is difficult to reconcile with the predictions of some hierarchical evolutionary models. We shortly discuss the possible implications of such discrepancies.

preprint2004arXiv

Near Infrared Spectroscopy of High Redshift Active Galactic Nuclei. I. A Metallicity-Accretion Rate Relationship

We present new near infrared spectroscopic measurements of the H_beta region for a sample of 29 luminous high redshift quasars. We have measured the width of H_beta in those sources, and added archival H_beta width measurements, to create a sample of 92 active galactic nuclei (AGNs) for which H_beta width and rest-frame UV measurements of N V λ1240 and C IV λ1549 emission-lines are available. Our sample spans six orders of magnitude in luminosity and includes 31 radio-loud AGNs. It also includes 10 narrow-line Seyfert 1 galaxies and one broad absorption-line quasar. We find that metallicity, indicated by the N V/C IV line ratio, is primarily correlated with accretion rate, which is a function of luminosity and H_beta line-width. This may imply an intimate relation between starburst, responsible for the metal enrichment of the nuclear gas, and AGN fueling, represented by the accretion rate. The correlation of metallicity with luminosity, or black hole (BH) mass, is weaker in contrast with recent results which were based on measurements of the width of C IV. We argue that using C IV as a proxy to H_beta in estimating M_BH might be problematic and lead to spurious BH mass and accretion rate estimates in individual sources. We discuss the potential implications of our new result in the framework of the starburst-AGN connection and theories of BH growth.

preprint2004arXiv

Near Infrared Spectroscopy of High Redshift Active Galactic Nuclei. II. Disappearing Narrow Line Regions and the Role of Accretion

We present new near infrared spectroscopic measurements for 29 luminous high-z quasars and use the data to discuss the size and other properties of the NLRs in those sources. The high resolution spectra have been used to carefully model the Fe II blends and to provide reliable [O III], Fe II and Hb measurements. We find that about 2/3 of all high luminosity sources show strong [O III] lines while the remaining objects show no or very weak such line. While weak [O III] emitters are also found among lower luminosity AGN, we argue that the implications for very high luminosity objects are different. In particular, we suggest that the averaging of these two populations in other works gave rise to claims of a Baldwin relationship in [O III] which is not confirmed by our data. We also argue that earlier proposed relations of the type R_NLR \propto L_[O III]^{1/2}, where R_NLR is the NLR radius, are theoretically sound yet they must break down for R_NLR exceeding a few kpc. This suggests that the NLR properties in luminous sources are different from those observed in nearby AGN. In particular, we suggest that some sources lost their very large, dynamically unbound NLR while others are in a phase of violent star-forming events that produce a large quantity of high density gas in the central kpc. This gas is ionized and excited by the central radiation source and its spectroscopic properties may be different from those observed in nearby, lower luminosity NLRs. We also discuss the dependence of EW(Hb) and Fe II/Hb on L, M_BH, and accretion rate for a large sample of AGNs. The strongest dependence of the two quantities is on the accretion rate and the Fe II/Hb correlation is probably due to the EW(Hb) dependence on accretion rate. We show the most extreme values measured so far of Fe II/Hb and address its correlation with EW([O III]).

preprint2003arXiv

The infrared supernova rate in starburst galaxies

We report the results of our ongoing search for extincted supernovae (SNe) at near-infrared wavelengths. We have monitored at 2.2 micron a sample of 46 Luminous Infrared Galaxies and detected 4 SNe. The number of detections is still small but sufficient to provide the first estimate of supernova rate at near-infrared wavelengths. We measure a SN rate ofv 7.6+/-3.8 SNu which is an order of magnitude larger than observed in quiescent galaxies. On the other hand, the observed near-infrared rate is still a factor 3-10 smaller than that estimated from the far-infrared luminosity of the galaxies. Among various possibilities, the most likely scenario is that dust extinction is so high (Av>30) to obscure most SNe even in the near-IR. The role of type Ia SNe is also discussed within this context. We derive the type Ia SN rate as a function of the stellar mass of the galaxy and find a sharp increase toward galaxies with higher activity of star formation. This suggests that a significant fraction of type Ia SNe are associated with young stellar populations. Finally, as a by-product, we give the average K-band light curve of core-collapse SNe based on all the existing data, and review the relation between SN rate and far-infrared luminosity.