Metal Content in Relativistically Jetted and Radio-Quiet Quasars in the Main Sequence Context

: The optical and UV properties of radio-quiet (RQ) and radio-loud (RL, relativistically “jetted”) active galactic nuclei (AGN) are known to differ markedly; however, it is still unclear what is due to a sample selection and what is associated with intrinsic differences in the inner workings of their emitting regions. Chemical composition is an important parameter related to the trends of the quasar main sequence. Recent works suggest that in addition to physical properties such as density, column density, and ionization level, strong Fe II emitters require very high metal content. Little is known, however, about the chemical composition of jetted radio-loud sources. In this short note, we present a pilot analysis of the chemical composition of low-z radio-loud and radio-quiet quasars. Optical and UV spectra from ground and space were combined to allow for precise measurements of metallicity-sensitive diagnostic ratios. The comparison between radio-quiet and radio-loud was carried out for sources in the same domain of the Eigenvector 1/main sequence parameter space. Arrays of dedicated photo-ionization simulations with the input of appropriate spectral energy distributions indicate that metallicity is sub-solar for RL AGN, and slightly sub-solar or around solar for RQ AGN. The metal content of the broad line-emitting region likely reﬂects a similar enrichment story for both classes of AGN not involving recent circum-nuclear or nuclear starbursts.


Introduction
Type-1 active galactic nuclei (AGN) are characterized by the presence of broad and narrow optical and UV lines (for an introduction, see e.g., [1][2][3]).AGN spectra show an overwhelming variety of broad emission line profiles not only among different objects but also among different lines in the spectrum of the same object.For stars, the identification of optical spectral types and luminosity classes allows for the knowledge of the most relevant physical parameters, including the star's evolutionary status [4].This feat is not yet possible for AGN.However, progress in empirical classification has yielded a coarse contextualization of the main accretion properties, such as black hole mass, Eddington ratio, outflow prominence, spectral energy distributions, emitting region size, the metal content of the line-emitting gas, etc. (e.g., [5][6][7][8][9][10][11]).The main set of correlations was derived from a principal component analysis (PCA) of a sample of several tens of quasars [12].The importance of Eigenvector 1 derived from the PCA has revealed itself over the years [13], leading to the definition of what has become known as the main sequence of quasars [14,15].The optical plane of this main sequence is identified by the line width of the HI Balmer line Hβ (FWHM Hβ) and the prominence of a singly ionized emission, defined as the flux ratio between the FeII blends centered at λ4570 and Hβ itself (hereafter R FeII ).The distribution of type-1 unobscured AGN takes the form of an elbow-shaped sequence in the optical plane (see e.g., Figure 1).The ranges of R FeII and FWHM Hβ define spectral types [15,16], and different classes may show different occupations in the plane (Figure 1).A case in point is provided by RQ and RL [17]: most powerful, relativistically jetted sources cluster in Population B, with moderate FeII emission and relatively broad Hβ line profiles (4000 km s −1  FWHM(Hβ) 8000 km s −1 ).The broader spectral type shows higher fractions of RL; however, these spectral bins have a low prevalence, and a sizeable fraction of all type-1 AGN falls only in bin B1.This spectral type is, therefore, well suited for an inter-comparison between RQ and RL properties.In the following, we will build composite spectra for jetted and non-jetted sources of spectral type B1 (Section 2).The analysis takes advantage of the main sequence (MS) correlations concerning line profiles in the spectral type B1 (Section 3).It is focused on the measurements of intensity ratios intended to be diagnostics of the metal content of the line-emitting gas.Intensity ratios are interpreted using arrays of photo-ionization simulations covering a broad range of metallicity.Results (Section 4) suggest slightly subsolar or solar metallicity and provide evidence of significant diversity in terms of chemical enrichment and evolutionary status along the main sequence (Section 5).
Figure 1.Sketch of the optical plane of the MS, FWHM Hβ vs. R FeII .Numbers in square brackets provide, from top to bottom, the complete sample fractional occupation, the fraction of jetted (coredominated + Fanaroff-Riley II) sources in each spectral bin.The shaded area identifies the region along the sequence with a high prevalence of jetted sources (c.f.Ref. [18]).The original sample is described in Ref. [19].Points mark the location of the composite spectra (red: RL; blue: RQ) and of PKS 0226-038 (black).

Composite Spectra for Spectral Type B1 Quasars
The obtaining of spectral data of good quality covering the full spectral ranges from 1000 Å to ≈6000 Å (i.e., from Lyα to Hβ included) for the same object is still a non-trivial feat for low-z quasars, as the UV coverage demands space-based observations.Here, we consider one sample of 20 RL and one sample of 16 RQ sources, all belonging to spectral type B1, in the redshift range ≈0.002-0.5 and ≈0.25-0.65,for RQ and RL, respectively.Based on the absolute magnitudes derived from the quick-look magnitudes in the NASA Extragalactic Database (NED), the range is between −21 and −27 (RQ), and be-tween −23.5 and −26.5 (RL), which correspond to bolometric luminosities in the range log L ∼ 45-47 [erg s −1 ].The UV data are HST/FOS observations analyzed in [20], and the optical spectra were obtained from [19].The spectral similarity ensures the consistency of black hole mass and Eddington ratio.We built median and average composites for the RQ and RL classes that look consistent.As the median combination is not improving the S/N, fitting analysis was performed only on the averages.The average S/N is high with S/N ≈ 90 for RQ composites, and ≈130 and ≈55 for the visual and UV ranges in the RL composite, respectively.The composite spectra are shown in Figure 2. The RQ and RL composite spectra covering the optical and the UV domain.The green shaded area marks the range used to measure the FeIIλ4570 blend.The grey-shaded region has not been considered in the analysis due to the limited number of available spectra in this region.The RL spectrum has been vertically shifted for clarity.The dot-dashed lines identify the main emission feature included in the non-linear χ 2 fitting procedures.

Line Profiles
Population B Hβ profiles show a redward asymmetry modeled with a broader redshifted (FWHM ∼10,000 km s −1 , shift at line base ∼1000-2000 km s −1 ) and a narrower Gaussian [16,21].The very broad Gaussian component is meant to represent the innermost part of the broad line region (BLR), providing a simple representation of the radial stratification of the BLR in Pop.B suggested by reverberation mapping [22].This component (hereafter the very broad component, VBC) has been associated with a physical region of high-ionization virialized gas and closest to the continuum source-the very broad line region, VBLR [23][24][25][26].Although the physical properties of the VBLR line-emitting gas are not well known, a decomposition of the full Hβ profile into a symmetric, unshifted Hβ component (Hβ BC ) and a Hβ VBC provide an excellent fit to most Hβ Pop.B profiles [16,27].
The multicomponent fits were performed using the SPECFIT routine from IRAF [28].This routine allows for the simultaneous minimum-χ 2 fit of the continuum approximated by a power law and the spectral line components yielding FWHM, peak wavelength, and intensity for all line components.In the optical range, we fit the Hβ profile and the [OIII]λλ4959,5007 emission lines, and the FeII multiplets for the composite objects.

Diagnostics of Metallicity and Photo-Ionization Modeling
Diagnostics from the rest-frame UV spectrum take advantage of the observations of strong resonance lines that are collisionally excited [29,30] and at least constrain density n H , ionization parameter U, and chemical abundance Z.For instance, C IVλ1549/Lyα, C IVλ1549/(Si IV + O IV])λ1400, C IVλ1549/He IIλ1640, N Vλ1240/He IIλ1640 are sensitive to metallicity; and Al IIIλ1860/Si III]λ1892, Si III]λ1892/C III]λ1909 are sensitive to density, since inter-combination lines have a well defined critical density [29,31].Ratios of lines involving different ionic stages of the same element are obviously sensitive to the ionization parameter.The lines emitted from ionic species of silicon and aluminum deserve special attention-they are two elements greatly enhanced in supernova ejecta [32].This approach has yielded tight constraints, especially for sources radiating at high Eddington ratio [29,33,34], where physical properties are consequently well constrained because they converge toward an extreme.
The photo-ionization code Cloudy [35] models the ionization, chemical, and thermal state of gas exposed to a radiation field and predicts its emission spectra and physical parameters.Cloudy simulations require inputs in terms of n H , U, Z, quasar spectral energy distribution (SED), and column density N c .The ionization parameter L ν /hν is the number of ionizing photons, provides the ratio between photon and hydrogen number density and is dependent on the spectral energy distribution of the ionizing continuum.The simulations were carried out assuming the RL and RQ SEDs from Laor et al. [36] representative of the SED of B1 objects (work in preparation).The geometry was assumed to be open and plane-parallel, meaning that a slab of emitting gas is exposed to a radiation field only on one side.Arrays of Cloudy photo-ionization models ( 1 ) for a given metallicity Z and N c , constant density n and U were evaluated at steps of 0.25 dex covering the ranges 7 ≤ log n H ≤ 13 [cm −3 ], −3 ≤ log U ≤ 1.The single-value metallicity arrays were computed for log Z at −3.0, −2.7, −2.3, −2.0, −1.7, −1.3, −1.0, −0.7, −0.3, 0, 0.3., 0.7, 1.0, and 1.3 in solar units, i.e., from 0.001 Z to 20 Z .No dust and no microturbulence broadening were included in the calculations.The Z calculations are based on a single zone assumption for the BLR.The lowest values of the density may bias the solutions toward cases where significant [OIII]λλ4959,5007 is expected; all cases with [OIII]λ5007/Hβ > 0.1 were excluded, as no broad [OIII]λλ4959,5007 is observed.A more refined analysis in the framework of the locally optimized emitting cloud model [37][38][39] is deferred to an eventual work.

Composite Spectra
The composite average spectra are shown in Figure 2. The spectral similarity between the two classes already evident in Figure 2 is confirmed by the profile comparison of Figure 3.Both spectra are FeII weak; show weak NV; AlIII is weak and the CIII] line is by far the most prominent in the 1900 Å blend encompassing SiIII] and CIII] along with AlIII.All lines can be successfully decomposed into a BC and a VBC.Apparently, there is no VBC in the AlIII and SiIII] lines and in FeII Ṅo prominent BC is observed in the HeII optical and UV lines, at variance with Hβ, Lyα, CIV, and the other lines.The absence of a prominent BC in helium lines is typical [40], albeit the implications for the BLR structure are not obvious even on a qualitative basis.Table 1 reports the line fluxes normalized to Hβ.In view of the heuristic decomposition approach, fluxes are reported for the BC, VBC, and sum of the two components, i.e., for the full broad profile (total flux).Each value in this table has been normalized by the corresponding component of Hβ.Error analysis has been performed empirically or on the basis of previous analyses (Hβ and FeII: [19,41]; CIV: [42]; blend at 1900 Å: [34,41]; NV: [33]).The undetected HeII BC was considered in the models to be HeII/Hβ 0.1.A more refined treatment of uncertainties should be considered in future work.Spectral differences between the two classes are not striking on a qualitative basis.We note a stronger redward asymmetry in RL sources.This is, however, a result known for a long time [43,44].The optical FeII emission is twice as strong in RQ than in RL; this is also consistent with previous results [18].The intensity ratios most commonly used in Z estimates are marginally higher for RQ than for RL.SiIV+OIV]λ1402/CIV ≈ 0.24 vs. 0.15 for RQ and RL, respectively.Similarly, AlIII/CIII] ≈ 0.2 vs. 0.12; NV/Hβ ≈ 1.32 vs. 1.05, if total fluxes are considered.
We computed the normalized χ 2 ν in the following form, to constrain the value of the metallicity from the values of the diagnostic ratios.For each spectrum k, and for each component c, we can write (following [33]): where the summation is made over the available diagnostic ratios.The number of degrees of freedom n f is 8 for the composites and 7 for the case of PKS 0226-038 presented in Section 4.2.The χ 2 ν is computed to the results of the Cloudy simulations as a function of U, n H , and Z (subscript 'mod').The following independent intensity ratios were considered in the summation: NV/CIV, CIV/HeIIλ1640, CIV/Hβ, SiIV+OIV]λ1402/CIV, AlIII/CIV, AlIII/SiIII], AlIII/CIII], HeII/Hβ, FeIIλ4570/Hβ.The ratios were computed for BC, VBC, and total line intensity.Non-detections were treated as upper limits.
Figure 4 shows the projections of the 3D space U, n H , Z.Each contour in the plane braces elements of the grid of Cloudy parameter space that is consistent with the minimum χ 2 ν within the uncertainties at 1σ, 2σ, and 3σ confidence levels.The distribution of the data points is constrained in a range of U, n H , Z, at very low density, relatively high ionization, and low metallicity.Within the limit in U, and n H , the distribution of Z is slightly subsolar for the RQ (Z ≈ 0.5Z ), and even more subsolar for the radio-loud (Z ≈ 0.1Z ).There is a relatively broad range of density and ionization parameters and metallicities that are consistent within 2σ confidence level from the minimum χ 2 .This is probably a reflection of the stratified nature of the emitting region in Population B and the intrinsic heterogeneity of the composites.
The values reported in Table 2 for the sub-regions confirm the validity of the virial scenario, with U increasing by a factor 10 between BLR and VBLR.The derived density is not changing strongly, as implied by the strong CIII] VBC.Since the distance from the center of gravity scales with the inverse of the velocity dispersion squared (i.e., r ∝ FWHM −2 ), we might expect that U/U BLR ∼ FWHM 4  VBLR /FWHM  The Z values all fall in ranges that are consistent.The results for the sub-regions should be read with some care, as the VBC measurements are difficult for HeII, HeIIλ1640, NV.The fairly high metallicity value for RQ VBLR is due to several line ratios being consistently higher than for the RL case: CIV/Hβ, (SiIV+OIV]λ1402)/CIV, (SiIV+OIV]λ1402)/HeIIλ1640, and NV/HeIIλ1640 are between 30% and a factor two times higher for the VBLR of RQ than of RL.This epitomizes the need for very accurate measurements on the line profiles.
These results can be compared with the ones derived from the analysis of the Sloan composite spectrum [45].Values of the intensity ratios reported in their Table 2, and the assumption of typical 10% uncertainties for most ratios do not provide any constraint on Z, and loose constraints on density (logn H ∼ 9.5 +1.5 −2.5 ) and ionization (log U ∼ −1.75-−0.5).Averaging the full dataset of a large color-based spectroscopic survey presents major hindrances.Spectral bins are likely to be associated with intrinsic differences in Z. Building a composite combines all survey spectra with a weight proportional to the relative prevalence of each spectral type (e.g., [13]).Since the most populated spectral bin along the sequence is B1, with A2 the second most populated [15,46], the appearance of the composite spectrum resembling the Pop.B RQ composite has some key features, and A2 in others.For example, the moderate FeII/Hβ and SiIV+OIV]λ1402)/CIV values are consistent with solar or supersolar Z.The composite yields line ratios that reflect a combination of physical properties from individual spectra; on the other hand, it does not reflect the spectral diversity of a large population of type-1 AGN, as extreme sources are just a minority ( 5%; Figure 1).

A Typical RL Source at High z
PKS 0226-038 is a luminous jetted source at the cosmic noon (z ≈ 2.06922; M B ≈ −28.0, and Kellermann's radio loudness parameter [47], R K ∼ 10 3 ).Its optical and UV spectrum is markedly different from the Population B composite at low z: CIII] is weaker, FeII stronger, and AlIII is almost as strong as CIV (Figure 5).No diagnostic based on NV and HeII is available in this case.The measured diagnostic ratios are reported in Table 3.The uncertainties were assumed to be ≈10% for the intensities normalized to Hβ save for CIV/Hβ which is affected by a re-scaling (≈20%) and for FeIIλ4570 for which δR FeII = 0.1 was assumed.In this case, the FeIIλ4570 (Figure 5) is partly missing because of a telluric absorption.The FeII optical emission is well-represented by a "solid" template with fixed multiplet ratios.The derived values of the metallicity are well constrained with minimum χ 2 obtained for log Z/Z ≈ 1, with 1σ uncertainty equal to 0.3 dex.The ionization parameter log U ≈ −2.75 and density logn H ∼ 11.75 [cm −3 ] are significantly lower and higher with respect to the values at low-z.The three parameters reflect the moderate R FeII , high AlIII/CIV (yielding low ionization and high Z), and low CIII]/SiIII] (implying a high density).

Discussion: Metal Enrichment along the Quasar Main Sequence
Several basic results emerge from the present analysis: 1.
RQ Pop.B sources show slightly subsolar or solar metallicities; 2.
in the same spectral bin along the MS, RL sources show definitely subsolar chemical abundances, lower with respect to RQ.The difference is not marginal and is supported by several diagnostic indicators consistently observed to be lower in RL sources.It is also consistent with the location of the RL template in the MS, displaced toward broader Hβ and lower R FeII with respect to the RQ one; 3.
there is an important difference in properties at the extremes of the E1 MS: here, we focused only on the low-R FeII extreme deriving subsolar metallicity.However, highly supersolar Z estimates were derived at the other end of the MS, with R FeII 1 [33,34].
The first result is not completely unexpected, since a similar analysis yielded slightly subsolar abundances for NGC 1275 [48] and Pop.B sources are known to have similar intensity ratios [49].
The origin of the RQ and RL Z difference is likely associated with a difference in sample/host galaxy evolutionary history.At the modest accretion rates required to sustain the low Eddington ratio of B1 sources, there is no need to invoke "cataclysmic" events such as recent, major wet mergers to explain the AGN luminosity via super-Eddington accretion rates [50], and accretion material could be provided via stellar mass loss in earlytype galaxies [51].The general impression of low-z Population B is that of fairly evolved systems, past their prime of accretion events [52].This impression is further reinforced in the case of radio-loud sources, where sub-parsec binary black holes at the stable end of their in-spiral phase are relatively frequent [53].
The gradient in chemical composition is also no surprise.Although it is possible to account for an increase in FeII along the sequence only based on ionization degree and density [54], a thorough quantitative assessment also requires a change in chemical composition [9].We pass from sources where singly ionized iron emission is barely detected to sources where it dominates the thermal balance of the emitting regions, overwhelming the entire Balmer line emission.
The case of the high-z RL source PKS 0226-038 signifies the co-existence of powerful relativistic ejections and BLR physical conditions typical of sources radiating at L/L Edd 0.1-0.2(Population A according to [14]).It has been known for a long time that the radio-loud fraction among AGN increases with redshift [55,56].However, the realization of the high prevalence of unambiguously jetted AGN at z 1 radiating at high L/L Edd is a recent result.The Eddington ratio of PKS 0226-038 is logL/L Edd ≈ −0.43.It belongs to the spectral type B2, which include sources in the R FeII range 0.5-1, and with moderate accretion rate.Strong FeII emission is a rather rare occurrence among RLs at low-z (especially in Fanaroff-Riley II sources, [17]), it is frequently found at high redshift (Ref.[57], and an article in preparation), along with higher L/L Edd .This is probably a consequence of the increase in average accretion rate at the epochs corresponding to the cosmic noon [58][59][60], and much beyond.Indeed, highly accreting jetted sources are being discovered up to redshift ≈ 6 [61,62].

Conclusions
This pilot work revealed a population of quasars with slightly subsolar abundances.We can infer a metallicity trend along the quasar MS on the basis of previous results and of the analysis presented in this paper, supporting the assumptions in the study of Panda et al. [9], although the MS-oriented analysis has been so far focused only on spectral types, B1 and extreme Population A. Full analysis of individual spectral types isolated along the MS [16] is still needed.Special attention should also be devoted to high-z high-L sources.At high-z, the sample selection effects could be even more important, excluding a population of sources on the basis of their low accretion rate [63].At high luminosities, wind effects are expected to be maximized, and may introduce another factor in the chemical enrichment of the broad line-emitting gas.

Figure 2 .
Figure 2.The RQ and RL composite spectra covering the optical and the UV domain.The green shaded area marks the range used to measure the FeIIλ4570 blend.The grey-shaded region has not been considered in the analysis due to the limited number of available spectra in this region.The RL spectrum has been vertically shifted for clarity.The dot-dashed lines identify the main emission feature included in the non-linear χ 2 fitting procedures.

Figure 3 .
Figure 3. Analysis of prominent blends helpful for metallicity diagnostics.Continuum-subtracted spectra are shown in the rest frame, with the blend ordered with an increasing wavelength from left to right: the Lyα + NV blend, SiIV+ OIV]λ1402; CIV + HeII; the 1900 Å blend made up mainly of AlIII, SiIII], CIII]; the Hβ + [OIII]λλ4959,5007 region.Dashed magenta line: model spectrum; thick solid line: broad components; red thick line: very broad components; thin orange lines: narrow components of Hβ and [OIII]λ5007; blue lines: blue-shifted components.Green lines trace the scaled and broadened FeII emission template.The lower panels show the observed minus model residuals in the radial velocity scale.

Figure 4 .
Figure 4. Projections of the 3D parameter space n H , U, Z. Isophotal contour traces the 1,2,3σ confidence levels; the green spot identifies the values of n H , U, Z yielding the minimum χ 2 ν .The upper panels are for RQ, lower panels are for RL.Regions of the parameter space where [OIII]λ5007/Hβ > 0.1 have been assigned a very high χ 2 ν .

Figure 5 .
Figure 5.The spectrum of PKS 0226-038 in the optical and the UV domain.The dot-dashed lines identify the main emission features included in the non-linear χ 2 fitting procedures.The straight line traces a range where data are unavailable for this source.

Author Contributions:
Investigation, writing, P.M.; computations, S.P.; data analysis, A.D.M.; validation, A.d.O.All authors have read and agreed to the published version of the manuscript.
4BLR ∼ (2.32 − 2.55)4∼ 29 − 42.The U values reported in Table2indicate ratios between one and two orders of magnitudes, consistent with the ones expected from a virial velocity field.

Table 2 .
Derived values of U, Z, n H and 1σ ranges a .H in units of cm −3 .The ranges are defined by the limiting elements of the model grid that are compatible with the minimum χ 2 a : ionization parameter U, abundance Z in solar units, and hydrogen particle density n ν within 1 σ confidence level.

Table 3 .
Diagnostic line ratios measured on the PKS 0226-038 spectrum a .Line ratios refer to the total flux of the full broad line profiles. a