Press Release
Free Access
Issue
A&A
Volume 635, March 2020
Article Number L5
Number of page(s) 10
Section Letters to the Editor
DOI https://doi.org/10.1051/0004-6361/201937292
Published online 10 March 2020

© ESO 2020

1. Introduction

The most luminous active galactic nuclei (AGN) are expected to exhibit the strongest and clear-cut manifestations of winds (Menci et al. 2008; Faucher-Giguère & Quataert 2012). Indeed, the fastest and most energetic winds have been reported in hyper-luminous quasars (i.e. with a bolometric luminosity of Lbol ≳ 1047 erg s−1; Wu et al. 2011; Fiore et al. 2017; Vietri et al. 2018; Meyer et al. 2019; Perrotta et al. 2019). Luminous quasars are typically characterised by their low coronal X-ray luminosity (LX) compared to the disc’s ultraviolet (UV) and larger-scale dust-reprocessed mid-infrared (MIR) luminosities, which are represented as LUV and LUV, respectively. This characterisation is usually performed by adopting the αOX1 and LMIR/LX parameters (e.g. Vignali et al. 2003; Just et al. 2007; Lusso & Risaliti 2016; Stern 2015; Martocchia et al. 2017; Chen et al. 2017). Past studies on large quasar samples over a wide luminosity range (Lbol ≈ 1045 − 1048 erg s−1) have also reported indications of a further weak dependence of the αOX parameter on the velocities of the broad emission line (BEL) winds, which are parametrised by the CIV blueshift, that is, negative velocity shifts (vCIV) of the CIV emission line (e.g. Richards et al. 2011; Kruczek et al. 2011; Vietri et al. 2018). However, once the αOX luminosity dependence is removed by adopting ΔαOX, that is, the difference between the measured αOX and the one predicted by the αOXL2500 Å relation (e.g. Just et al. 2007), the dependence with vCIV is much less significant (e.g. Gibson et al. 2008; Ni et al. 2018; Timlin et al. 2020). It is important to notice that these dependences on vCIV may be regarded as representative for the bulk of the quasar population since CIV BEL winds are a common feature in quasars (e.g. Shen et al. 2011). A similar dependence has also been reported for the maximum velocity measured in line-of-sight detected winds, such as those reported in the broad-absorption line (BAL) quasars (e.g. Gallagher et al. 2006; Gibson et al. 2009; Stalin et al. 2011). Nonetheless, a peculiar category of z = 2 − 3 optically luminous quasars selected to have weak broad emission lines (with a rest-frame equivalent width of the CIV, REWCIV < 15 Å) and fast BEL winds (vCIV ≲ −2000 km s−1 Richards et al. 2011), has mostly revealed 1–2 dex weaker X-ray emission (Gibson et al. 2008; Wu et al. 2011) compared to the αOXL2500 Å expectations.

Sources with vCIV ≲ −2000 km s−1 and REWCIV ≲ 20 Å are reported in our sample of 86 broad-line unlensed highly accreting (λEdd >  0.4) MIR and optically bright z = 2 − 4 hyper-luminous quasars. These sources were selected to be the MIR-brightest WISE/SDSS quasars with z >  1.5 and a flux density of S22μm >  3 mJy (WISSH; Bischetti et al. 2017). The WISSH quasars exhibit widespread evidence of winds at all scales from nuclear BAL (Bruni et al. 2019) to ionised [OIII]/Ly α galactic- and circumgalactic-scale outflows (Bischetti et al. 2017; Travascio et al. 2020). In particular, in Vietri et al. (2018) we report a surprisingly high fraction (∼70%) of sources with weak UV and optical BEL spectra (e.g. REWCIV ≲ 20 Å) as well as extreme CIV blueshifts (vCIV < −2000 km s−1) in a sub-sample of the WISSH quasars. Furthermore, Martocchia et al. (2017) found a large spread in the LX with a non-negligible fraction that have ∼0.5–1 dex fainter values than the average.

In this Letter, we explore the relation between the extreme radiative field of the hyper-luminous WISSH quasars and their X-ray coronal properties, and we interpret it in the context of accretion disc wind scenarios. We adopt a ΛCDM cosmology with ΩΛ = 0.73 and H0 = 70 km s−1 Mpc−1 throughout the paper. Errors as well as upper and lower limits are quoted at 68.3% and 90% confidence levels, respectively.

2. Sample presentation and X-ray data reduction and analysis

In this work we consider the radio-quiet hyper-luminous WISSH sources with (i) reported vCIV measures that are relative to their systemic redshift and (ii) available X-ray data. The selected sample of thirteen sources is reported in Table A.1 and was mainly drawn from the Vietri et al. (2018) WISSH sub-sample of 18 quasars for which CIV emission line properties were derived. We also complement our work by including five Type 1 radio-quiet hyper-luminous sources at similar z with published vCIV (Coatman et al. 2017, hereafter C17; see Table A.1) and available X-ray archived data. Further details on the sample selection are reported in Appendix A.

We consider both Chandra and XMM-Newton observations that are available for each source (see Table B.1). For each dataset, we performed standard data reduction as detailed in Appendix B. The X-ray spectral modelling was performed in the 0.2–10 keV and 0.3–8 keV bands for XMM-Newton and Chandra, respectively. We employed an intrinsically absorbed power-law model, which was further modified by the Galactic absorption. Further details on the modelling and the derived parameters are reported in Appendix C and Table C.1.

3. Results

Figure 1 shows the unabsorbed (i.e. intrinsic) 2–10 keV luminosity (L2 − 10; left panel) and the X-ray unabsorbed and UV de-extincted αOX (middle panel) as a function of vCIV for the WISSH quasars (black). Both quantities strongly correlate with vCIV2 with the Spearman correlation coefficient ρ ≈ 0.6 and the two-sided null-hypothesis probability of p ≈ 0.02 (see Table 1). This is a consequence of the lack of significant correlation between 2500Å and vCIV (right panel; ρ = 0.32 and p ≈ 0.29) and the limited dispersion of 2500Å. Hence the sources with the largest negative vCIV (i.e. larger blueshifts) are 0.5–1 dex X-ray weaker and exhibit steeper αOX than the sources with the lowest vCIV.

thumbnail Fig. 1.

Left panel: L2 − 10 as a function of vCIV for the WISSH quasars (black circles). We also report hyper-luminous quasars at z = 2 − 3 (C17 sample, blue triangles). Middle panel: intrinsic αOX as a function of vCIV for the WISSH quasars. Red circles indicate X-ray weak sources (see Sect. 3 for details). Right panel: L2500 Å as a function of vCIV. Solid and dashed lines report the best-fit linear relation and 1σ uncertainties on the normalisation for the WISSH-only relations. Dotted line in the left panel reports the best-fit relation with the addition of the C17 sample. Orange lines in the right panel report the best-fit L2 − 10vCIV relation normalised to the mean values of L2500 Å.

Table 1.

Correlation and linear regression coefficients.

No significant correlation between vCIV and Lbol and L6μm is found (see left panels in Fig. 2 and Table 1). Therefore their ratio with L2 − 10, that is, the X-ray bolometric correction kbol, X = Lbol/L2 − 10 and L6μm/L2 − 10, are strongly dependent on vCIV (see right panels in Fig. 2 and Table 1). The inclusion of the hyper-luminous quasars from C17 in the relations involving L2 − 10, Lbol, and kbol, X further confirms the strength and significance of the dependence with vCIV. We note that Vietri et al. (2018) and C17 adopted slightly different definitions of vCIV (see notes on Table A.1). This difference does not change our result. Indeed, for the two WISSH quasars reported in both samples (i.e. J1106+6400 and J1201+1206), the vCIV reported by C17 are 400 − 500 km s−1 larger than the values reported by Vietri et al. (2018). This small systematic offset makes little difference in our correlations and if we correct the vCIV of the C17 sub-sample by 500 km s−1, we obtain a slightly stronger L2−10vCIV correlation with ρ = 0.63 and p = 0.0053.

thumbnail Fig. 2.

Top panels: log L6μm (left) and L6μm/L2 − 10 (right) as a function of vCIV. Bottom panels: log Lbol (left) and kbol, X (right) as a function of vCIV. The meaning of the data points and lines is the same as in Fig. 1.

The ΔαOX values (based on the αOX − 2500Å relation from Just et al. 2007) appear to be strongly dependent on vCIV (see Table 1). Similar to past works, we also calculated , that is, we derived it by performing spectral fitting with an unabsorbed power-law model for the 2 keV luminosity estimation (see Appendix C) and the relative (see Table C.1). Both quantities exhibit a strong and significant correlation with vCIV, which is similar to αOX and ΔαOX. We find that four and five sources (∼30 − 40% of the considered WISSH sub-sample) are X-ray weak (red circles in Fig. 1), as they have ΔαOX and values that are less than −0.2 (a typically adopted threshold value, e.g. Luo et al. (2015), used to identify X-ray weak sources), respectively. We notice that, given the small number of sources and the sample selection function (i.e. only X-ray detected sources for which basic X-ray spectral analysis can be performed), the fraction of X-ray weak sources must be considered with caution, although see Nardini et al. (2019) for similar results.

For all of the significant relations (i.e. with p ≲ 0.02), we also computed ρ and account for the uncertainties in the measurements. We generated 10 000 random realisations of the sample by Gaussian distributing each quantity that was considered according to its best-fit value and error. We obtain ρsim and its uncertainty by adopting the mean value and standard deviation of the distribution of ρ for the simulated datasets. We also calculated psim, 90, the 90% percentile on the distribution of the p values. The derived ρsim and psim, 90 confirm the significance of the relations. For all quantities, we derived linear relations with vCIV by employing the BCES(Y|X) method (Akritas & Bershady 1996). We show them in Figs. 1 and 2. Table 1 reports ρ, p, ρsim, psim, 90, and the slope (α) and y-intercept (β) of the linear relation.

As for the derived X-ray spectral parameters, we do not report a significant dependence of the photon index (Γ) on vCIV (see Table 1), but it is important to bear in mind that we have weak constraints on Γ for the X-ray weaker sources (see top panel of Fig. 3). As for the column density (NH), we find the WISSH quasars to exhibit moderate values of absorption which are compatible with their Type 1 nature. Indeed seven sources have measured NH ≲ 1023 cm−2 with a mean log(NH/cm−2)≈22.3. The other sources have upper limits in the range of log(NH/cm−2) = 21 − 23.7, with the largest values mainly driven by poor statistics in the spectra (≲50 net-counts). The lower panel of Fig. 3 reports NH as a function of vCIV. We computed the Spearman’s rank correlation coefficient, which accounts for upper limits4. We find a weak but not significant anti-correlation both including or excluding the C17 data (see Table 1).

thumbnail Fig. 3.

Γ and NH as a function of vCIV (top and bottom panels). The meaning of the data points and linear fits is the same as in Fig. 1.

Finally we mention a lack of correlation between λEdd (as reported in Vietri et al. 2018) and vCIV for the WISSH quasars. Indeed, all of these sources shine at λEdd ≈ 1 and the reported large uncertainties on the λEdd, which are mainly driven by the constraints on the supermassive black-hole (SMBH) mass (see Vietri et al. 2018), prevent an accurate investigation of a possible trend.

4. Discussion

We report a relation between L2 − 10 and αOX with vCIV in a sample of MIR and optically-selected hyper-luminous quasars, and we also note the lack of a similar dependence in the UV, MIR, and in Lbol. The use of good-quality X-ray data and the well-defined quasar selection, resulting in a narrow Lbol (≈ LUV) range, allows one to measure a marked αOX-vCIVαOX-vCIV) dependence (i.e. stronger in terms of the correlation coefficient than those reported to date in samples probing the bulk of the quasar population; Richards et al. 2011; Vietri et al. 2018; Timlin et al. 2020) and reveals a clear-cut dependence on L2 − 10. This suggests that at these luminosity regimes, L2 − 10 is the main driver of the αOX-vCIV relation. We notice that the large average ⟨vCIV⟩ for our quasars of ∼ − 2900 km s−1 is in agreement with the increasing trend of the relation reported, for example in Vietri et al. (2018), for an extended sample spanning more than three decades in Lbol (see also Timlin et al. 2020 for a similar result). Accordingly, the lack of a significant correlation for log Lbol-vCIV in our sample is the result of the restricted luminosity range spanned by our quasars.

Theoretical and observational arguments suggest that BEL winds in AGN are produced at accretion disc scales (e.g. Elvis 2000), where the AGN luminosity output is large in the UV and X-ray bands. Accretion disc winds in AGN can be sustained by either magnetic or radiative forces; the observational results of this Letter imply that regardless of the driving mechanism of the disc wind, the fastest UV winds appear in the X-ray weakest sources.

In magnetohydrodynamic (MHD)-driven scenarios, the presence of the disc wind does not depend on the X-ray radiative output; however the existence of observable ions does. Indeed, a general correlation between αOX and, for example, vCIV is expected to exist from simple photoionisation arguments. A larger X-ray luminosity is generally effective in stripping the bound electrons off the CIV atoms up to large disc radii, allowing for the formation of low-velocity winds. Lower X-ray luminosities instead allow for the presence of CIV at distances closer to the SMBH where larger terminal velocities are expected. In particular, in order to not get over-ionised, the UV wind must be accompanied by an inner shield of partially ionised gas absorbing the X-ray flux, which is itself part of the wind driven by the magnetic forces (see e.g. Fukumura et al. 2010).

In radiation-driven scenarios, the wind is driven by radiation pressure on spectral lines, and the relative contribution of the UV and X-ray emission is instead very important in determining the existence of the wind itself (Murray et al. 1995; Proga et al. 2000; Proga & Kallman 2004). In particular, for a given UV luminosity, a weak X-ray emission (i.e. steep αOX) is crucial for the existence of fast radiation-driven winds. Indeed, in the inner disc regions, the gas opacity to UV transitions drops as the atoms are over-ionised by an intense X-ray flux. In this case, the wind cannot be efficiently accelerated beyond escape velocity and eventually it falls back to the disc as a failed wind (FW). The FW contains clumps of dense gas in the proximity of the X-ray emitting corona. The FW effectively shields the gas that is located farther away from the ionising X-ray photons, and it allows for the acceleration of disc material at all radii where the radiation pressure is high enough to overcome the gravitational pull of the SMBH (Proga et al. 2000; Proga & Kallman 2004; Risaliti & Elvis 2010).

A high X-ray flux would produce a vast inner zone of FW and would allow for the launch of disc winds only at large radii. Conversely, a lower X-ray flux would produce a reduced inner FW zone, and would allow for the formation of disc winds on scales that are closer to the SMBH. As the terminal velocity of the wind is inversely proportional to its launching distance from the central SMBH, a lower X-ray emission favours, in general, the launch of faster radiation-driven accretion disc winds, compared to a higher X-ray emission (see e.g. Giustini & Proga 2019, for a recent review).

Interestingly, because of their high density, FW clumps can further cause an efficient cooling of the corona via bremsstrahlung emission. This, therefore, leads to a weakening (quenching) of the inverse Compton X-ray emission (Proga 2005; Laor & Davis 2014).

We notice that recent post-processing radiative transfer calculations suggest that the FW is not able to efficiently prevent over-ionisation and, therefore, it may not be so crucial for wind acceleration. Indeed, the FW is found to have a higher ionisation state (than previous estimates) and an extremely limited X-ray shielding power (Higginbottom et al. 2014). Evidence against the X-ray shielding scenario comes from X-ray observations of semi-relativistic BAL winds, reporting weak and moderate X-ray absorption (∼0.5 − 5 × 1022 cm−2; Hamann et al. 2013). Similar evidence also seems to be present for hyper-luminous quasars. Indeed, our result would support a high degree of intrinsic X-ray quenching (Fig. 1, left panel) rather than a dependence on nuclear X-ray shielding (Fig. 3, lower panel). Further support for the hypothesis of coronal X-ray weakness comes from NuSTAR estimates that at least ∼1/3 of luminous BAL quasars may exhibit significant X-ray weakness (Luo et al. 2014). However, the Luo et al. (2014) sample consists of heavily obscured (NH ≈ 1024 cm−2) BAL quasars, which, in principle, may lead to part of the X-ray emission being further suppressed by nuclear shielding.

Studies focusing on samples of weak-line quasars explain their properties in the context of simple orientation-dependent nuclear shielding of the X-ray emission without invoking coronal quenching (Wu et al. 2011; Luo et al. 2015; Ni et al. 2018). In such a model, the shield may likely be produced by the geometrically thick inner accretion disc regions, and it would lead to a significant dependence on NH and a lack of correlation with the intrinsic LX. In our MIR and optically selected WISSH quasars, we instead find a marked L2 − 10-vCIV dependence and a lack of correlation with NH. However, it is important to notice that only simple cold absorbers have been adopted in modelling their X-ray spectrum thus far. Probably the inner shielding disc regions would require a more complex and ionised absorber. In this sense, the findings by Hamann et al. (2013) on semi-relativistic BAL winds in sources with low and moderate observed cold NH may imply the existence of such an ionised absorber which would act as an ionisation shield. The current X-ray data quality for our sample is not sufficient enough to add further free parameters to account for ionised absorption.

We mention that orientation may also play a role in producing the large scatter in the L2 − 10-vCIV relation. Indeed it may reflect the projection of the wind velocity field depending on the line-of-sight inclination of the disc-corona structure. For instance, in the scenario envisaged by Elvis (2000), which predicts an extreme funnel-shaped geometry of the wind for luminous quasars (see Fig. 7 in Elvis 2000), we qualitatively expect that at a given L2 − 10 the highest emission line blueshifts are reported by sources seen at relatively large inclinations (compatible with their Type 1 nature). Whereas, the lowest blueshifts are seen in sources viewed pole-on.

X-ray observations of well-monitored ultra-fast outflows (UFOs; e.g. Tombesi et al. 2010; Gofford et al. 2013) seem to support the important role of radiation in driving AGN disc winds. For example, the hyper-luminous local quasar PDS 456, known to display recurrent and variable UFOs (e.g. Reeves et al. 2009; Nardini et al. 2015), fits well in the L2 − 10-vCIV relation as it is reported to have a CIV blueshift of ∼5200 km s−1 (O’Brien et al. 2005) and L2 − 10 ≈ (0.3 − 1) × 1045 erg s−1. Interestingly, Matzeu et al. (2017) report a positively correlated variability on PDS 456 (over a period of 13 years) between the UFO velocity and the X-ray luminosity5. If PDS 456 is representative of the high-redshift hyper-luminous quasars, this level of X-ray variability could account for the scatter in the relation. As for the opposite signs of correlation for CIV BEL winds and UFO velocities with L2 − 10, if proven to be a common characteristics in luminous quasars, they may as well be qualitatively explained in the context of the radiation-driven accretion disc wind scenario. Indeed, the X-ray luminosity, which acts on X-ray line transitions, is responsible for the radiative acceleration of the UFO in one case. In the other case, it acts as an ionisation state regulator of UV line-driven BEL winds. We notice though that recent photoionisation and radiative transfer calculations by Dannen et al. (2019) suggest that the line driving mechanisms may not be relevant in plasma with high ionisation parameters typical of UFOs (i.e. ξ >  1000; e.g. Tombesi et al. 2011; Nardini et al. 2015). In this case, MHD-driven winds may be a viable mechanism to explain the positive correlation exhibited by UFOs as recently reported by Fukumura et al. (2018).

Dedicated deep X-ray observations, which will also be performed at lower luminosity regimes, will be crucial in testing and shedding light on the origin of the vCIV dependence by better constraining the spectral parameter for the X-ray weak sources. In this regard, ATHENA will further enable us to investigate the properties of the accretion disc-scale absorber (e.g. kinematics and ionisation state), discriminate between competing disc wind scenarios, and further investigate the role of UFOs (e.g. Martocchia et al. 2017; Barret & Cappi 2019).

Future studies on αOX, kbol, X, and L6μm/L2 − 10 will need to consider the reported marked dependences on vCIV (Figs. 1 and 2) at these high luminosity regimes. This will be a necessary step in order to obtain better constrained relations and remove possible systematics due to the inclusion of highly blueshifted and hence X-ray weak sources.


1

αOX is the X-ray to optical spectral index between the rest-frame luminosities at 2500 Å and 2 keV; i.e. αOX = 0.3838 log(L2 keV /L2500) Å.

2

See Table A.1 for its definition.

3

We mention that the relation between L2 − 10 and the CIV BEL wind velocity is still in place (ρ ≈ 0.77 and p = 0.003) even if we replace vCIV, the peak of the CIV emission line, with the velocity of the CIV outflow component as estimated in the two component CIV spectral modelling by Vietri et al. (2018).

4

We used the ASURV package v. 1.2 to account for censored data (Lavalley et al. 1992; Feigelson & Nelson 1985; Isobe et al. 1986).

5

We notice that variable UFOs have also been observed in other four high-z bright quasars (APM 08279+5255, PG 1115+080, H 1413+117 and HS 1700+641; Chartas et al. 2003, 2007; Saez & Chartas 2011; Lanzuisi et al. 2012). A similar correlation between the quasar X-ray luminosity and the UFO velocity has been reported for APM 08279+5255 (Saez & Chartas 2011).

6

In order to derive the K-corrected rest-frame 6 cm flux densities from the 1.4 GHz observed ones, we assumed a power-law synchrotron spectrum S(ν)∝να with α = 0.5 (e.g. Jiang et al. 2007).

7

Compared to Martocchia et al. (2017) who analysed pn-only spectra, we further improved the statistics of the XMM-Newton spectra by adding the MOS data.

8

The MOS detectors have been reported to have at most 7% higher fluxes than the pn as detailed in the XMM-SOC-CAL-TN-0052 Issue 6.0 available in http://xmm2.esac.esa.int/docs/documents/CAL-TN-0052.ps.gz

Acknowledgments

We thank the referee for his/her comments and suggestions. We thank Silvia Martocchia for useful discussions. LZ, EP, AB and MB acknowledge financial support under ASI/INAF contract 2017-14-H.0. CR acknowledges support from the CONICYT+PAI Convocatoria Nacional subvencion a instalacion en la academia convocatoria año 2017 PAI77170080. MG is supported by the “Programa de Atracción de Talento” of the Comunidad de Madrid, grant number 2018-T1/TIC-11733. GM is supported by the Spanish State Research Agency (AEI) Project No. ESP-2017-86582-C4-1-R. This research has been partially funded by the AEI Project No. MDM-2017-0737 Unidad de Excelencia “María de Maeztu” - Centro de Astrobiología (INTA-CSIC).

References

  1. Akritas, M. G., & Bershady, M. A. 1996, ApJ, 470, 706 [NASA ADS] [CrossRef] [Google Scholar]
  2. Barret, D., & Cappi, M. 2019, A&A, 628, A5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Bischetti, M., Piconcelli, E., Vietri, G., et al. 2017, A&A, 598, A122 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  4. Bruni, G., Piconcelli, E., Misawa, T., et al. 2019, A&A, 630, A111 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Cash, W. 1979, ApJ, 228, 939 [NASA ADS] [CrossRef] [Google Scholar]
  6. Chartas, G., Brandt, W. N., & Gallagher, S. C. 2003, ApJ, 595, 85 [NASA ADS] [CrossRef] [Google Scholar]
  7. Chartas, G., Eracleous, M., Dai, X., Agol, E., & Gallagher, S. 2007, ApJ, 661, 678 [NASA ADS] [CrossRef] [Google Scholar]
  8. Chen, C.-T. J., Hickox, R. C., Goulding, A. D., et al. 2017, ApJ, 837, 145 [NASA ADS] [CrossRef] [Google Scholar]
  9. Coatman, L., Hewett, P. C., Banerji, M., et al. 2017, MNRAS, 465, 2120 [NASA ADS] [CrossRef] [Google Scholar]
  10. Dannen, R. C., Proga, D., Kallman, T. R., & Waters, T. 2019, ApJ, 882, 99 [Google Scholar]
  11. Duras, F., Bongiorno, A., Piconcelli, E., et al. 2017, A&A, 604, A67 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Elvis, M. 2000, ApJ, 545, 63 [NASA ADS] [CrossRef] [Google Scholar]
  13. Faucher-Giguère, C.-A., & Quataert, E. 2012, MNRAS, 425, 605 [NASA ADS] [CrossRef] [Google Scholar]
  14. Feigelson, E. D., & Nelson, P. I. 1985, ApJ, 293, 192 [NASA ADS] [CrossRef] [Google Scholar]
  15. Fiore, F., Feruglio, C., Shankar, F., et al. 2017, A&A, 601, A143 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Fukumura, K., Kazanas, D., Contopoulos, I., & Behar, E. 2010, ApJ, 723, L228 [NASA ADS] [CrossRef] [Google Scholar]
  17. Fukumura, K., Kazanas, D., Shrader, C., et al. 2018, ApJ, 864, L27 [NASA ADS] [CrossRef] [Google Scholar]
  18. Gallagher, S. C., Brandt, W. N., Chartas, G., et al. 2006, ApJ, 644, 709 [NASA ADS] [CrossRef] [Google Scholar]
  19. Gibson, R. R., Brandt, W. N., & Schneider, D. P. 2008, ApJ, 685, 773 [NASA ADS] [CrossRef] [Google Scholar]
  20. Gibson, R. R., Brandt, W. N., Gallagher, S. C., & Schneider, D. P. 2009, ApJ, 696, 924 [NASA ADS] [CrossRef] [Google Scholar]
  21. Giustini, M., & Proga, D. 2019, A&A, 630, A94 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Gofford, J., Reeves, J. N., Tombesi, F., et al. 2013, MNRAS, 430, 60 [Google Scholar]
  23. Hamann, F., Chartas, G., McGraw, S., et al. 2013, MNRAS, 435, 133 [NASA ADS] [CrossRef] [Google Scholar]
  24. Higginbottom, N., Proga, D., Knigge, C., et al. 2014, ApJ, 789, 19 [Google Scholar]
  25. Isobe, T., Feigelson, E. D., & Nelson, P. I. 1986, ApJ, 306, 490 [NASA ADS] [CrossRef] [Google Scholar]
  26. Jiang, L., Fan, X., Ivezić, Ž., et al. 2007, ApJ, 656, 680 [NASA ADS] [CrossRef] [Google Scholar]
  27. Just, D. W., Brandt, W. N., Shemmer, O., et al. 2007, ApJ, 665, 1004 [Google Scholar]
  28. Kellermann, K. I., Sramek, R., Schmidt, M., Shaffer, D. B., & Green, R. 1989, AJ, 98, 1195 [NASA ADS] [CrossRef] [Google Scholar]
  29. Kruczek, N. E., Richards, G. T., Gallagher, S. C., et al. 2011, AJ, 142, 130 [NASA ADS] [CrossRef] [Google Scholar]
  30. Lanzuisi, G., Giustini, M., Cappi, M., et al. 2012, A&A, 544, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  31. Laor, A., & Davis, S. W. 2014, MNRAS, 438, 3024 [NASA ADS] [CrossRef] [Google Scholar]
  32. Lavalley, M. P., Isobe, T., & Feigelson, E. D. 1992, BAAS, 24, 839 [Google Scholar]
  33. Luo, B., Brandt, W. N., Alexander, D. M., et al. 2014, ApJ, 794, 70 [NASA ADS] [CrossRef] [Google Scholar]
  34. Luo, B., Brandt, W. N., Hall, P. B., et al. 2015, ApJ, 805, 122 [NASA ADS] [CrossRef] [Google Scholar]
  35. Lusso, E., & Risaliti, G. 2016, ApJ, 819, 154 [NASA ADS] [CrossRef] [Google Scholar]
  36. Martocchia, S., Piconcelli, E., Zappacosta, L., et al. 2017, A&A, 608, A51 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  37. Matzeu, G. A., Reeves, J. N., Braito, V., et al. 2017, MNRAS, 472, L15 [NASA ADS] [CrossRef] [Google Scholar]
  38. Menci, N., Fiore, F., Puccetti, S., & Cavaliere, A. 2008, ApJ, 686, 219 [NASA ADS] [CrossRef] [Google Scholar]
  39. Meyer, R. A., Bosman, S. E. I., & Ellis, R. S. 2019, MNRAS, 487, 3305 [NASA ADS] [CrossRef] [Google Scholar]
  40. Murray, N., Chiang, J., Grossman, S. A., & Voit, G. M. 1995, ApJ, 451, 498 [NASA ADS] [CrossRef] [Google Scholar]
  41. Nardini, E., Reeves, J. N., Gofford, J., et al. 2015, Science, 347, 860 [Google Scholar]
  42. Nardini, E., Lusso, E., Risaliti, G., et al. 2019, A&A, 632, A109 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Ni, Q., Brandt, W. N., Luo, B., et al. 2018, MNRAS, 480, 5184 [NASA ADS] [Google Scholar]
  44. O’Brien, P. T., Reeves, J. N., Simpson, C., & Ward, M. J. 2005, MNRAS, 360, L25 [NASA ADS] [CrossRef] [Google Scholar]
  45. Perrotta, S., Hamann, F., Zakamska, N. L., et al. 2019, MNRAS, 488, 4126 [NASA ADS] [CrossRef] [Google Scholar]
  46. Piconcelli, E., Vignali, C., Bianchi, S., et al. 2015, A&A, 574, L9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  47. Proga, D. 2005, ApJ, 630, L9 [NASA ADS] [CrossRef] [Google Scholar]
  48. Proga, D., & Kallman, T. R. 2004, ApJ, 616, 688 [Google Scholar]
  49. Proga, D., Stone, J. M., & Kallman, T. R. 2000, ApJ, 543, 686 [NASA ADS] [CrossRef] [Google Scholar]
  50. Reeves, J. N., O’Brien, P. T., Braito, V., et al. 2009, ApJ, 701, 493 [NASA ADS] [CrossRef] [Google Scholar]
  51. Richards, G. T., Kruczek, N. E., Gallagher, S. C., et al. 2011, AJ, 141, 167 [NASA ADS] [CrossRef] [Google Scholar]
  52. Risaliti, G., & Elvis, M. 2010, A&A, 516, A89 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  53. Runnoe, J. C., Brotherton, M. S., & Shang, Z. 2012, MNRAS, 422, 478 [NASA ADS] [CrossRef] [Google Scholar]
  54. Saez, C., & Chartas, G. 2011, ApJ, 737, 91 [Google Scholar]
  55. Shen, Y., Richards, G. T., Strauss, M. A., et al. 2011, ApJS, 194, 45 [NASA ADS] [CrossRef] [Google Scholar]
  56. Stalin, C. S., Srianand, R., & Petitjean, P. 2011, MNRAS, 413, 1013 [NASA ADS] [CrossRef] [Google Scholar]
  57. Stern, D. 2015, ApJ, 807, 129 [NASA ADS] [CrossRef] [Google Scholar]
  58. Timlin, J. D., Brandt, W. N., Ni, Q., et al. 2020, MNRAS, 492, 719 [NASA ADS] [CrossRef] [Google Scholar]
  59. Tombesi, F., Cappi, M., Reeves, J. N., et al. 2010, A&A, 521, A57 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  60. Tombesi, F., Cappi, M., Reeves, J. N., et al. 2011, ApJ, 742, 44 [Google Scholar]
  61. Travascio, A., Zappacosta, L., Cantalupo, S., et al. 2020, A&A, in press https://doi.org/10.1051/0004-6361/201936197 [Google Scholar]
  62. Vietri, G., Piconcelli, E., Bischetti, M., et al. 2018, A&A, 617, A81 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  63. Vignali, C., Brandt, W. N., & Schneider, D. P. 2003, AJ, 125, 433 [NASA ADS] [CrossRef] [Google Scholar]
  64. Wachter, K., Leach, R., & Kellogg, E. 1979, ApJ, 230, 274 [Google Scholar]
  65. Wu, J., Brandt, W. N., Hall, P. B., et al. 2011, ApJ, 736, 28 [NASA ADS] [CrossRef] [Google Scholar]

Appendix A: Sample selection

The sources considered in this work are reported in Table A.1. They were selected starting from the 18 WISSH quasars analysed by Vietri et al. (2018) for which CIV emission line properties were estimated. We only include the twelve radio-quiet sources with available archived X-ray Chandra and XMM-Newton data (see Table B.1). We added an additional WISSH source (J1441+0454) to this sample with available X-ray archived data for which constraints for the CIV emission line have already been obtained and will be reported in Vietri et al. in preparation. In Fig. A.1 we report the LUV − z for our sub-sample (red points) together with the entire WISSH sample. The sources in our sub-sample are mainly clustered at z ≈ 2.1 and z ≈ 3.4; therefore, they encompass the large z range of the WISSH sources and exhibit LUV values that are representative of the entire WISSH sample. X-ray data for eleven sources have already been analysed by Martocchia et al. (2017). In this work, additional and new proprietary and publicly available X-ray data are considered for seven sources (J0958+2827, J1201+0116, J1236+6554, J1326-0005, J0900+4215, J1106+6400, J2123-0050). Two of these sources (J0958+2827, J1326-0005) are not in the Martocchia et al. (2017) sample. The final WISSH sample studied here consists of thirteen sources. Chandra and XMM-Newton X-ray data are available for eleven and six sources, respectively. Three sources have both Chandra and XMM-Newton coverage. The log of the X-ray observations for these sources is reported in Table B.1.

thumbnail Fig. A.1.

L6μm versus z for our WISSH sub-sample (filled red circles) and the entire WISSH sample (empty black circles).

Table A.1.

Properties of the WISSH quasars and the additional sources from the C17 sample.

In order to increase the number of hyper-luminous sources with available X-ray data and CIV spectral properties, we consider an additional sample starting from the sources analysed in C17 with Lbol >  1047 erg s−1 and available pointed X-ray data. The C17 sample is a heterogeneous compilation of 230 quasars with available infrared and optical spectroscopy. Our final sample consists of the five sources reported in Table A.1. Chandra and XMM-Newton data are available for two and three sources, respectively (see Table B.1).

Finally, the sources that are considered are all undetected in the FIRST/NVSS radio surveys with the only exception of three sources (J0900+4215, J1441+0454, J1549+1245) with 1.4 GHz integrated flux densities of ∼1.7 mJy. From the integrated FIRST flux-densities at 1.4 GHz, we calculated the radio loudness parameter R defined as the ratio between the rest-frame 6 cm and 2500 Å flux densities6. We obtained R <  10 for all the sources and therefore we conclude that they are not radio-loud (Kellermann et al. 1989).

Appendix B: X-ray data reduction

Data from sources listed in Table B.1 whose Observation ID (OBSID) is marked by an asterisk have been also analysed in Martocchia et al. (2017). We performed the reduction of the new Chandra data with CIAO 4.9 (with version 4.7.3 of the Chandra Calibration Database). We reprocessed each data set with the script CHANDRA_REPRO with the default script settings in order to generate level 2 data with the most updated calibration. We inspected the light-curves on each observation in order to check for possible high background periods and did not find any observation to be affected by such contaminations. The spectral extraction and the creation of the relative response files were performed using the CIAO script SPECEXTRACT. We used source extraction circular regions with radii in the range of 1.5–3 arcsec in order to include all of the source counts. For the background, we adopted annular source-free regions centred on the quasar with inner and outer radii of 6 and 30 arcsec, respectively. For the source j0958, we reduced two observations and then added spectral and response files with the FTOOLS script ADDASCASPEC.

Table B.1.

Log of the X-ray observations of the sources from the WISSH and C17 samples.

The reduction of the XMM-Newton data (both pn and MOS cameras7) was performed with SAS v16.0.0. All of the observations were performed in full-window mode with the thin filter applied. The light curves for pn and MOS exposures were screened at energies > 10 keV (10 − 12 keV for the pn) for high background flaring periods. We adopted a count-rate threshold filtering criterion with typical values of 0.3–0.5 and 0.1–0.2 counts s−1 for pn and MOS, respectively. The resulting net-exposure times are reported in Table B.1. We selected X-ray events corresponding to patterns 0–4 and 0–12 for pn and MOS, respectively. The source extraction was performed using the same circular apertures for both pn and MOS detectors. In order to include all of the source counts and simultaneously minimise the background contribution for each source, we adopted different apertures ranging from 12 to 30 arcsec. The background spectrum was extracted in the chip including the source from circular (1–2 arcmin radius) and annular (inner and outer radii up to 0.7 and 4 arcmin) source-free regions for pn and MOS, respectively. The same data reduction steps were followed in case of the sources belonging to the C17 sample (see Table B.1 for details on their X-ray observations).

For Chandra spectra, we obtained background-subtracted counts ranging from 24 to ∼200 in the 0.3–8 keV band, reaching ∼800 counts in one case (J2123−0050). As for XMM-Newton, we obtained spectra with background-subtracted counts ranging from 70 (126) to 1700 (1600) for pn (MOS1+MOS2) detectors. We consistently grouped all Chandra and XMM-Newton spectra at a minimum of one background-subtracted count per bin. The X-ray modelling was performed with XSPEC v.12.9.0i by adopting the Cash statistic (C-stat) with the implemented direct background subtraction (Cash 1979; Wachter et al. 1979).

Appendix C: X-ray spectral analysis

The X-ray modelling, which was consistently carried out with the same model for all of the sources, was performed in the 0.2–10 keV and 0.3–8 keV for the XMM-Newton and Chandra detectors, respectively. The adopted model is a simple power-law model absorbed by the Galactic interstellar medium and by the obscuring medium (i.e. nuclear absorber+interstellar matter) at the redshift of the source. The obscuration is parametrised as uniform cold absorbers adopting the WABS and ZWABS multiplicative models for Galactic and intrinsic absorbers, respectively. In the cases of sources with spectra from more than one detector (i.e. either for XMM-Newton or XMM-Newton+Chandra), a joint fit is performed with the addition in the modelling of a constant term to account for cross-calibration and possible source flux variation (in the case of observations taken at different epochs). During the modelling, we tied the pn and MOS constants whenever the MOS one was exceeded by more than 7% the pn constant8. For all of the sources except J0900+4215, the calibration constants have been linked among the XMM-Newton detectors. For J0900+4215 (the source with the best quality spectra), the MOS spectra were found to have constants ∼3% higher than the PN. For three sources, a joint fit analysis with XMM-Newton and Chandra was performed (i.e. J0900+4215, J1106+6400, J2123−0050). We find J0900+4215 and J1106+6400 to exhibit Chandra calibration constants, which significantly differ from the XMM-Newton ones by a factor ∼0.6 ± 0.1 and 1.8 ± 0.3 (errors are 90% level), respectively. This indicates that the sources have varied their flux between the two observations. For J2123−0050, the Chandra constant is a factor > 2 higher than the XMM-Newton ones. This is likely due to flux loss in the XMM-Newton spectra due to the presence of a source at ∼20 arcsec from the quasar, which forced us to shrink the quasar spectral extraction region to a circle with a radius of 12 arcsec in order to minimise the contamination. Apart from the calibration constants, the modelling was performed by linking all of the other parameters and leaving three parameters free to vary, that is, the source column density (NH), photon index (Γ), and power-law normalisation. Error estimation for NH and Γ was performed by leaving all of the previously discussed parameters free to vary during the calculation. The uncertainty on the unabsorbed X-ray luminosities (i.e. 2 keV and 2–10 keV) were estimated by freezing Γ to its best-fit value. For each source, Table C.1 reports the fit statistic (C-stat) and degrees of freedom (d.o.f.) of the modelling along with Γ, NH, the observed fluxes at 0.5–2 keV and 2–10 keV (f0.5 − 2 and f2 − 10, respectively), L2 − 10, and the corrected (for extinction and absorption) αOX and ΔαOX (based on Just et al. 2007). Determinations of αOX from the literature are usually derived based on 2 keV luminosity estimated under the assumption of a power-law X-ray spectral model without accounting for intrinsic absorption. Hence, for consistency, we also derived and by estimating the 2 keV luminosity from simple unabsorbed power-law spectral modelling (we only accounted for Galactic absorption). As for the quantities regarding fluxes and luminosities (i.e. f0.5 − 2, f2 − 10, L2 − 10, αOX, and ΔαOX), we adopted and report the pn-derived values for all of the sources observed by XMM-Newton-only in Table C.1. For the sources for which we performed joint XMM-Newton and Chandra modellings, we adopted the XMM-Newton-derived quantities for J1106+6400 (the Chandra observations are of lower quality) and the Chandra ones for J0900+4215 and J2123−0050. Indeed, for J0900+4215, the Chandra observation was performed closest in time to the SDSS spectrum in which the relative vCIV was measured. For J2123−0050, the Chandra spectrum does not suffer from contamination from the nearby X-ray source.

Table C.1.

Properties of the WISSH and C17 samples derived from the X-ray spectral analysis.

We tested for possible X-ray model-dependent systematics on the derived L2 − 10. For each WISSH source with best-fit values or upper limits on NH larger than 1022 cm−2, we parametrised the intrinsic absorption with a partially ionised absorber by adopting the model ZXIPCF with the covering fraction fc = 1. We left both NH and the ionisation parameter (ξ) free to vary. We obtain an estimated L2 − 10, which is on average 0.15 dex larger than the luminosities derived by using a cold absorption model. In this case, we still have a significant L2 − 10-vCIV relation (ρ = 0.74 and p = 0.006).

We mention that restricting the energy range of the WISSH spectral modellings by limiting the low energy bound to 0.5 keV for all the instruments does not make a significant difference in our result. Indeed, L2 − 10 does not vary by more than ±0.05 dex. This level of variation does not change the strength of the reported correlation with vCIV (i.e. ρ = 0.75 and p = 0.004).

All Tables

Table 1.

Correlation and linear regression coefficients.

Table A.1.

Properties of the WISSH quasars and the additional sources from the C17 sample.

Table B.1.

Log of the X-ray observations of the sources from the WISSH and C17 samples.

Table C.1.

Properties of the WISSH and C17 samples derived from the X-ray spectral analysis.

All Figures

thumbnail Fig. 1.

Left panel: L2 − 10 as a function of vCIV for the WISSH quasars (black circles). We also report hyper-luminous quasars at z = 2 − 3 (C17 sample, blue triangles). Middle panel: intrinsic αOX as a function of vCIV for the WISSH quasars. Red circles indicate X-ray weak sources (see Sect. 3 for details). Right panel: L2500 Å as a function of vCIV. Solid and dashed lines report the best-fit linear relation and 1σ uncertainties on the normalisation for the WISSH-only relations. Dotted line in the left panel reports the best-fit relation with the addition of the C17 sample. Orange lines in the right panel report the best-fit L2 − 10vCIV relation normalised to the mean values of L2500 Å.

In the text
thumbnail Fig. 2.

Top panels: log L6μm (left) and L6μm/L2 − 10 (right) as a function of vCIV. Bottom panels: log Lbol (left) and kbol, X (right) as a function of vCIV. The meaning of the data points and lines is the same as in Fig. 1.

In the text
thumbnail Fig. 3.

Γ and NH as a function of vCIV (top and bottom panels). The meaning of the data points and linear fits is the same as in Fig. 1.

In the text
thumbnail Fig. A.1.

L6μm versus z for our WISSH sub-sample (filled red circles) and the entire WISSH sample (empty black circles).

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.