Skip to content

Raychaudhuri Equations

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

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

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

The Raychaudhuri equations evolve the expansion of a smooth timelike or null congruence. They are local geometric identities: focusing follows only after the parameterization, twist, curvature convergence, and initial expansion have been fixed. The complete geometric derivation is in The Raychaudhuri Identity; the physical assumptions are developed in Raychaudhuri Equations for Timelike and Null Congruences.

In DD dimensions let uaua=−1u^au_a=-1, and define

hab=gab+uaub,aa=ub∇bua,Bab=hachbd∇duc.h_{ab}=g_{ab}+u_au_b, \qquad a^a=u^b\nabla_bu^a, \qquad B_{ab}=h_a{}^ch_b{}^d\nabla_du_c.

The projected derivative has the unique decomposition

Bab=θD−1hab+σab+ωab,θ=habBab=∇aua,σab=B(ab)−θD−1hab,ωab=B[ab].\begin{aligned} B_{ab}&=\frac{\theta}{D-1}h_{ab}+\sigma_{ab}+\omega_{ab},\\ \theta&=h^{ab}B_{ab}=\nabla_au^a,\\ \sigma_{ab}&=B_{(ab)}-\frac{\theta}{D-1}h_{ab}, &\omega_{ab}&=B_{[ab]}. \end{aligned}

Thus σab\sigma_{ab} is spatial, symmetric, and trace free, while ωab\omega_{ab} is spatial and antisymmetric. Brackets have unit weight, and the norm convention has no extra factor of one half: σ2=σabσab\sigma^2=\sigma_{ab}\sigma^{ab} and ω2=ωabωab\omega^2=\omega_{ab}\omega^{ab}. Because the rest metric is positive definite, both norms are nonnegative.

With the site curvature convention, the full timelike equation is

ua∇aθ=−θ2D−1−σabσab+ωabωab−Rabuaub+∇aaa.u^a\nabla_a\theta =-\frac{\theta^2}{D-1} -\sigma_{ab}\sigma^{ab} +\omega_{ab}\omega^{ab} -R_{ab}u^au^b +\nabla_aa^a.

For geodesics aa=0a^a=0; for a locally hypersurface-orthogonal congruence ωab=0\omega_{ab}=0. These are independent restrictions. In particular, the nongeodesic contribution is the divergence ∇aaa\nabla_aa^a, not the acceleration magnitude aaaaa_aa^a. The geodesic specialization and its curvature-convergence sign are checked against Guth and Turner 2017, Lecture 14, printed pp. 59–61, Eqs. (14.10)–(14.18).

Let kak^a be null and choose an auxiliary null vector ℓa\ell^a with k⋅ℓ=−1k\cdot\ell=-1. For D≥3D\ge3,

qab=gab+kaℓb+ℓakb,Bab=qacqbd∇dkcq_{ab}=g_{ab}+k_a\ell_b+\ell_ak_b, \qquad \mathcal B_{ab}=q_a{}^cq_b{}^d\nabla_dk_c

projects onto the positive-definite (D−2)(D-2)-dimensional screen. Define

Bab=θD−2qab+σab+ωab,θ=qabBab,\mathcal B_{ab} =\frac{\theta}{D-2}q_{ab}+\sigma_{ab}+\omega_{ab}, \qquad \theta=q^{ab}\mathcal B_{ab},

with screen-trace-free symmetric shear and antisymmetric twist. If

kb∇bka=κka,k^b\nabla_bk^a=\kappa k^a,

then the nonaffine null equation is

ka∇aθ=κθ−θ2D−2−σabσab+ωabωab−Rabkakb.k^a\nabla_a\theta =\kappa\theta -\frac{\theta^2}{D-2} -\sigma_{ab}\sigma^{ab} +\omega_{ab}\omega^{ab} -R_{ab}k^ak^b.

An affine parameter has κ=0\kappa=0. Both θ\theta and κ\kappa change when the generator is rescaled, so an affine-distance bound is meaningful only after fixing an affine normalization. Null hypersurface generators are twist free before caustics. The affine specialization and curvature-convergence term agree with Kontou and Sanders 2020, § 2, p. 8, Eqs. (16)–(17); the +κθ+\kappa\theta term is retained before specializing.

For a geodesic, twist-free timelike congruence satisfying Rabuaub≥0R_{ab}u^au^b\ge0,

dθdτ≤−θ2D−1.\frac{d\theta}{d\tau}\le-\frac{\theta^2}{D-1}.

If θ0<0\theta_0\lt0, a focal point occurs within proper time (D−1)/∣θ0∣(D-1)/\lvert\theta_0\rvert, provided the smooth congruence exists that far. Equality throughout an interval requires σab=0\sigma_{ab}=0 and Rabuaub=0R_{ab}u^au^b=0 in addition to the stated geodesic and twist-free hypotheses. The equality solution is

θ(τ)=θ01+θ0(τ−τ0)/(D−1).\theta(\tau) =\frac{\theta_0} {1+\theta_0(\tau-\tau_0)/(D-1)}.

For an affinely parametrized, twist-free null congruence satisfying Rabkakb≥0R_{ab}k^ak^b\ge0, replace D−1D-1 by D−2D-2. Initial θ0<0\theta_0\lt0 then gives the affine bound (D−2)/∣θ0∣(D-2)/\lvert\theta_0\rvert; equality additionally requires zero shear and Rabkakb=0R_{ab}k^ak^b=0. No such bound follows from the displayed nonaffine equation without first controlling κ\kappa.

Curvature convergence is not itself an energy condition. In Einstein gravity with D>2D\gt2,

Rabkakb=8πGDTabkakb,Rabuaub=8πGD(Tabuaub+TD−2)−2ΛD−2.\begin{aligned} R_{ab}k^ak^b&=8\pi G_D T_{ab}k^ak^b,\\ R_{ab}u^au^b &=8\pi G_D\left(T_{ab}u^au^b+\frac{T}{D-2}\right) -\frac{2\Lambda}{D-2}. \end{aligned}

Thus the null energy condition implies null convergence, whereas the matter strong energy condition implies timelike convergence only when the cosmological term is absent or separately controlled. A focal point is degeneration of the congruence’s transverse Jacobi map; it is not by itself curvature blow-up or geodesic incompleteness Kontou and Sanders 2020, § 5.1, printed pp. 37–42.

For an ingoing shear-free light bundle in DD-dimensional Minkowski spacetime, let its area scale as A∝(λ∗−λ)D−2A\propto(\lambda_*-\lambda)^{D-2}. Using θ=dln⁡A/dλ\theta=d\ln A/d\lambda gives

θ=−D−2λ∗−λ,dθdλ=−θ2D−2,\theta=-\frac{D-2}{\lambda_*-\lambda}, \qquad \frac{d\theta}{d\lambda} =-\frac{\theta^2}{D-2},

so the null coefficient, sign, equality solution, and finite affine focusing time are all fixed while curvature remains zero.

For a comoving flow in a DD-dimensional FLRW spacetime with scale factor A(t)A(t), u=∂tu=\partial_t, H=A˙/AH=\dot A/A, and

θ=(D−1)H,σab=ωab=aa=0,Rtt=−(D−1)(H˙+H2).\theta=(D-1)H, \qquad \sigma_{ab}=\omega_{ab}=a^a=0, \qquad R_{tt}=-(D-1)(\dot H+H^2).

Substitution reduces both sides to (D−1)H˙(D-1)\dot H, checking the timelike coefficient and Ricci sign. As a nongeodesic flat-space check, Rindler observers have zero expansion, shear, and twist but nonzero acceleration; ∇aaa=0\nabla_aa^a=0, so every term still vanishes even though aaaaa_aa^a does not.

  • Guth, Alan H., lectures with notes by Andrew P. Turner. 8.962 General Relativity, Spring 2017. Massachusetts Institute of Technology, 2017. Open PDF.
  • Kontou, Eleni-Alexandra, and Ko Sanders. “Energy Conditions in General Relativity and Quantum Field Theory.” Classical and Quantum Gravity 37 (2020) 193001. DOI; arXiv record.