Search

The number of baryons hosted in the disks of galaxies (stars and interstellar gas and dust) is lower than expected based on the mass of their dark-matter halos and the fraction of baryon-to-total matter in the universe, giving rise to the so called galaxy missing-baryon problem. The presence of cool circum-galactic matter gravitationally bound to the galaxy up to distances of at least ten times the size of the galaxy’s disk, mitigates the problem but is far from being sufficient for its solution. It has instead been suggested that the galaxy missing baryon may hide in a much hotter gaseous phase of the circum-galactic medium, possibly at the halo virial temperature and co-existing with the cool phase. However, current mass estimates of this possible hot phase are only loose lower limits. Here we exploit the best available X-ray spectroscopic data of known cool circum-galactic systems to show that the galaxy missing baryonic mass is indeed found in a hot medium surrounding the disks of galaxies and permeating their entire halos, at temperatures consistent with their virial temperatures. Main tracers of this hot medium are highly ionized light elements that only retain one or two electrons and that are only visible in the X-rays. Our measurements solve the long-standing galaxy missing baryon problems and contribute importantly to the understanding of the continuous cycle of baryons in-and-out of galaxies throughout the life of the universe.

Undisturbed galaxy clusters are characterized by a massive and large elliptical galaxy at their center, i.e. the Brightest Cluster Galaxy (BCG). How these central galaxies form is still debated. According to most models, a typical epoch for their assembly is z~1-2. We have performed a detailed multi-wavelength analysis of the core of XDCPJ0044.0-2033 (XDCP0044), one of the most massive and densest galaxy clusters currently known at redshift z~1.6, whose central galaxy population shows high star formation compared to lower-z clusters and an X-ray AGN located close to its center. SINFONI J-, H- and KMOS YJ-, H- bands spectroscopic data have been analyzed, together with deep archival HST photometric data in F105W, F140W, and F160W bands, Chandra X-ray, radio JVLA data at 1-2 GHz, and ALMA band-6 observations. In the central region of the cluster (~70x70 kpc^2), 2 systems of interacting galaxies have been identified and studied (Complex A and B), with a total of 7 confirmed cluster members. These galaxies show perturbed morphologies and 3 of them show signs of AGN activity. In particular, 2 type-1 AGN with typical broad lines have been found at the center of each complex (both of them X-ray obscured and highly accreting; Eddington ratio ~0.5), while a type-2 AGN has been discovered in Complex A. The AGN at the center of Complex B is also detected in X-ray while the other 2 are spatially related to radio emission. The 3 AGN provide one of the closest AGN triple at z>1 revealed so far with a minimum(maximum) projected distance of ~10(40) kpc. The observation of high star formation, merger signatures and nuclear activity in the core of XDCP0044 suggests that all these processes are key ingredients in shaping the nascent BCG. According to our data, XDCP0044 could form a typical massive galaxy of 10^12 Msun, hosting a Black Hole of 2x10^8-10^9 Msun, in a time scale of the order of ~2.5 Gyrs., 15 pages, 10 figures, 5 tables. Accepted to MNRAS in 2020 August 14

Improving our understanding of the nuclear properties of high-Eddington-ratio (λEdd) active galactic nuclei (AGN) is necessary since at this regime the radiation pressure is expected to affect the structure and efficiency of the accretion disc-corona system. This may cause departures from the typical nuclear properties of low-λEdd AGN, which have been largely studied so far. We present here the X-ray spectral analysis of 14 radio-quiet, λEdd ≳ 1 AGN at 0.4 ≤ z ≤ 0.75, observed with XMM-Newton. Optical/UV data from simultaneous Optical Monitor observations have also been considered. These quasars were selected to have relatively high values of black hole mass (MBH ∼ 108 − 8.5 M⊙) and bolometric luminosity (Lbol ∼ 1046 erg s−1) in order to complement previous studies of high-λEdd AGN at lower MBH and Lbol. We studied the relation between λEdd and other key X-ray spectral parameters, such as the photon index (Γ) of the power-law continuum, the X-ray bolometric correction (kbol, X), and the optical/UV-to-X-ray spectral index (αox). Our analysis reveals that, despite the homogeneous optical and supermassive black hole accretion properties, the X-ray properties of these high-λEdd AGN are quite heterogeneous. We indeed measured values of Γ between 1.3 and 2.5, at odds with the expectations based on previously reported Γ − λEdd relations, for which Γ ≥ 2 would be a ubiquitous hallmark of AGN with λEdd ∼ 1. Interestingly, we found that ∼30% of the sources are X-ray weak, with an X-ray emission about a factor of ∼10 − 80 fainter than that of typical AGN at similar UV luminosities. The X-ray weakness seems to be intrinsic and not due to the presence of absorption along the line of sight to the nucleus. This result may indicate that high-λEdd AGN commonly undergo periods of intrinsic X-ray weakness. Furthermore, results from follow-up monitoring with Swift of one of these X-ray weak sources suggest that these periods can last for several years. © ESO 2022., This work is based on observations obtained with XMM-Newton, an ESA science mission with instruments and contributions directly funded by ESA Member States and NASA. Part of this work is based on archival data, software and online services provided by the Space Science Data Center – ASI. This work has been partially supported by the ASI-INAF program I/004/11/4. ML acknowledges financial support from the PhD programme in Astronomy, Astrophysics and Space Science supported by MIUR (Ministero dell’Istruzione, dell’Università e della Ricerca). EP, LZ, FT, SB, AL and GV acknowledge financial support under ASI/INAF contract 2017-14-H.0. EP, SB, MB and GV acknowledge support from PRIN MIUR project “Black Hole winds and the Baryon Life Cycle of Galaxies: the stone-guest at the galaxy evolution supper”, contract #2017PH3WAT. AdO acknowledges financial support from the State Agency for Research of the Spanish MCIU through the project PID2019-106027GB-C41 and the “Center of Excellence Severo Ochoa” award to the Instituto de Astrofísica de Andalucía (SEV-2017-0709). RM acknowledges the financial support of INAF (Istituto Nazionale di Astrofisica), Osservatorio Astronomico di Roma, ASI (Agenzia Spaziale Italiana) under contract to INAF: ASI 2014-049-R.0 dedicated to SSDC. CR acknowledges support from the Fondecyt Iniciacion grant 11190831.

Accreting supermassive black holes (SMBHs), also known as active galactic nuclei (AGN), are generally surrounded by large amounts of gas and dust. This surrounding material reprocesses the primary X-ray emission produced close to the SMBH and gives rise to several components in the broadband X-ray spectra of AGN, including a power-law possibly associated with Thomson-scattered radiation. In this work, we study the properties of this scattered component for a sample of 386 hard-X-ray-selected, nearby ($z \sim 0.03$) obscured AGN from the 70-month Swift/BAT catalog. We investigate how the fraction of Thomson-scattered radiation correlates with different physical properties of AGN, such as line-of-sight column density, X-ray luminosity, black hole mass, and Eddington ratio. We find a significant negative correlation between the scattering fraction and the column density. Based on a large number of spectral simulations, we exclude the possibility that this anti-correlation is due to degeneracies between the parameters. The negative correlation also persists when considering different ranges of luminosity, black hole mass, and Eddington ratio. We discuss how this correlation might be either due to the angle dependence of the Thomson cross-section or to more obscured sources having a higher covering factor of the torus. We also find a positive correlation between the scattering fraction and the ratio of [OIII] $\lambda$5007 to X-ray luminosity. This result is consistent with previous studies and suggests that the Thomson-scattered component is associated with the narrow-line region., Comment: Accepted for publication in MNRAS, 18 pages, 9 figures, 2 tables, 1 equation

We study the impact of AGN ionised outflows on star formation in high-redshift AGN hosts, by combining NIR IFS observations, mapping the H$\alpha$ emission and [OIII] outflows, with matched-resolution observations of the rest-frame FIR emission. We present high-resolution ALMA Band 7 observations of eight X-ray selected AGN at z~2 from the SUPER sample, targeting the rest-frame ~260 um continuum at ~2 kpc (0.2'') resolution. We detected 6 out of 8 targets with S/N>10 in the ALMA maps, with continuum flux densities F = 0.27-2.58 mJy and FIR half-light radii Re = 0.8-2.1 kpc. The FIR Re of our sample are comparable to other AGN and star-forming galaxies at a similar redshift from the literature. However, we find that the mean FIR size in X-ray AGN (Re = 1.16+/- 0.11 kpc) is slightly smaller than in non-AGN (Re = 1.69+/-0.13 kpc). From SED fitting, we find that the main contribution to the 260 um flux density is dust heated by star formation, with < 4% contribution from AGN-heated dust and < 1% from synchrotron emission. The majority of our sample show different morphologies for the FIR (mostly due to reprocessed stellar emission) and the ionised gas emission (H$\alpha$ and [OIII], mostly due to AGN emission). This could be due to the different locations of dust and ionised gas, the different sources of the emission (stars and AGN), or the effect of dust obscuration. We are unable to identify any residual H$\alpha$ emission, above that dominated by AGN, that could be attributed to star formation. Under the assumption that the FIR emission is a reliable tracer of obscured star formation, we find that the obscured star formation activity in these AGN host galaxies is not clearly affected by the ionised outflows. However, we cannot rule out that star formation suppression is happening on smaller spatial scales than the ones we probe with our observations (< 2 kpc) or on different timescales., Comment: 32 pages, 13 figures. Accepted for publication in A&A

Context. Sources at the brightest end of the quasi-stellar object (QSO) luminosity function, during the peak epoch in the history of star formation and black hole accretion (z ∼ 2-4, often referred to as "Cosmic noon") are privileged sites to study the cycle of feeding & feedback processes in massive galaxies. Aims. We aim to perform the first systematic study of cold gas properties in the most luminous QSOs, by characterising their host-galaxies and environment. These targets exhibit indeed widespread evidence of outflows at nuclear and galactic scales. Methods. We analyse ALMA, NOEMA and JVLA observations of the far-infrared continuum, CO and [CII] emission lines in eight QSOs (bolometric luminosity LBol ≳ 3 × 1047 erg s-1) from the WISE-SDSS selected hyper-luminous (WISSH) QSOs sample at z ∼ 2.4-4.7. Results. 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-millimetre 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), mostly associated with rotating disks with disturbed kinematics. WISSH QSOs typically show lower CO luminosity and higher star formation efficiency than infrared matched, z ∼ 0-3 main-sequence galaxies, implying that, given the observed SFR ∼170-1100 M⊙ yr-1, molecular gas is converted into stars in ≲ 50 Myr. Most targets show extreme dynamical to black-hole mass ratios Mdyn/MBH ∼ 3-10, two orders of magnitude smaller than local relations. The molecular gas fraction in the host-galaxies of WISSH is lower by a factor of ∼10-100 than in star forming galaxies with similar M∗. Conclusions. Our analysis reveals that hyper-luminous QSOs at Cosmic noon undergo an intense growth phase of both the central super-massive black hole and of the host-galaxy. These systems pinpoint the high-density sites where giant galaxies assemble, where we show that mergers play a major role in the build-up of the final host-galaxy mass. We suggest that the observed low molecular gas fraction and short depletion timescale are due to AGN feedback, whose presence is indicated by fast AGN-driven ionised outflows in all our targets. © ESO 2020., This paper makes use of the following ALMA data: ADS/JAO.ALMA#2013.1.00417.S, ADS/JAO.ALMA#2019.1.01070.S, ADS/JAO.ALMA#2015.1.01602.S, ADS/JAO.ALMA#2016.1.00718.S, ADS/JAO. ALMA#2016.1.01515.S, ADS/JAO.ALMA#2018.1.01806.S. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada), MOST and ASIAA (Taiwan), and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ. This work is based on observations carried out under project number S17BW and W17DT with the IRAM NOEMA Interferometer. IRAM is supported by INSU/CNRS (France), MPG (Germany) and IGN (Spain). We are grateful to the anonymous referee for useful feedback which helped us to improve the paper. We thank Dr. Y. Wang for providing us with SUBARU images of J0209-0005 and Dr. P. Santini for useful discussion about SFE in high-z galaxies. MBi, MBr, CF, FF, AM, EP acknowledge support from PRIN MIUR project "Black Hole winds and the Baryon Life Cycle of Galaxies: the stone-guest at the galaxy evolution supper", contract #2017PH3WAT. MBi, CF and FF acknowledge support from INAF under PRIN SKA/CTA FORECaST and PRIN MAINSTREAM 2018 "Black hole winds and the baryon cycle". EP, LZ and MBi acknowledge financial support under ASI-INAF contract 2017-14-H.0. GB acknowledges financial support under the INTEGRAL ASI-INAF agreement 2019-35-HH.0 MPT acknowledges financial support from the State Agency for Research of the Spanish MCIU through the Center of Excellence Severo Ochoa award to the Instituto de Astrofisica de Andalucia (SEV-2017-0709) and through grant PGC2018-098915-B-C21 (MCI/AEI/FEDER, UE). RM acknowledges ERC Advanced Grant 695671 "QUENCH" and support by the Science and Technology Facilities Council (STFC). GVe acknowledges support from CONICYT Basal-CATA AFB-170002 and FONDECYT Postdoctorado 3200802 grants. CC acknowledges support from the Royal Society. MP is supported by the Programa Atraccion de Talento de la Comunidad de Madrid via grant 2018-T2/TIC-11715. GVi acknowledges financial support from Premiale 2015 MITiC (PI B. Garilli).

We present $NuSTAR$ observations of four active galactic nuclei (AGN) located within 15 Mpc. These AGN, namely ESO 121-G6, NGC 660, NGC 3486 and NGC 5195, have observed X-ray luminosities of $L_{\rm 2-10\ keV, obs} \lesssim$ 10$^{39}$ erg s$^{-1}$, classifying them as low luminosity AGN (LLAGN). We perform broadband X-ray spectral analysis for the AGN by combining our $NuSTAR$ data with $Chandra$ or $XMM-Newton$ observations to directly measure their column densities ($N_{\rm H}$) and infer their intrinsic power. We complement our X-ray data with archival and new high angular resolution mid-infrared (mid-IR) data for all objects, except NGC 5195. Based on our X-ray spectral analysis, we found that both ESO 121-G6 and NGC 660 are heavily obscured ($N_{\rm H}$ > 10$^{23}$ cm$^{-2}$; $L_{\rm 2-10\ keV,\ int} \sim$ 10$^{41}$ erg s$^{-1}$), and NGC 660 may be Compton-thick. We also note that the X-ray flux and spectral slope for ESO 121-G6 have significantly changed over the last decade, indicating significant changes in the obscuration and potentially accretion rate. On the other hand, NGC 3486 and NGC 5195 appear to be unobscured and just mildly obscured, respectively, with $L_{\rm 2-10\ keV,\ int} $ 10$^{41}$ erg s$^{-1}$ and $\lambda_{\rm Edd} \gtrsim$ 10$^{-3}$, and are detected in high angular resolution mid-IR imaging, indicating the presence of obscuring dust on nuclear scale. NGC 3486 however, is undetected in high-resolution mid-IR imaging, and the current data do not provide stringent constraints on the presence or absence of obscuring nuclear dust in the AGN., Comment: 18 pages, 11 figures, 3 tables, accepted for publication in MNRAS

During the last years, Ly$\alpha$ nebulae have been routinely detected around high-z, radio-quiet quasars (RQQs) thanks to the advent of sensitive integral field spectrographs. Constraining the physical properties of the Ly$\alpha$ nebulae is crucial for a full understanding of the circum-galactic medium (CGM), which is a venue of feeding and feedback processes. The most luminous quasars are privileged test-beds to study these processes, given their large ionizing fluxes and dense CGM environments in which they are expected to be embedded. We aim at characterizing the rest-frame UV emission lines in the CGM around a hyper-luminous, broad emission line, RQQ at z~3.6, that exhibits powerful outflows at both nuclear and host galaxy scales. We analyze VLT/MUSE observations of the quasar J1538+08 and perform a search for extended UV emission lines to characterize its morphology, emissivity, kinematics and metal content. We report the discovery of a very luminous ($\sim2 \times10^{44}~erg~s^{-1}$), giant Ly$\alpha$ nebula and a likely associated extended CIV nebula. The Ly$\alpha$ nebula emission exhibits moderate blueshift compared with the quasar systemic redshift and large average velocity dispersion ($\sigma_{v}$ ~700 $km~s^{-1}$) across the nebula, while the CIV nebula shows $\sigma_{v}$~$350~km~s^{-1}$. The Ly$\alpha$ line profile exhibits a significant asymmetry towards negative velocity values at 20-30 kpc south of the quasar and is well parameterized by two Gaussian components: a narrow ($\sigma$~$470~km~s^{-1}$) systemic one plus a broad ($\sigma$~1200 $km~s^{-1}$), blueshifted (~1500 $km~s^{-1}$) one. Our analysis of the MUSE observation of J1538+08 reveals metal-enriched CGM around this hyper-luminous quasar and our detection of blueshifted emission in the line profile of the Ly$\alpha$ nebula suggests that powerful nuclear outflows can propagate through the CGM over tens of kpc., Comment: 14 pages, 16 figures, 2 tables. Accepted for publication in A&A Journal

Aims. The SINFONI survey for Unveiling the Physics and Effect of Radiative feedback (SUPER) was designed to conduct a blind search for AGN-driven outflows on X-ray-selected AGNs at redshift z∼2 with high (∼2 kpc) spatial resolution, and to correlate them with the properties of their host galaxy and central black hole. The main aims of this paper are: (a) to derive reliable estimates for the masses of the black holes and accretion rates for the Type-1 AGNs in this survey; and (b) to characterise the properties of the AGN-driven winds in the broad line region (BLR). Methods. We analysed rest-frame optical and UV spectra of 21 Type-1 AGNs. We used Hα, Hβ, and MgII line profiles to estimate the masses of the black holes. We used the blueshift of the CIV line profile to trace the presence of winds in the BLR. Results. We find that the Hα and Hβ line widths are strongly correlated, as is the line continuum luminosity at 5100 Å with Hα line luminosity, resulting in a well-defined correlation between black hole masses estimated from Hα and Hβ. Using these lines, we estimate that the black hole masses for our objects are in the range Log (MBH/M·) = 8.4-10.8 and are accreting at λEdd = 0.04-1.3. Furthermore, we confirm the well-known finding that the CIV line width does not correlate with the Balmer lines and the peak of the line profile is blueshifted with respect to the [OIII]-based systemic redshift. These findings support the idea that the CIV line is tracing outflowing gas in the BLR for which we estimated velocities up to ∼4700 km s-1. We confirm the strong dependence of the BLR wind velocity on the UV-to-X-ray continuum slope, the bolometric luminosity, and Eddington ratio. We infer BLR mass outflow rates in the range 0.005-3 M· yr-1, revealing a correlation with the bolometric luminosity consistent with that observed for ionised winds in the narrow line region (NLR), and X-ray winds detected in local AGNs, and kinetic power ∼10-7-10-4 × LBol. The coupling efficiencies predicted by AGN-feedback models are much higher than the values reported for the BLR winds in the SUPER sample; although it should be noted that only a fraction of the energy injected by the AGN into the surrounding medium is expected to become kinetic power in the outflow. Finally, we find an anti-correlation between the equivalent width of the [OIII] line and the CIV velocity shift, and a positive correlation between this latter parameter and [OIII] outflow velocity. These findings, for the first time in an unbiased sample of AGNs at z∼2, support a scenario where BLR winds are connected to galaxy-scale detected outflows, and are therefore capable of affecting the gas in the NLR located at kiloparsec scale distances., With funding from the Spanish government through the "María de Maeztu Unit of Excellence" accreditation (MDM-2017-0737)

Feedback from AGN is thought to be key in shaping the life cycle of their host galaxies by regulating star-formation activity. Therefore, to understand the impact of AGN on star formation, it is essential to trace the molecular gas out of which stars form. In this paper we present the first systematic study of the CO properties of AGN hosts at z~2 for a sample of 27 X-ray selected AGN spanning two orders of magnitude in AGN bolometric luminosity (Lbol= 10^44.7-10^46.9 erg/s) by using ALMA Band 3 observations of the CO(3-2) transition (~1" angular resolution). To search for evidence of AGN feedback on the CO properties of the host galaxies, we compared our AGN with a sample of inactive (i.e., non-AGN) galaxies from the PHIBSS survey with similar redshift, stellar masses, and SFRs. We used the same CO transition as a consistent proxy for the gas mass for the two samples in order to avoid systematics involved when assuming conversion factors. By adopting a Bayesian approach to take upper limits into account, we analyzed CO luminosities as a function of stellar masses and SFRs, as well as the ratio LCO(3-2)/M* (proxy for the gas fraction). The two samples show statistically consistent trends in the LCO(3-2)-Lfir and LCO(3-2)-M* planes. However, there are indications that AGN feature lower CO(3-2) luminosities (0.4-0.7 dex) than inactive galaxies at the 2-3sigma level when we focus on the subset of parameters where the results are better constrained and on the distribution of the mean LCO(3-2)/M*. Therefore, even by conservatively assuming the same excitation factor r31, we would find lower molecular gas masses in AGN, and assuming higher r31 would exacerbate this difference. We interpret our result as a hint of the potential effect of AGN activity (e.g., radiation and outflows), which may be able to heat, excite, dissociate, and/or deplete the gas reservoir of the host galaxies. (abridged), Comment: Key figures: 2, 3, 4. 22 pages, 9 figures (5 in appendices), 3 tables. Accepted for publication in A&A

The SINFONI survey for Unveiling the Physics and Effect of Radiative feedback (SUPER) aims at tracing and characterizing ionized gas outflows and their impact on star formation in a statistical sample of X-ray selected Active Galactic Nuclei (AGN) at z$\sim$2. We present the first SINFONI results for a sample of 21 Type-1 AGN spanning a wide range in bolometric luminosity (log $\mathrm{L_{bol}}$ = 45.4-47.9 erg/s). The main aims of this paper are determining the extension of the ionized gas, characterizing the occurrence of AGN-driven outflows, and linking the properties of such outflows with those of the AGN. We use Adaptive Optics-assisted SINFONI observations to trace ionized gas in the extended narrow line region using the [OIII]5007 line. We classify a target as hosting an outflow if its non-parametric velocity of the [OIII] line, $\mathrm{w_{80}}$, is larger than 600 km/s. We study the presence of extended emission using dedicated point-spread function (PSF) observations, after modelling the PSF from the Balmer lines originating from the Broad Line Region. We detect outflows in all the Type-1 AGN sample based on the $\mathrm{w_{80}}$ value from the integrated spectrum, which is in the range 650-2700 km/s. There is a clear positive correlation between $\mathrm{w_{80}}$ and the AGN bolometric luminosity (99% correlation probability), but a weaker correlation with the black hole mass (80% correlation probability). A comparison of the PSF and the [OIII] radial profile shows that the [OIII] emission is spatially resolved for $\sim$35% of the Type-1 sample and the outflows show an extension up to $\sim$6 kpc. The relation between maximum velocity and the bolometric luminosity is consistent with model predictions for shocks from an AGN driven outflow. The escape fraction of the outflowing gas increase with the AGN luminosity, although for most galaxies, this fraction is less than 10%., Accepted for publication in A&A, 34 pages, 14 figures and 5 tables

We present the discovery of copious molecular gas in the halo of cid_346, a $z=2.2$ quasar studied as part of the SINFONI survey for Unveiling the Physics and Effect of Radiative feedback (SUPER). New Atacama Compact Array (ACA) CO(3-2) observations detect a much higher flux (by a factor of $14\pm5$) than measured on kiloparsec scales ($r\lesssim8$ kpc) using previous snapshot Atacama Large Millimeter/submillimeter Array data. Such additional CO(3-2) emission traces a structure that extends out to $r\sim200$ kpc in projected size, as inferred through direct imaging and confirmed by an analysis of the uv visibilities. This is the most extended molecular circumgalactic medium (CGM) reservoir that has ever been mapped. It shows complex kinematics, with an overall broad line profile (FWHM $= 1000$ km/s) that is skewed towards redshifted velocities up to at least $v\sim1000$ km/s. Using the optically thin assumption, we estimate a strict lower limit for the total molecular CGM mass observed by ACA of $M_{mol}^{CGM}>10^{10}~M_{\odot}$. There is however room for up to $M^{CGM}_{mol}\sim 1.7\times 10^{12}$ $M_{\odot}$, once optically thick CO emission with $\alpha_{\rm CO}=3.6$ $\rm M_{\odot}~(K~km~s^{-1}~pc^2)^{-1}$ and $L^{\prime}_{CO(3-2)}/L^{\prime}_{CO(1-0)}=0.5$ are assumed. Since cid_346 hosts quasar-driven ionised outflows and since there is no evidence of merging companions or an overdensity, we suggest that outflows may have played a crucial rule in seeding metal-enriched, dense gas on halo scales. However, the origin of such an extended molecular CGM remains unclear., Comment: Accepted for publication in A&A Letters. 9 pages and 9 figures including Appendix. Version updated after language editing

The WISSH quasars project aims at studying the nuclear and host galaxy properties of the most luminous quasars ($L_{bol}>10^{47}$ erg/s, $1.8, Comment: Accepted for publication on A&A

We report on the first ALMA observation of the CO(3$-$2) and rest-frame ~340 GHz continuum emission in PDS 456, which is the most luminous, radio-quiet QSO in the local Universe ($z$~0.18), with a bolometric luminosity $L_{\rm Bol}\sim10^{47}$ erg s$^{-1}$. ALMA angular resolution allowed us to map scales as small as ~700 pc. The molecular gas reservoir, traced by the core of the very bright CO(3$-$2) emission line, is distributed in a compact rotating disk, with size of ~1.3 kpc, seen close to face-on ($i$~25 deg). Fast CO(3$-$2) emission in the velocity range $v\in[-1000,+500]$ km s$^{-1}$ is also present. Specifically, we detect several blue-shifted clumps out to ~5 kpc from the nucleus, in addition to a compact ($R\lesssim1.2$ kpc), broad emission component. These components reveal a galaxy-wide molecular outflow, with a total mass $M_{\rm mol}^{\rm out}\sim2.5\times10^8$ $M_{\odot}$ and a mass outflow rate $\dot{M}_{\rm mol}\sim290$ $M_{\odot}$ yr$^{-1}$. The corresponding depletion time is ~8 Myr, shorter than the rate at which the molecular gas is converted into stars, indicating that the detected outflow is potentially able to quench star-formation in the host. The momentum flux of the molecular outflow normalised to the radiative momentum output (i.e. $L_{\rm Bol}/c$) is $\lesssim1$, comparable to that of the X-ray ultra-fast outflow (UFO) detected in PDS 456. This is at odds with the expectations for an energy-conserving expansion suggested for most of the large-scale outflows detected in low-luminosity AGN so far. We suggest three possible scenarios that may explain this observation: (i) in very luminous AGN such as our target the molecular gas phase is tracing only a fraction of the total outflowing mass; (ii) a small coupling between the shocked gas by the UFO and the host-galaxy ISM (iii) AGN radiation pressure may play an important role in driving the outflow., Accepted in A&A

Ultra Fast Outflows (UFOs) are an established feature in X-ray spectra of AGNs. According to the standard picture, they are launched at accretion disc scales with relativistic velocities, up to 0.3-0.4 c. Their high kinetic power is enough to induce an efficient feedback on galactic-scale, possibly contributing to the co-evolution between the central supermassive black hole (SMBH) and the host galaxy. It is therefore of paramount importance to fully understand the UFO physics, in particular the forces driving their acceleration and the relation with the accretion flow they originate from. In this paper we investigate the impact of special relativity effects on the radiative pressure exerted onto the outflow. The radiation received by the wind decreases for increasing outflow velocity v, implying that the standard Eddington limit argument has to be corrected according to v. Due to the limited ability of the radiation to counteract the SMBH gravity, we expect to find lower typical velocities with respect to the non-relativistic scenario. We integrate the relativistic-corrected outflow equation of motion for a realistic set of starting conditions. We concentrate on UFO typical values of ionisation, column density and launching radius. We explore a one-dimensional, spherical geometry and a 3D setting with a rotating thin accretion disc. We find that the inclusion of relativistic effects leads to sizeable differences in the wind dynamics and that v is reduced up to 50% with respect to the non-relativistic treatment. We compare our results with a sample of UFO from the literature, and we find that the relativistic-corrected velocities are systematically lower than the reported ones, indicating the need for an additional mechanism, such as magnetic driving, to explain the highest velocity components. These conclusions, derived for AGN winds, have a general applicability., Comment: 16 pages, 15 figures. Shortened abstract. Accepted for publication in A&A

Outflows are observed in a variety of astrophysical sources. Remarkably, ultra-fast ($v\geq 0.1c$), outflows in the UV and X-ray bands are often seen in AGNs. Depending on their energy and mass outflow rate, respectively $\dot{E}_{out}, \dot{M}_{out}$, such outflows may play a key role in regulating the AGN-host galaxy co-evolution process through cosmic time. It is therefore crucial to provide accurate estimates of the wind properties. Here, we concentrate on special relativistic effects concerning the interaction of light with matter moving at relativistic speed relatively to the source of radiation. Our aim is to assess the impact of these effects on the observed properties of the outflows and implement a relativistic correction in the existing spectral modelling routines. We define a simple procedure to incorporate relativistic effects in radiative transfer codes. Following this procedure, we run a series of simulations to explore the impact of these effects on the simulated spectra, for different $v$ and column densities of the outflow. The observed optical depth is usually considered a proxy for the wind $N_H$, independently on its velocity. However, our simulations show that the observed optical depth of an outflow with a given column density $N_H$ decreases rapidly as the velocity of the wind approaches relativistic values. This, in turn, implies that when estimating $N_H$ from the optical depth, it is necessary to include a velocity-dependent correction, already for moderate velocities (e.g. $v \geq 0.05c$). This correction linearly propagates to the derived $\dot{M}_{out}, \dot{E}_{out}$. As an example of these effects, we calculate the relativistically corrected values of $\dot{M}_{out}$ and $\dot{E}_{out}$ for a sample of $\sim 30$ Ultra-Fast Outflows taken from the literature, and find correction factors of $20-120 \%$ within the observed range of outflowing velocities., Comment: 6 pages, 5 figures. Accepted for publication in A&A

We present an analysis of a sample of 69 local obscured Swift/Burst Alert Telescope active galactic nuclei (AGNs) with X-ray spectra from NuSTAR and infrared (IR) spectral energy distributions from Herschel and WISE. We combine this X-ray and IR phenomenological modeling and find a significant correlation between reflected hard X-ray emission and IR AGN emission, with suggestive indications that this correlation may be stronger than the one between intrinsic hard X-ray and IR emissions. This relation between the IR and reflected X-ray emission suggests that both are the result of the processing of intrinsic emission from the corona and accretion disk by the same structure. We explore the resulting implications on the underlying distribution of covering fraction for all AGNs, by generating mock observables for the reflection parameter and IR luminosity ratio using empirical relations found for the covering fraction with each quantity. We find that the observed distributions of the reflection parameter and IR-to-X-ray ratio are reproduced with broad distributions centered around covering fractions of at least ~40%-50%, whereas narrower distributions match our observations only when centered around covering fractions of ~70%-80%. Our results are consistent with both independent estimates of the covering fractions of individual objects and the typical covering fraction obtained on the basis of obscured fractions for samples of AGNs. These results suggest that the level of reprocessing in AGNs, including X-ray reflection, is related in a relatively straightforward way to the geometry of the obscuring material., Published in the Astrophysical Journal: Lanz, L. et al. 2019, ApJ, 870, 26; 19 pages; 9 figures

Outflows from active galactic nuclei (AGN) are often invoked to explain the co-evolution of AGN and their host galaxies, and the scaling relations between the central black hole mass and the bulge velocity dispersion. Nuclear winds are often seen in the X-ray spectra through Fe K shell transitions and some of them are called ultra fast outflows (UFOs) due to their high velocities, up to some fractions of the speed of light. If they were able to transfer some percentage of the AGN luminosity to the host galaxy, this might be enough to trigger an efficient feedback mechanism. We aim to establish new constraints on the covering fraction and on the kinematic properties of the UFO in the powerful (L(bol) ~ 10^(47) erg/s) quasar PDS 456, an established Rosetta stone for studying AGN feedback from disk winds. This will allow us to estimate the mass outflow rate and the energy transfer rate of the wind, which are key quantities to understand the potential impact on the host galaxy. We analyze two sets of simultaneous XMM-Newton and NuSTAR observations taken in September 2013 and reported in Nardini et al. (2015) as having similar broadband spectral properties. We fit the Fe K features with a P-Cygni profile between 5 and 14 keV, using a novel Monte Carlo model for the WINd Emission (WINE). We find an outflow velocity ranging from 0.17 to 0.28 c, with a mean value of 0.23 c. We obtain an opening angle of the wind of 71(+13,-8) deg and a covering fraction of 0.7(+0.2,-0.3), suggesting a wide-angle outflow. We check the reliability of the WINE model by performing extensive simulations of joint XMM-Newton and NuSTAR observations. Furthermore, we test the accuracy of the WINE model in recovering the geometrical properties of UFOs by simulating observations with the forthcoming X-ray observatory ATHENA., Comment: 7 pages, 4 figures, accepted for publication in A&A

Theoretical models of galaxy formation suggest that the presence of an active galactic nucleus (AGN) is required to regulate the growth of its host galaxy through feedback mechanisms, produced by e.g. AGN-driven outflows. Although such outflows are common both at low and high redshift, a comprehensive picture is still missing. The peak epoch of galaxy assembly (1, Comment: 28 pages, 8 figures, 5 tables. Accepted for publication in A&A

It has been known for decades that the observed number of baryons in the local universe falls about 30-40% short of the total number of baryons predicted by Big-Bang Nucleosynthesis, as inferred from density fluctuations of the Cosmic Microwave Background and seen during the first 2-3 billion years of the universe in the so called Lyman-alpha Forest. A theoretical solution to this paradox locates the missing baryons in the hot and tenuous filamentary gas between galaxies, known as the warm-hot intergalactic medium. However, it is difficult to detect them there because the largest by far constituent of this gas - hydrogen - is mostly ionized and therefore almost invisible in far-ultraviolet spectra with typical signal-to-noise ratios. Indeed, despite the large observational efforts, only a few marginal claims of detection have been made so far. Here we report observations of two absorbers of highly ionized oxygen (OVII) in the high signal-to-noise-ratio X-ray spectrum of a quasar at redshift >0.4. These absorbers show no variability over a 2-year timescale and have no associated cold absorption, making the assumption that they originate from the quasar's intrinsic outflow or the host galaxy's interstellar medium implausible. The OVII systems lie in regions characterized by large (x4 compared to average) galaxy over-densities and their number (down to the sensitivity threshold of our data), agrees well with numerical simulation predictions for the long-sought warm-hot intergalactic medium (WHIM). We conclude that the missing baryons have been found., Appeared in Nature (Volume 558, Issue 7710) on 21 June 2018. The posted PDF version is the pre-editorial-change version and includes the main paper, its Methods section and the Extended Data section. A link to the (view-only) PDF of the final published version of the paper, is available here: https://rdcu.be/1eak

Observing the light emitted by the first accreting black holes (BHs) would dramatically improve our understanding of the formation of quasars at z > 6, possibly unveiling the nature of their supermassive black hole (SMBH) seeds. In previous works we explored the relative role of the two main competing BH seed formation channels, Population III remnants (low-mass seeds) and direct collapse BHs (high-mass seeds), investigating the properties of their host galaxies in a cosmological context. Building on this analysis, we predict here the spectral energy distribution and observational features of low- and high-mass BH seeds selected among the progenitors of a z~6 SMBH. We derive the processed emission from both accreting BHs and stars by using the photo-ionization code Cloudy, accounting for the evolution of metallicity and dust-to-gas mass ratio in the interstellar medium of the host galaxies, as predicted by the cosmological data- constrained model GAMETE/QSOdust. We show how future missions like JWST and ATHENA will be able to detect the light coming from SMBH progenitors already at z~16. We build upon previous complementary studies and propose a method based on the combined analysis of near infrared (NIR) colors, IR excess (IRX) and UV continuum slopes (i.e. color-color and IRX-Beta diagrams) to distinguish growing seed BH host galaxies from starburst-dominated systems in JWST surveys. Sources selected through this criterion would be the best target for follow-up X-ray observations., Comment: accepted for publicaztion in MNRAS

We present ALMA observations of the CO(6-5) and [CII] emission lines and the sub-millimeter continuum of the $z\sim6$ quasi-stellar object (QSO) SDSS J231038.88+185519.7. Compared to previous studies, we have analyzed a synthetic beam that is ten times smaller in angular size, we have achieved ten times better sensitivity in the CO(6-5) line, and two and half times better sensitivity in the [CII] line, enabling us to resolve the molecular gas emission. We obtain a size of the dense molecular gas of $2.9\pm0.5$ kpc, and of $1.4\pm0.2$ kpc for the 91.5 GHz dust continuum. By assuming that CO(6-5) is thermalized, and by adopting a CO--to--$H_2$ conversion factor $\rm \alpha_{CO} = 0.8~ M_{\odot}~K^{-1}~ (km/s)^{-1} ~pc^{2}$, we infer a molecular gas mass of $\rm M(H_2)=(3.2 \pm0.2) \times 10^{10}\rm M_{\odot}$. Assuming that the observed CO velocity gradient is due to an inclined rotating disk, we derive a dynamical mass of $\rm M_{dyn}~sin^2(i) = (2.4\pm0.5) \times 10^{10}~ M_{\odot}$, which is a factor of approximately two smaller than the previously reported estimate based on [CII]. Regarding the central black hole, we provide a new estimate of the black hole mass based on the C~IV emission line detected in the X-SHOOTER/VLT spectrum: $\rm M_{BH}=(1.8\pm 0.5) \times 10^{9}~ M_{\odot}$. We find a molecular gas fraction of $\rm \mu=M(H_2)/M^*\sim4.4$, where $\rm M^*\approx M_{dyn} - M(H_2)-M(BH)$. We derive a ratio $v_{rot}/\sigma \approx 1-2$ suggesting high gas turbulence, outflows/inflows and/or complex kinematics due to a merger event. We estimate a global Toomre parameter $Q\sim 0.2-0.5$, indicating likely cloud fragmentation. We compare, at the same angular resolution, the CO(6-5) and [CII] distributions, finding that dense molecular gas is more centrally concentrated with respect to [CII]. We find that the current BH growth rate is similar to that of its host galaxy., Comment: A&A in press

We present an ALMA high-resolution observation of the 840 um continuum and [CII] line emission in the WISE-SDSS selected hyper-luminous (WISSH) QSO J1015+0020 at z~4.4. Our analysis reveals an exceptional overdensity of [CII]-emitting companions with a very small (, Accepted

Studying the coupling between the energy output produced by the central quasar and the host galaxy is fundamental to fully understand galaxy evolution. Quasar feedback is indeed supposed to dramatically affect the galaxy properties by depositing large amounts of energy and momentum into the ISM. In order to gain further insights on this process, we study the SEDs of sources at the brightest end of the quasar luminosity function, for which the feedback mechanism is supposed to be at its maximum. We model the rest-frame UV-to-FIR SEDs of 16 WISE-SDSS Selected Hyper-luminous (WISSH) quasars at 1.8 < z < 4.6 disentangling the different emission components and deriving physical parameters of both the nuclear component and the host galaxy. We also use a radiative transfer code to account for the contribution of the quasar-related emission to the FIR fluxes. Most SEDs are well described by a standard combination of accretion disk+torus and cold dust emission. However, about 30% of them require an additional emission component in the NIR, with temperatures peaking at 750K, which indicates the presence of a hotter dust component in these powerful quasars. We measure extreme values of both AGN bolometric luminosity (LBOL > 10^47 erg/s) and SFR (up to 2000 Msun/yr). A new relation between quasar and star-formation luminosity is derived (LSF propto LQSO^(0.73)) by combining several Herschel-detected quasar samples from z=0 to 4. Future observations will be crucial to measure the molecular gas content in these systems, probe the impact between quasar-driven outflows and on-going star-formation, and reveal the presence of merger signatures in their host galaxies., 19 pages, 12 figures; Accepted for publication in Astronomy & Astrophysics on June 13, 2017

We present the 78-ks Chandra observations of the $z=6.4$ quasar SDSS J1148+5251. The source is clearly detected in the energy range 0.3-7 keV with 42 counts (with a significance $\gtrsim9\sigma$). The X-ray spectrum is best-fitted by a power-law with photon index $\Gamma=1.9$ absorbed by a gas column density of $\rm N_{\rm H}=2.0^{+2.0}_{-1.5}\times10^{23}\,\rm cm^{-2}$. We measure an intrinsic luminosity at 2-10 keV and 10-40 keV equal to $\sim 1.5\times 10^{45}~\rm erg~s^{-1}$, comparable with luminous local and intermediate-redshift quasar properties. Moreover, the X-ray to optical power-law slope value ($\alpha_{\rm OX}=-1.76\pm 0.14$) of J1148 is consistent with the one found in quasars with similar rest-frame 2500 \AA ~luminosity ($L_{\rm 2500}\sim 10^{32}~\rm erg~s^{-1}$\AA$^{-1}$). Then we use Chandra data to test a physically motivated model that computes the intrinsic X-ray flux emitted by a quasar starting from the properties of the powering black hole and assuming that X-ray emission is attenuated by intervening, metal-rich ($Z\geq \rm Z_{\odot}$) molecular clouds distributed on $\sim$kpc scales in the host galaxy. Our analysis favors a black hole mass $M_{\rm BH} \sim 3\times 10^9 \rm M_\odot$ and a molecular hydrogen mass $M_{\rm H_2}\sim 2\times 10^{10} \rm M_\odot$, in good agreement with estimates obtained from previous studies. We finally discuss strengths and limits of our analysis., Comment: 9 pages, 3 figures, 1 table, MNRAS in press

We present a study of the average X-ray spectral properties of the sources detected by the NuSTAR extragalactic survey, comprising observations of the E-CDFS, EGS and COSMOS fields. The sample includes 182 NuSTAR sources (64 detected at 8-24 keV), with 3-24 keV fluxes ranging between $f_{\rm 3-24 keV}\approx10^{-14}$ and $6\times10^{-13}$ erg/cm$^2$/s ($f_{\rm 8-24 keV}\approx3\times10^{-14}-3\times10^{-13}$ erg/cm$^2$/s) and redshifts of $z=0.04-3.21$. We produce composite spectra from the Chandra+NuSTAR data ($E\approx2-40$ keV, rest frame) for all the sources with redshift identifications (95%) and investigate the intrinsic, average spectra of the sources, divided into broad-line (BL) and narrow-line (NL) AGN, and also in different bins of X-ray column density and luminosity. The average power-law photon index for the whole sample is $\Gamma=1.65_{-0.03}^{+0.03}$, flatter than $\Gamma\approx1.8$ typically found for AGN. While the spectral slope of BL and X-ray unabsorbed AGN is consistent with typical values ($\Gamma=1.79_{-0.01}^{+0.01}$), a significant flattening is seen in NL AGN and heavily-absorbed sources ($\Gamma=1.60_{-0.05}^{+0.08}$ and $\Gamma=1.38_{-0.12}^{+0.12}$, respectively), likely due to the effect of absorption and to the contribution from Compton reflection to the high-energy flux (E>10 keV). We find that the typical reflection fraction in our spectra is $R\approx0.5$ (for $\Gamma=1.8$), with a tentative indication of an increase of the reflection strength with column density. While there is no significant evidence for a dependence of the photon index with X-ray luminosity in our sample, we find that $R$ decreases with luminosity, with relatively high levels of reflection ($R\approx1.2$) for $L_{\rm 10-40 keV}10^{44}$ erg/s AGN, assuming $\Gamma=1.8$., Comment: 19 pages, 14 figures, 4 tables. Accepted for publication in ApJ

We have performed a high sensitivity observation of the UFO/BAL quasar APM 08279+5255 at z=3.912 with NOEMA at 3.2 mm, aimed at detecting fast moving molecular gas. We report the detection of blueshifted CO(4-3) with maximum velocity (v95\%) of $-1340$ km s$^{-1}$, with respect to the systemic peak emission, and a luminosity of $L' = 9.9\times 10^9 ~\mu^{-1}$ K km s$^{-1}$ pc$^{-2}$ (where $\mu$ is the lensing magnification factor). We discuss various scenarios for the nature of this emission, and conclude that this is the first detection of fast molecular gas at redshift $>3$. We derive a mass flow rate of molecular gas in the range $\rm \dot M=3-7.4\times 10^3$ M$_\odot$/yr, and momentum boost $\dot P_{OF} / \dot P_{AGN} \sim 2-6$, therefore consistent with a momentum conserving flow. For the largest $\dot P_{OF}$ the scaling is also consistent with a energy conserving flow with an efficiency of $\sim$10-20\%. The present data can hardly discriminate between the two expansion modes. The mass loading factor of the molecular outflow $\eta=\dot M_{OF}/SFR$ is $>>1$. We also detect a molecular emission line at a frequency of 94.83 GHz, corresponding to a rest frame frequency of 465.8 GHz, which we tentatively identified with the cation molecule $\rm N_2H^+$(5-4), which would be the first detection of this species at high redshift. We discuss the alternative possibility that this emission is due to a CO emission line from the, so far undetected, lens galaxy. Further observations of additional transitions of the same species with NOEMA can discriminate between the two scenarios., Comment: Accepted for publication on Astronomy & Astrophysics

Feedback from accreting SMBHs is often identified as the main mechanism responsible for regulating star-formation in AGN host galaxies. However, the relationships between AGN activity, radiation, winds, and star-formation are complex and still far from being understood. We study scaling relations between AGN properties, host galaxy properties and AGN winds. We then evaluate the wind mean impact on the global star-formation history, taking into account the short AGN duty cycle with respect to that of star-formation. We first collect AGN wind observations for 94 AGN with detected massive winds at sub-pc to kpc spatial scales. We then fold AGN wind scaling relations with AGN luminosity functions, to evaluate the average AGN wind mass-loading factor as a function of cosmic time. We find strong correlations between the AGN molecular and ionised wind mass outflow rates and the AGN bolometric luminosity. The power law scaling is steeper for ionised winds (slope 1.29+/-0.38) than for molecular winds (0.76+/-0.06), meaning that the two rates converge at high bolometric luminosities. The molecular gas depletion timescale and the molecular gas fraction of galaxies hosting powerful AGN winds are 3-10 times shorter and smaller than those of main-sequence galaxies with similar SFR, stellar mass and redshift. These findings suggest that, at high AGN bolometric luminosity, the reduced molecular gas fraction may be due to the destruction of molecules by the wind, leading to a larger fraction of gas in the atomic ionised phase. The AGN wind mass-loading factor $\eta=\dot M_{OF}/SFR$ is systematically higher than that of starburst driven winds. Our analysis shows that AGN winds are, on average, powerful enough to clean galaxies from their molecular gas only in massive systems at z, This work was supported by ASI/INAF contract I/009/10/0 and INAF PRIN 2011, 2012 and 2014. MB acknowledges support from the FP7 Career Integration Grant “eEASy” (CIG 321913). LZ acknowledges support from ASI/INAF grant I/037/12/0. CF acknowledges funding from the European Union Horizon 2020 research and innovation programme under the Marie SklodowskaCurie grant agreement No 664931. CC acknowledges funding from the European Union Horizon 2020 research and innovation programme under the Marie SklodowskaCurie grant agreement No 664931 and support from Swiss National Science Foundation Grants PP00P2 138979 and PP00P2 166159. RM acknowledges the ERC Advanced Grant 695671 QUENCH and support from the Science and Technology Facilities Council (STFC).

We present the first full catalog and science results for the NuSTAR serendipitous survey. The catalog incorporates data taken during the first 40 months of NuSTAR operation, which provide ~20Ms of effective exposure time over 331 fields, with an areal coverage of 13 sq deg, and 497 sources detected in total over the 3-24 keV energy range. There are 276 sources with spectroscopic redshifts and classifications, largely resulting from our extensive campaign of ground-based spectroscopic followup. We characterize the overall sample in terms of the X-ray, optical, and infrared source properties. The sample is primarily comprised of active galactic nuclei (AGNs), detected over a large range in redshift from z = 0.002 - 3.4 (median of = 0.56), but also includes 16 spectroscopically confirmed Galactic sources. There is a large range in X-ray flux, from log( f_3-24keV / erg s^-1 cm^-2 ) ~ -14 to -11, and in rest-frame 10-40 keV luminosity, from log( L_10-40keV / erg s^-1 ) ~ 39 to 46, with a median of 44.1. Approximately 79% of the NuSTAR sources have lower energy ( 10^44 erg s^-1) to ~80% at the lowest luminosities (Lx < 10^43 erg s^-1). Our optical spectroscopic analysis finds that the observed fraction of optically obscured AGNs (i.e., the Type 2 fraction) is F_Type2 = 53(+14-15)%, for a well-defined subset of the 8-24 keV selected sample. This is higher, albeit at a low significance level, than the Type 2 fraction measured for redshift- and luminosity-matched AGNs selected by, 49 pages, 23 figures, 8 tables. Accepted for publication in ApJ

We present the 3-8 keV and 8-24 keV number counts of active galactic nuclei (AGN) identified in the NuSTAR extragalactic surveys. NuSTAR has now resolved 33-39% of the X-ray background in the 8-24 keV band, directly identifying AGN with obscuring columns up to approximately 1e25 /cm2. In the softer 3-8 keV band the number counts are in general agreement with those measured by XMM-Newton and Chandra over the flux range 5e-15 < S(3 - 8 keV)/(erg/cm2/s) < 1e-12 probed by NuSTAR. In the hard 8-24 keV band NuSTAR probes fluxes over the range 2e-14 < S(8-24 keV)/(erg/cm2/s) < 1e-12, a factor of approximately 100 fainter than previous measurements. The 8-24 keV number counts match predictions from AGN population synthesis models, directly confirming the existence of a population of obscured and/or hard X-ray sources inferred from the shape of the integrated cosmic X-ray background. The measured NuSTAR counts lie significantly above simple extrapolation with a Euclidian slope to low flux of the Swift/BAT 15-55 keV number counts measured at higher fluxes S(15-55 keV) > 1e-11 erg/cm2/s, reflecting the evolution of the AGN population between the Swift/BAT local (z, 8 pages; accepted by the Astrophysical Journal

We compute the non-thermal emissions produced by relativistic particles accelerated by the AGN-driven shocks in NGC 1068, and we compare the model predictions with the observed gamma-ray and radio spectra . The former is contributed by pion decay, inverse Compton scattering, and bremsstrahlung, while the latter is produced by synchrotron radiation. We derive the gamma-ray and radio emissions by assuming the standard acceleration theory, and we discuss how our results compare with those corresponding to other commonly assumed sources of gamma-ray and radio emissions, like Supernova remnants (SNR) or AGN jets. We find that the AGN-driven shocks observed in the circumnuclear molecular disk of such a galaxy provide a contribution to the gamma-ray emission comparable to that provided by the starburst activity when standard particle acceleration efficiencies are assumed, while they can yield the whole gamma-ray emission only when the parameters describing the acceleration efficiency and the proton coupling with the molecular gas are tuned to values larger than those assumed in standard, SNR-driven shocks. We discuss the range of acceleration efficiencies (for protons and electrons) and of proton calorimetric fractions required to account for the observed gamma-ray emission in the AGN outflow model. We further compare the neutrino flux expected in our model with constraints from current experiments, and we provide predictions for the detections by the upcoming KM3NeT neutrino telescope. This analysis strongly motivates observations of NGC 1068 at >TeV energies with current and future Cherenkov telescopes in order to gain insight into the nature of the gamma-rays source., 12 pages, 4 figures, accepted for publication in A&A

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.

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

We present XMM RGS and Chandra LETG observations of the blazar, H 2356-309, located behind the Sculptor Wall, a large-scale galaxy structure expected to harbor high-density Warm-Hot Intergalactic Medium (WHIM). Our simultaneous analysis of the RGS and LETG spectra yields a 3-sigma detection of the crucial redshifted O vii K-alpha line with a column density (>~ 10^{16} cm^{-2}) consistent with similar large-scale structures produced in cosmological simulations. This represents the first detection of non-local WHIM from X-ray absorption studies where XMM and Chandra data are analyzed simultaneously and the absorber redshift is already known, thus providing robust evidence for the expected repository of the "missing baryons"., Accepted for Publication in The Astrophysical Journal, 7 pages, 2 figures

We present the concentration (c)-virial mass (M) relation of 39 galaxy systems ranging in mass from individual early-type galaxies up to the most massive galaxy clusters, (0.06-20) x 10^{14} M_sun. We selected for analysis the most relaxed systems possessing the highest quality data currently available in the Chandra and XMM public data archives. A power-law model fitted to the X-ray c-M relation requires at high significance (6.6 sigma) that c decreases with increasing M, which is a general feature of CDM models. The median and scatter of the c-M relation produced by the flat, concordance LCDM model (Omega_m=0.3, sigma_8=0.9) agrees with the X-ray data provided the sample is comprised of the most relaxed, early forming systems, which is consistent with our selection criteria. Holding the rest of the cosmological parameters fixed to those in the concordance model the c-M relation requires 0.76< sigma_8 99% conf.) both open CDM models and flat CDM models with Omega_m ~1. This result provides novel evidence for a flat, low-Omega_m universe with dark energy using observations only in the local (z << 1) universe. Possible systematic errors in the X-ray mass measurements of a magnitude ~10% suggested by CDM simulations do not change our conclusions., Accepted for Publication in ApJ; 13 pages, 4 figures; minor clarifications and updates; correlation coefficients corrected in Table 1 (correct values were used in the analysis in previous versions); conclusions unchanged

We present the first direct measurements of the rest-frame 10-40 keV X-ray luminosity function (XLF) of Active Galactic Nuclei (AGNs) based on a sample of 94 sources at 0.1 < z, Comment: 14 pages, accepted for publication in the Astrophysical Journal

We present a broad-band (~0.3-70 keV) spectral and temporal analysis of NuSTAR observations of the luminous infrared galaxy NGC 6240, combined with archival Chandra, XMM-Newton and BeppoSAX data. NGC 6240 is a galaxy in a relatively early merger state with two distinct nuclei separated by ~1."5. Previous Chandra observations have resolved the two nuclei, showing that they are both active and obscured by Compton-thick material. Although they cannot be resolved by NuSTAR, thanks to the unprecedented quality of the NuSTAR data at energies >10 keV, we clearly detect, for the first time, both the primary and the reflection continuum components. The NuSTAR hard X-ray spectrum is dominated by the primary continuum piercing through an absorbing column density which is mildly optically thick to Compton scattering (tau ~ 1.2, N_H ~ 1.5 x 10^(24) cm^-2). We detect moderate hard X-ray (> 10 keV) flux variability up to 20% on short (15-20 ksec) timescales. The amplitude of the variability is maximum at ~30 keV and is likely to originate from the primary continuum of the southern nucleus. Nevertheless, the mean hard X-ray flux on longer timescales (years) is relatively constant. Moreover, the two nuclei remain Compton-thick, although we find evidence of variability of the material along the line of sight with column densities N_H, Comment: 14 pages, 10 figures, accepted for publication in A&A

We present initial results and the source catalog from the NuSTAR survey of the Extended Chandra Deep Field South (hereafter, ECDFS) - currently the deepest contiguous component of the NuSTAR extragalactic survey program. The survey covers the full ~30 arcmin x 30 arcmin area of this field to a maximum depth of ~360 ks (~220 ks when corrected for vignetting at 3-24 keV), reaching sensitivity limits of ~1.3 x 10^-14 erg/cm2/s (3-8 keV), ~3.4 x 10^-14 erg/cm2/s (8-24 keV) and ~3.0 x 10^-14 erg/cm2/s (3-24 keV). Fifty four (54) sources are detected over the full field, although five of these are found to lie below our significance threshold once contaminating flux from neighboring (i.e., blended) sources is taken into account. Of the remaining 49 that are significant, 19 are detected in the 8-24 keV band. The 8-24 keV to 3-8 keV band ratios of the twelve sources that are detected in both bands span the range 0.39-1.7, corresponding to a photon index range of Gamma ~ 0.5-2.3, with a median photon index of 1.70 +/- 0.52. The redshifts of the 49 sources in our main sample span the range z = 0.21-2.7, and their rest-frame 10-40 keV luminosities (derived from the observed 8-24 keV fluxes) span the range L(10-40 keV) ~ (0.7-300) x 10^43 erg/s, sampling below the "knee" of the X-ray luminosity function out to z ~ 0.8-1. Finally, we identify one NuSTAR source that has neither a Chandra nor an XMM-Newton counterpart, but that shows evidence of nuclear activity at infrared wavelengths, and thus may represent a genuine, new X-ray source detected by NuSTAR in the ECDFS., Comment: 20 pages, published in the Astrophysical Journal

To provide the census of the sources contributing to the X-ray background peak above 10 keV, NuSTAR is performing extragalactic surveys using a three-tier "wedding cake" approach. We present the NuSTAR survey of the COSMOS field, the medium sensitivity and medium area tier, covering 1.7 deg2 and overlapping with both Chandra and XMM-Newton data. This survey consists of 121 observations for a total exposure of ~3 Ms. To fully exploit these data, we developed a new detection strategy, carefully tested through extensive simulations. The survey sensitivity at 20% completeness is 5.9, 2.9 and 6.4 x 10^-14 erg/cm2/s in the 3-24 keV, 3-8 keV and 8-24 keV bands, respectively. By combining detections in 3 bands, we have a sample of 91 NuSTAR sources with 10^42 -10^45.5 erg/s luminosities and redshift z=0.04-2.5. Thirty two sources are detected in the 8-24 keV band with fluxes ~100 times fainter than sources detected by Swift-BAT. Of the 91 detections, all but four are associated with a Chandra and/or XMM-Newton point-like counterpart. One source is associated with an extended lower energy X-ray source. We present the X-ray (hardness ratio and luminosity) and optical-to-X-ray properties. The observed fraction of candidate Compton-thick AGN measured from the hardness ratio is between 13-20%. We discuss the spectral properties of NuSTAR J100259+0220.6 (ID 330) at z=0.044, with the highest hardness ratio in the entire sample. The measured column density exceeds 10^24 /cm2, implying the source is Compton-thick. This source was not previously recognized as such without the >10 keV data., 20 pages, published in the Astrophysical Journal

Q0056-363 is the most powerful X-ray quasar known to exhibit a broad, likely relativistic iron line (Porquet & Reeves 2003). It has been observed twice by XMM-$Newton$, three and half years apart (July 2000 and December 2003). In the second observation, the UV and soft X-ray fluxes were fainter, the hard X-ray power law flatter, and the iron line equivalent width (EW) smaller than in the 2000 observation. These variations can all be explained, at least qualitatively, if the disc is truncated in the second observation. We report also on the possible detection of a transient, redshifted iron absorption line during the 2003 observation., Accepted for publication in A&A

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

We present a possible OVIII X-ray absorption line at $z=0.117 \pm 0.001$ which, if confirmed, will be the first one associated with a broad HI Ly$\beta$ (BLB: FWHM=$160^{+50}_{-30}$ km s$^{-1}$) absorber. The absorber lies along the line of sight to the nearby ($z=0.1372$) Seyfert 1 galaxy PKS~0558-504, consistent with being a WHIM filament. The X-ray absorber is marginally detected in two independent XMM-Newton spectra of PKS~0558-504, a long $\sim 600$ ks Guest-Observer observation and a shorter, $\sim 300$ ks total, calibration observation, with a combined single line statistical significance of 2.8$\sigma$ (2.7$\sigma$ and 1.2$\sigma$ in the two spectra, respectively). When fitted with our self-consistent hybrid-photoionization WHIM models, the combined XMM-{\em Newton} spectrum is consistent with the presence of OVIII K$\alpha$ at $z=(0.117 \pm 0.001)$. This model gives best fitting temperature and equivalent H column density of the absorber of log$T=6.56_{-0.17}^{+0.19}$ K, and logN$_H=(21.5 \pm 0.3) (Z/Z_{0.01\odot})^{-1}$ cm$^{-2}$. The statistical sigificance of this single X-ray detection is increased by the detection of broad and complex HI Ly$\beta$ absorption in archival FUSE spectra of PKS~0558-504, at redshifts $z=0.1183 \pm 0.0001$ consistent with the best-fitting redshift of the X-ray absorber. The single line statistical significance of this line is 4.1$\sigma$ (3.7$\sigma$ if systematics are considered), and thus the combined (HI+OVIII) statistical significance of the detection is of 5.0$\sigma$. The detection of both metal and H lines at a consistent redshift, in this hot absorbing system, allows us to speculate on its metallicity. By associating the bulk of the X-ray absorber with the BLB line detected in the FUSE spectrum at $z_{BLB}=0.1183 \pm 0.0001$, we obtain a metallicity of 1-4\% Solar., Comment: 37 pages, 12 figures, 3 Tables. Accepted for publication by the ApJ.

We report the discovery of a galaxy cluster serendipitously detected as an extended X-ray source in an offset observation of the group NGC 5044. The cluster redshift, z=0.281, determined from the optical spectrum of the brightest cluster galaxy, agrees with that inferred from the X-ray spectrum using the Fe K alpha complex of the hot ICM (z=0.27 +/- 0.01). Based on the 50 ks XMM observation, we find that within a radius of 383 kpc the cluster has an unabsorbed X-ray flux, f_X (0.5-2 keV) = 3.34 (+0.08, -0.13) x 10^{-13} erg/cm^2/s, a bolometric X-ray luminosity, L_X = 2.21 (+0.34, -0.19) x 10^{44} erg/s, kT = 3.57 +/- 0.12 keV, and metallicity, 0.60 +/- 0.09 solar. The cluster obeys the scaling relations for L_X and T observed at intermediate redshift. The mass derived from an isothermal NFW model fit is, M_vir = 3.89 +/- 0.35 x 10^{14} solar masses, with a concentration parameter, c = 6.7 +/- 0.4, consistent with the range of values expected in the concordance cosmological model for relaxed clusters. The optical properties suggest this could be a ``fossil cluster''., Comment: 5 pages, 4 colour figures, accepted for publication in ApJ

We present the results of spectroscopic observations of galaxies associated with the diffuse X-ray emitting structure discovered by Zappacosta et al. (2002). After measuring the redshifts of 161 galaxies, we confirm an overdensity of galaxies with projected dimensions of at least 2 Mpc, determine its spectroscopic redshift in z=0.401+/-0.002, and show that it is spatially coincident with the diffuse X-ray emission. This confirms the original claim that this X-ray emission has an extragalactic nature and is due to the Warm-Hot Intergalactic Medium (WHIM). We used this value of the redshift to compute the temperature of the emitting gas. The resulting value depends on the metallicity that is assumed for the IGM, and is constrained to be between 0.3 and 0.6 keV for metallicities between 0.05 and 0.3 solar, in good agreement with the expectations from the WHIM., 9 pages, A&A, in press, minor language changes

We present radial mass profiles within 0.3 r_vir for 16 relaxed galaxy groups-poor clusters (kT range 1-3 keV) selected for optimal mass constraints from the Chandra and XMM data archives. After accounting for the mass of hot gas, the resulting mass profiles are described well by a two-component model consisting of dark matter (DM), represented by an NFW model, and stars from the central galaxy. The stellar component is required only for 8 systems, for which reasonable stellar mass-to-light ratios (M/L_K) are obtained, assuming a Kroupa IMF. Modifying the NFW dark matter halo by adiabatic contraction does not improve the fit and yields systematically lower M/L_K. In contrast to previous results for massive clusters, we find that the NFW concentration parameter (c_vir) for groups decreases with increasing M_vir and is inconsistent with no variation at the 3 sigma level. The normalization and slope of the c_vir-M_vir relation are consistent with the standard LambdaCDM cosmological model with sigma_8 = 0.9. The small intrinsic scatter measured about the c_vir-M_vir relation implies the groups represent preferentially relaxed, early forming systems. The mean gas fraction (f =0.05 +/- 0.01) of the groups measured within an overdensity Delta=2500 is lower than for hot, massive clusters, but the fractional scatter (sigma_f/f=0.2) for groups is larger, implying a greater impact of feedback processes on groups, as expected., Comment: Accepted for Publication in ApJ; 30 pages, 9 figures. No changes from previous version

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., 16 pages, 17 figures. Accepted for publication in A&A