## Jet formation in GRBs: A semi-analytic model of MHD flow in Kerr geometry with realistic plasma injection

Noemie Globus 1 and Amir Levinson 1

## ABSTRACT

We construct a semi-analytic model for MHD flows in Kerr geometry, that incorporates energy loading via neutrino annihilation on magnetic field lines threading the horizon. We compute the double-flow structure for a wide range of energy injection rates, and identify the different operation regimes. At low injection rates the outflow is powered by the spinning black hole via the Blandford-Znajek mechanism, whereas at high injection rates it is driven by the pressure of the plasma deposited on magnetic field lines. In the intermediate regime both processes contribute to the outflow formation. The parameter that quantifies the load is the ratio of the net power injected below the stagnation radius and the maximum power that can be extracted magnetically from the black hole.

## 1. Introduction

An issue of considerable interest in the theory of gamma-ray bursts (GRBs) is the nature of the gammaray emitting jet. The conventional wisdom has been that the jet is produced by a hyper-accreting black hole that results from a neutron star merger in case of short GRBs (Eichler et al. 1989), or the core-collapse of a massive star in case of long GRBs (MacFadyen, Woosley 1999). The black hole is likely to be immersed in a strong magnetic field ( B H /similarequal 10 15 G) seeded by the progenitor and advected inwards during the formation of the central engine. Feedback from a rapidly rotating black hole may dominate the torque experienced by the surrounding torus, leading to a state of suspended accretion in long GRBs (van Putten &amp; Ostriker, 2001; van Putten &amp; Levinson 2003), provided that mass loading of magnetic field lines anchored to the disk is, somehow, strongly suppressed (Komissarov &amp; Barkov 2009; Globus &amp; Levinson 2013; hereafter GL13). Rapid heating of the inner regions of the hyper-accretion disk, or the torus in the suspended accretion state if established, leads to prodigious emission of MeV neutrinos (and anti-neutrinos), with luminosities in the range L ν = 10 51 -10 54 erg s -1 , depending on accretion rate and specific angular momentum a of the black hole (e.g., Popham et al. 1999; Chen &amp; Beloborodov 2007).

In the context of the picture outlined above, two competing jet formation mechanisms have been widely discussed in the literature; magnetic extraction of the spin down power of a Kerr black hole, and outflow formation via neutrino annihilation in the polar region, above the horizon (Paczy´ski 1990, Levinson &amp; n Eichler 1993, Levinson 2006). These two processes are commonly treated under idealized conditions: models of Blandford-Znajek jets usually invoke the force-free limit and ignore loading of magnetic field lines (but c.f., Komissarov &amp; Barkov 2009), whereas models of jets driven by νν ¯ annihilation (MacFadyen &amp; Woosley 1999, Fryer &amp; M´ esz´ros 2003) are usually constructed within the pure hydrodynamic limit. a In general, however, both processes might be at work, and it is desirable to characterize the interplay between them. The approach undertaken in this paper is to treat νν ¯ annihilation in the magnetosphere as external plasma load. It has been shown elsewhere (GL13) that injection of relativistically hot plasma on horizon threading field lines

1 School of Physics &amp; Astronomy, Tel Aviv University, Tel Aviv 69978, Israel

always leads to the formation of a double-flow structure in the magnetosphere. The plasma inflowing into the black hole carries positive energy that tends to counteract the BZ process. The plasma outflowing to infinity contributes to the total asymptotic power. The question addressed in this paper is how the structure of the MHD flow depends on the details of the plasma injection process.

In a preliminary investigation (GL13), we considered the effect of the load on the activation of the BZ mechanism, assuming that the plasma source is confined to an infinitely thin layer, outside which the MHD flow is ideal and adiabatic. We derived solutions for the inflow section only, and evaluated the critical load above which the BZ process switches off. We then argued that this critical value differentiates magnetically extracted from pressure-driven flows when the injected plasma is relativistically hot. In this paper we extend our analysis to more realistic situations, and obtain solutions of the MHD equations for the entire double-flow structure. We focus on the conditions anticipated in GRBs, and employ a realistic injection profile computed recently by Zalamea &amp; Beloborodov (2011, hereafter ZB11). We identify the different operation regimes, including the intermediate regime where the transition from magnetically extracted to pressure driven flows occurs. We also calculate overloaded solutions for which the BZ process is switch off and compare them with pure hydrodynamic flows derived in a previous study (Levinson &amp; Globus 2013, hereafter LG13).

## 2. Model

The double-flow structure established in the magnetosphere is illustrated in Figure 1: plasma inflow into the black hole and outflow to infinity are ejected from a stagnation radius, r st , located between the inner and outer light surfaces. The plasma consists of relativistically hot e ± pairs created via annihilation of MeV neutrinos emitted from the surrounding accretion disk. The exact location of the stagnation surface depends, quite generally, on the energy injection profile, and is treated as an eigenvalue of the MHD equations. The double flow possesses six critical surfaces, corresponding to the characteristic phase speeds of the three MHD waves propagating in the medium: two slow magnetosonic, two Alfv´ enic and two fast magnetosonic. We consider an infinitely conducting, stationary and axisymmetric flow. In general, the flow is characterized by a stream function Ψ( r, θ ) that defines the geometry of magnetic flux surfaces, and by the following functionals of Ψ: the angular velocity of magnetic field lines Ω F (Ψ), the ratio of mass and magnetic fluxes η (Ψ), and the energy, angular momentum and entropy per baryon, denoted by E (Ψ), L (Ψ) and s (Ψ), respectively. These quantities are given explicitly in Equations (A7)-(A10).

The ideal MHD condition implies that Ω F (Ψ) is conserved along magnetic flux tubes, as usual. All other quantities change along streamlines, owing to plasma injection by the external source, according to Equations (A11)-(A14). We assume that in the acceleration zone the plasma is relativistically hot with a negligible baryonic content. We can therefore adopt the equation of state w = ρc h 2 ¯ = 4 p . To simplify the analysis we invoke a split monopole configuration for the magnetic field lines, described by a stream function of the form Ψ( r, θ ) = Ψ (1 0 -cos θ ). The energy, angular momentum and entropy fluxes, Equation (A6), then have only a radial component: /epsilon1 r = ρ E u r , l r = ρ L u r , s r = ρsu r . Note that E , L and s diverge in the limit ρ → 0 (baryon-free flow), whereas the corresponding fluxes, /epsilon1 r , l r and s r , remain finite and are well defined also in the baryon-free case. For the relativistic equation of state adopted above the entropy per unit volume is given by S = ρs = 4 p/kT ∝ p 3 / 4 , whereby s r ∝ u p r 3 / 4 , hence the pressure p can be used as a free variable instead of S . Since our analysis encompases the force-free limit, we find it convenient to use /epsilon1 r , l r and p as our free variables. With the above simplifications, Equations (A11)-(A13) reduce to:

$$\frac { 1 } { \sqrt { - g } } \partial _ { r } ( \sqrt { - g } \epsilon ^ { r } ) = - q _ { t },$$

$$\frac { 1 } { \sqrt { - g } } \partial _ { r } ( \sqrt { - g } l ^ { r } ) & = q _ { \varphi }, \\ 3 _ { \infty } \quad, \quad \partial _ { r } \quad \omega \quad u _ { \alpha } q ^ { \alpha } & \quad \omega \quad.$$

$$\frac { \ddot { 3 } } { 4 } \partial _ { r } \ln p = - \partial _ { r } \ln ( \Sigma u ^ { r } ) - \frac { u _ { \alpha } q ^ { \alpha } } { 4 p u ^ { r } },$$

here √ -g = Σsin θ , and Σ = r 2 + a 2 cos 2 θ . This set needs to be augmented by an equation of motion for the velocity u r . Instead of u r we use the poloidal velocity u p = √ u u r r = √ g rr u r . Its rate of change along streamlines is derived in appendix A and can be written in the form,

$$\partial _ { r } \ln u _ { p } = \frac { N _ { a d } + N _ { q } } { D },$$

where N ad , N q and D are functionals of /epsilon1 r , l r , p, u p , Ω , given explicitly by Equations (A23)-(A25). Equations F (1)-(4) form a complete set that governs the structure of the double MHD flow. The solutions for the radial profiles of the free variables /epsilon1 r , l r , p , and u p , depend on the particular choice of the angle θ that characterizes magnetic flux surfaces. The angular velocity Ω F ( θ ) is given as an input. The energy and angular momentum flow rates per solid angle (along a particular flux surface) are defined, respectively, as:

$$\dot { \mathcal { E } } ( r, \theta ) \equiv \Sigma \epsilon ^ { r }, \quad \dot { \mathcal { L } } ( r, \theta ) \equiv \Sigma l ^ { r }.$$

## 2.1. Source terms for the process νν ¯ → e + e -

The energy-momentum deposition rate by the reaction νν ¯ → e + e -was computed in a number of works, under different simplifying assumptions (Popham et al. 1999, Chen &amp; Beloborodov 2007, ZB11). In what follows we use the recent analysis by ZB11 which includes general relativistic effects. Following ZB11 we denote by Q α νν ¯ the local energy-momentum deposition rate measured by a zero-angular-momentum observer (ZAMO). In general, those rates are functions of the Boyer-Lindquist coordinates r , θ and ϕ , as can be seen from Figures 2 and 3 in ZB11. In terms of the metric components defined in appendix A, g ϕϕ = /pi1 2 , g tϕ = -ωg ϕϕ , g tt = -α 2 + ω g 2 ϕϕ and g rr , we have the following relations between the ZAMO rates Q α νν ¯ and the source terms q α measured by a distant observer: αq t = Q t νν ¯ , /pi1q ϕ = Q ϕ νν ¯ + /pi1ωQ t νν ¯ /α , √ g rr q r = Q r νν ¯ . From that we obtain

$$- \, q _ { t } \ & = \ \alpha Q ^ { t } _ { \nu \bar { \nu } } + \varpi \omega Q ^ { \varphi } _ { \nu \bar { \nu } }, \\ \ \hat { \ } = \ \frac { } { } - \varpi \, & \quad \, \frac { } { } \, \nu \, \cdot$$

$$q _ { \varphi } \ = \ \varpi Q _ { \nu \bar { \nu } } ^ { \varphi },$$

$$u _ { \alpha } q ^ { \cdot ^ { \cdot } } \ = \ - \alpha u ^ { t } Q ^ { t } _ { \nu \bar { \nu } } + \varpi Q ^ { \varphi } _ { \nu \bar { \nu } } ( u ^ { \varphi } - \omega u ^ { t } ) + u _ { p } Q ^ { r } _ { \nu \bar { \nu } }.$$

The total power deposited in the magnetosphere can be expressed as,

$$\dot { E } ^ { t o t } _ { \nu \bar { \nu } } & = \int _ { r \geq r _ { H } } \left ( \alpha Q ^ { t } _ { \nu \bar { \nu } } + \varpi \omega Q ^ { \varphi } _ { \nu \bar { \nu } } \right ) \sqrt { - g } d r d \theta d \varphi \,. \\ \dots \dots \dots \dots \dots \pi \pi \pi \dots \dots \dots \dots \dots \pi \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \end{array}$$

A fit to the numerical results by ZB11 yields: ˙ E tot νν ¯ /similarequal 10 52 ˙ m 9 / 4 acc x -4 8 . mso erg s -1 for a black hole mass M BH = 3 M /circledot , and accretion rates (henceforth measured in units of M /circledot s -1 ) in the range 0 02 . &lt; m ˙ acc &lt; 1, where x mso is the radius of the marginally stable orbit in units of m = GM BH /c 2 .

Unfortunately, ZB11 do not exhibit results for the azimuthal term Q ϕ νν ¯ . It is also difficult to fit their result for Q r νν ¯ . We shall therefore set Q ϕ νν ¯ = Q r νν ¯ = 0. This should not alter much q t and u α q α , as the first term on the right hand side of Equations (6) and (8) dominates anyhow. However, for this choice q ϕ = 0,

implying that the angular momentum flow rate, ˙ ( L r, θ ) ≡ Σ l r , is conserved, as readily seen from Equation (2).

For our radial flow model it is sufficient to use the angle-averaged energy deposition rate. We adopt the form

$$Q ^ { t } _ { \nu \bar { \nu } } ( r ) = \dot { Q } _ { 0 } f ( x ),$$

where x = r/m is a fiducial radius, and f ( x ) is normalized such that f (1) = 1. From figures 2 and 3 in ZB11 we obtain the approximate profile f ( x ) /similarequal x -b , with b = 4 5 for a black hole spin parameter ˜ . a ≡ a/m = 0 95, . and b = 3 5 for ˜ = 0. The injected power per solid angle, from the horizon to a given radius . a x , is then given by

$$\dot { \mathcal { E } } _ { \nu \bar { D } } ( x, \theta ) = m ^ { 3 } \dot { Q } _ { 0 } \int _ { x _ { H } } ^ { x } \alpha \Sigma f ( x ^ { \prime } ) d x ^ { \prime }. \\ \text{information in to some hominom how to } \partial < \theta < \pi / \partial \text{i.s}$$

The cumulative power distribution in the upper hemisphere (0 ≤ θ ≤ π/ 2) is

$$\dot { E } _ { \nu \bar { \vartheta } } ( x ) & = 2 \pi \int _ { 0 } ^ { \pi / 2 } \dot { \mathcal { E } } _ { \nu \bar { \vartheta } } ( x, \theta ) \sin \theta d \theta. \\ \ddots \.$$

It is related to the total power through ˙ E tot νν ¯ = ˙ E νν ¯ ( x = ∞ ). The amount absorbed by the black hole along a particular field line equals the power per solid angle injected in the inflow section (between the horizon and the stagnation radius):

$$\i r a u s \} \colon & & \dot { \mathcal { E } } _ { \nu \bar { p } } ^ { i n } ( \theta ) = \dot { \mathcal { E } } _ { \nu \bar { p } } ( x _ { s t }, \theta ) = m ^ { 3 } \dot { Q } _ { 0 } \int _ { x _ { H } } ^ { x _ { s t } } \alpha \Sigma x ^ { - b } d x. & & ( 1 3 ) \\ \dot { \lambda } t o t ( \alpha ) \quad \dot { \lambda } i n ( \alpha ) \quad \dot { \lambda } t o t ( \alpha ) \ = \ \dot { \lambda } \ \ \partial _ { 0 } \quad \cdots \ e \ \cdots, \quad \tau \ \tau \ \cdots.$$

The rest, ˙ E out νν ¯ ( θ ) = ˙ E tot νν ¯ ( θ ) - E ˙ in νν ¯ ( θ ), where ˙ E tot νν ¯ ( θ ) ≡ E ˙ νν ¯ ( ∞ , θ ), emerges at infinity. The total power intercepted by the black hole in one hemisphere is

$$\dot { E } _ { \nu \bar { p } } ^ { i n } = 2 \pi \int _ { 0 } ^ { \pi / 2 } \dot { \mathcal { E } } _ { \nu \bar { p } } ^ { i n } ( \theta ) \sin \theta d \theta = \dot { E } _ { \nu \bar { p } } ( x _ { s t } ).$$

A plot of ˙ E νν ¯ ( x ) (Equation (12)), is exhibited in Figure 2. For our computations we use b = 4 5 for a spin . parameter ˜ = 0 95 (the solid line in Figure 2). a .

## 2.2. The load parameter κ θ ( )

Let us denote by ˙ E H ( θ ) = ˙ ( E r H , θ ), ˙ E st ( θ ) = ˙ ( E r st , θ ), and ˙ E ∞ ( θ ) = ˙ ( E ∞ , θ ) the angular distribution of the power at the horizon, stagnation radius and infinity, respectively, where ˙ ( E r, θ ) is defined in Equation (5). Integration of Equation (1) yields

$$\dot { \mathcal { E } } _ { H } ( \theta ) & = \dot { \mathcal { E } } _ { s t } ( \theta ) - \dot { \mathcal { E } } _ { \nu \bar { \mathcal { D } } } ^ { i n } ( \theta ), & ( 1 5 ) \\ \dot { \mathcal { C } } _ { \nu \lambda } & = \dot { \mathcal { C } } _ { \nu \lambda + \dot { \mathcal { C } } o u t } ( \theta ) & ( 1 \epsilon )$$

$$\dot { \mathcal { E } } _ { \infty } ( \theta ) = \dot { \mathcal { E } } _ { s t } ( \theta ) + \dot { \mathcal { E } } _ { \nu \bar { \mathcal { D } } } ^ { o u t } ( \theta ),$$

where Equations (6), (10) and (13) have been employed.

Now, the specific energy of the injected plasma is positive, hence ˙ E in νν ¯ ( θ ) ≥ 0, ˙ E out νν ¯ ( θ ) ≥ 0, as can be inferred from Equation (11). In the situations envisaged here ˙ E ∞ ( θ ) &gt; 0, but ˙ E H ( θ ) can be negative or positive, depending on the load. In the force-free limit ˙ E tot νν ¯ ( θ ) → 0, whereby Equations (15) and (16) yield

˙ E ∞ ( θ ) = ˙ E H ( θ ) &gt; 0. The extracted power per solid angle is given, in this limit, by (Blandford &amp; Znajek, 1977):

$$\dot { \mathcal { E } } _ { H } ( \theta ) = P _ { B Z } ( \theta ) \equiv \frac { c } { 6 4 \pi ^ { 3 } } \alpha _ { \Omega } ( 1 - \alpha _ { \Omega } ) \left ( \frac { \tilde { a } } { m } \right ) ^ { 2 } \frac { ( x _ { H } ^ { 2 } + \tilde { a } ^ { 2 } ) \sin ^ { 2 } \theta } { x _ { H } ^ { 2 } ( x _ { H } ^ { 2 } + \tilde { a } ^ { 2 } \cos ^ { 2 } \theta ) } \Psi _ { 0 } ^ { 2 } \,, \quad \, ( 1 7 ) \intertext { m o f t h o b l a b o l a m i n } \ e n v { \ v t h o b l a b o l a m i n } \tilde { \ s u m i n o i o n o c o n o r o n o c o n o c o n o c o n o c o n o c o n o c o n o c o n o c o n o c o n o c o n o c o n o c o n o c o n o c o n o c o n o c o n o c o n o c o n o c o n o c o n o c o n o c o n o c o n o c o n o c o n o c o n o c o n o c o n o c o n o c o n o c o n o c o n o c o n o c o n o c o n o c o n o c o n o c o n o c o n o c o n o c o n o c o n o c o n o c o n o c o n o c o n o c o n o c o n o c o n o c o n o } = \Omega _ { s } / / \ s.$$

in terms of the black hole spin ˜, magnetic flux Ψ , and the dimensionless parameter a 0 α Ω = Ω F /ω H . In general ˙ E H ( θ ) &lt; P BZ ( θ ), as readily seen from Equation (15).

Henceforth, we shall quantify the load on a specific streamline θ by the parameter

$$\kappa ( \theta ) \equiv \frac { \dot { \mathcal { E } } _ { \nu \bar { \mathcal { H } } } ^ { i n } ( \theta ) } { P _ { B Z } ( \theta ) }.$$

At κ θ ( ) &lt;&lt; 1 the flow along the streamline θ is nearly force-free. At κ θ ( ) &gt;&gt; 1 the flow is nearly hydrodynamic, whereby ˙ E st ( θ ) /similarequal 0 and ˙ E H ( θ ) /similarequal -E ˙ in νν ¯ ( θ ) &lt; 0, namely the energy injected below the stagnation radius on the flux surface θ is completely absorbed by the black hole. As shown below, the transition between the two regimes occurs, in general, at κ θ ( ) /similarequal 1.

## 3. Integration method

In general, the double flow must pass smoothly through 6 critical surfaces. Solutions that satisfy this requirement can be obtained, in principle, only if the trans-field equation is solved simultaneously with Equations (A11)-(A13), as the exact location of the critical surfaces is contingent upon the actual shape of the magnetic surfaces. Such an analysis is beyond the scope of this paper. Fixing the geometry of magnetic field lines renders the system of MHD equations, Eqs. (1)-(4), over constrained. The reason is that there are 4 regularity conditions (the regularity conditions at the inner and outer Alfv´n surfaces are automatically e satisfied), but only 3 adjustable parameters; the location of the stagnation point, r st , at which u p = 0, and the values of /epsilon1 r and p at r st , henceforth denoted by /epsilon1 r st , and p st , respectively. The value of the angular momentum flux at r st is related to /epsilon1 r st through the Bernoulli condition, Equation (A31): l r st = Ω -1 F /epsilon1 r st . Since we are interested in determining the energy flux on the horizon, we seek solutions that are regular on all 3 inner surfaces, and on the outer slow magnetosonic surface, but not necessarily on the outer fast magnetosonic surface. In practice we find that any solution that crosses the outer slow magnetosonic surface, is also regular on the outer Alfv´n surface. e

Our strategy is to start with some initial guess for the three adjustable parameters, r st , /epsilon1 r st , and p st , whereupon Equations (1)-(4) are integrated numerically from r st inwards to the horizon, and outwards to the outer Alfv´n surface, along a given field line e θ . The integration is repeated many times, where in each run the values r st , /epsilon1 r st , p st are readjusted until a solution that crosses the desired critical points smoothly is obtained. In all the examples presented below we neglected the change in linear and angular momentum owing to plasma injection, that is, we set Q ϕ νν ¯ = Q r νν ¯ = 0 in Equations (6)-(9). Equation (2) then readily yields ˙ L = Σ l r = const. The black hole spin ˜, angular velocity of the streamlines Ω a F , and load parameter κ are given as input parameters. In practice, however, κ cannot be determined a priori, since r st is unknown. We therefore use instead of κ the dimensionless parameter ˜ p B ≡ Ψ 2 0 c/ (32 π Q m 3 ˙ 0 5 ) as an indicator for the load. Once a solution is obtained and r st is determined, κ is computed by employing Equations (13), (17) and (18).

## 4. Results

The family of solutions can be divided into two classes that are distinguished by the sign of the energy flux on the horizon, /epsilon1 r H . This devision is dictated by the load parameter κ θ ( ), as discussed further below (Figure 4). We find that in case of underloaded solutions, defined as those for which κ θ ( ) &lt;&lt; 1, the specific energy is negative in the entire region encompassed by the plasma inflow (below the stagnation radius), including the horizon, whereby /epsilon1 r H &gt; 0. In case of overloaded solutions ( κ θ ( ) &gt;&gt; 1) we find /epsilon1 r H &lt; 0. Interestingly, the energy flux of overloaded solutions changes sign at some radius below the stagnation point, implying that there is still a region where the specific energy of the plasma inflow is negative. Typical examples are shown in Figure 3, where the velocity profiles (left panel) and the corresponding energy fluxes (right panel) of underloaded ( κ = 10 -5 ) and overloaded ( κ = 20) equatorial flows are exhibited, for a black hole spin parameter ˜ = 0 95, and an energy deposition profile a . f ( x ) = x -4 5 . . In the right panel we also exhibit solutions with κ = 0 4, 1, 6, and . κ = ∞ (a purely hydrodynamic flow), that are not shown in the left panel for clarity. As seen, the energy flux of the underloaded solutions is positive everywhere, whereas that of overloaded solutions changes sign below the stagnation radius. The location r 0 at which ˙ ( E r 0 , θ ) = 0 approaches r st as κ θ ( ) → ∞ . We think that this peculiar behavior stems from the fact that in the regime Ω F &lt; ω H the Poynting flux measured by a distant observer is always driven by the black hole (i.e., by frame dragging). To elucidate this point we employ Equation (A16) to obtain the Poynting flux on the horizon:

$$\left ( \frac { F ^ { r \theta } F _ { \theta t } } { 4 \pi } \right ) _ { H } = - \frac { \varpi _ { H } ^ { 2 } \Omega _ { F } ( \omega _ { H } - \Omega _ { F } ) } { M _ { H } ^ { 2 } + \varpi _ { H } ^ { 2 } ( \omega _ { H } - \Omega _ { F } ) ^ { 2 } } ( \epsilon _ { H } ^ { r } - \omega _ { H } l _ { H } ^ { r } ). \\ \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots$$

Now, /epsilon1 r H -ω H H l r = ρ H u r H ( E H -ω H L H ) is always negative, since the ZAMO energy is always positive, viz., E ZAMO = E ω L &gt; 0, and u r H &lt; 0. Thus, for any value of the load parameter, the electromagnetic flux on the horizon is positive if ω H &gt; Ω . Note also that in the force-free limit F M 2 H → 0, and (19) reduces to the familiar result, /epsilon1 r H = ( F rθ F θt ) H / 4 π , whereas in the pure hydrodynamic case M 2 H → ∞ and the Poynting flux vanishes, as expected. Figure 4 shows the electric current, I = α/pi1B ϕ , for two overloaded solutions, and it is seen that it never changes sign. We find that this is true in general in the regime 0 &lt; Ω F &lt; ω H , implying that for any value of κ there is a continuous flow of Poynting energy from the horizon outwards, against the inflow of injected plasma. In particular, at the stagnation radius B ϕ ( r st ) &lt; 0, and from Equation (A7) we obtain 1

$$\epsilon _ { s t } ^ { r } & = - \left ( \frac { \varpi \Omega _ { F } B _ { \varphi } B _ { r } } { 4 \pi \sqrt { g _ { r r } } } \right ) _ { s t } = \left ( \frac { F ^ { r \theta } F _ { \theta t } } { 4 \pi } \right ) _ { s t } > 0. \\ \intertext { a n d a t i n d } r & \leq \alpha \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \end{cases}$$

Since for overloaded solutions /epsilon1 r H &lt; 0, it is evident that the energy flux must vanish at some radius r &lt; r st , as seen in the right panel of Figure 3. Our interpretation is that the outward flow of electromagnetic energy driven by the black hole is counteracted by an inward flow of kinetic energy injected on magnetic flux tubes. When the latter exceeds the former the net energy flux becomes negative. For κ θ ( ) &lt;&lt; 1 this never happens, implying energy extraction from the black hole. For κ θ ( ) &gt;&gt; 1 this happens close to the stagnation radius, and since /epsilon1 r H &lt; 0 we infer that the black hole is being fed by the energy of the overloaded inflow (i.e., the energy injected below the stagnation radius).

Figure 5 displays the dependence of the outflow power on the load in the regime where rotational energy extraction is switched on, viz., /epsilon1 r H &gt; 0. For reference, powers are normalized to the equatorial BZ power,

1 This point was not properly understood in GL13. The claim made there, that for overloaded solutions B ϕ must vanish at r st is incorrect. However, the conclusion regarding the activation of the BZ process remains valid, as confirmed in the present analysis.

P BZ ( π/ 2). The horizontal axis gives values of the parameter ˙ E in νν ¯ ( θ /P ) BZ ( π/ 2), which for the equatorial flow ( θ = π/ 2) is just the load parameter κ θ ( ) defined in Equation (18). For other streamlines, the load parameter is obtained by multiplying values on the horizontal axis by the factor P BZ ( π/ 2) /P BZ ( θ ). The dashed lines delineate the normalized power extracted from the black hole, ˙ E H ( θ /P ) BZ ( π/ 2), and the solid lines the asymptotic power at infinity, ˙ E ∞ ( θ /P ) BZ ( π/ 2). These lines are essentially the locus of solutions obtained from the numerical integration of Equations (1)-(4). Specific cases are indicated by the symbols; the circles correspond to solutions for which Ω F = ω H / 2 and the triangles to solutions for which Ω F = ω H / 4. As seen, the effect of the load is highly insensitive to the value of Ω F , even though the structure of the flow does depend on this parameter (see Figure 6 below). This analysis confirms that the transition from underloaded to overloaded flows occurs at κ θ ( ) /similarequal 1. We also computed solutions for different black hole spins, and found the same behavior (see, for example, Figure 4 in GL13). One caveat is the possibility of a nonlinear feedback of the load on the magnetic flux in the vicinity of the horizon. Such a feedback may, in principle, change somewhat the activation condition, but not in a drastic way. It may be possible to test it using numerical simulations.

As explained above, in our model the angular velocity Ω F is given as an input. In reality it is determined by global conditions. For nearly force-free flows numerical simulations indicate that Ω F /similarequal ω H / 2. We therefore used this value for the underloaded solutions. However, when the inertia of the injected plasma becomes important, it is likely to affect Ω F . To study how the properties of the flow depend on this parameter, we sought solutions with different values of Ω F ( θ ), but the same value of κ θ ( ). An example is presented in Figure 6, and it is seen that while the velocity profile depends on Ω F , the power profile is insensitive to the choice of this parameter. In particular, it does not affect at all the activation condition. We find this trend is quite general, and therefore conclude that the result exhibited in Figure 5 is robust.

## 5. Conclusion

We constructed a semi-analytic model for the double-transonic flow established in the magnetosphere of a Kerr black hole under conditions anticipated in GRBs, incorporating plasma deposition on magnetic field lines via annihilation of MeV neutrinos emitted by the surrounding hyper-accretion flow. We examined the effect of energy loading on the properties of the flow, and identified the different operation regimes. We find that magnetic extraction of the black hole spin energy ensues, as long as the power deposited below the stagnation radius separating the inflow and outflow sections is smaller than the force-free BZ power. The transition from underloaded flows that are powered by the black hole spin energy, to overloaded flows that are powered by the neutrino source is continuous, as seen in Figure 5.

To relate the load parameter derived in Equation (18) to the accretion rate ˙ m acc (henceforth measured in units of M /circledot s -1 ), we employ the scaling relation derived in ZB11. As mentioned above, their analysis, that exploit an advanced disk model, yields a total energy deposition rate of

$$\dot { E } _ { \nu \bar { \nu } } ^ { t o t } \simeq 1 0 ^ { 5 2 } \left ( M _ { B H } / 3 M _ { \odot } \right ) ^ { - 3 / 2 } \dot { m } _ { a c c } ^ { 9 / 4 } x _ { m s o } ^ { - 4. 8 } \ \text{erg s} ^ { - 1 },$$

for accretion rates in the range 0 02 . &lt; m ˙ acc &lt; 1, where x mso is the radius of the marginally stable orbit in units of m . Combining the latter result with the activation condition derived from Figure 5, and using the angle averaged energy deposition rate, yields a rough estimate for the accretion rate at which a transition from underloaded to overloaded solutions occurs:

$$\dot { m } _ { c } \simeq 1 \left ( \frac { M _ { B H } } { 3 M _ { \odot } } \right ) ^ { - 2 / 9 } \left ( \frac { \Psi _ { 0 } } { 1 0 ^ { 2 8 } \, G \, \text{cm} ^ { 2 } } \right ) ^ { 8 / 9 } f ( \tilde { a } ).$$

Here Ψ 0 is the magnetic flux accumulated in the vicinity of the horizon, and the function f (˜) is displayed a in Figure 7. According to this relation, when ˙ m acc &lt; m ˙ c the outflow is powered by the BZ process, whereas for ˙ m acc &gt; m ˙ c it is driven by the neutrino source. In reality, the magnetic flux Ψ 0 should also depend on the accretion rate, however, the sensitivity of this relation to the assumptions underlying the specific disk model adopted for its calculation renders it highly uncertain. Furtheremore, the presence of sufficiently strong magnetic field in the inner disk regions may affect the neutrino luminosity. For illustration, we use the disk model of Chen &amp; Beloborodov (2007) to estimate the magnetic flux. Unfortunately, it is difficult to derive scaling relations from the results presented in this paper, but from Figures 1 and 2 there we obtained Ψ 0 ∼ 2 × 10 28 √ ξ B G cm 2 for a black hole mass M BH = 3 M /circledot , angular momentum ˜ = 0 95, viscosity a . parameter α vis = 0 1, and accretion rate . ˙ m acc = 0 2, assuming that the magnetic pressure in the inner . regions of the disk is a fraction ξ B of the total pressure. For this choice we infer that with ξ B on the order of a few percents, as naively expected, the transition from underloaded to overloaded flows may occur at accretion rates ˙ m acc &gt; 0 1 or so. . We emphasize that this estimate is highly uncertain.

The above results may also have some implications for the jet structure. To be concrete, the dependence of the activation condition on the inclination angle θ of magnetic surfaces (see Figure 5), and the approximate uniformity of the angular distribution of the energy deposition rate indicated in Figures 2 and 3 of ZB11, suggest that for accretion rates ˙ m acc &lt; ∼ ˙ m c , and unless the magnetic flux near horizon is extremely high, the outflow produced in the polar region may consist of an inner core inside which the power is dominated by the thermal energy of the hot plasma, and outside which it is dominated by the Poynting flux driven by frame dragging.

This research was supported by a grant from the Israel Science Foundation no. 1277/13

## A. Derivation of the flow equations

The stress-energy tensor of a magnetized fluid takes the form,

$$T ^ { \alpha \beta } & = \bar { h } \rho c ^ { 2 } u ^ { \alpha } u ^ { \beta } + p g ^ { \alpha \beta } + \frac { 1 } { 4 \pi } \left ( F ^ { \alpha \gamma } F ^ { \beta } _ { \gamma } - \frac { 1 } { 4 } g ^ { \alpha \beta } F ^ { 2 } \right ), \\ & \quad \quad \quad$$

here u α is the four-velocity measured in units of c, ¯ h = ( ρc 2 + e int + ) p /ρc 2 the dimensionless specific enthalpy, ρ the baryonic rest-mass density, p the pressure, and g µν the coefficients of the metric tensor of the Kerr spacetime. In the following we use geometrical units ( c = G = 1), unless otherwise stated, and express the Kerr metric in the regular Boyer-Lindquist coordinates, ds 2 ≡ g µν dx dx µ ν with the non-zero metric coefficients given by: g rr = Σ ∆, / g θθ = Σ, g ϕϕ ≡ /pi1 2 = A sin 2 θ/ Σ, g tt = -α 2 + ω g 2 φφ , g tφ = -ωg φφ , in terms of ∆ = r 2 + a 2 -2 mr , Σ = r 2 + a 2 cos 2 θ , A = ( r 2 + a 2 ) 2 -a 2 ∆sin 2 θ , α = √ Σ∆ /A , and ω = 2 mra/A . The parameters m and a are the mass and specific angular momentum per unit mass of the hole, α is the time lapse and ω the frame dragging potential between a zero-angular-momentum observer (ZAMO) and an observer at infinity. The angular velocity of the black hole is defined as the value of ω on the horizon, viz., ω H ≡ ω r ( = r H ) = a/ (2 mr H ), here r H = m + √ m 2 -a 2 is the radius of the horizon, obtained from the condition ∆ H = 0.

The dynamics of the flow is governed by the energy-momentum equations:

$$\frac { 1 } { \sqrt { - g } } ( \sqrt { - g } T ^ { \alpha \beta } ) _ {, \alpha } + \Gamma ^ { \beta } _ { \ \mu \nu } T ^ { \mu \nu } = q ^ { \beta },$$

mass conservation:

and Maxwell's equations:

$$\frac { 1 } { \sqrt { - g } } \partial _ { \alpha } ( \sqrt { - g } \rho u ^ { \alpha } ) = q _ { n },$$

$$F _ { ; \alpha } ^ { \beta \alpha } & = \frac { 1 } { \sqrt { - g } } ( \sqrt { - g } F ^ { \beta \alpha } ) _ {, \alpha } = 4 \pi j ^ { \beta }, \\ & F _ { \alpha \beta \cdot \gamma } + F _ { \beta \cdot \alpha } + F _ { \gamma \alpha \cdot \beta } = 0,$$

$$F _ { \alpha \beta, \gamma } + F _ { \beta \gamma, \alpha } + F _ { \gamma \alpha, \beta } = 0,$$

subject to the ideal MHD condition F µν u ν = 0. Here, q β denotes the source terms associated with energymomentum transfer by an external agent, q n is a particle source, and Γ β µν denotes the affine connection.

The energy, angular momentum and entropy fluxes, can be expressed explicitly as

$$\epsilon ^ { a } \equiv - T _ { t } ^ { a } = \rho u ^ { a } \mathcal { E }, \ \ l ^ { a } \equiv T _ { \varphi } ^ { a } = \rho u ^ { a } \mathcal { L }, \ \ s ^ { a } = ( \rho / m _ { N } ) u ^ { a } s,$$

in terms of the energy per baryon,

$$\mathcal { E } = - \bar { h } u _ { t } - \frac { \alpha \varpi \Omega _ { F } } { 4 \pi \eta } B _ { \varphi },$$

$$\mathcal { L } = \bar { h } u _ { \varphi } - \frac { \alpha \varpi B _ { \varphi } } { 4 \pi \eta },$$

and the entropy per baryon s , where

$$\alpha _ { F } = v ^ { \varphi } - \frac { v _ { p } B _ { \varphi } } { \varpi B _ { p } }$$

is the angular velocity of magnetic field lines,

$$\eta = \frac { \rho u _ { p } } { B _ { p } }$$

is the ratio of mass and magnetic fluxes, and the index a runs over r and θ . In the above equations u p = ± ( u u r r + u u θ θ ) 1 / 2 is the poloidal velocity, where the plus sign applies to outflow lines and the minus sign to inflow lines, v p = u /γ p , with γ = u α t being the Lorentz factor measured by a ZAMO, v ϕ = u ϕ /u t , B p = ( B 2 r + B 2 θ ) 1 / 2 /α is the redshifted poloidal magnetic field, and B r = F θϕ / √ A sin θ , B θ = √ ∆ F ϕr / √ A sin θ and B ϕ = √ ∆ F rθ / Σ the magnetic field components measured by a ZAMO (see van Putten and Levinson 2012 for details). Note that with our sign convention the value of η is defined to be positive on outflow lines and negative on inflow lines.

For a stationary and axisymmetric ideal MHD flow, Equations (A2)-(A5) can be reduced to (GL13)

$$\frac { 1 } { \sqrt { - g } } \partial _ { a } ( \sqrt { - g } \epsilon ^ { a } ) & = - q _ { t }, & ( \text{A11} ) \\ 1 \quad \sim \sim \frac { 1 } { 2 } \sim \frac { 1 } { 2 } \sim \frac { 1 } { 2 } \sim \cos$$

$$\stackrel { \kappa } { \frac { k T } { \sqrt { - g } } } \partial _ { a } ( \sqrt { - g } s ^ { a } ) & = - u _ { a } q ^ { \alpha }, & ( \text{A13} ) \\ \stackrel { w ^ { a } \lambda } { w ^ { a } \lambda } - \stackrel { u _ { p } q _ { n } } { \underline { u _ { p } q _ { n } } } & \text{$\ell \, \Lambda \, \Lambda$}$$

$$\frac { 1 } { \sqrt { - g } } \partial _ { a } ( \sqrt { - g } l ^ { a } ) & = q _ { \varphi }, & ( \text{A12} ) \\ k T \quad \sim \tilde { \sim } \longrightarrow \tilde { \sim } \tilde { \sim } \quad \tilde { \sim }$$

$$\begin{smallmatrix} \mathfrak { v } & \mathfrak { s } \\ u ^ { a } \partial _ { a } \eta = \frac { u _ { p } q _ { n } } { B _ { p } }, \\ \dots \dots \cdot \end{smallmatrix}$$

$$u ^ { a } \partial _ { a } \Omega ( \Psi ) & = 0.$$

angular momentum per baryon,

It can be readily seen that for q n = q µ = 0, the quantities Ω(Ψ), E (Ψ), L (Ψ), η (Ψ) and s (Ψ) are conserved on magnetic flux surfaces Ψ( r, θ ) =const, as is well known (e.g., Camenzind 1986).

From (A7)-(A10) we obtain the expressions:

$$B _ { \varphi } & = - \frac { 4 \pi \eta \mathcal { E } } { \alpha \varpi } \, \frac { \alpha ^ { 2 } \tilde { L } - \varpi ^ { 2 } ( \Omega _ { F } - \omega ) ( 1 - \tilde { L } \omega ) } { k _ { 0 } - M ^ { 2 } } \\, & \quad \mathcal { E } \, \alpha ^ { 2 } ( 1 - \Omega _ { F } \tilde { L } ) - M ^ { 2 } ( 1 - \omega \tilde { L } ) & \dots$$

$$u ^ { \varphi } = \frac { \cdot \cdot } { \overline { h } } \frac { \stackrel { \cdot \cdot } { \alpha } ^ { 2 } \Omega _ { F } ( 1 - \Omega _ { F } \tilde { L } ) - M ^ { 2 } \omega ( 1 - \omega \tilde { L } ) - M ^ { 2 } \tilde { L } \alpha ^ { 2 } \varpi ^ { - 2 } } { \alpha ^ { 2 } \left ( k _ { 0 } - M ^ { 2 } \right ) } \,, \quad \quad \quad$$

$$u ^ { t } & = \frac { \stackrel { \cdots } { \bar { h } } } { \bar { h } } \frac { \stackrel { \cdots } { \alpha ^ { 2 } ( 1 - \Omega _ { F } \tilde { L } ) - M ^ { 2 } ( 1 - \omega \tilde { L } ) } } { \alpha ^ { 2 } \left ( k _ { 0 } - M ^ { 2 } \right ) }, \\ \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

here ˜ = L L E / , M as the poloidal Alfv´nic Mach number, defined through e M 2 ≡ 4 πhη c /ρ ¯ 2 2 = u /u 2 p 2 A , and u 2 A = B / 2 p (4 πhρc ¯ 2 ). Combining the latter relations with the normalization condition u α u α = -1 yields the Bernoulli equation (Camenzind, 1986; Takahashi et al. 1990):

$$u _ { p } ^ { 2 } + 1 = \left ( \frac { \mathcal { E } } { \bar { h } } \right ) ^ { 2 } \frac { k _ { 0 } k _ { 2 } - 2 k _ { 2 } M ^ { 2 } - k _ { 4 } M ^ { 4 } } { ( k _ { 0 } - M ^ { 2 } ) ^ { 2 } } \,,$$

where

$$k _ { 2 } \ = \ \begin{array} { c c } ( 1 - \tilde { L } \Omega ) ^ { 2 }, \\ \tilde { r } \, 2 \quad ( 1 - \tilde { r }, \lambda ^ { 2 } ) \end{array}$$

$$k _ { 0 } \ & = \ \alpha ^ { 2 } - \varpi ^ { 2 } \left ( \Omega _ { F } - \omega \right ) ^ { 2 }, \\ 1. \ & = \ \bar { \ t } 1 \quad \tilde { \ t } \, \varpi ^ { 2 } \quad \bar { \ t } \, \varpi ^ { 1 } \, \searrow$$

$$k _ { 4 } \ = \ \frac { \lambda } { \varpi ^ { 2 } } - \frac { \lambda ^ { \cdot } } { \alpha ^ { 2 } } \,.$$

In terms of the free variables, /epsilon1 r , l r , p , used in our integration we have: η E = √ Σ /epsilon1 r / ( √ ∆ B p ), and E /h ¯ = /epsilon1 r / (4 pu r ), ˜ = L l r //epsilon1 r .

The equation of motion (4) is obtained upon differentiating Equation (A19) along a given streamline, using Equations (A11)-(A13) with the source terms q t = -αQ f x ˙ 0 ( ), q φ = 0 and u α q α = u q t t , which are derived in section 2.1, and noting that in the split monopole geometry, the redshifted poloidal field reduces to B p = Ψ 0 / (2 π √ Σ∆), the poloidal velocity is u p = √ Σ ∆ / u r , and the convective derivative reduces to u α ∂ α = u ∂ r r = √ ∆ Σ / u ∂ p r . This yields the following expressions for the functionals D N , ad and N q in Equation (4):

$$u \, \psi _ { \alpha } - u \, \psi _ { r } - \, \sqcup \, \psi _ { d } \, \psi _ { r }. \, \ i m s y e s u n e \, \text{inowing expressions} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{incase} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in class} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in cases} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in Case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{inCase} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{inclass} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{in case} \, \text{inclosing} \, \text{inclosing} \, \text{inclosing} \, \text{inclosing} \, \text{inclosing} \, \text{inclosing} \, \text{inclosing} \, \text{inclosing} \, \text{inclosing} \, \text{inclosing} \, \text{inclosing} \, \text{inclosing} \, \text{inclosing} \, \text{inclosing} \, \text{inclosing} \, \text{inclosing} \, \text{inclosing} \, \text{inclosing} \, \text{inclosing} \, \text{inclosing} \, \text{inclosing} \, \text{inclosing} \, \text{inclosing}$$

where the derivatives are defined by

$$\partial _ { x } k _ { 0 } \ = \ \underset { \mathfrak { g } _ { 0 } } { \partial _ { x } ( \alpha ^ { 2 } ) - ( \Omega _ { F } - \omega ) ^ { 2 } } \, \partial _ { x } ( \varpi ^ { 2 } ) + 2 \varpi ^ { 2 } \, ( \Omega _ { F } - \omega ) \, \partial _ { x } \omega, \\ \mathfrak { g } _ { 0 } \tilde { \tau } _ { 1 } \,,$$

$$\partial _ { x } k _ { a d } \ = \ \frac { 2 } { \alpha ^ { 2 } } ( 1 - \tilde { L } \omega ) ^ { 2 } \partial _ { x } \ln \alpha + \frac { 2 \tilde { L } \omega } { \alpha ^ { 2 } } ( 1 - \tilde { L } \omega ) \partial _ { x } \ln \omega - \frac { 2 \tilde { L } ^ { 2 } } { \varpi ^ { 2 } } \partial _ { x } \ln \varpi \,.$$

## A.1. The stagnation point

At the stagnation point x = x st , where u p = 0, Equations (1)-(4) with q ϕ = 0 yield:

$$\partial _ { x } u _ { p _ { | _ { x = x _ { s t } } } } \ = \ - \frac { ( \Sigma q _ { t } ) _ { s t } } { 4 p _ { s t } \alpha _ { s t } ( k _ { 0 s t } \, A _ { s t } ) ^ { 1 / 2 } },$$

$$= \ - \frac { ( \langle q _ { t } ) s t } { 4 p _ { s t } \alpha _ { s t } ( k _ { 0 s t } \, A _ { s t } ) ^ { 1 / 2 } }$$

$$( A 2 8 )$$

$$\partial _ { x } ( \Sigma \epsilon ^ { r } ) _ { | _ { x = x _ { s t } } } \ = \ - & ( \Sigma q _ { t } ) _ { s t } \,, \\ \intertext { a } \, \begin{matrix} \partial _ { x } ( \Sigma q _ { t } ) _ { s t } \,, \\ \end{matrix} \, \begin{matrix} \partial _ { x } ( \Sigma q _ { t } ) _ { s t } \,, \\ \end{matrix} \, \begin{matrix} \partial _ { x } ( \Sigma q _ { t } ) _ { s t } \,, \\ \end{matrix} \, \begin{matrix} \partial _ { x } ( \Sigma q _ { t } ) _ { s t } \,, \end{matrix} \, \begin{matrix} \partial _ { x } ( \Sigma q _ { t } ) _ { s t } \,, \end{matrix} \, \begin{matrix} \partial _ { x } ( \Sigma q _ { t } ) _ { s t } \,, \end{matrix} \, \begin{matrix} \partial _ { x } ( \Sigma q _ { t } ) _ { s t } \,, \end{matrix} \, \begin{matrix} \partial _ { x } ( \Sigma q _ { t } ) _ { s t } \,, \end{matrix} \, \begin{matrix} \partial _ { x } ( \Sigma q _ { t } ) _ { s t } \,, \end{matrix} \, \end{matrix}$$

$$\frac { 3 } { 4 } \partial _ { x } \ln ( p ) _ { | _ { x = x _ { s t } } } ^ { \mu ^ { \nu ^ { \nu ^ { \nu } } } } \ = \ - \frac { x _ { s t } - 1 } { \Delta _ { s t } } - \frac { x _ { s t } } { \Sigma _ { s t } } - \frac { \Sigma _ { s t } ^ { 2 } f ( x _ { s t } ) \epsilon _ { s t } ^ { r } ( 1 - \omega _ { s t } / \Omega _ { F } ) } { 8 p _ { s t } \alpha _ { s t } k _ { 0 s t } \tilde { p } _ { B } },$$

where -(Σ q t ) st = α st Σ st ˙ Q f x 0 ( st ), and the parameter ˜ p B ≡ Ψ 2 0 c/ (32 π Q m 3 ˙ 0 5 ) is an indicator for the load.

The Bernoulli condition can be rewritten

$$( 1 - \tilde { L } _ { s t } \Omega _ { F } ) \frac { \epsilon _ { s t } ^ { r } } { ( 4 p u ^ { r } ) _ { s t } } = \sqrt { k _ { 0 s t } } \,, & & ( A 3 1 ) \\ \otimes & \text{as we have shown. there is always extraction of angular momentum from the}$$

implying ˜ L st Ω F = 1, since as we have shown, there is always extraction of angular momentum from the black hole so that /epsilon1 r st &gt; 0.

## A.2. The Alfv´ en surfaces

The location of the Alfv´n surfaces is defined by the roots of the denominator in Equations (A16)-(A18), e M 2 A = k 0 A = α 2 A -/pi1 2 A (Ω F -ω A ) 2 , here the subscript A denotes values on this surface. The latter equation has two roots, corresponding to the inner and outer Alfv´ en surfaces. The requirement that B ϕ , u t and u ϕ in Equations (A16)-(A18) are continuous there imposes a condition on the ratio of the ZAMO energy, E ZAMO = E ω L , and the energy E of an observer at infinity:

$$\left ( \frac { \mathcal { E } } { \mathcal { E } ^ { Z A M O } } \right ) _ { A } = 1 - \alpha _ { A } ^ { - 2 } \varpi _ { A } ^ { 2 } \omega _ { A } ( \omega _ { A } - \Omega _ { F } ). & ( A 3 2 ) \\ \intertext { t h D 7 } \text{when to attributed } \tilde { r } \, = \, \partial r \, \mathcal { C } \, \circ \, \text{and} \, \, \text{and} \, \text{and} \, \text{and} \, \text{when} \, \mathcal { C } \, \circ \, \text{and} \, \text{when} \, \text{below}$$

In the regime where the BZ process is activated, /epsilon1 r H = ρu r E &gt; 0, and we must have E &lt; 0 anywhere below the stagnation radius where u r &lt; 0, and in particular at the inner Alfv´n point ( e IA ). Since the ZAMO energy is always positive, the latter condition, combined with Equation (A32), readily implies ω IA &gt; Ω , F and defines the range of Alfv´n radii that are permitted for energy extraction: e

$$\alpha _ { I A } \, \frac { \alpha _ { I A } } { \sqrt { \omega _ { I A } ( \omega _ { I A } - \Omega _ { F } ) } } & < \varpi _ { I A }. & ( A 3 3 ) \\ \vdots \text{marks this range for solutions with } \Omega _ { F } = \omega _ { H } / 2. \text{ In the outflow section all} \\ \tilde { \Phi } _ { \Phi } \,. & > \, \Omega _ { \Phi } \,. & \text{ and}$$

The shaded area in Figure 8 marks this range for solutions with Ω F = ω H / 2. In the outflow section all energies are positive, yielding ω OA &lt; Ω , and F

at the outer Alfv´n point ( e OA ).

$$\mathbb { m } { \omega } _ { O A } & \sim \mathbb { m } { F }, \, \Gamma _ { A } \\ & \quad \frac { \alpha _ { O A } } { \sqrt { \omega _ { O A } ( \omega _ { O A } - \Omega _ { F } ) } } > \varpi _ { O A }, & ( A 3 4 ) \\ \mathbb { m } { A } ).$$

## B. Pure hydrodynamic flows

The equations governing a purely hydrodynamic flow can be obtained formally from the MHD equations derived above upon taking the limit M 2 →∞ :

$$( u _ { p } ^ { 2 } - c _ { s } ^ { 2 } ) & \partial _ { x } \ln u _ { p } = N _ { a d } + N _ { q }, \\ / & \bullet \quad / \bullet 2 \circ \circ \bullet$$

$$N _ { q } = - \frac { 3 \, q _ { t } } { 2 \alpha ^ { 2 } \epsilon ^ { r } } ( 1 - \tilde { L } \omega ) \left ( \frac { \mathcal { E } } { \hbar } \right ) ^ { 2 } \left [ 1 + \left ( \frac { \mathcal { E } } { \hbar } \right ) ^ { 2 } \frac { 4 \, k _ { 4 } } { 3 } \right ] \,.$$

$$N _ { a d } = ( 1 + u _ { p } ^ { 2 } ) c _ { s } ^ { 2 } \left ( \frac { x } { \Sigma } + \frac { x - 1 } { \Delta } \right ) _ { \substack { \ell \\ \ell \lambda } } ^ { 2 } \left ( \frac { \mathcal { E } } { \bar { h } } \right ) ^ { 2 } \frac { 3 \partial _ { x } k _ { a d } } { 4 }, \\ \mathfrak { r } _ { \alpha }. \right ) _ { \substack { \ell \lambda } } ^ { 2 } \left [ \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \right ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( { \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \ \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \bottom ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left (  \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \top ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left (
 \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \bar { \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left (     \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left  \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( { { \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( {  \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( (  \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left { \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( ( { \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left ( \left                           \        \ {$$

$$\frac { 3 } { 4 } \partial _ { x } \ln \tilde { p } = - \partial _ { x } \ln u _ { p } - \frac { x } { \Sigma } - \frac { x - 1 } { \Delta } - \frac { q _ { t } ( 1 - \tilde { L } \omega ) } { \alpha ^ { 2 } \epsilon ^ { r } } \left ( \frac { \mathcal { E } } { \tilde { h } } \right ) ^ { 2 },$$

$$\partial _ { x } ( \Sigma \epsilon ^ { r } ) = - \Sigma q _ { t } \,, \quad \quad \quad$$

where now ˜ is a free parameter that describes the family of solutions, ˜ = L p p/Q t ˙ 0 d is the normalized pressure, with t d = GM BH /c 3 . The Bernoulli condition (A31) implies, in this limit, /epsilon1 r st = 0, as expected in the absence of magnetic fields. It can be readily shown that when M 2 → ∞ the slow-magnetosonic and Alfven speeds approach zero, whereas the fast-magnetosonic speed approaches the sound speed, c s = 1 / √ 2. Consequently, the above system of equations has critical points at u p = ± c s , as can be directly verified.

The regularity conditions at the sonic points, obtained from Equations (B2) and (B3), read:

$$2 \left ( \frac { x _ { c 1 } } { \Sigma _ { c 1 } } + \frac { x _ { c 1 } - 1 } { \Delta$$

$$2 \left ( \frac { x _ { c 2 } } { \Sigma _ { c 2 } } + \frac { x _ { c 2 } - 1 } { \Delta$$

/negationslash denoting the sonic point of the inflow (outflow) by x c 1 ( x c 2 ), and noting that u p = -1 / √ 2 at x c 1 , and u p = 1 / √ 2 at x c 2 . As seen, the existence of two sonic points, that is, x c 1 = x c 2 , is a consequence of energy injection. When f ( x ) = 0 the solutions of (B6) and (B7) merge, and the system has only one critical point, for either an inflow or an outflow, depending on the boundary conditions.

At the stagnation point, x = x st , the above equations yield

$$\tau \cdots \quad \partial _ { x } u _ { p _ { | x = \tilde { v } _ { s t } } } \ & = \$$

$$\partial _ { x } \ln ( \tilde { p } ) _ { | _ { x = x _ { s t } } } \ = \ \frac { 2 } {$$

$$\partial _ { x } u _ { p | _ { x = x _ { s t } } } & \ = \ \frac { \nu \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cd. } { 4 \tilde { p } _ { s t } \sqrt { \Delta _ { s t } } \alpha _ { s t } \sqrt { - k _ { 4, s t } } }, \\ \partial _ { x } ( \Sigma \epsilon ^ { r } ) _ { | _ { x = x _ { s t } } } & \ = \ \frac { - ( \Sigma q _ { t } ) _ { s t } \,, } { 2 \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdots }$$

where ˜ p st = ˜( p x st ) is the normalized stagnation pressure. Thus, for a given choice of f ( x ) the solution is fully determined once x st and ˜ p st are known, since /epsilon1 r st = 0.

## REFERENCES

Blandford, R. D. &amp; Znajek, R. L. 1977, MNRAS, 179, 433 Camenzind, M. 1986, ApJ, 162, 32 Chen, W-X. &amp; Beloborodov, A. 2007, ApJ, 657, 383 Eichler, D., Livio, M., Piran, T. &amp; Schramm, D. N. 1989, Nature, 340, 126 Fryer, C. L., &amp; M´sz´ros, P. 2003, ApJ, 588, L25 e a Globus, N. &amp; Levinson, A 2013, Phys. Rev. D, 88, 4046 Komissarov, S. &amp; Barkov, M. 2009, MNRAS, 397, 1153 Levinson, A. 2006, ApJ, 648, 510 Levinson, A. &amp; Eichler, D. 1993, ApJ, 418, 386 Levinson, A &amp; Globus, N. 2013, ApJ, 770, 159 MacFadyen, A. I. &amp; Woosley, S. E. 1999, ApJ, 524, 262 Paczy´ski, B. 1990, ApJ, 363, 218 n Popham, R., Woosley, S. E., &amp; Fryer, C. 1999, ApJ, 518, 356 Takahashi, M., Nitta, S., Tatematsu, Y. &amp; Tomimatsu, A. 1990, ApJ, 363, 206 van Putten, M. H. P. M. &amp; Levinson, A. 2003, ApJ, 584, 937 van Putten, M. H. P. M. &amp; Ostriker, E. C. 2001, ApJ, 552, 31

Thorne, K. S., MacDonald, D. A. &amp; Price, R. H. 1986, Black Holes: The Membrane Paradigm (Yale University Press)

van Putten, M. H. P. M. &amp; Levinson, A. 2012, Relativistic Astrophysics of the Transient Universe (Cambridge University Press)

Zalamea, I. &amp; Beloborodov, A. 2011, MNRAS, 410, 2302

This preprint was prepared with the AAS L T E X macros v5.2. A

Fig. 1.- Illustration of the double-transonic flow model.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0126v1-figure-499f1330-3efd-4770-8ed3-0a0323bfb111.jpg)

Fig. 2.- Net power deposited on magnetic field lines via neutrino annihilation as a function of radius x (Equation (12)), for two different injection profiles, f ( x ) = x -4 5 . and f ( x ) = x -3 5 . (˜ = 0 95). a .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0126v1-figure-64496ce9-5879-4548-bfb8-69e81b84a99c.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0126v1-figure-04e40c9f-5370-4b0c-ad4c-30cac543f933.jpg)

Fig. 3.- Left panel: velocity profiles of underloaded ( κ = 10 -5 ) and overloaded ( κ = 20) solutions, with f ( x ) = x -4 5 . , ˜ = 0 95, and a . θ = 90 . The region above (below) the horizontal dotted line ◦ u r = 0, corresponds to the outflow (inflow) sections. The inner and outer slow magnetosonic points (SMP), Alfv´n points (AP), e and fast magnetosonic points (FMP) are indicated. The vertical red line delineates the horizon. Right panel: profiles of the outflow power per solid angle for different values of the load parameter κ . The κ = 10 -5 and κ = 0 4 curves are rescaled for convenience. . The cross symbole on each curve marks the location of the stagnation radius. In the region of the inflow where the energy flux is positive (between the point of zero flux and the stagnation point) the specific energy is negative.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0126v1-figure-81086980-089c-4952-9c0b-9de9a10e2000.jpg)

Fig. 4.- Electric current distribution, I ( x ) = B /pi1α ϕ , for equatorial flow solutions.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0126v1-figure-98107f5a-6943-4aa2-aff5-1c04ceb9cd12.jpg)

Fig. 5.- Dependence of the outflow power on the load for two different streamlines, θ = 90 , and ◦ θ = 30 . ◦ For reference, powers are normalized by the equatorial BZ power, P BZ ( π/ 2), given in Equation (17). The dashed line in each case gives the normalized power per solid angle on the horizon, ˙ E H ( θ /P ) BZ ( π/ 2), and the solid line the total power per solid angle, ˙ E ∞ ( θ /P ) BZ ( π/ 2). Specific cases are indicated by the symbols, with the circles corresponding to solutions for which Ω F = ω H / 2, and the triangles to solutions with Ω F = ω H / 4.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0126v1-figure-4dbbbbc1-046c-4a9e-9948-388355e22518.jpg)

Fig. 6.- A comparison between two equatorial flow solutions with the same load parameter, κ = 0 4, and . different angular velocities Ω F , as indicated. The left panel displays the velocity profiles and the right panel the corresponding power profiles.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0126v1-figure-70237884-6139-4d2e-9b3b-f823511678c0.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0126v1-figure-880bde26-a039-42f9-856c-99323b8ba754.jpg)

Fig. 7.- A plot of the function f (˜) defined in Equation (22). a

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0126v1-figure-97807742-2855-48d5-b2a4-ec151a9a895e.jpg)

Fig. 8.- Effect of the loading efficiency on the position of the Alfv´n surfaces. The profile of e ˜ Ω L F is shown for 4 solutions corresponding to different values of κ , and ˜ = 0 95, Ω a . F = ω H / 2. The range of Alfv´n radii e that allows rotational energy extraction is delineated by the shaded area. In the force-free case ( κ → 0, ˙ E H = P BZ ) the Alfv´n points coincide with the light surfaces. e This result is in accord with that derived by Takahashi et al. (1990) for an adiabatic flow.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0126v1-figure-d1c0ecea-d7c9-4aac-bd87-b66d4176911e.jpg)
Noemie Globus Amir Levinson
Jet formation in GRBs: A semi-analytic model of MHD flow in Kerr geometry with realistic plasma injection
We construct a semi-analytic model for MHD flows in Kerr geometry, that
incorporates energy loading via neutrino annihilation on magnetic field lines
threading the horizon. We compute the double-flow structure for a wide range of
energy injection rates, and identify the different operation regimes. At low
injection rates the outflow is powered by the spinning black hole via the
Blandford-Znajek mechanism, whereas at high injection rates it is driven by the
pressure of the plasma deposited on magnetic field lines. In the intermediate
regime both processes contribute to the outflow formation. The parameter that
quantifies the load is the ratio of the net power injected below the stagnation
radius and the maximum power that can be extracted magnetically from the black
hole.