Search

We present spatially resolved VLT/SINFONI spectroscopy with adaptive optics of type-2 active galactic nuclei (AGN) from the SINFONI Survey for Unveiling the Physics and Effect of Radiative feedback (SUPER), which targeted X-ray bright ($L_{2-10 keV}\gtrsim10^{42}$ erg s$^{-1}$) AGN at Cosmic Noon ($z\sim2$). Our analysis of the rest-frame optical spectra unveils ionised outflows in all seven examined targets, as traced via [OIII]$\lambda$5007 line emission, moving at $v\gtrsim600$ km s$^{-1}$. In six objects these outflows are clearly spatially resolved and extend on 2-4 kpc scales, whereas marginally resolved in the remaining one. Interestingly, these SUPER type-2 AGN are all heavily obscured sources ($N_{H}\gtrsim10^{23}$ cm$^{-2}$) and host faster ionised outflows than their type-1 counterparts within the same range of bolometric luminosity ($L_{bol} \sim 10^{44.8-46.5}$ erg s$^{-1}$). SUPER has hence provided observational evidence that the type-1/type-2 dichotomy at $z\sim2$ might not be driven simply by projection effects, but might reflect two distinct obscuring life stages of active galaxies, as predicted by evolutionary models. Within this picture, SUPER type-2 AGN might be undergoing the 'blow-out' phase, where the large amount of obscuring material efficiently accelerates large-scale outflows via radiation pressure on dust, eventually unveiling the central active nucleus and signal the start of the bright, unobscured type-1 AGN phase. Moreover, the overall population of ionised outflows detected in SUPER has velocities comparable with the escape speed of their dark matter halos, and in general high enough to reach 30-50 kpc distances from the centre. These outflows are hence likely to sweep away the gas (at least) out of the baryonic disk and/or to heat the host gas reservoir, thus reducing and possibly quenching star formation., Comment: 21 pages, 13 figures. Key figure is 8. Accepted for publication in A&A

IRASF11119 is an ultra-luminous IR galaxy with post-merger morphology, hosting a type-1 QSO at z=0.189. Its 2013 Suzaku spectrum shows a prominent Ultra Fast Outflow (UFO) absorption feature (v_out~0.25c). In 2021, we obtained the first XMM-Newton long look of the target, coordinated with a simultaneous NuSTAR observation. The new high-quality data allow us to detect at P>99.8% c.l. multiple absorption features associated with the known UFO. Furthermore, an emission plus absorption feature at 1.1-1.3 keV reveals the presence of a blueshifted P-Cygni profile in the soft band. We associate the hard band features with blends of FeXXV and FeXXVI He$\alpha$-Ly$\alpha$ and He$\beta$-Ly$\beta$ line pairs and infer a large column (N$_H$~$10^{24}$ cm$^{-2}$) of highly ionized (log$\xi$~5) gas outflowing at v_out=0.27c. The 1 keV feature can be associated with a blend of Fe and Ne transitions, produced by a lower column (N$_H$~$10^{21}$ cm$^{-2}$) and ionization (log$\xi$~2.6) gas component outflowing at the same speed. Using a radiative-transfer disk wind model to fit the highly ionized UFO, we derive a large mass outflow rate, comparable with the mass accretion rate (M$_{out}$=4.25 M$_{Sun}$/yr, ~1.6 M$_{acc}$), and kinetic energy and momentum flux among the highest reported in the literature. We measure an extremely low high-energy cut-off (E$_c$~25 keV). Several other cases in the literature suggest that a steep X-ray continuum may be related to the formation of powerful winds. The lack of a significant momentum boost between the nuclear UFO and the different phases of the large-scale outflow, observed in IRASF11119 and in a growing number of sources with powerful UFOs, can be explained by (i) a momentum-driven expansion, (ii) an inefficient coupling of the UFO with the host ISM, or (iii) by repeated energy-driven expansion episodes with low duty-cycle, that average out on long time-scales., Comment: 21 pages, 13 figures, accepted for publication in A&A

The spectral energy distribution (SED) of low-luminosity active galactic nuclei (LLAGN) presents unique challenges due to their comparable radiation output to their host galaxies and complex accretion dynamics. We introduce a novel module within the CIGALE framework specifically designed for SED fitting of LLAGN, incorporating both empirical relationships like $L_\mathrm{X}$--$L_\mathrm{12\mu m}$ and physically-based accretion models such as advection-dominated accretion flows (ADAFs) and truncated accretion disks. This allows for more accurate depiction of LLAGN central emissions. Using this module, we analyzed a set of 52 X-ray-detected local galaxies, primarily LINERs and Seyferts, and compared its performance to higher-luminosity AGN from the COSMOS and SDSS datasets. Our results show that the module adeptly estimates bolometric luminosities with high precision, despite significant galaxy contamination. It also introduces a versatile X-ray bolometric correction formula covering a vast range of luminosities. Further, our study explored the $\alpha_\mathrm{ox}$ index, which measures the UV to X-ray emission slope, showing that unlike quasars, LLAGN display either stable or only slightly varying $\alpha_\mathrm{ox}$ values, indicating differing accretion and photon production processes in the low luminosity regime. Additionally, we observed a significant drop of 1.4 dex in specific star formation rates when moving from whole galaxies to a central 9-arcsecond aperture in LLAGN, suggesting potential feedback mechanisms at play. Overall, our findings underscore the importance of a multiwavelength approach in AGN studies, highlighting distinct behaviors of LLAGN compared to quasars, thus enhancing our understanding of LLAGN and providing a framework for future comprehensive AGN population studies., Comment: Submitted to A&A

PDS 456 is the most luminous RQQ at z<0.3 and can be regarded as a local counterpart of the powerful QSOs shining at Cosmic Noon. It hosts a strong nuclear X-ray ultra-fast outflow, and a massive and clumpy CO(3-2) molecular outflow extending up to 5 kpc from the nucleus. We analyzed the first MUSE WFM and AO-NFM optical integral field spectroscopic observations of PDS456. The AO-NFM observations provide an unprecedented spatial resolution, reaching up to 280 pc. Our findings reveal a complex circumgalactic medium around PDS 456, extending up to a maximum projected size of ~46 kpc. This includes a reservoir of gas with a mass of ~1e7-1e8 Modot, along with eight companion galaxies, and a multi-phase outflow. WFM and NFM MUSE data reveal an outflow on a large scale (~12 kpc from the quasar) in [OIII], and on smaller scales (within 3 kpc) with higher resolution (about 280 pc) in Halpha, respectively. The [OIII] outflow mass rate is 2.3 +/- 0.2 Modot/yr which is significantly lower than those typically found in other luminous quasars. Remarkably, the Ha outflow shows a similar scale, morphology, and kinematics to the CO(3-2) molecular outflow, with the latter dominating in terms of kinetic energy and mass outflow rate by two and one orders of magnitude, respectively. Our results therefore indicate that mergers, powerful AGN activity, and feedback through AGN-driven winds will collectively contribute to shaping the host galaxy evolution of PDS 456, and likely, that of similar objects at the brightest end of the AGN luminosity function across all redshifts. Moreover, the finding that the momentum boost of the total outflow deviates from the expected energy-conserving expansion for large-scale outflows highlights the need of novel AGN-driven outflow models to comprehensively interpret these phenomena., Comment: 15 pages, 13 figures accepted for publication in Astronomy & Astrophysics

The detection of blue-shifted absorption lines likely associated with ionized Iron K-shell transitions in the X-ray spectra of many Active Galactic Nuclei (AGN) suggests the presence of a highly ionized gas outflowing with mildly relativistic velocities (0.03c-0.6c), named Ultra-Fast Outflow (UFO). Within the SUBWAYS project we characterized these winds starting from a sample of 22 radio-quiet quasars at 0.1 < z < 0.4, and compared the results with similar studies in the literature on samples of 42 local radio-quiet Seyfert galaxies and 14 high redshift radio-quiet quasars. The scope of our work is a statistical study of UFO parameters and incidence, considering key physical properties of the sources, e.g. supermassive black hole (SMBH) mass, bolometric luminosity, accretion rates and Spectral Energy Distribution, with the aim of gaining new insights into the UFO launching mechanisms. We find indications that highly luminous AGN with steeper X-ray/UV ratio, are more likely to host UFO. The presence of UFO is not significantly related to any other AGN property in our sample. These findings suggest that the UFO phenomenon may be transient. Focusing on AGN with UFO, other important results are: (1) faster UFO have larger ionization parameters and column densities; (2) X-ray radiation plays a more crucial role in driving highly ionized winds compared to UV; (3) the correlation between outflow velocity and luminosity is significantly flatter than what expected for radiatively driven winds; (4) more massive BH experience higher wind mass-losses, suppressing accretion of matter onto the BH; (5) the UFO launching radius is positively correlated with the Eddington ratio. Furthermore, our analysis suggest the involvement of multiple launching mechanisms, including radiation pressure and magneto-hydrodynamic processes, rather than pointing to a single, universally applicable mechanism., Comment: 52 pages, 154 figures and 7 tables. Accepted in Astronomy & Astrophysics, in press

Powerful outflows from active galactic nuclei (AGN) can significantly impact the gas reservoirs of their host galaxies. However, it is still unclear how these outflows can propagate from the very central regions of galaxies to their outskirts, and whether nuclear winds can be driven by and/or be responsible for drastic spectral transitions. In this work we test feedback propagation models on the case test of 2MASS 0918+2117 (2M0918), a z=0.149 X-ray variable AGN, which showed tentative evidence for nuclear ultra-fast outflows (UFOs) in a 2005 XMM-Newton observation. We also investigate whether UFOs can be related to the observed X-ray variability. We observed 2M0918 with XMM-Newton and NuSTAR in 2020 to confirm the presence and characterize the UFOs. We perform a kinematic analysis of the 2005 SDSS optical spectrum to reveal and measure the properties of galaxy-scale ionized outflows. Furthermore, we construct 20-year-long lightcurves of observed flux, line-of-sight column density, and intrinsic accretion rate from the spectra of the first 4 SRG/eROSITA all-sky surveys and archival observations from Chandra and XMM-Newton.We significantly detect UFOs with v$\sim$0.16c and galaxy-scale ionized outflows with velocities of $\sim$ 700 km/s. We also find that the drastic X-ray variability (factors >10) can be explained both in terms of variable obscuration and variable intrinsic luminosity.Comparing the energetics of the two outflow phases, 2M0918 is consistent with momentum-driven wind propagation. 2M0918 expands the sample of AGN with both UFOs and ionized gas winds from 5 to 6, and brings the sample of AGN hosting multiscale outflows to 19, contributing to a clearer picture of feedback physics. From the variations in accretion rate, column density, and ionization level of the obscurer, we propose a scenario that connects obscurers, an accretion enhancement, and the emergence of UFOs, Comment: 18 pages, 18 figures

The eROSITA instrument aboard the Spectrum Roentgen Gamma (SRG) satellite has performed its first all-sky survey between December 2019 and June 2020. This paper presents the resulting hard X-ray (2.3-5 keV) sample, the first created from an all-sky imaging survey in the 2-8 keV band, for sources within western galactic sky. The 5466 hard X-ray selected sources detected with eROSITA are presented and discussed. The Bayesian statistics-based code NWAY is used to identify the counterparts for the X-ray sources. These sources are classified based on their multiwavelength properties, and the literature is searched to identify spectroscopic redshifts, which further inform the source classification. A total of 2547 sources are found to have good-quality counterparts, and 111 of these are detected only in the hard band. Comparing with other hard X-ray selected surveys, the eROSITA hard sample covers a larger redshift range and probes dimmer sources, providing a complementary and expanded sample as compared to Swift-BAT. Examining the column density distribution of missed and detected eROSITA sources present in the follow-up catalog of Swift BAT 70 month sources, it is demonstrated that eROSITA can detect obscured sources with column densities $>10^{24}$ cm$^{-2}$, but that the completeness drops rapidly after $10^{23}$ cm$^{-2}$. A sample of hard-only sources, many of which are likely to be heavily obscured AGN, is also presented and discussed. X-ray spectral fitting reveals that these sources have extremely faint soft X-ray emission and their optical images suggest that they are found in more edge-on galaxies with lower b/a. The resulting X-ray catalog is demonstrated to be a powerful tool for understanding AGN, in particular heavily obscured AGN found in the hard-only sample., Comment: 20 pages, 21 figures. Submitted to A&A. Associated with eROSITA data release. Associated catalogs will be released after journal publication

During its calibration and performance verification phase, the eROSITA instrument aboard the SRG satellite performed a uniform wide--area X-ray survey of approximately 140 deg$^{2}$ in a region of the sky known as the eROSITA Final Equatorial Depth Survey (eFEDS). The primary aim of eFEDS is to demonstrate the scientific performance to be expected at the end of the 8-pass eROSITA all sky survey. This will provide the first focussed image of the whole sky in the hard X-ray ($>2$~keV) bandpass. The expected source population in this energy range is thus of great interest, particularly for AGN studies. We use the 2.3--5 keV selection presented by Brunner et al. (2022) to construct a sample of 246 point-like hard X-ray sources for further study and characterization. These are classified as either extragalactic ($\sim 90$~\%) or Galactic ($\sim 10$~\%), with the former consisting overwhelmingly of AGN and the latter active stars. We concentrate our further analysis on the extragalactic/AGN sample, describing their X-ray and multiwavelength properties and comparing them to the eFEDS main AGN sample selected in the softer 0.2-2.3 keV band. The eROSITA hard band selects a subsample of sources that is a factor $>10$ brighter than the eFEDS main sample. The AGN within the hard population reach up to $z=3.2$ but on the whole are relatively nearby, with median $z$=0.34 compared to $z$=0.94 for the main sample. The hard survey probes typical luminosities in the range $\log L_{\rm X} = 43-46$. X-ray spectral analysis shows significant intrinsic absorption (with $\log N_{\rm H}>21$) in $\sim 20$~\% of the sources, with a hard X-ray power law continuum with mean $<\Gamma>=1.83\pm0.04$, typical of AGN, but slightly harder than the soft-selected eROSITA sample. (abridged), Comment: Submitted to A&A as part of eROSITA Data Release 1

The eROSITA telescope array aboard the Spektrum Roentgen Gamma (SRG) satellite began surveying the sky in December 2019, with the aim of producing all-sky X-ray source lists and sky maps of an unprecedented depth. Here we present catalogues of both point-like and extended sources using the data acquired in the first six months of survey operations (eRASS1; completed June 2020) over the half sky whose proprietary data rights lie with the German eROSITA Consortium. We describe the observation process, the data analysis pipelines, and the characteristics of the X-ray sources. With nearly 930000 entries detected in the most sensitive 0.2-2.3 keV energy range, the eRASS1 main catalogue presented here increases the number of known X-ray sources in the published literature by more than 60%, and provides a comprehensive inventory of all classes of X-ray celestial objects, covering a wide range of physical processes. A smaller catalogue of 5466 sources detected in the less sensitive but harder 2.3-5 keV band is the result of the first true imaging survey of the entire sky above 2 keV. We show that the number counts of X-ray sources in eRASS1 are consistent with those derived over narrower fields by past X-ray surveys of a similar depth, and we explore the number counts variation as a function of the location in the sky. Adopting a uniform all-sky flux limit (at 50% completeness) of F_{0.5-2 keV} > 5 \times 10^{-14}$ erg\,s$^{-1}$\,cm$^{-2}$, we estimate that the eROSITA all-sky survey resolves into individual sources about 20% of the cosmic X-ray background in the 1-2 keV range. The catalogues presented here form part of the first data release (DR1) of the SRG/eROSITA all-sky survey. Beyond the X-ray catalogues, DR1 contains all detected and calibrated event files, source products (light curves and spectra), and all-sky maps. Illustrative examples of these are provided., Comment: 39 pages, 23 figures. Accepted for publication in A&A. Accompanying eROSITA-DE Data Release 1

We present MEGARA (Multi-Espectr\'ografo en GTC de Alta Resoluci\'on para Astronom\'ia) Integral Field Unit (IFU) observations of 5 local type-2 quasars (QSO2s, z $\sim 0.1$) from the Quasar Feedback (QSOFEED) sample. These active galactic nuclei (AGN) have bolometric luminosities of 10$^{45.5-46}$ erg/s and stellar masses of $\sim$10$^{11}$ M$_{\odot}$. We explore the kinematics of the ionized gas through the [O~III]$\lambda$5007 $\r{A}$ emission line. The nuclear spectra of the 5 QSO2s, extracted in a circular aperture of $\sim$ 1.2" ($\sim$ 2.2 kpc) in diameter, show signatures of high velocity winds in the form of broad (full width at half maximum; 1300$\leq$FWHM$\leq$2240 km/s and blueshifted components. We find that 4 out of the 5 QSO2s present outflows that we can resolve with our seeing-limited data, and they have radii ranging from 3.1 to 12.6 kpc. In the case of the two QSO2s with extended radio emission, we find that it is well-aligned with the outflows, suggesting that low-power jets might be compressing and accelerating the ionized gas in these radio-quiet QSO2s. In the four QSO2s with spatially resolved outflows, we measure ionized mass outflow rates of 3.3-6.5 Msun/yr when we use [S~II]-based densities, and of 0.7-1.6 Msun/yr when trans-auroral line-based densities are considered instead. We compare them with the corresponding molecular mass outflow rates (8 - 16 Msun/yr), derived from CO(2-1) ALMA observations at 0.2" resolution. Both phases show lower outflow mass rates than those expected from observational scaling relations where uniform assumptions on the outflow properties were adopted. This might be indicating that the AGN luminosity is not the only driver of massive outflows and/or that these relations need to be re-scaled using accurate outflow properties. We do not find a significant impact of the outflows on the global star formation rates., Comment: 24 pages, 13 figures and 4 tables. Accepted for publication in A&A; A&A 665, A55 (2023); doi: 10.1051/0004-6361/202347715

The existence of luminous quasars (QSO) at the Epoch of Reionization (EoR; i.e. z>6) powered by supermassive black holes (SMBH) with masses $\gtrsim10^9~M_\odot$ challenges models of early SMBH formation. To shed light on the nature of these sources we started a multiwavelength programme based on a sample of 18 HYPerluminous quasars at the Epoch of ReionizatION (HYPERION). These are the luminous QSOs whose SMBH must have had the fastest mass growth during the Universe first Gyr. In this paper we present the HYPERION sample and report on the first of the 3 years planned observations of the 2.4 Ms XMM-Newton Multi-Year Heritage program on which HYPERION is based. The goal of this program is to accurately characterize the X-ray nuclear properties of QSOs at the EoR. Through a joint X-ray spectral analysis of 10 sources, in the rest-frame $\sim2-50$ keV range, we report a steep average photon index ($\Gamma\sim2.4\pm0.1$). Absorption is not required. The average $\Gamma$ is inconsistent at $\geq4\sigma$ level with the canonical 1.8-2 value measured in QSO at z<6. This spectral slope is also much steeper than that reported in lower-z QSOs with similar luminosity or accretion rate, thus suggesting a genuine redshift evolution. Alternatively, we can interpret this result as the presence of an unusually low-energy cutoff $E_{cut}\sim20$ keV on a standard $\Gamma=1.9$ power-law. We also report on mild indications that HYPERION QSOs show higher soft X-ray emission at 2 keV compared to the UV one at 2500A than expected by lower-z luminous AGN. We speculate that a redshift-dependent coupling between the corona and accretion disc or intrinsically different coronal properties may account for the steep spectral slopes, especially in the presence of powerful winds. The reported slopes, if confirmed at lower luminosities, may have an important impact on future X-ray AGN studies in the early Universe., Comment: 21 pages (including appendix), 12 figures, 4 tables. Accepted for pubblication in A&A

The paper presents the analysis of optically selected GAMA groups and clusters in the SRG/eROSITA X-ray map of eFEDS (eROSITA Final Equatorial Depth Survey), in the halo mass range $10^{13}-5{\times}10^{14}$ $M_{\odot}$ and at $z < 0.2$. All X-ray detections have a clear GAMA counterpart, but most of the GAMA groups in the halo mass range $10^{13}-10^{14}$ $M_{\odot}$ remain undetected. We compare the X-ray surface brightness profiles of the eROSITA detected groups with the mean stacked profile of the undetected low-mass halos at fixed halo mass. Overall, we find that the undetected groups exhibit less concentrated X-ray surface brightness, dark matter, and galaxy distributions with respect to the X-ray detected halos. The mean gas mass fraction profiles are consistent in the two samples within 1.5$\sigma$, indicating that the gas follows the dark matter profile. The low mass concentration and the magnitude gap indicate that these systems are young. They reside with a higher probability in filaments while X-ray detected groups favor the nodes of the Cosmic Web. Because of the lower central emission, the undetected systems tend to be X-ray under-luminous at fixed halo mass and to lie below the $L_X-M_{halo}$ relation. Interestingly, the X-ray detected systems inhabiting the nodes scatter the less around the relation, while those in filaments tend to lie below it. We do not observe any strong relationship between the system X-ray appearance and the AGN activity. We cannot exclude the role of the past AGN feedback in affecting the gas distribution over the halo lifetime. However, the data suggests that the observed differences might be related to the halo assembly bias., Comment: 16 pages, 14 figures, Accepted for publication in MNRAS

Studies indicate strong evidence of a scaling relation in the local Universe between the supermassive black hole mass ($M_\rm{BH}$) and the stellar mass of their host galaxies ($M_\star$). They even show similar histories across cosmic times of their differential terms: star formation rate (SFR) and black hole accretion rate (BHAR). However, a clear picture of this coevolution is far from being understood. We select an X-ray sample of active galactic nuclei (AGN) up to $z=2.5$ in the miniJPAS footprint. Their X-ray to infrared spectral energy distributions (SEDs) have been modeled with CIGALE, constraining the emission to 68 bands. For a final sample of 308 galaxies, we derive their physical properties (e.g., $M_\star$, $\rm{SFR}$, $\rm{SFH}$, and $L_\rm{AGN}$). We also fit their optical spectra for a subsample of 113 sources to estimate the $M_\rm{BH}$. We calculate the BHAR depending on two radiative efficiency regimes. We find that the Eddington ratios ($\lambda$) and its popular proxy ($L_\rm{X}$/$M_\star$) have 0.6 dex of difference, and a KS-test indicates that they come from different distributions. Our sources exhibit a considerable scatter on the $M_\rm{BH}$-$M_\star$ relation, which can explain the difference between $\lambda$ and its proxy. We also model three evolution scenarios to recover the integral properties at $z=0$. Using the SFR and BHAR, we show a notable diminution in the scattering between $M_\rm{BH}$-$M_\star$. For the last scenario, we consider the SFH and a simple energy budget for the AGN accretion, obtaining a relation similar to the local Universe. Our study covers $\sim 1$ deg$^2$ in the sky and is sensitive to biases in luminosity. Nevertheless, we show that, for bright sources, the link between SFR and BHAR, and their decoupling based on an energy limit is the key that leads to the local $M_\rm{BH}$-$M_\star$ scaling relation., Comment: Accepted for publication in A&A

We present a new X-ray spectroscopic study of $22$ luminous ($2\times10^{45}\lesssim L_{\rm bol}\rm /erg\,s^{-1} \lesssim 2\times10^{46}$) active galactic nuclei (AGNs) at intermediate-redshift ($0.1 \lesssim z \lesssim 0.4$), as part of the SUpermassive Black hole Winds in the x-rAYS (SUBWAYS) sample, mostly composed of quasars (QSOs) and type\,1 AGN. Here, 17 targets were observed with \textit{XMM-Newton} between 2019--2020 and the remaining 5 are from previous observations. The aim of this large campaign ($1.45\,\rm Ms$ duration) is to characterise the various manifestations of winds in the X-rays driven from supermassive black holes in AGN. In this paper we focus on the search and characterization of ultra-fast outflows (UFOs), which are typically detected through blueshifted absorption troughs in the Fe\,K band ($E>7\,\rm keV$). By following Monte Carlo procedures, we confirm the detection of absorption lines corresponding to highly ionised iron (e.g., Fe\,\textsc{xxv}\,H$\alpha$, Fe\,\textsc{xxvi}\,Ly$\alpha$) in 7/22 sources at the $\gtrsim95\%$ confidence level (for each individual line). The global combined probability of such absorption features in the sample is $>99.9\%$. The SUBWAYS campaign extends at higher luminosity and redshifts than previous local studies on Seyferts, obtained using \xmm and \suzaku observations. We find a UFO detection fraction of $\sim30\%$ on the total sample that is in agreement with the previous findings. This work independently provides further support for the existence of highly-ionised matter propagating at mildly relativistic speed ($\gtrsim0.1c$) in a considerable fraction of AGN over a broad range of luminosities, which is expected to play a key role in the self-regulated AGN feeding-feedback cycle, as also supported by hydrodynamical multiphase simulations., Comment: 30 pages, 10 figures, accepted for publication in A&A

We present a UV spectroscopic study of ionized outflows in 21 active galactic nuclei (AGN), observed with the HST. The targets of the SUBWAYS sample were selected with the aim to probe the parameter space of the underexplored AGN between the local Seyfert galaxies and the luminous quasars at high redshifts. Our targets, spanning redshifts of 0.1-0.4 and bolometric luminosities (L_bol) of 10^45-10^46 erg/s, have been observed with a large multi-wavelength campaign. Here, we model the UV spectra and look for different types of AGN outflows. We find that 60% of our targets show a presence of outflowing H I absorption, while 40% exhibit ionized outflows seen as absorption by either C IV, N V, or O VI. This is comparable to the occurrence of ionized outflows seen in the local Seyfert galaxies. All UV absorption lines in the sample are relatively narrow, with outflow velocities reaching up to -3300 km/s. We did not detect any UV counterparts to the X-ray ultra-fast outflows (UFOs), most likely due to their being too highly ionized. However, all SUBWAYS targets with an X-ray UFO demonstrate the presence of UV outflows at lower velocities. We find significant correlations between the column density (N) of the UV ions and L_bol of the AGN, with N of H I decreasing with L_bol, while N of O VI is increasing with L_bol. This is likely to be a photoionization effect, where toward higher AGN luminosities, the wind becomes more ionized, resulting in less absorption by neutral or low-ionization ions and more absorption by high-ionization ions. In addition, we find that N of the UV ions decreases as their outflow velocity increases. This may be explained by a mechanical power that is evacuating the UV-absorbing medium. Our observed relations are consistent with multiphase AGN feeding and feedback simulations indicating that a combination of both radiative and mechanical processes are in play., Comment: Accepted for publication in Astronomy & Astrophysics (A&A), 22 pages, 12 figures

Water vapor (H$_{2}$O) is one of the brightest molecular emitters after carbon monoxide (CO) in galaxies with high infrared (IR) luminosity, and allows us to investigate the warm dense phase of the interstellar medium (ISM) where star formation occurs. However, due to the complexity of its radiative spectrum, H$_{2}$O is not frequently exploited as an ISM tracer in distant galaxies. Therefore, H$_{2}$O studies of the warm and dense gas at high-$z$ remains largely unexplored. In this work we present observations conducted with the Northern Extended Millimeter Array (NOEMA) toward three $z>6$ IR-bright quasars J2310+1855, J1148+5251, and J0439+1634 targeted in their multiple para-/ortho-H$_{2}$O transitions ($3_{12}-3_{03}$, $1_{11}-0_{00}$, $2_{20}-2_{11}$, and $4_{22}-4_{13}$), as well as their far-IR (FIR) dust continuum. By combining our data with previous measurements from the literature we estimate dust masses and temperatures, continuum optical depths, IR luminosities, and the star-formation rates from the FIR continuum. We model the H$_{2}$O lines using the MOLPOP-CEP radiative transfer code and find that water vapor lines in our quasar host galaxies are primarily excited in warm dense (gas kinetic temperature and density of $T_{\rm kin} = 50\,{\rm K}$, $n_{\rm H_{2}}\sim 10^{4.5}-10^{5}\,{\rm cm^{-3}}$) molecular medium with water vapor column density of $N_{\rm H_{2}O}\sim 2\times10^{17}-3\times10^{18}\,{\rm cm^{-3}}$. High-$J$ H$_{2}$O lines are mainly radiatively pumped by the intense optically-thin far-IR radiation field associated with a warm dust component with temperatures of $T_{\rm dust}\sim 80-190\,{\rm K}$ that account for $<5-10\%$ of the total dust mass. Our results are in agreement with expectations based on the H$_{2}$O spectral line energy distribution of local and high-$z$ ultra-luminous IR galaxies and AGN. [abridged], Comment: Accepted for publication in A&A: 18 pages, 9 figures

We present a new X-Ray Accretion Disk-wind Emulator (\textsc{xrade}) based on the 2.5D Monte Carlo radiative transfer code which provides a physically-motivated, self-consistent treatment of both absorption and emission from a disk-wind by computing the local ionization state and velocity field within the flow. \textsc{xrade} is then implemented through a process that combines X-ray tracing with supervised machine learning. We develop a novel emulation method consisting in training, validating, and testing the simulated disk-wind spectra into a purposely built artificial neural network. The trained emulator can generate a single synthetic spectrum for a particular parameter set in a fraction of a second, in contrast to the few hours required by a standard Monte Carlo radiative transfer pipeline. The emulator does not suffer from interpolation issues with multi-dimensional spaces that are typically faced by traditional X-ray fitting packages such as \textsc{xspec}. \textsc{xrade} will be suitable to a wide number of sources across the black-hole mass, ionizing luminosity, and accretion rate scales. As an example, we demonstrate the applicability of \textsc{xrade} to the physical interpretation of the X-ray spectra of the bright quasar PDS 456, which hosts the best-established accretion-disk wind observed to date. We anticipate that our emulation method will be an indispensable tool for the development of high-resolution theoretical models, with the necessary flexibility to be optimized for the next generation micro-calorimeters on board future missions, like \textit{XRISM/resolve} and \textit{Athena/X-IFU}. This tool can also be implemented across a wide variety of X-ray spectral models and beyond., Comment: 20 pages, 17 figures, MNRAS accepted for publication

We analyze Near-Infrared Integral Field Spectrograph (NIFS) observations of the type-2 quasar (QSO2) SDSS J094521.33+173753.2 to investigate its warm molecular and ionized gas kinematics. This QSO2 has a bolometric luminosity of 10$^{45.7}$ erg s$^{-1}$ and a redshift of z = 0.128. The K-band spectra provided by NIFS cover a range of 1.99-2.40 $\mu$m where low-ionization (Pa$\alpha$ and Br$\delta$), high ionization ([S XI]$\lambda$1.920 $\mu$m and [Si~VI]$\lambda$1.963 $\mu$m) and warm molecular lines (from H$_2$ 1-0S(5) to 1-0S(1)) are detected, allowing us to study the multi-phase gas kinematics. Our analysis reveals gas in ordinary rotation in all the emission lines detected and also outflowing gas in the case of the low- and high-ionization emission lines. In the case of the nuclear spectrum, which corresponds to a circular aperture of 0.3\arcsec~(686 pc) in diameter, the warm molecular lines can be characterized using a single Gaussian component of full width at half maximum (FWHM)= 350-400 km s$^{-1}$, while Pa$\alpha$, Br$\delta$, and [Si~VI] are best fitted with two blue-shifted Gaussian components of FWHM$\sim$800 and 1700 km s$^{-1}$, in addition to a narrow component of $\sim$300 km s$^{-1}$. We interpret the blue-shifted broad components as outflowing gas, which reaches the highest velocities, of up to $-$840 km s$^{-1}$, in the south-east direction (PA$\sim$125$^{\circ}$), extending up to a distance of $\sim$3.4 kpc from the nucleus. The ionized outflow has a maximum mass outflow rate of $\dot{\text{{M}}}_{\text{{out, max}}}$=42-51 M$_\odot$ yr$^{-1}$, and its kinetic power represents 0.1$\%$ of the quasar bolometric luminosity., Comment: 24 pages, 20 figures and 5 tables. Accepted for publication in A&A

We combine results from deep ALMA observations of massive ($M_*>10^{10}\;M_{\odot}$) galaxies at different redshifts to show that the column density of their inter stellar medium (ISM) rapidly increases towards early cosmic epochs. Our analysis includes objects from the ASPECS and ALPINE large programs, as well as individual observations of $z\sim 6$ QSO hosts. When accounting for non-detections and correcting for selection effects, we find that the median surface density of the ISM of the massive galaxy population evolves as $\sim(1+z)^{3.3}$. This means that the ISM column density towards the nucleus of a $z>3$ galaxy is typically $>100$ times larger than locally, and it may reach values as high as Compton-thick at $z\gtrsim6$. Remarkably, the median ISM column density is of the same order of what is measured from X-ray observations of large AGN samples already at $z\gtrsim2$. We develop a simple analytic model for the spatial distribution of ISM clouds within galaxies, and estimate the total covering factor towards active nuclei when obscuration by ISM clouds on the host scale is added to that of pc-scale circumnuclear material (the so-called 'torus'). The model includes clouds with a distribution of sizes, masses, and surface densities, and also allows for an evolution of the characteristic cloud surface density with redshift, $\Sigma_{c,*}\propto(1+z)^\gamma$. We show that, for $\gamma=2$, such a model successfully reproduces the increase of the obscured AGN fraction with redshift that is commonly observed in deep X-ray surveys, both when different absorption thresholds and AGN luminosities are considered. Our results suggest that 80-90\% of supermassive black holes in the early Universe ($z>6-8$) are hidden to our view, primarily by the ISM in their hosts. [abridged], Comment: 22 pages, 15 figures, accepted for publication in A&A

Theoretical models of galaxy-AGN co-evolution ascribe an important role for the feedback process to a short, luminous, obscured, and dust-enshrouded phase during which the accretion rate of the SMBH is expected to be at its maximum and the associated AGN-driven winds are also predicted to be maximally developed. To test this scenario, we have isolated a text-book candidate from the eROSITA Final Equatorial-Depth Survey (eFEDS) obtained within the Performance and Verification program of the eROSITA telescope on board Spectrum Roentgen Gamma. From an initial catalog of 246 hard X-ray selected sources matched with the photometric and spectroscopic information available within the eROSITA and Hyper Suprime-Cam consortia, three candidates Quasars in the feedback phase have been isolated applying the diagnostic proposed in Brusa et al. (2015). Only one source (eFEDSU J091157.5+014327) has a spectrum already available (from SDSS-DR16, z=0.603) and it unambiguously shows the presence of a broad component (FWHM~1650 km/s) in the [OIII]5007 line. The associated observed L_[OIII] is ~2.6x10^{42} erg/s, one to two orders of magnitude larger than that observed in local Seyferts and comparable to those observed in a sample of z~0.5 Type 1 Quasars. From the multiwavelength data available we derive an Eddington Ratio (L_bol/L_Edd) of ~0.25, and a bolometric correction in the hard X-ray of k_bol~10, lower than those observed for objects at similar bolometric luminosity. The presence of an outflow, the high X-ray luminosity and moderate X-ray obscuration (L_X~10^44.8 erg/s, N_H~2.7x10^22 cm^-2) and the red optical color, all match the prediction of quasars in the feedback phase from merger driven models. Forecasting to the full eROSITA all-sky survey with its spectroscopic follow-up, we predict that by the end of 2024 we will have a sample of few hundreds such objects at z=0.5-2., Comment: 11 pages, 6 figures, 2 tables. Resubmitted to A&A for the Special Issue: "The Early Data Release of eROSITA and Mikhail Pavlinsky ART-XC on the SRG Mission"

In November 2019, eROSITA on board of SRG observatory started to map the entire sky in X-rays. After the 4-year survey program, it will reach flux limits about 25 times deeper than ROSAT. During the SRG Performance Verification phase, eROSITA observed a contiguous 140 deg$^2$ area of the sky down to the final depth of the eROSITA all-sky survey ("eROSITA Final Equatorial-Depth Survey": eFEDS), with the goal of obtaining a census of the X-ray emitting populations (stars, compact objects, galaxies, clusters of galaxies, AGN) that will be discovered over the entire sky. This paper presents the identification of the counterparts to the point-sources detected in eFEDS in the Main and Hard samples described in Brunner et al 2021, and their multi-wavelength properties, including redshift. For the identification of the counterparts we combined the results from two independent methods (NWAY and ASTROMATCH), trained on the multi-wavelength properties of a sample of 23k XMM-Newton sources detected in the DESI Legacy Imaging Survey DR8. Then spectroscopic redshifts and photometry from ancillary surveys are collated for the computation of photometric redshifts. The eFEDS sources with a reliable counterparts are 24774/27369 (90.5\%) in the Main sample and 231/246 (93.9\%) in the Hard sample, including 2514 (3) sources for which a second counterpart is equally likely. [abridged] This paper is accompanying the eROSITA early data release of all the observations performed during the performance and verification phase. Together with the catalogs of primary and secondary counterparts to the Main and Hard samples of the eFEDS survey this paper releases their multi-wavelength properties and redshifts., Comment: Accepted by A&A for the Special Issue: The Early Data Release of eROSITA and Mikhail Pavlinsky ART-XC on the SRG Mission. NOTE: Paper completely revised wrt the first version submitted to arXiv. Revised catalogs available via https://erosita.mpe.mpg.de/edr/eROSITAObservations/Catalogues/

Context. The eROSITA X-ray telescope on board the Spectrum-Roentgen-Gamma (SRG) observatory combines a large field of view and collecting area in the energy range $\sim$0.2 to $\sim$8.0 keV with the capability to perform uniform scanning observations of large sky areas. Aims. SRG/eROSITA performed scanning observations of the $\sim$140 square degrees eROSITA Final Equatorial Depth Survey (eFEDS) field as part of its performance verification phase. The observing time was chosen to slightly exceed the depth of equatorial fields after the completion of the eROSITA all-sky survey. We present a catalog of detected X-ray sources in the eFEDS field providing source positions and extent information, as well as fluxes in multiple energy bands and document a suite of tools and procedures developed for eROSITA data processing and analysis, validated and optimized by the eFEDS work. Methods. A multi-stage source detection procedure was optimized and calibrated by performing realistic simulations of the eROSITA eFEDS observations. We cross-matched the eROSITA eFEDS source catalog with previous XMM-ATLAS observations, confirming excellent agreement of the eROSITA and XMM-ATLAS source fluxes. Result. We present a primary catalog of 27910 X-ray sources, including 542 with significant spatial extent, detected in the 0.2-2.3 keV energy range with detection likelihoods $\ge 6$, corresponding to a (point source) flux limit of $\approx 6.5 \times 10^{-15}$ erg/cm$^2$/s in the 0.5-2.0 keV energy band (80% completeness). A supplementary catalog contains 4774 low-significance source candidates with detection likelihoods between 5 and 6. In addition, a hard band sample of 246 sources detected in the energy range 2.3-5.0 keV above a detection likelihood of 10 is provided. The dedicated data analysis software package, calibration database, and calibrated data products are described in an appendix., Comment: 25 pages, 15 figures, accepted, A&A Special Issue: The Early Data Release of eROSITA and Mikhail Pavlinsky ART-XC on the SRG Mission

Outflows driven by active galactic nuclei (AGN) are expected to have a significant impact on the host galaxy evolution, but it is still debated how they are accelerated and propagate on galaxy-wide scales. This work addresses these questions by studying the link between X-ray, nuclear ultra-fast outflows (UFOs) and extended ionised outflows, for the first time in two quasars close to the peak of AGN activity ($z\sim2$), where AGN feedback is expected to be more effective. As targets, we selected two multiple-lensed quasars at $z\sim1.5$, HS 0810+2554 and SDSS J1353+1138, known to host UFOs and observed with the near-IR integral field spectrometer SINFONI at the VLT. We performed a kinematical analysis of the [O III]$\lambda$5007 optical emission line, in order to trace the presence of ionised outflows. We detected spatially resolved ionised outflows in both galaxies, extended more than 8 kpc and moving up to $v>2000$ km/s. We derived mass outflow rates of $\sim$12 M$_{sun}$/yr and $\sim$2 M$_{sun}$/yr for HS 0810+2554 and SDSS J1353+1138. Comparing with the co-hosted UFO energetics, the ionised outflow energetics in HS 0810+2554 is broadly consistent with a momentum-driven regime of wind propagation, while in SDSS J1353+1138 it differs by a factor of $\sim$100 from theoretical predictions, requiring either a massive molecular outflow or a high variability of the AGN activity to account for such a discrepancy. By additionally considering our results with those from the small sample of well-studied objects (all local but one), with both UFO and extended (ionised/atomic/molecular) outflow detections, we found that in 10 out of 12 galaxies the large-scale outflow energetics is consistent with the theoretical predictions of either a momentum- or an energy-driven scenario. This suggests that such models explain relatively well the acceleration mechanism of AGN-driven winds on large scales., Comment: 19 pages, 9 figures. Accepted for publication in A&A

VANDELS is an ESO Public Spectroscopic Survey designed to build a sample of high signal to noise, medium resolution spectra of galaxies at redshift between 1 and 6.5. Here we present the final Public Data Release of the VANDELS Survey, comprising 2087 redshift measurements. We give a detailed description of sample selection, observations and data reduction procedures. The final catalogue reaches a target selection completeness of 40% at iAB = 25. The high Signal to Noise ratio of the spectra (above 7 in 80% of the spectra) and the dispersion of 2.5{\AA} allowed us to measure redshifts with high precision, the redshift measurement success rate reaching almost 100%. Together with the redshift catalogue and the reduced spectra, we also provide optical mid-IR photometry and physical parameters derived through SED fitting. The observed galaxy sample comprises both passive and star forming galaxies covering a stellar mass range 8.3< Log(M*/Msolar)<11.7. All catalogues and spectra are accessible through the survey database (http://vandels.inaf.it) where all information can be queried interactively, and via the ESO Archive (https://www.eso.org/qi/)., Comment: 16 pages, 13 figures, accepted for publication in Astronomy & Astrophysics

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) was designed to conduct a blind search for AGN-driven outflows on X-ray selected AGN at redshift z$\sim$2 with high ($\sim$2 kpc) spatial resolution, and correlate them to the properties of the host galaxy and central black hole. The main aims of this paper are: a) to derive reliable estimates for the BH mass and accretion rates for the Type-1 AGN in this survey; b) to characterize the properties of the AGN driven winds in the BLR. We analyzed rest-frame optical and UV spectra of 21 Type-1 AGN. We found that the BH masses estimated from H$\alpha$ and H$\beta$ lines are in agreement. We estimate BH masses in the range Log(M$\rm_{BH}/M_{\odot}$)=8.4-10.8 and Eddington ratios $\rm\lambda_{Edd}$ =0.04-1.3. We confirm that the CIV line width does not correlate with the Balmer lines and the peak of the line profile is blue-shifted with respect to the [OIII]-based systemic redshift. These findings support the idea that the CIV line is tracing outflowing gas in the BLR, with velocities up to $\sim$4700 km/s. We confirm the strong dependence of the BLR wind velocity with the UV-to-Xray continuum slope, L$\rm_{Bol}$ and $\rm\lambda_{Edd}$. We inferred BLR mass outflow rates in the range 0.005-3 M$_{\odot}$/yr, showing a correlation with the bolometric luminosity consistent with that observed for ionized winds in the NLR and X-ray winds detected in local AGN, and kinetic power $\sim$10$^{[-7:-4]}\times$ L$\rm_{Bol}$. Finally, we found an anti-correlation between the equivalent width of the [OIII] line with respect to the CIV shift, and a positive correlation with [OIII] outflow velocity. These findings, for the first time in an unbiased sample of AGN at z$\sim$2, support a scenario where BLR winds are connected to galaxy scale detected outflows, and are capable of affecting the gas in the NLR located at kpc scale., Comment: Accepted for publication in A&A

We report on ALMA observations of D1, a system at z~6.15 with stellar mass M_* ~ 10^7 M_sun containing globular cluster (GC) precursors, strongly magnified by the galaxy cluster MACS J0416.1-2403. Since the discovery of GC progenitors at high redshift, ours is the first attempt to probe directly the physical properties of their neutral gas through infrared observations. A careful analysis of our dataset, performed with a suitable procedure designed to identify faint narrow lines and which can test various possible values for the unknown linewidth value, allowed us to identify a 4-sigma tentative detection of [CII] emission with intrinsic luminosity L_CII=(2.9 +/- 1.4) 10^6 L_sun, one of the lowest values ever detected at high redshift. This study offers a first insight on previously uncharted regions of the L_CII-SFR relation. Despite large uncertainties affecting our measure of the star formation rate, if taken at face value our estimate lies more than 1 dex below the values observed in local and high redshift systems. Our weak detection indicates a deficiency of [CII] emission, possibly ascribed to various explanations, such as a low-density gas and/or a strong radiation field caused by intense stellar feedback, and a low metal content. From the non-detection in the continuum we derive constraints on the dust mass, with 3-sigma upper limit values as low as a few 10^4 M_sun, consistent with the values measured in local metal-poor galaxies., Comment: MNRAS, accepted. 15 pages, 5 figures

eROSITA (extended ROentgen Survey with an Imaging Telescope Array) is the primary instrument on the Spectrum-Roentgen-Gamma (SRG) mission, which was successfully launched on July 13, 2019, from the Baikonour cosmodrome. After the commissioning of the instrument and a subsequent calibration and performance verification phase, eROSITA started a survey of the entire sky on December 13, 2019. By the end of 2023, eight complete scans of the celestial sphere will have been performed, each lasting six months. At the end of this program, the eROSITA all-sky survey in the soft X-ray band (0.2--2.3\,keV) will be about 25 times more sensitive than the ROSAT All-Sky Survey, while in the hard band (2.3--8\,keV) it will provide the first ever true imaging survey of the sky. The eROSITA design driving science is the detection of large samples of galaxy clusters up to redshifts $z>1$ in order to study the large-scale structure of the universe and test cosmological models including Dark Energy. In addition, eROSITA is expected to yield a sample of a few million AGNs, including obscured objects, revolutionizing our view of the evolution of supermassive black holes. The survey will also provide new insights into a wide range of astrophysical phenomena, including X-ray binaries, active stars, and diffuse emission within the Galaxy. Results from early observations, some of which are presented here, confirm that the performance of the instrument is able to fulfil its scientific promise. With this paper, we aim to give a concise description of the instrument, its performance as measured on ground, its operation in space, and also the first results from in-orbit measurements., Comment: 16 pages, 19 figures, accepted by Astronomy & Astrophysics

We use MUSE adaptive optics (AO) data in Narrow Field Mode to study the properties of the ionised gas in MR 2251-178 and PG 1126-041, two nearby (z~0.06) bright quasars hosting sub-pc scale Ultra Fast Outflows (UFOs) detected in the X-ray band. We decompose the optical emission from diffuse gas into a low- and a high-velocity components. The former is characterised by a clean, regular velocity field and a low (~80 km/s) velocity dispersion. It traces regularly rotating gas in PG 1126-041, while in MR 2251-178 it is possibly associated to tidal debris from a recent merger or flyby. The other component is found to be extended up to a few kpc from the nuclei, and shows a high (~800 km/s) velocity dispersion and a blue-shifted mean velocity, as expected from AGN-driven outflows. We estimate mass outflow rates up to a few Mo/yr and kinetic efficiencies between 0.1-0.4 per cent, in line with those of galaxies hosting AGNs of similar luminosity. The momentum rates of these ionised outflows are comparable to those measured for the UFOs at sub-pc scales, consistent with a momentum-driven wind propagation. Pure energy-driven winds are excluded unless about 100x additional momentum is locked in massive molecular winds. By comparing the outflow properties of our sources with those of a small sample of well-studied QSOs hosting UFOs from the literature, we find that winds seem to systematically lie either in a momentum-driven or in an energy-driven regime, indicating that these two theoretical models bracket very well the physics of AGN-driven winds., Comment: 18 pages, 10 figures. Key figures are 4 and 9. Accepted for publication in A&A

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., Comment: Accepted for publication in A&A

Based on ALMA Band 3 observations of the CO(2-1) line transition, we report the discovery of three new gas-rich (M_H2 ~ 1.5-4.8 x 10^10 M_sun, SFRs in the range ~5-100 M_sun/yr) galaxies in an overdense region at z=1.7, that already contains eight spectroscopically confirmed members. This leads to a total of 11 confirmed overdensity members, within a projected distance of ~ 1.15 Mpc and in a redshift range of Dz = 0.012. Under simple assumptions, we estimate that the system has a total mass of >= 3-6 x 10^13 M_sun, and show that it will likely evolve into a >~ 10^14 M_sun cluster at z = 0. The overdensity includes a powerful Compton-thick Fanaroff-Riley type II (FRII) radio-galaxy, around which we discovered a large molecular gas reservoir (M_H2 ~ 2 x 10^11 M_sun). We fitted the FRII resolved CO emission with a 2-D Gaussian model with major (minor) axis of ~ 27 (~ 17) kpc, that is a factor of ~3 larger than the optical rest-frame emission. Under the assumption of a simple edge-on disk morphology, we find that the galaxy interstellar medium produces a column density towards the nucleus of ~ 5.5 x 10^23 cm^-2. Such a dense ISM may then contribute significantly to the total nuclear obscuration measured in the X-rays (N_(H,X) ~ 1.5 x 10^24 cm^-2) in addition to a small, pc-scale absorber around the central engine. The velocity map of this source unveils a rotational motion of the gas that is perpendicular to the radio-jets. The FRII is located at the center of the projected spatial distribution of the structure members, and its velocity offset from the peak of the redshift distribution is well within the structure's velocity dispersion. All this, coupled with the large amount of gas around the FRII, its stellar mass of ~ 3 x 10^11 M_sun, SFR of ~ 200-600 M_sun/yr, and powerful radio-to-X-ray emission, suggests that this source is the likely progenitor of the future brightest cluster galaxy., Comment: 7 pages, 4 figures, 2 tables. Accepted for publication as a Letter 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%., Comment: Accepted for publication in A&A, 34 pages, 14 figures and 5 tables

We study the co-evolution between the black hole accretion rate (BHAR) and the star formation rate (SFR) in different galaxy life phases: main sequence star-forming galaxies, quiescent and starburst galaxies at different cosmic epochs. We take advantage of the X-ray data from the Chandra COSMOS-Legacy survey and of the extensive multiwavelength ancillary observations in the COSMOS field presented in the COSMOS2015 catalog. We perform an X-ray stacking analysis and combine it with detected sources, in a broad redshift interval ($0.1

We present X-ray spectra of 185 bright sources detected in the XMM-Newton deep survey of the Chandra Deep Field South, combining the three EPIC cameras. The 2-10 keV flux limit of the sample is 2e-15 erg/s/cm2. The sources are distributed over a redshift range of z=0.1-3.8. Eleven new X-ray redshift measurements are included. A spectral analysis was performed using a simple model to obtain absorbing column densities, rest-frame 2-10 keV luminosities and Fe K line properties of 180 sources at z>0.4. Obscured AGN are found to be more abundant toward higher redshifts. Using the XMM-Newton data alone, seven Compton-thick AGN candidates are identified, which makes the Compton-thick AGN fraction to be ~4%. An exploratory spectral inspection method with two rest-frame X-ray colours and an Fe line strength indicator is introduced and tested against the results from spectral fitting. This method works reasonably well to characterise a spectral shape and can be useful for a pre-selection of Compton-thick AGN candidates. We found six objects exhibiting broad Fe K lines out of 21 unobscured AGN of best data quality, implying a detection rate of ~30%. Five redshift spikes, each of which has more than six sources, are identified in the redshift distribution of the X-ray sources. Contrary to the overall trend, the sources at the two higher-redshift spikes at z=1.61 and z=2.57 have puzzlingly low obscuration., Comment: Accepted for publication in Astronomy & Astrophysics. Table 1 and Fig. 2 are available in electronic form the CDS. Typos corrected and some language editing and references added

We present the X-ray source catalog for the 479 ks Chandra exposure of the SDSS J1030+0524 field, that is centered on a region that shows the best evidence to date of an overdensity around a z > 6 quasar, and also includes a galaxy overdensity around a Compton-thick Fanaroff-Riley type II radio galaxy at z = 1.7. Using wavdetect for initial source detection and ACIS Extract for source photometry and significance assessment, we create preliminary catalogs of sources that are detected in the full, soft, and hard bands, respectively. We produce X-ray simulations that mirror our Chandra observation to filter our preliminary catalogs and get a completeness level of > 91% and a reliability level of 95% in each band. The catalogs in the three bands are then matched into a final main catalog of 256 unique sources. Among them, 244, 193, and 208 are detected in the full, soft, and hard bands, respectively. This makes J1030 field the fifth deepest extragalactic X-ray survey to date. The field is part of the Multiwavelength Survey by Yale-Chile (MUSYC), and is also covered by optical imaging data from the Large Binocular Camera (LBC) at the Large Binocular Telescope, near-IR imaging data from the Canada France Hawaii Telescope WIRCam, and Spitzer IRAC. Thanks to its dense multi-wavelength coverage, J1030 represents a legacy field for the study of large-scale structures around distant accreting supermassive black holes. Using a likelihood ratio analysis, we associate multi-band counterparts for 252 (98.4%) of the 256 Chandra sources, with an estimated reliability of 95%. Finally, we compute the cumulative number of sources in each X-ray band, finding that they are in general agreement with the results from the Chandra Deep Fields., Comment: A&A in press, 20 pages, 15 figures. The figures A.2 and and A.3 (postage-stamp images of the catalog) can be found online at: http://www.oabo.inaf.it/~LBTz6/1030/chandra_1030

We present a comparative analysis of the properties of AGN emitting at radio and X-ray wavelengths. The study is performed on 907 X-ray AGN and 100 radio AGN selected on the CDFS and UDS fields and makes use of new and ancillary data available to the VANDELS collaboration. Our results indicate that the mass of the host galaxy is a fundamental quantity which determines the level of AGN activity at the various wavelengths. Indeed large stellar masses are found to be connected with AGN radio emission, as virtually all radio-active AGN reside within galaxies of M*>10^{10} Msun. Large stellar masses also seem to favour AGN activity in the X-ray, even though X-ray AGN present a mass distribution which is more spread out and with a non-negligible tail at M*<10^{9} Msun. Stellar mass alone is also observed to play a fundamental role in simultaneous radio and X-ray emission: the percentage of AGN active at both wavelengths increases from around 1% of all X-ray AGN residing within hosts of M*<10^{11} Msun to about 13% in more massive galaxies. In the case of radio-selected AGN, such a percentage moves from about 15% to about 45% (but up to 80% in the deepest fields). Neither cosmic epoch, nor radio luminosity, X-ray luminosity, Eddington ratio or star-formation rate of the hosts are found to be connected to an enhanced probability for joint radio+X-ray emission of AGN origin. Furthermore, only a loose relation is observed between X-ray and radio luminosity in those AGN which are simultaneously active at both frequencies., Comment: 18 pages, 18 figures, to appear on MNRAS

Hyperluminous quasars ($L_{\rm bol}\gtrsim 10^{47}$ erg s$^{-1}$) are ideal laboratories to study the interaction and impact of extreme radiative field and the most powerful winds in the AGN nuclear regions. They typically exhibit low coronal X-ray luminosity ($L_{\rm X}$) compared to the UV and MIR radiative outputs ($L_{\rm UV}$ and $L_{\rm MIR}$) with a non-negligible fraction of them reporting even $\sim$1 dex weaker $L_{\rm X}$ compared to the prediction of the well established $L_{\rm X}$-$L_{\rm UV}$ and $L_{\rm X}$-$L_{\rm MIR}$ relations followed by the bulk of the AGN population. We report in our WISE/SDSS-selected Hyperluminous (WISSH) $z=2-4$ broad-line quasar sample, the discovery of a dependence between the intrinsic 2-10 keV luminosity ($L_{\rm 2-10}$) and the blueshifted velocity of the CIV emission line ($v_{\rm CIV}$) indicative of accretion disc winds. In particular, sources with fastest winds ($v_{\rm CIV}\gtrsim 3000~\rm km s^{-1}$) possess $\sim$0.5-1 dex lower $L_{\rm 2-10}$ than sources with negligible $v_{\rm CIV}$. No similar dependence is found on $L_{\rm UV}$, $L_{\rm MIR}$, $L_{\rm bol}$, photon index and absorption column density. We interpret these findings in the context of accretion disc wind models. Both magnetohydrodynamic and line-driven models can qualitatively explain the reported relations as a consequence of X-ray shielding from the inner wind regions. In case of line-driven winds, the launch of fast winds is favoured by a reduced X-ray emission, and we speculate that these winds may play a role in directly limiting the coronal hard X-ray production., Comment: 10 pages (including Appendix), 4 figures. Astronomy & Astrophysics Letters in press. Title numbering changed to match the published version

We look to provide a detailed description of the SPectroscopic IDentification of ERosita Sources (SPIDERS) survey, an SDSS-IV programme aimed at obtaining spectroscopic classification and redshift measurements for complete samples of sufficiently bright X-ray sources. We describe the SPIDERS X-ray Point Source Spectroscopic Catalogue, considering its store of 11,092 observed spectra drawn from a parent sample of 14,759 ROSAT and XMM sources over an area of 5,129 deg$^2$ covered in SDSS-IV by the eBOSS survey. This programme represents the largest systematic spectroscopic observation of an X-ray selected sample. A total of 10,970 (98.9\%) of the observed objects are classified and 10,849 (97.8\%) have secure redshifts. The majority of the spectra (10,070 objects) are active galactic nuclei (AGN), 522 are cluster galaxies, and 294 are stars. The observed AGN redshift distribution is in good agreement with simulations based on empirical models for AGN activation and duty cycle. Forming composite spectra of type 1 AGN as a function of the mass and accretion rate of their black holes reveals systematic differences in the H-beta emission line profiles. This study paves the way for systematic spectroscopic observations of sources that are potentially to be discovered in the upcoming eROSITA survey over a large section of the sky., Comment: accepted in A&A Please visit our project page : http://www.mpe.mpg.de/XraySurveys/SPIDERS The catalogs are here : https://data.sdss.org/sas/dr16/eboss/spiders/analysis/VAC_SPIDERS_2RXS_DR16.fits and https://data.sdss.org/sas/dr16/eboss/spiders/analysis/VAC_SPIDERS_XMMSL2_DR16.fits

We report the discovery of a galaxy overdensity around a Compton-thick Fanaroff-Riley type II (FRII) radio galaxy at z=1.7 in the deep multiband survey around the z=6.3 QSO SDSS J1030+0524. Based on a 6hr VLT/MUSE and on a 4hr LBT/LUCI observation, we identify at least eight galaxy members in this structure with spectroscopic redshift z=1.687-1.699, including the FRII galaxy at z=1.699. Most of the identified overdensity members are blue, compact galaxies that are actively forming stars at rates of $\sim$8-60 $M_{\odot}$ yr$^{-1}$. Based on a 500ks Chandra ACIS-I observation we found that the FRII nucleus hosts a luminous QSO ($L_{2-10keV}=1.3\times10^{44}$ erg s$^{-1}$, intrinsic and rest-frame) that is obscured by Compton-thick absorption ($N_H=1.5\pm0.6\times 10^{24}$ cm$^{-2}$). Our Chandra observation, the deepest so far for a distant FRII within a galaxy overdensity, revealed significant diffuse X-ray emission within the region covered by the overdensity. In particular, X-ray emission extending for $\sim$240 kpc is found around the Eastern lobe of the FRII. Four out of the six MUSE star forming galaxies in the overdensity are distributed in an arc-like shape at the edge of this diffuse X-ray emission. The probability of observing by chance four out of the six $z=1.7$ sources at the edge of the diffuse emission is negligible. In addition, these four galaxies have the highest specific star formation rates of the MUSE galaxies in the overdensity and lie above the main sequence of field galaxies of equal stellar mass at z=1.7. We propose that the diffuse X-rays originate from an expanding bubble of gas that is shock-heated by the FRII jet, and that star formation is promoted by the compression of the cold interstellar medium of the galaxies around the bubble, which may be remarkable evidence of positive AGN feedback on cosmological scales. [shortened version], Comment: 17 pages, 13 figures. Accepted for publication in A&A

X-ray emission from AGN is believed to be produced via Comptonization of optical/UV seed photons emitted by the accretion disk, up-scattered by hot electrons in a corona surrounding the black hole. A critical compactness vs. temperature threshold is predicted above which any increase in the source luminosity, for a fixed size, would then generate positron-electron pairs rather than continue heating the coronal plasma. Current observations seem to confirm that all AGN populate the region below this critical line. These models, however, have never been probed by observations in the high-luminosity regime, where the critical line is expected to reach low temperatures. To fill this observational gap, we selected two luminous ($log(L_{bol})>47.5$ erg/s) quasars, 2MASSJ1614346+470420 ($z=1.86$) and B1422+231 ($z=3.62$), and obtained XMM-Newton and NuSTAR deep observations for them. We performed detailed spectral analysis of their quasi-simultaneous soft and hard X-ray data, in order to constrain the parameters of their coronae. Using a phenomenological cut-off power-law model, with the inclusion of a reflection component, we derived rest-frame values of the high energy cut-off of $E_{cut}=106^{+102}_{-37}$ keV and $E_{cut}=66^{+17}_{-12}$ keV, respectively. Comptonization models consistently give as best-fit parameters electron temperatures of $\sim45$ keV and $\sim28$ keV, respectively, and optically thick coronae ($\tau>1$). These low coronal temperatures fall in the limited region allowed at these luminosities to avoid runaway pair production., Comment: 8 pages, 2 figures, accepted for publication in ApJ Letter

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., Comment: Accepted in A&A

Powerful winds driven by supermassive black holes (SMBHs) are likely the main mechanism through which SMBHs regulate their own growth and influence the host galaxy evolution. However, their origin and their capability to impact the large-scale environment are still highly debated. Fundamental results will come from high-energy and spatial resolution X-ray observatories., Comment: Science White Paper Submitted for the Astro2020 Decadal Survey on Astronomy and Astrophysics

In the context of the upcoming SRG/eROSITA survey, we present an N-body simulation-based mock catalogue for X-ray selected AGN samples. The model reproduces the observed hard X-ray AGN luminosity function (XLF) and the soft X-ray logN-logS from redshift 0 to 6. The XLF is reproduced to within $\pm5\%$ and the logN-logS to within $\pm20\%$. We develop a joint X-ray -- optical extinction and classification model. We adopt a set of empirical spectral energy distributions to predict observed magnitudes in the UV, optical and NIR. With the latest eROSITA all sky survey sensitivity model, we create a high-fidelity full-sky mock catalogue of X-ray AGN. It predicts their distributions in right ascension, declination, redshift and fluxes. Using empirical medium resolution optical spectral templates and an exposure time calculator, we find that $1.1\times10^6$ ($4\times10^5$) fiber-hours are needed to follow-up spectroscopically from the ground the detected X-ray AGN with an optical magnitude $21

We investigated the interstellar medium properties of the disc and outflowing gas in the central regions of nine nearby Seyfert galaxies, all characterised by prominent outflows. These objects are part of the Measuring Active Galactic Nuclei Under MUSE Microscope survey, which aims to probe their physical conditions and ionisation mechanism by exploiting MUSE unprecedented sensitivity. We studied the different properties of the gas in the disc and outflow with spatially and kinematically resolved maps by dividing the strongest emission lines in velocity bins. We associated the core of the lines with the disc, consistent with the stellar velocity, and the redshifted and the blueshifted wings with the outflow. We find that the outflowing gas is characterised by higher values of density and ionisation parameter than the disc, which presents a higher dust extinction. Moreover, we distinguish high- and low-ionisation regions across the portion of spatially resolved narrow-line region traced by the outflowing gas. The high-ionisation regions characterised by the lowest [NII]/H{\alpha} and [SII]/H{\alpha} line ratios generally trace the innermost parts along the axis of the emitting cones where the [SIII]/[SII] line ratio is enhanced, while the low-ionisation regions follow the cone edges and/or the regions perpendicular to the axis of the outflows, also characterised by a higher [OIII] velocity dispersion. A possible scenario to explain these features relies on the presence of two distinct populations of line emitting clouds: one is optically thin to the radiation and is characterised by the highest excitation, while the other is optically thick and is impinged by a filtered, and thus harder, radiation field which generates strong low-excitation lines. The highest values of [NII]/H{\alpha} and [SII]/H{\alpha} line ratios may be due to shocks and/or a hard filtered radiation field from the AGN., Comment: 28 pages, 19 figures. Forthcoming article 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

Ionized outflows, revealed by broad asymmetric wings of the [OIII] line, are commonly observed in AGN but the low intrinsic spatial resolution of observations has generally prevented a detailed characterization of their properties. The MAGNUM survey aims at overcoming these limitations by focusing on the nearest AGN, including NGC 1365, a nearby Seyfert galaxy (D~17 Mpc), hosting a low-luminosity AGN (Lbol ~ 2x10^43 erg/s). We want to obtain a detailed picture of the ionized gas in the central ~5 kpc of NGC 1365 in terms of physical properties, kinematics, and ionization mechanisms. We also aim to characterize the warm ionized outflow as a function of distance from the nucleus and its relation with the nuclear X-ray wind. We employed VLT/MUSE optical integral field spectroscopic observations to investigate the warm ionized gas and Chandra ACIS-S X-ray data for the hot highly-ionized phase. We obtained flux, kinematic, and diagnostic maps of the optical emission lines, which we used to disentangle outflows from disk motions and measure the gas properties down to a spatial resolution of ~70 pc. [OIII] emission mostly traces an AGN-ionized kpc-scale biconical outflow with velocities up to ~200 km/s. H{\alpha} emission traces instead star formation in a circumnuclear ring and along the bar, where we detect non-circular motions. Soft X-rays are mostly due to thermal emission from the star-forming regions, but we could isolate the AGN photoionized component which matches the [OIII] emission. The mass outflow rate of the extended ionized outflow matches that of the nuclear X-ray wind and then decreases with radius. However, the hard X-ray emission from the circumnuclear ring suggests that star formation might contribute to the outflow. The integrated mass outflow rate, kinetic energy rate, and outflow velocity are broadly consistent with the typical relations observed in more luminous AGN., Comment: Accepted for publication on Astronomy and Astrophysics. 22 pages, 14 figures