Superluminal Motion in Jets

Point a very-long-baseline radio array at the heart of a quasar and, over a few years, you can watch it move. The parsec-scale jet resolves into a stationary "core" and bright "knots" that drift outward, and converting their angular proper motion into a transverse velocity using the source's distance gives a number routinely several times the speed of light — close to forty times cc in the most extreme sources. This is superluminal motion, and after its discovery in 1971 it briefly looked like a crisis for relativity. It is not. The resolution needs no new physics: it is a pure light-travel-time projection effect for material streaming toward us at nearly cc. Because an approaching blob almost keeps pace with the photons it emits, successive flashes arrive bunched together, and the sky-projected displacement divided by that compressed interval can exceed cc though nothing physical does. Remarkably, Martin Rees wrote the prediction down in 1966, five years before anyone saw it. What the observations added was a measurement: a single apparent speed places a hard lower bound on the bulk Lorentz factor, γ≥1+βapp2\gamma\ge\sqrt{1+\beta_{\rm app}^2}, turning "are these jets relativistic?" into a settled yes. The same geometry that fakes faster-than-light motion also Doppler-boosts the approaching jet, dims its receding twin below detection, and swells the variability — so superluminal sources are the pole-on blazars of the AGN unified model. Apparent kinematics and observed luminosity are two readings of one orientation, and this article derives the effect from scratch, extracts the Lorentz factor and viewing angle from data, and traces the line from Rees to the MOJAVE survey.

Key takeaways
  • Compact jet knots in quasars and radio galaxies cross the sky at apparent transverse speeds βapp\beta_{\rm app} well above 1, sometimes reaching βapp∼10\beta_{\rm app}\sim 10–4040 in units of cc.
  • It is a geometric illusion: plasma at bulk speed β→1\beta\to 1 at a small angle θ\theta to the line of sight nearly overtakes its own light, compressing photon arrival times by (1−βcos⁡θ)(1-\beta\cos\theta).
  • The apparent speed is βapp=βsin⁡θ1−βcos⁡θ\beta_{\rm app}=\dfrac{\beta\sin\theta}{1-\beta\cos\theta}, peaking at cos⁡θ=β\cos\theta=\beta at βapp,max=βγ\beta_{\rm app,max}=\beta\gamma, with γ=(1−β2)−1/2\gamma=(1-\beta^2)^{-1/2} the bulk Lorentz factor.
  • A single measured βapp\beta_{\rm app} forces a hard bound γ≥1+βapp2\gamma\ge\sqrt{1+\beta_{\rm app}^2} — kinematic proof that AGN jets carry highly relativistic plasma.
  • The identical geometry drives relativistic Doppler beaming, tying superluminal sources to blazars and the orientation-based AGN unified model.
  • Key idea: an emitter moving at β→1\beta\to 1 at a small angle θ\theta to the line of sight nearly catches its own emitted light, so the observed interval between two positions is compressed by (1−βcos⁡θ)(1-\beta\cos\theta); the sky-projected speed βapp=βsin⁡θ/(1−βcos⁡θ)\beta_{\rm app}=\beta\sin\theta/(1-\beta\cos\theta) can then exceed unity though the true speed never does.
The core insight

An emitter moving at β→1\beta\to 1 at a small angle θ\theta to the line of sight nearly overtakes its own emitted light, so the observed time between two positions is compressed by the factor (1−βcos⁡θ)(1-\beta\cos\theta); the sky-projected speed βapp=βsin⁡θ/(1−βcos⁡θ)\beta_{\rm app}=\beta\sin\theta/(1-\beta\cos\theta) can then exceed unity even though the true speed cannot.

Where this lives today

The modern engine of the field is systematic, multi-epoch VLBI monitoring. The MOJAVE program has tracked hundreds of AGN jets at 15 GHz for over two decades, measuring proper motions, accelerations, and knot bending. Its apparent-speed distribution — most sources modestly superluminal, a tail reaching βapp∼40\beta_{\rm app}\sim40 — is just what an orientation-selected population with intrinsic γ∼5\gamma\sim5–4040 should produce, and fitting it with Doppler factors yields the speed and orientation statistics that anchor the AGN unified model. Superluminal motion is now a routine diagnostic: with a Doppler factor it delivers γ\gamma and θ\theta (worked example 2), and its mere presence certifies a jet as relativistic. The Event Horizon Telescope connects the jet base in M87 near the black hole to the superluminal knots further out, while Fermi links the fastest blazar jets to the flaring gamma-ray sky. Open questions remain sharp — how and where jets accelerate to these Lorentz factors, and whether a structured "spine–sheath" flow shows up in the kinematics — but the backbone, Rees's 1966 projection effect, is settled physics.

The physical picture

Picture a blob that emits a flash at position A, then moves toward you and emits a second at B. Because it has travelled closer to you in between, the second flash has less distance to cover, so the two reach you separated by less time than the true flight time from A to B. Divide the real transverse displacement by that compressed interval and the quotient — the apparent sky speed — can exceed cc. Nothing outran light; the messenger simply ran alongside the source. The first figure fixes the geometry and reads the formula off directly.

Figure
A blob moves at speed v at angle θ to the line of sight. Between two flashes it advances Δz = vΔt cosθ toward the observer, so the second flash arrives only Δt_obs = Δt(1 − βcosθ) later, while it sweeps Δx⊥ = vΔt sinθ across the sky — hence β_app = βsinθ ⁄ (1 − βcosθ). Angle exaggerated for clarity.
The light-travel-time projection: the second flash is emitted Δz closer to you line of sight observer v = βc θ A flash 1 B flash 2 shorter photon path Δx⊥ = vΔt sinθ (across the sky) Δz = vΔt cosθ toward observer → arrival gap shrinks to Δt(1 − βcosθ)

Two competing factors are the whole story. The numerator's sin⁡θ\sin\theta vanishes head-on — a blob aimed straight at you shows no transverse motion. The denominator's 1/(1−βcos⁡θ)1/(1-\beta\cos\theta) blows up head-on, where compression is strongest. Their product peaks at the "sweet spot" cos⁡θ=β\cos\theta=\beta, where βapp=βγ≈γ\beta_{\rm app}=\beta\gamma\approx\gamma. The apparent speed is thus a Lorentz-factor detector: to fake a large βapp\beta_{\rm app} the flow must genuinely move at γ≳βapp\gamma\gtrsim\beta_{\rm app}. The second figure plots this and marks the peaks.

Figure
Apparent transverse speed β_app versus viewing angle θ for bulk Lorentz factors γ = 5 and γ = 10. Each curve climbs past the β_app = 1 threshold, peaks at cosθ = β where β_app = βγ (9.95 at θ = 5.7° for γ = 10, and 4.90 at θ = 11.5° for γ = 5), then falls back to β near θ = 90°. The peak apparent speed is essentially the Lorentz factor itself.
The punchline: β_app(θ) peaks at cosθ = β, where β_app = βγ ≈ γ β_app = 1 (the speed of light) 0 2 4 6 8 10 0 15 30 45 60 90 viewing angle θ (degrees) apparent speed β_app γ = 10: β_app,max = 9.95 at θ = 5.7° γ = 5: β_app,max = 4.90 at θ = 11.5°

The amber curve (γ=5\gamma=5) tops out near 4.9, the blue curve (γ=10\gamma=10) near 9.95: doubling the Lorentz factor doubles the peak apparent speed while shrinking the angle at which it occurs, since sin⁡θmax=1/γ\sin\theta_{\rm max}=1/\gamma. Fast jets beam their fastest motion into an ever-narrower cone — which is why the superluminal sky is dominated by the rare few pointed almost straight at us.

The core equations

The master relation follows from the first figure: transverse displacement over light-compressed observed time. Everything downstream — maximum, Lorentz bound, Doppler links — is algebra on this one expression.

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

where: βapp\beta_{\rm app} — apparent transverse speed in units of cc; β=v/c\beta=v/c — true bulk speed; θ\theta — angle between the velocity and the line of sight. Reading it: the denominator is the compression factor, always below 1 for an approaching source, and it pushes the ratio above unity even though the numerator βsin⁡θ\beta\sin\theta is itself below 1. Maximising equation (1) over θ\theta at fixed β\beta gives a clean optimum: the largest apparent speed a jet of Lorentz factor γ\gamma can display is βγ\beta\gamma, when it points a little off the line of sight.

βapp,max=βγ,atcos⁡θ=β,  sin⁡θ=1γ(2)\beta_{\rm app,max}=\beta\gamma,\qquad \text{at}\quad\cos\theta=\beta,\ \ \sin\theta=\frac{1}{\gamma} \qquad (2)

where: γ=(1−β2)−1/2\gamma=(1-\beta^2)^{-1/2} — bulk Lorentz factor. Reading it: for a relativistic jet β→1\beta\to1, so βapp,max≈γ\beta_{\rm app,max}\approx\gamma — the peak apparent speed reads off the Lorentz factor, at an optimal angle of roughly 1/γ1/\gamma radians. Because the observed value can never exceed that maximum, a single βapp\beta_{\rm app} — viewing angle unknown — already forces a floor on γ\gamma, the most model-independent statement in the subject.

γ ≥ 1+βapp2(3)\gamma\ \ge\ \sqrt{1+\beta_{\rm app}^2} \qquad (3)

where: βapp\beta_{\rm app} — measured apparent speed; γ\gamma — minimum Lorentz factor consistent with it. Reading it: βapp=10\beta_{\rm app}=10 demands γ≥101≈10.05\gamma\ge\sqrt{101}\approx10.05; no angle lets a slow jet fake a fast one. Turning a raw observation into βapp\beta_{\rm app} needs cosmology: an interferometer measures an angular proper motion μ\mu, not a velocity, and distance plus a redshift factor convert it.

μ=βapp cDA (1+z)(4)\mu=\frac{\beta_{\rm app}\,c}{D_A\,(1+z)} \qquad (4)

where: μ\mu — angular proper motion (rad s−1^{-1}); DAD_A — angular-diameter distance; zz — redshift, whose (1+z)(1+z) converts observed to emitted time. Reading it: two sources at equal true speed but different zz show different μ\mu, so a cosmological model is unavoidable when quoting βapp\beta_{\rm app}. The same geometry also boosts the flux, through the relativistic Doppler factor — the same denominator as equation (1), now controlling brightness.

δ=1γ (1−βcos⁡θ)(5)\delta=\frac{1}{\gamma\,(1-\beta\cos\theta)} \qquad (5)

where: δ\delta — relativistic Doppler factor. Reading it: near the optimal angle δ≈γ\delta\approx\gamma, so a γ=15\gamma=15 jet is boosted by δ∼15\delta\sim15 in frequency and time, and by a large power in flux. Both βapp\beta_{\rm app} and δ\delta are functions of the same (β,θ)(\beta,\theta), so they satisfy exact identities; measuring both pins down γ\gamma and θ\theta.

βapp=γ δ βsin⁡θ,βapp2+δ2=2γδ−1(6)\beta_{\rm app}=\gamma\,\delta\,\beta\sin\theta,\qquad \beta_{\rm app}^2+\delta^2=2\gamma\delta-1 \qquad (6)

where: δ\delta — Doppler factor; βapp\beta_{\rm app} — apparent speed. Reading it: the second identity inverts to γ=(βapp2+δ2+1)/(2δ)\gamma=(\beta_{\rm app}^2+\delta^2+1)/(2\delta), from which β\beta and θ\theta follow — the standard route to full jet kinematics. (This corrects a common misprint: the left identity is βapp=γδβsin⁡θ\beta_{\rm app}=\gamma\delta\beta\sin\theta, since γδ=1/(1−βcos⁡θ)\gamma\delta=1/(1-\beta\cos\theta).) The flux boost scales as a power of the Doppler factor, the exponent set by geometry and the synchrotron spectrum.

Sobs=δ p Sem,p=3−α  (discrete blob),p=2−α  (continuous jet)(7)S_{\rm obs}=\delta^{\,p}\,S_{\rm em},\qquad p=3-\alpha\ \ (\text{discrete blob}),\quad p=2-\alpha\ \ (\text{continuous jet}) \qquad (7)

where: Sobs,SemS_{\rm obs},S_{\rm em} — observed and rest-frame flux densities; α\alpha — spectral index, Sν∝ναS_\nu\propto\nu^{\alpha}. Reading it: a synchrotron jet with α≈−0.7\alpha\approx-0.7 has p≈2.7p\approx2.7–3.73.7, so even modest δ\delta boosts enormously. Applied to a symmetric jet's two sides, equation (7) predicts a colossal asymmetry — why parsec-scale jets look one-sided.

SjetScounterjet=(1+βcos⁡θ1−βcos⁡θ)p(8)\frac{S_{\rm jet}}{S_{\rm counterjet}}=\left(\frac{1+\beta\cos\theta}{1-\beta\cos\theta}\right)^{p} \qquad (8)

where: Sjet,ScounterjetS_{\rm jet},S_{\rm counterjet} — flux of approaching versus receding jet. Reading it: for β≈1\beta\approx1, small θ\theta this reaches millions, so the counterjet drops below detection — the geometry that fakes superluminal motion also hides the far side.

Derivation AThe superluminal formula from light-travel time

The goal is to derive the master relation βapp=βsin⁡θ/(1−βcos⁡θ)\beta_{\rm app}=\beta\sin\theta/(1-\beta\cos\theta) from first principles, using nothing but the geometry of two emission events and the finite speed of light — no velocity addition, no Lorentz transformation, just bookkeeping of when photons leave the blob and when they reach us. We compare the true transverse motion to the compressed interval over which we receive the light.

Step 1. The blob emits a flash at coordinate time t1=0t_1=0 and a second at t2=Δtt_2=\Delta t, moving a coordinate distance v Δtv\,\Delta t in between. Why: two events define a motion, and coordinate (source-frame) time keeps the light-travel effect separate from any intrinsic clock effect.

Step 2. Resolve onto the line of sight: the blob advances Δz=v Δtcos⁡θ\Delta z=v\,\Delta t\cos\theta toward the observer between flashes. Why: this approach shortens the second photon's path and so shrinks the interval over which the flashes are received.

Step 3. The second photon starts Δz\Delta z nearer, so it needs Δz/c\Delta z/c less travel time; the received interval is Δtobs=Δt−Δz/c=Δt (1−βcos⁡θ)\Delta t_{\rm obs}=\Delta t-\Delta z/c=\Delta t\,(1-\beta\cos\theta). Why: this single factor (1−βcos⁡θ)(1-\beta\cos\theta) is the whole physics of the illusion — as β→1\beta\to1, θ→0\theta\to0 it drives the observed interval toward zero.

Step 4. Resolve perpendicular to the line of sight: the transverse motion an interferometer measures is Δx⊥=v Δtsin⁡θ\Delta x_\perp=v\,\Delta t\sin\theta. Why: only motion across the sky gives the changing angular position we call proper motion; line-of-sight motion shows up only through the Step 3 timing effect.

Step 5. Divide the transverse displacement by the received interval and normalise by cc: βapp=Δx⊥/(c Δtobs)=βsin⁡θ/(1−βcos⁡θ)\beta_{\rm app}=\Delta x_\perp/(c\,\Delta t_{\rm obs})=\beta\sin\theta/(1-\beta\cos\theta). Why: this is the ratio an observer forms, and it exceeds 1 whenever βsin⁡θ>1−βcos⁡θ\beta\sin\theta>1-\beta\cos\theta — easily met by fast jets at small angles. *Note: βapp>1\beta_{\rm app}>1 requires only β>1/(sin⁡θ+cos⁡θ)\beta>1/(\sin\theta+\cos\theta) — for example, β>0.71\beta>0.71 at θ=45°\theta=45° — so even modestly relativistic jets can appear superluminal at favourable angles.*

Go deeper · Derivation B · 1 minThe optimal viewing angle and the peak apparent speed
Derivation BThe optimal viewing angle and the peak apparent speed

The goal is to find the viewing angle that maximises the apparent speed at fixed bulk speed, and to evaluate the maximum there. This matters because the largest apparent speeds across a jet population trace an envelope set by the true Lorentz factors, so knowing where and how high the curve peaks lets us read γ\gamma almost directly off the fastest sources. We differentiate βapp(θ)\beta_{\rm app}(\theta) at fixed β\beta, solve for the stationary point, and substitute back.

Step 1. Differentiate βapp(θ)=βsin⁡θ/(1−βcos⁡θ)\beta_{\rm app}(\theta)=\beta\sin\theta/(1-\beta\cos\theta) by the quotient rule, giving the numerator βcos⁡θ (1−βcos⁡θ)−β2sin⁡2θ\beta\cos\theta\,(1-\beta\cos\theta)-\beta^2\sin^2\theta over (1−βcos⁡θ)2(1-\beta\cos\theta)^2. Why: the maximum is a stationary point, so setting the numerator to zero locates it.

Step 2. Expand: βcos⁡θ−β2cos⁡2θ−β2sin⁡2θ=β(cos⁡θ−β)\beta\cos\theta-\beta^2\cos^2\theta-\beta^2\sin^2\theta=\beta(\cos\theta-\beta), which vanishes at cos⁡θmax=β\cos\theta_{\rm max}=\beta. Why: the Pythagorean identity collapses the trigonometric terms into a constant β2\beta^2, so the optimum is startlingly simple — the cosine of the best angle equals the speed itself.

Step 3. At the optimum, sin⁡θmax=1−β2=1/γ\sin\theta_{\rm max}=\sqrt{1-\beta^2}=1/\gamma and 1−βcos⁡θmax=1−β2=1/γ21-\beta\cos\theta_{\rm max}=1-\beta^2=1/\gamma^2. Why: expressing both through γ\gamma makes the substitution a one-line cancellation and shows the optimal angle shrinks as 1/γ1/\gamma.

Step 4. Substitute back: βapp,max=β(1/γ)/(1/γ2)=βγ\beta_{\rm app,max}=\beta(1/\gamma)/(1/\gamma^2)=\beta\gamma. Why: for γ≫1\gamma\gg1 this is ≈γ\approx\gamma, so the peak apparent speed is essentially the Lorentz factor — the key to interpreting the fastest sources. *The calculus confirms the textbook maximum: γ=10\gamma=10 gives βapp,max=9.95\beta_{\rm app,max}=9.95 at θ=5.7°\theta=5.7°, and γ=5\gamma=5 gives 4.904.90 at θ=11.5°\theta=11.5° (Figure 2).*

Go deeper · Derivation C · 1 minThe kinematic lower bound on the Lorentz factor
Derivation CThe kinematic lower bound on the Lorentz factor

The goal is to extract the most robust statement from a single apparent-speed measurement when the viewing angle is unknown: the smallest bulk Lorentz factor that could reproduce it. This is what makes superluminal motion evidence rather than curiosity — it converts one measured number into a hard floor on the jet's true speed, and historically it settled whether AGN jets are relativistic. We need only the peak relation of the previous derivation and one identity.

Step 1. The observed speed cannot exceed the maximum for the true Lorentz factor: βapp≤βapp,max=βγ\beta_{\rm app}\le\beta_{\rm app,max}=\beta\gamma. Why: whatever the unknown angle, equation (2) caps the apparent speed at βγ\beta\gamma, so the measurement constrains βγ\beta\gamma from below.

Step 2. Square and use β2γ2=(1−γ−2)γ2=γ2−1\beta^2\gamma^2=(1-\gamma^{-2})\gamma^2=\gamma^2-1, giving βapp2≤γ2−1\beta_{\rm app}^2\le\gamma^2-1. Why: the identity eliminates β\beta, leaving a clean inequality in γ\gamma alone — no need to know speed and angle separately.

Step 3. Solve: γ≥1+βapp2\gamma\ge\sqrt{1+\beta_{\rm app}^2}. Why: this bound means βapp≫1\beta_{\rm app}\gg1 forces γ≳βapp\gamma\gtrsim\beta_{\rm app}, proving the emitting matter carries relativistic bulk kinetic energy and is not a slow-moving pattern. *For example, βapp=8⇒γ≥65≈8.06\beta_{\rm app}=8\Rightarrow\gamma\ge\sqrt{65}\approx8.06 and βapp=10⇒γ≥10.05\beta_{\rm app}=10\Rightarrow\gamma\ge10.05, consistent with the γ∼5\gamma\sim5–4040 range seen across VLBI surveys.*

Limiting cases

Reading the master relation βapp=βsin⁡θ/(1−βcos⁡θ)\beta_{\rm app}=\beta\sin\theta/(1-\beta\cos\theta) in its extreme regimes shows both when the effect switches on and why it must vanish in the everyday limit.

CaseResultMeaning
Non-relativistic (β≪1\beta\ll1)βapp≈βsin⁡θ≪1\beta_{\rm app}\approx\beta\sin\theta\ll1Ordinary projection; the apparent speed is just the true transverse speed, always sub-luminal.
In the plane of the sky (θ=90∘\theta=90^\circ)βapp=β<1\beta_{\rm app}=\beta<1No line-of-sight motion, so no arrival-time compression and no enhancement.
Head-on (θ→0\theta\to0)βapp→0\beta_{\rm app}\to0Maximum time compression but zero transverse displacement; the knot barely shifts on the sky.
Optimal angle (cos⁡θ=β\cos\theta=\beta)βapp,max=βγ≈γ\beta_{\rm app,max}=\beta\gamma\approx\gammaLargest apparent speed for a given γ\gamma, at sin⁡θ=1/γ\sin\theta=1/\gamma.
Ultra-relativistic near optimumβapp≈γ≈δ\beta_{\rm app}\approx\gamma\approx\deltaApparent speed, Doppler factor, and Lorentz factor coincide — the blazar regime.

Worked examples

1 · From proper motion to Lorentz factor (fully worked)

A quasar at z=0.5z=0.5 has a jet knot drifting at μ=0.2\mu=0.2 mas yr−1^{-1}. In a standard flat cosmology (H0=70H_0=70, Ωm=0.3\Omega_m=0.3) the angular-diameter distance is DA≈1267D_A\approx1267 Mpc. Find the apparent speed and minimum Lorentz factor.

Step 1. Convert μ\mu to SI. One mas is 10−3/20626510^{-3}/206265 rad, one year 3.156×1073.156\times10^{7} s, so μ=0.2×10−3/206265/(3.156×107)=3.07×10−17\mu=0.2\times10^{-3}/206265/(3.156\times10^{7})=3.07\times10^{-17} rad s−1^{-1}.

Step 2. Distance: DA=1267×3.086×1022=3.91×1025D_A=1267\times3.086\times10^{22}=3.91\times10^{25} m.

Step 3. Equation (4): βapp=μDA(1+z)/c=(3.07×10−17)(3.91×1025)(1.5)/(2.998×108)=6.0\beta_{\rm app}=\mu D_A(1+z)/c=(3.07\times10^{-17})(3.91\times10^{25})(1.5)/(2.998\times10^{8})=\mathbf{6.0} — six times cc.

Step 4. Apply bound (3): γ≥1+36=37=6.08\gamma\ge\sqrt{1+36}=\sqrt{37}=\mathbf{6.08}. If the jet sits at the optimal angle the bound is tight, with sin⁡θ=1/γ=0.164\sin\theta=1/\gamma=0.164, θ≈9.4∘\theta\approx\mathbf{9.4^\circ}. Result: βapp=6.0\beta_{\rm app}=6.0, so the plasma moves at γ≥6.1\gamma\ge6.1 — genuinely relativistic.

2 · Solving for the jet's speed and angle (fully worked)

VLBI gives a blazar knot βapp=15\beta_{\rm app}=15, and variability plus the synchrotron self-Compton limit give a Doppler factor δ=12\delta=12. Find the bulk Lorentz factor, true speed, and viewing angle.

Step 1. Solve equation (6)'s second identity for γ\gamma: γ=(βapp2+δ2+1)/(2δ)=(225+144+1)/24=15.4\gamma=(\beta_{\rm app}^2+\delta^2+1)/(2\delta)=(225+144+1)/24=\mathbf{15.4}. One inversion delivers the Lorentz factor from the two observables.

Step 2. Recover the speed: β=1−1/γ2=0.9979\beta=\sqrt{1-1/\gamma^2}=\mathbf{0.9979}, i.e. 99.79%99.79\% of cc.

Step 3. Recover the angle from δ=1/[γ(1−βcos⁡θ)]\delta=1/[\gamma(1-\beta\cos\theta)]: 1−βcos⁡θ=1/(γδ)=1/185.01-\beta\cos\theta=1/(\gamma\delta)=1/185.0, so cos⁡θ=(1−1/185.0)/β=0.9967\cos\theta=(1-1/185.0)/\beta=0.9967, θ≈4.7∘\theta\approx\mathbf{4.7^\circ}. Just inside the optimal θmax=arcsin⁡(1/γ)≈3.7∘\theta_{\rm max}=\arcsin(1/\gamma)\approx3.7^\circ, so both observables are large at once.

Step 4. Check: 1−βcos⁡θ=0.005411-\beta\cos\theta=0.00541 gives βapp=15.0\beta_{\rm app}=15.0 ✓ and δ=1/(15.4×0.00541)=12.0\delta=1/(15.4\times0.00541)=12.0 ✓. Result: γ≈15.4\gamma\approx15.4, β≈0.998\beta\approx0.998, θ≈4.7∘\theta\approx4.7^\circ — a fast jet pointed almost at us.

3 · Why the counterjet disappears

Take the same jet, β=0.998\beta=0.998 at θ=5∘\theta=5^\circ, with a synchrotron spectrum α=−0.7\alpha=-0.7 so the continuous-jet exponent is p=2−α=2.7p=2-\alpha=2.7. Estimate the approaching-to-receding brightness ratio.

Step 1. With βcos⁡θ=0.998×0.99619=0.9942\beta\cos\theta=0.998\times0.99619=0.9942, the bracket is (1+βcos⁡θ)/(1−βcos⁡θ)=1.9942/0.00582=343(1+\beta\cos\theta)/(1-\beta\cos\theta)=1.9942/0.00582=\mathbf{343}.

Step 2. Raise to p=2.7p=2.7: Sjet/Scounterjet=3432.7≈7×106S_{\rm jet}/S_{\rm counterjet}=343^{2.7}\approx\mathbf{7\times10^{6}}. The approaching jet outshines its receding twin seven million to one — the counterjet is undetectable, which is why superluminal jets look one-sided. (For a discrete blob, p=3.7p=3.7 gives ∼2×109\sim2\times10^{9}.)

Physics of the portrait

Superluminal motion is one of the rare cases where theory arrived before data. Five years before any interferometer could catch a jet moving, a graduate student wrote down the projection effect, the formula, and the prediction that expanding radio sources should appear to fly apart faster than light.

Origins · 1 minMartin Rees and how this was discovered
MR
Martin Rees
1942– · relativistic astrophysics, black holes, and cosmology

Martin Rees was a 24-year-old research student at Cambridge when, in 1966, he published a short Nature paper titled "Appearance of Relativistically Expanding Radio Sources." The radio astronomy of the day had found quasars varying so fast that, taken naïvely, the emitting region seemed to be expanding at many times the speed of light. Rees showed no such thing was needed: a source moving relativistically toward the observer would look superluminal purely through light-travel-time compression, and he wrote down the geometry that gives βapp=βsin⁡θ/(1−βcos⁡θ)\beta_{\rm app}=\beta\sin\theta/(1-\beta\cos\theta) — a prediction with almost no supporting evidence, since the resolving power to test it did not yet exist. When very-long-baseline interferometry matured around 1971, the archetypal quasars 3C 273 and 3C 279 were resolved into components separating superluminally, precisely as his picture required. Rees went on to co-author, with Blandford and McKee, the relativistic-beaming framework that ties superluminal motion to Doppler boosting and jet one-sidedness, and became one of the architects of modern high-energy astrophysics — arguing early that quasars are powered by accretion onto supermassive black holes, and later ranging across gamma-ray bursts and cosmology as Astronomer Royal and President of the Royal Society. The superluminal prediction remains a textbook case of theory-led astrophysics: a back-of-envelope relativity argument that dissolved an apparent paradox before the instruments could even pose it.

History
1966
— Rees predicts that relativistically expanding radio sources will appear to move superluminally
1971
— early VLBI resolves 3C 279 and 3C 273 into components separating faster than light
1977
— Blandford, McKee and Rees set out the relativistic-beaming interpretation of the phenomenon
1979
— Blandford and Königl publish the relativistic jet model unifying superluminal motion with one-sidedness
1994
— the Very Long Baseline Array begins routine multi-epoch monitoring of AGN jet kinematics
2018
— the MOJAVE survey catalogues apparent speeds and accelerations for hundreds of AGN jets

Going further

Texts. Carroll & Ostlie, An Introduction to Modern Astrophysics, 2nd ed., §28 (Active Galaxies) — the standard treatment of superluminal motion and beaming · Rybicki & Lightman, Radiative Processes in Astrophysics, §4.8 (relativistic beaming) · Dodelson, Modern Cosmology (angular-diameter distance, equation 4) · Mo, van den Bosch & White, Galaxy Formation and Evolution, Ch. 14 (AGN and unification) · Binney & Tremaine, Galactic Dynamics · Ryden, Introduction to Cosmology, 2nd ed.

Primary sources. Rees 1966, Nature 211, 468 (the prediction) · Blandford, McKee & Rees 1977, Nature 267, 211 (beaming) · Blandford & Königl 1979, ApJ 232, 34 (jet model) · Lister et al. (MOJAVE) 2019, ApJ — VLBI apparent-speed statistics. The kinematics rest on pure special-relativistic light-travel-time geometry. The key relations are: (1) βapp=βsin⁡θ/(1−βcos⁡θ)\beta_{\rm app}=\beta\sin\theta/(1-\beta\cos\theta); (2) the maximisation giving cos⁡θ=β\cos\theta=\beta and βapp,max=βγ\beta_{\rm app,max}=\beta\gamma; (3) the bound γ≥1+βapp2\gamma\ge\sqrt{1+\beta_{\rm app}^2}, derived via β2γ2=γ2−1\beta^2\gamma^2=\gamma^2-1; (4) the identity βapp2+δ2=2γδ−1\beta_{\rm app}^2+\delta^2=2\gamma\delta-1 and its inversion γ=(βapp2+δ2+1)/(2δ)\gamma=(\beta_{\rm app}^2+\delta^2+1)/(2\delta); (5) the corrected left identity βapp=γδβsin⁡θ\beta_{\rm app}=\gamma\delta\beta\sin\theta, since γδ=1/(1−βcos⁡θ)\gamma\delta=1/(1-\beta\cos\theta). Numerical anchors: βapp=6.0\beta_{\rm app}=6.0 from μ=0.2\mu=0.2 mas yr−1^{-1}, z=0.5z=0.5, DA=1267D_A=1267 Mpc; the joint solve βapp=15, δ=12⇒γ=15.4, θ=4.7∘\beta_{\rm app}=15,\ \delta=12\Rightarrow\gamma=15.4,\ \theta=4.7^\circ; and the counterjet brightness ratio ∼7×106\sim7\times10^{6}.

Prerequisites. Special Relativity and the Lorentz Factor · Relativistic Jets · AGN: the Unified Model · Cosmological Distances and Redshift. Leads to. Blazars and Relativistic Beaming · Black Hole Spin · Quasars.

Ask the tutor: Give me a measured apparent speed (and a Doppler factor if you have one) and I will walk you through extracting the jet's Lorentz factor, true speed, and viewing angle step by step.

Does superluminal motion mean something actually moves faster than light?

No. The physical bulk speed of the jet plasma is always below cc — in the worked examples it is 99.8%99.8\% of cc, never more. What exceeds cc is the apparent transverse speed βapp\beta_{\rm app}: the sky-projected distance a knot seems to cover divided by the light-distorted time over which we receive its light. Because an approaching blob nearly keeps pace with its own photons, they arrive bunched together, and dividing a real transverse motion by an artificially short arrival interval gives a ratio that can top cc. No matter, energy, or signal travels superluminally, so relativity is untouched — the effect is a projection, like a searchlight spot sweeping a distant wall.

Why does the apparent speed peak at an intermediate angle instead of head-on?

Two effects fight each other. The transverse motion you can see scales as sin⁡θ\sin\theta, which vanishes head-on — a jet aimed straight at you shows no sideways drift. The arrival-time compression scales as 1/(1−βcos⁡θ)1/(1-\beta\cos\theta), strongest head-on and weakening as the angle grows. Their product βapp=βsin⁡θ/(1−βcos⁡θ)\beta_{\rm app}=\beta\sin\theta/(1-\beta\cos\theta) must therefore peak in between, and calculus pins the sweet spot at cos⁡θ=β\cos\theta=\beta, i.e. sin⁡θ=1/γ\sin\theta=1/\gamma. For a fast jet that angle is only a few degrees, and the peak value there is βγ≈γ\beta\gamma\approx\gamma. The largest apparent speeds come not from jets aimed exactly at us, but from those aimed just slightly off.

How does one measured β_app become a limit on the Lorentz factor?

The apparent speed can never exceed the maximum available to the true Lorentz factor, βapp,max=βγ\beta_{\rm app,max}=\beta\gamma, so a measurement forces βapp≤βγ\beta_{\rm app}\le\beta\gamma. Squaring and using β2γ2=γ2−1\beta^2\gamma^2=\gamma^2-1 eliminates the speed, leaving βapp2≤γ2−1\beta_{\rm app}^2\le\gamma^2-1, i.e. γ≥1+βapp2\gamma\ge\sqrt{1+\beta_{\rm app}^2}. No viewing angle is needed: whatever the geometry, the Lorentz factor must be at least this large. A knot at βapp=10\beta_{\rm app}=10 requires γ≥10.05\gamma\ge10.05, and no slow flow can counterfeit a fast apparent speed — which is why one clean proper-motion measurement certifies a jet as relativistic.

If a single β_app only gives a lower bound, how do we get the real γ and θ?

By measuring a second, independent quantity — usually the Doppler factor δ\delta, from the source's variability timescale and brightness temperature or from synchrotron self-Compton arguments. Since both βapp\beta_{\rm app} and δ\delta depend on the same (β,θ)(\beta,\theta), they satisfy exact identities; the cleanest, βapp2+δ2=2γδ−1\beta_{\rm app}^2+\delta^2=2\gamma\delta-1, inverts to γ=(βapp2+δ2+1)/(2δ)\gamma=(\beta_{\rm app}^2+\delta^2+1)/(2\delta), and then β\beta and θ\theta follow. Worked example 2 runs this through for βapp=15\beta_{\rm app}=15, δ=12\delta=12, giving γ≈15.4\gamma\approx15.4, θ≈4.7∘\theta\approx4.7^\circ. Two observables break the degeneracy and fix the full kinematics.

Why do these jets so often show only one side?

The same geometry that fakes superluminal motion boosts the approaching jet and suppresses the receding one. Flux scales as S∝δpS\propto\delta^{p} with p≈2p\approx2–44, and the two sides have Doppler factors set by 1∓βcos⁡θ1\mp\beta\cos\theta, so their brightness ratio is [(1+βcos⁡θ)/(1−βcos⁡θ)]p[(1+\beta\cos\theta)/(1-\beta\cos\theta)]^{p} — reaching millions to billions for β≈1\beta\approx1 and small θ\theta (worked example 3). The counterjet is not absent; it is de-boosted below any realistic detection threshold. This one-sidedness is itself evidence for relativistic beaming.

Why is a cosmological model needed to quote an apparent speed?

Because an interferometer measures an angular rate — a proper motion μ\mu in mas yr−1^{-1} — not a velocity. Converting it to a transverse velocity needs the physical distance, which comes from the redshift through a cosmological model: βapp=μDA(1+z)/c\beta_{\rm app}=\mu D_A(1+z)/c. The angular-diameter distance DAD_A depends on H0H_0, Ωm\Omega_m, ΩΛ\Omega_\Lambda, and the (1+z)(1+z) accounts for time dilation between emitted and observed frames. Two knots at the same true speed but different redshifts show different proper motions, and the wrong cosmology gives the wrong βapp\beta_{\rm app}. The bound γ≥1+βapp2\gamma\ge\sqrt{1+\beta_{\rm app}^2} inherits this dependence only through βapp\beta_{\rm app}.

Before this

Where this leads