Search

We aim at characterizing the X-ray-to-optical/near-infrared broad-band emission of luminous QSOs in the first Gyr of cosmic evolution to understand whether they exhibit differences compared to the lower-\textit{z} QSO population. Our goal is also to provide for these objects a reliable and uniform catalog of SED fitting derivable properties such as bolometric and monochromatic luminosities, Eddington ratios, dust extinction, strength of the hot dust emission. We characterize the X-ray/UV emission of each QSO using average SEDs from luminous Type 1 sources and calculate bolometric and monochromatic luminosities. Finally we construct a mean SED extending from the X-rays to the NIR bands. We find that the UV-optical emission of these QSOs can be modelled with templates of $z\sim$2 luminous QSOs. We observe that the bolometric luminosities derived adopting some bolometric corrections at 3000 \AA\ ($BC_{3000\text{\AA}}$) largely used in the literature are slightly overestimated by 0.13 dex as they also include reprocessed IR emission. We estimate a revised value, i.e. $BC_{3000\text{\AA}}=3.3 $ which can be used for deriving $L_\text{bol}$ in \textit{z} $\geq$ 6 QSOs. A sub-sample of 11 QSOs is provided with rest-frame NIR photometry, showing a broad range of hot dust emission strength, with two sources exhibiting low levels of emission. Despite potential observational biases arising from non-uniform photometric coverage and selection biases, we produce a X-ray-to-NIR mean SED for QSOs at \textit{z} $\gtrsim$ 6, revealing a good match with templates of lower-redshift, luminous QSOs up to the UV-optical range, with a slightly enhanced contribution from hot dust in the NIR.

One of the main open questions in the field of luminous ($L_{\rm bol}>10^{47}\,\rm erg\,s^{-1}$) quasars (QSOs) at $z \gtrsim 6$ is the rapid formation ($< 1\,$Gyr) of their supermassive black holes (SMBHs). For this work we analysed the relation between the X-ray properties and other properties describing the physics and growth of both the accretion disc and the SMBH in QSOs at the Epoch of Reionization (EoR). The sample consists of 21 $z>6$ QSOs, which includes 16 sources from the rapidly grown QSOs from the HYPERION sample and five other luminous QSOs with available high-quality archival X-ray data. We discovered a strong and statistically significant ($>3\sigma$) relation between the X-ray continuum photon index ($\Gamma$) and the $\rm C\,IV$ disc wind velocity ($v_{\rm C\,IV}$) in $z>6$ luminous QSOs, whereby the higher the $v_{\rm C\,IV}$, the steeper the $\Gamma$. This relation suggests a link between the disc-corona configuration and the kinematics of disc winds. Furthermore, we find evidence at $>2-3\sigma$ level that $\Gamma$ and $v_{\rm C\,IV}$ are correlated to the growth rate history of the SMBH. Although additional data are needed to confirm it, this result may suggest that, in luminous $z>6$ QSOs, the SMBH predominantly grows via fast accretion rather than via initial high seed BH mass., Comment: Accepted for publication in Astronomy & Astrophysics

The energy released by AGN has the potential to heat or remove the gas of the ISM, thus likely impacting the cold molecular gas reservoir of host galaxies at first, with star formation following on longer timescales. Previous works on high-z galaxies have yielded conflicting results, possibly due to selection biases and other systematics. To provide a reliable benchmark for galaxy evolution models at cosmic noon (z=1-3), two surveys were conceived: SUPER and KASHz, both targeting unbiased X-ray-selected AGN at z>1 that span a wide bolometric luminosity range. In this paper, we assess the effects of AGN feedback on the molecular gas content of host galaxies in a statistically robust, uniformly selected, coherently analyzed sample of AGN at z=1-2.6, drawn from the KASHz and SUPER surveys. By using ALMA data in combination with dedicated SED modeling, we retrieve CO and FIR luminosity as well as $M_*$ of SUPER and KASHz AGN. We selected non-active galaxies from PHIBBS, ASPECS and multiple ALMA/NOEMA surveys of sub-mm galaxies. By matching the samples in z, $M_*$ and $L_{FIR}$, we compared the properties of AGN and non-active galaxies within a Bayesian framework. We find that AGN hosts at given $L_{FIR}$ are on average CO depleted compared to non-active galaxies, confirming what was previously found in the SUPER survey. Moreover, the molecular gas fraction distributions of AGN and non-active galaxies are statistically different, with that of of AGN being skewed to lower values. Our results indicate that AGN can indeed reduce the total cold molecular gas reservoir of their host galaxies. Lastly, by comparing our results with predictions from three cosmological simulations (TNG, Eagle and Simba) filtered to match the observed properties, we confirm already known discrepancies and highlight new ones between observations and simulations.[Abridged], Comment: 26 pages (17 in the main body of the paper and 9 in the appendix), 14 figures (11 in the main body of the paper and 3 in the appendix), 8 tables (2 in the main body of the paper and 6 in the appendix). Accepted for publication in A&A, in press

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

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

We used low- to high-frequency ALMA observations to investigate the cold gas and dust in ten QSOs at $z\gtrsim 6$. Our analysis of the CO(6-5) and CO(7-6) emission lines in the selected QSOs provided insights into their molecular gas masses, which average around $10^{10}\ \rm M_\odot$, consistent with typical values for high-redshift QSOs. Proprietary and archival ALMA observations in bands 8 and 9 enabled precise constraints on the dust properties and star formation rate (SFR) of four QSOs in our sample for the first time. The examination of the redshift distribution of dust temperatures revealed a general trend of increasing $T_{\rm dust}$ with redshift, which agrees with theoretical expectations. We computed a mean cold dust spectral energy distribution considering all ten QSOs. This offers a comprehensive view of the dust properties of high-$z$ QSOs. The QSOs marked by a more intense growth of the supermassive black hole (HYPERION QSOs) showed lower dust masses and higher gas-to-dust ratios on average, but their $\rm H_2$ gas reservoirs are consistent with those of other QSOs at the same redshift. The observed high SFR in our sample yields high SF efficiencies and thus very short gas depletion timescales ($\tau_{\rm dep}\sim 10^{-2}$ Gyr). Beyond supporting the paradigm that high-$z$ QSOs reside in highly star-forming galaxies, our findings portrayed an interesting evolutionary path at $z>6$. Our study suggests that they are undergoing rapid galaxy growth that might be regulated by strong outflows. Their inferred evolutionary path shows a convergence toward the massive end of the local relation, which supports the idea that they are candidate progenitors of local massive galaxies. The observed pathway involves intense BH growth followed by substantial galaxy growth, in contrast with a symbiotic growth scenario. The abstract has been shortened (full version in the article)., Comment: 29 pages; 6 tables; 18 figures. Accepted by A&A. A section about SF efficiency has been added compared to the previous version

One of the key open questions in Cosmology is the nature of the sources that completed the cosmological hydrogen Reionization at z~5.2. High-z primeval galaxies have been long considered the main drivers for Reionization, with a minor role played by high-z AGN. However, in order to confirm this scenario, it is fundamental to measure the photo-ionization rate produced by active SMBHs close to the epoch of Reionization. Given the pivotal role played by spectroscopically complete observations of high-z QSOs, in this paper we present the first results of the RUBICON (Reionizing the Universe with BrIght COsmological Nuclei) survey. It consists of a color selected sample of bona-fide z~5 QSO candidates from the Hyper Suprime-Cam Subaru Strategic Survey. Our QSO candidates have been validated both by photometric redshifts based on SED fitting and by spectroscopic redshifts, confirming that they lie at 4.5

We present bolometric luminosities, black hole masses and Eddington ratios for 42 luminous quasars at z>6 using high signal-to-noise ratio VLT/X-Shooter spectra, acquired in the enlarged ESO Large Programme XQR-30. In particular, we derive bolometric luminosities from the rest-frame 3000 A, luminosities using a bolometric correction from the literature, and the black hole masses by modelling the spectral regions around the CIV 1549A and the MgII 2798A emission lines, with scaling relations calibrated in the local universe. We find that the black hole masses derived from both emission lines are in the same range, and the scatter of the measurements agrees with expectations from the scaling relations. The MgII-derived masses are between ~(0.8-12) x 10^9 Msun, and the derived Eddington ratios are within ~0.13-1.73, with a mean (median) of 0.84 (0.72). By comparing the total sample of quasars at z>5.8, from this work and from the literature, to a bolometric luminosity distribution-matched sample at z~1.5, we find that quasars at high redshift host slightly less massive black holes which accrete slightly more rapidly than at lower-z, with a difference in the mean Eddington ratios of the two samples of ~0.27, in agreement with recent literature work., Comment: 9 pages; 5 figures; accepted for publication in A&A

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

The final phase of the reionization process can be probed by rest-frame UV absorption spectra of quasars at z>6, shedding light on the properties of the diffuse intergalactic medium within the first Gyr of the Universe. The ESO Large Programme "XQR-30: the ultimate XSHOOTER legacy survey of quasars at z~5.8-6.6" dedicated ~250 hours of observations at the VLT to create a homogeneous and high-quality sample of spectra of 30 luminous quasars at z~6, covering the rest wavelength range from the Lyman limit to beyond the MgII emission. Twelve quasar spectra of similar quality from the XSHOOTER archive were added to form the enlarged XQR-30 sample, corresponding to a total of ~350 hours of on-source exposure time. The median effective resolving power of the 42 spectra is R~11400 and 9800 in the VIS and NIR arm, respectively. The signal-to-noise ratio per 10 km/s pixel ranges from ~11 to 114 at $\lambda \simeq 1285$ \AA rest frame, with a median value of ~29. We describe the observations, data reduction and analysis of the spectra, together with some first results based on the E-XQR-30 sample. New photometry in the H and K bands are provided for the XQR-30 quasars, together with composite spectra whose characteristics reflect the large absolute magnitudes of the sample. The composite and the reduced spectra are released to the community through a public repository, and will enable a range of studies addressing outstanding questions regarding the first Gyr of the Universe., Comment: 21 pages, 10 figures. Final version accepted by MNRAS

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

In order to investigate the impact of radio jets on the interstellar medium (ISM) of galaxies hosting active galactic nuclei (AGN), we present subarcsecond resolution Atacama Large Millimeter/submillimeter Array (ALMA) CO(2-1) and CO(3-2) observations of the Teacup galaxy. This is a nearby ($D_{\rm L}$=388 Mpc) radio-quiet type-2 quasar (QSO2) with a compact radio jet ($P_{\rm jet}\approx$10$^{43}$ erg s$^{-1}$) that subtends a small angle from the molecular gas disc. Enhanced emission line widths perpendicular to the jet orientation have been reported for several nearby AGN for the ionised gas. For the molecular gas in the Teacup, not only do we find this enhancement in the velocity dispersion but also a higher brightness temperature ratio (T32/T21) perpendicular to the radio jet compared to the ratios found in the galaxy disc. Our results and the comparison with simulations suggest that the radio jet is compressing and accelerating the molecular gas, and driving a lateral outflow that shows enhanced velocity dispersion and higher gas excitation. These results provide further evidence that the coupling between the jet and the ISM is relevant to AGN feedback even in the case of radio-quiet galaxies., Comment: 9 pages, 8 figures. Accepted for publication in A&A Letters

We present spatially resolved H$\alpha$ properties of 21 type 1 AGN host galaxies at z$\sim$2 derived from the SUPER survey. These targets were observed with the adaptive optics capabilities of the SINFONI spectrograph, a near-infrared integral field spectrograph, that provided a median spatial resolution of 0.3 arcsec ($\sim$2 kpc). We model the H$\alpha$ emission line profile in each pixel to investigate whether it traces gas in the narrow line region or if it is associated with star formation. To do this, we first investigate the presence of resolved H$\alpha$ emission by removing the contribution of the AGN PSF. We find extended H$\alpha$ emission in sixteen out of the 21 type 1 AGN host galaxies (76%). Based on the BPT diagnostics, optical line flux ratios and the line widths (FWHM), we show that the H$\alpha$ emission in five galaxies is ionised by the AGN (30%), in four galaxies by star formation (25%) and for the rest (45%), the ionisation source is unconstrained. Two galaxies show extended H$\alpha$ FWHM $>$600 km/s, which is interpreted as a part of an AGN-driven outflow. Morphological and kinematic maps of H$\alpha$ emission in targets with sufficient signal-to-noise ratio suggest the presence of rotationally supported disks in six galaxies and possible presence of companions in four galaxies. In two galaxies, we find an anti-correlation between the locations of extended H$\alpha$ emission and [OIII]-based ionised outflows, indicating possible negative feedback at play. However, in the majority of galaxies, we do not find evidence of outflows impacting H$\alpha$ based star formation., Comment: 16 pages, 11 figures, 2 tables, accepted for publication in MNRAS; Appendix = 2 tables and 3 figures

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

The outstanding mass growth of supermassive black holes (SMBHs) at the Reionisation Epoch and how it is related to the concurrent growth of their host galaxies, poses challenges to theoretical models aimed at explaining how these systems formed in short timescales (<1 Gyr). To trace the average evolutionary paths of quasi-stellar objects (QSOs) and their host galaxies in the BH mass-host mass ($M_{\rm dyn}$) plane, we compare the star formation rate (SFR), derived from the accurate estimate of the dust temperature and the dust mass ($T_{\rm dust}, M_{\rm dust}$), with the BH accretion rate. To this aim, we analysed a deep, $900$ pc resolution ALMA observation of the sub-mm continuum, [CII] and H$_2$O of the $z\sim 6$ QSO J2310+1855, enabling a detailed study of dust properties and cold gas kinematics. We performed an accurate SED analysis obtaining a dust temperature of $T_{\rm dust} = 71$ K and a dust mass of $M_{\rm dust}= 4.4 \times 10^8\ \rm M_{\odot}$. The implied AGN-corrected SFR is $1240 \ \rm M_{\odot}yr^{-1}$, a factor of 2 smaller than previously reported for this QSO. We derived the best estimate of the dynamical mass $M_{\rm dyn} = 5.2\times 10^{10}\ \rm M_{\odot}$ within $r = 1.7$ kpc, based on a dynamical model of the system. We found that ${\rm SFR}/M_{\rm dyn}>\dot M_{\rm BH}/M_{\rm BH}$, suggesting that AGN feedback might be efficiently acting to slow down the SMBH accretion, while the stellar mass assembly is still vigorously taking place in the host galaxy. In addition, we were also able to detect high-velocity emission on the red and blue sides of the [CII] emission line, that traces a gaseous outflow, and for the first time, we mapped a spatially-resolved water vapour disk through the H$_2$O v=0 $3_{(2,2)}-3_{(1,3)}$ emission line detected at $\nu_{\rm obs} = 274.074$ GHz, whose kinematic properties and size are broadly consistent with those of the [CII] disk., Comment: 15 pages, 10 figures, 4 tables. Accepted in A&A

We report on the variability of a multi-component broad absorption line (BAL) system observed in the hyper-luminous quasar J1538+0855 at z=3.6. Observations from SDSS, VLT, LBT and Subaru telescopes taken at five different epochs, spanning 17 yr in the observed frame, are presented. We detect three (A, B, C) CIV variable troughs exhibiting extreme velocities ($\sim$40,000-54,000 km s$^{-1}$) similar to the ultra-fast outflows (UFOs) typically observed in the X-ray spectra. The A component of the BAL UFO ($\rm v_{ufo}$ $\sim$0.17 c) shows strength variations, while B ($\rm v_{ufo}$ $\sim$0.15 c) and C ($\rm v_{ufo}$ $\sim$0.13 c) components show changes both in shape and strength, appearing and disappearing at different epochs. In addition, during the last observation on June 2021 the entire BAL system disappears. The variability trends observed during the first two epochs (1.30 yr rest-frame) in the CIV, SiIV, OVI and NV absorption spectral regions are the same for B and C troughs, while the A component of the BAL varies independently. This suggests a change in the ionization state of the absorbing gas for B and C components and tangential motion for the A component, as causes of this temporal behavior. Accordingly, it is possible to provide an upper limit for distance of the gas responsible for the A component of $R\rm_{out}^{A}$$\le$58 pc, and in turn, a kinetic power of $\dot{E}\rm_{K,ufo}$ $\le$5.2 $\times$ 10$^{44}$ erg s$\rm^{-1}$. We also obtain $R\rm_{out}^{B,C}$ $\le$2.7 kpc for B and C components, which implies an upper limit estimation of $\dot{E}\rm_{K,ufo}$ $\le$2.1$\times$10$^{46}$ erg s$\rm^{-1}$ and $\dot{E}\rm_{K,ufo}$ $\le$1.4$\times$10$^{46}$ erg s$\rm^{-1}$, respectively. Future spectral monitoring with high-resolution instruments is mandatory to accurately constrain physical properties of the BAL UFO discovered in the UV spectrum of J1538+0855., Comment: Accepted for publication in A&A

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

We present CO(2-1) and adjacent continuum observations of 7 nearby radio-quiet type-2 quasars (QSO2s) obtained with ALMA at ~0.2" resolution (370 pc at z~0.1). The CO morphologies are diverse, including disks and interacting systems. Two of the QSO2s are red early-type galaxies with no CO(2-1) detected. In the interacting galaxies, the central kpc contains 18-25% of the total cold molecular gas, whereas in the spirals it is only 5-12%. J1010+0612 and J1430+1339 show double-peaked CO morphologies which do not have optical counterparts. Based on our analysis of the ionized and molecular kinematics and mm continuum emission, these CO morphologies are most likely produced by AGN feedback in the form of outflows, jets, and/or shocks. The CO kinematics of the QSO2s are dominated by rotation but also reveal noncircular motions. According to our analysis of the kinematics, these noncircular motions correspond to molecular outflows mostly coplanar with the CO discs in four of the QSO2s, and either to a coplanar inflow or vertical outflow in the case of J1010+0612. These outflows represent 0.2-0.7% of the QSO2s' total molecular gas mass and have maximum velocities of 200-350 km/s, radii from 0.4 to 1.3 kpc, and outflow rates of 8-16 Msun/yr. These properties are intermediate between those of the mild molecular outflows measured for Seyferts, and the fast and energetic outflows of ULIRGs. This suggests that it is not only AGN luminosity that drives massive molecular outflows. Other factors such as jet power, coupling between winds, jets, and/or ionized outflows and the CO discs, and amount or geometry of dense gas in the nuclear regions might be also relevant. Thus, although we do not find evidence for a significant impact of quasar feedback on the total molecular gas reservoirs and SFRs, it appears to be modifying the distribution of cold molecular gas in the central kpc of the galaxies., Comment: 29 pages including 18 figures, 7 tables and appendices. A&A accepted version including minor edits

Improving our understanding of the nuclear properties of high-Eddington ratio ($\lambda_\mathrm{Edd}$) active galactic nuclei (AGN) is necessary since the bulk of X-ray spectroscopic studies have been focused on low-Eddington AGN. We present here the X-ray spectral analysis of 14 radio-quiet, $\lambda_\mathrm{Edd}\gtrsim1$ AGN at $0.4\leq z \leq 0.75$, observed with XMM-Newton. Optical/UV data from simultaneous Optical Monitor observations have been also considered. These AGN have been selected to have relatively high values of black hole mass ($M_\mathrm{BH}\sim10^{8-8.5}M_\odot$) and bolometric luminosity ($L_\mathrm{bol} \sim 10^{46}$ erg s$^{-1}$), in order to complement previous studies of high-Eddington AGN at lower $M_\mathrm{BH}$ and $L_\mathrm{bol}$. We studied the relation between $\lambda_\mathrm{Edd}$ and other key X-ray spectral parameters, such as the photon index of the power-law continuum $\Gamma$, the X-ray bolometric correction $k_\mathrm{bol,X}$ and $\alpha_\mathrm{ox}$. Despite the homogeneous optical and SMBH accretion properties, the X-ray properties of these high-Eddington AGN are quite heterogeneous. We measured values of $\Gamma$ comprised between 1.3 and 2.5, at odds with the expectations based on previously reported $\Gamma-\lambda_\mathrm{Edd}$ relations, by which $\Gamma\geq2$ would be an ubiquitous hallmark of AGN with $\lambda_\mathrm{Edd}\sim1$. We found that $\sim30\%$ of the sources are X-ray weak, with an X-ray emission about a factor of $\sim10-80$ fainter than that of typical AGN at similar UV luminosities. The X-ray weakness seems to be intrinsic and not due to intervening obscuration. This may indicate that high-Eddington AGN commonly undergo periods of intrinsic X-ray weakness. Furthermore, results from a follow-up monitoring with Swift of one of these X-ray weak AGN suggest that these periods can last for several years., Comment: 17 pages, 8 figures (with 2 more figures in Appendix A and B). Accepted for publication in A&A

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

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. The ALMA and SINFONI maps presented in this paper are available at: https://doi.org/10.25405/data.ncl.23694765

We present the analysis of the ALMA CO(2-1) emission line and the underlying 1.2 mm continuum of Mrk509 with spatial resolution of 270 pc. This local Seyfert 1.5 galaxy, optically classified as a spheroid, is known to host a ionised disc, a starburst ring, and ionised gas winds on both nuclear and galactic scales. From CO(2-1) we estimate a molecular gas mass $M_{H_2}=1.7\times 10^9\, \rm M_{\odot}$, located within a disc of size 5.2 kpc, with $M_{dyn}$=(2.0$\pm$1.1) $\times$ $10^{10}\, \rm M_{\odot}$ inclined at $44\pm10$ deg. The molecular gas fraction within the disc is $\mu_{gas}=5\%$. The gas kinematics in the nuclear region within r=700 pc suggests the presence of a warped nuclear disc. Both the presence of a molecular disc with ongoing star-formation in a starburst ring, and the signatures of a minor merger, are in agreement with the scenario where galaxy mergers produce gas destabilization, feeding both star-formation and AGN activity. The spatially-resolved Toomre Q-parameter across the molecular disc is in the range $Q_{gas}=0.5-10$, and shows that the disc is marginally unstable across the starburst ring, and stable at nucleus and in a lopsided ring-like structure located inside of the starburst ring. We find complex molecular gas kinematics and significant kinematics perturbations at two locations, one within 300 pc from the nucleus, and one 1.4 kpc away close to the region with high $Q_{gas}$, that we interpret as molecular winds. The total molecular outflow rate is in the range 6.4-17.0 $\rm M_\odot/yr$. The molecular wind total kinetic energy is consistent with a multiphase momentum-conserving wind driven by the AGN with $\dot{P}_{of}/\dot{P}_{rad}$ in the range 0.06-0.5. The spatial overlap of the inner molecular wind with the ionised wind, and their similar velocity suggest a cooling sequence within a multiphase AGN driven wind., Comment: 13 pages, 11 figures, 3 tables Accepted by A&A

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

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

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., Comment: 15 pages, 10 figures, 5 tables. Accepted to MNRAS in 2020 August 14

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

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

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

We present a detailed study of ionized outflows in a large sample of ~650 hard X-ray detected AGN. Using optical spectroscopy from the BAT AGN Spectroscopic Survey (BASS) we are able to reveal the faint wings of the [OIII] emission lines associated with outflows covering, for the first time, an unexplored range of low AGN bolometric luminosity at low redshift (z~0.05). We test if and how the incidence and velocity of ionized outflow is related to AGN physical parameters: black hole mass, gas column density, Eddington Ratio, [OIII], X-ray, and bolometric luminosities. We find a higher occurrence of ionized outflows in type 1.9 (55%) and type 1 AGN (46%) with respect to type 2 AGN (24%). While outflows in type 2 AGN are evenly balanced between blue and red velocity offsets with respect to the [OIII] narrow component, they are almost exclusively blueshifted in type 1 and type 1.9 AGN. We observe a significant dependence between the outflow occurrence and accretion rate, which becomes relevant at high Eddington ratios (> -1.7). We interpret such behaviour in the framework of covering factor-Eddington ratio dependence. We don't find strong trends of the outflow maximum velocity with AGN physical parameters, as an increase with bolometric luminosity can be only identified when including samples of AGN at high luminosity and high redshift taken from literature., Comment: Accepted for publication in MNRAS, 15 pages, 19 figures, 1 table

We present ALMA rest-frame 230 GHz continuum and CO(2-1) line observations of the nearby Compton-thick Seyfert galaxy ESO428-G14, with angular resolution 0.7 arcsec (78 pc). We detect CO(2-1) emission from spiral arms and a circum-nuclear ring with 200 pc radius, and from a transverse gas lane with size of $\sim100$ pc, which crosses the nucleus and connects the two portions the circumnuclear ring. The molecular gas in the host galaxy is distributed in a rotating disk with intrinsic circular velocity $v_{rot}=135$ km/s, inclination $i=57$ deg, and dynamical mass $M_{dyn }=5\times 10^9~\rm M_{\odot}$ within a radius of $\sim 1$ kpc. In the inner 100 pc region CO is distributed in a equatorial bar, whose kinematics is highly perturbed and consistent with an inflow of gas towards the AGN. This inner CO bar overlaps with the most obscured, Compton-thick region seen in X-rays. We derive a column density of $\rm N(H_2) \approx 2\times10^{23}~ cm^{-2}$ in this region, suggesting that molecular gas may contribute significantly to the AGN obscuration. We detect a molecular outflow with a total outflow rate $\rm \dot M_{of}\approx 0.8~M_{\odot}/yr$, distributed along a bi-conical structure with size of $700$ pc on both sides of the AGN. The bi-conical outflow is also detected in the $\rm H_2$ emission line at 2.12 $\mu$m, which traces a warmer nuclear outflow located within 170 pc from the AGN. This suggests that the outflow cools with increasing distance from the AGN. We find that the hard X-ray emitting nuclear region mapped with Chandra is CO-deprived, but filled with warm molecular gas traced by $\rm H_2$ - thus confirming that the hard (3-6 keV) continuum and Fe K$\alpha$ emission are due to scattering from dense neutral clouds in the ISM., Comment: Submitted to ApJ

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

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

We present the stacking analysis of a sample of 48 QSOs at 4.5

We present the analysis of the rest-frame optical-to-UV spectrum of APM 08279+5255, a well-known lensed broad absorption line (BAL) quasar at $z = 3.911$. The spectroscopic data are taken with the optical DOLoRes and near-IR NICS instruments at TNG, and include the previously unexplored range between C III] $\lambda$1910 and [O III] $\lambda\lambda$4959,5007. We investigate the possible presence of multiple BALs by computing "balnicity" and absorption indexes (i.e. BI, BI$_0$ and AI) for the transitions Si IV $\lambda$1400, C IV $\lambda$1549, Al III $\lambda$1860 and Mg II $\lambda$2800. No clear evidence for the presence of absorption features is found in addition to the already known, prominent BAL associated to C IV, which supports a high-ionization BAL classification for APM 08279+5255. We also study the properties of the [O III], H$\beta$ and Mg II emission lines. We find that [O III] is intrinsically weak ($F_{\rm [OIII]}/F_{\rm H\beta} \lesssim 0.04$), as it is typically found in luminous quasars with a strongly blueshifted C IV emission line ($\sim$2500 km s$^{-1}$ for APM 08279+5255). We compute the single-epoch black hole mass based on Mg II and H$\beta$ broad emission lines, finding $M_{\rm BH} = (2 \div 3) \times 10^{10}\mu^{-1}$ M$_\odot$, with the magnification factor $\mu$ that can vary between 4 and 100 according to CO and rest-frame UV-to-mid-IR imaging respectively. Using a Mg II equivalent width (EW)-to-Eddington ratio relation, the EW$_{\rm MgII} \sim 27$ \AA\ measured for APM 08279+5255 translates into an Eddington ratio of $\sim$0.4, which is more consistent with $\mu=4$. This magnification factor also provides a value of $M_{\rm BH}$ that is consistent with recent reverberation-mapping measurements derived from C IV and Si IV., Comment: 10 pages, 4 figures, 4 tables, accepted for publication in A&A

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 (<150 km/s) velocity shift with respect to the QSO redshift. We report the discovery of the closest companion observed so far in submillimetre observations of high-z QSOs. It is only 2.2 kpc distant and merging with J1015+0020, while two other [CII] emitters are found at 8 and 17 kpc. Two strong continuum emitters are also detected at <3.5 arcsec. They are likely associated to the same overdense structure of J1015+0020, as they exceed by a factor of 100 the number of expected sources, considering the Log(N)-Log(S) at 850 um. The host galaxy of J1015+0020 shows a SFR of about 100 Msun/yr while the total SFR of the QSO and its companion galaxies is a factor of 10 higher, indicating that substantial stellar mass assembly at early epochs may have taken place in the QSO satellites. For J1015+0020 we compute a SMBH mass MBH~6E9 Msun and a dynamical mass Mdyn~4E10 Msun . This translates into an extreme ratio Mdyn/MBH~7. The total stellar mass of the QSO host galaxy plus the [CII] emitters already exceeds 1E11 Msun at z~4.4. These sources will likely merge and develop into a giant galaxy of 1.3E12 Msun. Under the assumption of constant mass or Eddington accretion rate equal to the observed values, we find that the growth timescale of the host galaxy is comparable or even shorter than that inferred for the SMBH., Comment: Accepted

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 have undertaken a multi-band observing program aimed at obtaining a complete census of winds in a sample of WISE/SDSS selected hyper-luminous (WISSH) QSOs at z~2-4. We have analyzed the rest-frame optical (LBT/LUCI and VLT/SINFONI) and UV (SDSS) spectra of 18 randomly selected WISSH QSOs to measure the SMBH mass and study the properties of winds both in the NLR and BLR traced by blueshifted/skewed [OIII] and CIV emission lines, respectively. These WISSH QSOs are powered by SMBH with masses $\ge$10$^9$ Msun accreting at 0.4<$\lambda_{Edd}$<3.1. We have found the existence of two sub-populations characterized by the presence of outflows at different distances from the SMBH. One population ([OIII] sources) exhibits powerful [OIII] outflows, rest-frame EW (REW) of the CIV emission REW$_{CIV}\approx$20-40 A and modest CIV velocity shift (v$_{CIV}^{peak}$) with respect to the systemic redshift (<=2000 km/s). The second population (Weak [OIII] sources), representing ~70% of the analyzed WISSH QSOs, shows weak/absent [OIII] emission and an extremely large v$_{CIV}^{peak}$ (up to ~8000 km/s and REW$_{CIV}$<=20 A). We propose two explanations for the observed behavior of the strength of the [OIII] emission in terms of orientation effects of the line of sight and ionization cone. The dichotomy in the presence of BLR and NLR winds could be likely due to inclination effects considering a polar geometry scenario for the BLR winds. We find a strong correlation with L$_{Bol}$ and an anti-correlation with $\alpha_{ox}$, whereby the higher L$_{Bol}$, the steeper $\alpha_{ox}$ and the larger is the v$_{CIV}^{peak}$. Finally, the observed dependence v$_{CIV}^{peak}\propto L_{Bol}^{0.28\pm0.04}$ is consistent with radiatively driven winds scenario, where strong UV continuum is necessary to launch the wind and a weakness of the X-ray emission is fundamental to prevent overionization of the wind itself., Comment: Accepted for publication in A&A

We perform a survey of the X-ray properties of 41 objects from the WISE/SDSS selected Hyper-luminous (WISSH) quasars sample, composed by 86 broad-line quasars (QSOs) with bolometric luminosity $L_{Bol}\geq 2\times 10^{47}\,erg\, s^{-1}$, at z~2-4. All but 3 QSOs show unabsorbed 2-10 keV luminosities $L_{2-10}\geq10^{45} \,erg \,s^{-1}$. Thanks to their extreme radiative output across the Mid-IR-to-X-ray range, WISSH QSOs offer the opportunity to significantly extend and validate the existing relations involving $L_{2-10}$. We study $L_{2-10}$ as a function of (i) X-ray-to-Optical (X/O) flux ratio, (ii) mid-IR luminosity ($L_{MIR}$), (iii) $L_{Bol}$ as well as (iv) $\alpha_{OX}$ vs. the 2500$\mathring{A}$ luminosity. We find that WISSH QSOs show very low X/O(<0.1) compared to typical AGN values; $L_{2-10}/L_{MIR}$ ratios significantly smaller than those derived for AGN with lower luminosity; large X-ray bolometric corrections $k_{\rm Bol,X}\sim$ 100-1000; and steep $-2<\alpha_{OX}<-1.7$. These results lead to a scenario where the X-ray emission of hyper-luminous quasars is relatively weaker compared to lower-luminosity AGN. Models predict that such an X-ray weakness can be relevant for the acceleration of powerful high-ionization emission line-driven winds, commonly detected in the UV spectra of WISSH QSOs, which can in turn perturb the X-ray corona and weaken its emission. Accordingly, hyper-luminous QSOs represent the ideal laboratory to study the link between the AGN energy output and wind acceleration. Additionally, WISSH QSOs show very large BH masses ($\log[M_{\rm BH}/M_{\odot}]$>9.5). This enables a more robust modeling of the $\Gamma-M_{BH}$ relation by increasing the statistics at high masses. We derive a flatter $\Gamma$ dependence than previously found over the broad range 5 <$\log(M_{\rm BH}/M_{\odot})$ < 11., Comment: 20 pages, 14 Figures. Accepted for publication on A&A

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

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., Comment: 19 pages, 12 figures; Accepted for publication in Astronomy & Astrophysics on June 13, 2017

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<=2, i.e. a strong form of co-evolution between SMBHs and galaxies appears to break down for the least massive galaxies., Comment: A&A in press

Models and observations suggest that both power and effects of AGN feedback should be maximised in hyper-luminous (L_Bol>10^47 erg/s) quasars, i.e. objects at the brightest end of the AGN luminosity function. We present the first results of a multi-wavelength observing program, focusing on a sample of WISE/SDSS selected hyper-luminous (WISSH) broad-line quasars at z~1.5-5. The WISSH quasars project has been designed to reveal the most energetic AGN-driven outflows, estimate their occurrence at the peak of quasar activity and extend the study of correlations between outflows and nuclear properties up to poorly-investigated extreme AGN luminosities (L_Bol~10^47 -10^48 erg/s). We present NIR, long-slit LBT/LUCI1 spectroscopy of five WISSH quasars at z~2.3-3.5 showing prominent [OIII] emission lines with broad (FWHM~1200-2200 km/s) and skewed profiles. The luminosities of the broad [OIII] wings are the highest measured so far (L_[OIII]^broad >~5x10^44 erg/s) and reveal the presence of powerful ionised outflows with mass outflow rates Mdot >~1700 M_Sun/yr and kinetic powers Edot >~10^45 erg/s. Although these estimates are affected by large uncertainties, due to the use of [OIII] as tracer of ionized outflows and the very basic outflow model we assume, these results suggest that the AGN is highly efficient in pushing outwards large amounts of ionised gas in our targets. The mechanical outflow luminosities for WISSH quasars correspond to higher fractions (~1-3%) of L_Bol than those derived for lower L_Bol AGN. Our targets host very massive (M_BH>~2x10^9 M_Sun) black holes which are still accreting at a high rate (i.e. a factor of ~0.4-3 of the Eddington limit). These findings demonstrate that WISSH quasars offer the opportunity of probing the extreme end of both luminosity and SMBH mass functions and revealing powerful ionised outflows able to affect the evolution of their host galaxies., Comment: 18 pages, 11 figures, A&A accepted