Search

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

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

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

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., Comment: Published in MNRAS. The source code is temporarily available at https://gitlab.lam.fr/gyang/cigale/tree/xray, but will be eventually integrated into CIGALE (https://cigale.lam.fr/) in its next release

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

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

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., Comment: 11 pages, 5 figures, accepted for publication in Astronomy & Astrophysics, Main Journal

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 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 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., Comment: 15 pages, 11 figures (including Appendix). Version matches A&A published version (after proof)

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

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