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.

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

CIGALE is a powerful multiwavelength spectral energy distribution (SED) fitting code for extragalactic studies. However, the current version of CIGALE is not able to fit X-ray data, which often provide unique insights into AGN intrinsic power. We develop a new X-ray module for CIGALE, allowing it to fit SEDs from the X-ray to infrared (IR). We also improve the AGN fitting of CIGALE from UV-to-IR wavelengths. We implement a modern clumpy two-phase torus model, SKIRTOR. To account for moderately extincted type 1 AGNs, we implement polar-dust extinction. We publicly release the source code (named X-CIGALE). We test X-CIGALE with X-ray detected AGNs in SDSS, COSMOS, and AKARI-NEP. The fitting quality (as indicated by reduced $\chi^2$) is good in general, indicating that X-CIGALE is capable of modelling the observed SED from X-ray to IR. We discuss constrainability and degeneracy of model parameters in the fitting of AKARI-NEP, for which excellent mid-IR photometric coverage is available. We also test fitting a sample of AKARI-NEP galaxies for which only X-ray upper limits are available from Chandra observations, and find that the upper limit can effectively constrain the AGN SED contribution for some systems. Finally, using X-CIGALE, we assess the ability of Athena to constrain the AGN activity in future extragalactic studies.

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.

The masses of supermassive black holes at the centres of local galaxies appear to be tightly correlated with the mass and velocity dispersions of their galactic hosts. However, the local Mbh-Mstar relation inferred from dynamically measured inactive black holes is up to an order-of-magnitude higher than some estimates from active black holes, and recent work suggests that this discrepancy arises from selection bias on the sample of dynamical black hole mass measurements. In this work we combine X-ray measurements of the mean black hole accretion luminosity as a function of stellar mass and redshift with empirical models of galaxy stellar mass growth, integrating over time to predict the evolving Mbh-Mstar relation. The implied relation is nearly independent of redshift, indicating that stellar and black hole masses grow, on average, at similar rates. Matching the de-biased local Mbh-Mstar relation requires a mean radiative efficiency ~0.15, in line with theoretical expectations for accretion onto spinning black holes. However, matching the "raw" observed relation for inactive black holes requires a mean radiative efficiency around 0.02, far below theoretical expectations. This result provides independent evidence for selection bias in dynamically estimated black hole masses, a conclusion that is robust to uncertainties in bolometric corrections, obscured active black hole fractions, and kinetic accretion efficiency. For our fiducial assumptions, they favour moderate-to-rapid spins of typical supermassive black holes, to achieve a mean radiative efficiency ~0.12-0.20. Our approach has similarities to the classic Soltan analysis, but by using galaxy-based data instead of integrated quantities we are able to focus on regimes where observational uncertainties are minimized.

The WISSH quasars project aims at studying the nuclear and host galaxy properties of the most luminous quasars ($L_{bol}>10^{47}$ erg/s, $1.8<z<4.6$). Nuclear winds are manifested as UV broad ($\geq$ 2,000 km/s) absorption lines (BAL) in $\sim$ 15\% of quasars. We aim at studying the incidence and properties of such winds in the WISSH sample, to investigate possible differences with respect to lower luminosity AGN regimes. We collected optical spectra from the SDSS data release 12, and identified those showing absorption troughs in the region between \siiv and \civ emission lines. We find a higher observed fraction of \civ BAL quasars in the WISSH sample (24\%) with respect to previous catalogues (10-15\%). These WISSH BAL quasars are also characterized by a larger average BI ($\sim$4,000 km/s) and maximum velocity ($\sim$17,000 km/s). Moreover, for two objects we discovered BAL features bluewards of the \siiv peak, which can be associated to \civ absorption with velocity of 0.15c. Finally, we updated previous studies on the dependence of maximum outflow velocity upon bolometric luminosity, showing that BAL winds have intermediate properties compared to molecular/ionized winds and ultra fast outflows (UFOs). Finally, the radio properties of the WISSH BAL quasars as a whole are in line with those of samples at lower luminosities from previous studies. Our results suggest that the higher $L_{\rm bol}$ of the WISSH quasars likely favours the acceleration of BAL outflows and that their most likely driving mechanism is radiation pressure. Furthermore we estimate that the kinetic power associated to these winds in hyperluminous quasars is sufficient, for the highest column density and fastest winds, to provide efficient feedback onto the host galaxy.

We present a 155ks NuSTAR observation of the $z\sim2$ hot dust-obscured galaxy (hot DOG) W1835+4355. We extracted spectra from the two NuSTAR detectors and analyzed them jointly with the archival XMM PN and MOS spectra. We performed a spectroscopic analysis based on both phenomenological and physically motivated models employing toroidal and spherical geometry for the obscurer. In all the modelings, the source exhibits a Compton-thick column density $N_{\rm H} \gtrsim 10^{24}$ cm$^{-2}$, a 2-10 keV luminosity $L_{2-10}\approx2\times10^{45}$ erg s$^{-1}$ , and a prominent soft excess ($\sim5-10$ % of the primary radiative output), which translates into a luminosity $\sim10^{44}$ erg s$^{-1}$. We modeled the spectral energy distribution from 1.6 to 850 $\mu m$ using a clumpy two-phase dusty torus model plus a modified blackbody to account for emission powered by star formation in the far-infrared. We employed several geometrical configurations consistent with those applied in the X-ray analysis. In all cases we obtained a bolometric luminosity $L_{\rm bol}\approx3-5\times10^{47}$ erg s$^{-1}$, which confirms the hyperluminous nature of this active galactic nucleus. Finally, we estimate a prodigious star formation rate of $\sim$3000 $M_{\odot}\,yr^{-1}$, which is consistent with the rates inferred for $z\approx2-4$ hyperluminous type I quasars. The heavily obscured nature, together with $L_{\rm bol}$, the ratio of X-ray to mid-infrared luminosity, the rest-frame optical morphology, and the host star formation rate are indicative of its evolutionary stage. We can interpret this as a late-stage merger event in the transitional, dust-enshrouded, evolutionary phase eventually leading to an optically bright AGN.

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.

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.

We imaged with ALMA and ARGOS/LUCI the molecular gas and the dust and stellar continuum in XID2028, an obscured QSO at z=1.593, where the presence of a massive outflow in the ionized gas component traced by the [O III]5007 emission has been resolved up to 10 kpc. This target represents a unique test case to study QSO 'feedback in action' at the peak epoch of AGN-galaxy coevolution. The QSO has been detected in the CO(5-4) transition and in the 1.3mm continuum, at ~30 and ~20 {\sigma} significance respectively, with both emissions confined in the central (<4 kpc) radius area. Our analysis suggests the presence of a fast rotating molecular disc (v~400 km/s) on very compact scales, and well inside the galaxy extent seen in the rest-frame optical light (~10 kpc, as inferred from the LUCI data). Adding available measurements in additional two CO transitions, CO(2-1) and CO(3-2), we could derive a total gas mass of ~10$^{10}$ M$_\odot$, thanks to a critical assessment of CO excitation and the comparison with Rayleigh-Jeans continuum estimate. This translates into a very low gas fraction (<5%) and depletion time scales of 40-75 Myr, reinforcing the result of atypical gas consumption conditions in XID2028, possibly due to feedback effects on the host galaxy. Finally, we also detect at ~5{\sigma} the presence of high velocity CO gas, which we interpret as a signature of galaxy-scale molecular outflow, spatially coincident with the ionised gas outflow. XID2028 represents therefore a unique case where the measurement of total outflowing mass (~500-800 M$_\odot$/yr) including the molecular and atomic components, in both the ionised and neutral phases, has been attempted for a high-z QSO.

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.

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.

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.

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.