Force-Free Black-Hole Magnetospheres

A rotating black hole can be surrounded by an electromagnetic field whose stress greatly exceeds the inertia of the charges carrying its current. In that limit the plasma does not disappear: it supplies charge and current while experiencing essentially no net Lorentz force. The resulting force-free model is self-consistent only where the field is degenerate and magnetically dominated, its current closes globally, its field lines cross every light surface regularly, and the field is regular for an observer falling through the future horizon. These conditions turn a seemingly arbitrary electromagnetic configuration into a constrained magnetosphere and identify precisely where kinetic plasma physics must replace it.

Required background. Classical Field Theory supplies Maxwell theory, stress-energy conservation, and covariant field invariants.

Helpful background. Black-Hole Geometry and Mechanics supplies Kerr Killing fields and horizon regularity; Relativistic Matter and Compact Objects supplies relativistic plasma and equation-of-state context; and Fluid and Kinetic Theory supplies the continuum and kinetic limits behind charge supply and current closure.

Force-free electrodynamics and its domain

Let $F_{ab}$ be the electromagnetic two-form and $j^a$ the four-current. Maxwell’s equations are

$$\nabla_b F^{ab}=j^a, \qquad \nabla_{[a}F_{bc]}=0,$$

with electromagnetic units absorbed into the definition of $j^a$. The divergence of the electromagnetic stress tensor is the negative Lorentz-force density. If matter stress is negligible at leading order, conservation therefore requires

$$F_{ab}j^b=0.$$

This is the force-free equation, not the statement $j^a=0$. Indeed, a black-hole magnetosphere needs currents to produce its toroidal magnetic field and Poynting flux. The covariant formulation and the relation between force freedom, degeneracy, and magnetic domination are developed in Gralla and Jacobson 2014, §3, pp. 2502–2505.

For any local observer with four-velocity $u^a$, define electric and magnetic fields by $E_a=F_{ab}u^b$ and $B_a={}^*F_{ab}u^b$. A degenerate field obeys

$$F_{ab}{}^*F^{ab}=0 \quad\Longleftrightarrow\quad \boldsymbol E\mathbin{\cdot}\boldsymbol B=0.$$

Magnetic domination adds the invariant inequality

$$\frac12 F_{ab}F^{ab}=B^2-E^2>0.$$

Together these statements guarantee a local subluminal frame in which the electric field vanishes. In a frame where $\boldsymbol E\perp\boldsymbol B$, its velocity is the field-line drift

$$\boldsymbol v_D=\frac{\boldsymbol E\times\boldsymbol B}{B^2}, \qquad v_D^2=\frac{E^2}{B^2}<1.$$

The inequality is consequential. At $B^2-E^2=0$ the drift becomes null; if $E^2>B^2$, no subluminal magnetic frame exists and an ideal force-free evolution cannot represent the local dissipation or particle acceleration. Large magnetization, roughly $\sigma\sim B^2/(\rho h)\gg1$, motivates neglecting inertia, but it does not by itself guarantee $\boldsymbol E\cdot\boldsymbol B=0$, adequate pair supply, or a valid solution in a current sheet.

Flux surfaces, currents, and light surfaces

For a stationary, axisymmetric, degenerate magnetosphere, magnetic surfaces can be labelled by a poloidal flux function $\psi(r,\theta)$. With the force-free condition as well as symmetry and the homogeneous Maxwell equation, two quantities are constant on each connected flux surface:

• the field-line angular velocity $\Omega_F(\psi)$; and
• the polar current $I(\psi)$, which also fixes the toroidal magnetic field.

These are functions to be determined, not freely chosen constants. The reduction of the stationary axisymmetric field to $\psi$, $\Omega_F(\psi)$, and $I(\psi)$ is given in Gralla and Jacobson 2014, § 7.2. Force freedom is needed for $I$ to be a function of $\psi$; symmetry and degeneracy alone are insufficient. The current leaving along one set of flux surfaces must return along others or through current sheets. A local solution that never specifies this return path has not established a globally closed magnetosphere.

Let $t^a$ and $\phi^a$ be the stationary and axial Killing fields. A magnetic surface rotates with the helical vector

$$\chi_F^a=t^a+\Omega_F\phi^a.$$

A light surface occurs where this vector is null,

$$g_{tt}+2\Omega_F g_{t\phi}+\Omega_F^2g_{\phi\phi}=0.$$

It is the curved-spacetime generalization of the light cylinder. For prescribed $I(\psi)$ and $\Omega_F(\psi)$, the stationary stream equation is elliptic away from the light surfaces but degenerates there: its second-derivative coefficient vanishes. Smooth crossing imposes a first-order regularity relation among $\psi$, its transverse derivative, $I$, and $\Omega_F$. This is not an elliptic-to-hyperbolic transition of that stationary equation. A rotating black-hole magnetosphere generally has inner and outer light surfaces, and satisfying one does not guarantee regularity at the other Gralla and Jacobson 2014, §§ 7.4.2 and 9.3.

The future horizon is a regularity surface

The event horizon is generated by

$$\chi_H^a=t^a+\Omega_H\phi^a,$$

where $\Omega_H$ is the horizon angular velocity. Boyer–Lindquist components can look singular at the horizon even when the electromagnetic tensor is smooth in horizon-penetrating coordinates. Requiring finite fields for an ingoing observer gives the Znajek relation: the horizon value of $I$ is tied to $(\Omega_F-\Omega_H)\partial_\theta\psi$ times a positive metric factor. Its exact prefactor depends on whether $\psi$ denotes total, hemispheric, or $2\pi$-normalized flux, so a flux convention must accompany any displayed component formula.

This condition is often described informally as a horizon boundary condition. Causally and mathematically it is better understood as future-horizon regularity. The black hole is not a conducting material surface that prescribes an outgoing signal. A stationary exterior solution must instead be compatible with smooth ingoing fields. The Kerr regularity analysis and its relation to energy extraction are given in Gralla and Jacobson 2014, §9, pp. 2523–2526; the original construction appears in Blandford and Znajek 1977, §§2–3, pp. 435–440, PDF.

Two linked flux relations expose the physics. Along a stationary flux tube,

$$\frac{d\dot E}{d\psi}=\Omega_F\frac{d\dot J}{d\psi},$$

while horizon regularity makes the outward energy flux proportional to

$$\Omega_F(\Omega_H-\Omega_F)(\partial_\theta\psi)^2$$

times a nonnegative geometrical weight. Thus $0<\Omega_F<\Omega_H$ permits outward electromagnetic power. The next page derives this Blandford–Znajek extraction condition and its normalization.

Independent checks and failure boundaries

A proposed force-free magnetosphere should survive four checks. First, compute both invariants rather than inferring force freedom from a plotted field line: $\boldsymbol E\cdot\boldsymbol B$ must vanish and $B^2-E^2$ must remain positive in the claimed domain. Second, verify that $\Omega_F$ and $I$ are constant along each flux surface. Third, check finite, compatible crossings of every light surface and the future horizon in regular coordinates. Fourth, integrate the current over the full domain and identify the return path.

The limits are diagnostic. Setting $I=0$ removes the toroidal field and the stationary poloidal Poynting flux. Taking $\Omega_H\to0$ removes the rotational extraction window. Current sheets, charge-starved gaps, pair-creation zones, shocks, and heavily mass-loaded regions can violate the force-free assumptions even when the surrounding funnel is well described by them. Force-free electrodynamics then supplies outer boundary data for a kinetic or resistive calculation; it does not predict the particle distribution or emitted spectrum by itself.

Exercises

Magnetic-frame diagnostic. In a local orthonormal frame let $\boldsymbol E=(0.6,0,0)B_0$ and $\boldsymbol B=(0,0,1)B_0$. Decide whether a frame with $\boldsymbol E'=0$ exists, find its drift velocity and the magnetic-field magnitude in that frame, and repeat the existence test after replacing $0.6$ by $1.1$.

Solution

The fields are perpendicular, so $\boldsymbol E\cdot\boldsymbol B=0$. In the first case $B^2-E^2=(1-0.36)B_0^2=0.64B_0^2>0$. The magnetic frame therefore exists and moves with

$$\boldsymbol v_D=\frac{\boldsymbol E\times\boldsymbol B}{B^2}=(0,-0.6,0),$$

in units where the speed of light is one. Because $B^2-E^2$ is invariant and $E'=0$, its magnetic field has magnitude $B'=0.8B_0$. With $E=1.1B_0$, the invariant becomes $B^2-E^2=-0.21B_0^2$. The required drift would be superluminal, so no physical magnetic frame exists and the force-free model has crossed its validity boundary.

References

• Blandford, Roger D., and Roman L. Znajek. “Electromagnetic Extraction of Energy from Kerr Black Holes.” Monthly Notices of the Royal Astronomical Society 179 (1977): 433–456. DOI. Open PDF.
• Gralla, Samuel E., and Ted Jacobson. “Spacetime Approach to Force-Free Magnetospheres.” Monthly Notices of the Royal Astronomical Society 445 (2014): 2500–2534. DOI.