Issue 
A&A
Volume 612, April 2018



Article Number  A63  
Number of page(s)  22  
Section  Astrophysical processes  
DOI  https://doi.org/10.1051/00046361/201731793  
Published online  27 April 2018 
Nonradial and nonpolytropic astrophysical outflows
X. Relativistic MHD rotating spine jets in Kerr metric
^{1}
Laboratoire Univers et Théories, Observatoire de Paris  PSL, UMR 8102 du CNRS, Université Paris Diderot,
92190
Meudon,
France
email: loic.chantry@obspm.fr
^{2}
Section of Astrophysics, Astronomy and Mechanics, Department of Physics, University of Athens,
Panepistimiopolis Zografos, Athens
15783,
Greece
email: christophe.sauty@obspm.fr
Received:
17
August
2017
Accepted:
15
December
2017
Context. Highresolution radio imaging of active galactic nuclei (AGN) has revealed that the jets of some sources present superluminal knots and transverse stratification. Recent observational projects, such as ALMA and γray telescopes, such as HESS and HESS2 have provided new observational constraints on the central regions of rotating black holes in AGN, suggesting that there is an inner or spinejet surrounded by a disk wind. This relativistic spinejet is likely to be composed of electronpositron pairs extracting energy from the black hole and will be explored by the future γray telescope CTA.
Aims. In this article we present an extension to and generalization of relativistic jets in Kerr metric of the Newtonian meridional selfsimilar mechanism. We aim at modeling the inner spinejet of AGN as a relativistic light outflow emerging from a spherical corona surrounding a Kerr black hole and its inner accretion disk.
Methods. The model is built by expanding the metric and the forces with colatitude to first order in the magnetic flux function. As a result of the expansion, all colatitudinal variations of the physical quantities are quantified by a unique parameter. Unlike previous models, effects of the light cylinder are not neglected.
Results. Solutions with high Lorentz factors are obtained and provide spinejet models up to the polar axis. As in previous publications, we calculate the magnetic collimation efficiency parameter, which measures the variation of the available energy across the field lines. This collimation efficiency is an integral part of the model, generalizing the classical magnetic rotator efficiency criterion to Kerr metric. We study the variation of the magnetic efficiency and acceleration with the spin of the black hole and show their high sensitivity to this integral.
Conclusions. These new solutions model collimated or radial, relativistic or ultrarelativistic outflows in AGN or γray bursts. In particular, we discuss the relevance of our solutions to modeling the M 87 spinejet. We study the efficiency of the central black hole spin to collimate a spinejet and show that the jet power is of the same order as that determined by numerical simulations.
Key words: black hole physics / magnetohydrodynamics (MHD) / relativistic processes / galaxies: jets
© ESO 2018
1 Introduction
Active galactic nuclei (AGN) jets are now recognized as being multicomponent outflows related to accretion onto a supermassive black hole (SMBH). The outer kiloparsec/megaparsecscale lobes are fed by a powerful hadronic plasma, which most likely originates from the Keplerian accretion disk via the magnetocentrifugal launching mechanism, Blandford & Payne (1982). This hadronic population could be responsible for the second component peaking in the γray band of the spectral energy distribution (SED) of some blazars. As explained in Böttcher et al. (2013), this is possible if a relativistic jet of protons contributes significantly to the radiative output through proton synchrotron emission or photopion production. However the purely hadronic synchrotron models for blazars present a major problem due to the requirement for having very high powers in the jets. Indeed, hadronic processes are very inefficient (Sikora 2011). They are only energetically favorable for highfrequency synchrotron peak (HSP) blazars, where the jet is likely to be highly magnetized (Petropoulou & Dermer 2016).
The alternative approach to explaining the origin of the highenergy emission is to consider leptonic models. There, the radiative output throughout the electromagnetic spectrum is assumed to be dominated by leptons both at low and high frequencies. Inverse Compton scattering of soft photons (IC) by relativistic nonthermal electrons in the jet is supposed to be the most probable mechanism for γray production, at least for strongline and lowfrequency synchrotron peak – LSP – blazars. Until recently, leptonic models in quasi steadystate were very successful in modeling the SED for almost all classes of blazars (e.g., Celotti & Ghisellini 2008). However leptonic models with a onezone component where only one portion of the jet dominates the emission are now questioned. Mostly, they fail to produce the extremely high bulk Lorentz factors in a very compact emission region required for ultrafast variability of some TeV blazars (Begelman et al. 2008). Despite several proposed mechanisms mostly based on smallscale inhomogeneities in the jet (see references in Vovk & Babić 2015), it is not yet possible to know the location of the variable γray emission; it may be at the base of the jet or up to parsecscale distances from the central black hole. Moreover most of the sources in which ultrarapid variability events have been observed are BL Lac objects emitting veryhighenergy radiation.
For BL Lacs the application of the onezone synchrotronself Compton model to the SED implies that the jet is weakly magnetized and the emitting region is far from equipartition. Both conditions are required if magneticreconnectionpowered emission is at the origin of the ultrarapid variability (Tavecchio 2016). A most important constraint comes from the launching and the acceleration processes, which should lead to almost equipartition between the magnetic and the kinetic energy fluxes as can be seen in other blazars. In order to unify the BL Lac and radiogalaxy populations, Ghisellini et al. (2005) already proposed a structured jet model with two components; a faster core (the spine) surrounded by a slower sheath or layer. This two component model is able to explain peculiar features of TeVemitting BL Lacs, such as the absence of fast superluminal components and the presence of a limbbrightened radio structure (Giroletti et al. 2004; Nagai et al. 2014). Tavecchio & Ghisellini (2016) showed that this structured jet model gives an extra source of soft photons intervening in the IC emission for BL Lacs. It is due to the radiative interplay between the two components and allows to reproduce the emission in equipartition conditions. This spinesheathjet structure has also been explored in Sikora et al. (2016) for strongline blazars in order to explain why γray variations are often observed to have much larger amplitudes than the corresponding optical variations, as well as other observed features of γray flares. For the nearby radiogalaxy M 87, a direct detection of different motions was possible and an extended mildly relativistic flow is seen surrounding a relativistic central jet (Mertens et al. 2016).
The possible seed photon population at the origin of IC scattering in AGNs is questionable. The dominant seed source depends critically on the location of the primary emitting region in the jet (Finke 2016). Radio array observations of individual blazars are used to localize the γray emission site among the superluminal components seen along the jet. However, they cannot put strong constraints on localizing the emissionregions of the highenergy flares. More theoretical and numerical modelings and carefully analyzed polarization data are needed (see Kravchenko et al. 2016 and references therein). We note also the treatment used in Hada et al. (2011) for locating the central black hole relatively close to the jet baseradio core in M 87 using the dependence of the core position on frequency. There are cases where the spinesheath model should work. The leptons in superluminal components after leaving the black hole corona interact with the sheath where IC scattering of photons occurs and produce the observed γray activity (Marscher et al. 2010; Casadio et al. 2015).
For AGN jets, twocomponent models have been theoretically studied in Sol et al. (1989) and further developed by several authors (e.g., Fabian & Rees 1995). In those studies the spine jet is the inner leptonic component (electronpositron pairs), which is self collimated inside a hadronic disk wind, and should be separated from the wind by a forcefree, empty region. On the other hand, it has been shown numerically that the high collimation efficiency of the outerdisk wind is sufficient to prevent decollimation from the inertia of the inner relativistic plasma, in the case where the initial magnetic field lines going out from the launching region, that is, the black hole corona, are radial (e.g., Gracia et al. 2009). Moreover synthetic synchrotron maps built from such magnetohydrodynamic (MHD) models can reproduce the opening angle of M 87 up to a large distance from the black hole, and the change between centerbrightening near the core to limbbrightening further away.
The standard accretionjet model includes an advectiondominated accretion flow (ADAF) which takes place at least in the inner part of the accretion disk in AGNs, except for very high accretion efficiency (Narayan & Yi 1994). The presence of this ADAF structure in the disk center does not prevent the launching of the internal jet as can be seen from observations near the black hole for radiogalaxies (see the Faraday rotation measures from observations of M 87 in Feng et al. 2016). In the highaccretion regime for AGNs, where a standard thin disk is found and powerful outflows are observed, both an internal jet and an outer disk wind coexist.
Numerical simulations in the framework of general relativity have been performed to model the inner jet formation assuming a largescale magnetic field right from the beginning, mostly using monopolar configuration anchored in the ergosphere. Komissarov (2007) was able to ensure the viability of the Blandford–Znajek mechanism (Blandford & Znajek 1977) and the possibility of building a relativistic particle outflow. Similar relativistic simulations of jets were obtained by McKinney & Blandford (2009) where the central Poyntingfluxdominated forcefree jet is selfconsistently feeding the funnel of the accretiondisk wind created along the axis by the very high toroidal field in a magnetically arrested disk (MAD). McKinney et al. (2012) and Tchekhovskoy et al. (2011) have also derived a scaling law between the magnetic flux of the field threading the black hole and the massaccretion rate on the horizon. This has been confirmed observationally by Zamaninasab et al. (2014). In addition, the jet power, directly related to the magnetic flux on the horizon, was calculated in their simulations, allowing for the derivation of a net flow efficiency between the accretion and the jet.
Beyond the launching region, the internal structure of magnetized and relativistic jets can be probed by relativistic MHD (RMHD) simulations, as has been done by Martí et al. (2016). They characterized the internal jet structure of overpressured, steady jets in connection with the dominant energy (internal, restmass, and magnetic). They showed that transverse equilibrium with a significant toroidal magnetic component implies a structure with a central spine and a surrounding layer with lower thermal and total pressures. Thus, we have strong clues, both from observations and simulations, that relativistic jets have a transverse structure. On the other hand, the jet launching mechanism associated with magnetic flux threading the black hole seems to produce only very light spine jets with rarified gas such that the inner outflow is forcefree. The inner beam would be so light that it would be invisible at large distances. An alternative issue, in order to model the spine jet, is to extend the selfsimilar solutions for jet outflows in young stellar objects to general relativity. This model can be derived analytically and used as initial conditions in numerical simulations.
A selfsimilar model of the spine jet has been proposed by Meliani et al. (2006) for nonrotating black holes in Schwarzschild metric, extending the previous Newtonian selfsimilar model of Sauty & Tsinganos (1994). The Meliani et al. (2010) model aimed to explain the dichotomy between FRI and FRII sources in terms of magnetic collimation efficiency. Then, Globus et al. (2014), partially generalized this model to a Kerr metric. However they fail to give an exact generalization of the magnetic collimation criterion linked to the efficiency of the magnetic rotator. In both cases, the authors assume that the light cylinder is sufficiently far from the spine jet to have its effects neglected. Physically, this means that both models cannot produce solutions consistent across the light cylinder.
In the present paper, we present a new extension of the nonrelativistic meridional selfsimilar solutions of Sauty & Tsinganos (1994), both in Schwarzschild and Kerr metrics. This model can produce solutions for relativistic jets emerging from a spherical corona surrounding the central part of a Kerr black hole and its inner accretion disk.
Unlike previous models in a Schwarzschild or a Kerr metric, this new model for a nonrotating, as well as for a rotating black hole includes a selfconsistent crossing of the light cylinder in the solution. The Alfvén Mach surfaces are not spherically symmetric. However, the Alfvén transition surface is spherical and includes the light cylinder.
In Sect. 2 we summarize the ideal MHD equations in the 3+1 formalism and in Sect. 3 construct the model and derive its equations. In Sect. 4 we describe how to solve the equations and illustrate this procedure with various new solutions in a Schwarzschild metric. Then, four different solutions in a Kerr metric are presented in Sect. 5 with high Lorentz factor and different geometries potentially applicable to AGNs and γray bursts (GRB). Some characteristics of those solutions of the model are discussed as well as how the magnetic collimation evolves with the rotation of the black hole.
2 Steady axisymmetric relativistic MHD outflows
2.1 Kerr metric
The first step in building selfsimilar solutions in relativistic flows is to define the metric. In fact the central massive black hole dominates the gravitational field in the near regions and determines completely the metric field. Thus, in Kerr metric, the geodesics are defined as follows. (1)
We have the usual notations of the elements,
where and .
We note that is the angular momentum of the massive central object, is its mass, h is the lapse function, ω is the angular velocity of zero angular momentum observers (ZAMO) and we use a for the lengthscale related to the angular momentum of the black hole (Kerr scale). We can define the dimensionless spin of the black hole a_{H} in units of the gravitational radius r_{s}∕2 such that, a_{H} = 2a∕r_{s}. Furthermore, β is the shift vector.
The lapse function h, the angular velocity ω of zero angular momentum observers (ZAMO) and the shift vector coordinates can be written as:
with . The corresponding line elements for the Kerr metric are given in Appendix A.
2.2 Maxwell’s equations
The next step is to define the electromagnetic field in this metric. Using covariant derivatives, we can write Maxwell’s equations in Kerr space, assuming stationarity and axisymmetry,
where is the space orthonormal basis. We note that all quantities in the above equations are given in the ZAMO frame. We can split all vector fields in a poloidal component in the meridional plane and a toroidal one along the azimuthal direction. The poloidal magnetic field B_{p} can be expressed in terms of the magnetic flux function A, (11)
and using Faraday’s law, we get the electric field, (12)
We can introduce the electric potential, Φ, but the electric field E is not directly proportional to the gradient of the electric potential, (13)
The condition of ideal MHD for infinite electrical conductivity leads to (14)
We note that ∇ is the covariant derivative on a space hypersurface; see Appendix B for the expression of its coordinates.
2.3 Equations of motion
In the Kerr metric, the 3+1 formalism gives the equation for mass conservation, the Euler equation, and the energy conservation, respectively,
Here is the toroidal component of the bulk flow speed as seen by the ZAMO. The factor γ is the bulk Lorentz factor, ρ_{0} is the mass density, and ξc^{2} the specific enthalpy measured in the comoving frame of the outflow, that contains kinetic enthalpy of perfect relativistic gas ξ_{K} and some heating term Q∕c^{2}. (18)
For the kinetic enthalpy we use the Taub–Matthews approximation of ideal fluid equation of state; for more details see Taub (1948), Meliani et al. (2004) and Mignone et al. (2005). (19)
The energy conservation has been derived in the frame of the ZAMO and equivalently the first law of thermodynamics can be obtained by projecting the conservation of the energymomentum tensor along the fluid 4velocity but in the comoving frame. Assuming infinite conductivity, the contribution of the electromagnetic field is null and only the thermal energy affects the variation of the enthalpy of the fluid, giving, (20)
2.4 Constants of motion
Under the assumptions of steadiness and axisymmetry, the MHD equations in general relativity can be partially integrated to yield severalfield/streamline constants (Beskin 2010). We already deduced those constants including the magnetic field with the same formalism in Cayatte et al. (2014). Here we present the derivation of the equations, following the notations of Tsinganos (1982), in order to compare with previous selfsimilar models, for example, Meliani et al. (2006) and give the choice of the first integrals.
Steady and axisymmetric flows are characterized by a function A that defines the geometry of the magnetic flux surfaces. In the poloidal plane, field lines are lines of constant magnetic flux A and first integrals will be functions of A, among which the mass flux Ψ. The poloidal velocity can be expressedin terms of Ψ, (21)
The frozenin condition for ideal MHD flows, gives, in the toroidal direction, combined with Eq. (11),
where Ψ_{A} ≡dΨ∕dA is the magnetic to mass flux ratio.
The poloidal components of the law of flux freezing (Eq. (14)) give in turn the isorotation law, (22)
where Ω(A) ≡ cdΦ∕dA, which is the isorotation frequency, constant along each magnetic flux tube.
By integrating the Euler equation in the toroidal direction, we get the conservation of the angular momentum flux L(A), (23)
The last equation to integrate is the energy conservation. In other words, we may take the Euler equation projected along the time axis of the 3+1 decomposition, and integrate it under the hypothesis of steadiness, (24)
2.5 Toroidal fields
Using the three last integrals of motion, we may express the toroidal components of the velocity and the magnetic fields and the enthalpy density as functions of these first integrals and the poloidal components. Using the standard procedure of inversion we get
where we have defined the poloidal Alfvén Mach number, (28)
This definition of the poloidal Alfvén Mach number is consistent with the definition used by Meliani et al. (2006) and includes the lapse function. This is also the definition taken by Breitmoser & Camenzind (2000) because the velocity, hV_{p}, calculated with the universal time is continuous across the event horizon.
The numerator and denominator of Eq. (25) are zero at the Alfvén transition surface, if the following two equations are satisfied, (29) (30)
The denominators of Eqs. (26) and (27) are identical to the one of Eq. (25) and there numerators are a linear combination of Eqs. (29) and (30). Hence, the numerators and the denominators of Eqs. (26) and (27) are also zero at the Alfvén transition surface.
We can reformulate the above equations by changing variables. We rescale the cylindrical radius with ch∕(Ω − ω) leading to the dimensionless cylindrical radius x and introduce the parameter x_{MR}, (31)
The second condition, Eq. (30), at the Alfvén transition surface becomes, keeping the first one unchanged, (35)
In Kerr metric, the parameter is an extension of defined by Meliani et al. (2006). It measures the amount of energy carried by the electromagnetic field. This is the energy flux of the magnetic rotator (MR) divided by the total energy flux of the outflow in the corotating frame, . This new parameter x_{MR}, conversely to the previous x_{A}, is not any more constant along a field line, since ω is not an integral of the motion.
3 Model equations
3.1 Angular expansion
The MHD equations and the metric constitute a coupled set of highly nonlinear equations that cannot be solved analytically. The approach followed so far for Newtonian flows has been to look for solutions with separable variables in the frame of selfsimilarity. However, this technique cannot be applied in the frame of general relativity due to the complexity of the metric even for the simpler cases of a Schwarzschild or a Kerr metric. Instead, we may model the jet close to its symmetry axis, that is, to describe the spine jet, by expanding all variables with sin θ to second order.
Along the polar axis where ϖ and θ go to zero, we may define the spherical Alfvén radius to be the distance r_{⋆} from the center where the Alfvén transition surface condition, , applies. The subscript ⋆ denotes the value of a physical quantity at the Alfvén transition surface, along the polar axis. We shall use this location to write all our quantities in dimensionless form. Thus, the dimensionless spherical radius is (36)
At R = 1, the velocity is V_{⋆}, the magnetic field B_{⋆}, the density ρ_{⋆}, the enthalpy ξ_{⋆} and the lapse function h_{⋆}. Because of the Alfvén transition along the polar axis, we have, (37)
Thus the dimensionless magnetic flux function α is defined as, (38)
Moreover we can expand to the second order the metric of the system in dimensionless form using the characteristic dimensions of the system defined in the previous Section. This introduces the two following new parameters, (39)
which are respectively the Schwarzschild radius in units of the Alfvén radius and the dimensionless black hole spin.
Another dimensionless parameter is needed to describe the gravitational potential, as in the classical model. This parameter ν represents the escape speed at the Alfvén point along the polar axis in units of V_{⋆}. Then, the value of V_{⋆} is fixed by the following condition, (40)
Thus, tosecond order in sin θ the ZAMO angular velocity and the lapse function are written as,
In order to simplify our notation, we define the lapse function along the polar axis, (43)
and the polar shift of the metric, (44)
See Appendix A for details.
It will be useful to introduce the dimensionless polar shift function (see also Eq. (59)), (45)
We also expand the magnetic flux function to second order in sin θ. The magnetic flux is an even function which is zero along the polar axis due to axisymmetry and because of the symmetry around the equatorial plane. Thus all odd orders are zero and the first nonvanishing even order is the second order in colatitude. If we keep the lowest order in the expansion we get (46)
where f is the inverse of the classical expansion factor for solar coronal holes (see Tsinganos & Sauty 1992). This expansion, similarly to the classical self similar model of Sauty & Tsinganos (1994), is equivalent to a hypothesis of separation of the variables in the magnetic flux function.
Thus from Eq. (A.17), the cylindrical radius can also be seen as an expansion in the magnetic flux. This is physically more meaningful as the magnetic flux is constant on a given mass flux tube. Moreover, several free integrals depend solely on this magnetic flux. We define the dimensionless cylindrical radius G in units of the polar Alfvén radius as, (47)
The cylindrical radius can be written in the various following forms, (48)
We can also write the metric as an expansion in α (see also Appendix A),
Of course we can always reverse our point of view and go back to the expansion in θ. This would be the case if we wanted to use the steady analytical solution as initial conditions for numerical simulations.
We can parametrize the geometry of the flux tubes with the logarithm derivative of f denoted F, (51)
The angle χ of the magnetic poloidal field line with the radial direction (see Sauty et al. 1999) is given in our metric by, (52)
3.2 Choice of the Alfvén surface and pressure
We can expand all physical quantities to the first order in α. Thus the Alfvén number is given by, (53)
Contrary to previous selfsimilar models, the Alfvén number cannot be spherically symmetric because of the presence of the cylindrical radius in units of the “light cylinder” x, in the numerator and the denominator of Eqs. (32)–(34). This is induced by the regularity conditions, Eqs. (29) and (30), and the sphericity of the Alfvén surface. The surface x = 1 is the socalled outer “light cylinder”. Of course this surface may not be exactly cylindrical if x depends also on α, which may be the case for instance close to the black hole where ω has a strong dependence on α or if Ω is not constant with α. Therefore, this is rather a light surface, but for the sake of simplicity we refer to it as a “light cylinder” in the rest of the text.
Similarly, the pressure can be expanded to first order, (54)
where P_{0} is a constant.
In order to simplify and as a first step, we assume for both equations that the radial dependence of the nonpolar component of the Alfvén number and the pressure are simply constant, M_{1} (R) = m_{1} = cst, K(R) = κ = cst. Thus, (55)
We note that we have m_{1} = 0 in previous models; see Meliani et al. (2006) and Globus et al. (2014).
3.3 Choice of the free integrals
Free integrals are also expanded to the first order in the magnetic flux. The mass to magnetic flux ratio is similar to the one in the classical case, expanded as, (56)
where δ is a free parameter describing the deviations from spherical symmetry of the ratio number density/enthalpy as in Meliani et al. (2006) and not of the density itself, conversely to Sauty & Tsinganos (1994).
The total angular momentum loss flux density is given by (57)
Thus it is natural to expand the quantity LΨ_{A} rather than L itself. LΨ_{A} is also the poloidal current density along the polar axis and writes as (58)
The isorotation law can be expanded to first order as well as the total energy, (59)
where we see from Eq. (24) that .
Although we have some freedom with the choice of w_{1} and e_{1}, we could choose e_{1} = 0 and w_{1} = −δ∕2 to restrict ourselves to the values of the previous models, in particular in Schwarzschild metric; see Meliani et al. (2006) and Meliani et al. (2010). In fact, the isorotation function Ω does not need to be expanded beyond the zerothorder term because Ω always appears multiplied by another quantity as in (Ω − ω)ϖ or LΩ. Thus, the value of w_{1} is free and does not affect the solution. Conversely, the value of e_{1} affects the whole dynamics and we shall study the effects of its variation in a future publication. We already discussed the fact that taking a weak dependence of Ω on α has the advantage of minimizing the variation of the “light cylinder” near the base of the jet. Thus for the sake of simplicity, we shall study here the case where e_{1} = 0 and w_{1} = 0.
3.4 Constraints on the Alfvén Mach number, the isorotation law, and the angular momentum flux
The value of m_{1} is, in fact, determined by the prescription to cross the Alfvén transition surface. In order for the denominator in Eqs. (25)–(27) to vanish at the Alfvénic transition, the two following relations given in Eqs. (29) and (30) must be fulfilled. They can be expanded to first order. For the first regularity condition we get, (61)
The first term in the right part of Eq. (61) is due to the “light cylinder”, and the second one to the nonsphericity of the gravitational field in Kerr metric. m_{1} is negligible whenever the rotational speed λV_{⋆} is subrelativistic and either the μ parameter or the angular momentum of the black hole are negligible too. Since m_{1} < 0, there is a limiting field line where we have M_{Alf} = 0 since the magnetic flux increases going out from the polar axis.
To apply the second regularity condition we use the numerator of Eq. (25) and we get to the first order in α: (62)
The regularity conditions on the Alfvén surface fixes the value of m_{1}. Thus the critical Alfvén surface is a sphere, as in previous meridional selfsimilar models. We warn, however, that the Alfvén transition surface is a generalized or modified Alfvén surface as it takes into account the modification by the “light cylinder”.
Simultaneously, surfaces of constant Poloidal Alfvén Mach Number, M_{Alf} = const. (see Eq. (55)) are not spherical surfaces, conversely to the one defined by Meliani et al. (2006). Two effects modify it; first the “light cylinder” effect, which was neglected in Meliani et al. (2006) and Globus et al. (2014), and second the framedragging effect (LenseThirring).
3.5 Expansion of the velocity and magnetic fields
The model is obtained using an expansion to the second order for sin θ in the Euler equation. Due to axisymmetry, firstorder terms are zero along r and ϕ while the antisymmetry along θ gives the zeroth and second orders as null along the colatitude.
Then, for the poloidal velocity field, this gives (65)
And for the poloidal magnetic field, we get
Now from Eqs. (25) and (26), we can calculate to the first order in sin θ the toroidal components of fields
where the functions N_{V}, N_{B}, and D have been generalized,
3.6 Expansion of the enthalpy, densities, and electric field
We used Eq. (27) to deduce the enthalpy, (73)
and the mass density is given by (74)
In GRMHD, we also need the expressions of the electric field and the charge density. The electric field is a secondorder term for the radial component and a firstorder term for the θcomponent,
Using Maxwell–Gauss Eq. (7), we calculate the charge density from the divergence of the above electric field, to zeroth orderonly, (77)
With all these quantities we are able to expand the Euler equation. The radial component is expanded to the second order and the colatitude component to the first order. From the expansion of poloidal components in the Euler equation (Eq. (16)) and using Eq. (51), we can reverse the system to get the equations of the model (see Appendix C for details).
3.7 “Light cylinder”
The rescaling value ch∕(Ω − ω) for the cylindrical radius used in Eqs. (32)–(34) has been defined by Meliani et al. (2006) as the “light cylinder”. It is a surface of revolution Σ_{LC} where, (78)
On the “lightcylinder”, the electric field E is equal to the poloidal component of the magnetic field B_{p}. In the present publication, Σ_{LC} designed the external “light cylinder”, that is, x = +1, though, strictly speaking, this is not a cylinder as explained earlier, but a surface of revolution. This external “light cylinder” is outside the Alfvén surface since the denominator of Eqs. (32)–(34) is equal to on the “light cylinder”, is negative before crossing the Alfvén transition surface, and is positive after crossing it. At large distance in the jet, the lapse function h goes to unity and Ω − ω tends to Ω, which is assumed constant in our model. Thus, Σ_{LC} is located on a constant cylindrical radius along the z axis, becoming a real cylinder.
From the isorotation law, we get, (79)
As in special relativity, the second term of Eq. (79) cannot be neglected in the vicinity of the “light cylinder”. The sign of is such that always remains less than the speed of light (Vlahakis 2015). Moreover after crossing the “light cylinder” one of the two following conditions must be fulfilled. Either, we have » B_{p} or V_{p} » , or both.
The term x was neglected in the equation of the previous relativistic meridionalselfsimilar models, Meliani et al. (2006); Globus et al. (2014). Hence, these models could not produce jets crossing the “light cylinder”. Conversely, in this model this quantity is taken into account. We assume an expansion in sin(θ) of this quantity.
Contrary to the two previous models we can choose the dependence of the isorotation frequency with the magnetic flux (see discussion on Eq. (59)) and this choice will not affect the solution. Thus, if Ω does not depend strongly on the magnetic flux A, even at the base of the jet, the ratio h∕(Ω − ω) will be nearly constant. The reason for this is that ϖ is larger than the Alfvén radius, which is at least a few times the Schwarzschild radius. As a consequence, the departure of Σ_{LC} from a real cylinder is unnoticeable.
3.8 Domain of validity
The equations of the model are the result of an inversion of the expanded conservation equations. Therefore, it will be useful to quantify therelative error of the expansion we made in order to properly analyze our results and to obtain the domain of validity of these results. To have an idea of the domain of validity, we quantify the rest in the expansion of the momentum equation, for each force F^{i} (R, sin θ), (80)
where F is one of the following forces: gravitational, centrifugal, inertial, electric or magnetic pressure, and so on. We define a new function in order to map the relative error: (81)
For example, in the case of the electric force, in Schwarzschild metric, we get, assuming solid rotation (w_{1} = 0), (82)
The relative error on the electric force which tends to zero in the asymptotic regime of cylindrical jets is defined as (83)
Even at the base of the jet this error can be reduced, as can be seen in Fig. 1 for the solution in Kerr metric presented in Sect. 5.1 when the colatitude is less than 30^{°}. To get an estimate of the error in the expanded forces, we should add all relative error terms or take the largest one. This gives an estimate of the domain of validity of the solutions for a given set of parameters. We postpone the full error analysis for a future paper.
Fig. 1 Relative error on the electric force for a recollimating oscillating solution in Kerr metric (K1, see Sect. 5.1). Color isocontours correspond to the relative error in the electric force. Field lines anchored into the black hole magnetosphere and in the accretion disk are plotted in black solid lines. The limiting field line between the inner jet coming from black hole corona and the outflow outgoing from the accretion disk is plotted in red. The “light cylinder” is indicated by a green solid line. The cylindrical radius and the distance above the equatorial plane are in units of Schwarzschild radius. 
3.9 The magnetic collimation efficiency, ϵ
By writing the first law of thermodynamics in the frame of the fluid along streamlines of an axisymmetric flow, we can construct a constant of the motion, as in the classical case. The first law of thermodynamics reduces to the adiabatic law if the heating is included in some effective enthalpy (see Eq. (20)). Thus (ξc^{2}) is an effective specific enthalpy, like for polytropic flows where the enthalpy also hides the heating (cf. Sauty & Tsinganos 1994) but generalized for relativistic outflows (see Eq. (18)). Using Eq. (28), we can rewrite the first law of thermodynamics in the following form, (84)
As the magnetictomassflux ratio and the total energy flux are constant along each streamline, this is equivalent to (85)
We note that Ψ_{A}ξc^{2} is proportional to the thermal energy. If we write (86)
and using the expressions of the pressure and the Mach number, we get an equation of the form (87)
We see, as in the classical case, that the second term of the pressure is proportional to the first one such that, (88)
We deduce from the previous equation that the quantity ϵ, which is defined by (89)
is a dimensionless constant for all the field lines. To give explicitly and , it may be useful to write
Finally the calculation leads to (90)
This equation is similar to Eq. (71) in Meliani et al. (2006) and can be interpreted the same way. The parameter ϵ measures the efficiency of the magnetic rotator to collimate the flow. At the outflow base, ϵ is the relative difference of the transverse variation of internal energy that is simply the exchange of work done by the macroscopic forces. As this is perpendicular to the flow axis, this means that ϵ really measures the transverse force which collimates the flow and mainly its magnetic component.
We note that the quantity − 2m_{1} appears twice in Eq. (90). First, it is associated with κ, having a similar effect to the nonspherically symmetric pressure in the term which is given as a factor of M^{4}. Second, it is associated with 2e_{1} in the term corresponding to the excess or the deficit of the gravitational energy not compensated by the thermal driving at the base of the jet.
To conclude we can also derive the magnetic collimation efficiency in a different form; after some calculations, we can write (91)
This new relation brings a link between the total enthalpy on the axis and its logarithmic variation with α. In particular, the sign of ϵ seems to connect the balance between logarithmic variation of total enthalpy per unit of volume and the meridional increase of the pressure. The factor indicates that ϵ probably tends to decrease for solutions which reach ultrarelativistic speed.
4 Methodology for obtaining solutions
4.1 Numeral integration
Appendix C gives the coupled ordinary differential Eqs. (C.1), (C.5), (C.6) for the four quantities of the Alfvén number, dimensionless radius, expansion factor, and pressure, (M^{2}, G^{2}, F, Π). Details on the method and numerical techniques for the integration of these differential equations of the model system can be found in Sauty & Tsinganos (1994) and Meliani et al. (2006). Briefly, by using a Runge–Kutta scheme we start integrating from the Alfvén surface using the continuity relations. Integrating upwind by adjusting the value of the slope of the derivative of the Alfvén number at the Alfvén transition, the unique value of F_{⋆} can be found that allows the crossing of the modified slow magnetosonic surface, for a given value of the pressure at the Alfvén transition Π_{⋆}. Then, integrating downwind, the value of Π_{⋆} can be further adjusted. After several iterations of upstream and downstream integrations, the program automatically finds the solution that crosses all critical points, a proxy that the asymptotic pressure converges towards the demanded value. In particular, in the following Sections we have selected only solutions with the minimum possible value of Π_{⋆}, the socalled limiting solutions (see Sauty et al. 2004). We either have a collimated jet where, at infinity, Π_{∞} is minimum, or a conical wind where, at infinity, Π_{∞} is zero. If necessary, we can always add a constant value to the pressure P_{0} to ensure that the pressure is positive everywhere in the flow.
In the following, we only outline briefly the points which differ from previous related studies for building our model in the framework of general relativistic MHD in a Kerr metric. Some interesting spine jet solutions, both for their properties close to the black hole and at large distance, are presented in Sect. 5. The presented solutions depend on a number of parameters and a systematic parametric study of the model is postponed for a following paper.
4.2 Alfvén regularity conditions
The regularity conditions at the Alfvén transition (R = 1) for the azimuthal components and enthalpy have already been discussed in Sect. 2.5. In order for the field lines to not have a kink at the Alfvén transition (i.e., F is a continuous function across R = 1), we impose an extra regularity condition on the transfield equation giving the four physical quantities (M^{2}, G^{2}, F, Π) at R = 1 (Tsinganos & Trussoni 1991, Sauty & Tsinganos 1994). In other words, similarly to the classical model Sauty & Tsinganos (1994), we should take appropriately into account this regularity condition at R = 1, which gives the ratio τ which is involved in the magnetic toroidal component, (92)
where p is the slope of the square of the Alfvén number at the Alfvén transition, and (93)
More specifically, the expansion factor F is determined by Eqs. (C.1) and (C.3). Also in Eq. (C.3) are the ratios of N_{V} and N_{B} with D. Furthermore, D appears in the denominator with a power higher by one than the powers of N_{V} and N_{B} in the numerator. Hence, in order that dF∕dR does not diverge and the slope of F(R) is continuous across R = 1 we require that such that dF∕dR is finite at R = 1 wherein we have 0/0. Thus, near the Alfvén transition, we may define the function . Then, in order to avoid a singularity at the Alfvén transition, a necessary condition is to choose p, F_{⋆}, Π_{⋆} such that this function is zero, that is, , and thus . After some algebra we finally get a second degree polynomial for F_{⋆}, namely, (94)
In this way, the regularity condition at the Alfvén transition is automatically satisfied and no more constraints are needed at R = 1.
4.3 Effect of a nonspherical Alfvén number in a Schwarzschild metric
To illustrate our model, we present the following two solutions, which are built in the framework of the Schwarzschild metric. The first one corresponds to a solution of a model presented in Meliani et al. (2006), in which m_{1} = 0.
The chosen values of the other parameters are λ = 1.0, κ = 0.2, δ = 1.2, ν = 0.8, ℓ = 0, μ = 0.1, e_{1} = 0. We compare this solution to a solution with the same parameters but by keeping the value of m_{1} given by Eq. (61), m_{1} = −0.078. In both solutions, the value of Π_{⋆} is the minimum value of the limiting solution. Such solutions have the minimum amplitude of oscillations in the jet.
In Fig. 2 the radial velocity on the polar axis is compared for the two solutions, while field lines in the poloidal plane are plotted in Fig. 3. We note that in Eqs. (C.1) and (C.4) and Eqs. (C.1) and (C.3), giving the plasma acceleration and the variation of the expansion factor with the radius R, respectively, there are several terms proportionalto the factor (κ − 2m_{1}). As m_{1} is always negative, it is evident that − m_{1} effectively increases the transverse pressure gradient, which is proportional to κ. In other words, the effect of − m_{1} is similar to the effect of κ which enforces collimation for κ > 0. Hence, taking into account a nonsphericalAlfvén number (m_{1}≠0) introduces an extra collimation force which explains why the second solution with a nonzero m_{1} is more collimated. Indeed, the width of the jet at infinity over its value at the base, G_{∞} ∕G_{0}, decreases from 13.74 to 9.42, which is similar to an increase of κ and this fact can be checked directly by looking at the poloidal field line’s shape shown in Fig. 3.
Since more tightly collimated solutions have a smaller superAlfvénic acceleration, the second solution reaches a lower Lorentz factor asymptotically. Thus, similarly to κ, the introduction of a negative m_{1} ≠ 0 leads to a decrease of the velocity because of the tighter collimation. The higher collimation reduces the pressure gradient along the axis, which in turn decreases the acceleration due to pressure driving on large distances (see Sauty et al. 2004).
Additionally, m_{1} appears within the term (κ + 2m_{1} − δ − 2e_{1}) which appears in the plasma acceleration function in Eqs. (C.1) and (C.4) and the function determining the expansion factor F in Eqs. (C.1) and (C.3). The first three terms (κ + 2m_{1} − δ) arise from the variation across the field lines of the heat content P∕ρ with α, that is, , while the fourth term e_{1} is proportional to the variation of the total energy with α. The bigger this term is, the larger the initial acceleration (see Sauty & Tsinganos 1994), because it is linked with the distribution of the heating which opposes gravity to accelerate the outflow. As the weight of the plasma decreases with the latitude, then the pressure gradient increases along the axis resulting in a larger acceleration close to the base, as explained in Tsinganos & Sauty (1992). This term decreases rapidly as the Alfvén surface is reached. Thus, it is responsible only for the initial acceleration. With the parameter m_{1} being negative, the second solution is more accelerated between the base and the Alfvén surface. The velocity of the second solution reaches the velocity of the first one at the Alfvén surface. This effect disappears far from the source.
Fig. 2 Evolution of the radial velocity along the polar axis for solutions in a Schwarzschild metric, with m_{1} = 0 (blue) and m_{1} = −0.078 (red). The second case has a smaller terminal velocity. 
Fig. 3 Field lines for a solution in a Schwarzschild metric, with parameters λ = 1.0, κ = 0.2, δ = 1.2, ν = 0.8, ℓ = 0, μ = 0.1, e_{1} = 0 and m_{1}= 0 (left) and m_{1}= −0.078 (right). We note that the case m_{1} = −0.078 corresponds to a more tightly collimated jet. Lengths are in units of the Schwarzschild radius. The red lines are connected to the magnetosphere of the central object while the green lines are connected to the disk. The separating line is in blue and the light cylinder in black. 
Set of parameters used for the four selected solutions in the Kerr metric.
5 Solutions in a Kerr metric
In the following, we discuss four different solutions in a Kerr metric to illustrate the present model. A more detailed parametric study is postponed to a following paper. For the purposes of the present paper, we show three cylindrically collimated solutions with high asymptotic Lorentz factor, typical of AGNs and GRBs. Those solutions cross the “light cylinder” and are sorted with increasing magnetic collimation efficiency parameter ϵ. We also exhibit a conical solution crossing the “light cylinder” with high Lorentz factor and strongly negative ϵ, something that was not possible with the previous relativistic meridionally selfsimilar solutions.
In order to get a Lorentz factor as high as possible in the asymptotic part of the collimated part of the jet, we know from the study of the classical solutions that among all cylindrical solutions, the limiting solutions with the lowest value of Π_{⋆} reach the highest terminal velocity. These solutions are the so called limiting solutions in Sauty et al. (2002). As Π_{∞} is negative for the limiting solutions, we have to add a positive P_{0} value to thepressure. Of course, it is always possible for those cylindrical solutions to have a higher pressure P_{0}, but by doing so it also increases the effective temperature, in particular in the asymptotic part. For the same set of parameters, it is also possible to get cylindrical solutions by increasing Π_{⋆}. However, such solutions usually have a strong initial decollimation associated with a peak in both the Lorentz factor and temperature, while the asymptotic jet is decelerated to lower Lorentz factors and smaller radii, a result that we used to interpret the FRI/FRII dichotomy [cf. Meliani et al. (2010)].
Solutions K1 and K2 have been obtained for maximally rotating black holes, that is, a_{H} close to 1 (a ≃ r_{s}∕2). In solution K4 the value of a_{H} has been fixed to 0.73 (a = 0.73r_{s}∕2). We do not expect all black holes to be maximally rotating. For example, in M 87, the dimensionless spin should be above 0.65 (i.e., a > 0.65r_{s}∕2), (Li et al. 2009) but not too close to one. Other examples can be found and for K3 we will use the value a_{H} = 0.6 (a = 0.6r_{s}∕2) adopted in Mertens et al. (2016) for M 87.
Output parameters for the four solutions in the Kerr metric.
Fig. 4 Threedimensional (3D) representation of the field lines and streamlines for the thermally collimated solution K1 at the base of the jet and for two flux tubes. The blue lines correspond to streamlines, the red lines to magnetic field lines. The length is in units of the Alfvén radius, that is, ten times the Schwarzschild radius. 
5.1 A mildly relativistic collimated solution with oscillations (K1)
The collimated solution K1 corresponds to an overpressured outflow (κ ≥ 0) in a Kerr metric. As ϵ ≤ 0, the collimation of the jet is not fully magnetic but it has a significant contribution by the gas pressure, at least during the phase of strong acceleration up to ≃30r_{s}. In this solution, field lines are strongly oscillating compared to the two previous solutions in the Schwarzschild metric. The outflow undergoes a series of strong oscillations connected to the balance between the toroidal magnetic tension and the decollimation forces (centrifugal and electric forces and transverse pressure gradient) of the plasma.
The significant contribution of the transverse pressure gradient also explains these strong oscillations in the flow, as is also shown in the classical solutions. The parameters of this solution are displayed in the first row K1 of Table 1and the output values of m_{1} and ϵ in Table 2. We note that m_{1} = −0.062 is a relatively small value, which clearly indicates that the Alfvén surface is almost spherically symmetric in this case. As a consequence, the light cylinder is relatively far from the jet axis. Most of the central field lines (see the inner 5 to 7 central red lines in Fig. 5) remain within the “light cylinder”, which means that despite the important role of the magnetic field in the collimation, the jet is pressure or enthalpy driven in the relativistic case. However, unlike in the classical solutions, the electric field is the dominant decollimating force for the lines that cross the “light cylinder”. This decollimation and expansion after the Alfvén surface is associated with a strong pressure gradient yielding a strong acceleration of the jet in the superAlfvénic regime. More details on this will be given in the following solution. The pressure gradient is the gas pressure gradient close to the axis but assisted by the toroidal magnetic pressure outside the “light cylinder”. This is similar to superfast flows in radially selfsimilar models for disk winds (see Vlahakis & Königl 2003a and Vlahakis & Königl 2003b). Moreover, in the relativistic regime the inertia increases faster when the flow is accelerated such that the collimation from the magnetic field is delayed to larger distances.
In Fig. 5 we see that the expansion after the Alfvén surface is strong and leads to a late acceleration of the flow. After the large expansion, the jet recollimates smoothly and consequently decelerates slightly because of the compression.
In Fig. 4, we clearly see that there is a strong azimuthal magnetic field, although the scale of the figure tends to exaggerate this phenomenon.
The Lorentz factor γ of this solution reaches a relatively small value around 3.7, typical of less powerful AGN jets like some of the FRI radiogalaxies.
Fig. 5 Poloidal field lines and “light cylinder” for the thermally collimated solution K1, for λ = 1.0, κ = 0.2, δ = 2.3, ν = 0.9, μ = 0.1, ℓ = 0.05, e_{1} = 0. The length unit is the Schwarzschild radius. 
5.2 A highly relativistic collimated solution with oscillations (K2)
The solution K2 is collimated and has an extremely high Lorentz factor, which may be typical of GRBs. This K2 model corresponds to the values of the parameters given in the second line of Table 1, that is, λ = 1.0, κ = 0.2, δ = 1.35, ν = 0.46223, ℓ = 0.5, μ = 0.1, e_{1} = 0 and the second line of Table 2 for the output parameters, m_{1} = −0.234 and ϵ = −0.04.
For this model, the outflow starts very close to the black hole horizon at r_{0} = 1.57 r_{s}, (see Fig. 7), and thus at the base of the jet the effects of general relativity play an important role. The final velocity is highly relativistic, as is shown in Fig. 6.
The parameter ν of the K2 solution is accurately adjusted (to the fifth digit), so as to obtain a rather high Lorentz factor (larger than 100). This proves the versatility of the model which handles any magnitude of Lorentz factors. In order to obtain such high Lorentz factors, we must carefully tune the parameter directly linked to gravitation, ν, as mentioned above. The same parameter is also responsible for the thermal acceleration in the classical model (Sauty & Tsinganos 1994).
In Figs. 8 and 9 we plot the forces along and perpendicular to a field line defined by α =0.01α_{lim} where α_{lim} is the dimensionless magnetic flux between the inner jet and an external accretion disk wind.
The strong decollimation associated with the slow acceleration enhances the electric force as in the previous solution K1. This can be seen in Fig. 9. However due to the higher rotation here, more field lines cross the light cylinder, which is very close to the axis, such that the decollimation from the electric field is much stronger in this solution.
The feedback of this strong electric field is to further increase the decollimation beyond the Alfvén surface at large distances. Again, the large expansion increases the pressure and enthalpy gradient as seen in Fig. 8. Thus, the pressure force increases, resulting in a very long acceleration phase up to 10^{6} r_{s}. Thus, the thermal acceleration becomes very efficient and the plasma asymptotically reaches an extremely high Lorentz factor. Indeed, the jet radius increase of this solution is also very high with an expansion factor G_{∞} ∕G_{0} ≃ 1300, as can be seen in Fig. 7.
The value of ϵ is still negative but very close to zero (ϵ = −0.04). As in the classical case, this means the magnetic efficiency to collimate the flow is higher in this model at larges distances. However the decollimating force that ensures the equilibrium is no longer the centrifugal force or the pressure gradient but the electric force on the lines that cross the “light cylinder”. This is a specific feature of relativistic jets.
To further analyze the jet acceleration we may calculate the contribution of the different components of the total energy and the conversion of the magnitude of each component to another form along the streamlines. First, we want to obtain a physically acceptable heating term which goes to zero at infinity along the axis of the flow. By defining the enthalpy in Eqs. (18) and (19) and its analytical expression in our model in Eq. (34), we may fix a streamline limiting the zone where the pressure is positive. Once this field line has been chosen, for example, at α =0.4α_{lim}, the dimensionless pressure P_{0} ∕ρ_{⋆}c^{2} can be calculated. In the particular case presented here, we have taken P_{0}∕ρ_{⋆}c^{2} = 4.3 × 10^{−7}. The pressure P_{0} in Eq. (54) is chosen such that the gas pressure is equal to zero, when Π (R) reaches its minimum value. For some value of α, the term will be too large to ensure that Q goes to zero at infinity along the axis. Indeed, α can take any value below a maximum α_{max}. We can consider that the solution can be valid only in the region wherein the pressure is positive. When α is fixed and P_{0} ∕ρ_{⋆}c^{2} is deduced, we are able to calculate ξ_{⋆}. In this particular case, a value ξ_{⋆} ≃ 78 is chosen.
Figure 10 shows the normalized total energy on the axis , the kinetic component hγξ_{K} and the external heating hγQ∕c^{2}. We find a decrease of the external heating and a related increase of the kinetic part. Thus, the kinetic enthalpy represents the major component up to r = 10^{2}r_{s}.
Figure 11 shows the same energetic distribution, but on a streamline with α = 0.05α_{lim}. Out of the polar axis, there are extra energetic components, such as the framedragging and the Poynting fluxes. Both these energetic contributions are very small on this field line, as compared to the total energy. Hence, the jet is enthalpydriven from the axis right up to the limiting line. We also note that at the base of the jet, the framedragging energy is of the same order as the Poynting flux. Contrary to the energetic distribution along the axis, the external heating constitutes the larger part of energy at infinity. While along the axis a high value of γ_{∞} ≃ 100 is obtained, the Lorentz factor at infinity on this particular line is γ_{∞}≃ 3.6. Hence, since the acceleration of the plasma and the resulting final flow speed at infinity on this particular field line are small, the external heating is not consumed when accelerating the flow and therefore it is left unused at infinity, contrary to what happens along the axis; Fig. 10.
Fig. 6 Lorentz factor γ for the K1 (blue line) and K2 (red line) solutions. Distances are given in Schwarzschild radius units. 
Fig. 7 Poloidal field lines and “light cylinder” for the K2 solution, that is, for λ = 1.0, κ = 0.2, δ = 1.35, ν = 0.46223, μ = 0.1, ℓ = 0.05, e_{1} = 0. Distances are given in Schwarzschild radius units. 
Fig. 8 Longitudinal forces, that is, along the field line, for the K2 solution, along the line α = 0.01α_{lim}. Distances are given in Schwarzschild radius units. 
Fig. 9 Transverse forces, that is, perpendicular to the field line, for the K2 solution, along the line α = 0.01α_{lim}. We see that the Lorentz force is collimating and is balanced by the electric force that decollimates. Distances are given in Schwarzschild radius units. 
Fig. 10 Relative normalized contribution to the total energy of the kinetic enthalpy hγξ_{k} and the external heating distribution hγQ∕c^{2} along the axis. 
Fig. 11 Relative normalized contribution to the total energy of the kinetic enthalpy hγξ_{k} and the external heating distribution hγQ∕c^{2} along the streamline α = 0.05α_{lim} 
5.3 A mildly relativistic collimated solution without oscillations (K3)
From VLBI imaging, Mertens et al. (2016) recovered a detailed twodimensional velocity field in the jet of M 87 at subparsec scales. They confirmed the stratification of the flow from the very beginning of the jet and identified a relativistic sheath, that is, an accelerating layer, which is launched from the inner part of the accretion disk at a cylindrical distance of around 5r_{s}. Mertens et al. (2016) interpret this outer sheath layer as the internal part of an external disk wind. They also interpret the inner spine jet as a component coming from the internal accretion disk. However, it is not clear that the inner spine necessarily has to be one of the disk wind components. Instead we propose that the inner spine jet originates from the black hole corona and has a higher Lorentz factor. The authors note that this fast inner spine jet cannot be detected in their data because of its lower emissivity compared to that of the sheath layer or because its speed is too high.
We propose here, as an alternative scenario, that the spine beam may originate from the magnetosphere of the black hole, either connected to the black hole itself, as in the Blandford–Znajek–Penrose mechanism, or connected to the inner part of the accretion disk. In the first case, the spine jet would be a leptonic plasma, and in the second it would be a hadronic one. In such a case, we can model the spine jet with our meridionally selfsimilar solutions. We note that unlike Poyntingfluxdominated models, our model is valid on the jet axis.
The Lorentz factor profile inferred by Mertens et al. (2016, see their Fig. 19) supposes that the velocity structure observed in the 43 GHz VLBA maps at a deprojected distance of z ≃ 10 mas ≃ 1300r_{s} is due to the sheath layer. The Lorentz factor of the sheath has a value γ ≃ 2.4 at R ≃ 1000r_{s}.
From the curve deduced for the spine jet with an approximated MHD disk wind solution, those authors, following Anderson et al. (2003), model the acceleration and collimation of the flow with a Lorentz factor γ ≈ 7. We propose that it is possible to construct a MHD solution along the lines of the model analysed in this paper, by choosing a positive value of ϵ ≈ 0.5 and a similar velocity profile. The solution K3 presented here has a Lorentz factor γ ≃ 5 at R ≃ 1000r_{s}. By appropriately tuning the parameters, we could obtain even higher Lorentz factors from 7 to 10, as is observed at the distance of the knot HST1, where this velocity is observed in the optical band. However, the spine jet is deboosted relative to the sheath and a precise measurement of its Lorentz factor is difficult.
Solution K3 with a maximum Lorentz factor on the axis γ ≃ 5.5 has an asymptotic spine jet radius G_{∞}≃ 20r_{s}. At this distance from the axis, the Lorentz factor has dropped to a value γ ≃ 2.4 consistent with the radius and the Lorentz factor at the inner observed distance of the outer sheath jet of Mertens et al. (2016). Thus, our solution may model the initial spine jet inside the sheath layer. However, to confirm that, we need to use this initial solution in simulations similar to those in Hervet et al. (2017). The spine jet/sheath jet interaction will probably produce shocks and rarefaction waves that may further accelerate the jet.
Therefore, although we can obtain such a type of collimated solution for different sets of the parameters, we focus here on the specific solution K3, where λ = 1.2, κ = 0.005, δ = 2.3, ν = 0.409, ℓ = 0.024 and μ = 0.08. Compared to the other Kerr solutions studied in this paper, λ is higher and κ is very small, leading to a solution with a positive value of the magnetic collimation efficiency, ϵ = 0.55. Another advantage of this solution is the fact that the pressure depends only very weakly on the magnetic flux function, that is, on a particular field line.
Interestingly, in Mertens et al. (2016), the radius of the base of this spine jet is equal to r_{0} = 2.4r_{s}, which is consistent with their model of a disk wind solution, by fixing the jet shape and solving the Bernoulli equation. In our selfsimilar solution, we have a similar radius for the magnetospheric polar cup, where our jet solution starts. Clearly, this is an alternative scenario.
Moreover, the angular velocity of the field lines anchored in this polar cup above the black hole can be calculated from our parameters, (98)
By taking the value of the M 87 distance and the black hole mass, as in Mertens et al. (2016), we can calculate Ω_{⋆} in the context of this K3 solution. This value may directly be compared to the values they deduced in two jet regions from the conservationof total energy and angular momentum fluxes in the approximation of special relativity.
Hence, we find Ω_{⋆} ≃ 1.03 × 10^{−6}s^{−1}, a value which is for the spine jet almost the same value as the isorotation frequency of a Keplerian speed at the launching location of the sheath layer. It also corresponds to the initial toroidal velocity of the Blandford & Payne (1982) mechanism.
Our model corresponds to an alternative configuration, because the spine jet may either originate from the Keplerian disk, in which case it would be hadronic, or form via the generalized Penrose–Blandford–Znajek mechanism. In this second alternative, the jet would be a leptonic beam with an angular frequency proportional to the spin of the central black hole. The Blandford–Znajek mechanism allows us to extract energy from the black hole when 0 ≤ Ω ≤ ω_{BH}, with a maximum value for 0.5ω_{BH}. We note that ω_{BH} is, by definition, the angular velocity of ZAMO at the location of the outer event horizon and is given by, (99)
where a_{H} is the dimensionless spin of the black hole in units of the gravitational radius r_{s} ∕2. Indeed simulations of such Poyntingdominated and forcefree jets (Tchekhovskoy 2015) have shown that the angular speed of a field line anchored in the magnetosphere is about half the black hole angular speed ω_{BH}. Our value of Ω_{⋆} is one third of 0.5ω_{BH} ((Nathanail & Contopoulos 2014)). In order to determine if the spine jet of our solutions originates from a Keplerian disk or froma black hole via a generalized Penrose–Blandford–Znajek mechanism, we need to solve the MHD equations up stream up to the black hole horizon. Moreover, in order to model the full jet of M 87, a complete MHD simulation including a disk wind and a spine jet must be developed. Something that should be done in the future.
We notethat as is already known, at the interface of the spine jet and the sheath layer, a recollimation shock may occur, producing compression and rarefaction waves which may accelerate the flow (Hervet et al. 2017). For those reasons we have chosen here to adjust the value of ν to ν = 0.42, in order to get the solution K3, with a lowerLorentz factor but suppressing completely the oscillations of the field lines. The value calculated above for Ω_{⋆} is not changed, but radio emission maps, which are obtained for the M 87 jet, will be produced by recollimation shocks due to the interaction between the fast spine jet and the sheath layer.
As can be seen in Fig. 12 for the K3 model, the Lorentz factor reaches a nearly constant value at the distance of the B structure observed in the 43 GHz VLBA maps and γ is larger than that deduced for the sheath layer by Mertens et al. (2016). As explained above, the interaction between the sheath layerand the spine jet can induce a bulk flow acceleration. However, as we mention, the Lorentz factor γ is maximum along the axis and decreases with latitude such that at its outer boundary it matches the sheath layer value.
The values of the parameters of the K3 model are given in the third line of Tables 1 and 2. The radius at the base of this spine jet is equal to r_{0} = 2.55r_{s} and the magnetic efficiency to collimate the flow is larger compared to the other two Kerr solutions discussed in this paper, as ϵ is now positive, which means that the jet is fully magnetically collimated. The field lines displayed for K3 in Fig. 13 are nearly cylindrical above the equatorial plane, at distances 10^{4} r_{s}, with a smooth flaring occurring after the Alfvén surface.
Fig. 12 Lorentz factor for the nonoscillating collimated solution K3. 
Fig. 13 Poloidal field lines and “light cylinder” for the nonoscillating collimated solution K3, that is, for λ = 1.2, κ = 0.005, δ = 2.3, ν = 0.42, μ = 0.08, ℓ = 0.024, e_{1} = 0. Distances are given in Schwarzschild radius units. 
Fig. 14 Poloidal field lines and “light cylinder” for the conical solution K4, that is, for λ = 0.0143, κ = 1.451, δ = 3.14, ν = 0.8, μ = 0.41, ℓ = 0.15, e_{1} = 0. Distances are given in Schwarzschild radius units. 
5.4 Conical solution
The parameters of the conical solution K4 are given in Tables 1 and 2. Such radial solutions could be useful to describe relativistic noncollimated outflows, such as those seen in association with radioquiet galaxies, such as Seyfert galaxies. Conical solutions could also be useful to model GRBs, wherein an unstable noncollimated relativistic wind may fragment into small pieces under some instabilities. In such cases, the apparent collimation of the GRB would be due to fragmentation – see Meliani & Keppens (2010); van Eerten et al. (2011).
A conical solution is a solution in which the limiting value of F at infinity is 0. In this solution, the spherical part of the Mach number diverges, M →∞. The same effect occurs for the cylindrical radius of the flow, G. To obtain this conical solution, we started with the parameters used for modeling the solar wind, as in Sauty et al. (2005). The magnetic collimation parameter must be strongly negative, ϵ∕(2λ^{2}) ≤ 0. We adjust the solution to the relativistic case and increase the velocity by taking a larger value for δ. For the radial solution K4, the parameters λ and κ are adjustedas well in order to obtain a terminal Lorentz factor larger than 8. Thus, we find a conical solution K4 for the following set of parameters: λ = 0.0143, κ = 1.451, δ = 3.14, ν = 0.8, ℓ = 0.15, μ = 0.41, m_{1} = −0.004 and e_{1} = 0.
For thisset of parameters we get a much more negative value for the magnetic collimation efficiency parameter, ϵ = −5.68. The solution quickly reaches the conical regime (see Fig. 14), and most of the field lines cross the “light cylinder” plotted as a black solid line. The axial radial velocity profile is plotted in Fig. 15. As can be seen, high Lorentz factors are obtained.
5.5 Magnetic collimation efficiency versus black hole spin
As already mentioned, the constant ϵ is a measure of the efficiency of the magnetorotational forces to collimate the flow. We have studied how ϵ relates to the black hole spin a_{H} = 2l∕μ, keeping unchanged all the other parameters: λ, κ, δ, ν, and μ.
In Fig. 16, we plot ϵ versus the black hole spin a_{H}, for several cases studied in the context of our model. First we note that for solutions with parameters similar to the K3 solution, referred to here as K3type solutions, the value of ϵ is positive, it increases with the black hole spin a_{H} and shows the largest variation in relative magnitude with an absolute total variation equal to 0.05. On the other hand, for K1type and K4type solutions, the values of epsilon are negative. The total variation for K1type solutions is equal to 0.10 and for the K4type solutions is equal to 0.25. We see that ϵ is increasing with the black hole spin a_{H} for K1type and K3type solutions. It also presents a minimum value for K4 at a black hole spin slightly smaller than a_{H} = 0. This means that the magnetic collimation efficiency is lower for counterrotating black holes in relation to their accretion disk for K1type and K3type solutions and depends weakly on the spin direction for K4type ones. However, K4type solutions are conical, contrary to all other solutions, which are cylindrically collimated. In all cases ϵ does not vary linearly with the black hole spin a_{H}.
This nonlinear variation of epsilon with the black hole spin a_{H} is explainedif we try to derive it at the base of the jet where the Alfvénic number is equal to 0, since some terms of the second order in l cannot be neglected in the following equation. (100)
An increase in ϵ coincides with a decrease in the maximal Lorentz factor for the collimating Kerr solutions, as we will see for K1type and K3type solutions. This explains why it was not possible to obtain physical solutions by decreasing l for K2type solutions. For this second Kerr solution, we performed a fine tuning of the parameters to obtain the largest possible Lorentz factor at large distances from the maximally rotating black hole. Then, decreasing l leads automatically to exceed the value of the speed of light, c, for the polar velocity at some distance in the jet. The acceleration phase does not vary for K1type and K3type solutions with l except immediately before R = 500 where the Lorentz factor reaches a plateau. The value of the maximum γ increases when l decreases.
For the two collimated solutions K1 and K3, there is a clear effect of collimation induced by the rotation of the black hole. In Fig. 17 we plot for five values of the black hole spin the evolution of the factor G along the jet, that is„ the ratio of the jet cylindrical radius divided by its value at the Alfvén surface. From this plot it can be seen that the maximum of the radius of the jet is reached at different distances, as the spin varies, and decreases when a_{H} increases. The same trend is observed for the terminal jet radius but the ratio of G_{∞} /G_{0} decreases only by a factor of 0.96 between a nonrotating and a maximally rotating black hole (a_{H} = 0.99). Hence, the faster the black hole rotates, the smaller the maximal jet radius. This result is expected because for a fast black hole rotation, the magnetic collimation efficiency parameter (ϵ) is higher.
The case of K3type solutions is simple, as the factor G increases with distance until it reaches a constant value. The ratio G_{∞} /G_{0} gives directly the expansion factor which is decreasing when a_{H} increases from − 0.99 up to 0.99.
As can be seen from Fig. 18, the Lorentz factor maximum follows the opposite trend to ϵ for the conical K4type solutions: the minimum value of ϵ is obtained for a_{H} ≃ 0 and for a nonrotating black hole the Lorentz factor curve reaches a maximum before decreasing up to a plateau at large distances. This type of curve is no longer observed when the absolute value of a_{H} goes above some threshold. The increase of the plateau value for the Lorentz factor is much more pronounced for a_{H} > 0 but can be seen also for negative values of the spin. At a distance of r = 1000r_{s} γ increases from a value of the order of 6 for a_{H} = 0 up to 14 for a_{H} = 0.99.
For the conical solutions wherein the asymptotic geometry is fixed, their geometry is less affected. The effect of varying ϵ is different in this case, as compared to this effect for the cylindrically collimated solutions. To analyze the effect of collimation induced by the black hole spin, we plot the ratio of the cylindrical to spherical radius of the field lines for K4type solutions and five different values of the black hole spin. We note that the ratio of the cylindrical to spherical radius of the field lines gives their opening angle with respect to the axis. The field lines start at the base of the jet with an opening angle which increases with the radial distance and expands away from the polar axis more rapidly for higher black hole rotation. The effect is more pronounced for negative spin parameters. At some distances above the equatorial plane (z ≃ 100r_{s}) the opening angle becomes constant as the field lines are finally becoming radial. It is found that for a given ratio of α∕α_{lim}, where α_{lim} is the last open field line, the asymptotic opening angle is constant for all K4type solutions, regardless of the spin of the black hole.
As seen in Fig. 19, the collimation of the jet increases with ∣ a_{H} ∣ which is consistent with the increase of ϵ. We note also that the base of the jet becomes slightly closer to the black hole horizon as ∣ a_{H} ∣ increases.
Globally we see that the total geometry of the solution, and in particular the expansion of the jet radius, is extremely sensitive to the black hole speed.
Fig. 15 Lorentz factor for the conical solution K4. Distances are given in Schwarzschild radius units. 
Fig. 16 Variation of the magnetic collimation efficiency parameter ϵ vs. the black hole spin parameter a_{H} for K1type, K3type, and K4type solutions. The value of ϵ has been normalized by ϵ(0) when it is negative and by ϵ(0) when it is positive. 
Fig. 17 Cylindrical jet radius normalized to its value at the Alfvén surface, G, for K1type solutions, as a function of the distance along the polar axis, for five different values of the black hole spin a_{H}. The functionG is equal to 1 at the Alfvén distance r = 10r_{s}. 
Fig. 18 Lorentz factor γ for K4type solutions when the black hole spin parameter a_{H} varies between − 0.99 and 0.99. 
Fig. 19 Ratio of the cylindrical radius to the spherical radius for the conical K4type solutions with a spin parameter a_{H} varying between − 0.99 and 0.99 vs. the distance z above the equatorial plane in units of the Schwarzschild radius. 
5.6 The magnetic flux and power of the jets
Intuitively, by physical arguments of magnetic flux conservation, it is expected that magnetic fields not only play a dominant role in collimating largescale AGN jets, but also critically affect the origin of the jets in accretion disks of black hole systems, which are accordingly termed magnetically arrested disks (Narayan et al. 2003). Indeed, theoretical modeling concludes that magnetic fields at the base of AGN jets are related to the corresponding accretion rate (Tchekhovskoy & McKinney 2012). Zamaninasab et al. (2014) reported that the measured magnetic flux of the jet and the accretion disk luminosity are tightly correlated over several orders of magnitude for a sample of many radioloud AGN, concluding thus that the jetlaunching region is threaded by a dynamically important magnetic field. The magnetic fields of AGN can be measured either by the effect of a frequencydependent shift of the VLBI core position (known as the coreshift effect), or by Faraday rotation (e.g., MartíVidal et al. 2015, who reported magnetic fields of at least tens of Gauss on scales of the order of several light days – 0.01 parsecs – from the black hole). Furthermore, MHD simulations in the frame of general relativity allow us to calculate the saturation or equilibrium value for the poloidal magnetic flux Φ_{BH} threading the black hole (McKinney 2005, Tchekhovskoy et al. 2011, McKinney et al. 2012). In Zamaninasab et al. (2014) this poloidal flux Φ_{BH} is calculated as a function of the massaccretion rate Ṁ as, (101)
We consider that since such a strong magnetic flux can thread the black hole, we can use this formula to link our value of poloidal magnetic field at the Alfvén point B_{⋆} with the massaccretion rate Ṁ. In our model the magnetic flux from each hemisphere is given by . Then, the magnitude of the magnetic field is calculated from the expression (102)
In orderto compare the jet power for the K2 and K3 solutions with the one obtained by general relativistic MHD simulations, we calculate this power in terms of the parameters of our model. Similarly to the way we deduced the angular momentum flux density in Eqs. (57) and (58), we may calculate the jet+counterjet power by substituting Ψ_{A} from Eq. (56) and with from Eq. (60), and with the help of Eq. (40), we obtain, in terms of the constants of our model, (103)
Therefore, for the efficiency η_{jet} ≡ P_{jet}∕Ṁc^{2} for our K2 model we obtain a value η_{jet} ≃ 0.52, while for our K3 model we obtain η_{jet} ≃ 0.40. On the other hand, McKinney (2005) determined selfconsistently the jet power in Blandford–Znajek numerical models and deduced an efficiency η_{jet} between 0.01 and 0.1 for ultrarelativistic Poyntingdominated jets with a_{H} larger than 0.8. Later, Tchekhovskoy et al. (2011) and McKinney et al. (2012) increased the magnetic flux which can be pushed near the black hole leading to magnetically arrested accretion and obtained values of the net flow efficiency larger than 1 for rapidly spinning black holes with a_{H} larger than 0.9. Their models that develop a highly nonaxisymmetric magnetically choked accretion flow, initially have the poloidal component of the magnetic field dominant and the wind has an efficiency always smaller than that of the jet. We note that the net flow efficiency for the jet is equal to where Ṁ_{H} is the black hole massaccretion rate and is the timeaveraged value of accretion power. This black hole massaccretion rate could be smaller than the massaccretion rate measured by Zamaninasab et al. (2014). In fact, they deduced the accretion power by dividing the bolometric luminosity with a radiative efficiency of 0.4. Larger values of the inflow rates Ṁ_{in, i} and Ṁ_{in, o} have been obtained by McKinney et al. (2012) at radii 5r_{s} and 25r_{s}, respectively.Zdziarski et al. (2015) found also that the jet power moderately exceeds the accretion power Ṁ c^{2} for blazars estimating the magnetic flux from the radiojet coreshift effect and the selfabsorbed flux evaluation. However, there is a large scatter around the mean value for blazars η_{jet} ≃ 1.3 and the jet power for radio galaxies is smaller, especially for M 87. Therefore, our estimations of jet power from our K2 and K3 solutions with mass loading could perfectly match that for lessefficient blazars and radio galaxies. Moreover, they do not overly depend on the spin parameter a_{H}, since α_{lim} keeps a value slightly smaller than 1 when the spin parameter varies for the K3type solutions, even for retrograde black holes.
At this point, we prefer to postpone a further discussion of the jet power, until we have completed our study, which also includes inflow solutions and leads to a spinenergy extraction or addition from the black hole (work in progress).
6 Summary and conclusions
As was pointed already in 1957 by Parker (see also Parker 1963), for the driving of the solar wind and similar enthalpydriven astrophysical outflows, some energy/momentum addition is required. The original isothermal and polytropic models with a heat conduction, have shown that effectively energy and/or momentum are necessary for producing supersonic/superAlfvenic outflows at large distances, to also meet the respective causality requirements. Quasiradial windtype astrophysical outflows with shock transitions (Habal & Tsinganos 1983) have been applied to explain the appearance of emission knots in galactic (Silvestro et al. 1987) and extragalactic objects (Ferrari et al. 1984, Ferrari et al. 1986), also in the framework of special relativity (Ferrari et al. 1985).
When deviations of the outflow geometry from radial expansion exist and the problem is fully twodimensional, these outflows can be collimated in the form of jets (Sauty & Tsinganos 1994) mainly by magnetic fields with a suitable external gas pressure distribution. Along these lines, in Vlahakis & Tsinganos (1998), the original Parker model was extended to include general MHD effects, in the context of meridional selfsimilarity. The present paper takes the extra step of using the framework of a Kerr metric to explore analogous enthalpy or generalized pressuredriven outflows from the environment of a rotating black hole.
Specifically, in this paper we present an exact MHD solution for an outflow in a Kerr metric, constructed by using the assumption of self similarity and the mechanism for driving the outflow which is developed in Sauty & Tsinganos (1994). Additionally, the model is based on a firstorder expansion of the governing general relativistic equations in the magnetic flux function around the symmetry axis of the system. It yields four nonlinear and coupled differential equations as a function of the radius, for the Alfvén number, the gas pressure, the expansion function and the radius of the jet. The model depends on seven parameters. Two of them are the meridional increase of the gas pressure and the mass to magnetic flux ratio, κ and δ, respectively.There is also the meridional increase of the total energy with the magnetic flux function, e_{1}, the poloidal current density flowing along the system axis, λ, the escape speed in units of the Alfvén speed, ν, the Schwarzschild radius in units of the Alfvén radius, μ, and the dimensionless black hole spin l. In addition to those seven parameters, we have to adjust the pressure at the Alfvénic transition. We chose toadjust it such as to minimize the oscillations of the magnitude of the flow speed along the axis and taking the limiting solution. We also fix the magnetic field at the Alfvén transition, B_{⋆}, and a uniform pressure constant P_{0} to ensure a zero external heating at infinity along the axis where the Lorentz factor is maximum.
The model takes into account the light cylinder effects and the meridional increase of the Alfvén number with the magnetic flux function, m_{1}. This parameter is deduced from the regularity conditions at the Alfvén transition surface.
The classical energetic criterion for the transition from conical winds to cylindrical jets is generalized in general relativity and it amounts to say that if the total available energy along a nonpolar streamline exceeds the corresponding energy along the axis, then the outflow collimates in a jet.
In the framework of a Kerr metric, we illustrate the model with four different enthalpydriven solutions wherein the contribution of the Poynting flux is rather small. The first three solutions are cylindrically collimated, while the fourth represents a conical outflow at infinity. The flow collimation is induced by electromagnetic forces. In all four models, relativistic speeds are obtained, while in one of them the Lorentz factor γ obtains ultrarelativistic values. A preliminary application of one of our Kerr solutions (K3) was explored to model the spine jet in M 87, yielding encouraging results. A more complete modeling for the M 87 jet including an external diskwind component will be explored in another study.
Our analytical solutions of the full general relativistic MHD equations in a Kerr metric may contribute to a better understanding of relativistic AGN jets and are complementary to sophisticated numerical simulations of such jets (e.g., McKinney 2005, Tchekhovskoy et al. 2011, McKinney et al. 2012). In both the analytical and the numerical approach, the outflows are electromagnetically confined. However, while in the above numerical simulations the outflow is driven electromagnetically (e.g., via the Blandford–Znajek mechanism), in the present analytical solutions the outflow from the hot corona surrounding the black hole is enthalpy or generalized pressuredriven (e.g., via the Sauty–Tsinganos mechanism). Nevertheless, it is interesting to note that the jet powers for the two representative analytical solutions we present in this paper are similar to those determined by the numerical simulations.
The present model can also serve to construct an inflow solution, in order to link it with an outflow solution and the physical creation of leptonic pairs to determine the energy balance of the black hole, via a generalized Penrose process as compared to the Blandford–Znajek mechanism. This undertaking is in progress and will be presented in another publication.
Acknowledgements
Part of this work was supported by the French Action Fédératrice CTA project and the “Programme Blanc” of Paris Observatory. The authors thank the referee for his constructive suggestions which lead to a better presentation of the results of the paper, in particular for a comparison with the results of numerical simulations.
Appendix A Kerr metric elements
We can rewrite the Kerrmetric in a simpler form: (A.1)
The Kerr metric elements are defined as
Thus, to second order in sin θ the metric is
In orderto simplify our notation, we define along the polar axis, the axial lapse function, (A.11)
and the axial shift of the metric, (A.12)
We can also write the metric as an expansion in α:
Appendix B Vectorial operators in Boyer–Lindquist coordinates of the Kerr metric.
In this Appendix we summarize the expressions of the vectorial operators under the assumption of axisymmetry (∂_{ϕ} = 0).
1. The gradient vector ∇ in the ZAMO orthonormal bases, (B.1)
2. The divergence of a vector V, (B.2)
3. The scalar Laplace operator,
4. The curvature operator on a vector V, (B.3)
After some algebra we get for the poloidal component of the advection term: (B.4)
It can be useful to get the nonsymmetric advection term, (B.5)
Appendix C Final differential equations of the MHD problem
The final ordinary differential equations of our model can be written as: (C.1)
where X_{+} = R^{2} + ℓ^{2} X_{−} = R^{2} − ℓ^{2}. (C.4) (C.5) (C.6)
References
 Anderson, J. M., Li, Z.Y., Krasnopolsky, R., & Blandford, R. D. 2003, ApJ, 590, L107 [NASA ADS] [CrossRef] [Google Scholar]
 Begelman, M. C., Fabian, A. C., & Rees, M. J. 2008, MNRAS, 384, L19 [NASA ADS] [CrossRef] [Google Scholar]
 Beskin, V. S. 2010, MHD Flows in Compact Astrophysical Objects: Accretion, Winds and Jets (Berlin, Heidelberg: SpringerVerlag) [CrossRef] [Google Scholar]
 Blandford, R. D., & Znajek, R. L. 1977, MNRAS, 179, 433 [NASA ADS] [CrossRef] [Google Scholar]
 Blandford, R. D., & Payne, D. G. 1982, MNRAS, 199, 883 [NASA ADS] [CrossRef] [Google Scholar]
 Böttcher, M., Reimer, A., Sweeney, K., & Prakash, A. 2013, ApJ, 768, 54 [NASA ADS] [CrossRef] [Google Scholar]
 Breitmoser, E., & Camenzind, M. 2000, A&A, 361, 207 [NASA ADS] [Google Scholar]
 Cayatte, V., Vlahakis, N., Matsakos, T., et al. 2014, ApJ, 788, L19 [NASA ADS] [CrossRef] [Google Scholar]
 Casadio, C., et al. 2015, ApJ, 813, 51 [NASA ADS] [CrossRef] [Google Scholar]
 Celotti, A., & Ghisellini, G. 2008, MNRAS, 385, 283 [NASA ADS] [CrossRef] [Google Scholar]
 Fabian, A. C., & Rees, M. J. 1995, MNRAS, 277, L55 [NASA ADS] [Google Scholar]
 Feng, J., Wu, Q., & Lu, R.S. 2016, ApJ, 830, 6 [NASA ADS] [CrossRef] [Google Scholar]
 Ferrari, A., Habbal, S. R., Rosner, R., & Tsinganos, K. 1984, ApJ, 277, L35 [NASA ADS] [CrossRef] [Google Scholar]
 Ferrari, A., Trussoni, E., Rosner, R., & Tsinganos, K. 1985, ApJ, 294, 397 [NASA ADS] [CrossRef] [Google Scholar]
 Ferrari, A., Trussoni, E., Rosner, R., & Tsinganos, K. 1986, ApJ, 300, 577 [NASA ADS] [CrossRef] [Google Scholar]
 Finke, J. D. 2016, ApJ, 830, 94 [Google Scholar]
 Giroletti, M., et al. 2004, ApJ, 600, 127 [NASA ADS] [CrossRef] [Google Scholar]
 Ghisellini, G., Tavecchio, F., & Chiaberge, M. 2005, A&A, 432, 401 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Globus, N., Sauty, C., Cayatte, V., & Celnikier, L.M. 2014, Phys. Rev. D, 89, 124015 [NASA ADS] [CrossRef] [Google Scholar]
 Gracia, J., Vlahakis, N., Agudo, I., Tsinganos, K., & Bogovalov, S. V. 2009, ApJ, 695, 503 [NASA ADS] [CrossRef] [Google Scholar]
 Hada, K., Doi,A., Kino, M., et al. 2011, Nature, 477, 185 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Habal, S. R., & Tsinganos, K., JGR, 88, 1965 [Google Scholar]
 Hervet, O., Meliani, Z., Zech, A., et al. 2017, A&A, 606, A103 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Komissarov, S. S. 2007, MNRAS, 382, 995 [NASA ADS] [CrossRef] [Google Scholar]
 Kravchenko, E. V., Kovalev, Y. Y., Hovatta, T., & Ramakrishnan, V. 2016, MNRAS, 462, 2747 [NASA ADS] [CrossRef] [Google Scholar]
 Li, Y.R., Yuan, Y.F., Wang, J.M., Wang, J.C., & Zhang, S. 2009, ApJ, 699, 513 [NASA ADS] [CrossRef] [Google Scholar]
 Marscher, A. P., et al. 2010, ApJL, 710, 126 [Google Scholar]
 Martí, J. M., Perucho, M., & Gómez, J. L. 2016, ApJ, 831, 163 [NASA ADS] [CrossRef] [Google Scholar]
 MartíVidal, I., Muller, S., Vlemmings, W., Horellou, C., & Aalto, S. 2015, Science, 348, 311 [NASA ADS] [CrossRef] [Google Scholar]
 McKinney, J. C. 2005, ApJ, 630, L5 [NASA ADS] [CrossRef] [Google Scholar]
 McKinney, J. C., & Blandford, R. D. 2009, MNRAS, 394, L126 [NASA ADS] [CrossRef] [Google Scholar]
 McKinney, J. C., Tchekhovskoy, A., & Blandford, R. D. 2012, MNRAS, 423, 3083 [NASA ADS] [CrossRef] [Google Scholar]
 Meliani, Z., & Keppens, R. 2010, A&A, 520, L3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Meliani, Z., Sauty, C., Tsinganos, K., & Vlahakis, N. 2004, A&A, 425, 773 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Meliani, Z., Sauty, C., Vlahakis, N., Tsinganos, K., & Trussoni, E. 2006, A&A, 447, 797 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Meliani, Z., Sauty, C., Tsinganos, K., Trussoni, E., & Cayatte, V. 2010, A&A, 521, A67 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mertens, F., Lobanov, A. P., Walker, R. C., & Hardee, P. E. 2016, A&A, 595, A54 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mignone, A., Plewa, T., & Bodo, G. 2005, ApJS, 160, 199 [NASA ADS] [CrossRef] [Google Scholar]
 Nagai, H., et al. 2014, ApJ, 785, 53 [NASA ADS] [CrossRef] [Google Scholar]
 Narayan, R., & Yi, I., 1994, ApJ, 428, L13 [NASA ADS] [CrossRef] [Google Scholar]
 Narayan, R., Igumenshchev, I. V., & Abramowicz, M. A. 2003, PASJ, 55, L69 [NASA ADS] [Google Scholar]
 Nathanail, A., & Contopoulos, I. 2014, ApJ, 788, 186 [NASA ADS] [CrossRef] [Google Scholar]
 Parker, E. N. 1957, Phys. Rev., 107, 924 [NASA ADS] [CrossRef] [Google Scholar]
 Parker, E. N. 1963, Interplanetary Dynamical Processes (New York: Interscience Publishers Division, John Wiley and Sons) [Google Scholar]
 Petropoulou, M., & Dermer, C. 2016, ApJ, 825, L11 [NASA ADS] [CrossRef] [Google Scholar]
 Prieto, M. A., FernándezOntiveros, J. A., Markoff, S., Espada, D., & GonzálezMartín, O. 2016, MNRAS, 457, 3801 [NASA ADS] [CrossRef] [Google Scholar]
 Sauty, C., & Tsinganos, K. 1994, A&A, 287, 893 [NASA ADS] [Google Scholar]
 Sauty, C., Tsinganos, K., & Trussoni, E. 1999, A&A, 348, 327 [NASA ADS] [Google Scholar]
 Sauty, C., Trussoni, E., & Tsinganos, K. 2002, A&A, 389, 1068 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Sauty, C., Trussoni, E., & Tsinganos, K. 2004, A&A, 421, 797 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Sauty, C., Lima, J. J. G., Iro, N., & Tsinganos, K. 2005, A&A, 432, 687 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Sikora, M. 2011, in Jets at all Scales, eds. G. Romero, R. Sunyaev, & T. Bellon (Cambridge: Cambridge Univ. Press), IAU Symp. 275, 59 [Google Scholar]
 Sikora, M., Rutkowski, M., & Begelman, M. C. 2016, MNRAS, 457, 1352 [NASA ADS] [CrossRef] [Google Scholar]
 Silvestro, G., Ferrari, A., Rosner, R., Trussoni, E., & Tsinganos, K. 1987, MNRAS, 325, 228 [Google Scholar]
 Sol, H., Pelletier, G., & Asseo, E. 1989, MNRAS, 237, 411 [NASA ADS] [Google Scholar]
 Taub, A. H. 1948, PhRv, 74, 328 [Google Scholar]
 Tavecchio, F. 2016, AIP Conf. Proc., 1792, 020007 [CrossRef] [Google Scholar]
 Tavecchio, F., & Ghisellini, G. 2016, MNRAS, 456, 2374 [NASA ADS] [CrossRef] [Google Scholar]
 Tchekhovskoy, A. 2015, in The Formation and Disruption of Black Hole Jets, eds. I. Contopoulos, D. Gabuzda, & N. Kylafis, Astrophys. Space Sci. Lib., 414, 45 [NASA ADS] [CrossRef] [Google Scholar]
 Tchekhovskoy, A., Narayan, R., & McKinney, J. C. 2011, MNRAS, 418, L79 [NASA ADS] [CrossRef] [Google Scholar]
 Tchekhovskoy, A., & McKinney, J. C. 2012, MNRAS, 423, L55 [NASA ADS] [CrossRef] [Google Scholar]
 Tsinganos, K. 1982, ApJ, 252, 775 [NASA ADS] [CrossRef] [Google Scholar]
 Tsinganos, K., & Trussoni, E. 1991, A&A, 249, 156 [NASA ADS] [Google Scholar]
 Tsinganos, K., & Sauty, C., 1992, A&A, 255, 405 [Google Scholar]
 Urry, C. M., & Padovani, P. 1995, PASP, 107, 803 [NASA ADS] [CrossRef] [Google Scholar]
 van Eerten, H. J., Meliani, Z., Wijers, R. A. M. J., & Keppens, R. 2011, MNRAS, 410, 2016 [NASA ADS] [Google Scholar]
 Vlahakis, N. 2015, in The Formation and Disruption of Black Hole Jets, eds I. Contopoulos, D. Gabuzda, & N. Kylafis Astrophys. Space Sci. Lib. 414, 177 [Google Scholar]
 Vlahakis, N., & Tsinganos, A. 1998, MNRAS, 298, 777 [NASA ADS] [CrossRef] [Google Scholar]
 Vlahakis, N., & Königl, A. 2003a, ApJ, 596, 1104 [NASA ADS] [CrossRef] [Google Scholar]
 Vlahakis, N., & Königl, A. 2003b, ApJ, 596, 1080 [NASA ADS] [CrossRef] [Google Scholar]
 Vlahakis, N., & Königl, A. 2004, ApJ, 605, 656 [NASA ADS] [CrossRef] [Google Scholar]
 Vovk, I., & Babić, A. 2015, A&A, 578, A92 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Zamaninasab, M., ClausenBrown, E., Savolainen, T., & Tchekhovskoy, A. 2014, Nature, 510, 126 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Zdziarski, A. A., Sikora, M., Pjanka, P., & Tchekhovskoy A. 2015, MNRAS, 451, 927 [NASA ADS] [CrossRef] [Google Scholar]
All Tables
All Figures
Fig. 1 Relative error on the electric force for a recollimating oscillating solution in Kerr metric (K1, see Sect. 5.1). Color isocontours correspond to the relative error in the electric force. Field lines anchored into the black hole magnetosphere and in the accretion disk are plotted in black solid lines. The limiting field line between the inner jet coming from black hole corona and the outflow outgoing from the accretion disk is plotted in red. The “light cylinder” is indicated by a green solid line. The cylindrical radius and the distance above the equatorial plane are in units of Schwarzschild radius. 

In the text 
Fig. 2 Evolution of the radial velocity along the polar axis for solutions in a Schwarzschild metric, with m_{1} = 0 (blue) and m_{1} = −0.078 (red). The second case has a smaller terminal velocity. 

In the text 
Fig. 3 Field lines for a solution in a Schwarzschild metric, with parameters λ = 1.0, κ = 0.2, δ = 1.2, ν = 0.8, ℓ = 0, μ = 0.1, e_{1} = 0 and m_{1}= 0 (left) and m_{1}= −0.078 (right). We note that the case m_{1} = −0.078 corresponds to a more tightly collimated jet. Lengths are in units of the Schwarzschild radius. The red lines are connected to the magnetosphere of the central object while the green lines are connected to the disk. The separating line is in blue and the light cylinder in black. 

In the text 
Fig. 4 Threedimensional (3D) representation of the field lines and streamlines for the thermally collimated solution K1 at the base of the jet and for two flux tubes. The blue lines correspond to streamlines, the red lines to magnetic field lines. The length is in units of the Alfvén radius, that is, ten times the Schwarzschild radius. 

In the text 
Fig. 5 Poloidal field lines and “light cylinder” for the thermally collimated solution K1, for λ = 1.0, κ = 0.2, δ = 2.3, ν = 0.9, μ = 0.1, ℓ = 0.05, e_{1} = 0. The length unit is the Schwarzschild radius. 

In the text 
Fig. 6 Lorentz factor γ for the K1 (blue line) and K2 (red line) solutions. Distances are given in Schwarzschild radius units. 

In the text 
Fig. 7 Poloidal field lines and “light cylinder” for the K2 solution, that is, for λ = 1.0, κ = 0.2, δ = 1.35, ν = 0.46223, μ = 0.1, ℓ = 0.05, e_{1} = 0. Distances are given in Schwarzschild radius units. 

In the text 
Fig. 8 Longitudinal forces, that is, along the field line, for the K2 solution, along the line α = 0.01α_{lim}. Distances are given in Schwarzschild radius units. 

In the text 
Fig. 9 Transverse forces, that is, perpendicular to the field line, for the K2 solution, along the line α = 0.01α_{lim}. We see that the Lorentz force is collimating and is balanced by the electric force that decollimates. Distances are given in Schwarzschild radius units. 

In the text 
Fig. 10 Relative normalized contribution to the total energy of the kinetic enthalpy hγξ_{k} and the external heating distribution hγQ∕c^{2} along the axis. 

In the text 
Fig. 11 Relative normalized contribution to the total energy of the kinetic enthalpy hγξ_{k} and the external heating distribution hγQ∕c^{2} along the streamline α = 0.05α_{lim} 

In the text 
Fig. 12 Lorentz factor for the nonoscillating collimated solution K3. 

In the text 
Fig. 13 Poloidal field lines and “light cylinder” for the nonoscillating collimated solution K3, that is, for λ = 1.2, κ = 0.005, δ = 2.3, ν = 0.42, μ = 0.08, ℓ = 0.024, e_{1} = 0. Distances are given in Schwarzschild radius units. 

In the text 
Fig. 14 Poloidal field lines and “light cylinder” for the conical solution K4, that is, for λ = 0.0143, κ = 1.451, δ = 3.14, ν = 0.8, μ = 0.41, ℓ = 0.15, e_{1} = 0. Distances are given in Schwarzschild radius units. 

In the text 
Fig. 15 Lorentz factor for the conical solution K4. Distances are given in Schwarzschild radius units. 

In the text 
Fig. 16 Variation of the magnetic collimation efficiency parameter ϵ vs. the black hole spin parameter a_{H} for K1type, K3type, and K4type solutions. The value of ϵ has been normalized by ϵ(0) when it is negative and by ϵ(0) when it is positive. 

In the text 
Fig. 17 Cylindrical jet radius normalized to its value at the Alfvén surface, G, for K1type solutions, as a function of the distance along the polar axis, for five different values of the black hole spin a_{H}. The functionG is equal to 1 at the Alfvén distance r = 10r_{s}. 

In the text 
Fig. 18 Lorentz factor γ for K4type solutions when the black hole spin parameter a_{H} varies between − 0.99 and 0.99. 

In the text 
Fig. 19 Ratio of the cylindrical radius to the spherical radius for the conical K4type solutions with a spin parameter a_{H} varying between − 0.99 and 0.99 vs. the distance z above the equatorial plane in units of the Schwarzschild radius. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.