Abstract

The analysis of the frequency dependence of the observed shift of the cores of relativistic jets in active galactic nuclei (AGNs) allows us to evaluate the number density of the outflowing plasma ne and, hence, the multiplicity parameter λ = ne/nGJ, where nGJ is the Goldreich–Julian number density. We have obtained the median value for λmed = 3 × 1013 and the median value for the Michel magnetization parameter σM, med = 8 from an analysis of 97 sources. Since the magnetization parameter can be interpreted as the maximum possible Lorentz factor Γ of the bulk motion which can be obtained for relativistic magnetohydrodynamic (MHD) flow, this estimate is in agreement with the observed superluminal motion of bright features in AGN jets. Moreover, knowing these key parameters, one can determine the transverse structure of the flow. We show that the poloidal magnetic field and particle number density are much larger in the centre of the jet than near the jet boundary. The MHD model can also explain the typical observed level of jet acceleration. Finally, casual connectivity of strongly collimated jets is discussed.

1 INTRODUCTION

Strongly collimated jets represent one of the most visible signs of the activity of compact astrophysical sources. They are observed both in relativistic objects such as active galactic nuclei (AGNs) and microquasars, and in young stars where the motion of matter is definitely non-relativistic. This implies that we are dealing with some universal and extremely efficient mechanism of energy release.

At present the magnetohydrodynamic (MHD) model of activity of compact objects is accepted by most astrophysicists (Mestel 1999; Krolik 1999). At the heart of the MHD approach lies the model of the unipolar inductor, i.e. a rotating source of direct current. It is believed that the electromagnetic energy flux – the Poynting flux – plays the main role in the energy transfer from the ‘central engine’ to active regions. The conditions for the existence of such a ‘central engine’ are satisfied in all the compact sources mentioned above. Indeed, all compact sources are assumed to harbour a rapidly spinning central body (black hole, neutron star or young star) and some regular magnetic field, which leads to the emergence of strong induction electric fields. The electric fields, in turn, lead to the appearance of longitudinal electric currents resulting in effective energy losses and particle acceleration.

The first studies of the electromagnetic model of compact sources (namely radio pulsars) were carried out as early as the end of the 1960s (Goldreich & Julian 1969; Michel 1969). It was evidenced that there are objects in the Universe in which electrodynamical processes can play the decisive role in the energy release. Then, Blandford (1976) and Lovelace (1976) independently suggested that the same mechanism can also operate in AGNs, and for nearly 40 years this model has remained the leading one.

Remember that within the MHD approach the total energy losses Pjet can be easily evaluated as Pjet ∼ IU, where I is the total electric current flowing along the jet, and U ∼ ER0 is the electro-motive force exerted on the black hole on the scale R0. If the central engine (black hole for AGNs) rotates with the angular velocity Ω in the external magnetic field B0, one can evaluate the electric field as E ∼ (ΩR0/c)B0. On the other hand, assuming that for AGNs the current density j is fully determined by the relativistic outflow of the Goldreich–Julian (GJ) charge density ρGJ = ΩB0/(2πc) (i.e. the minimum charge density required for the screening of the longitudinal electric field in the magnetosphere), one can write down j ≈ cρGJ. It gives
(1)
As for AGNs, one can set R0 ≈ rg = 2GM/c2 to obtain the well-known evaluation (Blandford & Znajek 1977)
(2)
In particular, comparing expressions (1) and (2), one can straightforwardly obtain
(3)
For AGNs it corresponds to 1019–1021 A. Certainly, the question as to whether it is possible to consider a black hole immersed in an external magnetic field as a unipolar inductor turned out to be also rather nontrivial (Punsly 2001; Beskin 2009; Okamoto 2009).

As a result, the MHD model was successfully used to describe a lot of processes in active nuclei including the problem of the stability of jets (Benford 1981; Hardee & Norman 1988; Appl & Camenzind 1992; Istomin & Pariev 1994; Bisnovatyi-Kogan 2007; Lyubarsky 2009) and their synchrotron radiation (Blandford & Königl 1979; Pariev, Istomin & Beresnyak 2003; Lyutikov, Pariev & Gabuzda 2005). In particular, it was shown both analytically (Bogovalov 1995; Heyvaerts & Norman 2003; Beskin & Nokhrina 2009) and numerically (Komissarov et al. 2006; Tchekhovskoy, McKinney & Narayan 2009; Porth et al. 2011; McKinney, Tchekhovskoy & Blanford 2012) that for sufficiently small ambient pressure the dense core can be formed. This is related both to advances in the theory which have at last formulated sufficiently simple analytical relations (Blandford & Znajek 1977; Beskin 2010), and to the breakthrough in numerical simulations (Komissarov et al. 2006; Tchekhovskoy et al. 2009; Porth et al. 2011; McKinney et al. 2012) which confirmed theoretical predictions.

Moreover, recently Kronberg et al. (2011) demonstrated that in QSO 3C 303 the jet does possess large enough toroidal magnetic field, the appropriate longitudinal electric current along the jet I ≈ 1.7 × 1019 A being as large as the electric current IGJ (equation 1) which is necessary to support the Poynting energy flux. Besides, the lack of γ-ray radiation as probed by the Fermi Observatory for AGN jets observed at small enough viewing angles ϑ < 5° (Savolainen et al. 2010) can be easily explained as well. Indeed, as was found (see e.g. Beskin 2010), within well-collimated magnetically dominated MHD jets the Lorentz factors of the particle bulk motion can be evaluated as
(4)
where r is the distance from the jet axis, and RL = c/Ω is the light cylinder radius. Thus, the energy of particles radiating in small enough angles ϑ with respect to the jet axis is to be much smaller than that corresponding to peripheral parts of a jet.
The most important MHD parameters describing relativistic flows (which was originally introduced for radio pulsars) are the Michel magnetization parameter σM and the multiplicity parameter λ. The first one determines the maximum possible bulk Lorentz factor Γ of the flow when all the energy transported by the Poynting flux is transmitted to particles. The second one is the dimensionless multiplicity parameter λ = ne/nGJ, which is defined as the ratio of the number density ne to the GJ number density nGJ = ΩB/2πce. It is important that these two parameters are connected by the simple relation (Beskin 2010)
(5)
Here, PA = m2c5/e2 ≈ 1017 erg s−1 is the minimum energy losses of the central engine which can accelerate particles to relativistic energies, and Pjet is the total energy losses of the compact object.

Unfortunately, up to now neither the magnetization nor the multiplicity parameters were actually known as the observations could not give us the direct information about the number density and bulk energy of particles. The core-shift method has been applied to obtain the concentration ne, magnetic field B (Lobanov 1998; O'Sullivan & Gabuzda 2009; Pushkarev et al. 2012; Zdziarski et al. 2014) and the jet composition (Hirotani 2005) in AGN jets. However, evaluation of multiplicity and Michel magnetization parameters, which needs to estimate the total jet power, has not been done. From a theoretical point of view if the inner parts of the accretion disc are hot enough, then electron–positron pairs can be produced by two-photon collisions, where photons with sufficient energy originate from the inner parts of the accretion disc (Blandford & Znajek 1977; Moscibrodzka et al. 2011). In this case λ ∼ 1010–1013, and the Michel magnetization parameter σM ∼ 10–103. The second model takes into account the appearance of the region where the GJ plasma density is equal to zero due to general relativity effects that corresponds to the outer gap in the pulsar magnetosphere (Beskin, Istomin & Pariev 1992; Hirotani & Okamoto 1998). This model gives λ ∼ 102–103, and σM ∼ 1010–1013.

This large difference in the estimates for the magnetization parameter σM leads to two completely different pictures of the flow structure in jets. In particular, it determines whether the flow is magnetically or particle dominated. The point is that for ordinary jets r/RL ∼ 104–105. As a result, using the universal asymptotic solution Γ ≈ r/RL (4), one can obtain that the values σM ∼ 10–103 correspond to the saturation regime when almost all the Poynting flux Pem is transmitted to the particle kinetic energy flux Ppart. On the other hand, for σM ∼ 1012 the jet remains magnetically dominated (Ppart ≪ Pem). Thus, the determination of the Michel magnetization parameter σM is the key point in the analysis of the internal structure of relativistic jets.

The paper is organized as follows. In Section 2, it is shown that VLBI (very long baseline interferometry) observations of synchrotron self-absorption in AGN jets allow us to evaluate the number density of the outflowing plasma ne and, hence, the multiplicity parameter λ. We discuss the source sample and present the result for multiplicity and Michel magnetization parameters in Section 3. The values λ ∼ 1013 obtained from the analysis of 97 sources shows that for most jets the magnetization parameter σM ≲ 30. Since the magnetization parameter is the maximum possible value of the Lorentz factor of the relativistic bulk flow, this estimate is consistent with observed superluminal motion. In Section 4, it is shown that for physical parameters determined above, the poloidal magnetic field and particle number density are much larger in the centre of the jet than near its boundary. Finally, in Section 5 the casual connectivity of strongly collimated supersonic jets is discussed. Throughout the paper, we use the Λ cold dark matter cosmological model with H0 = 71 km s−1 Mpc−1, Ωm = 0.27 and |$\Omega _\Lambda =0.73$| (Komatsu et al. 2009).

2 THE METHOD

2.1 General relations

To determine the multiplicity parameter λ and the Michel magnetization parameter σM one can use the dependence on the visible position of the core of the jet from the observation frequency (Blandford & Königl 1979; Gould 1979; Marscher 1983; Lobanov 1998; Hirotani 2005; Kovalev et al. 2008; O'Sullivan & Gabuzda 2009; Sokolovsky et al. 2011; Pushkarev et al. 2012). This effect is associated with the absorption of the synchrotron photon gas by relativistic electrons (positrons) in a jet.

Typically, the parsec-scale radio morphology of a bright AGN manifests a one-sided jet structure due to Doppler boosting that enhances the emission of the approaching jet. The apparent base of the jet is commonly called the ‘core’, and it is often the brightest and most compact feature in VLBI images of AGN. The VLBI core is thought to represent the jet region where the optical depth is equal to unity.

We will employ the following model to connect the physical parameters at the jet launching region with the observable core-shift. There is a MHD relativistic outflow of non-emitting plasma moving with bulk Lorentz factor Γ and number density ne in the observer rest frame. On the latter we superimpose the flow of emitting particles with distribution |${\rm d}n_{{\rm syn}*} = k_{{\rm e}*}\gamma _{*}^{-1+2\alpha }{\rm d}\gamma _{*}$|⁠, γ* ∈ [γmin*; γmax*]. Here, nsyn* is number density of emitting plasma, ke* is number density amplitude and γ* is the emitting particles’ Lorentz factor. All the parameters with subscript ‘*’ are taken in the non-emitting plasma rest frame, i.e. in the frame which locally moves with the bulk Lorentz factor Γ.

We suppose that the emitting particles radiate synchrotron photons in the jet's magnetic field, and these photons scatter off the same electrons, which lead to the photon absorption (Gould 1979; Lobanov 1998; Hirotani 2005). The corresponding turn-over frequency νm*, the frequency at which the flux density Sν has a maximum, can be evaluated using expressions from Gould (1979) as
(6)
The function cα ( · ) is a composition of gamma-functions defined by Gould (1979), and for α = −1/2 we have cα (2) = 1.396. Constants e, m and c are the electron charge, electron mass and the speed of light correspondingly. Finally, B* is the magnitude of disordered magnetic field in an emitting region with a characteristic size R* along the line of sight.
Although we assume that the toroidal magnetic field dominates in the jet, an assumption of disordered magnetic field in our opinion can be retained, because, for an optically thin jet, the photon meets both directions of field. Thus, the mean magnetic field along the photon path is almost zero, which mimics the behaviour of a disordered field. As a result, the parameters in the observer rest frame and plasma rest frame are connected by the following equations:
(7)
(8)
(9)
(10)
where z is the redshift,
(11)
is the Doppler factor, χ is the jet half-opening angle and φ is a viewing angle.
Further, the number density of emitting electrons nsyn is connected with the amplitude ke as
(12)
where γ = γ*Γ. For α = −1/2, we get
(13)
We also put nsyn = ξne. Here ξ is a ratio of the number density of emitting particles to the MHD flow number density. The portion of particles effectively accelerated by the internal shocks was found by Sironi, Spitkovsky & Arons (2013) to be about 1 per cent, so we take ξ ≈ 0.01.
Finally, we assume (Lobanov 1998; Hirotani 2005) the following power-law behaviour for the magnetic field and particle density dependence on distance:
(14)
(15)
where B1 is the magnetic field and n1 is the number density at r1 = 1 pc, respectively. For these scalings of particle density and magnetic field with the distance the turn-over frequency νm as a function of r does not depend on α and can be written as
(16)
This scaling has been confirmed by Sokolovsky et al. (2011) in measurements of core shifts for 20 AGNs made for nine frequencies each. Using these dependences of magnetic field and particle number density of distance r, we obtain in the observer rest frame
(17)
where |$C=c_{\alpha }^2(1-2\alpha )/5(5-2\alpha )$|⁠.
On the other hand, the values B1 and n1 can be related through introducing the flow magnetization parameter σ – the ratio of Poynting vector to particle kinetic energy flux at a given distance along the flow (see Appendix A). Let us define the magnetization σξ as a ratio of Poynting vector to the total kinetic energy flux of emitting and non-emitting particles:
(18)
Here, the kinetic energy flux of emitting electrons is
(19)
and function F ( · ) for α = −1/2 is defined by the following expression:
(20)
Estimating now the Poynting vector as
(21)
and particle kinetic energy fluxes as
(22)
we obtain the following relationship between magnetic field and particle number density:
(23)
In what follows, we neglect the term ξFγmin in comparison with Γ. Further on we omit the index ξ. Using equation (23), we get
(24)
As to the number density n1, it can be defined through the multiplicity parameter λ and total jet energy losses Pjet as (see Appendix B for more detail)
(25)

2.2 The saturation regime

To determine the intrinsic parameters of relativistic jets, let us consider two cases for the different magnetization at 1 pc. In what follows, we assume that the flow at its base is highly, or at least mildly, magnetized, i.e. σM > 1.

First, we assume that up to the distance r = 1 pc the plasma has been effectively accelerated so that the Poynting flux is smaller in comparison with the particle kinetic energy flux, i.e. σ1 ≲ 1. In other words, the acceleration reached the saturation regime (Beskin 2010). Combining now equations (4) and (A8), it is easy to obtain that this case corresponds to σM < 105. Accordingly, the bulk Lorentz factor at r = 1 pc can be evaluated as
(26)
In this case equation (24) can be rewritten as
(27)
Using now the relationship between the angular distance θd and the distance from the jet base r
(28)
where DL is the luminosity distance, we obtain
(29)
This expression can be rewritten as a following relationship between the core position and the observation frequency:
(30)
Having the measured core-shift Δrmas in milliarcseconds for two frequencies νm, 1 and νm, 2, we obtain for α = −1/2:
(31)
Accordingly, using equation (5), we obtain
(32)
As we see, this value is in agreement with our assumption σM < 105.

2.3 Highly magnetized outflow

Let us now assume that the flow is still highly magnetized at a distance of the observable core. This implies that the Michel magnetization parameter σM > 105. Using now relation (A8), one can obtain
(33)
On the other hand, equation (24) can be rewritten as
(34)
This gives the following expression for the Michel magnetization parameter
(35)
As we see, these values are in contradiction with our assumption σM > 105. Thus, one can conclude that it is the saturation limit that corresponds to parsec-scale relativistic jets under consideration.

3 THE STATISTICS FOR MULTIPLICITY PARAMETER

Several methods can be applied to measure the apparent shift of the core position as discussed by Kovalev et al. (2008). As a result, a magnitude of the shift, designated by η, can be measured and presented in units [mas GHz] or [pc GHz]. Knowing this quantity, one can use the expressions (31)–(35) to estimate the multiplicity and magnetization parameters.

3.1 The sample of objects

In our analysis, we use the results of two surveys of the apparent core shift in AGN jets: Sokolovsky et al. (2011) show results for 20 objects obtained from nine frequencies between 1.4 and 15.3 GHz (S-sample) and Pushkarev et al. (2012) have results for 163 AGN from four frequencies covering 8.1–15.3 GHz (P-sample). Of these we use only those sources for which the apparent opening angle is known from Pushkarev et al. (2009). As a result, 97 sources are left from the P-sample and 5 from the S-sample. Although all of S-sample sources are in P-sample, we have included them as an independent measurement of core shift. Moreover, for objects 0215+015 and 1219+285, the two measurements of core shift have been made for two different epochs, and we included them too. This leaves us with 97 sources and 104 measurements of core shift.

The distance to the objects is determined from the redshift and accepted cosmology model. For a Doppler factor we use the estimate δ ≈ βapp, where measured apparent velocity βapp is a ratio of apparent speed of a bright feature in a jet to the speed of light. We believe this to be a good estimate because Cohen et al. (2007) have showed using Monte Carlo simulations that the probability density p(δ| βapp) to observe a Doppler factor for a given apparent velocity is peaked around unity. This is done under an assumption that the measured βapp does represent the underlying jet flow. The redshifts z and the apparent velocities βapp are taken from Lister et al. (2013).

The value of observation angle φ we obtain from the set of equations for Doppler factor δ and apparent velocity
(36)
Taking δ = βapp, we obtain from equations (11) and (36) for the observation angle the relation
(37)
The half-opening angle related to the observed opening angle χapp as
(38)
We use the values for χapp derived by Pushkarev et al. (2009) with typical errors of 1| $_{.}^{\circ}$|5. We also have chosen parameter γmin = 1.
We evaluate the total jet power Pjet through the relationship (Cavagnolo et al. 2010) between the luminosities of jets in radio band and mechanical jet power, needed to form the cavities in surrounding gas. The power law, found by Cavagnolo et al. (2010) for a range of frequencies 200–400 MHz is
(39)
In order to find flux density measurements at the 92 cm band for each source we use the CATS data base (Verkhodanov et al. 1997) which accumulates measurements at different epochs and from the different catalogues. The data which we use in this paper were originally reported by De Breuck et al. (2002), Douglas et al. (1996), Ghosh, Gopal-Krishna & Rao (1994), Kühr et al. (1979), Kühr et al. (1981), Gregory & Condon (1991), Mitchell et al. (1994) and Rengelink et al. (1997) with a typical flux density accuracy of about 10 per cent.

The typical error for core-shift measurements in Pushkarev et al. (2012) and Sokolovsky et al. (2011) is 0.05 mas. There are 23 objects in our sample that have the core-shift values less than 0.05 mas. For them we have replaced the core-shift values by 0.05 mas for our calculations for convenience of the λ and σM analysis.

3.2 Results and discussion

Using the formula (31), we obtain the following result for the equipartition regime. The obtained values for the multiplication parameter λ and magnetization parameter σM are presented in Table 1. Their distributions are shown in Fig. 1 and Fig. 2, respectively. In cases when more than one estimate is determined per source (e.g. for 0215+015), an average value is used in the histograms. The resultant median value for the multiplicity parameter λmed = 3 × 1013, and median value for magnetization parameter σM, med = 8. The multiplicity parameter for our sample lies in the interval (3 × 1012; 4 × 1014), and the Michel magnetization parameter σM lies correspondingly in the (0.4; 61) interval.

Distributions of the multiplicity parameter λ for the sample of 97 sources. Two objects with λ = 2.8 × 1014 and 3.6 × 1014 lie out of the shown range of values.
Figure 1.

Distributions of the multiplicity parameter λ for the sample of 97 sources. Two objects with λ = 2.8 × 1014 and 3.6 × 1014 lie out of the shown range of values.

Distributions of the Michel magnetization parameter σM for the sample of 97 sources.
Figure 2.

Distributions of the Michel magnetization parameter σM for the sample of 97 sources.

Table 1.

Jet parameters and derived multiplication and magnetization parameters.

SourcezβappχappS0.3PjetReferenceΔrcoreEpochλσ
(c)(°)(Jy)(1045 erg s−1)for S0.3(mas)for Δrcore(1013)
(1)(2)(3)(4)(5)(6)(7)(8)(9)(10)(11)
0003−0660.3478.4016.32.171.0720.0352006-07-071.219.69
0106+0132.09924.3723.62.8510.5060.0052006-07-072.0223.67
0119+1150.57018.5715.62.241.8660.3472006-06-153.844.27
0133+4760.85915.3621.71.632.5480.1312006-08-093.525.80
0202+1490.40515.8916.46.252.3930.1222006-09-061.6310.90
0202+3191.46610.1513.40.762.9980.0132006-08-092.1711.13
0212+7352.3676.5516.41.545.1780.1492006-07-079.824.41
0215+0151.71525.0636.70.883.9050.0882006-04-283.757.54
0.2412006-12-017.973.54
0234+2851.20621.9919.81.453.5250.2752006-09-065.294.79
0333+3211.25913.078.03.686.7270.2792006-07-074.108.62
0336−0190.85224.4526.82.493.2660.1172006-08-092.668.68
0403−1320.57120.8016.47.624.4520.3462006-05-243.666.95
0420−0140.9165.7422.70.871.8430.2672006-10-0613.491.30
0458−0202.28613.5723.13.5411.2030.0062006-11-103.2017.26
0528+1342.07017.3416.11.025.8520.1672006-10-064.817.40
0529+0751.25418.0356.41.754.5420.1102006-08-093.837.58
0605−0850.87019.1914.01.432.3930.0962006-11-101.6611.95
0607−1570.3231.91835.12.310.9630.2402006-09-0628.460.38
0642+4493.3968.5323.40.707.6880.1102006-10-065.278.29
0730+5040.72014.0714.80.711.2080.2622006-05-244.273.21
0735+1780.4505.0421.01.811.2380.0392006-04-282.565.07
0736+0170.18913.7917.91.790.4230.0792006-06-150.828.38
0738+3130.63110.7210.51.261.4830.1832006-09-062.855.23
0748+1260.88914.5816.21.452.6520.0982006-08-092.418.69
0754+1000.26614.4013.70.740.3920.2662006-04-282.063.32
0804+4991.4361.1535.30.602.4980.0942006-10-0635.880.61
0805−0771.83741.7618.82.609.2620.2072006-05-243.1214.09
0823+0330.50512.8813.40.630.7150.1412006-06-152.124.71
0827+2430.94219.8114.60.711.8020.1502006-05-242.526.92
0829+0460.17410.1318.70.670.2120.1092006-07-071.283.82
0836+7102.21821.0812.45.0617.7520.1862006-09-063.8116.43
0851+2020.30615.1428.51.110.5630.0282006-04-281.087.69
0906+0151.02622.0817.51.543.0530.1682006-10-063.047.56
0917+6241.45312.0715.91.274.0780.1122006-08-093.957.14
0923+3920.6952.7610.83.283.0350.0422006-07-073.236.69
0945+4081.24920.2014.02.946.3020.0832006-06-151.8018.92
1036+0540.4735.726.50.750.8020.1952006-05-242.773.80
1038+0641.26510.696.71.594.3220.1062006-10-062.0214.01
1045−1880.59510.518.02.792.4620.1562006-09-062.019.45
1127−1451.18414.8916.15.638.9460.0962006-08-092.7314.77
1150+8121.25010.1115.01.393.8180.0872006-06-153.318.03
1156+2950.72524.5916.74.333.8980.1622006-09-062.1511.40
1219+0440.9660.8213.01.142.5240.1332006-05-2422.110.94
1219+2850.1039.1213.91.770.1930.1822006-02-121.114.04
0.1422007-04-300.934.87
0.1992006-11-101.193.78
1222+2160.43426.6010.83.981.9050.1802006-04-281.1414.13
1226+0230.15814.861063.723.2530.0202006-03-090.3160.77
1253−0550.53620.5814.416.566.3130.0482006-04-050.7539.84
1308+3260.99727.4818.51.422.7980.1432006-07-072.359.33
1334−1270.53916.3312.61.911.7120.2372006-10-062.615.98
1413+1350.2471.788.82.740.8120.2302006-08-096.021.64
1458+7180.9046.614.519.6413.3080.0812006-09-061.4632.22
0.1362007-03-012.1621.84
1502+1061.83917.5337.91.084.9230.0522006-07-073.598.92
1504−1660.8763.9418.41.802.7930.1482006-12-019.562.05
1510−0890.36028.0015.22.751.2230.1222006-04-280.9313.47
1514−2410.0496.397.82.060.0830.1882006-04-280.565.15
1538+1490.6068.7416.12.822.3630.0322006-06-151.6811.09
1546+0270.41412.0812.90.700.6130.0102006-08-090.8710.32
1606+1061.23219.0924.02.675.3070.0572006-07-072.1114.84
1611+3431.40029.1526.94.208.4430.0572006-06-151.7922.54
1633+3821.81329.2222.62.518.2880.1192006-09-063.0713.50
1637+5740.75113.5910.71.321.8880.1172006-05-241.928.94
1638+3981.66615.8553.80.643.1180.0072006-08-094.685.36
1641+3990.59319.2712.99.935.1380.2112006-06-152.2911.99
1655+0770.62114.775.52.362.3360.0802006-11-100.7325.39
0.0862007-06-010.7724.05
1726+4550.7172.3016.50.490.9580.0092006-09-065.182.34
1730−1300.90227.3510.46.466.5430.1742006-07-071.6719.72
1749+0960.3227.9016.81.200.6160.0612006-06-151.436.15
1751+2881.1183.8712.10.401.5520.0072006-10-063.604.62
1758+3882.0922.2117.90.181.8280.0792006-11-1013.981.42
1803+7840.68010.7918.41.922.2380.0292006-09-061.7110.80
0.0612007-05-031.989.31
1823+5680.66426.176.82.632.5280.0522006-07-070.4246.21
1828+4870.69213.077.147.7815.6030.1172006-08-091.3935.35
1849+6700.65723.0816.60.861.2280.0242006-05-240.8815.50
1908−2011.1194.3923.92.705.2120.2462006-03-0918.031.69
1928+7380.3028.179.84.811.4080.1472006-04-281.727.66
1936−1551.6575.3435.20.673.4520.2152006-07-0722.921.15
2008−1591.1804.859.70.732.4120.0082006-11-102.657.89
2022−0771.38823.2319.62.636.6720.0062006-04-051.5123.67
2121+0531.94111.6634.00.633.9920.1522006-06-1510.292.83
2128−1230.5015.995.01.471.2330.2232006-10-062.525.20
2131−0211.28419.9618.42.666.1160.0892006-08-092.3914.14
2134+0041.9325.0415.20.994.8560.1882006-07-0712.352.60
2136+1412.4274.1532.50.946.1660.0082006-09-0610.283.64
2145+0670.9992.8323.23.765.1830.0082006-10-066.974.31
2155−1520.67218.1217.62.412.4330.4052006-12-015.343.60
2200+4200.0699.9526.21.820.1280.0322006-04-050.477.30
2201+1711.07617.6613.61.002.6320.3802006-05-245.643.82
2201+3150.2958.2712.81.820.8830.3472006-10-063.902.67
0.1922007-04-302.503.80
2209+2361.1252.2914.20.391.5120.0382006-12-016.032.73
2216−0380.9016.7315.62.253.5760.0112006-08-092.559.57
2223−0521.40420.3411.713.5918.0030.1992006-10-063.2118.33
2227−0881.5602.0015.81.415.1420.1862006-07-0722.851.40
2230+1141.0378.6213.38.519.2530.2782006-02-127.365.45
2243−1230.6325.2414.81.451.7110.1612006-09-065.732.79
2251+1580.85913.7740.912.479.3930.1242006-06-158.314.72
2345−1670.57611.4715.82.812.2130.1672006-11-103.245.54
2351+4561.98621.5620.12.238.5480.1962006-05-245.357.99
SourcezβappχappS0.3PjetReferenceΔrcoreEpochλσ
(c)(°)(Jy)(1045 erg s−1)for S0.3(mas)for Δrcore(1013)
(1)(2)(3)(4)(5)(6)(7)(8)(9)(10)(11)
0003−0660.3478.4016.32.171.0720.0352006-07-071.219.69
0106+0132.09924.3723.62.8510.5060.0052006-07-072.0223.67
0119+1150.57018.5715.62.241.8660.3472006-06-153.844.27
0133+4760.85915.3621.71.632.5480.1312006-08-093.525.80
0202+1490.40515.8916.46.252.3930.1222006-09-061.6310.90
0202+3191.46610.1513.40.762.9980.0132006-08-092.1711.13
0212+7352.3676.5516.41.545.1780.1492006-07-079.824.41
0215+0151.71525.0636.70.883.9050.0882006-04-283.757.54
0.2412006-12-017.973.54
0234+2851.20621.9919.81.453.5250.2752006-09-065.294.79
0333+3211.25913.078.03.686.7270.2792006-07-074.108.62
0336−0190.85224.4526.82.493.2660.1172006-08-092.668.68
0403−1320.57120.8016.47.624.4520.3462006-05-243.666.95
0420−0140.9165.7422.70.871.8430.2672006-10-0613.491.30
0458−0202.28613.5723.13.5411.2030.0062006-11-103.2017.26
0528+1342.07017.3416.11.025.8520.1672006-10-064.817.40
0529+0751.25418.0356.41.754.5420.1102006-08-093.837.58
0605−0850.87019.1914.01.432.3930.0962006-11-101.6611.95
0607−1570.3231.91835.12.310.9630.2402006-09-0628.460.38
0642+4493.3968.5323.40.707.6880.1102006-10-065.278.29
0730+5040.72014.0714.80.711.2080.2622006-05-244.273.21
0735+1780.4505.0421.01.811.2380.0392006-04-282.565.07
0736+0170.18913.7917.91.790.4230.0792006-06-150.828.38
0738+3130.63110.7210.51.261.4830.1832006-09-062.855.23
0748+1260.88914.5816.21.452.6520.0982006-08-092.418.69
0754+1000.26614.4013.70.740.3920.2662006-04-282.063.32
0804+4991.4361.1535.30.602.4980.0942006-10-0635.880.61
0805−0771.83741.7618.82.609.2620.2072006-05-243.1214.09
0823+0330.50512.8813.40.630.7150.1412006-06-152.124.71
0827+2430.94219.8114.60.711.8020.1502006-05-242.526.92
0829+0460.17410.1318.70.670.2120.1092006-07-071.283.82
0836+7102.21821.0812.45.0617.7520.1862006-09-063.8116.43
0851+2020.30615.1428.51.110.5630.0282006-04-281.087.69
0906+0151.02622.0817.51.543.0530.1682006-10-063.047.56
0917+6241.45312.0715.91.274.0780.1122006-08-093.957.14
0923+3920.6952.7610.83.283.0350.0422006-07-073.236.69
0945+4081.24920.2014.02.946.3020.0832006-06-151.8018.92
1036+0540.4735.726.50.750.8020.1952006-05-242.773.80
1038+0641.26510.696.71.594.3220.1062006-10-062.0214.01
1045−1880.59510.518.02.792.4620.1562006-09-062.019.45
1127−1451.18414.8916.15.638.9460.0962006-08-092.7314.77
1150+8121.25010.1115.01.393.8180.0872006-06-153.318.03
1156+2950.72524.5916.74.333.8980.1622006-09-062.1511.40
1219+0440.9660.8213.01.142.5240.1332006-05-2422.110.94
1219+2850.1039.1213.91.770.1930.1822006-02-121.114.04
0.1422007-04-300.934.87
0.1992006-11-101.193.78
1222+2160.43426.6010.83.981.9050.1802006-04-281.1414.13
1226+0230.15814.861063.723.2530.0202006-03-090.3160.77
1253−0550.53620.5814.416.566.3130.0482006-04-050.7539.84
1308+3260.99727.4818.51.422.7980.1432006-07-072.359.33
1334−1270.53916.3312.61.911.7120.2372006-10-062.615.98
1413+1350.2471.788.82.740.8120.2302006-08-096.021.64
1458+7180.9046.614.519.6413.3080.0812006-09-061.4632.22
0.1362007-03-012.1621.84
1502+1061.83917.5337.91.084.9230.0522006-07-073.598.92
1504−1660.8763.9418.41.802.7930.1482006-12-019.562.05
1510−0890.36028.0015.22.751.2230.1222006-04-280.9313.47
1514−2410.0496.397.82.060.0830.1882006-04-280.565.15
1538+1490.6068.7416.12.822.3630.0322006-06-151.6811.09
1546+0270.41412.0812.90.700.6130.0102006-08-090.8710.32
1606+1061.23219.0924.02.675.3070.0572006-07-072.1114.84
1611+3431.40029.1526.94.208.4430.0572006-06-151.7922.54
1633+3821.81329.2222.62.518.2880.1192006-09-063.0713.50
1637+5740.75113.5910.71.321.8880.1172006-05-241.928.94
1638+3981.66615.8553.80.643.1180.0072006-08-094.685.36
1641+3990.59319.2712.99.935.1380.2112006-06-152.2911.99
1655+0770.62114.775.52.362.3360.0802006-11-100.7325.39
0.0862007-06-010.7724.05
1726+4550.7172.3016.50.490.9580.0092006-09-065.182.34
1730−1300.90227.3510.46.466.5430.1742006-07-071.6719.72
1749+0960.3227.9016.81.200.6160.0612006-06-151.436.15
1751+2881.1183.8712.10.401.5520.0072006-10-063.604.62
1758+3882.0922.2117.90.181.8280.0792006-11-1013.981.42
1803+7840.68010.7918.41.922.2380.0292006-09-061.7110.80
0.0612007-05-031.989.31
1823+5680.66426.176.82.632.5280.0522006-07-070.4246.21
1828+4870.69213.077.147.7815.6030.1172006-08-091.3935.35
1849+6700.65723.0816.60.861.2280.0242006-05-240.8815.50
1908−2011.1194.3923.92.705.2120.2462006-03-0918.031.69
1928+7380.3028.179.84.811.4080.1472006-04-281.727.66
1936−1551.6575.3435.20.673.4520.2152006-07-0722.921.15
2008−1591.1804.859.70.732.4120.0082006-11-102.657.89
2022−0771.38823.2319.62.636.6720.0062006-04-051.5123.67
2121+0531.94111.6634.00.633.9920.1522006-06-1510.292.83
2128−1230.5015.995.01.471.2330.2232006-10-062.525.20
2131−0211.28419.9618.42.666.1160.0892006-08-092.3914.14
2134+0041.9325.0415.20.994.8560.1882006-07-0712.352.60
2136+1412.4274.1532.50.946.1660.0082006-09-0610.283.64
2145+0670.9992.8323.23.765.1830.0082006-10-066.974.31
2155−1520.67218.1217.62.412.4330.4052006-12-015.343.60
2200+4200.0699.9526.21.820.1280.0322006-04-050.477.30
2201+1711.07617.6613.61.002.6320.3802006-05-245.643.82
2201+3150.2958.2712.81.820.8830.3472006-10-063.902.67
0.1922007-04-302.503.80
2209+2361.1252.2914.20.391.5120.0382006-12-016.032.73
2216−0380.9016.7315.62.253.5760.0112006-08-092.559.57
2223−0521.40420.3411.713.5918.0030.1992006-10-063.2118.33
2227−0881.5602.0015.81.415.1420.1862006-07-0722.851.40
2230+1141.0378.6213.38.519.2530.2782006-02-127.365.45
2243−1230.6325.2414.81.451.7110.1612006-09-065.732.79
2251+1580.85913.7740.912.479.3930.1242006-06-158.314.72
2345−1670.57611.4715.82.812.2130.1672006-11-103.245.54
2351+4561.98621.5620.12.238.5480.1962006-05-245.357.99

Notes. Columns are as follows: (1) source name (B1950); (2) redshift z as collected by (Lister et al. 2013); (3) apparent velocity measured by (Lister et al. 2013); (4) apparent opening angle measured by Pushkarev et al. (2009); (5) flux density at the 92 cm band; (6) derived total jet power; (7) 92 cm flux density reference: 1 – De Breuck et al. (2002), 2 – Douglas et al. (1996), 3 – Ghosh et al. (1994), 4 – Gregory & Condon (1991), 5 – Kühr et al. (1979), 6 – Kühr et al. (1981), 7 – Mitchell et al. (1994), 8 – Rengelink et al. (1997); (8) core shift for frequencies 8.1 − 15.3 GHz, measured in mas (Sokolovsky et al. 2011; Pushkarev et al. 2012); (9) an epoch of the core-shift measurements by Pushkarev et al. (2012) for the year 2006 and by Sokolovsky et al. (2011) for the year 2007; (10) derived multiplicity parameter; (11) derived the Michel magnetization parameter.

Table 1.

Jet parameters and derived multiplication and magnetization parameters.

SourcezβappχappS0.3PjetReferenceΔrcoreEpochλσ
(c)(°)(Jy)(1045 erg s−1)for S0.3(mas)for Δrcore(1013)
(1)(2)(3)(4)(5)(6)(7)(8)(9)(10)(11)
0003−0660.3478.4016.32.171.0720.0352006-07-071.219.69
0106+0132.09924.3723.62.8510.5060.0052006-07-072.0223.67
0119+1150.57018.5715.62.241.8660.3472006-06-153.844.27
0133+4760.85915.3621.71.632.5480.1312006-08-093.525.80
0202+1490.40515.8916.46.252.3930.1222006-09-061.6310.90
0202+3191.46610.1513.40.762.9980.0132006-08-092.1711.13
0212+7352.3676.5516.41.545.1780.1492006-07-079.824.41
0215+0151.71525.0636.70.883.9050.0882006-04-283.757.54
0.2412006-12-017.973.54
0234+2851.20621.9919.81.453.5250.2752006-09-065.294.79
0333+3211.25913.078.03.686.7270.2792006-07-074.108.62
0336−0190.85224.4526.82.493.2660.1172006-08-092.668.68
0403−1320.57120.8016.47.624.4520.3462006-05-243.666.95
0420−0140.9165.7422.70.871.8430.2672006-10-0613.491.30
0458−0202.28613.5723.13.5411.2030.0062006-11-103.2017.26
0528+1342.07017.3416.11.025.8520.1672006-10-064.817.40
0529+0751.25418.0356.41.754.5420.1102006-08-093.837.58
0605−0850.87019.1914.01.432.3930.0962006-11-101.6611.95
0607−1570.3231.91835.12.310.9630.2402006-09-0628.460.38
0642+4493.3968.5323.40.707.6880.1102006-10-065.278.29
0730+5040.72014.0714.80.711.2080.2622006-05-244.273.21
0735+1780.4505.0421.01.811.2380.0392006-04-282.565.07
0736+0170.18913.7917.91.790.4230.0792006-06-150.828.38
0738+3130.63110.7210.51.261.4830.1832006-09-062.855.23
0748+1260.88914.5816.21.452.6520.0982006-08-092.418.69
0754+1000.26614.4013.70.740.3920.2662006-04-282.063.32
0804+4991.4361.1535.30.602.4980.0942006-10-0635.880.61
0805−0771.83741.7618.82.609.2620.2072006-05-243.1214.09
0823+0330.50512.8813.40.630.7150.1412006-06-152.124.71
0827+2430.94219.8114.60.711.8020.1502006-05-242.526.92
0829+0460.17410.1318.70.670.2120.1092006-07-071.283.82
0836+7102.21821.0812.45.0617.7520.1862006-09-063.8116.43
0851+2020.30615.1428.51.110.5630.0282006-04-281.087.69
0906+0151.02622.0817.51.543.0530.1682006-10-063.047.56
0917+6241.45312.0715.91.274.0780.1122006-08-093.957.14
0923+3920.6952.7610.83.283.0350.0422006-07-073.236.69
0945+4081.24920.2014.02.946.3020.0832006-06-151.8018.92
1036+0540.4735.726.50.750.8020.1952006-05-242.773.80
1038+0641.26510.696.71.594.3220.1062006-10-062.0214.01
1045−1880.59510.518.02.792.4620.1562006-09-062.019.45
1127−1451.18414.8916.15.638.9460.0962006-08-092.7314.77
1150+8121.25010.1115.01.393.8180.0872006-06-153.318.03
1156+2950.72524.5916.74.333.8980.1622006-09-062.1511.40
1219+0440.9660.8213.01.142.5240.1332006-05-2422.110.94
1219+2850.1039.1213.91.770.1930.1822006-02-121.114.04
0.1422007-04-300.934.87
0.1992006-11-101.193.78
1222+2160.43426.6010.83.981.9050.1802006-04-281.1414.13
1226+0230.15814.861063.723.2530.0202006-03-090.3160.77
1253−0550.53620.5814.416.566.3130.0482006-04-050.7539.84
1308+3260.99727.4818.51.422.7980.1432006-07-072.359.33
1334−1270.53916.3312.61.911.7120.2372006-10-062.615.98
1413+1350.2471.788.82.740.8120.2302006-08-096.021.64
1458+7180.9046.614.519.6413.3080.0812006-09-061.4632.22
0.1362007-03-012.1621.84
1502+1061.83917.5337.91.084.9230.0522006-07-073.598.92
1504−1660.8763.9418.41.802.7930.1482006-12-019.562.05
1510−0890.36028.0015.22.751.2230.1222006-04-280.9313.47
1514−2410.0496.397.82.060.0830.1882006-04-280.565.15
1538+1490.6068.7416.12.822.3630.0322006-06-151.6811.09
1546+0270.41412.0812.90.700.6130.0102006-08-090.8710.32
1606+1061.23219.0924.02.675.3070.0572006-07-072.1114.84
1611+3431.40029.1526.94.208.4430.0572006-06-151.7922.54
1633+3821.81329.2222.62.518.2880.1192006-09-063.0713.50
1637+5740.75113.5910.71.321.8880.1172006-05-241.928.94
1638+3981.66615.8553.80.643.1180.0072006-08-094.685.36
1641+3990.59319.2712.99.935.1380.2112006-06-152.2911.99
1655+0770.62114.775.52.362.3360.0802006-11-100.7325.39
0.0862007-06-010.7724.05
1726+4550.7172.3016.50.490.9580.0092006-09-065.182.34
1730−1300.90227.3510.46.466.5430.1742006-07-071.6719.72
1749+0960.3227.9016.81.200.6160.0612006-06-151.436.15
1751+2881.1183.8712.10.401.5520.0072006-10-063.604.62
1758+3882.0922.2117.90.181.8280.0792006-11-1013.981.42
1803+7840.68010.7918.41.922.2380.0292006-09-061.7110.80
0.0612007-05-031.989.31
1823+5680.66426.176.82.632.5280.0522006-07-070.4246.21
1828+4870.69213.077.147.7815.6030.1172006-08-091.3935.35
1849+6700.65723.0816.60.861.2280.0242006-05-240.8815.50
1908−2011.1194.3923.92.705.2120.2462006-03-0918.031.69
1928+7380.3028.179.84.811.4080.1472006-04-281.727.66
1936−1551.6575.3435.20.673.4520.2152006-07-0722.921.15
2008−1591.1804.859.70.732.4120.0082006-11-102.657.89
2022−0771.38823.2319.62.636.6720.0062006-04-051.5123.67
2121+0531.94111.6634.00.633.9920.1522006-06-1510.292.83
2128−1230.5015.995.01.471.2330.2232006-10-062.525.20
2131−0211.28419.9618.42.666.1160.0892006-08-092.3914.14
2134+0041.9325.0415.20.994.8560.1882006-07-0712.352.60
2136+1412.4274.1532.50.946.1660.0082006-09-0610.283.64
2145+0670.9992.8323.23.765.1830.0082006-10-066.974.31
2155−1520.67218.1217.62.412.4330.4052006-12-015.343.60
2200+4200.0699.9526.21.820.1280.0322006-04-050.477.30
2201+1711.07617.6613.61.002.6320.3802006-05-245.643.82
2201+3150.2958.2712.81.820.8830.3472006-10-063.902.67
0.1922007-04-302.503.80
2209+2361.1252.2914.20.391.5120.0382006-12-016.032.73
2216−0380.9016.7315.62.253.5760.0112006-08-092.559.57
2223−0521.40420.3411.713.5918.0030.1992006-10-063.2118.33
2227−0881.5602.0015.81.415.1420.1862006-07-0722.851.40
2230+1141.0378.6213.38.519.2530.2782006-02-127.365.45
2243−1230.6325.2414.81.451.7110.1612006-09-065.732.79
2251+1580.85913.7740.912.479.3930.1242006-06-158.314.72
2345−1670.57611.4715.82.812.2130.1672006-11-103.245.54
2351+4561.98621.5620.12.238.5480.1962006-05-245.357.99
SourcezβappχappS0.3PjetReferenceΔrcoreEpochλσ
(c)(°)(Jy)(1045 erg s−1)for S0.3(mas)for Δrcore(1013)
(1)(2)(3)(4)(5)(6)(7)(8)(9)(10)(11)
0003−0660.3478.4016.32.171.0720.0352006-07-071.219.69
0106+0132.09924.3723.62.8510.5060.0052006-07-072.0223.67
0119+1150.57018.5715.62.241.8660.3472006-06-153.844.27
0133+4760.85915.3621.71.632.5480.1312006-08-093.525.80
0202+1490.40515.8916.46.252.3930.1222006-09-061.6310.90
0202+3191.46610.1513.40.762.9980.0132006-08-092.1711.13
0212+7352.3676.5516.41.545.1780.1492006-07-079.824.41
0215+0151.71525.0636.70.883.9050.0882006-04-283.757.54
0.2412006-12-017.973.54
0234+2851.20621.9919.81.453.5250.2752006-09-065.294.79
0333+3211.25913.078.03.686.7270.2792006-07-074.108.62
0336−0190.85224.4526.82.493.2660.1172006-08-092.668.68
0403−1320.57120.8016.47.624.4520.3462006-05-243.666.95
0420−0140.9165.7422.70.871.8430.2672006-10-0613.491.30
0458−0202.28613.5723.13.5411.2030.0062006-11-103.2017.26
0528+1342.07017.3416.11.025.8520.1672006-10-064.817.40
0529+0751.25418.0356.41.754.5420.1102006-08-093.837.58
0605−0850.87019.1914.01.432.3930.0962006-11-101.6611.95
0607−1570.3231.91835.12.310.9630.2402006-09-0628.460.38
0642+4493.3968.5323.40.707.6880.1102006-10-065.278.29
0730+5040.72014.0714.80.711.2080.2622006-05-244.273.21
0735+1780.4505.0421.01.811.2380.0392006-04-282.565.07
0736+0170.18913.7917.91.790.4230.0792006-06-150.828.38
0738+3130.63110.7210.51.261.4830.1832006-09-062.855.23
0748+1260.88914.5816.21.452.6520.0982006-08-092.418.69
0754+1000.26614.4013.70.740.3920.2662006-04-282.063.32
0804+4991.4361.1535.30.602.4980.0942006-10-0635.880.61
0805−0771.83741.7618.82.609.2620.2072006-05-243.1214.09
0823+0330.50512.8813.40.630.7150.1412006-06-152.124.71
0827+2430.94219.8114.60.711.8020.1502006-05-242.526.92
0829+0460.17410.1318.70.670.2120.1092006-07-071.283.82
0836+7102.21821.0812.45.0617.7520.1862006-09-063.8116.43
0851+2020.30615.1428.51.110.5630.0282006-04-281.087.69
0906+0151.02622.0817.51.543.0530.1682006-10-063.047.56
0917+6241.45312.0715.91.274.0780.1122006-08-093.957.14
0923+3920.6952.7610.83.283.0350.0422006-07-073.236.69
0945+4081.24920.2014.02.946.3020.0832006-06-151.8018.92
1036+0540.4735.726.50.750.8020.1952006-05-242.773.80
1038+0641.26510.696.71.594.3220.1062006-10-062.0214.01
1045−1880.59510.518.02.792.4620.1562006-09-062.019.45
1127−1451.18414.8916.15.638.9460.0962006-08-092.7314.77
1150+8121.25010.1115.01.393.8180.0872006-06-153.318.03
1156+2950.72524.5916.74.333.8980.1622006-09-062.1511.40
1219+0440.9660.8213.01.142.5240.1332006-05-2422.110.94
1219+2850.1039.1213.91.770.1930.1822006-02-121.114.04
0.1422007-04-300.934.87
0.1992006-11-101.193.78
1222+2160.43426.6010.83.981.9050.1802006-04-281.1414.13
1226+0230.15814.861063.723.2530.0202006-03-090.3160.77
1253−0550.53620.5814.416.566.3130.0482006-04-050.7539.84
1308+3260.99727.4818.51.422.7980.1432006-07-072.359.33
1334−1270.53916.3312.61.911.7120.2372006-10-062.615.98
1413+1350.2471.788.82.740.8120.2302006-08-096.021.64
1458+7180.9046.614.519.6413.3080.0812006-09-061.4632.22
0.1362007-03-012.1621.84
1502+1061.83917.5337.91.084.9230.0522006-07-073.598.92
1504−1660.8763.9418.41.802.7930.1482006-12-019.562.05
1510−0890.36028.0015.22.751.2230.1222006-04-280.9313.47
1514−2410.0496.397.82.060.0830.1882006-04-280.565.15
1538+1490.6068.7416.12.822.3630.0322006-06-151.6811.09
1546+0270.41412.0812.90.700.6130.0102006-08-090.8710.32
1606+1061.23219.0924.02.675.3070.0572006-07-072.1114.84
1611+3431.40029.1526.94.208.4430.0572006-06-151.7922.54
1633+3821.81329.2222.62.518.2880.1192006-09-063.0713.50
1637+5740.75113.5910.71.321.8880.1172006-05-241.928.94
1638+3981.66615.8553.80.643.1180.0072006-08-094.685.36
1641+3990.59319.2712.99.935.1380.2112006-06-152.2911.99
1655+0770.62114.775.52.362.3360.0802006-11-100.7325.39
0.0862007-06-010.7724.05
1726+4550.7172.3016.50.490.9580.0092006-09-065.182.34
1730−1300.90227.3510.46.466.5430.1742006-07-071.6719.72
1749+0960.3227.9016.81.200.6160.0612006-06-151.436.15
1751+2881.1183.8712.10.401.5520.0072006-10-063.604.62
1758+3882.0922.2117.90.181.8280.0792006-11-1013.981.42
1803+7840.68010.7918.41.922.2380.0292006-09-061.7110.80
0.0612007-05-031.989.31
1823+5680.66426.176.82.632.5280.0522006-07-070.4246.21
1828+4870.69213.077.147.7815.6030.1172006-08-091.3935.35
1849+6700.65723.0816.60.861.2280.0242006-05-240.8815.50
1908−2011.1194.3923.92.705.2120.2462006-03-0918.031.69
1928+7380.3028.179.84.811.4080.1472006-04-281.727.66
1936−1551.6575.3435.20.673.4520.2152006-07-0722.921.15
2008−1591.1804.859.70.732.4120.0082006-11-102.657.89
2022−0771.38823.2319.62.636.6720.0062006-04-051.5123.67
2121+0531.94111.6634.00.633.9920.1522006-06-1510.292.83
2128−1230.5015.995.01.471.2330.2232006-10-062.525.20
2131−0211.28419.9618.42.666.1160.0892006-08-092.3914.14
2134+0041.9325.0415.20.994.8560.1882006-07-0712.352.60
2136+1412.4274.1532.50.946.1660.0082006-09-0610.283.64
2145+0670.9992.8323.23.765.1830.0082006-10-066.974.31
2155−1520.67218.1217.62.412.4330.4052006-12-015.343.60
2200+4200.0699.9526.21.820.1280.0322006-04-050.477.30
2201+1711.07617.6613.61.002.6320.3802006-05-245.643.82
2201+3150.2958.2712.81.820.8830.3472006-10-063.902.67
0.1922007-04-302.503.80
2209+2361.1252.2914.20.391.5120.0382006-12-016.032.73
2216−0380.9016.7315.62.253.5760.0112006-08-092.559.57
2223−0521.40420.3411.713.5918.0030.1992006-10-063.2118.33
2227−0881.5602.0015.81.415.1420.1862006-07-0722.851.40
2230+1141.0378.6213.38.519.2530.2782006-02-127.365.45
2243−1230.6325.2414.81.451.7110.1612006-09-065.732.79
2251+1580.85913.7740.912.479.3930.1242006-06-158.314.72
2345−1670.57611.4715.82.812.2130.1672006-11-103.245.54
2351+4561.98621.5620.12.238.5480.1962006-05-245.357.99

Notes. Columns are as follows: (1) source name (B1950); (2) redshift z as collected by (Lister et al. 2013); (3) apparent velocity measured by (Lister et al. 2013); (4) apparent opening angle measured by Pushkarev et al. (2009); (5) flux density at the 92 cm band; (6) derived total jet power; (7) 92 cm flux density reference: 1 – De Breuck et al. (2002), 2 – Douglas et al. (1996), 3 – Ghosh et al. (1994), 4 – Gregory & Condon (1991), 5 – Kühr et al. (1979), 6 – Kühr et al. (1981), 7 – Mitchell et al. (1994), 8 – Rengelink et al. (1997); (8) core shift for frequencies 8.1 − 15.3 GHz, measured in mas (Sokolovsky et al. 2011; Pushkarev et al. 2012); (9) an epoch of the core-shift measurements by Pushkarev et al. (2012) for the year 2006 and by Sokolovsky et al. (2011) for the year 2007; (10) derived multiplicity parameter; (11) derived the Michel magnetization parameter.

The Doppler factor of a flow can be also obtained through the variability method by measuring the amplitude and duration of a flare (Hovatta et al. 2009). Making an assumption that the latter corresponds to the time needed for light to cross the emitting region, and assuming that the intrinsic brightness temperature is known (from the equipartition argument), one can derive the beaming Doppler factor. We have used the variability Doppler factors obtained by Hovatta et al. (2009) for 50 objects with measured core-shifts (Pushkarev et al. 2012) instead of our original assumption for Doppler factor δ = βapp and have found that our estimates for λ and σM stay the same within a factor of 2.

We estimate the total typical accuracy of λ and σM values in Table 1 to be of a factor of a few. It is mostly due to the assumptions and simplifications introduced and, to a less of an extent, due to accuracy of observational parameters of the jets. We note that while an estimate for every source is not highly accurate, the distributions in Figs 1 and 2 should represent the sample properties well.

There are three objects in our sample that have the Michel magnetization parameter σM < 1, which means that the flow is not magnetically dominated at its base. And we have overall nine sources with σM < 2, which is in contradiction with our assumption of at least a mildly magnetized flow. This is a small fraction (9 per cent) of all 97 sources, so we feel that for the majority sources there is no contradiction of our assumptions and the resultant value for the Michel magnetization parameter.

For the highly magnetized regime, we come to a contradiction. Indeed, taking, for example, a source 0215+015, which has the Michel magnetization parameter σ0215 + 015 = σM, med, we obtain from equation (35) for a highly magnetized regime the following value:
(40)
In a highly magnetized regime, the scaling (4) holds, and for Γ = r/RL ≈ 104–105 we come to σM, mag ≈ 10−2–10−1. This is in contradiction with our assumption for a magnetized regime with initial magnetization σM > 105.

We see that the magnetization parameter λ obtained from the observed core-shift has the order of magnitude 1012–1014 which agrees with the two-photon conversion model of plasma production in a black hole magnetosphere (Blandford & Znajek 1977; Moscibrodzka et al. 2011). Thus, we obtain the key physical parameters of the jets being σM ∼ 10 and λ ∼ 1013. As a result, knowing these parameters and using a rather simple one-dimensional MHD approach, we can determine the internal structure of jets.

4 ON THE INTERNAL STRUCTURE OF JETS

As was shown by Beskin & Malyshkin (2000), Beskin & Nokhrina (2009) and Lyubarsky (2009), for well-collimated jets the one-dimensional cylindrical MHD approximation (when the problem is reduced to the system of two ordinary differential equations, see above mentioned papers for more detail) allows us to reproduce main results obtained later by two-dimensional numerical simulation (Komissarov et al. 2006; Tchekhovskoy et al. 2009; Porth et al. 2011; McKinney et al. 2012). In particular, both analytical and numerical consideration predict the existence of a dense core in the centre of a jet for low enough ambient pressure pext. Thus, knowing main parameters obtained above we can determine the transverse structure of jets using rather simple 1D analytical approximation. The only parameters we need are the Michel magnetization σM and the transverse dimension of a jet Rjet (or, the ambient pressure pext). In particular, transversal profiles of the Lorentz factor Γ, number density ne and magnetic field B can be well reproduced. In this section, we apply this approach to clarify the real structure of relativistic jets.

In Fig. 3, we present logarithmic plots of Lorentz factor and number density across the jet for λ = 1013, jet radius Rjet = 1 pc and σ = 5, 15 and 30. Fig. 4 shows logarithmic plots of poloidal and toroidal components of magnetic field across the jet with the same parameters as in Fig. 3. As we see, these results point to the existence of more dense central core in the centre of a jet. Indeed, for our parameters the number density in the centre of a jet is greater by a factor of a thousand than at the edge. However the Lorentz factor in the central core is small (see Fig. 3b). Thus, these results are in qualitative agreement with previous studies.

Transversal profile of the number density ne (a) and Lorentz factor Γ (b) in logarithmical scale for λ = 1013, jet radius Rjet = 1 pc and three different values of σ: 5 (solid line), 15 (dashed line) and 30 (dotted line).
Figure 3.

Transversal profile of the number density ne (a) and Lorentz factor Γ (b) in logarithmical scale for λ = 1013, jet radius Rjet = 1 pc and three different values of σ: 5 (solid line), 15 (dashed line) and 30 (dotted line).

Transversal profile of poloidal (a) and toroidal (b) components of magnetic field in logarithmical scale for the same parameters and line types as in Fig. 3.
Figure 4.

Transversal profile of poloidal (a) and toroidal (b) components of magnetic field in logarithmical scale for the same parameters and line types as in Fig. 3.

Knowing how the Lorentz factor on the edge of jet depends on its radius and making a simple assumption about the form of the jet, we can calculate the dependence of the Lorentz factor on the coordinate along the jet. The result is presented in Fig. 5 for the cases of parabolic |$\zeta \propto r_{\perp }^2$| and |$\zeta \propto r_{\perp }^3$| form of the jet. Here, ζ is the distance along the axis. We also assume that the jet has a radius of about 10 pc at the distance 100 pc in both cases, which corresponds to a half-opening angle of the jet θjet ≈ 0.1. According to Fig. 5, particle acceleration in the frame of the AGN host galaxy on the scales 60–100 pc has values about |$\dot{\Gamma }/\Gamma =10^{-3}$| yr−1 with very little dependence of this value on the particular form of a jet boundary. This agrees nicely with results of the VLBI acceleration study in AGN jets by Homan et al. (2009, 2014).

Dependence of Lorentz factor on coordinate along the jet in assumption of $\zeta \propto r_{\perp }^3$ (solid line) and $\zeta \propto r_{\perp }^2$ (dashed line) form of the jet.
Figure 5.

Dependence of Lorentz factor on coordinate along the jet in assumption of |$\zeta \propto r_{\perp }^3$| (solid line) and |$\zeta \propto r_{\perp }^2$| (dashed line) form of the jet.

5 STUDYING THE CAUSAL CONNECTIVITY OF THE CYLINDRICAL JET MODEL

The calculated multiplication parameter λ with the Michel magnetization parameter σM as well as the observed half-opening angle of a jet χ allow us to test causal connectivity across a jet for the cylindrical model. Every spacial point of a supermagnetosonic outflow has its own ‘Mach cone’ of causal influence. In case of a uniform flow the cone originating at the given point with its surface formed by the characteristics of a flow is a domain, where any signal from the point is known. For a non-uniform flow the cone becomes some vortex-like shape, depending on the flow property, but sustaining the property of a causal domain for a given point.

In a jet, if the characteristic inlet from any point of a set of boundary points reaches the jet axis, we say that the axis is causally connected with the boundary. On the contrary, if there is a characteristic that does not reach the axis, we have a causally disconnected flow. In the latter case, a question arises about the self-consistency of an MHD solution of the flow, since the inner parts of such a flow do not have any information about the properties of the confining medium. The examples of importance of causal connectivity in a flow and its connection with the effective plasma acceleration has been pointed out by Komissarov et al. (2009) and Tchekhovskoy et al. (2009).

In the case of the cylindrical jet model, the question of causality is even more severe. For a cylindrical model we take into account the force balance across a jet only, so the trans-field equation governing the flow becomes one-dimensional. For every initial condition at the axis its solution gives the flow profile and the position of a boundary, defined so as to contain the whole magnetic flux. Any physical value at the boundary such as, for example, the pressure, may be calculated from this solution. Or, we in fact reverse the problem, and for a given outer pressure at the boundary we find the initial conditions at the jet axis. Thus, we use the dependence of the jet properties at the axis from the conditions at the boundary. In this case, the boundary and the axis must be causally connected. In other words, for a strictly cylindrical flow the conditions at the boundary at the distance ζ0 from the jet origin must be ‘known’ to the point at the axis at the same ζ0.

In the cylindrical model, the dynamics of a flow along the jet is achieved by ‘piling up’ the described above cross cuts so as to either make the needed boundary form, or to model the variable outer pressure. In this case, the jet boundaries should be constrained by the ‘Mach cone’ following causal connection for the model to be self-consistent. Thus, we come to the following criteria: we may assert that we can neglect the jet-long derivatives in a trans-field equation if any characteristic, outlet from a boundary at ζ0, not only reaches the axis, but does it at ζ: |ζ − ζ0| ≪ ζ0.

For an axisymmetric flow, the condition of a causal connectivity across the flow may be written (Tchekhovskoy et al. 2009) in the simplest case as
(41)
where θF is a half-opening angle of a fast Mach cone at the boundary. This condition means that the characteristic from the jet boundary, locally having its half-opening angle θF with regard to the local poloidal flow velocity, reaches the axis. For an ultrarelativistic flow, θF may be defined as (Tchekhovskoy et al. 2009)
(42)
In the cylindrical approach, we can check the causal connection across the jet both by applying condition (41), and by tracking the net of characteristics, outlet from the boundary. This can be done for a different jet boundary shapes. Let us introduce the causality function
(43)
It follows from (41) and (42) that for f > 1 causal connectivity holds, and for f < 1 it does not. If a jet boundary form is given by a function ζ = ζ(r), where r is an axial radius, the half-opening jet angle is defined by
(44)
Fig. 6 shows the causality function for a paraboloidal flow (see Beskin & Nokhrina 2009) |$\zeta \propto r_{\perp }^{2}$|⁠, for a jet with a boundary shaped as |$\zeta \propto r_{\perp }^{3}$|⁠, and |$\zeta \propto r_{\perp }^{3/2}$|⁠. For the latter flow shape f < 1 for every distance. Thus, the first two outflows are causally connected, and the last one may be causally disconnected.
The causality function f for different magnetic surface shapes for Γ > Γmax/2, i.e. further the equipartition. Solid lines are shown for σM = 100 while dotted lines for σM = 10. The upper curves correspond to $\zeta \propto r_{\perp }^{3}$, the curves in the middle – to $\zeta \propto r_{\perp }^{2}$, and the lower curves – to $\zeta \propto r_{\perp }^{3/2}$.
Figure 6.

The causality function f for different magnetic surface shapes for Γ > Γmax/2, i.e. further the equipartition. Solid lines are shown for σM = 100 while dotted lines for σM = 10. The upper curves correspond to |$\zeta \propto r_{\perp }^{3}$|⁠, the curves in the middle – to |$\zeta \propto r_{\perp }^{2}$|⁠, and the lower curves – to |$\zeta \propto r_{\perp }^{3/2}$|⁠.

The cylindrical approach allows us to investigate the set of characteristics to check the causality of a flow. Let us discuss the paraboloidal flow first. We can calculate the Mach half-opening angle at each point, starting from the boundary and so tracking the exact characteristics. This half-opening angle is defined with regard to the flow velocity direction. Although in the cylindrical approach all the velocities have only the ζ-component, we may introduce the r-component by taking into account the given form of each magnetic surface. The latter is defined by the function
(45)
and for inner parts of a flow Ψ = Ψ0X. Thus, we define the angle θ of a field-line tangent to a vertical direction as |${\rm tg}\theta _{\ell } = \left|B_{r_{\perp }}\right|/\left|B_{\zeta }\right| = \mathrm{\partial} r_{\perp }/\mathrm{\partial} \zeta$|⁠, and at the given point of a flow we outlet the fast characteristic with regard to thus defined flow direction. We present in Fig. 7 (centre panel) the net of characteristics for a paraboloidal flow depicted as described above. The characteristics are calculated starting from the jet boundary towards the axis. They are parametrized by the square of fast magnetosonic Mach number MF2 at the axis at the same ζ, where the characteristic starts. This is done uniformly regarding the Mach number, thus the characteristics are plotted at different distances from each other. It can be seen that all the characteristics for parabolic jet boundary reach the axis at ζ not much greater than ζ0. The same result holds for a flow |$\zeta \propto r_{\perp }^{3}$| either (see Fig. 7, left-hand panel). On the contrary, each characteristic, depicted in Fig. 7 (right-hand panel) for a flow |$\zeta \propto r_{\perp }^{3/2}$|⁠, reaches the axis at a distance along a flow much greater than the ζ-coordinate of its origin. This suggests that the cylindrical approach is definitely not valid in this case.
The net of inbound characteristics for a power-law form of a jet boundary $\zeta \propto r_{\perp }^{3}$ (left-hand panel), $\zeta \propto r_{\perp }^{2}$ (centre panel), $\zeta \propto r_{\perp }^{3/2}$ (right-hand panel).
Figure 7.

The net of inbound characteristics for a power-law form of a jet boundary |$\zeta \propto r_{\perp }^{3}$| (left-hand panel), |$\zeta \propto r_{\perp }^{2}$| (centre panel), |$\zeta \propto r_{\perp }^{3/2}$| (right-hand panel).

6 DISCUSSION

We show that the multiplicity parameter λ, which is the ratio of number density ne of outflowing plasma to GJ number density nGJ, can be obtained from the direct observations of core shift, apparent opening angle and radio power of a jet. The formula (31) uses the following assumptions, taken from the theoretical model: (i) the acceleration process of plasma effectively stops (saturates) when there is an equipartition regime, i.e. the Poynting flux is equal to the plasma kinetic energy flux; (ii) we assume the certain power-law scalings for magnetic field B(r) and number density ne(r) as a functions of distance r (Lobanov 1998). These scalings are confirmed by Sokolovsky et al. (2011). We also see that these power laws are a good approximation from modelling the internal jet structure in Section 4.

In contrast with Lobanov (1998) and Hirotani (2005) we do not assume the equipartition regime of radiating particles with magnetic field, but the relation between the particles (radiating and non-radiating) kinetic energy and Poynting flux. We assume that only the small fraction of particles ∼1 per cent radiates (Sironi et al. 2013) and introduce the correlation between particle number density and magnetic field through the flow magnetization σ. Although for σ ∼ 1 both approaches give effectively the same relation between particle number density and magnetic field, for highly magnetized regime σ ≫ 1 our approach yields the different result. Probing both the equipartition regime σ ∼ 1 and highly magnetized regime σ ≫ 1 at parsec scales we conclude that the latter does not hold.

Using the obtained Michel magnetization parameter σM, one can easily explain the observationally derived values 2Γχ ≈ 0.1 (Clausen-Brown et al. 2013; Zamaninasab et al. 2014), where Γ is Lorentz factor of bulk plasma motion and χ is a jet half-opening angle. Indeed, as was found by Tchekhovskoy, McKinney & Narayan (2008), Beskin (2010), 2Γχ ≈ 1 in the whole domain where Γ ≪ σ, independent of the collimation geometry. This implies that 2Γχ ≈ 1 up to the distance from the origin whence the transverse dimension of a jet Rjet/RL = σM. At larger distances Γ remains practically constant, but for a parabolic geometry the opening angle decreases with the distance ζ as ζ−1/2 ≈ r−1/2. As a result, one can write down
(46)
This result is in agreement with the criteria of casual connectivity across a jet. Indeed, for an outflow with an equipartition between the Poynting and particle energy flux, we can write down
(47)
for a boundary causally connected with an axis. In Section 5, we have shown that for flows collimated better than a parabola, casuality connectivity across the jet holds further, i.e. for Γ > Γmax/2.

7 SUMMARY

The analysis of the frequency dependence of the observed shift of the core of relativistic AGN jets allows us to determine physical parameters of the jets such as the plasma number density and the magnetic field inside the flow. We have estimated the multiplicity parameter λ to be of the order of 1012–1015. It is consistent with the Blandford–Znajek model (Blandford & Znajek 1977) of the electron–positron generation in the magnetosphere of the black hole (see Moscibrodzka et al. 2011, as well). These values are in agreement with the particle number density ne which was found independently by Lobanov (1998).

As the transverse jet structure depends strongly on the flow regime, whether it is in equipartition or magnetically dominated, it is important to know the relation between the observed and maximum Lorentz factor. The Michel magnetization parameter σM is equal to the maximum Lorentz factor of plasma bulk motion. Typical derived values of σM ≲ 30, in agreement with the Lorentz factor estimated from VLBI jet kinematics (e.g. Cohen et al. 2007; Lister et al. 2009b, 2013) and radio variability (Jorstad et al. 2005; Hovatta et al. 2009; Savolainen et al. 2010). This implies that a flow is in the saturation regime. Since for strongly collimated flow the condition of causal connection is fulfilled (see e.g. Komissarov et al. 2009; Tchekhovskoy et al. 2009), the internal structure of an outflow can be modelled within the cylindrical approach (Beskin & Malyshkin 2000; Beskin & Nokhrina 2009). It has been shown that the results of the modelling, such as Lorentz factor dependence on the jet distance, are in good agreement with the observations. In particular, the relative growth of Lorentz factor |$\dot{\Gamma }/\Gamma$| with the distance along the axis is slow for the jets in saturation regime, having the magnitude ∼10−3 yr−1. This result may account for the recent measurements of acceleration in AGN jets (Homan et al. 2014).

We plan to address the following points in a separate paper: (i) the role of the inhomogeneity of the magnetic field and particle number density in a core, (ii) the action of the radiation drag force (Li, Begelman & Chiueh 1992; Beskin, Zakamska & Sol 2004; Russo & Thompson 2013), (iii) the possible influence of mass loading (Komissarov 1994; Derishev et al. 2003; Stern & Poutanen 2006) on the jet magnetization and dynamics.

We would like to acknowledge E. Clausen-Brown, D. Gabuzda, M. Sikora, A. Lobanov, T. Savolainen, M. Barkov and the anonymous referee for useful comments. We thank the anonymous referee for suggestions which helped to improve the paper. This work was supported in part by the Russian Foundation for Basic Research grant 13-02-12103. YYK was also supported in part by the Dynasty Foundation. This research has made use of data from the MOJAVE data base that is maintained by the MOJAVE team (Lister et al. 2009a), and data accumulated by the CATS data base (Verkhodanov et al. 1997). This research has made use of NASA's Astrophysics Data System.

REFERENCES

Appl
S.
Camenzind
M.
A&A
,
1992
, vol.
256
pg.
354
Benford
G.
ApJ
,
1981
, vol.
247
pg.
792
Beskin
V. S.
MHD Flows in Compact Astrophysical Objects
,
2009
Berlin
Springer-Verlag
Beskin
V. S.
Phys.-Usp.
,
2010
, vol.
53
pg.
1199
Beskin
V. S.
Malyshkin
L. M.
Astron. Lett.
,
2000
, vol.
26
pg.
4
Beskin
V. S.
Nokhrina
E. E.
MNRAS
,
2009
, vol.
397
pg.
1486
Beskin
V. S.
Istomin
Ya. N.
Pariev
V. I.
SvA
,
1992
, vol.
36
pg.
642
Beskin
V. S.
Zakamska
N. L.
Sol
H.
MNRAS
,
2004
, vol.
347
pg.
587
Bisnovatyi-Kogan
G. S.
Ap&SS
,
2007
, vol.
311
pg.
287
Blandford
R. D.
MNRAS
,
1976
, vol.
176
pg.
465
Blandford
R. D.
Königl
A.
ApJ
,
1979
, vol.
232
pg.
34
Blandford
R. D.
Znajek
R. L.
MNRAS
,
1977
, vol.
179
pg.
433
Bogovalov
S. V.
Astron. Lett.
,
1995
, vol.
21
pg.
565
Cavagnolo
K. W.
McNamara
B. R.
Nulsen
P. E. J.
Carilli
C. L.
Jones
C.
Bîrzan
L.
ApJ
,
2010
, vol.
720
pg.
1066
Clausen-Brown
E.
Savolainen
T.
Pushkarev
A. B.
Kovalev
Y. Y.
Zensus
J. A.
A&A
,
2013
, vol.
558
pg.
A144
Cohen
M. H.
Lister
M. L.
Homan
D. C.
Kadler
M.
Kellermann
K. I.
Kovalev
Y. Y.
Vermeulen
R. C.
ApJ
,
2007
, vol.
658
pg.
232
De Breuck
C.
Tang
Y.
de Bruyn
A. G.
Röttgering
H.
van Breugel
W.
A&A
,
2002
, vol.
394
pg.
59
Derishev
E. V.
Aharonian
F. A.
Kocharovsky
V. V.
Kocharovsky
Vl. V.
Phys. Rev. D
,
2003
, vol.
68
pg.
043003
Douglas
J. N.
Bash
F. N.
Bozyan
F. A.
Torrence
G. W.
Wolfe
C.
AJ
,
1996
, vol.
111
pg.
1945
Ghosh
T.
Gopal-Krishna Rao
A. P.
A&AS
,
1994
, vol.
106
pg.
29
Goldreich
P.
Julian
W. H.
ApJ
,
1969
, vol.
157
pg.
869
Gould
R. J.
A&A
,
1979
, vol.
76
pg.
306
Gregory
P. C.
Condon
J. J.
ApJS
,
1991
, vol.
75
pg.
1011
Hardee
P. E.
Norman
M. L.
ApJ
,
1988
, vol.
334
pg.
70
Heyvaerts
J.
Norman
C.
ApJ
,
2003
, vol.
596
pg.
1240
Hirotani
K.
ApJ
,
2005
, vol.
619
pg.
73
Hirotani
K.
Okamoto
I.
ApJ
,
1998
, vol.
497
pg.
563
Homan
D. C.
Kadler
M.
Kellermann
K. I.
Kovalev
Y. Y.
Lister
M. L.
Ros
E.
Savolainen
T.
Zensus
J. A.
ApJ
,
2009
, vol.
706
pg.
1253
Homan
D. C.
Lister
M. L.
Kovalev
Y. Y.
Pushkarev
A. B.
Savolainen
T.
Kellermann
K. I.
Richards
J. L.
Ros
E.
ApJ
,
2014
 
preprint (arXiv:1410.8502)
Hovatta
T.
Valtaoja
E.
Tornikoski
M.
Lähteenmäki
A.
A&A
,
2009
, vol.
498
pg.
723
Istomin
Ya. N.
Pariev
V. I.
MNRAS
,
1994
, vol.
267
pg.
629
Jorstad
S. G.
et al.
AJ
,
2005
, vol.
130
pg.
1418
Komatsu
E.
et al.
ApJS
,
2009
, vol.
180
pg.
330
Komissarov
S.
MNRAS
,
1994
, vol.
269
pg.
394
Komissarov
S.
Barkov
M.
Vlahakis
N.
Königl
A.
MNRAS
,
2006
, vol.
380
pg.
51
Komissarov
S. S.
Vlahakis
N.
Königl
A.
Barkov
M. V.
MNRAS
,
2009
, vol.
394
pg.
1182
Kovalev
Y. Y.
Lobanov
A. P.
Pushkarev
A. B.
Zensus
J. A.
A&A
,
2008
, vol.
483
pg.
759
Krolik
J.
Active Galactic Nuclei: From the Central Black Hole to the Galactic Environment
,
1999
Princeton, NJ
Princeton Univ. Press
Kronberg
P. P.
Lovelace
R. V. E.
Lapenta
G.
Colgate
S. A.
ApJ
,
2011
, vol.
741
pg.
L15
Kühr
H.
Nauber
U.
Pauliny-Toth
I. I. K.
Witzel
A.
MPIfR Preprint no. 55, A Catalogue of Radio Sources
,
1979
Bonn
Max-Planck Institut für Radioastronomie
pg.
1
Kühr
H.
Witzel
A.
Pauliny-Toth
I. I. K.
Nauber
U.
A&AS
,
1981
, vol.
45
pg.
367
Li
Z.-Y.
Begelman
M.
Chiueh
T.
ApJ
,
1992
, vol.
384
pg.
567
Lister
M. L.
et al.
AJ
,
2009a
, vol.
137
pg.
3718
Lister
M. L.
et al.
AJ
,
2009b
, vol.
138
pg.
1874
Lister
M. L.
et al.
AJ
,
2013
, vol.
146
pg.
120
Lobanov
A. P.
A&A
,
1998
, vol.
330
pg.
79
Lovelace
R. W. E.
Nature
,
1976
, vol.
262
pg.
649
Lyubarsky
Yu. E.
ApJ
,
2009
, vol.
308
pg.
1006
Lyutikov
M.
Pariev
V. I.
Gabuzda
D. C.
MNRAS
,
2005
, vol.
360
pg.
869
McKinney
J. C.
Tchekhovskoy
A.
Blanford
R. D.
MNRAS
,
2012
, vol.
423
pg.
2083
Marscher
A. P.
ApJ
,
1983
, vol.
264
pg.
296
Mestel
L.
Stellar Magnetism
,
1999
Oxford
Clarendon Press
Michel
F. C.
ApJ
,
1969
, vol.
158
pg.
727
Mitchell
K. J.
et al.
ApJS
,
1994
, vol.
93
pg.
441
Moscibrodzka
M.
Gammie
C. F.
Dolence
J. C.
Shiokawa
H.
ApJ
,
2011
, vol.
735
pg.
9
O'Sullivan
S. P.
Gabuzda
D. C.
MNRAS
,
2009
, vol.
400
pg.
26
Okamoto
I.
PASJ
,
2009
, vol.
61
pg.
971
Pariev
V. I.
Istomin
Ya. N.
Beresnyak
A. R.
A&A
,
2003
, vol.
403
pg.
805
Porth
O.
Fendt
Ch.
Meliani
Z.
Vaidya
B.
ApJ
,
2011
, vol.
737
pg.
42
Punsly
B.
Black Hole Gravitohydromagnetics
,
2001
Berlin
Springer-Verlag
Pushkarev
A. B.
Kovalev
Y. Y.
Lister
M. L.
Savolainen
T.
A&A
,
2009
, vol.
507
pg.
L33
Pushkarev
A. B.
Hovatta
T.
Kovalev
Y. Y.
Lister
M. L.
Lobanov
A. P.
Savolainen
T.
Zensus
J. A.
A&A
,
2012
, vol.
545
pg.
A113
Rengelink
R. B.
Tang
Y.
de Bruyn
A. G.
Miley
G. K.
Bremer
M. N.
Roettgering
H. J. A.
Bremer
M. A. R.
A&AS
,
1997
, vol.
124
pg.
259
Russo
M.
Thompson
Ch.
ApJ
,
2013
, vol.
767
pg.
142
Savolainen
T.
Homan
D. C.
Hovatta
T.
Kadler
M.
Kovalev
Y. Y.
Lister
M. L.
Ros
E.
Zensus
J. A.
A&A
,
2010
, vol.
512
pg.
A24
Sironi
L.
Spitkovsky
A.
Arons
J.
ApJ
,
2013
, vol.
771
pg.
54
Sokolovsky
K. V.
Kovalev
Y. Y.
Pushkarev
A. B.
Lobanov
A. P.
A&A
,
2011
, vol.
532
pg.
A38
Stern
B. E.
Poutanen
J.
MNRAS
,
2006
, vol.
372
pg.
1217
Tchekhovskoy
A.
McKinney
J.
Narayan
R.
MNRAS
,
2008
, vol.
388
pg.
551
Tchekhovskoy
A.
McKinney
J.
Narayan
R.
ApJ
,
2009
, vol.
699
pg.
1789
Verkhodanov
O. V.
Trushkin
S. A.
Andernach
H.
Cherenkov
V. N.
Hunt
G.
Payne
H.
ASP Conf. Ser. Vol. 125, Astronomical Data Analysis Software and Systems VI
,
1997
San Fransisco
Astron. Soc. Pac.
pg.
322
Zamaninasab
M.
Clausen-Brown
E.
Savolainen
T.
Tchekhovskoy
A.
Nature
,
2014
, vol.
510
pg.
126
Zdziarski
A. A.
Sikora
M.
Pjanka
P.
Tchekhovskoy
A.
MNRAS
,
2014
 
preprint (arXiv:1410.7310)

APPENDIX A: MAGNETIZATION PARAMETER

The standard Grad–Shafranov approach for MHD flows uses the following energy integral E, conserved at the magnetic flux surface Ψ:
(A1)
Here magnetic field in spherical coordinates {r; θ; φ} is defined by
(A2)
electric field
(A3)
ΩF is a rotational velocity and a function of magnetic flux Ψ, I is a current and
(A4)
Here we assume the flow to be cold, so the particle enthalpy is simply its rest mass. The Poynting vector
(A5)
The particle kinetic energy flux
(A6)
Thus, we can introduce the magnetization parameter, variable along the flow, as
(A7)
On the other hand, there is Michel's magnetization parameter σM, which has a meaning of magnetization at the base of an outflow, is constant, and is defined by
(A8)

APPENDIX B: MULTIPLICITY PARAMETER

The continuity of relativistic plasma number density flux through two cuts, perpendicular to bulk flow, can be written as
(B1)
Here we assume the flow velocity equal to c. Using the expression for the GJ concentration
(B2)
and expression for the full losses due to currents flowing in a magnetosphere
(B3)
we get
(B4)