Skip to content

Raychaudhuri Evolution, Null Focusing, and Renormalized Stress

Raychaudhuri’s equation converts local curvature along a null congruence into the evolution of its cross-sectional area. With quantum matter, the equation itself remains geometric; what changes is that the renormalized TkkT_{kk} need not be pointwise nonnegative. A focusing conclusion therefore requires the initial expansion, shear, twist, affine normalization, stress prescription, and approximation order—not only the sign of one energy pulse.

Required background. Classical energy conditions and quantum violations supplies TkkT_{kk}; Levi–Civita connection and curvature supplies geodesic deviation; and globally hyperbolic spacetimes supplies the causal setting.

Helpful background. Null-smeared stress observables fixes the quantum source, while relativistic causality clarifies which later global inference requires more than local focusing.

Let kμ=dxμ/dλk^\mu=dx^\mu/d\lambda be tangent to an affinely parametrized null congruence in dd dimensions. The screen-space deformation decomposes into expansion θ\theta, shear σμν\sigma_{\mu\nu}, and twist ωμν\omega_{\mu\nu}. With the site Riemann convention,

dθdλ=θ2d2σμνσμν+ωμνωμνRμνkμkν.\frac{d\theta}{d\lambda} =-\frac{\theta^2}{d-2} -\sigma_{\mu\nu}\sigma^{\mu\nu} +\omega_{\mu\nu}\omega^{\mu\nu} -R_{\mu\nu}k^\mu k^\nu.

For generators orthogonal to a smooth null hypersurface, ωμν=0\omega_{\mu\nu}=0. With

Gμν+Λgμν=8πGTμν,G_{\mu\nu}+\Lambda g_{\mu\nu} =8\pi G\,T_{\mu\nu},

null contraction removes both the trace and cosmological terms:

Rkk=8πGTkk.R_{kk}=8\pi G\,T_{kk}.

Thus NEC and nonzero shear make θ\theta nonincreasing. If θ(λ0)<0\theta(\lambda_0)<0 and NEC holds thereafter, dropping the other nonpositive terms gives

dθdλθ2d2,\frac{d\theta}{d\lambda} \le-\frac{\theta^2}{d-2},

so θ\theta diverges to -\infty no later than

λλ0=d2θ(λ0).\lambda-\lambda_0 =\frac{d-2}{\lvert\theta(\lambda_0)\rvert}.

This is a conjugate-point estimate, not by itself a singularity theorem. Geodesic incompleteness additionally needs global causal and genericity hypotheses; the separation between the local focusing lemma and the global causal argument is explicit in Wald 1984, Theorem 9.3.5 and § 9.5.

The structure map places Raychaudhuri between a renormalized null stress and any later causal consequence. Sampling or entropy conditions enter only when pointwise NEC is unavailable.

Initial expansion, shear, affine generator, and renormalized null stress determine Raychaudhuri evolution before any qualified focusing or global causal inference

Null focusing with quantum matter. The map is schematic and not to scale; Raychaudhuri evolution is exact geometrically, while replacing curvature by a mean renormalized stress and deriving a global conclusion each require separate hypotheses.

Take a twist-free, initially parallel congruence with negligible shear and

θ()=0,Tkk(λ)ren=Aπτeλ2/τ2,dλTkk=A.\theta(-\infty)=0, \qquad \langle T_{kk}(\lambda)\rangle_{\mathrm{ren}} =\frac{\mathcal A}{\sqrt{\pi}\tau} e^{-\lambda^2/\tau^2}, \qquad \int_{-\infty}^{\infty}d\lambda\,T_{kk}=\mathcal A.

Assume 8πGAτ18\pi G\lvert\mathcal A\rvert\tau\ll1, so the expansion generated by the pulse changes little over its own width and θ2\theta^2 is perturbative. Integrating Raychaudhuri gives

Δθ=8πGA+O ⁣((GA)2τ,dλσμνσμν),\Delta\theta =-8\pi G\mathcal A +O\!\left((G\mathcal A)^2\tau, \int d\lambda\,\sigma_{\mu\nu}\sigma^{\mu\nu}\right),

A positive pulse produces negative expansion; a negative pulse defocuses at this order. This is an initial-value result for the mean geometry. It does not say that an isolated negative pulse is an allowed complete quantum history; a QEI or ANEC statement may constrain the state and accompanying stress elsewhere.

After a positive pulse, suppose Tkk=0T_{kk}=0 and shear remains zero. The exact vacuum evolution is

θ(λ)=θ+1+θ+(λλ+)/(d2).\theta(\lambda) =\frac{\theta_+} {1+\theta_+(\lambda-\lambda_+)/(d-2)}.

For θ+<0\theta_+<0, a caustic occurs after affine interval (d2)/θ+(d-2)/\lvert\theta_+\rvert. If shear is generated by the pulse, σ20\sigma^2\ge0 moves the caustic earlier. Omitting shear can therefore overestimate, not underestimate, the available affine distance.

Under λ=aλ+b\lambda'=a\lambda+b with a>0a>0,

kμ=kμa,θ=θa,A=dλTkk=Aa.k'^\mu=\frac{k^\mu}{a}, \qquad \theta'=\frac{\theta}{a}, \qquad \mathcal A' =\int d\lambda'\,T_{k'k'} =\frac{\mathcal A}{a}.

The jump law becomes Δθ=8πGA\Delta\theta'=-8\pi G\mathcal A', and the affine distance to the caustic scales by aa. The spacetime event is unchanged. Quoting A\mathcal A or a focusing length without the normalization of kμk^\mu is therefore meaningless.

A second adversarial run sets the initial shear to a localized nonzero profile but solves the shear-free equation. Comparing with the full equation exposes a missed negative contribution. The correction is physical and cannot be repaired by rescaling λ\lambda.

The failure map highlights these two errors: hidden affine normalization changes numerical bounds, whereas omitted shear changes the actual congruence.

A null focusing estimate fails when affine normalization is hidden, shear is omitted, the stress prescription is unspecified, or a local caustic estimate is promoted to global incompleteness

Failure tests for a Raychaudhuri calculation. The diagram is schematic and not to scale; affine rescaling must transform tangent, expansion, stress projection, integration measure, and endpoints together.

See the chapter domain and failure-conditions table. The pulse calculation assumes an affine, twist-free congruence, a fixed renormalized mean stress, negligible shear at the retained order, and weak backreaction during the pulse. It fails for nonaffine generators unless the nonaffinity term is restored, near caustics where linearization fails, or when fluctuation effects invalidate a mean-stress geometry.

Integrate dθ/dλ=θ2/(d2)d\theta/d\lambda=-\theta^2/(d-2) from θ(0)=θ0<0\theta(0)=\theta_0<0 and locate the caustic.

Solution

Separation gives d(1/θ)/dλ=1/(d2)d(1/\theta)/d\lambda=1/(d-2), hence

1θ(λ)=1θ0+λd2.\frac1{\theta(\lambda)} =\frac1{\theta_0}+\frac{\lambda}{d-2}.

The denominator vanishes at λc=(d2)/θ0\lambda_c=(d-2)/\lvert\theta_0\rvert.

  • Wald, R. M. General Relativity. University of Chicago Press, 1984. Publisher.