Launching and Quenching of Black Hole Relativistic Jets
At Low Accretion Rate
Abstract
Relativistic jets are launched from black hole (BH) Xray binaries and active galactic nuclei when the disk accretion rate is below a certain limit (i.e., when the ratio of the accretion rate to the Eddingtion accretion rate, , is below about ) but quenched when above. We propose a new paradigm to explain this observed coupling between the jet and the accretion disk by investigating the extraction of the rotational energy of a BH when it is surrounded by different types of accretion disk. At low accretion rates (e.g., when ), the accretion near the event horizon is quasispherical. The accreting plasmas fall onto the event horizon in a wide range of latitudes, breaking down the forcefree approximation near the horizon. To incorporate the plasma inertia effect, we consider the magnetohydrodynamical (MHD) extraction of the rotational energy from BHs by the accreting MHD fluid, as described by the MHD Penrose process. It is found that the energy extraction operates, and hence a relativistic jet is launched, preferentially when the accretion disk consists of an outer ShakuraSunyaev disk (SSD) and an inner advectiondominated accretion flow. When the entire accretion disk type changes into an SSD, the jet is quenched because the plasmas brings more restmass energy than what is extracted from the hole electromagnetically to stop the extraction. Several other observed BH diskjet couplings, such as why the radio luminosity increases with increasing Xray luminosity until the radio emission drops, are also explained.
Subject headings:
accretion, accretion disks — black hole physics — Galaxies: active — magnetic fields — MHD — Xrays: binaries1. Introduction
Black hole (BH) relativistic jets are observed from black hole Xray binaries (BHXBs) and active galactic nuclei(AGNs). A puzzling coupling between the jet and the accretion disk has been recognized. Observationally, the radio luminosity (which is likely due to the radiation from the relativistic jet or may also be from the semirelativistic disk wind) increases with increasing Xray luminosity (which is likely due to the radiation from the disk), until it shows a sudden drop when the Xray luminosity exceeds a certain limit (Gallo et al., 2003; Maccarone et al., 2003; Fender et al., 2004; Trump et al., 2011). In terms of the dimensionless accretion rate , this limit has a value about , where , is the accretion rate, is the Eddingtion accretion rate, is the Eddingtion luminosity, is the efficiency converting the energy of the accreting mass into the radiation energy, and is the speed of light. In addition, when the Xray luminosity changes, the Xray spectral state also changes, indicating that the accretion disk changes its type with the accretion rate (Fender et al., 2004; Trump et al., 2011; Esin et al., 1997; Remillard & McClintock, 2006). Since it is believed that BH relativistic jets are originated from the rotating BHs that are threaded by largescale magnetic lines and are powered by the rotational energy of the BH (Blandford & Znajek, 1977), an investigation of whether the BH energy can be extracted by holethreading magnetic field lines when the BH is surrounded by different type of accretion disks, may offers the key to an understanding of the diskjet coupling of BHXBs and AGNs.
How the disk type varies is described below. When the accretion rate is extremely low, the energy transfer due to Coulomb collisions between ions and electrons is very inefficient. As a result, the ions cannot efficiently transfer the heat to the electrons, and the disk becomes an advectiondominated accretion flow or ADAF for short (Narayan & Yi, 1995; Narayan et al., 1997, 1998), which corresponds to a hot, radiatively inefficient, optically thin and geometrically thick disk. As the accretion rate increases, the outer part of the disk cools down due to efficient Coulomb collisions, and becomes a cold, optically thick and geometrically thin disk (Shakura & Sunyaev, 1973; Novikov & Thorne, 1973), hereafter we use the term ”SSD ” (ShakuraSunyaev disk) to refer to such disk; consequently, a combined disk is formed with an inner ADAF and an outer SSD. In a combined disk, the radius where the transition from thin disk to ADAF occurs is called the transition radius, . In general, decreases with increasing (Honma, 1996; Manmoto et al., 2000). When the accretion rate further increases and exceeds a critical value with , the Coulomb interaction becomes so efficient at every radius of the disk that the whole disk cools down, changing into the SSD. The plasma rotational motion is nearly Keplerian in a SSD, whereas it becomes subKeplerian in an ADAF, because the pressure gradient force contributed in the radial force balance in the latter type of accretion disk.
It is interesting to note that the accretion flow geometries near the horizon are quasispherical for all the above disks (see Section 2 for more details). As a result, there are considerable plasma flow along the holethreading line, reducing the electromagnetic extraction of BH energy by contributing their restmass energy. This consideration mainly differentiates our model from previous explanations of the diskjet couplings (e.g., Meier 2005; Ferreira et al. 2006). We emphasize that this effect is important at low accretion rate, e.g., when , due to the quasispherical geometry near the event horizon; at higher accretion rate , e.g., when as large as (Abramowicz et al., 2010), the accretion disk type further varies to a slim disk (Abramowicz et al., 1988, 2010; Sa̧dowski, 2009) and a disktype geometry near the horizon is realized. Therefore, a pure electromagnetic extraction of BH energy at most of the latitudes becomes a good approximation.
To consider the contribution of the quasispherical accreting plasma near the BH, we adopt the magnetohydrodynamics (MHD) theory and consider the “MHD Penrose process” (Takahashi et al., 1990), which includes both the electromagnetic and plasma inertia effects, to examine the extraction of the BH rotational energy^{1}^{1}1 In comparison, in the magnetically dominated limit, the BlandfordZnajek process (Blandford & Znajek, 1977) becomes a good approximation of the MHD Penrose process.. By constraining the magnetic field strength and the plasma density from the theory of accretion disks in a strong gravitational field, we investigate the importance of plasma inertia for different accretion disk types at low accretion rates. It is found that the MHD extraction of BH energy takes place preferentially in a combined disk and therefore enables the launch of a relativistic jet; however, such relativistic jets are quenched if a magnetic dominance breaks down when the entire accretion disk become a thin disk. We also propose a new paradigm of the coupling between the accretion disk and the jet for both BHXBs and AGNs. The observed launching and quenching of relativistic jets at low accretion rates and several other puzzling observational results, are naturally explained.
We introduce our model in Section 2. In Section 3, we summarize necessary formulae and concepts of the general relativistic ideal MHD flow in a stationary and axisymmetric BH magnetosphere. Next, the result is given in Section 4. Finally, discussion and conclusion are respectively presented in Section 5 and 6.
2. The Model
To qualitatively investigate whether a relativistic jet can be launched when the BH is surrounded by different type of accretion disks, we consider a model with following assumptions. Comments on the assumptions are provided in Section 5.3.
1. Accretion disks are described by the analytical solutions adopted from the accretion disk theories. The accretion disk type is important since both the largescale field strength and the mass flux per magnetic flux tube (hereafter, mass loading) are related to the disk. For simplicity, we adopt the analytical solutions of the ADAF and the SSD to represent the disk properties. For a combined disk, these two types of disk solutions are connected at the transition radius. The analytical solution we adopted is provided in the Appendix B. The middle region’ solution of a SSD, in which the pressure and the opacity is respectively dominated by the gas pressure and the electronscattering opacity, is used.
2. Parabolic largescale holethreading magnetic field lines are considered. We consider largescale magnetic field lines that are dragged in from a distant region by the accretion flow. It is expected that the field lines can be finally accumulated near the BH with a parabolic geometry. Thus, instead of solving for the field geometry, we adopt the paraboloidal field solution given in Blandford & Znajek (1977) for the magnetic flux function ,
(1) 
where the field strength is to be determined (see Appendix C for more details). Here, we focus on magnetic field lines that thread the horizon, since they are responsible for the extraction of the rotational energy from a rotating BH.
3. The accretion flow has a quasispherical geometry near the event horizon. An accretion flow will have a quasispherical geometry near the horizon if the flow become transonic before the radius of the innermost stable circular orbit, , (Abramowicz & Zurek, 1981). Compare to a disklike geometry, in which the accreting plasma enters the BH through a plunging region near the equatorial plane, the accreting plasma in a quasispherical geometry can fall onto the horizon at a wider range of latitudes. Since the sonic point of the accretion flow is well outside of for a typical ADAF, which has (Narayan et al., 1997), and the sonic point is close to for an SSD, the plasma distribution near the horizon should have a quasispherical geometry for all the disk types we considered here, ADAF, combined disk, and SSD. Note that, for the SSD, the plasmas are conventionally supposed to plunge into the horizon in a disklike geometry, which forms a striking contrast with our current picture.
4. The magnetic field strength is parameterized by the gravitational binding energy. Providing that the diffuse effect of the largescale magnetic field is relatively unimportant, largescale magnetic field lines can be ‘arrested’ by the accretion disk, being dragged from the outer region toward the central BH (Narayan et al., 2003; Spruit & Uzdensky, 2005; Rothstein & Lovelace, 2008). It is convenient to parameterize its strength at radius by (see Appendix A)
(2) 
where is the gravitational constant, is the mass of the central BH , is the surface density of the disk, and is the ratio of the gravitational binding energy of the disk at to the largescale magnetic field energy inside radius . The determination of is not yet clear; however, it can be related to the ionization degree of the accretion flow and the magnetic Prandtl number.
5. The accretion disk type is a function of the accretion rate. As described in Section 1, from low to high , the disk type varies as follows: (I) ADAF, (II) combined disk that consists of an inner ADAF and outer SSD, and (III) SSD. Hereafter, we call the corresponding range of as “Range I”, “Range II”, and “Range III”, respectively. Many studies on the transition radius between inner ADAF and outer SSD give similar results (Honma, 1996; Manmoto et al., 2000). To qualitatively describe how the disk configuration changes with , we adopt the analytical solution of from Honma (1996)
(3) 
where is the viscosity of ADAF in the prescription (Shakura & Sunyaev, 1973), is the critical accretion rate with a typical value , beyond which there is no ADAF solution and hence the entire disk becomes the SSD type. We recognize that the entire disk becomes thin disk if is met (Honma, 1996; Manmoto et al., 2000), and define the critical accretion rate, , such that is obtained at . By further assuming that the accretion disk has a finite outer radius , different ranges can be defined by (see also Figure 1) Range I: ; Range II: ; and Range III: .
Our strategy to investigate the formation of relativistic jet at different disk types is described below. In contrast to that there is only outflow along the largescale, diskthreading field lines for disk winds (Blandford & Payne, 1982; Sa̧dowski & Sikora, 2010), there are both inflow and outflow along the holethreading field lines. The inflow (or outflow) is launched when the gravitational force acting on the plasma is larger (or smaller) than the magnetocentrifugal forces do. Such a difference allows us to monitor the outflow properties by the inflow behavior; for an MHD flow along a holethreading line, a powerful outflow is realized only when the rotational energy is extracted by the inflow. Therefore, the condition for jet formation is to extract the rotational energy of BH, which can be examined by solving the equation of motion, namely, that the relativistic Bernoulli equation (BE), of the inflow along the field line (see Section 3 for more details).
The influences of the disk type lie in how they affect the largescale field strength () and the mass loading. At the injection point of the inflow, , by denoting the superscript “ADAF” and “SSD” as the analytical solutions of ADAF and SSD, respectively, and assuming that the same value of can be applied to both ADAF and SSD, these two physical properties can be estimated (see Table 1) in different region of the computation domain of accordingly.
Range I: ADAF type.— The largescale field at is estimated by the ADAF surface density, that is, by Equation (2). The mass flux per field line at is estimated by the ADAF solution as well, which gives .
Range II: Combined disk type.— The largescale magnetic field at , which contains an additional term due to the magnetic field advected from the thin disk to the ADAF, is estimated by . The mass flux per magnetic flux tube at is computed by , since the injection point is located inside the ADAF.
Range III: Thin disk type.— The largescale field at is estimated by and the mass flux per field line at is estimated by .
Range  Disk Type  Field Strength^{a}^{a}footnotemark: at the Injection Point  Mass Loading^{b}^{b}footnotemark: at the Injection Point 

I  ADAF  
II  ADAF+SSD (Combined)  
III  SSD 
See Equation (2).
The mass flux per magnetic flux tube.
In next section, we summarize necessary concepts of the ideal MHD flow in a BH magnetosphere, including the relativistic BE and the MHD Penrose process. The result of our model is provided in Section 4.
3. MHD Flows around a Rotating Black Hole
Stationary and axisymmetric ideal MHD flows in a rotating BH magnetosphere have been investigated in several studies (Camenzind, 1986a, b, 1987; Takahashi et al., 1990; Hirotani et al., 1992). In this paper, we adopt the Kerr metric as the background geometry with signature (,+,+,+) and use the geometrized units such that , that is, the length scale is in the unit of ‘’. Most of the derivations in this section follows the calculations in Takahashi et al. (1990). Note that a different signature (+,,,) is used in their paper.
3.1. The Relativistic Bernoulli Equation
In the BoyerLindquist coordinate, the metric is
(4) 
where , , and . The ideal MHD condition requires that the electric field vanish in the fluid’s rest frame,
(5) 
where is the electromagnetic field tensor satisfying the Maxwell’s equation, is the four velocity of the fluid, is the electromagnetic vector potential, and the comma refers to the derivative.
The fluid equation of motion can be split into the two components: the ‘relativistic BE’ (the poloidal equation) that describes the flow along the field line, and the ‘GradShafronov equation’ (the transfield equation) that describes the force balance perpendicular to the field lines, where the semicolon refers to the covariant derivative. The BE can be obtained by projecting the equation of the motion along the field line; however, it can alternatively derived from the definition of proper time,
(6) 
The energymomentum tensor consists of two terms. One is the electromagnetic term,
(7) 
and the other is the plasma term,
(8) 
where is the pressure and is the energy density. The proper number density obeys the continuity equation
(9) 
and the relativistic specific enthalpy satisfies
(10) 
where denotes the rest mass energy of the proton, is the adiabatic index, and the subscript “” is the quantity evaluated at the particle injection point. The assumption of an adiabatic flow requires that the entropy along the field line be constant, which ensures that the term remains constant along the field line.
There are four more conserved quantities along the field line: the angular velocity of the field line, , the particle flux per unit flux tube, , the total energy of the flow per particle, , and the total angular momentum of the flow per particle, . They can defined by (cf. Equations(2.3)(2.6) of Hirotani et al. 1992)
(11) 
(12) 
(13) 
(14) 
where , denotes the fluid angular velocity, and the toroidal magnetic field seen by a distant static observer with four velocity . One should note that the sign of is defined by the signs of and . Although the magnetic field is defined by
(15) 
however, it is convenient to introduce the rescaled poloidal field
(16) 
where runs over the poloidal coordinates, and . Combining Equations (12)(14), we have
(17) 
where
(18) 
and the poloidal Alfven Mach number is
(19) 
The poloidal velocity of the plasma is defined by
(20) 
Substituting Equation (17) into Equations (13) and (14), we can express the energy of the fluid, , and the angular momentum of the fluid, , as
(21) 
(22) 
Finally, combining Equations (4), (20), (21), and (22), we obtain the BE (cf. Equation (17) of Takahashi et al. 1990) ,
(23) 
where and are defined as
(24) 
(25) 
The BE, Equation (23), can be rewritten as a polynomial of with the aid of Equations (10) and (19). Note that the order of the polynomial equation depends on the second term of Equation (10), because particle number conservation gives .
3.2. The Cold Limit
The cold limit can be adopted when the second term in Equation (10) is relatively unimportant to its first term, i.e., (in c.g.s unit). Equation (10) therefore reduces to and the BE becomes a fourthorder polynomial equation of . By introducing a parameter , we have
The cold BE, therefore, can be expressed in terms of
(26) 
where
Even for a disk whose temperature is comparable with the virial temperature and the internal energy is comparable to the released gravitational energy (e.g., ADAF), the cold limit can still be modestly satisfied. This is because the released gravitational binding energy is still a small fraction of the restmass energy. Therefore, we solve the BE in cold limit, Equation (26), for simplicity and expect the result is at least qualitatively correct. In order to specify the coefficients and , we need to specify the mass loading (mass flux per magnetic flux tube), , and the strength and geometry of the magnetic field. It is provided in the Appendix C how we compute these coefficients.
3.3. Light Surfaces and the Separation Point
There are two light surfaces, which are defined by , in a BH magnetosphere. The outer light surface is formed by the centrifugal force in the same manner as the light cylinder in a pulsar magnetosphere, and the inner light surface is, on the other hand, formed by the strong gravity of the BH and the rotation of the magnetic field.
If a plasma starts with a negligible poloidal velocity, the inflow and outflow along the field line separate at the point where the gravitational force balances with the magnetocentrifugal forces. This point is called the “separation point” (or the “stagnation point”) and defined by , which leads to , where the prime denotes the derivative along the flow line. That is, if a plasma is injected at with (or with ), it will be accelerated inward (or outward), subsequently passing through the Alfven point and the inner (or outer) light surface. Figure 3 shows the contours of around a BH with spin parameter . Note that, at the injection point, the BE reduces to
(27) 
3.4. Critical Points
The critical points for an MHD flow can be found by differentiating Equation (23) along the flow line, , where explicit expressions of and are given in Camenzind (1986b) and Takahashi et al. (1990). For a general MHD inflow, there are three critical points at which both and vanish: the ‘slowmagnetosonic point’, the ‘Alfven point’, and the ‘fastmagnetosonic point’; the poloidal velocity there matches the slowmagnetosonic velocity, the Alfven velocity, and the fastmagnetosonic velocity, respectively. In the cold limit, the slowmagnetosonic velocity reduces to zero. Due to the causality requirement at the event horizon, any MHD inflow must pass through the fastmagnetosonic point, after passing through the Alfven point.
At the Alfven point , the Mach number equals ,
(28) 
Note that the requirement of implies that must be located between these two light surfaces. Besides, it is convenient to express and in terms of the quantities evaluated at the Alfven point. Since both and automatically vanish at the Alfven point, no additional constraints are imposed on the MHD flow at the Alfven point. Therefore, by requiring that the numerators in Equations (21)(23) vanish at , we obtain
(29) 
(30) 
and
(31) 
Combining Equations (27), (29), and (30), we can express the conserved quantities and in terms of and .
After the MHD inflow pass through the Alfven point, vanishes again in at the fastmagnetosonic point; thus, should also vanish there. We should notice here that both and can be satisfied at the fastmagnetosonic point only for a specific combination of the conserved quantities, , , , and . In previous works, e.g., Camenzind (1986a) and Takahashi et al. (1990), the authors gave , , and (or equivalently, , , and ) and search to select the trajectory that pass through the fastmagnetosonic point. In this work, we instead give , , and , considering plasma injection from the innermost region of the accretion disk, and solve for (or equivalently, ). As explained next, whether the rotational energy of the BH is extracted outward can be determined by the location of .
3.5. Negative Energy Flow
The sign of can be solely determined by via Equation (29) because it follows form Equation (27) that holds for any MHD flow starting from a position between the two light surfaces with a vanishing poloidal velocity. Note that becomes negative if . Thus, near the BH, a spatial region called the “negative energy region”, which satisfies , can be defined (Takahashi et al., 1990). If the Alfven point resides in this region, the MHD inflow has (which also implies because of ; whereas does not imply ).
A spinning BH differs from a static BH by an additional structure called “ergosphere” (which is defined by ), a region inside which any particle must orbit in the same direction of the hole spin. Being locating outside the horizon, the ergosphere can be viewed as the region that stores the rotation energy of a BH. Note that the negative energy region exists only when , and it is always located inside the ergosphere (see also Figure 3). It is also noteworthy that a negative energy MHD inflow can be launched from a region outside the ergosphere.
3.6. Outward Energy Flux and the MHD Penrose Process
Inflow Property  Infow Energy^{b}^{b}footnotemark:  Inflow Energy Flux^{c}^{c}footnotemark:  Jet Launch?^{d}^{d}footnotemark: 

plasma dominated ()^{a}^{a}footnotemark:  NO  
magnetically dominated ()  YES 
and are usually satisfied. See the text for more details.
The total energy , Equation (13), contains both the electromagnetic and the plasma contribution.
, Equation (32), where holds for an inflow.
A relativistic jet, which has , , and hence , can be launched only when it connects to an inflow with
The outward energy flux of the flow can be split into the electromagnetic part and the fluid part ,
(32) 
where,
(33) 
(34) 
Therefore, a negative energy inflow ( and ) results in a positive outward energy flux (), that is, the BH rotational energy is extracted. Note that usually holds because its sufficient condition, , is satisfied under normal conditions (MacDonald & Thorne, 1982), where refer to the rotational angular velocities of the BH. Note also that the plasmas usually cannot extract BH’s energy (i.e., usually holds) because of their restmass contribution, unless they are in a negative energy orbit, which is rare.
In a magnetically dominated magnetosphere ( ), because is continuous across the separation point, a stationary outflow solution (, and ) is possible only when it connects a negativeenergy inflow ( and , such that ). Consequently, a stationary MHD jet from a BH ergosphere ( a “ergospheric jet”) can be launched in a magnetically dominated magnetosphere only when the MHD Penrose process, , operates. If a magnetic dominance breaks down, can be discontinuous across the separation point, that is, an outflow () can connect to a positiveenergy inflow (), owing to the energy supplied by the accreting plasma near the separation point. However, this kind of outflow, similar to the disk wind, will not be accelerated into relativistic energies and therefore is not of interest in the present paper. To summarize, for an MHD flow along the holethreading field line, a stationary outflow with positive outward energy flux () should connect an inflow solution with at the separation point, provided that the injected energy flux by accretion is small compared to the Poynting flux. Thus, we can investigate the activity of a Poyntingfluxdominated jet (or an outflow) by examining the inflow energy, (see also Table 2).
It should be noted here that the term “MHD Penrose process” refers to the case of (remember is the ‘total’ energy, consisting both the electromagnetic plus plasma energy of the flow), as initially proposed by Takahashi et al. (1990); however, it was used to indicate (i.e., a negative energy orbit of the particle) in some previous studies (Hirotani et al., 1992; Koide et al., 2002; Semenov et al., 2004; Komissarov, 2005). The latter case is possible, but can only be transient (Koide et al., 2002; Komissarov, 2005). In comparison, in the original idea of the ‘Penrose process’ (Penrose, 1969), only the particle contribution () was considered; whereas in the BlandfordZnajek process (Blandford & Znajek, 1977), only the electromagnetic contribution () was considered.
4. Result
We choose representative parameters to demonstrate the jet behavior at different disk type qualitatively. Because the extraction of the rotational energy of BH is mainly related to whether the BH magnetosphere is magnetically dominated or not, the change of the parameters will result in only minor modification of the result in our model.
We adopt the dimensionless BH spin parameter to be and the angular velocity of the field, , to be half of the angular velocity of the hole, (i.e., ). The MHD inflow along the field line that intersects the event horizon at a modest latitude (, i.e., about from the pole; see Figures 2 and 3 for this field line) is considered, and the inflow is assumed to be launched at the separation point (that is, the injection point is located at the separation point). In addition, to investigate the importance of the BH mass, we consider separately the BHXB and AGN cases. For BHXB, we assume , whereas for AGN, . Furthermore, to explore the dependence of the ability of the disk in storing the largescale magnetic field, we adopt different values of (, and ).
After starting from the injection point, the accreting MHD plasma is accelerated inward, passing thorough the Alfven point, , and the fastmagnetosonic point, and finally enter the BH. The inflow energy can be computed by imposing the fastmagnetosonic condition once the injection point, the angular velocity of the field and the mass loading are specified. However, to discuss the conditions for jet launching () and quenching (), it is adequate to examine the sign of , which is uniquely specified by the location of . Specifically, if resides in the negative energy region, the MHD inflow has a negative energy () and the MHD Penrose process operates. Therefore, by solving from the BE for an MHD inflow onto a BH, we can investigate whether the energy extraction takes place as a function of . Detailed computation steps are provided in the Appendix.
Let us first consider the BHXB case with the BH mass, (Figures 4(a)(d)). Figure 4(a) shows how the largescale field strength at the injection point varies with the accretion rate. It is important to note that in a combined disk (Region II), the surface density of the disk does not monotonically increase inward but has a local maximum at the transition radius, . Since an SSD can confine more largescale field than an ADAF because of its greater surface density for a fixed , the largescale field arrested in the inner ADAF for a combined disk (solid line) is relatively stronger than what the ADAF alone would bring into the BH (dashed line). The locations of for four discrete values of are shown in Figures 4(b) and 4(c); the solid, dashed, dashdotted, and dashdotdot lines correspond to , , , and , respectively. It follows from the zoomin figure (Figure 4(c)) that enters the negative energy region (indicated by the shaded region) in Range II for , , and , indicating that the MHD Penrose process turns on in a combined disk. This is because, near the horizon, the relatively strong magnetic field (which is mainly provided by the outer SSD) and the relatively small mass loading (which is provided by the inner ADAF) result in a magnetically dominated environment. In comparison, a magnetically dominated environment near the horizon is relatively harder to achieved when the disk is an ADAF (Range I) or an SSD (Range III), because it is the same type of disk responsible for both the confinement of the largescale field and the loading of the plasma. Such diskjet coupling provides the reason why observed relativistic jet is launched when is near (but smaller) the value , but quenched when exceeds . In addition, since increases with increasing , increases negatively (i.e., the rotational energy of BH is more efficiently extracted outwards) with increasing once enters the negative energy region in Range II (Figure 4(d)). This is because the BH magnetosphere becomes more magnetically dominated when the transition radius decreases with , resulting from that the difference between the surface density of the outer SSD and inner ADAF becomes more significant.
Next, let us consider the AGN case with (Figures 4(e)(h)). Although the solutions look qualitatively similar, there are three major differences from the BHXB case. First, comparing Figures 4(b), (c), (f) and (g), we find that for an AGN, is realized in a wider range of for a fixed or in a wider range of for a fixed . Second, comparing Figures 4(d) and (h), we find that an AGN can extract more BH rotational energy per plasma particle than a BHXB can. Third, in an AGN, the MHD Penrose process works also in an SSD (i.e., Range III) for a relatively greater (e.g., ). However, the jet power decreases with increasing in an SSD (i.e., Range III), in contrast to the case of a combined disk (i.e., Range II). The above comparison can help to understand the observed properties of accreting BH in different masses (see below).
5. Discussion
5.1. Universal Paradigm of Black Hole DiskJet Coupling at Low Accretion Rate
The results in Section 4 can offer the key to an understanding of the observed connections between the disks and the jets for BHXBs and AGNs. A universal paradigm of diskjet couplings on the plane (c.f. Figures 4d & 4h) is diagrammatically presented in Figure 5. Transitions of diskjet states are indicated by the arrows for three representative values of . Since the terminal speed of an MHD outflow is roughly determined by the magnetic dominance at the launching point, the more magnetically dominated the disk innermost region becomes (i.e., the more a system moves downward on the plane), the greater speed an MHD outflow achieves. This discussion is valid for both the disk wind near the BH and the ergospheric jet. The two downward arrows in the right edge of this figure indicate how the terminal speed changes for a disk wind and a relativistic jet. Note that the ergospheric jet can only be launched when (i.e., in the lower half of the diagram).
The transitions of diskjet states for BHXBs, which are indicated by the solid arrows, explain the observed features of the BHXB diskjet couplings. Fender et al. (2004) show that two type of outflows associated with the Xray spectral types are required to explain the systematic properties of the radio emission of BHXBs. The first type is a steady outflow observed when the BHXB is in a hard state (HS), while the second type is a transient, ballistic, powerful and relativistic jet when the source is in a steep powerlaw (SPL) state and being switching from the HS to the soft state (SS). We interpret that the steady outflow (first type) is associated with the disk wind and the transient relativistic jet (second type) is associated with the ergospheric jet ejected for a modest or small . Note that, as shown in the paradigm, when the entire disk become an SSD, not only that the ergospheric jet is turned off but also that the disk wind weakens, which explains the sudden drop of observed radio luminosity when the state transits to the SS. In addition, the energy of the jet per particle, , as well as the ergospheric jet power, , increase with increasing in a combined disk, where the surface integral is carried out on a closed surface (e.g., the horizon); this likely corresponds to the observational fact that both the Lorentz factor of a jet and its radio luminosity rapidly increase with increasing before the jet is quenched. When the BH is surrounded by a combined disk and the MHD Penrose process takes place, the relatively strong field for the inner ADAF can result in a configuration which has been studied in previous simulations (Narayan et al., 2003; Igumenshchev et al., 2003; Igumenshchev, 2008; Punsly et al., 2009; McKinney et al., 2012), here we use the name magnetcially choked accretion flow (MCAF)^{2}^{2}2An MCAF is formed when the accumulated magnetic field near the BH reaches the equipartition level. The strong magnetic field is able to support (or disrupt) the disk at certain radius , where . Although slowed down by the magnetospheric barrier near , accreting plasma can still accrete inward in the form of highly non axisymmetry, irregular streams due to magnetic RayleighTaylor instability (Igumenshchev et al., 2003; Punsly et al., 2009; McKinney et al., 2012). We stress that the magnetically dominated magnetosphere developed when the BH surrounded by a combined disk can result in an MCAF or other conceptually similar accretion flow. according to McKinney et al. (2012). Accompanying by the jet formation, the formation of an MCAF can also explain some observation features when the source is in SPL. For an MCAF, because the slowed down flow can earn more time to release their heat energy via radiation, the efficiency for the flow to convert mass to radiation energy increases (Narayan et al., 2003). This can be the reason why the disk luminosity in SPL state is higher than that in HS or SS (Fender et al., 2004). In addition, the highfrequency quasiperiodic oscillation (QPO) feature observed at SPL state may be related to the jetdisk QPO mechanisms found by McKinney et al. (2012).
For AGNs, the transitions are indicated by the dashed arrows, which imply the following diskjet couplings. First, an efficient MHD Penrose process (i.e., a greater ) is viable at a relatively smaller in an AGN compared to a BHXB. Since the SSD exists only in the outer region when is small, it is consistent with the observational fact that radioloud AGNs favor a spectrum without a UV bump (Kawakatu et al., 2009), which possibly arises in the inner edge of an SSD. Second, the MHD Penrose process works more efficiently with increasing BH mass. This conclusion explains why AGN jets show higher Lorentz factors than BHXB jets. Moreover, it is also consistent with the observational fact (Liu et al., 2006; Chiaberge & Marconi, 2011) that most radioloud AGN host BHs with larger masses (), while radioquiet AGNs host BHs with no obvious mass constraints. Third, although a BHXB jet disappears when the disk transits from a combined disk to a SSD, an AGN jet could still be launched from the ergosphere, because the SSD in an AGN tends to be more magnetically dominated than in a BHXB owing to a smaller mass loading, , around a more massive BH. We interpret that these ergospheric jets ejected from the inner edge of the thin disks in AGNs result in (at least part of) the radioloud quasistellar objects; this forms a striking contrast to BHXBs from which no strong or steady radio emission has been detected in their SS. Estimated from the dynamical timescale, the characteristic time scale for an accreting BH system to sustain the relativistic jet, , is scaled by the central BH mass with the relation . Therefore, with equals to days for the relativistic jet from a BHXB with , relativistic jets from an AGN with may last for years.
Note also that the relation between the accretion rate, , and the ergospheric jet power, , inverses in Ranges II and III. That is, increases (or decreases) with increasing , when is relatively low (or high) and the disk is a combined disk (or an SSD). Provided that the Xray luminosity indicates the disk accretion rate, these conclusions naturally explain the empirical relation for both BHXBs and AGNs: the radio luminosity increases with increasing Xray luminosity when the bolometric luminosity, , is small, i.e., when (Merloni et al., 2003; Falcke et al., 2004) , while the radio luminosity decreases with increasing Xray luminosity when some BH systems are luminous, i.e., when (King et al., 2011) .
5.2. The Estimated Speed of Ergospheric Jet
At large distances, the spacetime is described by the Minkowski metric; thus, the Lorentz factor of the jet, , is defined by . Because is conserved quantitatively along the field lines and because the second term of Equation (13) becomes negligibly small compared to the first term at large distances, we obtain
(35) 
where denotes the total energy of the MHD outflow.
We can further estimate the terminal jet speed, by linking with the total energy of the MHD inflow, , as:
(36) 
where the subscript ”in” and ”out”, respectively denotes the quantities for the inflow and the outflow. The above equation is obtained by considering the conservation of the outward energy flux near the separation point, provided that the magnetic field configuration does not change significantly there. In general, the ratio of the inflow and outflow particle flux, , is greater than unity because only a small portion of the plasmas will escape as an outflow. As a result, Equations (35) and (36) give
(37) 
Hence, the inverse values of the total energy shown in Figures 4(d) and (h), as shown in Figure 6, give an estimation of the lower limit of the ergospheric jet Lorentz factor for BHXBs and AGNs. The BHXB and AGN relativistic jet speeds inferred from the observations, which has the typical value for BHXBs (e.g., Table 1 of Fender et al. 2004) and for AGNs, can be consistently explained by our model.
5.3. Validity of the Model
To investigate the reason of the launching and quenching of BH relativistic jets when varies, we solve for stationary negative energy inflow solutions (which corresponds to the formation of relativistic jets) by adopting analytical solution of the ADAF and the SSD, assuming a spherical accretion geometry near the horizon, and estimating the field strength of the largescale holethreading magnetic field by the parameter . The major drawback of our approach is that we ignore the possible disk instability when the SSD becomes radiation pressure dominated. In addition, our axialsymmetry and stationary solution miss the temporal and the nonaxisymmetry properties of the jet formation. Also, the disk dynamo effect and the magnetic reconnection process can further modify the largescale fields, although their effects can be in general included in the parameter .
For our purpose, we ignore several aspects which we think are minor factors in determining whether the jet is on or off for simplicity. Several discussions on the assumptions and the robustness of the result are provided in the following.
We ignore the BlandfordZnajek extraction in the forcefree region near the pole of the horizon in a quasispherical accretion flow. However, such extraction is expected to be a relatively minor factor to determine the jet launching. The reason is that most of the holethreading field lines have considerable plasma loading on them, and that the power for magnetic fields extracting the energy from a rotating BH is proportional to (Blandford & Znajek, 1977).
The most important feature of a combined disk is that the significant surface density difference of the SSD and the ADAF at the transition radius. Although here we describe the combined disk solution by patching the analytical solutions of the SSD and the ADAF at the transition radius, the above feature is expected to be true even if a more detailed disk solution is considered.
There are uncertainties in determining the largescale magnetic fields. For example, the parameter can vary with the radius , and the total flux of the largescale diskthreading magnetic fields is not unclear. However, the field lines that threads the disk are expected to be redistributed when the accretion disk types vary. The diskjet coupling features described by our model can be qualitatively preserved as long as the diffuse effect is relatively unimportant than the advection effect for the largescale diskthreading field lines. Especially, a magnetically dominated magnetosphere near the horizon is still most likely realized when the BH is surrounded by a combined disk, provided that the largescale magnetic field advected inward from the outer SSD is large enough for the inner ADAF. Such requirement is not difficult to be satisfied as the transition radius gradually decreases to a small enough radius (Figures 4(a) and (e)).
The result presented in Figure 4 shows that a negative energy flow may form when the BH is surrounded by a combined disk (Range II). The conclusion can actually still be valid when more detailed pictures are considered. Below we discuss two possible concerns. First, we consider whether the result would change when numerical ADAF solutions, e.g., Narayan et al. (1997) and Gammie & Popham (1998), instead of the selfsimilar solutions, Equations (B1)(B3), are used. It is true that the selfsimilar analytical solutions are not a good approximation near the hole. For example, the causality requires a numerical solution of be greater than Equation (B1) near the BH, while the surface density be less than the selfsimilar solution, Equation (B3). Nevertheless, the numerator in Equation (C2) will not change because of the mass conservation, , provided that the disk height remains unchanged. Since the magnetic field is predominantly supported by the outer thin disk, is essentially given by . As a result, neither nor contains significant errors if we adopt the selfsimilar solution of ADAF. Thus, the discussion given in this paper is not vulnerable for details of the ADAF solutions. Second, we discuss whether the estimation of the mass loading is fair when . In this case, an MCAF is formed. The highly nonaxisymmetry, irregular MCAF will try to avoid the strong field ranges by the magnetic RayleighTaylor instability and result in strongfield, lowdensity ‘magnetic islands’ (Punsly et al., 2009). As a result, both and decrease in a magnetic island. Thus, even when , the mass loading from the magnetic island located at is lower than the estimate with the selfsimilar solutions, Equations (B1) and (B2). This means that the MHD inflow launched from the magnetic island should be even more magnetically dominated and that the energy of the inflow is still negative.
6. Conclusions
In this paper we focus on the fundamental question: Why relativistic jets from accreting BHs can be selectively launched under a certain accretion condition? Here we solve this problem by incorporating the plasma inertia effect, which have been neglected in previous studies, in conjunction with the electromagnetic effect on the extraction of BH energy. At low accretion rate (e.g., ), the accretion flow has a quasispherical geometry near the BH event horizon. As a result, whether the rotational energy of the BH can be effectively extracted outward is determined by the inward restmass energy of the accreting plasmas. When the outward electromagnetic energy flux dominates, a relativistic, ergospheric jet can be therefore launched; on the contrary, when the inward restmass energy dominates, there is no powerful relativistic jet. By examining the MHD Penrose process for BHs surrounded by different types of accretion disk at low accretion rates, we conclude that the MHD inflow becomes magnetically dominated preferentially in a combined disk, which consists of an outer SSD and an inner ADAF, to enable the launch of an ergospheric jet. Such relativistic jets are quenched if a magnetic dominance breaks down when the entire disk becomes an SSD. We also provide a universal paradigm of the coupling between the BH accretion disks and their relativistic jets at low accretion rates. Several observed diskjet connections are naturally explained. Future study on the relation between the accretion geometry near the event horizon and the extraction of the BH energy will be important to the understanding of the BH jet formation mechanisms.
Appendix A Parameterized the LargeScale Magnetic Field Strength
The largescale field stored in a disk can be estimated from the point of view of energies. The largescale field energy inside the volume enclosed by a radius should not exceeds the gravitational binding energy of the disk near , otherwise the field would escape by buoyancy. By considering and , we define a ratio between this two quantities to obtain
where is the gravitational constant, is the mass of the central BH, and is the surface density of the disk.
Alternatively, we can estimate the largescale field stored in a disk by a force balance. If a disk accretion is supported by a largescale field at , as described in Narayan et al. (2003), the largescale field can be estimated by the force balance to give where is assumed. Similarly, a parameter can be defined by the ratio between the gravity force and the field force; we thus obtain
Since the two estimations above are consistent in order of magnitude, we simply parameterize in terms of , , and , such that
(A1) 
Appendix B Analytical Solution of the Accretion Disks
With the scaled units,
we adopt the following radial velocity , plasma density , and surface density solutions of an ADAF and a thin disk. Note that what we need is (because ) to solve the BE, and that can be more or less correctly evaluated even in the innermost region by an ADAF selfsimilar solution without invoking a numerical solution. We therefore adopt ADAF selfsimilar solution (Narayan & Yi, 1995; Narayan et al., 1998) and obtain
(B1) 
(B2) 
(B4) 
(B5) 
(B6) 
where
and is the radius of the innermost stable circular orbit. Note the dimensionless coordinate, , is related to by . The parameter describes the viscosity in the prescription (Shakura & Sunyaev, 1973). The superscripts “ADAF” and “SSD” denote the Advectiondominated accretion flow and the ShakuraSunyaev disk analytical solutions, respectively. To elucidate the general feature, we adopt typical values of the viscosity parameter such that and .
Appendix C Solving the Relativistic Bernoulli Equation for the Location of the Alfven Point
Once the field strength, field geometry, and the mass loading per field line at different are specified, we can calculate all the coefficients in the cold BE, Equation (26), and determine the location of the Alfven point, .
Since Equation (26) is solved in the geometrized unit (in length unit of ), extra transformations of units from c.g.s. unit are necessary, namely,
(C1) 
(C2) 
These unit transformations can be done by two steps. First, we convert the unit from the c.g.s into the geometrized unit (Table 3). Second, we rescale the length unit in geometrized unit from “” into “” (Table 4).
Next, we find the location of the Alfven point at various accretion rate by the following steps:

Calculate the coefficients and in Equation (26).

Guess to compute and by Equations (29) and (30).

Using and , calculate the coefficients and in Equation (26).

Using the calculated , , , , solve the cold BE equation(26).

Find so that the fluid can smoothly pass through the fastmagnetosonic point.
A correct will give a physical inflow solution that passes through the fast point before reaching the horizon with a nonzero speed (Figure 7). If is incorrect, the resulting inflow solution becomes unphysical, being either subfastmagnetosonic (Figure 8), or discontinuous (Figure 9).
Quantity  c.g.s Unit  Geometrized Unit  Dimension of the Geometrized Unit 

Restmass energy 