Quasars: from the Physics of Line Formation to Cosmology

TL;DR

This paper explores extreme Population A quasars with high accretion rates as stable, reliable distance indicators across cosmic history, analyzing their spectral properties and physical conditions for cosmological applications.

Contribution

It identifies spectral features of extreme Population A quasars that enable their use as redshift-independent distance estimators based on virial broadening.

Findings

Spectral properties allow easy identification of xA quasars in large surveys.

Physical conditions support the use of low-ionization line widths as virial estimators.

xA quasars can serve as stable cosmological distance indicators.

Abstract

Quasars accreting matter at very high rates (known as extreme Population A [xA] or super-Eddington accreting massive black holes) provide a new class of distance indicators covering cosmic epochs from the present-day Universe up to less than 1 Gyr from the Big Bang. The very high accretion rate makes it possible that massive black holes hosted in xA quasars radiate at a stable, extreme luminosity-to-mass ratio. This in turns translates into stable physical and dynamical conditions of the mildly ionized gas in the quasar low-ionization line emitting region. In this contribution, we analyze the main optical and UV spectral properties of extreme Population A quasars that make them easily identifiable in large spectroscopic surveys at low-z (z < 1) and intermediate-z (2 < z < 2.6), and the physical conditions that are derived for the formation of their emission lines. Ultimately, the analysis supports the possibility of identifying a virial broadening estimator from low-ionization line widths, and the conceptual validity of the redshift-independent luminosity estimates based on virial broadening for a known luminosity-to-mass ratio.

Figures

Equations

$$U=\frac{{\int }_{{\nu }_0}^{\infty }\frac{L_{\nu }}{h\nu }}{4\pi r_{\mathrm{BLR}}^{2}c{n}_{\mathrm{H}}}=\frac{Q(H)}{4\pi r_{\mathrm{BLR}}^{2}c{n}_{\mathrm{H}}},$$

$$L\propto \text{lamed}{M}_{\mathrm{BH}}\propto \text{lamed}{r}_{\mathrm{BLR}}(\delta v{)}^2$$

$$r_{\mathrm{BLR}}\propto {\left(\frac{L}{n_{\mathrm{H}}U}\right)}^{\frac{1}{2}}.$$

$$L\propto \text{lamed}{\left(\frac{L}{n_{\mathrm{H}}U}\right)}^{\frac{1}{2}}(\delta v{)}^2.$$

$$L\approx 7.8\cdot 10^{44}\frac{{}_1^2{\kappa }_{0.5}{f}_2^2}{{\tilde{\nu }}_{2.42\cdot 10^{16}}}\frac{1}{{\left(n_{\mathrm{H}}U\right)}_{9.6}}(\delta v_{1000}{)}^4$$

Full text

\pubvolume

xx \issuenum1 \articlenumber5

\historyReceived: date; Accepted: date; Published: date \TitleQuasars: from the Physics of Line Formation to Cosmology \AuthorP. Marziani1\orcidA, E. Bon2, N. Bon2, A. del Olmo3, M. L. Martínez-Aldama3, M. D’Onofrio4, D. Dultzin5, C. A. Negrete5, G. M. Stirpe6 \AuthorNamesPaola Marziani, Edi Bon and Natasa Bon

1 Introduction

1.1 Quasar spectra: emission from mildly ionized gas

The spectra of quasars can be easily recognized by the presence of broad and narrow optical and UV lines emitted by mildly- ionized species over a wide range of ionization potential. The type-1 composite quasar spectrum from the SDSS Vanden Berk et al. (2001) reveals Broad (FWHM $\gtrsim$ 1000 km s*-1*) and Narrow High Ionization lines (HILs,   50eV) and Low Ionization lines (LILs, < 20eV). Broad HILs encompass CIV$\lambda$1549, HeII$\lambda$1640 and HeII$\lambda$4686 as representative specimens. Broad LILs include HI Balmer lines (H$\beta$, H$\alpha$), MgII$\lambda$2800, the CaII IR Triplet, and FeII features. The FeII emission deserves a particular mention, as it is extended over a broad range of wavelengths (Fig. 6 of Vanden Berk et al. (2001)), and is especially prominent around MgII$\lambda$2800 and H$\beta$. The FeII emission is one of the dominant coolants in the broad line region (BLR) and therefore a main factor in its energetic balance ( the FeII emission extends from the UV to the FIR Marziani et al. (2006), and can reach the luminosity of Ly$\alpha$, Netzer (1990); Marinello et al. (2016)). So it may not appear surprising that an estimator of its strength plays an important role in the systematic organization of quasar properties (§2).

This paper reviews results obtained in the course of two decades (§2 and §3), attempting to explain how the spectral properties of a class of type-1 quasars and their physical interpretation can lead to the definition of “Eddington standard candles” (ESC, §4). In the following, we will restrict the presentation to type-1 quasars which are considered mainly “unobscured” sources with an unimpeded view of the BLR, and exclude type-2 AGN of quasars in which the broad lines are not detected in natural light (see Antonucci (1993) for an exhaustive review). We describe the physical basis of the method in §3 and § 4.1. We then introduce ESC selection criteria (§5) and preliminary cosmology results (§6) .

1.2 Quasars for cosmology: an open issue

The distribution of quasars in space and the intervening absorptions along the line of sight (i.e., the so called Ly$\alpha$ forest) have been since long considered as tracer of matter in the distant Universe (see D’Onofrio and Burigana (2009) and references therein). However, a relevant question may be why intrinsic properties of quasars have never been successfully used as cosmological probes. On the one hand, (1) quasars are easily recognizable and plentiful ($\gtrsim$500,000 in the data release 14 of the SDSS, Pâris et al. (2018)). (2) They are very luminous, can reach bolometric luminosity $L\gtrsim 10^{48}$ erg s*-1*; (3) they are observed in an extremely broad range of redshift $0\lesssim z\lesssim 7$, and (4) they are stable compared to transients that are employed as distance indicators in cosmology, such as type Ia supernovæ  (Sect. 2, Marziani et al. (2018) for a review). On the other hand, (1) quasars are anisotropic sources even if the degree of anisotropy is expected to be associated with the viewing angle of the accretion disk in radio-quiet quasars Negrete et al. (2018), and not large compared to radio-loud quasars whose optical continuum is in part beamed (see, for example Liu et al. (2006)); (2) quasars have an open-ended luminosity function ( i.e., without a clearly defined minimum, as the quasar highest spatial density occurs at the lowest luminosity); in other words, they are the “opposite” of a cosmological standard candle. In addition, (3) the long-term variability of radio-quiet quasars is poorly understood (see e.g. Bon et al. (2016, 2017) and references therein) (4) and the internal structure of the active nucleus ($\lesssim$ 1000 $r_{\mathrm{g}}$) is still a matter of debate (see e.g. a summary of open issues Netzer (2018) in (Marziani et al., 2018). Correlations with luminosity have been proved to be rather weak (see Sulentic et al. (2000), for a synopsis up to mid-1999). Selection effect may even cancel out the “Baldwin effect” Baldwin et al. (1978), a significant but weak anti-correlation between rest-frame equivalent width and continuum luminosity of CIV$\lambda$1549 that has been the most widely discussed luminosity correlation in the past decades.

2 Definition of a class of type-1 quasars with properties of Eddington standard candles

Nonetheless, new developments in the past decades have paved the road to the possibility of exploiting quasars as cosmological distance indicators in a novel way that would make them literal “Eddington standard candles” (ESC) (Teerikorpi (2011); Wang et al. (2013, 2014); Marziani and Sulentic (2014); see also Czerny et al. (2018) for a comprehensive review of secondary distance indicators including several techniques based on quasars). This possibility is based in the development of the concept of a quasar main sequence (MS), intended to provide a sort of H-R diagram for quasars Sulentic et al. (2008). The quasar MS can be traced in the plane defined by the prominence of optical FeII emission, $R_{\mathrm{FeII}}$ = I(FeII$\lambda$ 4570)/I(H$\beta$) (see Boroson and Green (1992); Sulentic et al. (2000, 2002); Zamfir et al. (2010); Shen and Ho (2014)). Fig. 1 provides a sketch of the MS in the optical plane FWHM(H$\beta$) vs. $R_{\mathrm{FeII}}$. It is possible to isolate spectral types in the optical plane of the MS as a function of $R_{\mathrm{FeII}}$ and FWHM H$\beta$ and, at a coarser level, two populations: Population A (FWHM H$\beta<$4000 km/s) and Population B of broader sources. Pop. A is rather heterogeneous, and encompasses a range of $R_{\mathrm{FeII}}$ from almost 0 to the highest values observed ($R_{\mathrm{FeII}}\gtrsim 2$ are very rare, $\lesssim$ 1% in optically-selected samples, Zamfir et al. (2010)). Along the quasar main sequence, the extreme Population A (xA) sources satisfying the condition $R_{\mathrm{FeII}}>1$  (about $10$ % of all quasars in optically-selected sample, green area in Fig. 1) show remarkably low optical variability, so low that it is even difficult to estimate the BLR radius via reverberation mapping Du et al. (2018). This is at variance with Pop. B sources that show more pronounced variability Dultzin-Hacyan et al. (1992); Giveon et al. (1999), the most extreme cases being observed among blazars which are low-accretors, at the opposite end in the quasar MS. Of the many multi-frequency trends along the main sequence (from the sources whose spectra show the broadest LILs [extreme Pop. B], and the weakest FeII emission, to sources with the narrowest LIL profiles and strongest FeII emission [extreme Pop. A]), we recall a systematic decrease of the CIV equivalent width, an increase in metallicity, and amplitude of HIL blueshifts (a more exhaustive list is provided by Table 1 of Sulentic et al. (2011)). Eddington ratio is believed to increase along with $R_{\mathrm{FeII}}$ (Boroson and Green, 1992; Kuraszkiewicz et al., 2004; Shen and Ho, 2014; Sun and Shen, 2015). The FWHM H$\beta$ is strongly affected by the viewing angle (i.e., the angle between the line of sight and the accretion disk axis), so that at least most narrow-line Seyfert 1s (NLSy1s) can be interpreted as Pop. A sources seen with the accretion disk oriented face-on or almost so Dultzin et al. (2011). At low-$z$ ($\lesssim$0.7), Pop. A implies low black hole mass $M_{\mathrm{BH}}$, and high Eddington ratio; on the converse, Pop. B is associated with high $M_{\mathrm{BH}}$ and low $L/L_{\mathrm{Edd}}$. This trend follows from the “downsizing” of nuclear activity at low-$z$ that helps give an elbow shape to the MS Fraix-Burnet et al. (2017): at low-$z$, very massive quasars ($M_{\mathrm{BH}}\gtrsim 10^{9}$ M*⊙*) do not radiate close to their Eddington limit but are, on the converse, low-radiators ($L/L_{\mathrm{Edd}}\lesssim 0.1$).

The inter-comparison between CIV$\lambda$1549 and H$\beta$ supports low-ionization lines virial broadening (in a system of dense clouds or in the accretion disk) + high-ionization lines (HILs) radial or vertical outflows, at least in Pop. A sources Leighly (2004); Sulentic et al. (2017). There is now a wide consensus on an accretion disk + wind system model Elvis (2000), and therefore on the existence of a “virialized” low-ionization subregion + higher ionization, outflowing subregion up to the highest quasar luminosities Bisogni et al. (2017); Sulentic et al. (2017); Vietri et al. (2018).

The most extreme examples at high accretion rate are a population of sources with distinguishing properties. They have been called extreme Pop. A or extreme quasars (xA), and are also known as super-Eddington accreting massive black holes (SEAMBHs) Negrete et al. (2018); Martínez-Aldama et al. (2018); Wang et al. (2013); Du et al. (2016). Observationally, xA quasars satisfy $R_{\mathrm{FeII}}\geq 1$ and still show LIL H$\beta$ profile consistent with emission from a virialized system. xA quasars may be well represent an early stage in the evolution of quasars and galaxies. In the hierarchical growth scenario for the evolution of galaxies, merging and strong interaction lead to accumulation of gas in the galaxy central regions, inducing enhanced star formation. Strong winds from massive stars and eventual Supernova explosions may ultimately provide enriched accretion fuel for the massive black hole at the galaxy Heller and Shlosman (1994); Collin and Zahn (1999); Williams et al. (1999). The active nucleus radiation force and the mechanical thrust of the accretion disk wind can then sweep the dust surrounding the black hole, at least within a cone coaxial with the accretion disk axis (see Fig. 7 of D’Onofrio and Marziani (2018)). The fraction of mass that is accreted by the black hole and the fraction that is instead ejected in the wind are highly uncertain; the outflow kinetic power can become comparable to the radiative output King and Pounds (2015); Marziani et al. (2017), especially in sources accreting at very high rate Nardini et al. (2015); interestingly, this seems to be true also for stellar-mass black holes Tetarenko et al. (2018). Feedback effects on the host galaxies are maximized by the high kinetic power of the wind, presumably made of gas much enriched in metals Baskin and Laor (2012).

3 Diagnostics of mildly-ionized gases

Diagnostics from the rest-frame UV spectrum takes advantage of the observations of strong resonance lines that are collisionally excited Negrete et al. (2012, 2013). The point is that the rest-frame UV spectrum offers a rich diagnostics that constrains at least gas density $n_{\mathrm{H}}$, ionization parameter $U$, chemical abundance $Z$. For instance, Si II$\lambda$ 1814/Si III]$\lambda$ 1892 sensitive to ionization CIV$\lambda$1549/Ly$\alpha$, CIV$\lambda$1549/(Si IV+OIV])$\lambda$1400, CIV$\lambda$1549/HeII$\lambda$1640, NV$\lambda$1240/HeII$\lambda$1640 are sensitive to metallicity; Al III$\lambda$1860/Si III]$\lambda$1892, Si III]$\lambda$1892/CIII]$\lambda$1909 are sensitive to density, since inter-combination lines have a well defined critical density (Negrete et al., 2012).

The photoionization code Cloudy models the ionization, chemical, and thermal state of gas exposed to a radiation field, and predicts its emission spectra and physical parameters Ferland et al. (2013, 2017). In Cloudy, collisional excitation and radiative processes typical of mildly ionized gases are included. Cloudy simulation require inputs in terms of $n_{\mathrm{H}}$, $U$, $Z$, quasar spectral energy distribution (SED), column density $N_{\mathrm{c}}$. The ionization parameter

[TABLE]

where $Q(H)$ is the number of ionizing photons, provides the ratio between photon and hydrogen number density. More importantly, the inversion of equation provides a measure of the emitting region radius $r_{\mathrm{BLR}}$ once the ionizing photon flux i.e., the product $Un_{H}$ is known. As we will see, the photon flux can be estimated with good precision from diagnostic line intensity ratios.

Maps built on an array of 551 Cloudy 08.00 – 13.00 photoionization models for a given metallicity $Z$ and $N_{\mathrm{c}}$, constant $n$ and $U$ evaluated at steps of 0.25 dex covering the ranges $7\leq\log n_{\mathrm{H}}\leq 14$ [cm*-3*], $-4.5\leq\log U\leq 0$. Given the measured intensity ratios for xA quasars, Cloudy simulations show convergence toward a well-defined value of $\log$ ($n_{\mathrm{H}}U$) Negrete et al. (2012); Martínez-Aldama et al. (2018). UV diagnostic ratios in the plane ionization parameter versus density indicate extremely high $n_{\mathrm{H}}~{}10^{12.5-13}$ cm*-3*, extremely low $\log U\sim-2.5-3$ (Fig. 3). Note the orthogonal information provided by the AlIII$\lambda$1860/SiIII]$\lambda$1892 that mainly depends on density. The left and right panel differ because of chemical abundances: the case with five times solar metallicity plus overabundance of Si and Al produces better agreement, displacing the solution toward lower density and higher ionization. Nonetheless, the product $Un_{\mathrm{H}}$ remains fairly constant. Diagnostic ratios sensible to chemical composition suggest high metallicity. The metallicities in the quasar BLR gas are a function of the ST along the MS: relatively low (solar or slightly sub-solar in extreme Pop. B sources (as estimated recently for NGC 1275 Punsly et al. (2018)), and relatively high for typical Pop. A quasars with moderate FeII emission ($Z\sim 5-10Z_{\odot}$, Shin et al. (2013); Sulentic et al. (2014)). If the diagnostic ratios are interpreted in terms of scaled $Z_{\odot}$, they may reach $Z\gtrsim 20Z_{\odot}$, even $Z\sim 100Z_{\odot}$ for xA quasars Martínez-Aldama et al. (2018). $Z$  values as high as $Z\sim 100Z_{\odot}$ are likely to be unphysical, and suggest relative abundances of elements deviating from solar values, as assumed in the previous example. The analysis of the gas chemical composition in the BLR of xA source has just begun. However, high or non-solar $Z$ are in line with the idea of xA sources being high accretors surrounded by huge amount of gas and a circum-nuclear star forming system, possibly with a top-heavy initial mass function Negrete et al. (2012). The high $n_{\mathrm{H}}$ is consistent with the low CIII]$\lambda$1909 emission that becomes undetectable in some cases. While in Pop. A and B we find evidence of ionization stratification within the low-ionization part of the BLR (Peterson and Wandel (2000, 1999), Gaskell (2009) and references therein), xA sources show intensity ratios that are consistent with a very dense “remnant” of the BLR, perhaps after lower density gas has been ablated away by radiation forces.

4 xA quasars as Eddington standard candles

There are several key elements that make it possible to exploit xA quasars as Eddington standard candles.

The first is the similarity of their spectra and hence of the physical condition in the mildly-ionized gas that is emitting the LILs. Line intensity ratios are similar (they scatter around a constant average with small dispersion). Since the line emitting gas is photoionized, intensity line ratios depend strongly on the ionizing continuum SED. Thus, also the ionizing SED is constrained within a small scatter. We remark that this is not true for the general population of quasars that show differences in line equivalent width and intensity ratios larger than an order of magnitude along the MS.

The mass reservoir in all xA sources is sufficient to ensure a very high accretion rate (possibly super-Eddington) that yields a radiative output close to the Eddington limit. The similarity of the SED and the presence of high rates of circumnuclear and galactic star formation as revealed by Spitzer Sani et al. (2010), has led to the conjecture that xA sources may be in particular stage of a quasar development, as mentioned above.

The second key element is the existence of a virialized low-ionization sub-region (possibly the accretion disk itself). This region coexists with outflowing gas even at extreme $L\gtrsim 10^{48}$ erg s*-1* and highest Eddington ratios, but is kinematically distinguishable on the basis of inter-line shifts between LILs and HILs, for exampleH$\beta$ and CIV$\lambda$1549.

In addition, xA quasars show extreme $L/L_{\mathrm{Edd}}$ along the MS with small dispersion. If the Eddington ratio is known, and constant, then $\lamed=L/L_{\mathrm{Edd}}\propto L/M_{\mathrm{BH}}$. Accretion disk theory teaches low radiative efficiency at high accretion rate, and that saturates toward a limiting values (Abramowicz et al. (1988); Mineshige et al. (2000); Abramowicz and Straub (2014) and references therein). Therefore, empirical evidence (the xA class of sources easily identified by their self-similar properties, scatters around a well-defined, extremal ) and theoretical support (the saturation of the radiative output per unit $M_{\mathrm{BH}}$) justified the consideration of xA sources potential ESCs.

4.1 Virial luminosity

The use of xA sources as Eddington standard candles requires several steps which should considered carefully.

The first step is the actual estimate of the accretion luminosity via a virial broadening estimator (VBE). The luminosity can be written as

[TABLE]

assuming virial motions of the low-ionionization part of the broad-line region (BLR). The $\delta v$ stands for a suitable VBE, usually the width of a convenient LIL (in practice, the FWHM of H$\beta$ or even Pa$\alpha$, La Franca et al. (2014)). 2. 2.

The $r_{\mathrm{BLR}}$ can be estimated from the inversion of Eq. 1 Negrete et al. (2012, 2013), taking again advantage that the ionizing photon flux shows a small scatter around a well defined value. In addition, another key assumption is that

[TABLE]

Eq. 3 implies that $r_{\mathrm{BLR}}$ scales with the square root of the luminosity. This is needed to preserve the $U$ parameter. If $U$ were going to change then the spectrum would also change as a function of luminosity. This is not evident comparing spectra over a wide luminosity range (4.5 dex), although some second order effects are possible.111The maximum temperature of the accretion disk is $\propto M_{\mathrm{BH}}^{-\frac{1}{2}}$; the SED is expected to become softer at high $M_{\mathrm{BH}}$, but this effect has not been detected yet at a high confidence level. 3. 3.

We can therefore write the virial luminosity as

[TABLE]

Making explicit the dependence of the number of ionizing photons on the SED, the virial luminosity becomes:

[TABLE]

where $\kappa$ is the fraction of ionizing luminosity scaled to 0.5, $\tilde{\nu}$ the average frequency of ionizing photons scaled to $2.42\cdot 10^{16}$ Hz, and $\left({n_{\mathrm{H}}U}\right)$ to $10^{9.6}$.

Eq. 5 is analogous to the Tully-Fisher Tully and Fisher (1977) and the early formulation of the Faber-Jackson Faber and Jackson (1976) laws for galaxies. Eq. 5 is applicable to xA quasars with Eddington ratio $\lamed\sim 1$ and dispersion $\delta\lamed\ll 1$, but in principle could be used to for every sample of quasars whose  is in a very restricted range.

5 Selection of Eddington standard candles

Selection criteria are based on emission line intensity ratios which are extreme along the quasar MS:

$R_{\mathrm{FeII}}$  $>1.0$; 2. 2.

UV AlIII$\lambda$1860/ SiIII]$\lambda$1892$>0.5$; 3. 3.

SiIII]$\lambda$1892/ CIII]$\lambda$1909$>1$

Marziani and Sulentic (2014). The first criterion can be easily applied to optical spectra of a large survey such as the SDSS for sources at $z\lesssim 1$. The second and third criterion can be applied to sources at $1\lesssim z\lesssim 4.5$ for which the 1900 blend lines are shifted into the optical and near IR domains. UV and optical selection criterions are believed to be equivalent. Due to a small sample size at low $z$ for which rest-frame optical and UV spectra are available, further testing is needed.

6 Tentative applications to cosmology and the future perspectives

Preliminary results were collected from 3 quasar samples (62 sources in total), unevenly covering the redshift range $0.4\lesssim z\lesssim 2.6$. For redshift $z\gtrsim 2.$ the UV AlIII$\lambda$1860 FWHM was used as a VBE for the rest-frame UV range, save a few cases for which H$\beta$ was availabele. This explorative application to cosmology yielded results consistent with concordance cosmology, and allowed the exclusion of some extreme cosmologies Marziani and Sulentic (2014). A more recent application involved the Marziani and Sulentic (2014) sample, along with the H$\beta$ sample of Negrete et al. (2018) and preliminary measurements from Martínez-Aldama et al. (2018). The resulting Hubble diagram is shown in Fig. 4. The plots in Fig. 4 involve $\approx 220$ sources and indicate a scatter $\delta\mu\approx$ 1.2 mag. The slope of the residuals ($b\approx-0.002\pm 0.104$) is not significantly different from 0, indicating good statistical agreement between luminosities derived from concordance cosmological parameters and from the virial equation. The Hubble diagram of Fig. 4 confirms the conceptual validity of the virial luminosity relation, Eq. 5.

Mock samples of several hundreds of objects, even with significant dispersion in luminosity with rms($\log L$) = 0.2 – 0.3, indicate that quasars covering the redshift range between 0 and 3 (i.e., a range of cosmic epochs from now to 2 Gyr since the Big Bang) could yield significant constraints on the cosmological parameters. A synthetic sample of 200 sources uniformly distributed in the redshift range 0 – 3 with a scatter of 0.2 dex yields $\Omega_{\mathrm{M}}\approx 0.28\pm 0.02$ at 1$\sigma$ confidence level, assuming $H_{0}=70$ km s*-1* Mpc*-1*, and flatness ($\Omega_{M}+\Omega_{\Lambda}$=1). If $\Omega_{M}+\Omega_{\Lambda}$ is unconstrained, $\Omega_{\mathrm{M}}\approx 0.30^{+0.12}_{-0.09}$ at 1$\sigma$ confidence level Marziani and Sulentic (2014). The comparison between the constraints set by supernova surveys and by a mock sample of 400 quasars with rms = 0.3 dex in $\log L$ shows the potential ability of the quasar sample to better constrain $\Omega_{M}$ Marziani and Sulentic (2014). The scheme of Fig. 5 illustrates the difference in sensitivity to cosmological parameters over the redshift range 0 – 4: supernovæ are sensitive to $\Omega_{\Lambda}$ since the effect of $\Omega_{\Lambda}$, in a concordance cosmology scenario, became appreciable only at relatively recent cosmic epochs. High redshift quasars provide information on a redshift range where the expansion of the Universe was still being decelerated by the effect of $\Omega_{M}$, a range that is not yet covered by any standard ruler or candle.

6.1 Error budget

The large scatter in the luminosity estimates is apparently daunting in the epoch of precision cosmology. Statistical errors can be reduced to rms $\approx 0.2$ dex in $L$ by increasing numbers, collecting large samples ($\sim$ 500 quasars).

In the case of xA quasars, the SED cannot vary much since spectra of xA quasars are almost identical in terms of line ratios (a second order effect Shemmer and Lieber (2015) not yet detected as significant in the data considered by Marziani and Sulentic (2014); Martínez-Aldama et al. (2018) may become significant with larger samples). The scaling $r_{\mathrm{BLR}}\propto L^{0.5}$ should hold strictly: a small deviation would imply a systematic change in the ionization parameter and hence of ST with luminosity.

A simplified error budget for statistical errors Marziani and Sulentic (2014) clearly indicates that the virial luminosity estimate is affected by the VBE uncertainties (i.e., FWHM measurement uncertainties) which enter with the fourth power in Eq. 5. Orientation effects are expected to be determinant in the FWHM uncertainties, as they can contribute 0.3 dex if H$\beta$ or any other line used as a VBE is emitted in a highly-flattened configuration.222 The $M_{\mathrm{BH}}$ is not computed explicitly for the estimate of $L$ following Eq. 5. However, the VBE uncertainty associated with orientation is the main source of uncertainty on $M_{\mathrm{BH}}$ for the xA sample. This is most likely the case also for the general quasar population Jarvis and McLure (2006); Marziani et al. (2017). Modeling the effect of orientation by computing the difference between the $L$ from concordance cosmology and the virial luminosity indeed reduces the sample standard deviation by a factor $\approx$ 5 to $\approx$ 0.2 mag, and accounts for most of the rms $\approx$ 0.4 dex in the virial luminosity estimates of the sample shown in Fig. 4 Negrete et al. (2018). The rms $\approx$ 0.2 mag value is comparable to the uncertainty in supernova magnitude measurements. Work is in progress in order to make viewing angle estimates of xA quasars usable for cosmology.

7 Conclusions

This paper provided an overview of the physical conditions in the broad line emitting region of extreme spectral types of type-1 quasars (the extreme Pop. A). There is strong evidence that xA sources are radiating close to their Eddington limit (i.e., with Eddington ratio scattering around a well-defined value), at high accretion rates. Their physical properties appear to be very stable across a very wide range of luminosity, 4 – 5 dex. The assumption of a constant   makes it possible to write a relation between luminosity and virial broadening, analogous to the one expressed by the Tully-Fisher and the early formulation of Faber-Jackson laws.

The scatter in the Hubble diagram obtained from virial luminosity estimates is still very high, about 1 mag (although comparable to the scatter from a method based on the non-linear relation between the X-ray and the UV emission of quasars Risaliti and Lusso (2015)). Very large samples are needed for reduction of scatter (and statistical error). In addition, the inter-calibration of rest-frame visual and UV properties and their dependence on $L$ (expect systematic errors!) needs to be extended by dedicated observations of xA sources covering the rest frame UV and visual range. Simulations of statistical and systematic effects which influence the estimates of the cosmic parameters are also needed.

In principle, Eddington standard candles can cover a range of distances where the metric of the Universe has not been “charted” as yet to retrieve an independent estimate of $\Omega_{\mathrm{M}}$. If samples with uniform coverage over a wide range of redshift would become available, xA sources could also address the physics of accelerated expansion (i.e., provide measurements of the dark energy equation of state).

\authorcontributions

All authors significantly contributed to the papers on which this review is based.

\funding

PM wishes to thank the SOC of the SPIG 2018 meeting for inviting the topical lecture on which this paper is based, and acknowledges the Programa de Estancias de Investigación (PREI) No. DGAP/DFA/2192/2018 of UNAM, where this paper was written. The relevant research is part of the project 176001 ”Astrophysical spectroscopy of extragalactic objects“ and 176003 “Gravitation and the large scale structure of the Universe" supported by the Ministry of Education, Science and Technological Development of the Republic of Serbia. MLMA acknowledges a CONACyT postdoctoral fellowship. AdO, MLMA acknowledge financial support from the Spanish Ministry for Economy and Competitiveness through grants AYA2013-42227-P and AYA2016-76682-C3-1-P. MLMA, PM and MDO acknowledge funding from the INAF PRIN-SKA 2017 program 1.05.01.88.04. DD and AN acknowledge support from CONACyT through grant CB221398. DD and AN thank also for support from grant IN108716 53 PAPIIT, UNAM.

\conflictsofinterest

The authors declare no conflict of interest.

\abbreviations

The following abbreviations are used in this manuscript:

[TABLE]

\reftitle

References