Powers and Magnetization of Blazar Jets

In this work I review the observational constraints imposed on the energetics and magnetisation of quasar jets, in the context of theoretical expectations. The discussion is focused on issues regarding the jet production efficiency, matter content, and particle acceleration. I show that if the ratio of electron-positron-pairs to protons is of order $15$, as is required to achieve agreement between jet powers computed using blazar spectral fits and those computed using radio-lobe calorimetry, the magnetization of blazar jets in flat-spectrum-radio-quasars (FSRQ) must be significant. This result favors the reconnection mechanism for particle acceleration and explains the large Compton-dominance of blazar spectra that is often observed, without the need to postulate very low jet magnetization.


Introduction
One of the biggest unresolved puzzles in the theory of active galactic nuclei (AGN) regards their ability to launch very powerful, relativistic jets. According to several studies, the most powerful jets reach or even exceed the associated accretion disk luminosities [1][2][3][4][5]. Obviously, production of such jets cannot be treated as a marginal byproduct of the accretion flow and most likely is governed by the Blandford-Znajek mechanism [6], with the black holes possessing very large spins and magnetic fluxes. The magnetic flux required to explain the production of the most powerful jets has recently been found to agree with the maximum magnetic flux that can be confined on black holes by the ram pressure of 'magnetically-arrested-disks' (MAD) [7]. In recent years, the MAD scenario has been thoroughly investigated, and is now considered to be the likely remedy for the production of very powerful jets [8][9][10]. However, as numerical simulations suggest, powers of jets launched in the MAD scenario depend not only on the spin and magnetic flux, but also on the disk's geometrical thickness [11], with the jet power scaling approximately quadratically with all these quantities.
According to standard accretion disk models, accretion disks become geometrically thick only when they are advection dominated, i.e., for the Eddington-ratio λ Edd ≡ L d /L Edd larger than ∼0.3 [12][13][14] or smaller than ∼0.003, where L d is the disk luminosity and L Edd is the Eddington luminosity [15][16][17][18]. Hence, one should not expect to observe powerful jets in AGN with 0.003 < λ Edd < 0.3. This theoretical expectation contradicts with the observational fact that there is no deficit of radio-loud AGN in this range of λ Edd . On the contrary, some studies show a trend of increasing jet production efficiency with decreasing Eddington-ratio [19,20] (see Figure 1). A possible reason for the aforementioned contradiction is modulation of the accretion disk luminosity and jet production. Noting that the power of a jet calculated using the energy content of the radio lobes is actually the time-averaged jet power, averaged over the source lifetime, and that according to the MAD scenario modulation of a jet power at the base is driven by modulation of the accretion rate, modulation of the jet production on time scales shorter than the lifetime of radio lobes will lead to modulation of the "apparent" jet production efficiency and modulation of the Eddington ratio.
For example, modulation of the accretion power by a factor 10 will cause the object to have 10 times lower Eddington ratio and 10 times larger apparent jet production efficiency during its accretion rate minimum relative to that at its maximum, and for a duty cycle 1/2 its apparent jet production efficiency will be about 5 times larger than the real one. A natural driver of variability in the accretion rate and jet production is viscous instabilities in accretion disks [21,22]. Observational support for this hypothesis seems to come from the spatial modulation of the radio brightness distributions seen in some large scale jets [23]. Figure 1. P j /L d ratio as a function of Eddington ratio λ Edd ≡ L d /L Edd . Green diamonds-FRII Narrow-Line-Radio-Galaxies [5]; red squares-Broad-Line-Radio-Galaxies plus radio-loud quasars [20]; grey dots-FRII quasars [24]; and the ellipse marks the approximate location of γ-ray selected FSRQs assumed to have jets with zero pair content [2]. The vertical dashed lines mark the range of the Eddington ratio within which the standard accretion disks are predicted to be radiatively efficient and geometrically thin. (The figure is adopted from Rusinek et al. 2016, in preparation [19]).
However, the situation is complicated if we consider the jet power in γ-ray detected FSRQs (flat-spectrum-radio-quasars) calculated by fitting their broad-band spectra assuming the ERC (external-radiation-Compton) model for γ-ray production [25]. The median of P j /L d calculated by Ghisellini et al. (2014) for jets in FSRQs is ∼10, which is 50 times larger than the median of P j /L d calculated by van Velzen & Falcke (2013) for radio-selected FRII quasars using the radio-lobe calorimetry (Willott et al. 1999 [26]) (see Figure 1). Studies of jet powers in blazars using calorimetry of their extended radio sources [27,28] allowed to verify whether such a difference is associated with selection of two different populations of quasars. As it has been found, the P j /L d median calculated using radio-lobe calorimetry is still much smaller than the P j /L d median calculated using the blazar spectral fits, but now they differ by a factor ∼16 [29]. Whilst the larger averaged values of P j /L d for jets in FSRQs than in FR II quasars, in both cases calculated using the radio-lobe calorimetry, can be explained by selection procedures, the difference between medians of P j /L d measured by the blazar spectral fits and using radio-lobe calorimetry for these same samples must have a different explanation. Following possibilities have been considered: -jet energy losses during propagation between the blazar zone and radio lobes (e.g., by the work done against the external medium by reconfinement shocks which may change their sizes following a jet power modulation by the central engine); -overestimation of jet power using the blazar models (e.g., by assuming zero pair content and/or a one zone model); -a significant fraction of blazars may be hosted by young or short lived compact double radio sources. In this case, the methods of spectral decomposition and core subtraction adopted to use radio-lobe calorimetry in blazars may lead to underestimation of the lobe radio luminosities, and therefore underestimation of the jet power.
In this presentation I will assume that the main reason for the discrepancy between jet power estimates in blazars is overestimation of the jet power by the blazar models by assuming zero pair content. It will be shown that ∼15 pairs per proton is enough to reconcile the difference between the jet power calculated by the two methods and that with such a pair content the magnetic reconnection mechanism is strongly favored as the energy source of blazar activity.

Jet Powers
I will assume that significant jet energy losses (relative to the initial jet power) are taking place only via blazar radiation. In this case the power of a jet leaving the blazar zone, P j , is equal to the rate at which energy is delivered to the radio lobes, P (rl) j , while the power of the jet entering the blazar zone is Then the jet production efficiency is where P rad is the rate of the jet energy losses dominated by the blazar radiation,Ṁ d is the accretion rate, and all powers denoted by P include contribution from both a jet and a counter-jet (assumed to be equal). Now, noting that P j = P B + P p + P e = P (rl) j , and that for FSRQs typically P e P B [2], we have where P B , P p , and P e are the magnetic, proton, and electron energy fluxes, respectively. Since for leptonic radiative models P p = (n p /n e )P p,n p =n e , we can find n e n p = P p,n p =n e P (rl) j where P p,n p =n e is the proton energy flux assuming proton-electron plasma (zero pair content), n e ≡ n e + + n e − , and the number of pairs per proton is n pairs /n p ≡ n e + /n p = (n e /n p − 1)/2. For P p,n p =n e 2.5 × 10 46 ergs · s −1 , P B κ B 10 45 ergs · s −1 , P rad 2.0 × 10 45 ergs · s −1 , andṀ d c 2 = L d / d = 3.16 × 10 46 /( d /0.1) ergs · s −1 which are averages obtained by Ghisellini et al. (2014) [2] for a sample of 191 γ-ray detected FSRQs, and for P (rl) j ∼ 0.1 × P p,n p =n e ∼ 2.5 × 10 45 ergs · s −1 found by Pjanka et al. [29], Equations (1)-(4) give P p ∼ 0.83 × 10 45 ergs · s −1 (→ σ P B /P p ∼ 2) P j,0 ∼ 4.5 × 10 45 ergs · s −1 η j ∼ 0.13( d /0.1) n e /n p ∼ 30 (→ n pairs /n p ∼ 15) where d is the accretion disk radiation efficiency, κ B ≡ (p B + u B )/u B , u B is the magnetic energy density, and p B is the magnetic pressure (not included in Ghisellini et al. [2]). For purely toroidal magnetic field κ B = 2, while for turbulent, isotropic magnetic field κ B = 4/3. We assumed in our estimations κ B = 5/3, which corresponds to equal energy densities of toroidal and turbulent, isotropic magnetic field components.

Jet Magnetization
Jets which are launched by the Blandford-Znajek mechanism are initially strongly Poynting-flux dominated. Following the conversion of the magnetic energy flux to the kinetic energy flux the jets are accelerated [30][31][32][33] reaching bulk Lorentz factors on sub-parsec/parsec scales ranging from a few up to tens. Whether their magnetization parameter σ defined to be the magnetic-to-kinetic energy flux ratio drops to unity or much smaller values is still debated. Whereas blazar models allow a determination of the magnetic energy flux, the kinetic energy flux depends on the unknown proton content (fits of blazar spectra by leptonic models provide us information about a number of electrons plus positrons, not about the number of protons). As was discussed in the preceding sections, pure electron-proton models predict jet powers which, for a large fraction of gamma-ray detected blazars, are not reachable even by the MAD models with maximal jet production efficiency, and have problems accounting for jet powers calculated using the lobe energetics which are more than 10 times smaller than the jet powers calculated using the blazar spectral fits. On the other hand, reducing the jet power by postulating n p n e implies much larger efficiency of radiative processes in blazars than in case of e-p models. Specifically, such efficiency is j ≡ P rad P B + P p + P rad (5) which for P p = P p (n p = n e ) and averaged values of powers taken from Ghisellini et al. [2] gives j 0.07, while for P p P (rl) j − P B gives j 0.45. Such high radiative efficiency and magnetization σ κ B P B /P p 2 favors the reconnection of magnetic fields as the mechanism powering the blazar radiation [34]. Assuming that a change of the jet speed within the blazar zone is negligible, one can find that at the entrance to the blazar zone the jet magnetization is σ 0 ∼ P B,0 /P K (P B + P rad )/P p ∼ 4.4. The ERC model of γ-ray production predicts the ERC-to-synchrotron luminosity peak ratio q = L ERC,peak L syn,peak where ζ is the fraction of disk radiation reprocessed into broad emission lines, θ j is the half-opening angle of a jet, and u ext and u B are energy densities of external radiation field and of magnetic field, respectively [35]. For typical blazar parameters and n p = n e this gives the observed values of the Compton-dominance q 1 only for σ 1 [35,36]. For n p n e the jet production efficiency η j is much lower and therefore q 1 can be reached even for σ > 1.

Particle Acceleration and Spectral Peaks
The shape of the electron/positron energy distribution required to reproduce the observed spectra of high Compton-dominance blazars naturally results from injection of electrons with the break at energy corresponding roughly with the average electron injection energy [35]. Assuming that each electron is accelerated within the blazar zone only once, that energy is given bȳ γ e diss P j,0 m e c 2Ṅ e Γ whereṄ e = (n e /n p )Ṅ p = (n e /n p )Ṁ p /m p ,Ṁ p is the proton mass flux in a jet, and e is the fraction of dissipated energy channeled to relativistic electrons and positrons. For the jet power dominated by magnetic and proton energy fluxes, and then in case of Γ = Γ 0 , For Γ 0 1 these equations giveγ Electrons with such energy are expected to produce synchrotron and ERC luminosity peaks. In order to reproduce their location in the electromagnetic spectrum, with the synchrotron emission peaking at around 10 13 Hz and the ERC emission peaking at around 10 22 Hz [37], the electron energy should be of the order of 100 [38]. For the parameters derived above (σ 0 4.4, σ 2, and n e /n p 30), and noting that the fraction of dissipated energy channeled to electrons/positrons is expected to be 0.5 < e < 1, this condition is satisfied.

Discussion and Conclusions
The jet powers in FR II quasars calculated using radio-lobe calorimetry [24] are a factor ∼50 smaller than jet powers calculated using spectral fits of γ-ray detected FSRQs when assuming n p = n e jet plasma [2] (see Figure 1). Measurements of radio luminosities of extended radio-structures in some FSRQs [27,28] has enabled a comparison of the results of the two jet power measurement techniques by applying both techniques to the same objects. As Pjanka et al. [29] showed, the average difference for the cross-matched samples is smaller, but still significant, by a factor ∼16. Noting that the power calculated by Ghisellini et al. includes the blazar radiation, whilst the power calculated using the radio lobe energetics does not, the difference is reduced to a factor ∼8. Assuming that this remaining difference is due to overestimation the blazar power by Ghisellini et al. by ignoring pair plasma in a jet, we found that on average: • the pair content: n pairs /n p ∼ 15; • the jet production efficiency: η j 0.14 ( d /0.1); • the magnetization at entrance and exit of blazar zone is σ 0 4.4 and σ 2, respectively; • the average energies of accelerated electrons/positrons correspond to typical observed locations of synchrotron and γ-ray luminosity peaks provided the acceleration is powered by magnetic reconnection operating in a jet at a distance corresponding to spatial extension of the BLR.
The pair content obtained above is larger by a factor 3/2 than the maximal one, (n pairs /n p ) max ∼ 10, above which the jet is expected to be efficiently decelerated due to the Compton-rocket effect [39]. However this limit was obtained by assuming an isotropic distribution of the external diffuse radiation field and neglecting Klein-Nishina effects and therefore can be relaxed somewhat if some level of flattening of the BLR is considered, and taking into account Klein-Nishina recoil effects [40]. Noting that the MAD scenario predicts production of jets with powers P j,0 a 2 (H/R) 2Ṁ d c 2 (11) and assuming that spins a in γ-ray selected quasars are close to their maximum value, ∼1, the value of η j ∼ 0.14 obtained above may simply result from suppression of jet production efficiency due to having geometrically thin disks in quasars, i.e., with H/R 1. This requires H/R ∼ 1/3, which is much larger than predicted by standard accretion disk models [41,42]. However, as some studies indicate, real disks can be thicker than the standard disks [43][44][45]. Then the ∼7 times lower jet production efficiency in radio selected FRII quasars than in γ-ray detected FSRQs may result from on average ∼ √ 7 times lower BH spins in the entire FRII quasar population than in the γ-ray selected one. And this can be explained by noting, that in γ-ray loud quasars jets are more powerful and therefore are expected to be more relativistic. This interpretation is supported by the fact that only about 9% of the FSRQs are found to be γ-ray loud [46].