Null Geodesics and Light Cones

A Lorentzian metric does one thing before it does anything else: at every event it sorts the tangent vectors into three classes by the sign of gμνVμVνg_{\mu\nu}V^\mu V^\nu — timelike (negative), spacelike (positive), and null (zero). Massive bodies thread the timelike interior; light and every other massless field ride the null boundary between what can be reached and what cannot. The curves that trace that boundary as economically as the geometry allows — everywhere null, and parallel-transporting their own tangent — are the null geodesics, and the surface they sweep out from an event is its light cone. Everything causal in general relativity is written in this one object. Two conditions define a light ray at once. It is null, kμkμ=0k_\mu k^\mu=0 for the tangent kμ=dxμ/dλk^\mu=\mathrm{d}x^\mu/\mathrm{d}\lambda, so no proper time elapses along it and it carries zero interval between any two of its points. And it is a geodesic, kν∇νkμ=0k^\nu\nabla_\nu k^\mu=0, its tangent carried forward unbent. Because there is no arc length to lean on, the only natural clock a photon has is an affine parameter, fixed by the geometry up to λ→aλ+b\lambda\to a\lambda+b; the physical observables are never λ\lambda itself but ratios built from it — frequencies and redshifts. Three deep facts organize the rest. Null geodesics generate the boundaries ∂J±\partial J^{\pm} of the causal future and past, so they encode which events can signal which. They are conformally invariant, unchanged as point sets under gμν→Ω2gμνg_{\mu\nu}\to\Omega^2 g_{\mu\nu} — the reason light cones survive being drawn on a finite Penrose diagram. And their neighbouring rays focus under the Raychaudhuri equation, the engine behind trapped surfaces and the singularity theorems.

Key takeaways
  • A null geodesic is a light ray: a curve whose tangent kμk^\mu is everywhere null, gμνkμkν=0g_{\mu\nu}k^\mu k^\nu=0, and parallel-transported along itself, kν∇νkμ=0k^\nu\nabla_\nu k^\mu=0.
  • The null geodesics through an event sweep out its light cone; extended to the future they generate ∂J+(p)\partial J^+(p), the boundary of the causal future, so they carry the entire causal structure of a spacetime.
  • Proper time is useless as a label (dτ=0\mathrm{d}\tau=0 along the ray); an affine parameter λ\lambda is supplied intrinsically by parallel transport, fixed only up to λ→aλ+b\lambda\to a\lambda+b.
  • Null geodesics are conformally invariant — under gμν→Ω2gμνg_{\mu\nu}\to\Omega^2 g_{\mu\nu} they remain null geodesics (affine parameter rescaled) — which is exactly why light cones stay at 45° on a conformal diagram.
  • Observed frequency is the projection ω=−kμuμ\omega=-k_\mu u^\mu onto an observer's 4-velocity; only ratios of these — redshifts — are affine-independent physics.
  • Congruences of null geodesics obey the Raychaudhuri equation, and gravity always focuses light for matter obeying the null energy condition, driving beams to caustics and seeding the singularity theorems.
  • Key idea: A light ray is the unique curve that is both null and self-parallel; affineness fixes its parametrization, conformal invariance fixes that it depends only on the causal structure, and the Raychaudhuri equation governs how neighbouring rays focus.
The core insight

A light ray is the unique curve that is both null and self-parallel; affineness fixes its parametrization, conformal invariance fixes that the curve depends only on the causal structure, and the Raychaudhuri equation governs how neighbouring rays focus into caustics.

Where this lives today

Null geodesics are the working instrument of strong-gravity astronomy. Every image of a black hole — the Event Horizon Telescope's reconstructions of M87\ and Sagittarius A\, and every synthetic frame of Interstellar's Gargantua — is made by integrating null geodesics backward from a virtual camera, with the photon-sphere result of Worked Example 1 setting the shadow size at bcrit=33 Mb_{\rm crit}=3\sqrt3\,M. Gravitational lensing, from cluster arcs to weak-lensing shear maps, is null-geodesic propagation through inhomogeneous spacetime. In gravitational-wave astronomy, the multi-messenger event GW170817 sent a neutron-star merger's gravitational and electromagnetic signals along the same null cone; their arrival within seconds after 130 million years bounded the speed of gravity to that of light at the 10−1510^{-15} level — a direct test that gravitational waves ride null geodesics. At the frontier, the focusing of null congruences underwrites the black-hole area theorem, the conformal boundary I\mathscr{I} hosts the BMS symmetries and the memory effect, and cosmic censorship — whether trapped-surface focusing always hides its singularity — remains one of the deepest open problems in classical relativity.

The physical picture

Set off a flash of light at an event pp. The spherical wavefront expanding at cc traces, in spacetime, a cone — the future light cone. Each individual photon worldline is one straight ruling of that cone, and each ruling is a null geodesic. The cone is not an abstraction laid over the physics; it is literally the union of the light rays, and its interior is the set of events a massive body launched from pp could still reach. The first figure makes this anatomy concrete: null generators on the surface, the timelike causal future inside, the spacelike elsewhere outside, and the mirror-image past cone below.

Figure
A flash at event p and its light cone. The future cone is the union of null-geodesic generators, each a photon worldline with ds² = 0. Its timelike interior is the causal future J⁺(p) reachable by massive bodies at event q; the spacelike exterior is the elsewhere that no signal from p can reach, such as event s. In Minkowski coordinates every generator runs at the speed of light, a 45° cone
The light cone is the union of null-geodesic generators — ds² = 0 on every ray t p q timelike from p s spacelike: no contact future light cone ∂J⁺(p) past light cone ∂J⁻(p) null geodesic timelike (inside): g(V,V) negative null (on cone): g(V,V) = 0 spacelike (outside): g(V,V) positive

Because light has no rest frame and no proper time elapses along its path, "how far along the ray" cannot be a length — it is an affine parameter, fixed only up to scale and origin. The observable is never the parameter but a ratio: project the tangent kμk^\mu onto an observer's 4-velocity and ω=−kμuμ\omega=-k_\mu u^\mu is that observer's frequency; compare two observers and the ratio is a redshift, blind to how λ\lambda was scaled. Conformal invariance is the deeper secret: stretching the metric by any positive factor cannot turn a null vector timelike or spacelike, so cones — and hence causality — survive rescaling, letting an infinite spacetime be bent onto a finite page. The last idea is focusing. Neighbouring rays in a beam converge or diverge at a rate called the expansion θ\theta, and gravity, entering the Raychaudhuri equation through the Ricci tensor, always bends a bundle of light toward itself when matter obeys the null energy condition. Push a converging beam and it must reach a caustic where rays cross — the seed of trapped surfaces. The quantitative face of this lives in a black hole, where light orbits at a fixed radius; the second figure plots the effective potential that pins that orbit.

Figure
The Schwarzschild photon effective potential V(r) = (1/r²)(1 − 2M/r) for a light ray in geometric units M = 1, L = 1. It has a single maximum at the photon sphere r = 3M where V = 1/27 ≈ 0.0370, an unstable circular light orbit. A ray with E²/L² above the peak is captured; the critical value E²/L² = 1/27 gives the capture impact parameter b = 3√3 M ≈ 5.196 M, which sets the edge of a black-hole shadow
Photon effective potential: peak at r = 3M fixes b = 3√3 M ≈ 5.196 M 0 0.01 0.02 0.03 0.04 2 5 10 15 20 radius r / M effective potential V(r) photon sphere r = 3M, V = 1/27 unstable circular light orbit E²/L² above 1/27 ⇒ captured (b under 3√3 M) V(r) = (1/r²)(1 − 2M/r)

The core equations

A light ray obeys two statements simultaneously, and neither implies the other. The first is that the tangent has zero length everywhere — the ray lies on the local light cone at each point it passes.

gμν kμkν=0,kμ=dxμdλ(1)g_{\mu\nu}\,k^\mu k^\nu = 0,\qquad k^\mu=\frac{\mathrm{d}x^\mu}{\mathrm{d}\lambda} \qquad (1)

where: gμνg_{\mu\nu} — the spacetime metric (signature −+++-+++); kμk^\mu — the tangent to the ray, the wave 4-vector in optics; λ\lambda — an affine parameter along the ray. Read it as: the ray never leaves the null cone, so no observer ever measures it moving slower or faster than light locally. The second statement is that the tangent is carried forward by parallel transport with nothing left over — the geodesic condition, in the affine form that has zero right-hand side. Expanding the covariant derivative gives a second-order ODE for the coordinates.

kν∇νkμ=d2xμdλ2+Γμαβdxαdλdxβdλ=0(2)k^\nu\nabla_\nu k^\mu = \frac{\mathrm{d}^2 x^\mu}{\mathrm{d}\lambda^2}+\Gamma^\mu{}_{\alpha\beta}\frac{\mathrm{d}x^\alpha}{\mathrm{d}\lambda}\frac{\mathrm{d}x^\beta}{\mathrm{d}\lambda}=0 \qquad (2)

where: ∇ν\nabla_\nu — the Levi-Civita covariant derivative; Γμαβ\Gamma^\mu{}_{\alpha\beta} — the Christoffel symbols of gg. The clean zero on the right holds only for an affine λ\lambda; any other label puts a term ∝kμ\propto k^\mu there, which bends nothing but spoils the bookkeeping. Metric compatibility makes the null value kμkμ=0k_\mu k^\mu=0 propagate: it is set once and conserved for all λ\lambda. Frequency is not intrinsic to the ray — it belongs to whoever measures it. Contracting the tangent with an observer's 4-velocity gives that observer's angular frequency; the ratio at two ends of the ray is the physical, affine-independent redshift.

ω=− kμuμ,ωrecωemit=(kμuμ)rec(kμuμ)emit(3)\omega = -\,k_\mu u^\mu,\qquad \frac{\omega_{\rm rec}}{\omega_{\rm emit}}=\frac{(k_\mu u^\mu)_{\rm rec}}{(k_\mu u^\mu)_{\rm emit}} \qquad (3)

where: uμu^\mu — the observer's 4-velocity, uμuμ=−1u_\mu u^\mu=-1 (c=1c=1); ω\omega — the frequency that observer sees (its absolute value depends on the arbitrary scaling of λ\lambda, but the ratio does not). Read it as: the affine ambiguity cancels in every observable, so only shifts, never bare frequencies, are physics. The single most powerful structural property is that the whole construction depends only on the metric up to a positive scale. Under a conformal rescaling the null cones are untouched, and null geodesics map to null geodesics.

g~μν=Ω2gμν ⇒ k~ν∇~νkμ=2(kν∂νln⁡Ω) kμ ∝ kμ(4)\tilde g_{\mu\nu}=\Omega^2 g_{\mu\nu}\ \Rightarrow\ \tilde k^\nu\tilde\nabla_\nu k^\mu = 2\big(k^\nu\partial_\nu\ln\Omega\big)\,k^\mu\ \propto\ k^\mu \qquad (4)

where: Ω>0\Omega>0 — the conformal factor; the right-hand side is purely tangential, so the curve is still a geodesic, merely re-labelled. Read it as: causality is a conformal invariant — which events can signal which is decided before any scale is chosen, and this is what makes Penrose–Carter diagrams faithful. For a bundle of neighbouring rays, the expansion θ=∇μkμ\theta=\nabla_\mu k^\mu measures the fractional growth rate of the beam's cross-section, and its evolution is governed by shear, twist, and curvature — the null Raychaudhuri equation.

dθdλ=−12θ2−σμνσμν+ωμνωμν−Rμνkμkν(5)\frac{\mathrm{d}\theta}{\mathrm{d}\lambda}=-\tfrac{1}{2}\theta^2-\sigma_{\mu\nu}\sigma^{\mu\nu}+\omega_{\mu\nu}\omega^{\mu\nu}-R_{\mu\nu}k^\mu k^\nu \qquad (5)

where: σμν\sigma_{\mu\nu} — shear (trace-free distortion); ωμν\omega_{\mu\nu} — twist (rotation) of the congruence; RμνR_{\mu\nu} — the Ricci tensor, with Rμνkμkν≥0R_{\mu\nu}k^\mu k^\nu\ge0 the null energy condition via Einstein's equations. Read it as: every term except the twist drives θ\theta downward, so gravity focuses light — the mathematical root of the singularity theorems. In the high-frequency limit a wave ψ∼A eiS\psi\sim A\,e^{iS} has its phase surfaces governed by the eikonal equation, and the null geodesics are exactly the rays orthogonal to those surfaces — the bridge from wave optics to light rays.

gμν∂μS ∂νS=0,kμ=∂μS(6)g^{\mu\nu}\partial_\mu S\,\partial_\nu S=0,\qquad k_\mu=\partial_\mu S \qquad (6)

where: SS — the rapidly-varying phase; kμ=∂μSk_\mu=\partial_\mu S — the wave covector, automatically null and (being a gradient) twist-free. Read it as: geometric optics is the leading term of wave optics, and its "rays" are precisely null geodesics of the metric.

Derivation ANull geodesic and null condition from the energy action

The goal is to derive the affine geodesic equation and the null condition together from a single variational principle, and to see why the usual arc-length action is helpless for light while a quadratic "energy" action succeeds. The square-root Lagrangian ∫−gμνx˙μx˙ν dλ\int\sqrt{-g_{\mu\nu}\dot x^\mu\dot x^\nu}\,\mathrm{d}\lambda integrates to zero along a null curve and its variation is ill-defined there; the quadratic action sidesteps this, gives the same geodesics for every signature, and hands us an affine parameter for free — the one in which the tangent norm is conserved, so the null value can be imposed once and holds everywhere.

Step 1. Adopt the energy action S=12∫gμνx˙μx˙ν dλS=\tfrac12\int g_{\mu\nu}\dot x^\mu\dot x^\nu\,\mathrm{d}\lambda with x˙μ≡dxμ/dλ\dot x^\mu\equiv\mathrm{d}x^\mu/\mathrm{d}\lambda, and vary xμx^\mu to get the Euler–Lagrange equation ddλ(gμνx˙ν)−12∂μgαβx˙αx˙β=0\frac{\mathrm{d}}{\mathrm{d}\lambda}(g_{\mu\nu}\dot x^\nu)-\tfrac12\partial_\mu g_{\alpha\beta}\dot x^\alpha\dot x^\beta=0. Why: unlike the arc length, this Lagrangian is smooth and nonzero on null curves, so its variation is legitimate for light — the price is that it is no longer reparametrization-invariant, which is exactly what pins down an affine parameter.

Step 2. Expand the total derivative and symmetrize: gμνx¨ν+12(∂αgμβ+∂βgμα−∂μgαβ)x˙αx˙β=0g_{\mu\nu}\ddot x^\nu+\tfrac12(\partial_\alpha g_{\mu\beta}+\partial_\beta g_{\mu\alpha}-\partial_\mu g_{\alpha\beta})\dot x^\alpha\dot x^\beta=0, where the bracket is gμλΓλαβg_{\mu\lambda}\Gamma^\lambda{}_{\alpha\beta}. Why: the particular antisymmetric-looking combination of metric derivatives is engineered to reassemble the Christoffel symbol, so the field equation is secretly the covariant one.

Step 3. Contract with the inverse metric to reach x¨μ+Γμαβx˙αx˙β=0\ddot x^\mu+\Gamma^\mu{}_{\alpha\beta}\dot x^\alpha\dot x^\beta=0, which is Eq. (2). Why: this is the affine geodesic equation with no right-hand side — it holds identically for timelike, spacelike, and null geodesics, so light is not a special case but the null member of one family.

Step 4. Note that the Lagrangian has no explicit λ\lambda-dependence, so its value gμνx˙μx˙νg_{\mu\nu}\dot x^\mu\dot x^\nu is conserved along the solution; choose that constant to be 00. Why: the null condition gμνkμkν=0g_{\mu\nu}k^\mu k^\nu=0 is therefore a first integral, imposed as an initial condition and preserved automatically — this is the precise sense in which "null" and "geodesic" are compatible but logically independent statements. *In flat space Γ=0\Gamma=0 gives x¨μ=0\ddot x^\mu=0, so the solutions are straight lines xμ=x0μ+λkμx^\mu=x_0^\mu+\lambda k^\mu with ημνkμkν=0\eta_{\mu\nu}k^\mu k^\nu=0; the conserved norm is the affine hallmark.*

Go deeper · Derivation B · 1 minNull geodesics are conformally invariant
Derivation BNull geodesics are conformally invariant

The goal is to prove the claim that makes conformal diagrams honest: that multiplying the metric by any positive function Ω2\Omega^2 leaves the set of null geodesics unchanged, merely rescaling their affine parameter — and to isolate the single term whose vanishing (thanks to nullity) makes the result hold for light but fail for massive particles. We track how the Christoffel symbols transform under g~μν=Ω2gμν\tilde g_{\mu\nu}=\Omega^2 g_{\mu\nu}, feed the change into the geodesic operator, and watch the metric-dependent piece disappear precisely because kμkμ=0k_\mu k^\mu=0.

Step 1. Under g~μν=Ω2gμν\tilde g_{\mu\nu}=\Omega^2 g_{\mu\nu} the connection shifts by Γ~μαβ=Γμαβ+δαμ∂βln⁡Ω+δβμ∂αln⁡Ω−gαβ∂μln⁡Ω\tilde\Gamma^\mu{}_{\alpha\beta}=\Gamma^\mu{}_{\alpha\beta}+\delta^\mu_\alpha\partial_\beta\ln\Omega+\delta^\mu_\beta\partial_\alpha\ln\Omega-g_{\alpha\beta}\partial^\mu\ln\Omega. Why: a conformal factor is not a mere coordinate change, so the connection genuinely changes; the three extra terms are the entire content of that change, and only the last one carries the metric gαβg_{\alpha\beta}.

Step 2. Insert this into k~ν∇~νkμ\tilde k^\nu\tilde\nabla_\nu k^\mu with the same tangent kμk^\mu, giving kν∇νkμ+2(kν∂νln⁡Ω)kμ−(gαβkαkβ)∂μln⁡Ωk^\nu\nabla_\nu k^\mu+2(k^\nu\partial_\nu\ln\Omega)k^\mu-(g_{\alpha\beta}k^\alpha k^\beta)\partial^\mu\ln\Omega. Why: the two δ\delta-terms collapse into a single tangential piece ∝kμ\propto k^\mu, while the gαβg_{\alpha\beta}-term is the only one pointing off the ray — it is the potential troublemaker.

Step 3. Impose nullity gαβkαkβ=0g_{\alpha\beta}k^\alpha k^\beta=0 so the off-ray term vanishes, leaving k~ν∇~νkμ=2(kν∂νln⁡Ω) kμ∝kμ\tilde k^\nu\tilde\nabla_\nu k^\mu=2(k^\nu\partial_\nu\ln\Omega)\,k^\mu\propto k^\mu. Why: a right-hand side proportional to kμk^\mu is still a geodesic — just non-affinely parametrized — so the curve survives intact; for a timelike geodesic gαβkαkβ=−1≠0g_{\alpha\beta}k^\alpha k^\beta=-1\ne0 and the extra term bends the path, which is why only light cones are conformally invariant.

Step 4. Restore affineness by rescaling λ→λ~\lambda\to\tilde\lambda with dλ~/dλ=Ω2\mathrm{d}\tilde\lambda/\mathrm{d}\lambda=\Omega^2 (up to constants). Why: the residual ∝kμ\propto k^\mu source is removable by a one-line reparametrization, so in the rescaled metric the very same points form an affinely-parametrized null geodesic — the light ray is a conformal invariant as a curve. *As a check, for Ω=\Omega= const the result reduces to the identity; the surviving proportionality-to-kμk^\mu is exactly the removable non-affineness.*

Go deeper · Derivation C · 1 minRaychaudhuri focusing and the inevitability of caustics
Derivation CRaychaudhuri focusing and the inevitability of caustics

The goal is to show that a bundle of light rays which is momentarily converging must, under gravity plus the null energy condition, focus to a caustic within a bounded affine distance — the local mechanism behind trapped surfaces and the Penrose singularity theorem. We start from the null Raychaudhuri equation for the expansion, specialize to the hypersurface-orthogonal (twist-free) congruences that generate null surfaces, drop everything that only helps focusing, and integrate the resulting inequality to bound the distance to θ→−∞\theta\to-\infty.

Step 1. Take Eq. (5) and set the twist to zero, ωμν=0\omega_{\mu\nu}=0, which holds for the generators of any null hypersurface such as ∂J+(p)\partial J^+(p) or an event horizon. Why: a gradient wave covector kμ=∂μSk_\mu=\partial_\mu S is automatically twist-free, so the surface-forming congruences that carry all the causal information are exactly the ones this simplification applies to.

Step 2. Discard the shear and curvature terms as manifestly non-positive: with the null energy condition Rμνkμkν≥0R_{\mu\nu}k^\mu k^\nu\ge0 and σμνσμν≥0\sigma_{\mu\nu}\sigma^{\mu\nu}\ge0, the equation becomes the inequality dθ/dλ≤−12θ2\mathrm{d}\theta/\mathrm{d}\lambda\le-\tfrac12\theta^2. Why: both dropped terms subtract from dθ/dλ\mathrm{d}\theta/\mathrm{d}\lambda, so ignoring them only weakens the bound — whatever focusing the inequality forces, the full equation forces at least as strongly.

Step 3. Rewrite as d(θ−1)/dλ≥12\mathrm{d}(\theta^{-1})/\mathrm{d}\lambda\ge\tfrac12 and integrate to θ−1(λ)≥θ0−1+12λ\theta^{-1}(\lambda)\ge\theta_0^{-1}+\tfrac12\lambda. Why: dividing by −θ2-\theta^2 turns the nonlinear inequality into a linear one for the reciprocal, which integrates immediately and makes the blow-up transparent.

Step 4. Suppose the beam starts converging, θ0<0\theta_0<0. Then θ−1\theta^{-1} climbs from a negative value through zero, so θ→−∞\theta\to-\infty at an affine distance no greater than Δλ≤2/∣θ0∣\Delta\lambda\le 2/|\theta_0|. Why: θ→−∞\theta\to-\infty is a caustic — neighbouring rays cross and the cross-sectional area collapses — and this is precisely the ingredient Penrose combined with a non-compact Cauchy surface to prove geodesic incompleteness inside a trapped surface. *As a concrete illustration, a Schwarzschild ingoing null shell with initial expansion θ0<0\theta_0<0 caustics within affine distance Δλ≤2/∣θ0∣\Delta\lambda\le2/|\theta_0|, in agreement with the reciprocal bound derived above.*

Limiting cases

Special regimes make the behaviour of null geodesics and their cones concrete:

CaseResultMeaning
Flat spacetime (Minkowski)xμ(λ)=x0μ+λ kμ, ημνkμkν=0x^\mu(\lambda)=x_0^\mu+\lambda\,k^\mu,\ \eta_{\mu\nu}k^\mu k^\nu=0Null geodesics are straight lines along null vectors; light cones are exact 45° cones in every inertial frame.
Conformal rescalinggμν→Ω2gμνg_{\mu\nu}\to\Omega^2 g_{\mu\nu}: same null geodesics, λ~\tilde\lambda rescaledThe rays are unchanged as curves; the new affine parameter obeys dλ~/dλ=Ω2\mathrm{d}\tilde\lambda/\mathrm{d}\lambda=\Omega^2. Causality is conformally invariant.
Geometric-optics (eikonal) limitkμ=∂μS, gμν∂μS ∂νS=0k_\mu=\partial_\mu S,\ g^{\mu\nu}\partial_\mu S\,\partial_\nu S=0Surfaces of constant phase of a high-frequency wave obey the eikonal equation; their normals are null geodesics.
Schwarzschild photon spherer=3M, bcrit=33 M≈5.196 Mr=3M,\ b_{\rm crit}=3\sqrt3\,M\approx5.196\,MUnstable circular null geodesics orbit at r=3Mr=3M and set the black-hole shadow radius 33 M3\sqrt3\,M.
FLRW cosmology1+z=a(trec)a(temit)1+z=\dfrac{a(t_{\rm rec})}{a(t_{\rm emit})}Light on a radial null geodesic in an expanding universe is redshifted exactly by the ratio of scale factors.

Worked examples

1 · The photon sphere and the shadow of a black hole (fully worked)

A light ray moves in the equatorial plane of a Schwarzschild black hole of mass MM (geometric units G=c=1G=c=1). Find the radius of the unstable circular photon orbit, the critical impact parameter that separates capture from escape, and the angular radius of the resulting shadow for Sagittarius A\ and M87\.

Step 1. The Killing vectors ∂t\partial_t and ∂ϕ\partial_\phi give two conserved charges along the affinely-parametrized ray, E=(1−2M/r)t˙E=(1-2M/r)\dot t and L=r2ϕ˙L=r^2\dot\phi (dots are d/dλ\mathrm{d}/\mathrm{d}\lambda). These are constant only because λ\lambda is affine — a non-affine label would let each drift. The null condition gμνx˙μx˙ν=0g_{\mu\nu}\dot x^\mu\dot x^\nu=0 then collapses to a one-dimensional radial problem, r˙2=E2−V(r)\dot r^2=E^2-V(r) with effective potential V(r)=L2r2(1−2Mr)V(r)=\dfrac{L^2}{r^2}\Big(1-\dfrac{2M}{r}\Big).

Step 2. A circular orbit needs r˙=0\dot r=0 and r¨=0\ddot r=0, i.e. V′(r)=0V'(r)=0. Differentiating, V′(r)=L2(−2r−3+6Mr−4)=0⇒−2r+6M=0V'(r)=L^2\big(-2r^{-3}+6Mr^{-4}\big)=0\Rightarrow -2r+6M=0, so rph=3Mr_{\rm ph}=3M — the photon sphere. It is a maximum of VV, hence unstable: a ray perturbed inward plunges, one perturbed outward escapes.

Step 3. At r=3Mr=3M the potential is V(3M)=L29M2(1−23)=L227M2V(3M)=\dfrac{L^2}{9M^2}\big(1-\tfrac23\big)=\dfrac{L^2}{27M^2}. A circular orbit has E2=VE^2=V, so the critical impact parameter is bcrit=LE=27 M=33 M≈5.196 Mb_{\rm crit}=\dfrac{L}{E}=\sqrt{27}\,M=\mathbf{3\sqrt3\,M\approx5.196\,M}. Rays with b<bcritb<b_{\rm crit} are captured; b>bcritb>b_{\rm crit} escape; b=bcritb=b_{\rm crit} spirals asymptotically onto the photon sphere. This bcritb_{\rm crit}, not the horizon at 2M2M, is the apparent edge of the dark shadow.

Step 4. Project onto the sky. For Sagittarius A\, M≈4.15×106 M⊙M\approx4.15\times10^6\,M_\odot gives GM/c2≈6.13×106 kmGM/c^2\approx6.13\times10^6\ \mathrm{km}, so bcrit≈3.18×107 kmb_{\rm crit}\approx3.18\times10^7\ \mathrm{km}; at D≈8.15 kpcD\approx8.15\ \mathrm{kpc} the angular radius is bcrit/D≈26 μasb_{\rm crit}/D\approx\mathbf{26\ \mu as} (diameter ≈52 μas\approx52\ \mu\mathrm{as}). For M87\, M≈6.5×109 M⊙M\approx6.5\times10^9\,M_\odot at D≈16.8 MpcD\approx16.8\ \mathrm{Mpc} gives an angular radius ≈20 μas\approx\mathbf{20\ \mu as} (diameter ≈40 μas\approx40\ \mu\mathrm{as}). Both match the bright rings imaged by the Event Horizon Telescope — the shadow is drawn by null geodesics, not by luminous matter.

2 · Gravitational redshift from a compact star (numbers)

Light of frequency ωemit\omega_{\rm emit} leaves the surface of a static star at areal radius RR and is received by a distant static observer. Using the conserved quantity along the null geodesic, find the redshift for the Sun and for a 1.4 M⊙1.4\,M_\odot neutron star of radius 10 km10\ \mathrm{km}.

Step 1. In a static spacetime the timelike Killing vector ξμ=∂t\xi^\mu=\partial_t makes E=−kμξμE=-k_\mu\xi^\mu constant along the ray. A static observer's 4-velocity is the normalized Killing vector, uμ=ξμ/−gttu^\mu=\xi^\mu/\sqrt{-g_{tt}}, so the measured frequency is ω=−kμuμ=E/−gtt\omega=-k_\mu u^\mu=E/\sqrt{-g_{tt}} with −gtt=1−2M/r-g_{tt}=1-2M/r in Schwarzschild.

Step 2. Taking the ratio between emission at RR and reception at r→∞r\to\infty (where −gtt→1-g_{tt}\to1), the affine scale and EE cancel: ωrecωemit=1−2M/R1\dfrac{\omega_{\rm rec}}{\omega_{\rm emit}}=\sqrt{\dfrac{1-2M/R}{1}}, hence 1+z=ωemitωrec=11−2GM/Rc21+z=\dfrac{\omega_{\rm emit}}{\omega_{\rm rec}}=\dfrac{1}{\sqrt{1-2GM/Rc^2}}. The light climbs out of the potential well and loses frequency.

Step 3. For the Sun, 2GM/Rc2=Rs/R=2.95/6.96×105=4.24×10−62GM/Rc^2=R_s/R=2.95/6.96\times10^{5}=4.24\times10^{-6}, so z=(1−4.24×10−6)−1/2−1≈2.1×10−6z=(1-4.24\times10^{-6})^{-1/2}-1\approx\mathbf{2.1\times10^{-6}} — the classic solar redshift, equivalent to a 636 m/s636\ \mathrm{m/s} Doppler shift and confirmed spectroscopically.

Step 4. For the neutron star, 2GM/Rc2=(2×1.4×1.477 km)/10 km=0.4132GM/Rc^2=(2\times1.4\times1.477\ \mathrm{km})/10\ \mathrm{km}=0.413, giving z=(1−0.413)−1/2−1≈0.31z=(1-0.413)^{-1/2}-1\approx\mathbf{0.31}. A spectral line at 500 nm arrives at 654 nm — a 31% shift, and the reason neutron-star surface spectra are a probe of the mass-to-radius ratio. The same formula diverges as R→2MR\to2M: light from the horizon is infinitely redshifted.

Physics of the portrait

Null geodesics became the central objects of modern relativity when a generation realized that the global organization of light cones — not the local metric — is where black holes, singularities, and infinity actually live. No one did more to make that shift than Roger Penrose, who turned the light cone from a picture into a calculus.

Origins · 1 minRoger Penrose and how this was discovered
RP
Roger Penrose
b. 1931 · mathematical physicist who rebuilt causal structure from null geodesics and conformal methods

Penrose recognized that the causal skeleton of a spacetime is carried entirely by its null geodesics, and that this skeleton is conformally invariant — facts he forged into two of the most consequential tools in general relativity. In the early 1960s he introduced the conformal treatment of infinity, compactifying an infinite spacetime by a rescaling g~=Ω2g\tilde g=\Omega^2 g that keeps all light cones at 45°, so that outgoing radiation could be read off at a sharp boundary he named null infinity, or "scri" (I\mathscr{I}); the same construction gives the Penrose–Carter diagrams that put a black hole's global structure on a single page. Then in 1965, in a three-page paper of astonishing economy, he combined the focusing of null-geodesic congruences with the topology of a trapped surface to prove that gravitational collapse must produce a singularity, with no symmetry assumed. This was the first of the singularity theorems; extended with Stephen Hawking, it showed the Big Bang and black-hole singularities to be generic features of Einstein's theory. Penrose also introduced spinor and twistor methods, the Penrose process, and the still-open cosmic-censorship conjecture, and shared the 2020 Nobel Prize in Physics for showing that black-hole formation is a robust prediction of general relativity. Nearly every idea in this article — that light cones generate causal structure, survive conformal rescaling, and focus into singularities — is one he set in its modern form.

History
1908
— Minkowski introduces spacetime and the light cone in his Cologne lecture "Raum und Zeit"
1955
— Raychaudhuri derives the focusing equation for geodesic congruences
1963
— Penrose brings infinity to a finite conformal boundary, defining null infinity
1965
— Penrose's singularity theorem uses null-geodesic focusing and trapped surfaces
2019
— the Event Horizon Telescope images the photon ring around M87\*

Going further

Texts. Misner, Thorne & Wheeler, Gravitation, §22 (geodesics) and §34 (causal structure) · Wald, General Relativity, §3.3, §8, §9 (causal structure, Raychaudhuri, singularity theorems) · Carroll, Spacetime and Geometry, §3.3 (geodesics), §5.4 (Schwarzschild null orbits) · Hartle, Gravity, §8–9 (geodesics, light bending) · Schutz, A First Course in General Relativity, §6 · Hawking & Ellis, The Large Scale Structure of Space-Time, §4 (null congruences). Key results collected here for reference: the affine geodesic equation and null first integral follow from the energy action; the conformal transformation of the connection leaves a surviving ∝kμ\propto k^\mu source that nullity eliminates; the null Raychaudhuri equation yields the reciprocal focusing bound θ−1≥θ0−1+12λ\theta^{-1}\ge\theta_0^{-1}+\tfrac12\lambda; the Schwarzschild photon sphere sits at r=3Mr=3M with bcrit=33 M≈5.196 Mb_{\rm crit}=3\sqrt3\,M\approx5.196\,M, giving shadow angular radii ≈26 μas\approx26\ \mu\mathrm{as} (Sgr A\, M=4.15×106 M⊙M=4.15\times10^6\,M_\odot, D=8.15 kpcD=8.15\ \mathrm{kpc}) and ≈20 μas\approx20\ \mu\mathrm{as} (M87\, M=6.5×109 M⊙M=6.5\times10^9\,M_\odot, D=16.8 MpcD=16.8\ \mathrm{Mpc}); and the gravitational redshift 1+z=(1−2GM/Rc2)−1/21+z=(1-2GM/Rc^2)^{-1/2} gives z≈2.1×10−6z\approx2.1\times10^{-6} (Sun) and z≈0.31z\approx0.31 (1.4 M⊙M_\odot, 10 km neutron star).

Prerequisites. Lorentzian Metric and the Geodesic Equation · Christoffel Symbols and the Covariant Derivative · Causal Structure of Spacetime. Leads to. Chronological and Causal Futures and Pasts · Achronal Sets and Edges · Conformal Transformations and Causal Invariance · Penrose–Carter Conformal Diagrams · Null Infinity (Scri) · Global Hyperbolicity.

Ask the tutor: Want to see how a converging beam of light is forced to a caustic, or why a black hole's shadow is set by rays that never touch it — ask, and we will trace a null geodesic through a real spacetime together.

Why can't proper time parametrize a light ray?

Along a null curve the line element vanishes identically: ds2=gμνdxμdxν=0\mathrm{d}s^2=g_{\mu\nu}\mathrm{d}x^\mu\mathrm{d}x^\nu=0 between any two of its points, so proper time dτ\mathrm{d}\tau never advances — a photon carries no clock of its own. You must use an affine parameter λ\lambda, which is not borrowed from any ruler but supplied by the geometry through the demand that the tangent parallel-transport, kν∇νkμ=0k^\nu\nabla_\nu k^\mu=0. It is fixed only up to λ→aλ+b\lambda\to a\lambda+b, since there is no normalization gμνkμkν=±1g_{\mu\nu}k^\mu k^\nu=\pm1 to pin the scale as there is for massive or spacelike geodesics. Physically the affine parameter is the carrier of frequency: comparing the tangent against observers at emission and reception yields the redshift.

What does it mean that null geodesics are conformally invariant, and why does it matter?

It means that if you multiply the metric by any positive function, gμν→Ω2gμνg_{\mu\nu}\to\Omega^2 g_{\mu\nu}, the very same curves remain null geodesics — only their affine parameter rescales. The reason is that the one term in the transformed geodesic equation that could bend the ray is multiplied by gμνkμkνg_{\mu\nu}k^\mu k^\nu, which is zero for light. It matters because it means causality — which events can send signals to which — depends only on the metric up to scale. That is what lets a Penrose–Carter diagram compress an infinite spacetime onto a finite page while keeping every light cone at exactly 45°: the rescaling changes distances and times drastically but leaves the causal skeleton, carried by the null geodesics, perfectly intact.

Is the light cone the same thing as the event horizon?

No, though they are intimately related. A light cone is a local object: at every event, the null vectors there define which nearby events can be causally influenced. An event horizon is a global object: the boundary of the region from which null geodesics can ever escape to infinity — you cannot know where it is without following light rays all the way out. The horizon is in fact a special null surface, generated by null geodesics that neither escape nor fall in, so it is built from light rays. But whereas a light cone exists at every event in every spacetime, a horizon exists only when the global arrangement of cones traps a region. Confusing the two is a common source of error in reading black-hole diagrams.

Do gravitational waves travel on null geodesics too?

In the high-frequency (geometric-optics) limit, yes. A gravitational wave packet of wavelength short compared to the background curvature propagates along null geodesics of that background, exactly like light, with its own eikonal phase obeying gμν∂μS ∂νS=0g^{\mu\nu}\partial_\mu S\,\partial_\nu S=0. This is why the multi-messenger event GW170817 could be used to compare the arrival of gravitational and electromagnetic signals: both traversed essentially the same null cone across 130 million light-years, arriving within seconds and bounding the fractional difference in their speeds to about 10−1510^{-15}. Outside the short-wavelength limit gravitational waves also backscatter off curvature (tails), so they are not purely null, but the leading propagation is along the same light cones that carry electromagnetic radiation.

Why does gravity always focus light rather than defocus it?

Because of the sign structure of the Raychaudhuri equation. Every term on the right except the twist is non-positive: −12θ2-\tfrac12\theta^2, the shear −σμνσμν-\sigma_{\mu\nu}\sigma^{\mu\nu}, and the curvature −Rμνkμkν-R_{\mu\nu}k^\mu k^\nu. The last is non-negative precisely when matter obeys the null energy condition, Rμνkμkν≥0R_{\mu\nu}k^\mu k^\nu\ge0, which via Einstein's equations means energy densities as seen along a light ray are not negative — true for all ordinary matter and radiation. So an initially parallel or converging surface-forming beam can only be driven toward convergence; a passing mass acts like a converging lens, never a diverging one. Defocusing would require negative energy (a violation of the condition), which is why exotic energy is the ingredient every traversable-wormhole and warp-drive proposal must smuggle in.

How is a photon's frequency defined if it has no clock of its own?

Frequency is not a property of the ray alone — it belongs to whoever measures it. Take the ray's tangent kμk^\mu and an observer's 4-velocity uμu^\mu; the projection ω=−kμuμ\omega=-k_\mu u^\mu is the angular frequency that observer records. Because the affine parameter is arbitrary up to a scale, the absolute value of ω\omega is not physical, but the ratio between two observers is: ωrec/ωemit=(kμuμ)rec/(kμuμ)emit\omega_{\rm rec}/\omega_{\rm emit}=(k_\mu u^\mu)_{\rm rec}/(k_\mu u^\mu)_{\rm emit}, the redshift, is completely well-defined. In a static spacetime a conserved E=−kμξμE=-k_\mu\xi^\mu (from the timelike Killing vector) makes this ratio computable as (−gtt)emit/(−gtt)rec\sqrt{(-g_{tt})_{\rm emit}/(-g_{tt})_{\rm rec}} — the gravitational redshift. So a photon needs no clock; the observers supply the ticking.

Before this

Where this leads