Original Research ARTICLE
A Hubble Diagram for Quasars
- 1Harvard-Smithsonian Center for Astrophysics, Cambridge, MA, United States
- 2INAF, Osservatorio Astrofisico di Arcetri, Firenze, Italy
- 3Dipartimento di Fisica e Astronomia, Università degli Studi di Firenze, Firenze, Italy
- 4Centre for Extragalactic Astronomy, Department of Physics, Durham University, Durham, United Kingdom
The cosmological model is at present not tested between the redshift of the farthest observed supernovae (z ~ 1.4) and that of the Cosmic Microwave Background (z ~ 1,100). Here we introduce a new method to measure the cosmological parameters: we show that quasars can be used as “standard candles” by employing the non-linear relation between their intrinsic UV and X-ray emission as an absolute distance indicator. We built a sample of ~1,900 quasars with available UV and X-ray observations, and produced a Hubble Diagram up to z ~ 5. The analysis of the quasar Hubble Diagram, when used in combination with supernovae, provides robust constraints on the matter and energy content in the cosmos. The application of this method to forthcoming, larger quasar samples, will also provide tight constraints on the dark energy equation of state and its possible evolution with time.
Quasars are among the brightest sources in the universe, by now observable up to redshift z ~ 7 (Mortlock et al., 2011). For this reason they have always been considered as potential candidates for extending the distance ladder in a redshift range well beyond the limit imposed by the supernovae (z ~ 1.4). Quasars, however, are known to be extremely variable, anisotropic sources and characterized by a wide range in luminosity. Unlike supernovae, an “easy” connection between a spectral (or time-dependent) feature and the luminosity is not available. In short, the use of quasars as standard candles is not obvious. The fundamental requirement needed to employ these sources for a cosmological purpose is being able to measure a “standard luminosity” from which infer the distance.
Several attempts have been made in this sense, using different relations involving quasars luminosity, such as the Baldwin effect (Baldwin, 1977; Korista et al., 1998), the Broad Line Region radius—luminosity relation (Watson et al., 2011; Kilerci Eser et al., 2015), the wavelength-dependent time delays in the emission variability of the accretion disc (Collier et al., 1999), the distribution of the linewidths as a function of the inclination of the BLR and of the quasar luminosity function (Rudge and Raine, 1999), the properties of highly accreting quasars (Wang et al., 2013; Marziani and Sulentic, 2014), the relation between the mass of the SuperMassive Black Hole and the X-ray variability (La Franca et al., 2014) among the others, or invoking the use of classic geometrical arguments, such as the measurement of the apparent size of the Broad Line Region (Elvis and Karovska, 2002), whose intrinsic dimension can be inferred from reverberation mapping monitoring (Blandford and McKee, 1982; Peterson, 1993). Most of these methods suffer from the high scatter in the observed relation or are limited by a poor statistics—generally imposed by the long times required by the observations. In order to make cosmology with quasars we need to (1) have a measurement of their distance with a high precision, i.e., we need a luminosity-related relation with the smallest possible scatter; (2) make the most out of their large statistics, i.e., being able to apply the method to a large number of sources.
We recently proposed the relation between the luminosities in the X-rays (2 keV) and UV (2,500 Å) bands (Tananbaum et al., 1979; Zamorani et al., 1981; Steffen et al., 2006; Lusso et al., 2010; Young et al., 2010; Lusso and Risaliti, 2016, 2017) as a method to estimate quasars distances (Risaliti and Lusso, 2015). The method relies on the simple fact that the LX − LUV relation is not linear. If we express it in terms of the fluxes, it becomes a function of the slope, the intercept, the intrinsic dispersion and the luminosity distance, i.e., a function of the cosmological parameters ΩM and ΩΛ. We can then use the fit to the relation or, equivalently, we can build the Hubble diagram for quasars, to infer information on the cosmological parameters. The obvious advantage of this method is that it only requires a measure of the two fluxes to be employed, therefore avoiding long-monitoring programs and allowing its application on a large number of sources up to high redshifts.
In Risaliti and Lusso (2015) we presented the first Hubble diagram for a sample of ~800 quasars for which 2 keV and 2,500 Å measures were available from the literature. We extended the Hubble diagram in a redshift range that was never explored before by any other cosmological probe (up to z ~ 6) and verified the excellent agreement between the distance moduli inferred with supernovae and those inferred with quasars in the common redshift range (z = 0.01–1.4).
The reliability and effectiveness of the method strongly depend, respectively, on the non-evolution of the relation with reshift and the dispersion in the relation, that, as already mentioned, directly affects the precision in distance estimates. This first sample, although properly treated to suit cosmological studies, can be improved in terms of homogeneity, quality and statistics of the data. In order to do that we selected a new sample of ~8,000 objects for which X-rays and UV observations are available, crossmatching the SDSS Data Release 7 (Shen et al., 2011) and Data Release 12 (Pâris et al., 2017) catalogs with the 3XMM-DR5 (Rosen et al., 2016) catalog. The much higher statistics of this sample with respect to the previous one allows us to apply stronger quality cuts in order to significantly decrease the dispersion in the LX − LUV relation and eventually to verify to a much higher extent the non-evolution of the relation with redshift.
2.1. The Dispersion in the LX − LUV Relation
The LX − LUV relation has been known to be characterized by a high dispersion (~ 0.3–0.4 dex) that has deterred from using it as an absolute distance indicator. The observed dispersion, however, is the result of two distinct contributions: an intrinsic scatter in the relation, related to the still unknown physics involved, and an additional scatter due to observational issues. In Lusso and Risaliti (2016) and Lusso and Risaliti (2017) we demonstrated that a large part of the observed dispersion is ascribable to the latter contribution and, through the analysis of a smaller sample of high quality data and of sources with multiple observations available, that the magnitude of the instrinsic dispersion can be attested to be <0.20 dex. Removing the part of the observed dispersion related to observational issues makes the relation much tighter and our distance estimates much more precise. Significant efforts has then been made in this direction, once the major contributions to this “observational” scatter has been recognisee in the following issues:
1. Uncertainties in the measurement of the (2 keV) X-ray flux, which are, on average, a factor of 2 (unlike UV measurements that can be constrained with uncertainties lower than 10%),
2. Absorption in the spectrum in the UV and in the X-ray wavelength ranges (with the effect on the UV band being more severe and also less recognizable with respect to the one on the X-ray spectrum),
3. Variability of the source and non-simultaneity of the observation in the UV and X-ray bands,
4. Inclination effects affecting the intrinsic emission of the accretion disc,
5. Selection effects due to the flux limit in the surveys (the Eddington Bias, i.e., sources with an average flux below the detection threshold can be detected only while on a positive fluctuation).
By filtering the sample according to the quality of the X-ray measurements, the amount of intrinsic absorption in the source and the systematic effects due to the Eddington bias, we obtain a “clean” sample with a scatter in the relation reduced by ~0.3 dex1 (Figure 1). These criteria select a sample of ~1,900 objects out of the initial ~8,000 with a dispersion of ~0.27 dex. Considering that the dispersion in the LX − LUV relation is propagated to the distance modulus and that the dispersion in the Hubble diagram for supernovae 1A is ~0.07 at z ~ 1, this means that ~15 quasars can provide the same cosmological information as one supernova at this redshift.
Figure 1. Comparison between the relation obtained for the whole sample of quasars for which UV and X-ray observations are available (~8,000) and for the sample cleaned by the observational issues that contribute to the observed dispersion (~1,900). After the quality cuts are applied, the observed dispersion is decreased by ~0.3 dex.
2.2. The Non-evolution of the LX − LUV Relation with Redshift
For the relation to be used as a reliable method for distance estimates, we need to check its non-evolution with the redshift. We repeated the analysis performed in Risaliti and Lusso (2015) on the new, larger sample: we divided the sample in redshift bins and evaluated the relation within each one of them. The size of the redshift bin must fullfill two requirements: it has to be small enough so that differences in the luminosity distance within the same bin are negligible with respect to the intrinsic dispersion in the relation (the observed dispersion in each redshift bin is ~0.25 dex on average), but large enough so that the number of sources within each bin makes this analysis meaningful. Moreover, if the differences in the luminosity distances are negligible, we can consider the relation between fluxes instead of luminosities, making this check independent from the assumed cosmological model.
The two requirements are met when the width of the redshift bin is Δlogz = 0.08. The whole redshift range spanned by quasars (with the exclusion of the range z ≅ 0.01–0.3, for which a significant number of objects is not available) yields 12 redshift bins.
For each subsample in the 12 redshift bins, we fitted the X-tay to UV fluxes with the same log-linear relation adopted for the luminosities of the whole sample. The results of these fits for the slope parameter are shown in Figure 2. All the fits are consistent with α = 0.6, with an average < α > = 0.56 ± 0.08.
Figure 2. Fits of the FX − FUV in narrow redshift intervals. The distance differences among the objects in the same bin are small with respect to the dispersion of the relation, so we can use fluxes as proxies of the luminosities. In this way, the analysis is independent from the choice of a cosmological model. In each panel we report the median redshift, the number of objects N, the dispersion σ, and the best fit slope α of the relation.
3. The Hubble Diagram for Quasars
The LX − LUV relation has been used to estimate quasar distances and build a Hubble diagram of quasars. While slope of the relation can be obtained in a cosmology-independent way as described in the previous Section, its absolute calibration requires the comparison with other standard candles (analogously to the calibration of SN1a based on Cepheid stars). We obtained such a calibration by requiring an overlap of the Hubble diagram of SN1a and quasars in the overlapping redshift range, z ~ 0.2–1.4. The result is shown in Figure 3. The main interesting aspects of this work are the following:
1. The Hubble diagram of quasars perfectly overlap with the one of SN1a in the redshift range z = 0.2–1.4: the same cosmological model fits both the supernovae and quasars data, with no significant residuals beyond the expected Gaussian fluctuations (see inset of Figure 3). While the absolute calibration of quasars has been chosen in order to have such an agreement, the shape of the diagram has no free parameter once the slope of the relation has been fixed. Therefore, the match over the whole common redshift range is a strong confirmation of the reliability of our method.
2. There are a few hundred quasars at redshifts higher than z = 1.4, where no supernova has been observed. The Hubble diagram of these objects is a powerful check of the cosmological model, tracing the evolution of the Universe in the poorly investigated z = 1.4–6 redshift interval, the only other cosmological probes testing the expansion of the Universe above z = 1.4 being Gamma Ray Bursts (Ghirlanda et al., 2006) and Ly-α BAO at z ~ 2.4 (du Mas des Bourboux et al., 2017).
3. The average dispersion in the log(DL)-redshift relation is σ ~ 0.27 dex. This is a large improvement with respect to our first work based on literature samples (σ ~ 0.35) and also on our first Hubble diagram based on SDSS quasars (σ ~ 0.30, Risaliti and Lusso, 2017). This is due to our on-going refinement of both UV and X-ray flux measurements and sample selection. Our primary goal in the early phases of this new branch of observational cosmology is to obtain a clean sample, where biases and systematic effects are greatly reduced (a complete discussion of this point will be presented in a forthcoming paper, together with the results of the fits of the Hubble Diagram with several cosmological models). This comes at the expense of sample statistics: the final sample in Figure 3 consists of “only” 1,900 objects, out of a parent sample of about 8,000 quasars. The main rejection criteria are: X-ray absorption, optical-UV reddening, poor quality of X-ray observations. In the future a more detailed spectral analysis of both the X-ray and optical/UV spectra will allow us to recover a large fraction of the rejected quasars, greatly enhancing the power of this method in constraining the cosmological parameters.
Figure 3. Hubble diagram of quasars obtained with our “clean” sample and supernovae 1A. The orange points are single measurements for quasars, while the red points are quasar averages in small redshift bins. Type 1A supernovae are also plotted with cyan points (JLA sample, Betoule et al., 2014). The inset plot shows a zoom of the same diagram in the redshift range where supernova 1A and quasars overlap. In this case both red and cyan points are averages in small redshift bins for quasars and supernovae 1A respectively.
All authors listed have made a substantial, direct and intellectual contribution to the work, and approved it for publication.
Support for this work was provided by the National Aeronautics and Space Administration through Chandra Award Number AR7-18013 X issued by the Chandra X-ray Observatory Center, which is operated by the Smithsonian Astrophysical Observatory for and on behalf of the National Aeronautics Space Administration under contract NAS8-03060. EL is supported by a European Union COFUND/Durham Junior Research Fellowship (under EU grant agreement no. 609412).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
1. ^Since from the initial sample selection, jetted quasars and Broad Absorption Line quasars (BAL) were excluded. These sources are known, respectively, to have an additional contributions to the X-ray emission (with respect to the one coming from the hot corona) and to be heavily obscured.
Betoule, M., Kessler, R., Guy, J., Mosher, J., Hardin, D., Biswas, R., et al. (2014). Improved cosmological constraints from a joint analysis of the SDSS-II and SNLS supernova samples. Astron. Astrophys. 568:A22. doi: 10.1051/0004-6361/201423413
Collier, S., Horne, K., Wanders, I., and Peterson, B. M. (1999). A new direct method for measuring the Hubble constant from reverberating accretion discs in active galaxies. Month. Notices R. Astron. Soc. 302, L24–L28. doi: 10.1046/j.1365-8711.1999.02250.x
du Mas des Bourboux, H., Le Goff, J.-M., Blomqvist, M., Busca, N. G., Guy, J., Rich, J., et al. (2017). Baryon acoustic oscillations from the complete SDSS-III Lyα-quasar cross-correlation function at z = 2.4. ArXiv e-prints.
Kilerci Eser, E., Vestergaard, M., Peterson, B. M., Denney, K. D., and Bentz, M. C. (2015). On the scatter in the radius-luminosity relationship for active galactic nuclei. Astrophys. J. 801:8. doi: 10.1088/0004-637X/801/1/8
La Franca, F., Bianchi, S., Ponti, G., Branchini, E., and Matt, G. (2014). A new cosmological distance measure using active galactic nucleus X-ray variability. Astrophys. J. Lett. 787:L12. doi: 10.1088/2041-8205/787/1/L12
Lusso, E., Comastri, A., Vignali, C., Zamorani, G., Brusa, M., Gilli, R., et al. (2010). The X-ray to optical-UV luminosity ratio of X-ray selected type 1 AGN in XMM-COSMOS. Astron. Astrophys. 512:A34. doi: 10.1051/0004-6361/200913298
Pâris, I., Petitjean, P., Ross, N. P., Myers, A. D., Aubourg, É., Streblyanska, A., et al. (2017). The sloan digital sky survey quasar catalog: twelfth data release. Astron. Astrophys. 597:A79. doi: 10.1051/0004-6361/201527999
Rosen, S. R., Webb, N. A., Watson, M. G., Ballet, J., Barret, D., Braito, V., et al. (2016). The XMM-Newton serendipitous survey. VII. The third XMM-Newton serendipitous source catalogue. Astron. Astrophys. 590:A1. doi: 10.1051/0004-6361/201526416
Rudge, C. M., and Raine, D. J. (1999). Determining the cosmological parameters from the linewidths of active galaxies. Month. Notices R. Astron. Soc. 308, 1150–1158. doi: 10.1046/j.1365-8711.1999.02802.x
Shen, Y., Richards, G. T., Strauss, M. A., Hall, P. B., Schneider, D. P., Snedden, S., et al. (2011). A catalog of quasar properties from sloan digital sky survey data release 7. Astrophys. J. Suppl. 194:45. doi: 10.1088/0067-0049/194/2/45
Steffen, A. T., Strateva, I., Brandt, W. N., Alexander, D. M., Koekemoer, A. M., Lehmer, B. D., et al. (2006). The X-ray-to-optical properties of optically selected active galaxies over wide luminosity and redshift ranges. Astron. J. 131, 2826–2842. doi: 10.1086/503627
Tananbaum, H., Avni, Y., Branduardi, G., Elvis, M., Fabbiano, G., Feigelson, E., et al. (1979). X-ray studies of quasars with the Einstein Observatory. Astrophys. J. Lett. 234, L9–L13. doi: 10.1086/183100
Wang, J. M., Du, P., Valls-Gabaud, D., Hu, C., and Netzer, H. (2013). Super-eddington accreting massive black holes as long-lived cosmological standards. Phys. Rev. Lett. 110:081301. doi: 10.1103/PhysRevLett.110.081301
Young, M., Elvis, M., and Risaliti, G. (2010). The X-ray energy dependence of the relation between optical and X-ray emission in quasars. Astrophys. J. 708, 1388–1397. doi: 10.1088/0004-637X/708/2/1388
Keywords: cosmology: distance scale–cosmological parameters–observations, galaxies: active, quasars: general, X-ray: general, ultraviolet: general
Citation: Bisogni S, Risaliti G and Lusso E (2018) A Hubble Diagram for Quasars. Front. Astron. Space Sci. 4:68. doi: 10.3389/fspas.2017.00068
Received: 12 November 2017; Accepted: 19 December 2017;
Published: 08 January 2018.
Edited by:Mauro D'Onofrio, Università degli Studi di Padova, Italy
Reviewed by:Fabio La Franca, Università degli Studi Roma Tre, Italy
Alenka Negrete, Universidad Nacional Autónoma de México, Mexico
Copyright © 2018 Bisogni, Risaliti and Lusso. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Susanna Bisogni, email@example.com