Raychaudhuri equation

In general relativity, the Raychaudhuri equation, or Landau–Raychaudhuri equation, is a fundamental result describing the motion of nearby bits of matter.

The equation is important as a fundamental lemma for the Penrose–Hawking singularity theorems and for the study of exact solutions in general relativity, but has independent interest, since it offers a simple and general validation of our intuitive expectation that gravitation should be a universal attractive force between any two bits of mass–energy in general relativity, as it is in Newton's theory of gravitation.

The equation was discovered independently by the Indian physicist Amal Kumar Raychaudhuri and the Soviet physicist Lev Landau.

Mathematical statement
Given a timelike unit vector field $$\vec{X}$$ (which can be interpreted as a family or congruence of nonintersecting world lines via the integral curve, not necessarily geodesics), Raychaudhuri's equation can be written
 * $$\dot{\theta} = - \frac{\theta^2}{3} - 2 \sigma^2 + 2 \omega^2 - {E[\vec{X}]^a}_a + _{;a}$$

where
 * $$2\sigma^2 = \sigma_{mn} \, \sigma^{mn}, \; 2\omega^2 = \omega_{mn} \, \omega^{mn}$$

are (non-negative) quadratic invariants of the shear tensor
 * $$\sigma_{ab} = \theta_{ab} - \frac{1}{3} \, \theta \, h_{ab}$$

and the vorticity tensor
 * $$\omega_{ab} = {h^m}_a \, {h^n}_b X_{[m;n]}$$

respectively. Here,
 * $$\theta_{ab} = {h^m}_a \, {h^n}_b X_{(m;n)}$$

is the expansion tensor, $$\theta$$ is its trace, called the expansion scalar, and
 * $$h_{ab} = g_{ab} + X_a \, X_b$$

is the projection tensor onto the hyperplanes orthogonal to $$\vec{X}$$. Also, dot denotes differentiation with respect to proper time counted along the world lines in the congruence. Finally, the trace of the tidal tensor $$E[\vec{X}]_{ab}$$ can also be written as
 * $${E[\vec{X}]^a}_{a} = R_{mn} \, X^m \, X^n$$

This quantity is sometimes called the Raychaudhuri scalar.

Intuitive significance
The expansion scalar measures the fractional rate at which the volume of a small ball of matter changes with respect to time as measured by a central comoving observer (and so it may take negative values). In other words, the above equation gives us the evolution equation for the expansion of the timelike congruence. If the derivative (with respect to proper time) of this quantity turns out to be negative along some world line (after a certain event), then any expansion of a small ball of matter (whose center of mass follows the world line in question) must be followed by recollapse. If not, continued expansion is possible.

The shear tensor measures any tendency of an initially spherical ball of matter to become distorted into an ellipsoidal shape. The vorticity tensor measures any tendency of nearby world lines to twist about one another (if this happens, our small blob of matter is rotating, as happens to fluid elements in an ordinary fluid flow which exhibits nonzero vorticity).

The right hand side of Raychaudhuri's equation consists of two types of terms:
 * 1) terms which promote (re)-collapse
 * 2) * initially nonzero expansion scalar,
 * 3) * nonzero shearing,
 * 4) * positive trace of the tidal tensor; this is precisely the condition guaranteed by assuming the strong energy condition, which holds for the most important types of solutions, such as physically reasonable fluid solutions,
 * 5) terms which oppose (re)-collapse
 * 6) * nonzero vorticity, corresponding to Newtonian centrifugal forces,
 * 7) * positive divergence of the acceleration vector (e.g., outward pointing acceleration due to a spherically symmetric explosion, or more prosaically, due to body forces on fluid elements in a ball of fluid held together by its own self-gravitation).

Usually one term will win out. However, there are situations in which a balance can be achieved. This balance may be:
 * stable: in the case of hydrostatic equilibrium of a ball of perfect fluid (e.g. in a model of a stellar interior), the expansion, shear, and vorticity all vanish, and a radial divergence in the acceleration vector (the necessary body force on each blob of fluid being provided by the pressure of surrounding fluid) counteracts the Raychaudhuri scalar, which for a perfect fluid is $$E[\vec{X}]^a{}_a = 4 \pi ( \mu + 3 p )$$ in geometrized units. In Newtonian gravitation, the trace of the tidal tensor is $$4 \pi \mu$$; in general relativity, the tendency of pressure to oppose gravity is partially offset by this term, which under certain circumstances can become important.
 * unstable: for example, the world lines of the dust particles in the Gödel solution have vanishing shear, expansion, and acceleration, but constant vorticity just balancing a constant Raychuadhuri scalar due to nonzero vacuum energy ("cosmological constant").

Focusing theorem
Suppose the strong energy condition holds in some region of our spacetime, and let $$\vec{X}$$ be a timelike geodesic unit vector field with vanishing vorticity, or equivalently, which is hypersurface orthogonal. For example, this situation can arise in studying the world lines of the dust particles in cosmological models which are exact dust solutions of the Einstein field equation (provided that these world lines are not twisting about one another, in which case the congruence would have nonzero vorticity).

Then Raychaudhuri's equation becomes
 * $$\dot{\theta} = - \frac{\theta^2}{3} - 2 \sigma^2 - {E[\vec{X}]^a}_a$$

Now the right hand side is always negative or zero, so the expansion scalar never increases in time.

Since the last two terms are non-negative, we have
 * $$\dot{\theta} \leq - \frac{\theta^2}{3}$$

Integrating this inequality with respect to proper time $$\tau$$ gives
 * $$\frac{1}{\theta} \geq \frac{1}{\theta_0} + \frac{\tau}{3}$$

If the initial value $$\theta_0$$ of the expansion scalar is negative, this means that our geodesics must converge in a caustic ($$\theta$$ goes to minus infinity) within a proper time of at most $$3/|\theta_0|$$ after the measurement of the initial value $$\theta_0$$ of the expansion scalar. This need not signal an encounter with a curvature singularity, but it does signal a breakdown in our mathematical description of the motion of the dust.

Optical equations
There is also an optical (or null) version of Raychaudhuri's equation for null geodesic congruences.
 * $$\dot{\widehat{\theta}} = - \frac{1}{2}\widehat{\theta}^2 - 2 \widehat{\sigma}^2 + 2 \widehat{\omega}^2 - T_{\mu\nu} U^\mu U^\nu$$.

Here, the hats indicate that the expansion, shear and vorticity are only with respect to the transverse directions. When the vorticity is zero, then assuming the null energy condition, caustics will form before the affine parameter reaches $$2/\widehat{\theta}_0$$.

Applications
The event horizon is defined as the boundary of the causal past of null infinity. Such boundaries are generated by null geodesics. The affine parameter goes to infinity as we approach null infinity, and no caustics form until then. So, the expansion of the event horizon has to be nonnegative. As the expansion gives the rate of change of the logarithm of the area density, this means the event horizon area can never go down, at least classically, assuming the null energy condition.