Skip to content

The Blandford–Znajek Process

Scrollable formula. Use Left and Right arrow keys to reveal hidden content.

Scrollable table. Use Left and Right arrow keys to reveal additional columns.

Scrollable figure. Use Left and Right arrow keys to inspect hidden content.

The Blandford–Znajek process converts a rotating black hole’s angular momentum into outward electromagnetic power. Three ingredients must be kept distinct: magnetic flux threads the future horizon, field lines rotate at an angular velocity set by the global load, and the horizon field is regular in ingoing coordinates. Their combination makes the sign of the power transparent: for an aligned system, extraction occurs when 0<ΩF<ΩH0<\Omega_F<\Omega_H. The familiar scaling PBZ∝ΦBH2ΩH2P_{\rm BZ}\propto\Phi_{\rm BH}^2\Omega_H^2 follows only after a magnetic geometry and flux convention are declared. It describes horizon-scale extraction, not by itself the delivery, acceleration, collimation, or radiation of an observed jet.

Required background. Classical Field Theory supplies electromagnetic stress-energy, Killing-energy currents, and flux conservation.

Helpful background. Black-Hole Geometry and Mechanics supplies Kerr horizons and rotational energy; Relativistic Matter and Compact Objects supplies mass loading and relativistic plasma context; and Fluid and Kinetic Theory supplies the continuum limits needed to connect a force-free funnel to matter.

With c=1c=1, GG explicit, and dimensionless spin a∗=J/(GM2)a_*=J/(GM^2), the Kerr horizon radius and angular velocity are

rH=GM(1+1−a∗2),ΩH=a∗2rH.r_H=GM\left(1+\sqrt{1-a_*^2}\right), \qquad \Omega_H=\frac{a_*}{2r_H}.

The horizon generator is χHa=ta+ΩHϕa\chi_H^a=t^a+\Omega_H\phi^a. A stationary axisymmetric degenerate electromagnetic field is organized by a flux function ψ\psi; every connected flux surface has a field-line angular velocity ΩF(ψ)\Omega_F(\psi) and polar current I(ψ)I(\psi). The disk, remote load, current sheets, and plasma supply determine those functions globally. The horizon does not simply force ΩF=ΩH\Omega_F=\Omega_H.

Regularity at the future horizon ties II to (ΩF−ΩH)∂θψ(\Omega_F-\Omega_H)\partial_\theta\psi. Substituting that Znajek relation into the Killing-energy flux gives, on each horizon latitude,

dPoutdθ=W(θ) ΩF(ΩH−ΩF)(∂θψ)2,W(θ)>0,\frac{dP_{\rm out}}{d\theta} =\mathcal W(\theta)\, \Omega_F(\Omega_H-\Omega_F) \left(\partial_\theta\psi\right)^2, \qquad \mathcal W(\theta)>0,

where the positive geometrical weight W\mathcal W depends on the definition of ψ\psi and on the Kerr metric. This form is more robust than memorizing a convention-dependent prefactor. For an aligned field with ΩH>0\Omega_H>0, outward power requires

0<ΩF<ΩH.0<\Omega_F<\Omega_H.

The original slow-rotation construction and its horizon condition are given in Blandford and Znajek 1977, §§2–3, pp. 435–440, PDF. A coordinate-regular derivation and sign analysis appear in Gralla and Jacobson 2014, §9, pp. 2523–2526.

Stationarity relates the energy and angular-momentum fluxes carried outward along one magnetic surface:

dPout=ΩF dJ˙out.dP_{\rm out}=\Omega_F\,d\dot J_{\rm out}.

For a bundle with a single ΩF\Omega_F, its electromagnetic contributions to the black hole’s mass and angular-momentum rates are

M˙BHEM=−Pout,J˙BHEM=−PoutΩF.\dot M_{\rm BH}^{\rm EM}=-P_{\rm out}, \qquad \dot J_{\rm BH}^{\rm EM}=-\frac{P_{\rm out}}{\Omega_F}.

Here MBHM_{\rm BH} is gravitational mass, not the positive horizon rest-mass accretion rate M˙H\dot M_H used in disk budgets. Matter accretion supplies additional terms. In the quasistationary, subextremal limit the electromagnetic part of the first law gives

κH8πGA˙HEM=M˙BHEM−ΩHJ˙BHEM=Pout(ΩHΩF−1)>0\frac{\kappa_H}{8\pi G}\dot A_H^{\rm EM} =\dot M_{\rm BH}^{\rm EM}-\Omega_H\dot J_{\rm BH}^{\rm EM} =P_{\rm out}\left(\frac{\Omega_H}{\Omega_F}-1\right)>0

throughout the extraction interval. When ΩF\Omega_F differs between surfaces, integrate dPoutdP_{\rm out} and dPout/ΩF(ψ)dP_{\rm out}/\Omega_F(\psi) separately; a single unweighted angular velocity cannot replace that integral. Part of the rotational energy becomes outward Poynting flux; the remainder is irreversibly associated with horizon-area increase. At either endpoint, ΩF=0\Omega_F=0 or ΩF=ΩH\Omega_F=\Omega_H, the stationary power vanishes. The former can still exert a torque and increase area in the limiting expression; the latter has no slip between the field and horizon generator.

For a simple load, write x=ΩF/ΩHx=\Omega_F/\Omega_H. At fixed flux distribution the load-dependent factor is

4x(1−x),4x(1-x),

normalized to unity at its maximum. Differentiation gives x=1/2x=1/2. This is why the split-monopole solution has ΩF≃ΩH/2\Omega_F\simeq\Omega_H/2 at leading order Blandford and Znajek 1977, §6, pp. 447–449, PDF. It is an impedance-matched benchmark, not a universal law for every flux surface in every time-dependent magnetosphere.

Use Gaussian electromagnetic units, with c=1c=1, and define the unsigned flux through one effective horizon hemisphere by

ΦBH=12∫H∣Br∣−g dθ dϕ.\Phi_{\rm BH}=\frac12\int_H |B^r|\sqrt{-g}\,d\theta\,d\phi.

Here Br=∗FrtB^r={}^*F^{rt} is the coordinate magnetic-flux density in horizon-penetrating coordinates. With this convention, and a load close to the usual matched solutions, a widely used total-power form is

PBZ=κ4πΦBH2ΩH2f(xH),xH=GMΩH,P_{\rm BZ}=\frac{\kappa}{4\pi} \Phi_{\rm BH}^2\Omega_H^2 f(x_H), \qquad x_H=GM\Omega_H,

where κ\kappa encodes magnetic geometry, the argument xHx_H is dimensionless, and f→1f\to1 at low spin. Restoring units gives xH=rgΩH/cx_H=r_g\Omega_H/c and PBZ=κΦBH2ΩH2f(xH)/(4πc)P_{\rm BZ}=\kappa\Phi_{\rm BH}^2\Omega_H^2f(x_H)/(4\pi c). Under this convention, κ=0.053\kappa=0.053 for the split monopole and 0.0440.044 for a parabolic field Tchekhovskoy, Narayan, and McKinney 2011, equation (1), pp. L79–L80. These coefficients do not independently specify an arbitrary load: departing from the matched state requires the corresponding load factor and field solution. A different flux or field-unit normalization moves factors of two or 4π4\pi into κ\kappa.

Two independent limits check the scaling. If a∗=0a_*=0, then ΩH=0\Omega_H=0 and rotational extraction vanishes. If ΦBH=0\Phi_{\rm BH}=0, there is no horizon-threading magnetic lever arm and the power also vanishes. At small spin, ΩH≃a∗/(4GM)\Omega_H\simeq a_*/(4GM), so fixed-flux power is quadratic in a∗a_*. Near extremality, higher-order geometry and the spin dependence of the supplied flux make a pure a∗2a_*^2 extrapolation unreliable.

The derivation assumes a stationary, axisymmetric, force-free region and a regular future horizon. Real accretion flows provide the magnetic flux, close the current, load the funnel with pairs or baryons, and can reconnect or interchange on dynamical times. The next page treats that separate magnetic-flux supply problem. A large formal PBZP_{\rm BZ} is not realized unless the accretion flow can retain the required coherent flux.

Nor does horizon extraction identify every observed jet with this mechanism. A measured radio luminosity is separated from horizon power by acceleration, dissipation, radiative efficiency, beaming, duty cycle, and environmental work. Disk winds and disk-powered magnetic torques can coexist with a black-hole-powered spine. Demonstrating Blandford–Znajek operation in a model requires a closed energy and angular-momentum budget with consistent orientations; inferring it in a source requires excluding credible alternatives rather than matching one power scaling.

Load matching and irreversible area growth. Let a flux tube have fixed ψ\psi and power P(x)=4P0x(1−x)P(x)=4P_0x(1-x), where x=ΩF/ΩHx=\Omega_F/\Omega_H and 0<x<10<x<1. Find the maximizing xx. At x=1/2x=1/2, express the horizon-area production term (κH/8πG)A˙H(\kappa_H/8\pi G)\dot A_H in terms of the outward power.

Solution

Differentiating gives dP/dx=4P0(1−2x)dP/dx=4P_0(1-2x), so the stationary point is x=1/2x=1/2. The second derivative is −8P0<0-8P_0<0, hence it is the maximum. From the first-law relation,

κH8πGA˙H=Pout(1x−1).\frac{\kappa_H}{8\pi G}\dot A_H =P_{\rm out}\left(\frac1x-1\right).

At x=1/2x=1/2 the bracket equals one, so the positive area-production term equals PoutP_{\rm out}. This equality is specific to the matched value; it does not mean that the area itself carries an outward energy flux.

  • 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.
  • Tchekhovskoy, Alexander, Ramesh Narayan, and Jonathan C. McKinney. “Efficient Generation of Jets from Magnetically Arrested Accretion on a Rapidly Spinning Black Hole.” Monthly Notices of the Royal Astronomical Society: Letters 418 (2011): L79–L83. DOI.