We propose a modified version of the observed non-linear relation between the X-ray (2 keV) and the ultraviolet (2500 Å) emission in quasars (i.e. ) which involves the full width at half-maximum, , of the broad emission line, i.e. . By analysing a sample of 550 optically selected non-jetted quasars in the redshift range of 0.36–2.23 from the Sloan Digital Sky Survey cross matched with the XMM-Newton catalogue 3XMM-DR6, we found that the additional dependence of the observed correlation on the of the MgII broad emission line is statistically significant. Our statistical analysis leads to a much tighter relation with respect to the one neglecting , and it does not evolve with redshift. We interpret this new relation within an accretion disc corona scenario where reconnection and magnetic loops above the accretion disc can account for the production of the primary X-ray radiation. For a broad line region size depending on the disc luminosity as , we find that , which is in very good agreement with the observed correlation. \helveticabold
1 Keywords:active galactic nuclei, quasar, supermassive black holes, accretion disc, X-ray
The physical relation between disc and coronal emission]The physical relation between disc and coronal emission in quasars
Lusso & Risaliti]Lusso Elisabeta and Risaliti Guido \correspondance
One of the observational evidences for the link between the accretion disc and the X-ray corona in active galactic nuclei (AGN) is given by the observed non-linear correlation between the monochromatic ultraviolet luminosity at 2500 Å () and the one in the X–rays at 2 keV (). Such relationship (parameterised as ) exhibits a slope, , around 0.6, implying that optically bright AGN emit relatively less X-rays than optically faint AGN . The value of the slope in this relation does not depend on the sample selection, the correlation is also very tight ( dex, Lusso and Risaliti 2), and independent on redshift [3, 4, 5, 6, 7, 8, 9, 10, 11]. Recently, the relationship (or, more precisely, its version with fluxes) has also been employed as a distance indicator to estimate cosmological parameters such as and by building a Hubble diagram in a similar way as for Type Ia supernovae [12, 13, and references therein], but extending it up to [14, 15]. Yet, the challenge in interpreting such relation on physical grounds is that the origin of the X–ray emission in quasars is still a matter of debate.
Some attempts at explaining the quasar X-ray spectra were based upon the reprocessing of radiation from a non-thermal electron-positron pair cascade (e.g. Svensson 16, 17, Zdziarski et al. 18). Another possibility rests on a two-phase accretion disc model, where a fraction of gravitational power is dissipated via buoyancy and reconnection of magnetic fields in a uniform, hot ( keV K) plasma close to the cold opaque disc [19, 20, 21, 22, 23]. The scenario for this model is the following. Phase 1 is the optically thick “cold” (tens of eV) disc, whilst phase 2 is a hot optically thin plasma located above (and below) the disc. The seed disc photons illuminate the hot tenuous plasma, and a fraction of them is up-scattered to hard X-rays via inverse Compton scattering, providing the main source of cooling of the plasma. About half of these hard X-ray photons is irradiated back to the disc, contributing to its energy balance, whilst the rest escape and are observed. A stable disc-corona system is then in place only if there is a strong coupling between ultraviolet and X-ray photons. This model retrieves the photon index slope of the coronal X-ray spectrum (i.e. ), in close agreement with real data. Nonetheless, it also predicts a nearly equal amount of ultraviolet and X-ray radiation (i.e. ), which is not in agreement with the observed correlation. If the corona is not uniform, but a rather patchy medium, and only a fraction () of the accretion power is released in the hot phase, the resulting value is , but one needs to consider a rather arbitrary number of active blobs . Magnetic field turbulence has been recognised not only as an additional heating mechanism in the formation of the corona (e.g. Galeev et al. 25, Merloni and Fabian 26, Liu et al. 27), but also as an efficient means for the transport of the disc angular momentum (e.g. Balbus 28). Yet, the value of the fraction of gravitational power dissipated in the hot corona (needed for keeping the plasma at high temperatures), and how the coronal physical state depends on the black hole mass and disc accretion rate still remain a matter of debate.
In Lusso and Risaliti  we outlined a simple but physically motivated, ad-hoc model to interpret the observed correlation between the ultraviolet and the X-ray emission in terms of physical parameters such as the black hole mass (, here we considered the normalized value ), the accretion rate (, where is the accretion rate at Eddington), and the distance to the black hole (, where is the Schwarzschild radius). Our main aim was to link such relations with observable quantities (i.e. the observed ultraviolet and X-ray luminosities), thus obtaining a relation that can be then compared with the data.
Here, we further analyse the relation with the goal of understanding its physical origin in the context of black hole accretion physics.
3 The observed plane
The quasar sample considered by Lusso and Risaliti  is obtained by cross-matching the quasar SDSS catalogue published by Shen et al.  with the serendipitous X-ray source catalogue 3XMM–DR6 . Filters are also applied in order to select a clean quasar sample where biases and contaminants are minimised, namely: (i) jetted, broad absorption line quasars and sources with high levels of absorption in the optical () are removed from the sample; (ii) we selected those quasars with a full-width half-maximum, , for the MgII 2800 Å higher than 2000 km s; (iii) only sources with good X-ray data (i.e. S/N5 in the 0.2–12 keV EPIC band) and low levels of X-ray absorption (i.e. with an X-ray photon index ) are considered; and finally (iii) the Eddington bias (see also Risaliti and Lusso 14 and Lusso and Risaliti 2 for further details) is minimised by including only quasars whose minimum detectable X-ray flux is lower than the expected one in each observation. The FWHM values for the MgII emission line as well as the ultraviolet 2500Å luminosities are taken from the Shen et al.  catalogue, which have been measured from the continuum/emission line fitting of the SDSS spectra (see their Section 3 for further details). The final clean sample is composed of 550 quasars and it is only 25% of the initial sample (2,100 quasars with both soft 0.5–2 keV and hard 2–12 keV fluxes), resulting from the stringent filters mentioned above. In the future we will refine the treatment of the systematics by considering additional instrumental systematics, and developing large mock simulations to better understand the impact of the Eddington bias on the results. The main aim is to reduce the rejection fraction, thus maximizing the statistics of the final sample.
Figure 1 shows an edge-on view of the plane. The best-fit regression relation for this quasar sample is
with an observed dispersion of 0.22 dex, which is a very tight relation. To fit the data, we adopted emcee , which is a pure-Python implementation of Goodman & Weare’s affine invariant Markov chain Monte Carlo ensemble sampler.
We then divided the quasar sample into equally spaced narrow redshift intervals in with a to minimise the scatter due to the different luminosity distances within each interval. We split the sample into 10 intervals with and, for each redshift interval, we performed a fit of the relation, , with free and . The results for the best-fit slopes and as a function redshift are shown in Figure 1. Despite the large scatter, both and slopes do not show any significant evolution with time and they are consistent with a constant value in the redshift interval 0.3–2.
4 The toy model
Previous works in the literature suggest that the amount of gravitational energy from the accretion disc released in the hot plasma surrounding the disc itself, likely depends on the black hole mass , the mass accretion rate , and the spin of the black hole. In Lusso and Risaliti , the correlation was reproduced if the energy transfer from an optically thick, geometrically thin accretion disc to the corona is confined to the gas-pressure dominated region of the disc. In such simple model, the monochromatic ultraviolet and X–ray luminosities show an extra dependence on the , thus on the black hole mass and accretion rate as and , respectively. Assuming a broad line region size function of the disc luminosity as we have that
where is the standard disc viscosity parameter, is a calibration constant in the broad line region radius-luminosity relation (i.e. Trippe 33), and .
In this model the bulk of the radiation budget of the X-ray corona is measured at the transition radius () where the gas pressure in the accretion disc equates the radiation pressure, and it is defined as
where , and . The transition radius can vary from a few gravitational radii () to several hundreds for a black hole mass , with an accretion rate in the range , and . This parameter also depends on the value of (i.e. the higher is , the smaller is ). Yet, the effective location of the corona may also be placed at small radii of less than tens of . From a qualitative perspective, the torque (caused by the disc rotation) of magnetic field lines (originated at ) that connects two opposite sides of the accretion disc (e.g. magnetic reconnection, a schematic representation of our model is provided in Figure 2) could cause particle acceleration at the magnetic reconnection site, located closer to the black hole (at ), where particles lose their energy radiatively via the inverse Compton process.
5 Constraining the fraction of accretion power released in the corona
Our simple (but physically motivated) model also provides constraints on the fraction of accretion power released in the hot plasma in the vicinity of the accretion disc, , as a function the broad line region size, from the observed normalization of plane.
Figure 3 shows how the normalization of equation (2) changes as a function of , , and . We fixed the factor to , which corresponds to a broad line region size of about cm ( light days) and a bolometric luminosity of erg s, typical of quasars (e.g. Kaspi et al. 34). The viscosity parameter varies from 0.1 to 0.4 , whilst ranges in the interval 0–0.99. Although the normalization of the relation dsiplays a large scatter, which can be a factor of in logarithm (i.e. orders of magnitude), the data suggest a value of in the range 0.9–0.95.
In the future, by studying larger sample of quasars with and accretion rate values (which are all related to ), we will compare our results on with the expectations from numerical simulations (e.g. Jiang et al. 36) and other theoretical models (e.g. Merloni and Fabian 26, Cao 37, Liu et al. 38), especially those concerning the physical state and extent of the broad line region (e.g. Failed Radiatively Accelerated Dusty Outflows, FRADO, Czerny et al. 39).
Our modified relationship in quasars, which takes into account the full-width half-maximum of the quasar emission line (i.e. ), has an observed dispersion of 0.2 dex over 3 orders of magnitude in luminosity and indicates that there is a good “coupling” between the disc, emitting the primary radiation, and the hot-electron corona, emitting X–rays.
We interpreted such relation through a simple (but physically motivated) model based on the ones presented by Svensson and Zdziarski  and Merloni and Fabian , where a geometrically thin, optically thick accretion disc is coupled with a uniform hot plasma. We assumed that the bulk of the corona emission is mainly powered by the accretion disc and it is located at the transition radius () where the gas pressure equates the radiation pressure in the disc. Assuming a broad line region size function of the bolometric luminosity as we have that , which leads to the final relation . Such a relation is remarkably consistent with the fit obtained from a sample of 550 optically selected quasars from SDSS DR7 cross matched with the XMM–Newton catalogue 3XMM-DR6. The toy model we presented, although simplistic, is capable of making robust predictions on the X–ray luminosities (at a given ultraviolet emission and ) of unobscured/blue quasars, and puts observational constraints on the fraction of accretion power released in the hot plasma in the vicinity of the accretion disc, , as a function the broad line region size, from the observed normalization of plane. The latter result will provide a new vantage point for the next generation of semi-analytical and magneto-hydrodinamical simulations investigating on the physical link between the accretion disc and the X-ray corona.
The proposed relation does not show significant evolution with time in the redshift range covered by our data, , and thus it can be employed as a cosmological indicator to estimate cosmological parameters (e.g. , ).
EL is supported by a European Union COFUND/Durham Junior Research Fellowship (under EU grant agreement no. 609412). This work has been supported by the grants PRIN-INAF 2012 and ASI INAF NuSTAR I/037/12/0.
For all catalogue correlations we have used the Virtual Observatory software TOPCAT  available online (http://www.star.bris.ac.uk/mbt/topcat/). This research has made use of data obtained from the 3XMM XMM–Newton serendipitous source catalogue compiled by the ten institutes of the XMM–Newton Survey Science Centre selected by ESA. This research made use of matplotlib, a Python library for publication quality graphics .
- Tananbaum et al.  Tananbaum H, Avni Y, Branduardi G, Elvis M, Fabbiano G, Feigelson E, et al. X-ray studies of quasars with the Einstein Observatory. ApJ 234 (1979) L9–L13. doi:10.1086/183100.
- Lusso and Risaliti  Lusso E, Risaliti G. The Tight Relation between X-Ray and Ultraviolet Luminosity of Quasars. ApJ 819 (2016) 154. doi:10.3847/0004-637X/819/2/154.
- Vignali et al.  Vignali C, Brandt WN, Schneider DP. X-Ray Emission from Radio-Quiet Quasars in the Sloan Digital Sky Survey Early Data Release: The _ox Dependence upon Ultraviolet Luminosity. AJ 125 (2003) 433–443. doi:10.1086/345973.
- Strateva et al.  Strateva IV, Brandt WN, Schneider DP, Vanden Berk DG, Vignali C. Soft X-Ray and Ultraviolet Emission Relations in Optically Selected AGN Samples. AJ 130 (2005) 387–405. doi:10.1086/431247.
- Steffen et al.  Steffen AT, Strateva I, Brandt WN, Alexander DM, Koekemoer AM, Lehmer BD, et al. The X-Ray-to-Optical Properties of Optically Selected Active Galaxies over Wide Luminosity and Redshift Ranges. AJ 131 (2006) 2826–2842. doi:10.1086/503627.
- Just et al.  Just DW, Brandt WN, Shemmer O, Steffen AT, Schneider DP, Chartas G, et al. The X-Ray Properties of the Most Luminous Quasars from the Sloan Digital Sky Survey. ApJ 665 (2007) 1004–1022. doi:10.1086/519990.
- Green et al.  Green PJ, Aldcroft TL, Richards GT, Barkhouse WA, Constantin A, Haggard D, et al. A Full Year’s Chandra Exposure on Sloan Digital Sky Survey Quasars from the Chandra Multiwavelength Project. ApJ 690 (2009) 644–669. doi:10.1088/0004-637X/690/1/644.
- Lusso et al.  Lusso E, et al. The X-ray to optical-UV luminosity ratio of X-ray selected type 1 AGN in XMM-COSMOS. A&A 512 (2010) A34. doi:10.1051/0004-6361/200913298.
- Young et al.  Young M, Elvis M, Risaliti G. The X-ray Energy Dependence of the Relation Between Optical and X-ray Emission in Quasars. ApJ 708 (2010) 1388–1397. doi:10.1088/0004-637X/708/2/1388.
- Marchese et al.  Marchese E, Della Ceca R, Caccianiga A, Severgnini P, Corral A, Fanali R. The optical-UV spectral energy distribution of the unabsorbed AGN population in the XMM-Newton Bright Serendipitous Survey. A&A 539 (2012) A48. doi:10.1051/0004-6361/201117562.
- Jin et al.  Jin C, Ward M, Done C. A combined optical and X-ray study of unobscured type 1 active galactic nuclei - II. Relation between X-ray emission and optical spectra. MNRAS 422 (2012) 3268–3284. doi:10.1111/j.1365-2966.2012.20847.x.
- Suzuki et al.  Suzuki N, Rubin D, Lidman C, Aldering G, Amanullah R, Barbary K, et al. The Hubble Space Telescope Cluster Supernova Survey. V. Improving the Dark-energy Constraints above z ¿ 1 and Building an Early-type-hosted Supernova Sample. ApJ 746 (2012) 85. doi:10.1088/0004-637X/746/1/85.
- Betoule et al.  Betoule M, Kessler R, Guy J, Mosher J, Hardin D, Biswas R, et al. Improved cosmological constraints from a joint analysis of the SDSS-II and SNLS supernova samples. A&A 568 (2014) A22. doi:10.1051/0004-6361/201423413.
- Risaliti and Lusso  Risaliti G, Lusso E. A Hubble Diagram for Quasars. ApJ 815 (2015) 33. doi:10.1088/0004-637X/815/1/33.
- Risaliti and Lusso  Risaliti G, Lusso E. Cosmology with AGN: can we use quasars as standard candles? Astronomische Nachrichten 338 (2017) 329–333. doi:10.1002/asna.201713351.
- Svensson  Svensson R. The pair annihilation process in relativistic plasmas. ApJ 258 (1982) 321–334. doi:10.1086/160081.
- Svensson  Svensson R. Steady mildly relativistic thermal plasmas - Processes and properties. MNRAS 209 (1984) 175–208. doi:10.1093/mnras/209.2.175.
- Zdziarski et al.  Zdziarski AA, Ghisellini G, George IM, Fabian AC, Svensson R, Done C. Electron-positron pairs, Compton reflection, and the X-ray spectra of active galactic nuclei. ApJ 363 (1990) L1–L4. doi:10.1086/185851.
- Haardt and Maraschi  Haardt F, Maraschi L. A two-phase model for the X-ray emission from Seyfert galaxies. ApJ 380 (1991) L51–L54. doi:10.1086/186171.
- Haardt and Maraschi  Haardt F, Maraschi L. X-ray spectra from two-phase accretion disks. ApJ 413 (1993) 507–517. doi:10.1086/173020.
- Svensson and Zdziarski  Svensson R, Zdziarski AA. Black hole accretion disks with coronae. ApJ 436 (1994) 599–606. doi:10.1086/174934.
- Di Matteo  Di Matteo T. Magnetic reconnection: flares and coronal heating in active galactic nuclei. MNRAS 299 (1998) L15–l20. doi:10.1046/j.1365-8711.1998.01950.x.
- Różańska and Czerny  Różańska A, Czerny B. Vertical structure of the accreting two-temperature corona and the transition to an ADAF. A&A 360 (2000) 1170–1186.
- Haardt et al.  Haardt F, Maraschi L, Ghisellini G. A model for the X-ray and ultraviolet emission from Seyfert galaxies and galactic black holes. ApJ 432 (1994) L95–L99. doi:10.1086/187520.
- Galeev et al.  Galeev AA, Rosner R, Vaiana GS. Structured coronae of accretion disks. ApJ 229 (1979) 318–326. doi:10.1086/156957.
- Merloni and Fabian  Merloni A, Fabian AC. Coronal outflow dominated accretion discs: a new possibility for low-luminosity black holes? MNRAS 332 (2002) 165–175. doi:10.1046/j.1365-8711.2002.05288.x.
- Liu et al.  Liu BF, Mineshige S, Shibata K. A Simple Model for a Magnetic Reconnection-heated Corona. ApJ 572 (2002) L173–L176. doi:10.1086/341877.
- Balbus  Balbus SA. Enhanced Angular Momentum Transport in Accretion Disks. ARAA 41 (2003) 555–597. doi:10.1146/annurev.astro.41.081401.155207.
- Lusso and Risaliti  Lusso E, Risaliti G. Quasars as standard candles. I. The physical relation between disc and coronal emission. A&A 602 (2017) A79. doi:10.1051/0004-6361/201630079.
- Shen et al.  Shen Y, Richards GT, Strauss MA, Hall PB, Schneider DP, Snedden S, et al. A Catalog of Quasar Properties from Sloan Digital Sky Survey Data Release 7. ApJS 194 (2011) 45. doi:10.1088/0067-0049/194/2/45.
- Rosen et al.  Rosen SR, Webb NA, Watson MG, Ballet J, Barret D, Braito V, et al. The XMM-Newton serendipitous survey. VII. The third XMM-Newton serendipitous source catalogue. A&A 590 (2016) A1. doi:10.1051/0004-6361/201526416.
- Foreman-Mackey et al.  Foreman-Mackey D, Hogg DW, Lang D, Goodman J. emcee: The MCMC Hammer. PASP 125 (2013) 306–312. doi:10.1086/670067.
- Trippe  Trippe S. AGN Broad Line Regions Scale with Bolometric Luminosity. Journal of Korean Astronomical Society 48 (2015) 203–206. doi:10.5303/JKAS.2015.48.3.203.
- Kaspi et al.  Kaspi S, Maoz D, Netzer H, Peterson BM, Vestergaard M, Jannuzi BT. The Relationship between Luminosity and Broad-Line Region Size in Active Galactic Nuclei. ApJ 629 (2005) 61–71. doi:10.1086/431275.
- King et al.  King AR, Pringle JE, Livio M. Accretion disc viscosity: how big is alpha? MNRAS 376 (2007) 1740–1746. doi:10.1111/j.1365-2966.2007.11556.x.
- Jiang et al.  Jiang YF, Stone JM, Davis SW. A Global Three-dimensional Radiation Magneto-hydrodynamic Simulation of Super-Eddington Accretion Disks. ApJ 796 (2014) 106. doi:10.1088/0004-637X/796/2/106.
- Cao  Cao X. An accretion disc-corona model for X-ray spectra of active galactic nuclei. MNRAS 394 (2009) 207–213. doi:10.1111/j.1365-2966.2008.14347.x.
- Liu et al.  Liu JY, Qiao EL, Liu BF. Revisiting the Structure and Spectrum of the Magnetic-reconnection-heated Corona in Luminous AGNs. ApJ 833 (2016) 35. doi:10.3847/1538-4357/833/1/35.
- Czerny et al.  Czerny B, Li YR, Hryniewicz K, Panda S, Wildy C, Sniegowska M, et al. Failed Radiatively Accelerated Dusty Outflow Model of the Broad Line Region in Active Galactic Nuclei. I. Analytical Solution. ApJ 846 (2017) 154. doi:10.3847/1538-4357/aa8810.
- Taylor  Taylor MB. TOPCAT & STIL: Starlink Table/VOTable Processing Software. Shopbell P, Britton M, Ebert R, editors, Astronomical Data Analysis Software and Systems XIV (2005), Astronomical Society of the Pacific Conference Series, vol. 347, 29.
- Hunter  Hunter JD. Matplotlib: A 2D Graphics Environment. Computing in Science and Engineering 9 (2007) 90–95. doi:10.1109/MCSE.2007.55.