Quasar Black Hole Mass Estimates from High-Ionization Lines : Breaking a Taboo ?

Can high ionization lines such as CIVλ1549 provide useful virial broadening estimators for computing the mass of the supermassive black holes that power the quasar phenomenon? The question has been dismissed by several workers as a rhetorical one because blue-shifted, non-virial emission associated with gas outflows is often prominent in CIVλ1549 line profiles. In this contribution, we first summarize the evidence suggesting that the FWHM of low-ionization lines like Hβ and MgIIλ2800 provide reliable virial broadening estimators over a broad range of luminosity. We confirm that the line widths of CIVλ1549 is not immediately offering a virial broadening estimator equivalent to the width of low-ionization lines. However, capitalizing on the results of Coatman et al. (2016) and Sulentic et al. (2017), we suggest a correction to FWHM CIVλ1549 for Eddington ratio and luminosity effects that, however, remains cumbersome to apply in practice. Intermediate ionization lines (IP ∼ 20–30 eV; AlIIIλ1860 and SiIII]λ1892) may provide a better virial broadening estimator for high redshift quasars, but larger samples are needed to assess their reliability. Ultimately, they may be associated with the broad-line region radius estimated from the photoionization method introduced by Negrete et al. (2013) to obtain black hole mass estimates independent from scaling laws.


Introduction: Statement of the Problem
A defining property of type-1 quasars is the presence of broad and narrow optical and UV lines emitted by ionic species over a wide range of ionization potentials (IPs, [1]): high-ionization lines (HILs) involving IP > 50 eV, and low-ionization lines (LILs) from ionic species with IP < 20 eV.Over the years, it has turned expedient to consider the UV resonance line CIVλ1549 as a representative of broad HILs.Lines do not all show the same profiles, and redshifts measured on different lines often show significant differences.Internal line shifts involving both broad lines in the optical and UV spectra of quasars have offered a powerful diagnostic tool of the quasar innermost structure since a few years after the discovery of quasars [2].
The CIVλ1549 is a resonant doublet ( 2 P o The occurrence of CIVλ1549 large shifts constrains the suitability of the CIVλ1549 profile broadening as a virial black hole mass (M BH ) estimator (see, e.g., [4,5], for reviews).Results at low-redshift obtained in the mid-2010s suggest that the CIVλ1549 line is unsuitable for, at least, part of the quasar Population A sources (Section 4 [6]).A similar conclusion was reached at z ≈ 2 on a sample of 15 high-luminosity quasars [7].More recent work tends to confirm that the CIVλ1549 line width is not straightforwardly related to virial broadening (e.g., [8]).However, the CIVλ1549 line is strong and observable up to z ≈ 6 with optical spectrometers.It is so highly desirable to have a consistent virial broadening estimate up to the highest redshifts that various attempts (e.g., [9]) have been done at "rehabilitating" CIVλ1549 line width estimators to bring them in agreement with the width of LILs such as Hβ and MgIIλ2800 [10].
In this contribution, we briefly stress the importance of black hole mass estimates (Section 2) and recapitulate the basic method of M BH estimates applied to large samples of quasars under the virial assumption (Section 3).To improve estimates that have been known to be plagued by systematic and statistical uncertainties as large as a factor ∼100, we adopt the "rehabilitating" power of the quasar eigenvector 1 (E1; Section 4).Our analysis is then focused on virial broadening estimators (VBEs; Section 5) from: (a) LILs (IP < 15 eV: Hβ, MgII 2800); (b) HILs (IP > 40 eV: CIVλ1549); (c) intermediate-ionization lines (IILs: SiIII]λ1892, Al IIIλ1860), for which we provide a brief summary of preliminary results.Reported results from our group were published during the last decade [11][12][13][14].We finally suggest that "photoioionization" computations of M BH (Section 6) may offer a solution to some of the problems associated with the use of an average scaling law.

Importance of Black Hole Mass Determination
The relevance of supermassive black hole M BH estimates is not limited to the sake of a better understanding of quasars interpreted as accreting system.Black hole masses are a key parameter in the evolution of the galaxies and in cosmology as well.Massive, fast outflows are affected by the ratio of radiation to gravitational forces.They provide feedback effects to the host galaxy, and are even invoked to account for the M BH -bulge velocity dispersion correlation [15][16][17].The ratio between radiation and gravitation forces also influences broad-line region dynamics; lower column density material may flow out of the emitting region [14,18,19].Black hole masses of high redshift quasars provide constraints on primordial black hole collapse (and references therein [20,21]).Overestimates by a factor as large as ∼100 for supermassive black holes at high z may even pose a spurious challenge to concordance cosmology.Present-day estimates constrain the relation between the formation of the seed black holes and the collapse of the protogalaxy.Collapse of dark matter clumps yielding massive seeds that then rapidly grows by super-Eddington accretion appears necessary to explain the occurrence of supermassive black holes at very high redshifts.

Virial Black Hole Mass Estimates
The virial expression for the mass can be written as where f is a factor of order unity reflecting geometry and dynamics of the broad line emitting regions, and r BLR is a representative distance of the broad-line emitting region (BLR). 1 The virial broadening δv is provided by a measurement of the line profile width, which can be velocity dispersion σ, FWHM, or FWZI.The FWHM is by far the most handy measure employed as a VBE although some authors claim that a better estimator may be offered by the velocity dispersion [22].The latter measure is, however, not defined in the case of Lorentzian profiles and not of easy interpretation in the case of shifted profiles.The virial assumption had a spectacular confirmation (albeit limited to few cases) by type-1 sources for which several lines were monitored: the response time of the BLR line and the line width were found to be correlated exactly as expected for Keplerian motion around a massive central object, with δv ∝ r − 1 2 [23].We can subdivide the estimates of the radius of the BLR (r BLR ) as primary and secondary.Primary determinations are obtained as the peak or centroid of the cross-correlation function between continuum and line light curves.Secondary determinations stem from the correlation between r BLR and luminosity that has been derived from reverberation-mapped sources (e.g., [24][25][26][27]): r BLR ∝ L a , a ≈ 0.5-0.65 [24,28].The relation takes different form for different lines, and has been defined for Hβ and CIV.It then becomes possible to derive scaling laws for the black hole mass: the relation is consistent with the virial assumption (e.g., [29,30]), and the r BLR scaling laws.If a = 0.5, b = 2 (e.g., [31][32][33]), the mass scaling law still provides M BH estimates, but the accuracy of these estimates is likely to be sample-dependent (see Section 5.3).
The M BH scaling laws provide a simple recipe usable with single-epoch spectra of large sample of quasars.Estimates of the Eddington ratio are derived by applying a bolometric correction to the observed luminosity.The correction is typically assumed a factor 10-13 from the flux λ f λ at 5100 Å (measured in erg s −1 ) and 3.4-5 from λ f λ measured at 1450 Å in the UV, [34,35]).For the sake of the present review, we will stay with these constant corrections without forgetting that different bolometric corrections should be defined along the quasar "main sequence" (introduced in Section 4), and that the correction is most likely luminosity as well as orientation dependent [36].

Caveats
The estimate of the M BH is based on several assumptions underlying reverberation mapping studies: among them, a compact region wetting the quasar continuum, and a fairly monotonic response of the line emitting gas.This latter assumption has been apparently challenged by the unpredicted behavior of NGC 5548 in 2014, where a time delay τ much shorter than expected was found when the source was in a high-luminosity state [37][38][39].The physical origin of this behavior is not yet clear: in principle, shielding, optically thin gas, changing size of the continuum source could also give rise to a shorter (or a lack of) response to continuum change in the emitting line regions.In addition, periodic signals (in photometric and spectroscopic measurements) have been detected in a number of sources [40], including NGC 5548 (e.g., [41][42][43][44][45]).The origin of the periodicity is as yet unclear as well.A second, supermassive binary black hole has been suggested in some cases [41,46] but, in principle, black holes of masses much smaller than the primary could produce periodic photometric properties without leaving a detectable trace of their gravitational pull on the line profiles.As a matter of fact, the r BLR -L scaling relation has a non-negligible intrinsic dispersion, which may be, at least in part, accounted for by a dependence of r BLR estimates on the dimensionless accretion rate found in recent work [27].
The basic assumption underlying the search of virial broadening estimators is that a line (or line component) is symmetric and unshifted with respect to the quasar rest frame.Most line profiles are asymmetric, and shifted, but this has been ignored until recent years.If CIV FWHM is employed as a virial broadening estimator, the loss of information is so severe that the M BH distribution as a function of redshift cannot be distinguished if masses obtained from random values of the FWHM (e.g., [47]) are used in place of the actual estimates.
A single value of the structure factor is obtained by scaling the M BH to agree with the dynamical masses [48][49][50][51][52].We expect that the geometry and dynamics of the BLR is related to the accretion mode, although it is as yet not unclear in which way.The structure factor is therefore likely to be different for different type-1 quasar populations.For Populations A and B (defined in Section 4), Collin et al. [53] find f (FWHM) ≈ 2 and 0.5, respectively.

The Rehabilitating Power of Eigenvector 1
The eigenvector 1 (E1) was originally defined by a principal component analysis of ≈80 Palomar-Green quasars, and is associated with an anti-correlation between the strength of the FeIIλ4570 blend, measured by the flux ratio R FeII = F(FeIIλ4570)/F(Hβ) and the FWHM of Hβ.The E1 is (at the very least) an useful tool to organize quasar diversity through a sequence (the main sequence, MS) in the so-called optical plane of the E1, defined by FWHM(Hβ) vs. R FeII .The 4D E1 parameters space includes two more parameters: the soft X-ray photon index Γ soft associated with the accretion status, and the centroid at half maximum of CIVλ1549 c(1/2) CIV measuring the amplitude of the high ionization outflow originating from the BLR [54].More parameters, however, correlate with E1 (see Table 1 of Sulentic et al. [55] and Fraix-Burnet et al. [56]).Since the analysis of Boroson and Green [57], MS trends have been found in works dealing with a number of phenomenological and physical correlates [58][59][60][61][62][63][64][65].Recently, large Sloan Digital Sky Survey (SDSS) data samples were involved [9,[66][67][68][69].The E1 "rehabilitating" power stems, in this context, by the ability to identify systematic trends that may be missed and confused as statistical errors if structurally different sources are dumped together in a single sample.
The trends along the MS allow for the definition of spectral types [70] in bins of ∆R FeII = 0.5 (A1, A2, A3, A4 from R FeII = 0 to R FeII = 2, tracing a change of L/L Edd convolved with orientation, [71]) and FWHM (A for FWHM Hβ < 4000 km s −1 , B1, B1 + , B1 ++ for broader sources in steps of 4000 km s −1 , tracing mainly a change in orientation).Quasars can be classified as belonging to two quasar Populations, A and B [72].Population A (FWHM Hβ < 4000 km/s) includes NLSy1s.Population A and B(roader) sources are most likely associated with a different accretion mode, since Population A are at L/L Edd > 0.2-0.3 which is the theoretical limit above and is an optically and geometrically thick advection-dominated accretion disk forms [73].Population B of low L/L Edd may be consistent with a flat α disk.The distinction between Populations A and B separates sources that are also called wind-and disk-dominated by Richards et al. [66], or Populations 1 and 2 by Collin et al. [53].It supersedes the distinction between NLSy1s and rest of type-1 AGNs, which is based on a FWHM limit that has played an important historical role but is physically arbitrary.At low z (<0.7),Population A show low M BH , high L/L Edd Population B, high M BH , low L/L Edd , a reflection of the "downsizing" of nuclear activity: practically no black hole with very large mass (M BH ∼ 10 9 -10 10 M ) is accreting super-Eddington in the local Universe.

LILs: Hβ and MgIIλ2800
The profiles of LILs like Hβ and MgIIλ2800 change along the quasar MS.If we want to extract a VBE, we have to consider the behavior of Hβ and MgIIλ2800 from extreme Population B to extreme Population A.
In Population B, the broad profiles of Hβ and MgIIλ2800 are most frequently redward asymmetric.Composite spectra of individual spectral types can be modeled by a broad Gaussian component (the BC, symmetric and unshifted, assumed to be the virialized component), and a redshifted very-broad Gaussian component.This latter component presumably associated with "perturbed" virialized motions or gravitational redshift in the inner BLR [74] systematically increases the line width with respect to the BC that provides the VBE.
At the other end of the MS extreme, Population A sources show narrower, Lorentzian-like profiles, slightly blueward asymmetric.The Hβ and MgII profiles can be modeled by Lorentzian functions plus a blueshifted excess modeled with a skewed Gaussian [12,13].Both the LIL resonance line MgIIλ2800 and Hβ are affected by low-ionization outflows detected in the extreme Population A corresponding to spectral types A3-A4 [13].However, in most of Population A, the LILs are dominated by a symmetric, virialized broad component [8,13,30,75].
We can define a parameter ξ yielding a correction to the observed profile to obtain a VBE, as follows: where the VBE FWHM can be considered best estimated by the FWHM of the broad component FWHM BC , whose proxy can be obtained by symmetrizing the observed FWHM by various corrections described below.Even if the non-virial broadening mechanism is different along the E1 sequence, 0.75 ≤ ξ ≤ 1.0 for both Hβ and MgII.This implies that, even if LILs are not always asymmetric, a modest correction is, on average, sufficient to retrieve a VBE from the observed FWHM.Table 1 reports the current estimates along with bibliographic references.

LIL VBE at High L
At low z, the quasar population reaches bolometric luminosity ∼10 47 erg s −1 at most.It is interesting to consider the LIL VBE behavior if we add to a local sample sources more luminous than this limit.Composite spectra are helpful to outline the LIL Hβ behavior over a wide luminosity range.Marziani et al. [11] computed composites in step of 1 dex for 6 dex in luminosity, joining an HE/ISAAC high-L sample (52 sources) and an SDSS sample from Zamfir et al. [77].The Hβ line becomes broader with increasing L (over 43 < log L < 48.5 [erg/s]), but shapes are similar: the Population A/Population B differences are preserved at high L [11].At the same time, the minimum FWHM(Hβ) increases with luminosity.This result is consistent with the virial assumption and with the r BLR scaling law: The minimum, luminosity-dependent FWHM is obtained for a limiting Eddington ratio ≈1.As a consequence, the Population A limit is also luminosity dependent.
Extracting a VBE from LIL Hβ at high-L is possible by employing different symmetrization methods: • substitution of the BC in place of the full Hβ profile.This method requires a multicomponent maximum likelihood fitting, and is fairly reliable for Population A (even for extreme Population A) and Population B spectral types where the VBC is creating a profile inflection; • symmetrization of the profile: FWHM symm = FWHM−2 c(1/2); • correction based on spectral type.In practice, this means to correct Hβ for Population B sources by a factor ξ as reported in Table 1.
All symmetrization methods were found to be equivalent at low-and high-L [45].In other words, the Hβ profile shapes at high L are consistent with those at low-z, lower L (with some caveats, [11]).

CIV
Virial broadening estimators extracted from the HIL CIVλ1549 along the E1 sequence are unfortunately not easy to define.If we scale the Hβ emission and overlay it on CIVλ1549 (Figure 1), we obtain a strong excess of blueshifted emission.The profile can be interpreted as an almost symmetric, unshifted "virialized" emitting region + an outflow/wind component that dominates in A3/A4 spectral types (e.g., [78][79][80][81]) and at high luminosity.The largest shifts of CIVλ1549 centroid at 1/2 along the MS are found in Population A [4,6,14].The Vestergaard and Peterson [29] scaling laws assumed that the width of CIV and Hβ are equivalent, with a ≈ 0.5, b = 2, and c an average constant difference between the 5100 and the 1450 luminosity.Considering the FWHM CIV as a proxy of FWHM Hβ causes a bias along the E1 sequence.Especially for extreme Population A, errors can be as large as 2 dex. Figure 2 is the statement of the ensuing BH mass taboo for CIV. Figure 2 (based on the data of Sulentic et al. [6]) emphasizes a systematic trend for which a corrective could be estimated.To this aim, we defined a sample of RQ quasars with CIVλ1549 observations including 70 sources at z < 1 [6], and 25 high-luminosity sources at z > 1.4 [14].Matching Hβ data are available from our previous observations [11,82] or from published spectra.At L > 10 47 erg s −1 , high amplitude CIV 1549 blueshifts in both Population A and B are observed with median ≈ 3000 km s −1 for Population A; two extreme cases involve CIV c(1/2) blueshift amplitude larger than 5000 km s −1 .Widespread powerful outflows are affecting both Population A and B sources [14,83] with worrisome implications for M BH estimates from CIVλ1549 FWHM.
Population A sources are more frequently selected at high z, high L. Figure 2 of [84] depicts what may be called an Eddington ratio bias: for a fixed mass, higher Eddington ratio sources are preferentially better sampled.In this way, the FWHM CIV (if left uncorrected) may lead to systematic over-estimates of M BH by even one/two orders of magnitude, with the potential of creating the (spurious) result of a population of extremely massive black holes at high redshift (z > 2).

A CIVλ1549 VBE
The centroid shift at half-maximum c(1/2) is correlated with the FWHM in the CIV line, implying that the the FWHM excess with respect to Hβ is associated with a blueshifted component, as shown in Figure 1. 2 The comparison of Figure 1 has been possible because we have Hβ observations for this source.However, Hβ observations are available only for a few hundred high-redshift quasars since the line is shifted in the IR.
Nonetheless, there have been several attempts in the last few years to define scaling laws that corrected for the non/virial contamination of CIV emission.A scaling law that assumes M BH ∝ FWHM 0.5 [32] accounts for the over-broadening of Population A sources, but overcorrects for Population B [45].In general, corrections dependent on L/L Edd (which is correlated with FWHM of CIV, [85]) or an L/L Edd proxy such as the SiIV + OIV]1400 blend/CIV 1549 ratio are promising [9] but tend not to work well for Population B. Empirical corrections based on CIV blueshift work fairly well for high-L sources only [83,86].
It is important to consider that there is a threshold in CIV shift amplitude (c(1/2)) at L/L Edd ≈ 0.2 [14].There is a strong correlation with L/L Edd if blueshifts are significant.There is also a weak but significant correlation with luminosity in the sample of Sulentic et al. [14]: the partial 2 This happens also for MgII, but with much lower blueward displacements (few hundreds vs. few thousands km s −1 ).correlation coefficient between c(1/2) and L (L/L Edd hidden) is significant at about a 2σ confidence level.A multivariate analysis confirms the blueshift dependence on both L and L/L Edd in that sample.The blueshift dependences are consistent with a radiation-driven outflow, with a slope ≈0.15 for log L, and slope ≈0.5 for L/L Edd .A pure dependence on L arises for L/L Edd in a small range.A strong dependence on L/L Edd and a weak dependence on L can be obtained under a variety of physical scenarios.We expect that strength and form of L and L/L Edd correlations are sample dependent [14].
Considering the threshold and the dependence on both L and L/L Edd , using c(1/2) as a proxy for L/L Edd , HIL CIV corrections based on c(1/2) and L reduce scatter to 0.33 dex with respect to Hβ estimates, providing an unbiased M BH estimator.Unfortunately, the correction coefficients are different for Pops.A and B. For Population B, the correction is highly uncertain in the sample of Sulentic et al. [14], and a larger sample of sources is needed.In addition, (1) corrections derived by Coatman et al. [83] may be sample dependent because of the Eddington ratio bias; and (2) any correction based on c(1/2) requires a precise estimate of the quasar rest frame of reference, which is not easily obtained from UV data.A theoretical correction requires that c(1/2) CIV and ionization conditions in the BLR are accounted for, and has not been competed as yet.It is therefore not obvious whether M BH estimates based on CIV can be reliable for large samples of quasars.

Intermediate Ionization Lines (IILs)
An alternative to the use of CIV is to resort to the emission lines of the 1900 Å blend, whose constituents are mainly the resonant doublet of AlIII at λ1860 Å, and the intercombination lines SiIII] and CIII] at λ 1892 and λ1909 Å, respectively (Table 1 of Negrete et al. [87] provides detailed information of the atomic transitions involved).We base our analysis on ≈80 objects of the CIV sample described in Section 5.3 for which the 1900 blend has been covered.CIII]λ1909 measurements have serious problems: in Population A, CIII] is faint, and blended with strong and diffuse FeIII emission.In Population B, the CIII]λ1909 line is affected by VBC.However, considering AlIII and SiIII] it is possible to obtain a VBE consistent with Hβ.To this aim is necessary to assume that FWHM AlII 1860 = FWHM SiIII] 1892 i.e., to anchor the FWHM AlII to FWHM SiIII].A source of concern is that AlIII is a resonance doublet ( 2 P o ) and part of its emission may originate in an outflow, where the AlIII, as does CIV, acts as a resonant scatterer of continuum photons.However, measured AlIII shifts are < 0.2 CIV shifts and the AlIII and SiIII] profiles looks fairly symmetric in the wide majority of cases (Figure 1 shows a typical case).
If we compare the FWHM of the IILs to Hβ, we find that it is in very good agreement with Population A. Population B IILs are narrower than Hβ, but a modest correction is needed to enforce an unbiased agreement: FWHM Hβ BC = 1.35 FWHM AlIII BC .In other words, ξ ≈ 1.35 (Table 1) is needed to rescale the line widths to the one of Hβ BC assumed as a reference VBE.

Photoionization Masses
The ionization parameter U can be written as where Q(H) is the number of ionizing photons, n the hydrogen density, and c the speed of light.The expression for U can be inverted to obtain the BLR radius The radius r BLR estimates from photo-ionization agree with cτ from reverberation mapping for a sample of 12 sources [76].The photon flux n • U is estimated using diagnostic ratios involving AlIII 1860, SiIII] 1892, SiII 1816, CIV 1549, SiIV + OIV] 1400.The photoionization method provides an unbiased estimator of r BLR in the sample of Negrete et al. [76], but M BH estimates at high L remain largely untested [88].The photoionization method, in principle, could avoid the dispersion intrinsic to the scaling laws with L. To obtain accurate individual estimates of M BH is then necessary to have an accurate knowledge of f and of orientation effects that influence the FWHM.Both are still poorly known at the time of writing.

Conclusions
Retrieving a VBE that is representative of the broadening due virial motions in the low-ionization part of the BLR is a major goal in order to reduce systematic and statistical effects in single-epoch M BH determinations applied to large samples of quasars.The following remarks emerge from the present review.

•
Low-ionization lines (Hβ, MgII 2800) provide reliable virial broadening estimators by applying corrections to the observed line width.The corrections depend on the spectral type along the E1 MS, but they are relatively small (less than 30%), and work up to the highest L of quasars.

•
The HIL CIVλ1549 is not immediately providing a reliable virial broadening estimator.The profile is broadened by an excess emission on its blue side.The shift amplitude depends on both L/L Edd and L. Large shifts are observed in Population A, with Eddington ratio above a critical L/L Edd ≈ 0.2.

•
It is possible to apply corrections to the observed CIV broadening, but they remain cumbersome even for Population A. Population B sources at low Eddington ratio require a different correction (still ill-defined by the analysis of Marziani et al. [45] as Population B sources are most affected by the Eddington ratio bias mentioned in Section 5.3).

•
Preliminary results on the 1900 blend indicate that the IILs lines could provide a better choice than CIV; IIL FWHM measurements appear intrinsically more robust than those of CIV since they do not require corrections based on shift measurement with respect to rest frame.However, more data are needed to assess their reliability.
The ultimate solution may be to abandon scaling laws altogether and to attempt M BH estimates on an individual basis, considering r BLR from photoionization and f = f (L/L Edd , L, . ..) (where . . .indicates any relevant physical parameter, for example the spin parameter of the black hole) as well as orientation effects on the VBE along the E1 sequence.

Figure 1 .
Figure 1.Left: the CIV profile of HE0512, a Population A quasar in the sample of Sulentic et al. [14].The profile has been decomposed into a symmetric, unshifted component (the virialized component, thick black line) and a blue shifted component modeled by a skewed Gaussian (thick blue line).The flat-topped emission on the red side of CIV is mainly HeIIλ1640; Right: the 1900 blend for HE0512.The AlIII and SiIII] lines at 1860 Å and 1892 Å are unshifted and symmetric, and their FWHMs are in agreement with the Hβ and of the CIV virialized component profile.The thick black components can be thought as scaled Hβ profiles, as in the case of CIV.Note that AlIIIλ1860 lacks the blue-shifted excess of CIV.The thin line shows the CIII] λ1909 profile whose intensity remains highly uncertain because of the heavy blending with FeIII emission.

Figure 2 .
Figure 2. Bias of the M BH estimates in the optical E1 plane FWHM(Hβ) vs. R FeII , for the full sample of Sulentic et al. [6].Grey symbols represent differences |∆log M| = | log M(CIV) − log M(Hβ)| ≤ 0.3, blue dots ∆log M > 0.3, red circles ∆log M < −0.3.The size of the circles is proportional to the amplitude of the difference, with the largest blue circles showing CIV masses in excess by a factor 100. Mid-sized circles indicate a disagreement of a factor ≈10.The labels identify the spectral types as described in Section 4 and the thick dot-dashed line at 4000 km s −1 shows the boundary between Population A and Population B.

Table 1 .
The LIL and IIL ξ (Equation2) factor for the spectral types of the quasar MS, listed for the spectral types (SpT) of highest occupation.