Relativistic Jets and Launching Mechanisms

Relativistic jets are the most far-reaching structures black holes make. From a launching region only a few gravitational radii across (rg=GM/c2r_g = GM/c^2), a narrow beam of magnetized plasma stays coherent out to megaparsecs in giant radio galaxies — nine orders of magnitude in scale. The same phenomenon appears at every black-hole mass: stellar holes in X-ray binaries flicker jets on and off with accretion state, collapsing stellar cores fire gamma-ray-burst jets, and supermassive holes drive the quasar and blazar zoo. The puzzle that unites them is stark: how does a hot, chaotic, inflowing disk produce an outflow collimated to a few degrees, moving at Lorentz factors of tens, and almost free of baryons? The answer, settled over the past half-century, is that jets are launched electromagnetically: energy leaves the engine not as hot gas but as Poynting flux — stored in a large-scale, ordered magnetic field — and only converts to bulk kinetic energy far downstream. Two mechanisms, both from Roger Blandford's group, define the field. Blandford–Znajek taps the black hole's own rotation: field lines threading the spinning horizon are wound up by frame-dragging in the ergosphere, driving a Poynting outflow like a spinning magnet acting as a dynamo. Blandford–Payne instead flings disk gas outward along inclined field lines like beads on rotating wires. General-relativistic magnetohydrodynamic (GRMHD) simulations added the decisive piece — the Magnetically Arrested Disk, in which accumulated flux chokes the inflow and pushes jet efficiency past 100% of the accreted rest-mass energy, the fingerprint that the hole's spin is being drained. Jet launching thus links black-hole spacetime to galaxy-scale feedback, the M–sigma relation, and the blazar sky.

Key takeaways
  • A relativistic jet is a collimated, magnetically driven outflow launched near an accreting black hole; blazar jets reach Γ∼5\Gamma \sim 5–5050 and stay coherent from rgr_g to megaparsecs.
  • Blandford–Znajek extracts the black hole's spin energy through field lines threading the horizon, powered by ergospheric frame-dragging; the power scales as PBZ∝ΦBH2 ΩH2P_{\rm BZ} \propto \Phi_{\rm BH}^2\,\Omega_H^2.
  • Blandford–Payne extracts the disk's rotational energy, centrifugally slinging gas along field lines inclined more than 30∘30^\circ from the rotation axis.
  • Both deliver energy as Poynting flux first; the jet starts magnetically dominated (σ≫1\sigma \gg 1, slow) and accelerates to Γ∞≈σ0\Gamma_\infty \approx \sigma_0 only as the field energy converts to motion.
  • In the Magnetically Arrested Disk state the horizon flux saturates and jet efficiency can exceed 100%100\% of M˙c2\dot{M}c^2 — direct proof that spin, not accretion, powers the jet; beamed toward us the jet becomes a blazar, its flux boosted by δ3+α\delta^{3+\alpha} with apparent superluminal motion up to tens of cc.
  • Key idea: a jet is a Poynting-flux-dominated outflow — energy stored in an ordered magnetic field threading either the spinning horizon (Blandford–Znajek) or the rotating disk (Blandford–Payne) is extracted electromagnetically; rotation twists the field into a toroidal component that carries energy and angular momentum outward and, through hoop stress, collimates the flow into a beam.
The core insight

A jet is a Poynting-flux-dominated outflow: energy is stored in a large-scale, ordered magnetic field threading either the spinning horizon (Blandford–Znajek) or the rotating disk (Blandford–Payne) and is extracted electromagnetically. Rotation twists the field lines, the resulting toroidal component carries angular momentum and energy outward as an electromagnetic torque, and hoop stress from that toroidal field collimates the flow into a narrow beam.

Where this lives today

Jet launching is now a precision, multi-messenger subject. The Event Horizon Telescope resolved the shadow of M87* in 2019 and, in 2021, mapped an ordered, spiral horizon-scale field — exactly the geometry a magnetically arrested, Blandford–Znajek jet requires — while decades of VLBI track the jet's parabolic-to-conical collimation out to the Bondi radius. Global GRMHD codes (H-AMR, BHAC, KORAL) run three-dimensional, radiative, two-temperature models that reproduce over-unity efficiencies and feed EHT image libraries to constrain spin. The physics scales down: long gamma-ray bursts are modeled as Blandford–Znajek jets from newborn spinning stellar holes; tidal disruption events like Swift J1644+57 launched jets in real time; and X-ray binaries switch jets on and off across the hard-to-soft transition within weeks. Open problems stay sharp — how jets resist kink instabilities over megaparsecs, where the magnetization drops to power the radiation, how much flux builds a MAD, and what loads the baryons that set σ0\sigma_0. Each is a launching question in disguise.

The physical picture

Picture the field lines as stiff wires anchored in the event horizon or the disk. As the anchor rotates, each wire winds up, growing a toroidal (azimuthal) component that does two jobs at once: its tension along the axis pinches the flow inward — hoop-stress collimation — while its outward Poynting flux carries energy and angular momentum away. It is the torque of an electric motor run in reverse, as a generator: rotation of the anchor becomes an outgoing electromagnetic energy flux. The first figure shows the geometry and the governing relation.

Figure
Magnetic launching of a jet from a spinning black hole: poloidal flux Φ_BH threading the horizon is wound up by frame-dragging at rate Ω_H into a toroidal field whose hoop stress collimates the flow, while the Poynting power carried outward scales as P_BZ ∝ Φ_BH² Ω_H²
Electromagnetic launching: spin winds Φ_BH into a toroidal field that collimates and powers the jet spin Ω_H poloidal flux Φ_BH accretion disk ergosphere Poynting flux out toroidal field → hoop stress collimates P_BZ ∝ Φ_BH² Ω_H²

The two mechanisms differ only in where the anchor sits. In Blandford–Payne the wires root in the disk, the rotation is Keplerian, and disk plasma is flung out once a field line leans far enough from vertical. In Blandford–Znajek the wires thread the horizon, the rotation is ΩH\Omega_H of the spacetime, and the process works with no accreting gas at all, draining the hole's rotation through the vacuum field. Which dominates is set by how much ordered flux the flow can pin near the hole; when that flux saturates, the disk becomes magnetically arrested and the jet over-efficient. The second figure is the punchline the simulations deliver.

Figure
Jet efficiency η = P_jet/⟨Ṁ⟩c² versus black-hole spin a in the magnetically arrested state: because η ∝ Ω_H²(a) the curve climbs steeply, crosses 100% near a ≈ 0.93, and reaches η ≈ 1.4 at a ≈ 0.98 — far above the maximum accretion efficiency η_acc = 0.42, so the surplus must come from the hole's spin (schematic, GRMHD-calibrated to Tchekhovskoy, Narayan and McKinney 2011)
The punchline: jet efficiency crosses 100% — the hole's spin is being tapped 0 0.5 1.0 1.5 0 0.25 0.50 0.75 1.0 black-hole spin a jet efficiency η = P_jet / ⟨Ṁ⟩c² η_acc = 0.42 (max thin-disk accretion, a = 1) η = 100% spin-extraction threshold crosses 100% at a ≈ 0.93 MAD, a ≈ 0.98: η ≈ 1.4

The core equations

A small set of relations connects spin, flux, magnetization, and boosting. The central one is the Blandford–Znajek power: a spinning, threaded horizon radiates energy at a rate fixed by the flux it carries and how fast it turns.

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

where: PBZP_{\rm BZ} — power extracted from the hole; ΦBH\Phi_{\rm BH} — poloidal flux through one horizon hemisphere; ΩH\Omega_H — horizon angular velocity; κ≈0.053\kappa \approx 0.053 — geometric constant (split-monopole); f(ΩH)≈1f(\Omega_H) \approx 1 — higher-spin correction. Read it as: power grows with the square of both flux and spin, so a jet needs a rapidly rotating hole and a strong ordered field. The horizon's angular velocity is a pure spacetime quantity, set by the mass and spin of the Kerr solution; it vanishes for a non-rotating hole and saturates near maximal spin.

ΩH=a c32 GM 11+1−a2(2)\Omega_H = \frac{a\,c^3}{2\,G M}\,\frac{1}{1+\sqrt{1-a^2}} \qquad (2)

where: a=Jc/(GM2)a = Jc/(GM^2) — dimensionless spin, 0≤a≤10 \le a \le 1; MM — mass; JJ — angular momentum. For a≪1a \ll 1, ΩH≈ac/(2rg)\Omega_H \approx a c/(2 r_g) so PBZ∝a2P_{\rm BZ} \propto a^2; as a→1a \to 1, ΩH→c/(2rg)\Omega_H \to c/(2 r_g) — a Schwarzschild hole drives no BZ jet. To compare jets across many decades of mass, power is written as an efficiency — jet power per unit rest-mass energy supplied by accretion.

ηjet≡Pjet⟨M˙⟩ c2(3)\eta_{\rm jet} \equiv \frac{P_{\rm jet}}{\langle\dot{M}\rangle\,c^2} \qquad (3)

where: ηjet\eta_{\rm jet} — jet efficiency; ⟨M˙⟩\langle\dot{M}\rangle — time-averaged accretion rate. Read it as a budget: accretion caps radiative efficiency at ηacc≤0.42\eta_{\rm acc} \le 0.42 (extreme Kerr), so any ηjet>1\eta_{\rm jet} > 1 is energy the disk cannot supply — it is spin. The efficiency saturates because the flux does. Once enough poloidal flux piles up, its pressure pushes back on the inflow and the disk becomes magnetically arrested at a near-universal flux.

ΦMAD=ϕMAD M˙ rg2 c ,ϕMAD≈50(4)\Phi_{\rm MAD} = \phi_{\rm MAD}\,\sqrt{\dot{M}\,r_g^2\,c}\,,\qquad \phi_{\rm MAD}\approx 50 \qquad (4)

where: ϕMAD≈50\phi_{\rm MAD} \approx 50 — dimensionless saturation flux (Gaussian-cgs GRMHD units); rg=GM/c2r_g = GM/c^2. The flux the hole can hold is set by M˙\dot{M}, which is why PBZ∝Φ2∝M˙P_{\rm BZ} \propto \Phi^2 \propto \dot{M} in the arrested state and efficiency becomes a function of spin alone. Whether the jet is a battery or a bullet is set by the magnetization — field energy over rest-mass energy.

σ≡B24π ρ c2=Poynting fluxrest-mass flux(5)\sigma \equiv \frac{B^2}{4\pi\,\rho\,c^2} = \frac{\text{Poynting flux}}{\text{rest-mass flux}} \qquad (5)

where: σ\sigma — magnetization; BB — comoving field; ρ\rho — comoving rest-mass density. At the base σ≫1\sigma \gg 1 (magnetically dominated, slow); acceleration is the trade of σ\sigma for bulk motion. That trade is governed by a quantity conserved along each field line — total energy per rest-mass flux — which fixes the terminal speed.

μ=Γ (1+σ)=const    ⟹    Γ∞≈μ≈σ0(6)\mu = \Gamma\,(1+\sigma) = \text{const} \;\;\Longrightarrow\;\; \Gamma_\infty \approx \mu \approx \sigma_0 \qquad (6)

where: μ\mu — conserved energy-to-mass flux; Γ\Gamma — local Lorentz factor; Γ∞\Gamma_\infty — terminal Lorentz factor; σ0\sigma_0 — base magnetization. Read it as a ceiling: the jet reaches at most Γ∞=μ\Gamma_\infty = \mu, and only if the field energy fully converts (σ→0\sigma \to 0). Finally, what we observe is not Γ\Gamma but the Doppler-boosted appearance, set by the beaming factor for a flow of speed β\beta seen at θ\theta.

δ=1Γ (1−βcos⁡θ)(7)\delta = \frac{1}{\Gamma\,(1-\beta\cos\theta)} \qquad (7)

where: δ\delta — Doppler factor; β=v/c\beta = v/c; θ\theta — angle to the line of sight. Head-on, δ≈2Γ\delta \approx 2\Gamma; at cos⁡θ=β\cos\theta = \beta (θ≈1/Γ\theta \approx 1/\Gamma), δ=Γ\delta = \Gamma. Observed flux scales as δ3+α\delta^{3+\alpha}, so the appearance is dominated by δ\delta, not Γ\Gamma.

Derivation AMagnetocentrifugal launching — the bead-on-a-wire criterion

The goal is to find the geometric condition under which a cold parcel of disk gas, frozen to a rigidly rotating field line, is flung outward with no energy barrier — the Blandford–Payne criterion. We treat the field line as a stiff wire corotating at the Keplerian rate, write the effective potential a bead feels along it, and ask when the footpoint is a potential maximum rather than a minimum, so the slightest outward nudge slides the bead to infinity. The answer is a single critical inclination angle.

Step 1. In the frame corotating at the footpoint rate Ω0\Omega_0, a unit-mass parcel feels gravity plus a centrifugal potential, Φeff(r,z)=−GMr2+z2−12Ω02r2\Phi_{\rm eff}(r,z) = -\dfrac{GM}{\sqrt{r^2+z^2}} - \tfrac{1}{2}\Omega_0^2 r^2, with Ω02=GM/r03\Omega_0^2 = GM/r_0^3 (Keplerian footpoint at r0r_0). Why: corotation makes the field "rigid," and centrifugal support is the whole point, so it must live in the potential.

Step 2. Parametrize a straight field line by arc length ss from the footpoint, inclined at angle θ\theta above the disk plane: r=r0+scos⁡θr = r_0 + s\cos\theta, z=ssin⁡θz = s\sin\theta. At the footpoint the net force along the wire vanishes, dΦeff/ds∣0=0\mathrm{d}\Phi_{\rm eff}/\mathrm{d}s\big|_0 = 0. Why: the launch is decided not by the slope (zero) but by the curvature — whether the footpoint is a hilltop or a valley.

Step 3. Expand to second order, Φeff(s)≈Φeff(0)+12Φ′′(0) s2\Phi_{\rm eff}(s) \approx \Phi_{\rm eff}(0) + \tfrac{1}{2}\Phi''(0)\,s^2. Setting GM=r0=1GM=r_0=1, differentiation gives Φ′′(0)=−3cos⁡2θ+sin⁡2θ=1−4cos⁡2θ\Phi''(0) = -3\cos^2\theta + \sin^2\theta = 1 - 4\cos^2\theta. Why: steep gravitational curvature (−3-3) along the radial direction and gentle centrifugal curvature combine into one expression whose sign is all that matters.

Step 4. Launching requires Φ′′(0)<0\Phi''(0) < 0, i.e. 1−4cos⁡2θ<0⇒cos⁡2θ>1/4⇒θ<60∘1 - 4\cos^2\theta < 0 \Rightarrow \cos^2\theta > 1/4 \Rightarrow \theta < 60^\circ from the disk plane — equivalently, more than 30∘30^\circ from the rotation axis. Why: below 60∘60^\circ the centrifugal barrier vanishes, so cold gas accelerates from rest with no thermal push, which is why jets are cold and fast rather than hot and slow. *The marginal inclination is arccos⁡(1/2)=60∘\arccos(1/2)=60^\circ from the disk plane (30° from the axis), confirmed by the sign change in Φ′′(0)=1−4cos⁡2θ\Phi''(0)=1-4\cos^2\theta. Beyond the Alfvén point the flow decouples from rigid corotation and is collimated by toroidal hoop stress.*

Go deeper · Derivation B · 1 minBlandford–Znajek power from horizon impedance matching
Derivation BBlandford–Znajek power from horizon impedance matching

The goal is to recover the PBZ∝Φ2ΩH2P_{\rm BZ} \propto \Phi^2\Omega_H^2 scaling from a circuit picture, treating the rotating horizon as a battery of finite internal resistance driving a load. The membrane paradigm gives the horizon a real surface resistivity; frame-dragging sets up an EMF, and the power delivered to the load (the jet) is maximized by impedance matching, as for any generator. This turns a general-relativistic result into Ohm's law and explains why power is quadratic in flux and spin.

Step 1. Frame-dragging rotates the field lines with the horizon, generating a unipolar-inductor EMF, E∼ΩHΦBH/c\mathcal{E} \sim \Omega_H \Phi_{\rm BH}/c — the physics of a Faraday disc. Why: relative rotation between the field's anchor and infinity does the work; ΩH\Omega_H and Φ\Phi enter linearly, so their product sets the driving voltage.

Step 2. Assign the horizon its membrane surface impedance ZH=4π/c≈377 ΩZ_H = 4\pi/c \approx 377\,\Omega, the impedance of free space — the "internal resistance" of the battery. Why: the horizon behaves as a resistive membrane in the 3+1 split; without finite impedance there is no dissipation and no steady power transfer.

Step 3. Maximum power is delivered when the external field geometry matches this impedance: P∼E2ZH∼c4π(ΩHΦBHc)2P \sim \dfrac{\mathcal{E}^2}{Z_H} \sim \dfrac{c}{4\pi}\left(\dfrac{\Omega_H \Phi_{\rm BH}}{c}\right)^2. Why: this is the maximum-power-transfer theorem; the factor c/4πc/4\pi becomes κ≈0.053\kappa \approx 0.053 once the split-monopole geometry is worked out.

Step 4. Collecting factors gives PBZ=κ4πc ΦBH2 ΩH2P_{\rm BZ} = \dfrac{\kappa}{4\pi c}\,\Phi_{\rm BH}^2\,\Omega_H^2. Why: the two squares are distinct — Φ2\Phi^2 because power is voltage-squared over a fixed impedance, ΩH2\Omega_H^2 because the EMF is itself proportional to the spin rate. *The impedance-matching estimate reproduces the Φ2ΩH2\Phi^2\Omega_H^2 scaling and the κ≈0.053\kappa\approx0.053 geometric prefactor of the split-monopole result.*

Go deeper · Derivation C · 1 minPoynting-to-kinetic conversion and the terminal Lorentz factor
Derivation CPoynting-to-kinetic conversion and the terminal Lorentz factor

The goal is to show why a jet carrying enormous energy starts out moving slowly, and what ultimate speed it reaches. For a steady, ideal, axisymmetric magnetized outflow the total energy flux per rest-mass flux is conserved along each field line; split between kinetic energy and Poynting flux, this invariant forces the terminal Lorentz factor to equal the launch-point magnetization once the field energy is fully mined. We track it from base to infinity.

Step 1. The relativistic Bernoulli invariant for a cold magnetized wind is μ=Γ(1+σ)=const\mu = \Gamma(1+\sigma) = \text{const} along a field line. Why: energy is conserved in ideal MHD, so what the jet gains in Γ\Gamma it loses in σ\sigma — acceleration is literally magnetic energy being spent.

Step 2. At the base the flow is nearly static and field-dominated, Γ0≈1\Gamma_0 \approx 1, σ0≫1\sigma_0 \gg 1, so μ≈σ0\mu \approx \sigma_0. Why: this resolves the apparent paradox — the jet carries a huge flux μ\mu, all stored magnetically, so the plasma barely moves.

Step 3. As magnetic energy converts, σ\sigma falls and Γ\Gamma rises; complete conversion (σ→0\sigma \to 0) gives Γ∞→μ≈σ0\Gamma_\infty \to \mu \approx \sigma_0. Why: the field sets a hard speed ceiling — you cannot accelerate past the energy you started with — and real jets convert incompletely, retaining σ∼0.1\sigma \sim 0.1–11.

Step 4. The equipartition milestone σ=1\sigma = 1 gives Γ≈μ/2\Gamma \approx \mu/2: half the energy is now kinetic, and this fast-magnetosonic point is where acceleration slows and dissipation (shocks, reconnection) lights up the blazar emission. Why: it marks the handoff from a launching problem to a radiation problem — what telescopes actually detect. *Here Γ∞=μ=Γ(1+σ)\Gamma_\infty = \mu = \Gamma(1+\sigma) is the conserved Bernoulli flux of ideal MHD for relativistic magnetized winds.*

Limiting cases

The mechanisms and regimes, read off the equations above:

CaseResultMeaning
No spin (BZ)PBZ∝ΩH2∝a2→0P_{\rm BZ} \propto \Omega_H^2 \propto a^2 \to 0No ergosphere, no frame-dragging EMF, no Blandford–Znajek jet — rotation is essential.
Near-maximal spin (BZ)ΩH→c/(2rg)\Omega_H \to c/(2 r_g); MAD ηjet≳1\eta_{\rm jet} \gtrsim 1Horizon spin saturates; with arrested flux the jet outpowers accretion — spin dominates.
Launching angle (BP)θ<60∘\theta < 60^\circ from disk plane (>30∘>30^\circ from axis)Field lines must lean far enough from vertical; steeper lines keep a centrifugal barrier.
Magnetically dominated baseσ≫1, Γ≈1\sigma \gg 1,\ \Gamma \approx 1Energy stored in the field; powerful but slow, and μ≈σ0\mu \approx \sigma_0 sets the ultimate speed.
Fully converted jetσ→0, Γ∞→μ\sigma \to 0,\ \Gamma_\infty \to \muAll magnetic energy became bulk motion; real jets stall near σ∼0.1\sigma \sim 0.1–11.
Beaming coneθ≲1/Γ, δ∼Γ\theta \lesssim 1/\Gamma,\ \delta \sim \GammaInside the cone the jet is a blazar; outside, the same source looks like a faint radio galaxy.

Worked examples

1 · The spin reservoir and spin-down time of M87*

M87 hosts a black hole of M=6.5×109 M⊙M = 6.5\times10^9\,M_\odot with spin a≈0.9a \approx 0.9. How much rotational energy is available to a Blandford–Znajek jet, and how long could a 1044 erg s−110^{44}\,\mathrm{erg\,s^{-1}} jet run on it?

Step 1. The rest energy is Mc2=(6.5×109)(1.989×1033 g)(2.998×1010)2=1.16×1064 ergMc^2 = (6.5\times10^9)(1.989\times10^{33}\,\mathrm{g})(2.998\times10^{10})^2 = \mathbf{1.16\times10^{64}\ erg}. Only the reducible part above the irreducible mass can be extracted.

Step 2. The extractable fraction for spin aa is 1−12(1+1−a2)1 - \sqrt{\tfrac{1}{2}(1+\sqrt{1-a^2})}. For a=0.9a=0.9 this is 1−12(1.436)=0.1531 - \sqrt{\tfrac12(1.436)} = \mathbf{0.153}, so Erot=0.153 Mc2=1.8×1063 ergE_{\rm rot} = 0.153\,Mc^2 = \mathbf{1.8\times10^{63}\ erg}. (At extremal a=1a=1 it rises to 0.29 Mc2≈3.4×10630.29\,Mc^2 \approx 3.4\times10^{63} erg.)

Step 3. The spin-down time is tspin∼Erot/Pjet=(1.8×1063)/(1044)=1.8×1019 s≈6×1011 yrt_{\rm spin} \sim E_{\rm rot}/P_{\rm jet} = (1.8\times10^{63})/(10^{44}) = \mathbf{1.8\times10^{19}\ s \approx 6\times10^{11}\ yr} — roughly 40 times the age of the universe. The spin is effectively a permanent battery; jet activity is limited by fuel supply and duty cycle, not by draining the reservoir. (Nuclear fusion liberates only 0.7%0.7\% of rest mass, so spin at ≈15%\approx 15\% here, up to 29%29\%, is a vastly larger tank.)

2 · Doppler beaming and superluminal motion in a blazar

A jet with Γ=10\Gamma = 10 points θ=5∘\theta = 5^\circ from our line of sight. Find its Doppler factor, flux boost, and apparent transverse speed.

Step 1. The bulk speed is β=1−1/Γ2=0.99499\beta = \sqrt{1-1/\Gamma^2} = \mathbf{0.99499}. With cos⁡5∘=0.99619\cos 5^\circ = 0.99619, the beaming denominator is 1−βcos⁡θ=1−(0.99499)(0.99619)=0.008801 - \beta\cos\theta = 1 - (0.99499)(0.99619) = \mathbf{0.00880}.

Step 2. The Doppler factor is δ=1/[Γ(1−βcos⁡θ)]=1/[10×0.00880]=11.4\delta = 1/[\Gamma(1-\beta\cos\theta)] = 1/[10\times0.00880] = \mathbf{11.4} — larger than Γ\Gamma itself because we sit inside the beaming cone (θ<1/Γ=5.7∘\theta < 1/\Gamma = 5.7^\circ); head-on the maximum would be δmax⁡≈2Γ=20\delta_{\max}\approx2\Gamma=20.

Step 3. The flux is boosted by δ3+α\delta^{3+\alpha}; a moving-blob δ4\delta^4 gives 11.44≈1.7×10411.4^4 \approx \mathbf{1.7\times10^{4}}. An intrinsically ordinary jet is amplified ten-thousand-fold — which is why blazars dominate the extragalactic gamma-ray sky.

Step 4. The apparent transverse speed is βapp=βsin⁡θ1−βcos⁡θ=(0.99499)(0.0872)0.00880=9.9\beta_{\rm app} = \dfrac{\beta\sin\theta}{1-\beta\cos\theta} = \dfrac{(0.99499)(0.0872)}{0.00880} = \mathbf{9.9}. The blob appears to move at 9.9c9.9c — superluminal motion, and no violation of relativity: a light-travel-time illusion of a sub-luminal flow aimed nearly at us.

3 · Blandford–Znajek power of M87* in the arrested state

With accretion rate M˙≈10−3 M⊙ yr−1\dot{M} \approx 10^{-3}\,M_\odot\,\mathrm{yr^{-1}} and MAD efficiency η≈0.8\eta \approx 0.8 at a=0.9a=0.9, estimate the jet power and compare with observations.

Step 1. Convert the rate: M˙=10−3×(1.989×1033)/(3.156×107)=6.3×1022 g s−1\dot{M} = 10^{-3}\times(1.989\times10^{33})/(3.156\times10^7) = \mathbf{6.3\times10^{22}\ g\,s^{-1}}, so the supply is M˙c2=(6.3×1022)(8.99×1020)=5.7×1043 erg s−1\dot{M}c^2 = (6.3\times10^{22})(8.99\times10^{20}) = \mathbf{5.7\times10^{43}\ erg\,s^{-1}}.

Step 2. The jet power is Pjet=η M˙c2=0.8×5.7×1043=4.5×1043 erg s−1P_{\rm jet} = \eta\,\dot{M}c^2 = 0.8\times5.7\times10^{43} = \mathbf{4.5\times10^{43}\ erg\,s^{-1}}, squarely in the observed range 104310^{43}–1044 erg s−110^{44}\,\mathrm{erg\,s^{-1}} from M87's radio lobes and X-ray cavities.

Step 3. Push the spin to a→1a \to 1, where the same MAD flux gives η≈1.4\eta \approx 1.4: then Pjet=1.4×5.7×1043=8×1043 erg s−1>M˙c2P_{\rm jet} = 1.4\times5.7\times10^{43} = \mathbf{8\times10^{43}\ erg\,s^{-1}} > \dot{M}c^2. This exceeds everything accretion can supply — the surplus is spin energy, the smoking gun of Blandford–Znajek.

Physics of the portrait

Relativistic jets have a single dominant author. Both launching mechanisms carry Roger Blandford's name, and the flat-spectrum theory of their emission is his too — the man who turned a "curious straight ray" in a nearby galaxy into black-hole electrodynamics.

Origins · 1 minRoger Blandford and how this was discovered
RB
Roger Blandford
b. 1949 · relativistic astrophysics and black-hole electrodynamics

Roger Blandford, born in 1949 and trained at Cambridge, built the theoretical scaffolding of nearly everything in this article. In 1977, with his student Roman Znajek, he showed that field lines threading a rotating Kerr horizon set up an electromagnetic circuit that extracts the hole's spin as a Poynting jet — the Blandford–Znajek process, still the leading explanation for the most powerful jets in the universe. In 1979, with Ari Königl, he explained the flat spectra of compact radio cores as a self-absorbed conical relativistic jet, the model still used to read VLBI images. In 1982, with David Payne, he derived the magnetocentrifugal disk-wind mechanism and its bead-on-a-wire criterion, the standard picture for disk-driven outflows from protostars to AGN. Blandford spent decades at Caltech before founding the Kavli Institute for Particle Astrophysics and Cosmology at Stanford, and chaired the 2010 US decadal survey. His hallmark is physical clarity — recasting curved-spacetime electrodynamics as circuits, impedances, and torques — which is why the field's foundations are so durable.

History
1918
— Curtis records a "curious straight ray" — the optical jet of M87
1969
— Penrose shows rotational energy can be extracted from a spinning black hole
1977
— Blandford and Znajek: electromagnetic energy extraction from Kerr holes
1982
— Blandford and Payne: the magnetocentrifugal disk-wind launching criterion
2011
— Tchekhovskoy, Narayan and McKinney: magnetically arrested jets exceed 100% efficiency

Going further

Texts. Carroll & Ostlie, An Introduction to Modern Astrophysics, 2nd ed., Ch. 17 (Kerr) and Ch. 28 (active galaxies) · Binney & Tremaine, Galactic Dynamics, 2nd ed. · Mo, van den Bosch & White, Galaxy Formation and Evolution · Blandford & Znajek (1977), MNRAS 179, 433 · Blandford & Payne (1982), MNRAS 199, 883 · Tchekhovskoy, Narayan & McKinney (2011), MNRAS 418, L79 · Event Horizon Telescope Collaboration (2019, 2021). Key numerical anchors: the PBZ∝Φ2ΩH2P_{\rm BZ}\propto\Phi^2\Omega_H^2 scaling carries κ≈0.053\kappa\approx0.053. Horizon angular velocity and irreducible mass are standard Kerr results: a=1⇒Mirr=M/2a=1\Rightarrow M_{\rm irr}=M/\sqrt{2}, Erot≈0.29 Mc2E_{\rm rot}\approx0.29\,Mc^2; a=0.9⇒Erot=0.153 Mc2a=0.9\Rightarrow E_{\rm rot}=0.153\,Mc^2. The Blandford–Payne launching angle follows from Φ′′(0)=1−4cos⁡2θ\Phi''(0)=1-4\cos^2\theta, giving arccos⁡(1/2)=60∘\arccos(1/2)=60^\circ from the disk plane (30° from the axis) — note this corrects a common misstatement of 54.7∘54.7^\circ. The terminal Lorentz factor Γ∞=μ=Γ(1+σ)\Gamma_\infty=\mu=\Gamma(1+\sigma) is the ideal-MHD Bernoulli invariant; ϕMAD≈50\phi_{\rm MAD}\approx50 and η∼1.4\eta\sim1.4 at high spin. Worked-example numbers: Mc2=1.16×1064Mc^2=1.16\times10^{64} erg, δ=11.4\delta=11.4, βapp=9.9\beta_{\rm app}=9.9, M˙c2=5.7×1043 erg s−1\dot{M}c^2=5.7\times10^{43}\,\mathrm{erg\,s^{-1}}.

Builds on: Black Hole Spin and the Ergosphere · Thin Accretion Disks and the Shakura–Sunyaev Model · Magnetorotational Instability and Accretion Turbulence · Advection-Dominated Accretion Flows. Leads to: Blazars and Relativistic Beaming · AGN Feedback · The M–sigma Relation and Black Hole–Galaxy Coevolution · Active Galactic Nuclei: the Unified Model.

Ask the tutor: Want to see how far a magnetically arrested disk can push jet efficiency at your chosen spin, or how the beaming cone reshapes a jet into a blazar — pick a black hole and let's launch its jet together.

How can a jet be more powerful than everything the accretion disk supplies?

It cannot — if accretion were the only source. Converting infalling rest mass to radiation caps out near 42%42\% even for a maximally spinning hole, so an efficiency above 100%100\% of M˙c2\dot{M}c^2 is energy the disk does not have. The resolution is a second reservoir: the rotational energy of the black hole itself, extracted electromagnetically by the Blandford–Znajek mechanism. GRMHD simulations of magnetically arrested disks routinely reach ηjet∼1.4\eta_{\rm jet}\sim1.4 at high spin, and this over-unity result is the cleanest dynamical proof that a jet drains a black hole's spin, spinning it down over cosmological times.

What is the difference between the Blandford–Znajek and Blandford–Payne mechanisms?

Both convert rotation into a Poynting-flux jet through twisted field lines, but they tap different flywheels. Blandford–Znajek threads field lines through the event horizon and extracts the black hole's spin; it needs an ergosphere and frame-dragging, and works even for a magnetically threaded hole with no gas accreting. Blandford–Payne anchors field lines in the accretion disk and extracts the disk's orbital energy, flinging gas centrifugally along inclined lines. In practice the two coexist: BZ powers a fast, low-density, highly magnetized spine along the axis, while BP drives a slower, denser, mass-loaded sheath around it — a structure now seen in both simulations and observations.

Why does the black hole's spin, not just its mass, matter for jet power?

Because the Blandford–Znajek power depends on the horizon's angular velocity, PBZ∝ΩH2P_{\rm BZ}\propto\Omega_H^2, and ΩH\Omega_H is set by spin: zero for a non-rotating (Schwarzschild) hole, rising toward c/2rgc/2r_g near maximal spin. A non-spinning hole has no ergosphere, no frame-dragging, no induced EMF, and no BZ jet — however massive it is or however fast it accretes. Mass sets the overall scale (through rgr_g and the sustainable accretion rate), but spin is the on/off switch and the throttle. That is why jet power, with the accretion rate, is one of the few observational handles on black-hole spin.

What is a magnetically arrested disk, and why does it matter?

A magnetically arrested disk (MAD) is a state in which so much ordered poloidal flux has piled up near the horizon that its magnetic pressure pushes back on the inflow, halting it — gas then leaks inward in discrete, bubble-like intervals rather than smoothly. The flux saturates at a near-universal value (ϕMAD≈50\phi_{\rm MAD}\approx50), so MAD jets reach a maximum efficiency that depends almost entirely on spin. It matters because it is the regime where jets are most powerful and where Blandford–Znajek efficiency exceeds unity. The 2021 EHT polarization images of M87* show exactly the ordered, spiral horizon-scale field a MAD state predicts, making it the leading model for the most luminous AGN jets.

Why do jets start out barely moving if they are called "relativistic"?

At the launch point almost all of a jet's energy is stored in its magnetic field, not in bulk motion — the magnetization σ=B2/(4πρc2)\sigma = B^2/(4\pi\rho c^2) is much greater than one. The conserved flux μ=Γ(1+σ)\mu = \Gamma(1+\sigma) can be enormous while Γ≈1\Gamma\approx1, because the large factor lives in σ\sigma. As the flow expands, magnetic energy converts to kinetic: σ\sigma falls, Γ\Gamma climbs, and the jet reaches its terminal Γ∞≈σ0\Gamma_\infty\approx\sigma_0 only far downstream. So "relativistic" describes the terminal state; the base is a slow, magnetically dominated launcher. Incomplete conversion leaves residual magnetization that later powers the radiation through shocks and reconnection.

What actually collimates a jet into such a narrow beam?

The same toroidal field that carries the energy also collimates it. As the anchor rotates, it winds the field into tight azimuthal loops, and the tension of those loops — the magnetic "hoop stress" — squeezes inward toward the axis, like stretched rubber bands pinching a column of fluid. This self-collimation focuses the flow into an opening angle of a few degrees, far narrower than any thermal outflow, and works over enormous distances because the ordered field is preserved along the jet. External pressure from the disk wind or interstellar medium tunes the profile, but the core mechanism is the jet's own toroidal field.

Before this

Where this leads