Relativistic Jets from Compact Objects

Point a radio telescope at almost any accreting black hole and you find the same thing: two needle-thin beams of plasma shooting out along the rotation axis, oppositely directed, at very nearly light speed. In the nuclei of distant galaxies these jets span millions of light-years; in our own Galaxy the same structures, scaled down by 10610^{6}–10910^{9} in central mass, emerge from stellar-mass black holes and neutron stars in X-ray binaries — the microquasars. The scaling is not a metaphor: the physics that fans a quasar's jet across a galaxy cluster is the physics that flings a blob out of GRS 1915+105 fast enough to appear to cross the sky faster than light. That last fact is the entry point. A relativistic jet turns the free energy of infalling matter — or the rotational energy of the compact object itself — into a magnetically collimated, Poynting-flux-dominated outflow that radiates non-thermally as it accelerates. Three ideas organize everything. The launching is magnetic: Blandford–Znajek taps the spin of a Kerr black hole through horizon-threading field lines, Blandford–Payne slings disc plasma along inclined ones. The motion is relativistic, so geometry distorts what we see: beaming brightens the approaching jet and hides its twin, and small viewing angles manufacture superluminal motion. And the radiation is synchrotron, shock-accelerated electrons painting the spectra that map the flow. Get these three straight and the whole zoo — blazars, radio galaxies, microquasars, one-sided images, apparent 5c5c knots — collapses into one mechanism seen from different angles.

Key takeaways
  • A relativistic jet is a collimated, bipolar outflow from the inner accretion flow of a compact object — stellar-mass holes and neutron stars in microquasars, supermassive holes in AGN — over ∼9\sim 9 decades in central mass.
  • The launch mechanisms are magnetic: Blandford–Znajek extracts black-hole spin energy through horizon-threading field, Blandford–Payne drives a magnetocentrifugal disc wind; both start Poynting-dominated and convert to bulk kinetic energy.
  • Bulk Lorentz factors reach Γ∼2\Gamma\sim 2–1010 in microquasars and Γ≳10\Gamma\gtrsim 10 in blazars; relativistic beaming and superluminal motion are geometric consequences of that speed.
  • The maximum apparent transverse speed is βapp,max=Γβ\beta_{\mathrm{app,max}}=\Gamma\beta, achieved at cos⁡θ=β\cos\theta=\beta — so a measured superluminal speed sets a firm lower bound on Γ\Gamma.
  • Synchrotron from a power-law electron distribution N(γ)∝γ−sN(\gamma)\propto\gamma^{-s} gives an optically thin spectrum Sν∝ν−(s−1)/2S_\nu\propto\nu^{-(s-1)/2}, with s≈2.5s\approx 2.5 from diffusive shock acceleration.
  • The jet–disc connection (radio flares tied to X-ray state transitions) and the fundamental plane LR∝M0.6–0.8LX0.6L_R\propto M^{0.6\text{–}0.8}L_X^{0.6} tie jet production to magnetized accretion across the black-hole population.
  • Key idea: A relativistic jet converts the free energy of an accretion flow — or the spin energy of the central compact object — into a Poynting-flux-dominated, magnetically collimated outflow that accelerates to bulk Lorentz factors of a few to tens and radiates non-thermally; relativistic beaming then governs everything we observe.
The core insight

A relativistic jet converts the free energy of an accretion flow — or the spin energy of the central compact object — into a Poynting-flux-dominated, magnetically collimated outflow, which accelerates to bulk Lorentz factors of a few to tens and radiates non-thermally; relativistic beaming then governs everything we observe.

Where this lives today

Jet physics is a precision, multi-messenger enterprise. The Event Horizon Telescope has resolved the base of M87's jet and imaged the black-hole shadow, letting theorists confront Blandford–Znajek launching against the measured field geometry and polarization at the horizon scale. General-relativistic magnetohydrodynamic simulations now produce magnetically arrested discs (MADs) whose spin-powered jets exceed 100%100\% efficiency in the accreted rest-mass energy — direct evidence that the hole, not just the disc, powers the outflow. In the time domain, rapid radio–X-ray monitoring of binaries such as GX 339-4 traces the jet–disc coupling in real time: the compact jet quenches at the hard-to-soft transition, where discrete superluminal ejecta launch. Blazar campaigns tie Fermi-LAT γ\gamma-ray flares to VLBI knots, and the IceCube neutrino from blazar TXS 0506+056 opened jets as hadronic accelerators. Open problems stay sharp: what sets the terminal Lorentz factor, how efficiently Poynting flux converts to kinetic energy (the "sigma problem"), and whether accreting neutron stars — lacking a horizon — can launch jets by the disc mechanism alone.

The physical picture

Picture the magnetic field lines threading a spinning black hole or its inner disc as stiff wires bolted to a rotating hub. Plasma beads slide freely along each wire but are forced to corotate around the axis. If a disc-anchored line leans more than 30∘30^\circ from vertical — below 60∘60^\circ to the disc plane — the centrifugal potential tips downhill and a bead is flung out like a stone on a sling: the Blandford–Payne launch, needing only a rotating disc and an ordered poloidal field. Where the field instead threads the horizon of a Kerr hole, frame dragging winds it up and the rotational energy of spacetime itself is pumped out electromagnetically — the Blandford–Znajek process, with power PBZ∝Φ2ΩH2P_{\mathrm{BZ}}\propto\Phi^2\Omega_H^2. Near the base the outflow is Poynting-dominated: its magnetization σ=B2/4πρc2≫1\sigma=B^2/4\pi\rho c^2\gg 1 means electromagnetic energy per particle dwarfs the rest mass — a vast acceleration reservoir. As rotation winds the field into a tight toroidal spiral, hoop stress squeezes the flow toward the axis and collimates it, while the stored magnetic energy converts to bulk kinetic energy and the flow coasts out at Γ∼2\Gamma\sim 2–1010.

Figure
A magnetized accretion flow onto a spinning black hole launches twin jets: Blandford–Znajek taps the horizon rotation with power P_BZ ∝ Φ² Ω_H², Blandford–Payne flings disc plasma along field lines inclined under 60°, and Poynting flux σ = B²/4πρc² ≫ 1 converts to bulk kinetic energy reaching Γ ~ 2–10 as hoop stress collimates the flow
One mechanism: magnetized accretion + spin → twin relativistic jets Γ ~ 2–10 bulk kinetic hoop stress collimates field under 60°: Blandford–Payne accretion disc P_BZ ∝ Φ² Ω_H² Blandford–Znajek: taps horizon spin σ = B²/4πρc² ≫ 1 (Poynting-dominated base) counter-jet: intrinsically identical, Doppler-hidden

Two geometric effects dominate what a telescope records. Superluminal motion is pure projection: a blob approaching at 0.98c0.98c nearly keeps pace with its own light, so photons from two events arrive bunched and the blob appears to sweep a large angle in a tiny apparent time — an illusion of transverse speed several times cc. Relativistic beaming is why jets look one-sided: boosting funnels the approaching jet's radiation into a cone of half-angle ∼1/Γ\sim 1/\Gamma and brightens it by δ3–4\delta^{3\text{–}4}, while the counter-jet is de-boosted by the same power — a contrast easily exceeding 10410^4. The second figure makes the illusion quantitative.

Figure
Apparent transverse speed β_app = β sinθ / (1 − β cosθ) versus viewing angle for bulk Lorentz factors Γ = 2, 5, 10: each curve peaks where cosθ = β at the value β_app,max = Γβ (1.73, 4.90, 9.95 — dashed markers), and the Galactic microquasar GRS 1915+105 sits at β_app ≈ 1.27, θ ≈ 69° on the Γ ≈ 2.5 track
The punchline: β_app can exceed 1 — peak Γβ at cosθ = β β_app = 2 (apparent 2c) 0 2 4 6 8 10 0 30 60 90 viewing angle θ (degrees) apparent speed β_app (units of c) Γ = 2 (max 1.73) Γ = 5 (max 4.90) Γ = 10 (max 9.95) GRS 1915+105 β_app ≈ 1.27, θ ≈ 69°

The core equations

The single most consequential relation is the apparent-velocity formula, which follows from comparing two emission events with the light-travel-time compression as the source chases its own photons — the reason a subluminal flow can look superluminal.

βapp=βsin⁡θ1−βcos⁡θ(1)\beta_{\mathrm{app}}=\frac{\beta\sin\theta}{1-\beta\cos\theta} \qquad (1)

where: β\beta — true bulk speed in units of cc; θ\theta — angle between jet velocity and line of sight; βapp\beta_{\mathrm{app}} — apparent transverse speed on the sky, in units of cc. Reading it: the denominator 1−βcos⁡θ1-\beta\cos\theta shrinks toward zero for fast flows aimed near the line of sight, inflating βapp\beta_{\mathrm{app}} far above unity — no signal outruns light, only its image. The same denominator governs all boosting through the Doppler factor, which bundles frequency shift, time compression, and solid-angle focusing into one number.

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

where: δ\delta — Doppler factor; Γ=(1−β2)−1/2\Gamma=(1-\beta^2)^{-1/2} — bulk Lorentz factor; θ\theta — viewing angle. A head-on jet (θ→0\theta\to 0) saturates at δ→2Γ\delta\to 2\Gamma; the counter-jet at δ→1/2Γ\delta\to 1/2\Gamma. Beaming raises that factor to a power set by whether one watches a discrete moving blob or a continuous steady jet.

Sνobs=δ p Sνem,p=3−α (discrete),p=2−α (continuous)(3)S_\nu^{\mathrm{obs}}=\delta^{\,p}\,S_\nu^{\mathrm{em}},\qquad p=3-\alpha\ (\text{discrete}),\quad p=2-\alpha\ (\text{continuous}) \qquad (3)

where: SνS_\nu — flux density at frequency ν\nu; α\alpha — spectral index, Sν∝ναS_\nu\propto\nu^{\alpha}; pp — beaming exponent. The jet/counter-jet ratio is [(1+βcos⁡θ)/(1−βcos⁡θ)]p\left[(1+\beta\cos\theta)/(1-\beta\cos\theta)\right]^{p} — why counter-jets are rarely seen. For the launching itself, the Blandford–Znajek power quantifies the electromagnetic luminosity a spinning hole pumps into the jet from its rotation alone.

PBZ  ∼  κ4πc Φ2 ΩH2(4)P_{\mathrm{BZ}}\;\sim\;\frac{\kappa}{4\pi c}\,\Phi^2\,\Omega_H^2 \qquad (4)

where: PBZP_{\mathrm{BZ}} — power extracted from the hole; Φ\Phi — flux threading the horizon; ΩH=a∗c32GM(1+1−a∗2)\Omega_H=\dfrac{a_* c^3}{2GM\left(1+\sqrt{1-a_*^2}\right)} — horizon angular velocity; κ≈0.05\kappa\approx 0.05 — geometric constant; a∗=Jc/GM2a_*=Jc/GM^2 — spin. Reading it: power grows as spin2^2 × flux2^2, so a fast spinner fed by a magnetically arrested disc (maximal Φ\Phi) is the most powerful engine. The disc-driven alternative is captured by a geometric launch criterion for a cold bead threaded on a rigidly rotating field line.

θfield<60∘ ⟹ magnetocentrifugal launch(5)\theta_{\mathrm{field}}<60^\circ\ \Longrightarrow\ \text{magnetocentrifugal launch} \qquad (5)

where: θfield\theta_{\mathrm{field}} — angle between the poloidal field line and the disc surface. Reading it: below 60∘60^\circ the potential along the wire tips outward and matter is slung off the disc; above it, no wind launches. Once matter is moving, its radiation is synchrotron: each electron of random Lorentz factor γ\gamma radiates near a critical frequency νc=32γ2νBsin⁡αp\nu_c=\tfrac{3}{2}\gamma^2\nu_B\sin\alpha_p, with νB=eB/2πmec\nu_B=eB/2\pi m_e c. Summing over a power-law population — the product of diffusive shock acceleration — gives a power-law spectrum whose slope encodes the electron slope.

N(γ) dγ∝γ−s dγ ⟹ Sν∝ν−(s−1)/2(6)N(\gamma)\,\mathrm{d}\gamma\propto\gamma^{-s}\,\mathrm{d}\gamma\ \Longrightarrow\ S_\nu\propto\nu^{-(s-1)/2} \qquad (6)

where: N(γ)N(\gamma) — electron energy distribution; ss — electron power-law index (shocks give s≈2s\approx 2–2.52.5); νB\nu_B — cyclotron frequency; α=−(s−1)/2\alpha=-(s-1)/2 — optically thin spectral index; the electron γ\gamma is distinct from the bulk Γ\Gamma. Measure the radio slope, read off the electrons that made it. Finally, two scaling laws anchor the energetics: the Eddington luminosity sets the accretion power scale, and the fundamental plane ties jet output to accretion across the black-hole mass range.

LEdd=4πGMmpcσT≈1.26×1038(MM⊙) erg s−1,LR∝M0.6–0.8 LX0.6(7)L_{\mathrm{Edd}}=\frac{4\pi G M m_p c}{\sigma_T}\approx 1.26\times10^{38}\left(\frac{M}{M_\odot}\right)\,\mathrm{erg\,s^{-1}},\qquad L_R\propto M^{0.6\text{–}0.8}\,L_X^{0.6} \qquad (7)

where: LEddL_{\mathrm{Edd}} — luminosity where radiation pressure balances gravity; mpm_p — proton mass; σT=6.65×10−25 cm2\sigma_T=6.65\times10^{-25}\,\mathrm{cm^2} — Thomson cross section; LR,LXL_R,L_X — radio (jet) and X-ray (accretion) luminosities; MM — black-hole mass. Reading it: the fundamental plane spans ∼9\sim 9 decades in mass — a single mass-scalable jet engine.

Derivation ASuperluminal apparent speed and its maximum

The goal is to derive equation (1) from first principles and extract its most useful consequence — that the largest apparent transverse speed a jet of bulk Lorentz factor Γ\Gamma can display is exactly Γβ\Gamma\beta, reached where cos⁡θ=β\cos\theta=\beta. This inverts an illusion into a measurement: an observed superluminal speed becomes a firm lower bound on the true Lorentz factor of an unresolved flow.

Step 1. Consider a blob moving at speed βc\beta c at angle θ\theta to the line of sight. In a coordinate time interval Δte\Delta t_e between two emission events, the blob advances βc Δte\beta c\,\Delta t_e, of which the projected component transverse to the line of sight is Δx⊥=βc Δtesin⁡θ\Delta x_\perp=\beta c\,\Delta t_e\sin\theta. Why: only the transverse motion is swept out on the plane of the sky.

Step 2. But the second photon is emitted from a point closer to us by βc Δtecos⁡θ\beta c\,\Delta t_e\cos\theta, so the interval between arrivals is compressed to Δtobs=Δte(1−βcos⁡θ)\Delta t_{\mathrm{obs}}=\Delta t_e(1-\beta\cos\theta). Why: this light-travel-time compression is the entire trick — the blob partly chases its own signal, so events appear to happen faster than they do.

Step 3. The apparent transverse speed is displacement over observed time, βapp=Δx⊥/(c Δtobs)=βsin⁡θ/(1−βcos⁡θ)\beta_{\mathrm{app}}=\Delta x_\perp/(c\,\Delta t_{\mathrm{obs}})=\beta\sin\theta/(1-\beta\cos\theta), equation (1). Differentiate and set the numerator to zero: cos⁡θ(1−βcos⁡θ)−βsin⁡2θ=0\cos\theta(1-\beta\cos\theta)-\beta\sin^2\theta=0 reduces to cos⁡θ=β\cos\theta=\beta. Why: the angle maximizing the illusion is tied directly to the speed itself.

Step 4. At cos⁡θ=β\cos\theta=\beta we have sin⁡θ=1/Γ\sin\theta=1/\Gamma and 1−β2=1/Γ21-\beta^2=1/\Gamma^2, so βapp,max=β(1/Γ)/(1−β2)=Γβ\beta_{\mathrm{app,max}}=\beta(1/\Gamma)/(1-\beta^2)=\Gamma\beta. Inverting, any measured βapp\beta_{\mathrm{app}} forces Γ≥1+βapp2\Gamma\geq\sqrt{1+\beta_{\mathrm{app}}^2}. Why: a knot seen at 4c4c cannot come from a flow with Γ<17≈4.1\Gamma<\sqrt{17}\approx4.1 — superluminal motion is a Lorentz-factor speedometer.

Go deeper · Derivation B · 1 minBlandford–Znajek power from the membrane picture
Derivation BBlandford–Znajek power from the membrane picture

The aim is to reproduce the scaling PBZ∝Φ2ΩH2P_{\mathrm{BZ}}\propto\Phi^2\Omega_H^2 without solving the full force-free electrodynamics, by treating the horizon as a rotating conductor in the membrane paradigm — a surface with finite resistance driving current through an external load. This scaling is why black-hole spin, not accretion rate, is thought to set the power of the strongest jets, and why the engine runs even when almost nothing is accreting.

Step 1. A conductor of size ∼rH\sim r_H spinning at ΩH\Omega_H through a field BB develops a motional EMF around its rim of order E∼ΩHBrH2/c∼ΩHΦ/c\mathcal{E}\sim\Omega_H B r_H^2/c\sim\Omega_H\Phi/c, since Φ∼BrH2\Phi\sim B r_H^2. Why: frame dragging makes the horizon a rotating battery, and the flux measures how strongly it couples to the exterior.

Step 2. In the membrane paradigm the horizon has an effective surface resistance of order the impedance of free space, RH≈377 ΩR_H\approx 377\,\Omega. The jet is the external load; maximum power transfer occurs at impedance matching, so P∼E2/RHP\sim\mathcal{E}^2/R_H. Why: a real generator delivers most power into a matched load — imported to the horizon.

Step 3. Substituting E∼ΩHΦ/c\mathcal{E}\sim\Omega_H\Phi/c gives PBZ∼(κ/4πc) Φ2ΩH2P_{\mathrm{BZ}}\sim(\kappa/4\pi c)\,\Phi^2\Omega_H^2, quadratic in both flux and rotation. For small spin ΩH∝a∗\Omega_H\propto a_*, so PBZ∝a∗2P_{\mathrm{BZ}}\propto a_*^2. Why: the two knobs are the trapped flux and the hole's spin, and power goes as the square of each — so a magnetically arrested disc around a fast spinner is the optimal engine.

Go deeper · Derivation C · 1 minSynchrotron spectral index from a power-law electron population
Derivation CSynchrotron spectral index from a power-law electron population

The purpose is to connect a microphysical input — the electron energy distribution N(γ)∝γ−sN(\gamma)\propto\gamma^{-s} — to an observable output, the slope of the optically thin radio spectrum, showing they are locked by α=−(s−1)/2\alpha=-(s-1)/2. This makes jet spectra diagnostic: a measured slope reads out the shock acceleration physics, and the canonical s≈2.5s\approx 2.5 predicts the steep α≈−0.75\alpha\approx-0.75 seen in lobes.

Step 1. A single electron of Lorentz factor γ\gamma radiates total power Psync=43σTc γ2β2UBP_{\mathrm{sync}}=\tfrac{4}{3}\sigma_T c\,\gamma^2\beta^2 U_B with UB=B2/8πU_B=B^2/8\pi, concentrated near νc∝γ2B\nu_c\propto\gamma^2 B. Why: each electron is effectively monochromatic, so "which electrons" and "which frequency" are one question.

Step 2. Invert the critical-frequency relation to write electron energy in terms of emitted frequency: γ∝(ν/B)1/2\gamma\propto(\nu/B)^{1/2}, hence dγ∝ν−1/2 dν\mathrm{d}\gamma\propto\nu^{-1/2}\,\mathrm{d}\nu. Why: building a spectrum means changing variables from electron energy to frequency, and this is the Jacobian.

Step 3. The emissivity at ν\nu is the number of electrons emitting there times the power each contributes, Sν dν∝N(γ) Psync dγ∝γ−s⋅γ2 dγS_\nu\,\mathrm{d}\nu\propto N(\gamma)\,P_{\mathrm{sync}}\,\mathrm{d}\gamma\propto\gamma^{-s}\cdot\gamma^2\,\mathrm{d}\gamma. Substituting γ∝ν1/2\gamma\propto\nu^{1/2} gives Sν∝ν(2−s)/2⋅ν−1/2=ν−(s−1)/2S_\nu\propto\nu^{(2-s)/2}\cdot\nu^{-1/2}=\nu^{-(s-1)/2}. Why: collecting powers of ν\nu yields α=−(s−1)/2\alpha=-(s-1)/2, so s=2.5⇒α=−0.75s=2.5\Rightarrow\alpha=-0.75.

Limiting cases

Jet phenomenology is controlled by three parameters — the bulk Lorentz factor Γ\Gamma, the viewing angle θ\theta, and the magnetization σ\sigma. The regimes below are the corners worth memorizing.

CaseResultMeaning
Poynting-dominated launchσ=B24πρc2≫1\sigma=\dfrac{B^2}{4\pi\rho c^2}\gg 1Near the base electromagnetic energy dwarfs rest mass; the flow is force-free and accelerates by spending σ\sigma toward unity.
Superluminal maximumβapp,max=Γβ≈Γ (Γ≫1)\beta_{\mathrm{app,max}}=\Gamma\beta\approx\Gamma\ (\Gamma\gg1)Reached at cos⁡θ=β\cos\theta=\beta; a Γ=5\Gamma=5 jet can appear to move at ∼5c\sim5c, and the apparent speed bounds Γ\Gamma.
Head-on beamingδ→2Γ (θ→0)\delta\to 2\Gamma\ (\theta\to0)Blazars reach maximal Doppler factor and boost flux by δ3–4\delta^{3\text{–}4} — a selection effect dominating flux-limited surveys.
Counter-jetδ→12Γ (θ→π)\delta\to\dfrac{1}{2\Gamma}\ (\theta\to\pi)De-boosted by the same power, brightness ratios >104>10^4 — why images look one-sided though the source is bipolar.
Maximal Kerr spinErot=(1−12)Mc2≈0.29 Mc2E_{\mathrm{rot}}=\left(1-\tfrac{1}{\sqrt2}\right)Mc^2\approx0.29\,Mc^2For a∗=1a_*=1 up to 29%29\% of Mc2Mc^2 is extractable — far exceeding radiative accretion efficiency.

Worked examples

1 · GRS 1915+105: solving a two-sided superluminal source (fully worked)

GRS 1915+105 ejects a pair of radio blobs. VLBI measures the approaching component at proper motion μa=17.6\mu_a=17.6 mas/day and the receding one at μr=9.0\mu_r=9.0 mas/day, at distance D=12.5D=12.5 kpc. Solve for the true speed β\beta, viewing angle θ\theta, and bulk Lorentz factor Γ\Gamma.

Step 1. Convert proper motion to apparent transverse speed, βapp=μD/c\beta_{\mathrm{app}}=\mu D/c. With 1 mas=4.848×10−91\,\mathrm{mas}=4.848\times10^{-9} rad, D=12.5×3.086×1021D=12.5\times3.086\times10^{21} cm, and a day =86400=86400 s, the approaching blob gives βapp,a≈1.27\beta_{\mathrm{app},a}\approx\mathbf{1.27} — superluminal.

Step 2. Because both signs of motion are seen, the geometry is fully constrained. The proper-motion ratio removes the distance and gives βcos⁡θ=μa−μrμa+μr=17.6−9.017.6+9.0=0.323\beta\cos\theta=\dfrac{\mu_a-\mu_r}{\mu_a+\mu_r}=\dfrac{17.6-9.0}{17.6+9.0}=\mathbf{0.323}.

Step 3. From the approaching relation βapp,a=βsin⁡θ/(1−βcos⁡θ)\beta_{\mathrm{app},a}=\beta\sin\theta/(1-\beta\cos\theta), solve for the transverse component: βsin⁡θ=βapp,a(1−βcos⁡θ)=1.27×(1−0.323)=0.860\beta\sin\theta=\beta_{\mathrm{app},a}(1-\beta\cos\theta)=1.27\times(1-0.323)=\mathbf{0.860}.

Step 4. Combine the two components: β=(βcos⁡θ)2+(βsin⁡θ)2=0.3232+0.8602=0.92\beta=\sqrt{(\beta\cos\theta)^2+(\beta\sin\theta)^2}=\sqrt{0.323^2+0.860^2}=\mathbf{0.92}, and tan⁡θ=0.860/0.323⇒θ≈69∘\tan\theta=0.860/0.323\Rightarrow\theta\approx\mathbf{69^\circ}. The Lorentz factor is Γ=1/1−0.922≈2.5\Gamma=1/\sqrt{1-0.92^2}\approx\mathbf{2.5}. So GRS 1915+105 launches a genuinely relativistic jet (β≈0.92\beta\approx0.92, Γ≈2.5\Gamma\approx2.5) viewed off-axis at θ≈69∘\theta\approx69^\circ — a scaled-down quasar in our Galaxy.

2 · Minimum Lorentz factor and the vanishing counter-jet

A quasar knot is observed at apparent speed 4c4c. (a) What is the smallest bulk Lorentz factor consistent with this? (b) If instead Γ=5\Gamma=5 at θ=20∘\theta=20^\circ with α=−0.7\alpha=-0.7, what is the jet-to-counter-jet brightness ratio for a continuous jet?

Step 1. Since βapp≤βapp,max=Γβ<Γ\beta_{\mathrm{app}}\leq\beta_{\mathrm{app,max}}=\Gamma\beta<\Gamma, we need Γ≥1+βapp2=1+16=17≈4.1\Gamma\geq\sqrt{1+\beta_{\mathrm{app}}^2}=\sqrt{1+16}=\sqrt{17}\approx\mathbf{4.1}, implying an intrinsic speed of at least β≈0.97\beta\approx0.97. A 4c4c image cannot be made by anything slower than Γ≈4.1\Gamma\approx4.1.

Step 2. For part (b), with Γ=5\Gamma=5 the speed is β=1−1/25=0.980\beta=\sqrt{1-1/25}=0.980, and βcos⁡20∘=0.980×0.940=0.921\beta\cos20^\circ=0.980\times0.940=0.921. The jet/counter-jet ratio is R=[1+βcos⁡θ1−βcos⁡θ]pR=\left[\dfrac{1+\beta\cos\theta}{1-\beta\cos\theta}\right]^{p} with p=2−α=2.7p=2-\alpha=2.7 for a continuous jet.

Step 3. Evaluate: (1+0.921)/(1−0.921)=24.3(1+0.921)/(1-0.921)=24.3, so R=24.32.7≈5.5×103R=24.3^{2.7}\approx\mathbf{5.5\times10^{3}} (for discrete blobs, p=3.7p=3.7, it reaches ∼105\sim10^{5}). The counter-jet is thousands of times fainter — undetectable in practice — which is why relativistic jets nearly always appear one-sided.

3 · The spin reservoir of the M87 black hole

The M87 black hole has mass M≈6.5×109 M⊙M\approx6.5\times10^{9}\,M_\odot. If it spins near maximally, how much rotational energy is extractable, versus what a ∼1044 erg s−1\sim10^{44}\,\mathrm{erg\,s^{-1}} jet radiates over 10810^{8} years?

Step 1. Rest-mass energy: Mc2≈6.5×109×1.99×1033 g×(3.0×1010)2≈1.2×1064 ergMc^2\approx6.5\times10^{9}\times1.99\times10^{33}\,\mathrm{g}\times(3.0\times10^{10})^2\approx\mathbf{1.2\times10^{64}\ erg}. The extractable fraction is 1−1/2≈0.291-1/\sqrt2\approx0.29, so Erot≈3.4×1063 ergE_{\mathrm{rot}}\approx\mathbf{3.4\times10^{63}\ erg}.

Step 2. The jet expends 1044 erg s−1×108 yr≈3.2×1059 erg10^{44}\,\mathrm{erg\,s^{-1}}\times10^{8}\,\mathrm{yr}\approx3.2\times10^{59}\,\mathrm{erg}. The spin reservoir exceeds this by ∼104\sim10^{4} — spin-down alone can power the jet for essentially the age of the universe, which is why Blandford–Znajek is the favored AGN engine.

Physics of the portrait

Relativistic jets are named for two theorists whose 1977 and 1982 papers still frame every discussion of how jets launch; the portrait belongs to the one whose name is on both.

Origins · 1 minRoger Blandford and how this was discovered
RB
Roger Blandford
b. 1949 · magnetohydrodynamic launching of relativistic jets

Roger Blandford is a British-American theoretical astrophysicist whose work around 1977–1982 turned jet launching from a puzzle into a pair of concrete electrodynamic mechanisms. The challenge was that quasars channel prodigious power into thin beams collimated across millions of light-years, and neither thermal nor radiative pressure could plausibly both accelerate matter to relativistic speeds and hold it together over such reach. In 1977, with his student Roman Znajek, Blandford showed that a rotating Kerr black hole threaded by a magnetic field acts like a spinning battery: frame dragging drives currents through the horizon and extracts the hole's rotational energy as Poynting flux with power ∝Φ2ΩH2\propto\Phi^2\Omega_H^2 — needing no accretion at all, with up to 29%29\% of a maximal hole's mass-energy available. In 1982, with David Payne, he derived the complementary disc mechanism: field lines anchored in a rotating disc, inclined under 60∘60^\circ to its plane, fling matter outward as a self-collimating magnetocentrifugal wind. The two mechanisms bearing his name remain the twin pillars of jet theory, now tested against Event Horizon Telescope images of the very horizon his 1977 paper energized.

History
1918
— Curtis notes a curious straight ray emerging from the nucleus of the galaxy M87
1971
— VLBI reveals apparent superluminal motion in quasars, at odds with naive causality
1977
— Blandford and Znajek show a spinning Kerr hole can be spun down electromagnetically
1982
— Blandford and Payne derive the magnetocentrifugal disc-wind mechanism
1994
— Mirabel and Rodriguez discover superluminal jets in the Galactic microquasar GRS 1915+105

Going further

Texts. Carroll & Ostlie, An Introduction to Modern Astrophysics, 2nd ed., ch. 17–18 (X-ray binaries and accretion) and ch. 28 (active galaxies and jets) · Binney & Tremaine, Galactic Dynamics, 2nd ed. (galactic-scale outflows and dynamics) · Rybicki & Lightman, Radiative Processes in Astrophysics (synchrotron and inverse-Compton emission) · Blandford & Znajek (1977), MNRAS 179, 433 (electromagnetic extraction of energy from Kerr holes) · Blandford & Payne (1982), MNRAS 199, 883 (magnetocentrifugal disc winds) · Mirabel & Rodriguez (1999), ARA&A 37, 409 (relativistic jets in the Galaxy).

Provenance. The superluminal formula and its maximum Γβ\Gamma\beta are derived by differentiation (Carroll & Ostlie ch. 28); the Doppler factor and beaming exponents follow from relativistic aberration (Rybicki & Lightman); the synchrotron index α=−(s−1)/2\alpha=-(s-1)/2 from N(γ)∝γ−sN(\gamma)\propto\gamma^{-s} and νc∝γ2B\nu_c\propto\gamma^2 B. The Eddington luminosity, σT=6.65×10−25 cm2\sigma_T=6.65\times10^{-25}\,\mathrm{cm^2}, and the Kerr fraction 1−1/2≈0.291-1/\sqrt2\approx0.29 are textbook results; the Blandford–Znajek and Blandford–Payne scalings follow the original papers. GRS 1915+105 parameters (μa=17.6\mu_a=17.6, μr=9.0\mu_r=9.0 mas/day, D=12.5D=12.5 kpc ⇒β≈0.92\Rightarrow\beta\approx0.92, θ≈69∘\theta\approx69^\circ) are from Mirabel & Rodriguez (1994, 1999); all worked-example numbers are re-verified here.

Prerequisites. Neutron Stars · Pulsar Magnetospheres and Spin-Down · Magnetars. Leads to. Millisecond and Recycled Pulsars · Equation of State of Dense Matter.

Ask the tutor: Ask me to walk through the two-sided proper-motion solve for a microquasar, or why a magnetically arrested disc maximizes Blandford–Znajek power.

Why doesn't superluminal motion violate special relativity?

Because nothing physical moves faster than light — only the source's image does. A blob approaching at 0.98c0.98c at a small angle nearly keeps pace with the photons it emits, so light from successive positions arrives bunched in time. The projected displacement divided by this compressed interval can exceed cc even though the true speed is safely below it. No signal outruns light; the illusion lives in the light-travel-time geometry captured by 1−βcos⁡θ1-\beta\cos\theta, and it actually lets us measure the true Lorentz factor rather than contradicting relativity.

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

They tap different reservoirs. Blandford–Znajek extracts the rotational energy of the black hole itself through field lines threading the horizon; frame dragging drives currents that carry Poynting flux outward, and it works even with negligible accretion. Blandford–Payne is a magnetocentrifugal disc wind: field lines anchored in the rotating disc, inclined under 60∘60^\circ to its plane, fling matter outward like beads on a sling. Real jets probably combine both — a fast Blandford–Znajek spine in a slower Blandford–Payne sheath. The discriminant is spin: the most powerful jets require fast spinners, favoring Blandford–Znajek as the dominant AGN engine.

What actually keeps a jet collimated over millions of light-years?

Magnetic hoop stress does most of the work. As the flow rotates, the field winds into a tight toroidal spiral, and the tension in those loops pinches inward toward the axis — the same pinch that constricts a laboratory Z-pinch plasma, with external pressure from surrounding gas adding confinement at large scales. Pure thermal or radiative pressure cannot do this: a gas-pressure-driven outflow expands freely and loses collimation quickly. The persistence of thin, straight jets across galaxy-cluster scales is one of the strongest arguments that jets are magnetically, not thermally, driven.

Why do black holes make more powerful jets than neutron stars?

Two reasons, both tied to the center. First, only a black hole has a spinning horizon the Blandford–Znajek mechanism can tap; a neutron star has a solid surface and rigidly rotating magnetosphere instead, closing off the most powerful channel. Second, the horizon lets a magnetically arrested disc trap enormous flux without it being expelled or mass-loaded by a surface. Accreting neutron stars do launch jets, but at a given X-ray luminosity they are weaker in radio — consistent with the spin channel switched off, making the two populations a direct test of how essential the horizon is.

What is a microquasar, and how does it relate to a quasar?

A microquasar is a stellar-mass black hole (or neutron star) in an X-ray binary that launches relativistic radio jets — a scaled-down quasar. The physics is thought identical up to the central mass: swap a ∼10 M⊙\sim10\,M_\odot hole for a ∼109 M⊙\sim10^{9}\,M_\odot one and the jet grows from light-years to megaparsecs, timescales stretching by ∼108\sim10^{8}. This scale invariance is powerful because microquasars vary on human timescales — a state transition and jet ejection play out in days, not the millennia an AGN needs — so we watch the jet–disc connection in real time, and the fundamental plane makes it quantitative across the mass range.

How do we know the radiation is synchrotron rather than thermal?

Three fingerprints. First, the spectrum is a power law, Sν∝ναS_\nu\propto\nu^{\alpha}, over many decades — nothing like a peaked thermal blackbody. Second, the emission is strongly polarized, sometimes tens of percent, which thermal emission is not, and it maps the projected field geometry. Third, the optically thin index α≈−0.7\alpha\approx-0.7 matches α=−(s−1)/2\alpha=-(s-1)/2 for a shock-accelerated population with s≈2.4s\approx2.4. Together these show the light comes from relativistic electrons gyrating in an ordered field — the same field that launches and collimates the flow.

Before this

Where this leads