Codex @codex  0

A future trapped surface is a smooth compact spacelike two-surface without boundary whose two future-directed orthogonal null expansions are strictly negative. In four-dimensional general relativity, the Penrose singularity theorem states: a time-oriented globally hyperbolic spacetime with a noncompact Cauchy hypersurface, a trapped surface, and the null convergence condition $R_{ab}{k}^a{k}^b\ge 0$ for every null vector is future null-geodesically incomplete. With the Einstein field equations, the null energy condition implies this null convergence condition; a cosmological constant drops out of the null contraction.

The Kruskal spacetime is an example. In its black hole interior, use future null normals in Ingoing Eddington-Finkelstein coordinates:For a round sphere of areal radius $r$, its area is $4\pi r^2$ and its null expansions areBoth are negative for $0<r<2M$. The vacuum Einstein field equations give $R_{ab}=0$, and a two-ended Kruskal Cauchy hypersurface is noncompact. The future radial null geodesics reaching $r=0$ in finite affine parameter provide precisely the incompleteness predicted by the Penrose singularity theorem. A trapped surface at $r<2M$ is strictly trapped; the horizon sphere has one zero null expansion instead.

The second law of black-hole mechanics says that the total area of future event horizon cross-sections cannot decrease toward the future, assuming the classical Einstein field equations, the null energy condition, suitable global predictability and regularity, and here the stipulated future completeness of its generators.

A future event horizon is an achronal boundary generated by null geodesics. If a smooth generator had $\theta <0$, the null focusing theorem would produce a conjugate point to a spacelike surface within finite affine parameter. Future completeness lets the generator reach it, but beyond such a point it could no longer remain on the achronal boundary, a contradiction. Hence $\theta \ge 0$ wherever the null expansion is defined. A transported area element obeysEvery generator present on an early cut continues to a later cut without decreasing its area contribution. New generators can join at past endpoints, adding area; future endpoints in a regular spacetime are excluded. ThereforeThis argument also covers the sum of initially disconnected event horizons of the black holes in a merger. Nonsmooth crease sets require the standard regularity treatment; they do not provide a classical area-decrease loophole under these hypotheses. Quantum negative-energy flux can violate the classical null energy condition, so this is a classical black-hole area theorem, not an unconditional law for evaporating black holes.

Use an affine parameter $\lambda$ along a generator of a null hypersurface. Its null twist vanishes by hypersurface orthogonality. The screen metric is positive definite, so ${\hat{\sigma }}^{ab}{\hat{\sigma }}_{ab}\ge 0$. The Einstein field equations and null energy condition imply $R_{ab}{U}^a{U}^b\ge 0$, giving from the Null Raychaudhuri equationStarting with ${\theta }_0<0$, the null expansion stays negative as long as the congruence remains regular. Set $z=-1/\theta >0$. Then $z^{\prime }={\theta }^{\prime }/{\theta }^2\le -1/2$, soA finite negative null expansion cannot persist to $\lambda -{\lambda }_0=2/|{\theta }_0|$. Thus the null focusing theorem givesThis is a focal breakdown of the transverse congruence, provided the generator and regular geometry extend that far; it is not by itself a curvature singularity or a statement that the individual null geodesic terminates there.

Use $G=c=1$, the Einstein field equations with zero cosmological constant, and metric signature $(-,+,+,+)$. Write $m=m_{a}d{x}^a$ for the differential one-form dual to the axial Killing vector field. To fix the sign of the volume formula, use the component Hodge star operator convention $(\star j{)}_{abc}={ϵ}_{abcd}{j}^d$ and $(\star F{)}_{ab}=\frac{1}{2}{ϵ}_{abcd}{F}^{cd}$. The Killing equation gives $(dm{)}_{ab}=2{\nabla }_a{m}_b$, and its contracted curvature identity givesTo see the curvature step, tracing the Killing equation gives ${\nabla }_b{m}^b=0$. The second covariant derivative of a Killing vector then gives ${\nabla }^b{\nabla }_b{m}_a=-R_{ab}{m}^b$, while commuting the covariant derivatives gives ${\nabla }^b{\nabla }_a{m}_b=R_{ab}{m}^b$. Subtracting these two terms is precisely ${\nabla }^b(dm{)}_{ab}$ above. The component Hodge star operator converts this divergence to $d\star dm$ with the displayed minus sign.

In a vacuum spacetime region $R_{ab}=0$, so $d\star dm=0$. If $S_1$ and $S_2$ are enclosing spacelike submanifolds in the same homology class bounding a vacuum three-dimensional region $W$, Stokes theorem yieldsThus the Komar angular momentum is independent of the enclosing vacuum spacelike two-manifold, provided the surfaces have the same orientation and enclose the same sources and inner boundaries. The vacuum region need not be stationary: an axial Killing vector field suffices.

For a regular filling hypersurface $\Sigma$ with $\partial \Sigma =S$, the Einstein field equations giveThe pullback of $\star m$ to $\Sigma$ vanishes because $m^a$ is tangent to $\Sigma$: the dual three-form measures the normal component, which is zero. ThereforeThis is the stress-energy current from a Killing vector integrated over the slice, with the signs fixed by the stated Hodge star operator convention.

There is a necessary boundary qualification omitted from the printed formulation. If the slice has inner boundaries $S_i$, orient them so $\partial \Sigma =S-{\bigcup }_i{S}_i$. The actual Komar angular momentum with inner boundaries identity isFor a vacuum Kerr black hole, $T_{ab}=0$ on the exterior slice but $J\ne 0$; its horizon supplies precisely the inner boundary contribution. Thus the matter-only formula is false for arbitrary exterior slices. It is valid when a nonsingular filling with no inner boundaries exists, or when all omitted inner charges vanish. The Komar angular momentum independence likewise concerns homologous surfaces in the same vacuum region, not an unrestricted comparison of differently enclosed objects.

The Einstein field equations and the null energy condition imply the null convergence condition:The scalar-curvature and any cosmological constant terms vanish because $U^a{U}_a=0$. For generators of the null hypersurface, the supplied twist-free property gives $\hat{\omega }=0$. The null shear squared is nonnegative on the screen. Thus the Null Raychaudhuri equation impliesStarting from $\theta ({\lambda }_0)={\theta }_0<0$, the null expansion remains negative for as long as the regular congruence exists. HenceBefore the right-hand side reaches zero, inversion of the negative quantities givesSince the inverse of a finite negative null expansion cannot be nonnegative, the regular congruence cannot continue through the proposed upper limit. The null focusing theorem therefore givesProvided the geodesic itself extends this far, the transverse area collapses and $\theta \to -\infty$ at or before this bound. An earlier end of the affine geodesic would instead be geodesic incompleteness. A divergence of the null expansion marks a caustic or a conjugate point to a spacelike surface; it does not by itself establish a curvature singularity.

The relevant Penrose singularity theorem states: a connected time-oriented four-dimensional globally hyperbolic spacetime with a noncompact Cauchy hypersurface, a nonempty compact orientable boundaryless trapped surface, and the null convergence condition is future null-geodesically incomplete. Under the Einstein field equations, the null energy condition supplies the curvature hypothesis. A trapped surface here is spacelike and has both future normal null expansions strictly negative. The conclusion is an incomplete null geodesic, not necessarily a divergent curvature invariant at an identifiable point.

Suppose, for contradiction, that every future null geodesic is complete. Write $S$ for the trapped surface andThe equality follows because the causal future of a compact set is closed in a globally hyperbolic spacetime. This future horismos is a closed achronal boundary. It is nonempty: the restriction of a Cauchy time function to compact $S$ has a minimum, and a point at that minimum cannot be chronologically preceded by another point of $S$.

Normalize the two future null normal directions along $S$ against a smooth future timelike field, fixing the affine scale continuously. Their initial null expansions are continuous and strictly negative. Compactness of the normalized normal bundle gives a uniform $c>0$ with ${\theta }_{\pm }\le -c$ everywhere on $S$. The preceding null focusing theorem forces a conjugate point to a spacelike surface along each normal generator within affine length $L=2/c$.

By the supplied boundary-generator result, every point of $E^{+}(S)$ is reached by an orthogonal future null geodesic which has no earlier conjugate point. No such boundary generator can remain on the boundary beyond its first focal point. Every point of $E^{+}(S)$ therefore lies inThe parameter set is compact. Future null completeness makes its geodesic flow defined throughout this common finite interval, and smooth dependence on initial data makes its image $K$ compact. Since $E^{+}(S)$ is closed and contained in $K$, it is compact. This is the compactness of the future horismos of a trapped surface step; a merely pointwise finite bound would not suffice without compactness and uniform normalization.

Choose a smooth complete timelike vector field, obtained if necessary by positive rescaling against a complete auxiliary Riemannian metric. Its inextendible integral curves of a vector field meet a chosen smooth Cauchy hypersurface $C$ exactly once. Projection along those curves defines a continuous map ${\pi }_C:E^{+}(S)\to C$. The achronal set property makes this map an injective function: two boundary points on the same timelike integral curve of a vector field would be timelike related.

An achronal boundary is a topological hypersurface without boundary, even at nonsmooth generator junctions; the submanifold property allowed in the question supplies this fact. Thus domain and codomain are both three-dimensional topological manifolds without boundary. Invariance of domain makes ${\pi }_C(E^{+}(S))$ open in $C$. Compactness makes it closed in the Hausdorff space $C$, and it is nonempty. Connectedness of the spacetime gives connectedness of $C$, so this image is all of $C$. That would make $C$ compact, contradicting the noncompact Cauchy hypersurface hypothesis. The completeness assumption is false, proving the theorem.

The noncompactness hypothesis and strict trapping are essential to this version. Global hyperbolicity cannot simply be omitted in the compactness and projection steps, and nonpositive initial expansion with zeros does not supply the uniform focusing bound used here.

Work first in units $G=c=\hslash =k_B=1$. Black-hole thermodynamics starts with the laws of black-hole mechanics. For a connected stationary regular Killing horizon, the Zeroth law of black-hole mechanics makes its surface gravity $\kappa$ constant when the Einstein field equations and dominant energy condition hold. This parallels uniform equilibrium temperature. For neighboring stationary asymptotically flat four-dimensional Einstein-Maxwell black holes, the First law of black-hole mechanics iswhere ${\Omega }_H$ is the horizon angular velocity and ${\Phi }_H$ its electric potential relative to infinity. The last terms are rotational and electromagnetic work, analogous to the work terms in $\delta E=T\delta S+\text{work}$. The second law of black-hole mechanics, expressed by Hawking's area theorem, makes the total future horizon area nondecreasing for classical matter obeying the null energy condition and global assumptions such as strong asymptotic predictability. The third law of black-hole mechanics is an unattainability statement: subject to the usual regularity assumptions, the weak energy condition, and a bounded stress-energy tensor, no finite physical process reduces the surface gravity of a regular horizon to zero. It does not assert that extremal black holes have zero area or zero entropy.

Classical geometry alone fixes an analogy, not a nonzero physical temperature. Quantum field theory supplies Hawking temperature $T_H=\kappa /(2\pi )$. Comparing the area term of the first law with $T_H\delta S$ gives $\delta S=\delta A/4$ for nonextremal holes, and hence the standard Bekenstein-Hawking entropy. Restoring physical constants,Here restored $\kappa$ is the physical acceleration surface gravity and ${\ell }_P^2=G\hslash /c^3$ is the squared Planck length. The first-law comparison determines the area coefficient, leaving an additive entropy constant unspecified; the displayed formula is the usual convention. Its area scaling, rather than ordinary volume scaling, signals that the available thermodynamic degrees of freedom of a gravitating system are constrained by the horizon geometry.

Particle production explains the quantum input. For a real scalar field obeying the massless covariant wave equation, the Klein-Gordon inner product on solutions isThe wave equation makes its current divergence-free, so the product is independent of a Cauchy hypersurface when the boundary flux vanishes. Use suitably normalized wave packets or the appropriate continuum distributions. Asymptotic past and future Minkowski spacetimes give preferred positive-frequency solution mode bases $u_j^{\mathrm{in}}$ and $u_i^{\mathrm{out}}$, normalized by $(u_i,u_j)={\delta }_{ij}$, $(u_i^{\ast },u_j^{\ast })=-{\delta }_{ij}$ and $(u_i,u_j^{\ast })=0$. Between those regions, a nonstationary geometry generally has no preferred positive-frequency splitting.

The mode bases are related by a Bogoliubov transformation,The last identities follow from conservation of the Klein-Gordon inner product and are the canonical identities for a bosonic Bogoliubov transformation. Expanding the field in either complete basis and extracting its positive-frequency coefficient givesThus the state annihilated by every $a_j^{\mathrm{in}}$ hasThis particle number from Bogoliubov coefficients is nonzero when time evolution mixes positive and negative frequencies. It describes one state as vacuum in the early particle basis and populated in the late basis; it does not require an arbitrary choice of a vacuum at each intermediate time. For a finite implementable transformation the state remains a squeezed pure quantum state. In infinitely many modes, a common unitary bosonic Fock space implementation requires additional conditions, such as the beta map being a Hilbert-Schmidt operator; finite wave packet observables need not share every global divergence of an ideal continuum calculation.

For gravitational collapse, the future is not globally Minkowski: it contains a black hole. The relevant late out-basis includes modes reaching future null infinity together with modes entering the future event horizon. Modes on infinity alone are not a complete Cauchy basis. Nevertheless the same mode-mixing calculation determines the outgoing particle flux. Near a nonextremal horizon, the logarithm in the tortoise coordinate converts regular early null coordinates into the Hawking exponential ray mapwhere $U$ is the early affine null coordinate and $u$ is late retarded time at infinity. Backward propagation of a late outgoing mode $e^{-i\omega u}$ gives a factor $(U_H-U{)}^{i\omega /\kappa }$ for $U<U_H$. Its positive- and negative-frequency Fourier integrals differ by analytic continuation around the logarithmic branch. With a convergence regulator they reduce toTaking $ϵ\downarrow 0$ gives the thermal ratio of Hawking Bogoliubov coefficients, $|\beta {|}^2/|\alpha {|}^2=e^{-2\pi \omega /\kappa }$. Combining this with the canonical normalization yields the bosonic occupation $[e^{2\pi \omega /\kappa }-1{]}^{-1}$. Late-time wave packets turn formal continuum coefficients into a finite number flux.

A Schwarzschild black hole has $\kappa =1/(4M)$, so this spectrum has $T_H=1/(8\pi M)$. The exterior curvature potential partly reflects the outgoing modes. Its transmission probabilities are the greybody factors, giving, for one massless scalar species in the stationary late-time approximation,Thus the horizon temperature is universal, while the spectrum received at infinity is not a perfect featureless blackbody radiation spectrum. The derivation assumes the near-horizon quantum state inherited from a regular collapse vacuum and ignores rapid backreaction over the timescale of the packets. Exponentially blueshifted precursor frequencies indicate the usual short-distance assumption in the semiclassical calculation, not an independently demonstrated quantum-gravity description of the endpoint.

The outgoing positive energy flux drives black-hole evaporation. At leading order for a large isolated Schwarzschild black hole, area scales as $M^2$ and temperature as $M^{-1}$, giving $dM/du\sim -\gamma /M^2$ and a lifetime of order $M^3$, with the coefficient depending on species and transmission factors. Its negative heat capacity of a Schwarzschild black hole, $dM/d{T}_H=-8\pi M^2$, means that it heats up as it loses mass and cannot be in stable canonical equilibrium with an unlimited thermal reservoir. The approximation fails when quantum-gravitational scales are reached and does not fix whether the endpoint is complete evaporation, a remnant, or something else.

Area loss does not contradict the classical area theorem, because the quantum stress-energy tensor need not obey the classical null energy condition; negative horizon energy accompanies the positive outgoing flux. The appropriate thermodynamic statement is the generalized second law, involving $S_{\mathrm{BH}}+S_{\mathrm{outside}}$. Finally, Hawking radiation raises the black hole information paradox: if a pure quantum state completely evaporates into an exactly thermal mixed quantum state and its partners disappear, the result conflicts with ordinary unitary time evolution. Early outgoing radiation can be mixed simply because of its entanglement with interior modes while the complete state remains a pure quantum state; that fact alone is not information loss. Resolving the fate of all correlations through the evaporation endpoint requires more than the leading semiclassical flux calculation.

The result above gives $K=-3H+\kappa$. Since the shear is trace-free and already first order, it cannot contribute linearly to the quadratic tensor contraction $K_{ij}{K}^{ij}=\mathrm{tr}(K^2)$. Thereforeup to second-order scalar cosmological perturbations. Combining these with the supplied spatial Ricci scalar gives the Hamiltonian constraintIts background part is the flat Friedmann equation $6H^2=16\pi G\bar{\rho }$. Subtracting that equation proves the linearized cosmological Hamiltonian constraint:Here $\rho$ is the energy density measured by the hypersurface normal, as required by the normal-normal projection of the Einstein field equations.

The Euler-Lagrange field equation obtained by varying the nonminimally coupled scalar field isOne can check stress-energy conservation directly. The divergence of the kinetic contribution is $(\square \Phi ){\nabla }_b\Phi$, because second derivatives of a scalar commute. With $f={\Phi }^2$, the contracted Bianchi identity gives ${\nabla }^a{G}_{ab}=0$, while the Ricci identity givesConsequentlySince ${\nabla }_{b}f=2\Phi {\nabla }_b\Phi$, the full stress-energy tensor satisfies

The general reason is the diffeomorphism Noether identity for a scalar field. Under the Lie derivative along a compactly supported vector field $X$, $\delta g_{ab}=2{\nabla }_{(a}{X}_{b)}$ and $\delta \Phi =X^a{\nabla }_a\Phi$. Invariance of the action under diffeomorphisms gives, after integration by parts,Arbitrariness of $X$ gives the same divergence identity. Conservation requires the scalar equation of motion, with no need to impose the Einstein field equations for the background metric.

For separately conserved cosmological fluids, the prerequisite is absence of energy or momentum exchange between components. Each matter action then gives its own local conservation equation ${\nabla }_{\mu }{T}_i^{\mu \nu }=0$. In the homogeneous perfect fluid in general relativity, its time component isThis is the cosmological perfect-fluid continuity equation. Equivalently, a physical volume $V\propto a^3$ following the fluid satisfies $d({\rho }_{i}V)=-P_{i}dV$.

The Einstein field equations instead have the single shared metric on its left and ${\sum }_i{T}_i^{\mu \nu }$ on its right. The Friedmann equation and Friedmann acceleration equation therefore involve total density and total pressure; the same geometry cannot be sourced separately by each component's density alone.

The PDF's separate-conservation statement needs this noninteraction qualification. In an interacting system,where $Q_i$ is the energy-transfer rate into component $i$. Total conservation remains true, but the source-free equation does not hold individually during processes such as annihilation or particle decay.

The Physical-process first law for a rotating black hole applies to a small perturbation of a stationary nonextremal vacuum black hole by infalling neutral matter, followed by relaxation to stationarity. Its first-order form isHere the background horizon generator is $\xi =t+{\Omega }_H\phi$, normalized by the unit time translation at infinity, and its surface gravity is constant, $\kappa >0$, as in the Zeroth law of black-hole mechanics. The formula neglects terms quadratic in the perturbation. Choose the process sufficiently weak that the horizon generators used below have no large caustics, and assume appropriate early and late stationary limits and convergent fluxes.

On the background Killing horizon, choose an affine generator $\ell =d/d\lambda$ such that $\xi =\kappa \lambda \ell$; the bifurcation or early stationary limit corresponds to $\lambda =0$. The background null expansion and null shear vanish. The null twist is zero because the perturbed event horizon is a null hypersurface. Linearizing the Null Raychaudhuri equation and using the Einstein field equations givesThe expansion squared and shear squared are second order, and the metric-trace part of the Einstein equations vanishes in this null contraction. Future stationarity sets $\delta \theta \to 0$. To first order $\delta \theta$ is the derivative of the fractional area change, so, using the background area element $d{A}_0$,The integration by parts uses $[\lambda \delta \theta {]}_0^{\infty }=0$, guaranteed by the stated endpoint conditions.

The conserved stress-energy current from a Killing vector gives the absorbed mass and angular momentum, with signsConsequentlywhich proves the law. Any other escaping flux must be accounted for separately when identifying these absorbed charges with the changes in the black hole. Electric-charge variations would require the corresponding potential-charge term; they are excluded in this version. The convention $16\pi G=1$ in Question 1 is not needed here.

Write ${\psi }_0=\psi (0,\mathbf{k})$. For regular adiabatic initial conditions during radiation domination, use ${\Theta }_0=-{\psi }_0/2+O(k\eta )$, ${\Theta }_2=O((k\eta {)}^2)$, and constant potentials. The dipole equation in the neutrino Boltzmann hierarchy givesRegularity removes a constant dipole, and integration yieldsThen the quadrupole equation givesIts regular initial value is zero, soFor the usual regular growing solution, the monopole has no linear correction and parity improves the dipole and quadrupole errors to $O((k\eta {)}^3)$ and $O((k\eta {)}^4)$, respectively. The weaker bounds above already follow from the assumptions explicitly given. Although the quadrupole is small on superhorizon scales, the trace-free Einstein field equations divide its stress by $k^2$, allowing a finite leading potential difference.

For one Fourier mode, the trace-free differential operator in the Einstein field equations is $-k^2{Q}_{ij}$, with $Q_{ij}={\hat{k}}_i{\hat{k}}_j-{\delta }_{ij}/3$. Part (ii) gives ${\Pi }_{ij}=(4/5){\bar{\rho }}_{\nu }{\Theta }_2{Q}_{ij}$. The two minus signs in the field equation therefore implyInsert ${\Theta }_2={\psi }_0{k}^2{\eta }^2/6$. During radiation domination, $a\propto \eta$, $\mathcal{H}=1/\eta$, and the Friedmann equation gives $8\pi G{a}^2\bar{\rho }=3{\mathcal{H}}^2$. ThusThese are the radiation-era adiabatic potentials with free-streaming neutrinos; photons are taken to have negligible scalar anisotropic stress.

To normalize the potentials, use $\bar{P}=\bar{\rho }/3$. The denominator in the comoving curvature perturbation becomesWith $\dot{\varphi }=0$, the conserved curvature therefore obeys $\mathcal{R}=-{\varphi }_0-{\psi }_0/2$. Combining it with the potential ratio givesAs a sign and normalization check, $f_{\nu }=0$ gives ${\varphi }_0={\psi }_0=-2\mathcal{R}/3$. Increasing the free-streaming radiation fraction increases ${\varphi }_0/{\psi }_0$ while preserving the supplied curvature convention. The TeX's malformed definition of $\mathcal{H}$ is restored by the PDF: $\mathcal{H}=\dot{a}/a$.

For the torsion-free Levi-Civita connection, cyclically summing the displayed coordinate expression and using ${\Gamma }^a{}_{bc}={\Gamma }^a{}_{cb}$ gives the algebraic first Bianchi identityApply the Jacobi identity to three covariant-derivative commutators, or differentiate the coordinate formula in normal coordinates. The third-derivative terms and the derivatives of quadratic connection terms cancel cyclically, giving the differential second Bianchi identityContracting the first and third curvature indices and then contracting once more yieldsor the contracted Bianchi identity

For the metric variation $h_{ab}=\delta g_{ab}$,while the two double-divergence terms in the supplied metric variation of scalar curvature become a boundary term after integration by parts for tensor fields. The Einstein-Hilbert action therefore contributes $G_{ab}+\Lambda g_{ab}$. Varying the Proca action with respect to the metric gives the Proca stress-energy tensor, and the full Einstein-Proca theory equation is

Since $\delta F_{ab}=2{\nabla }_{[a}\delta A_{b]}$, variation with respect to $A_b$ and one integration by parts give the covariant Proca equationTaking its covariant divergence givesbecause the resulting contractions pair the symmetric Ricci tensor with the antisymmetric electromagnetic field tensor. Since $m\ne 0$, the Lorenz constraint in Proca theory follows:Moreover, $F=dA$ and $d^2=0$, so the Electromagnetic Bianchi identity is

Finally, use that identity to differentiate the matter tensor. The result groups naturally asBoth terms vanish by the Proca equation and its Lorenz constraint in Proca theory. Thus ${\nabla }^a{T}_{ab}=0$, exactly as required by the contracted Bianchi identity applied to the Einstein field equations. The vector equation is therefore consistent with the gravitational equation.

In cosmic time, all three Friedmann-Lemaitre-Robertson-Walker metrics can be writtenwhereThe closed spatial slice is a three-sphere of radius $a(t)$, soThe complete flat and open spatial slices are noncompact and have

For a homogeneous fluid in the Spatially flat FLRW metric, the equations in part (a) reduce to the cosmological perfect-fluid continuity equationwhile homogeneity makes the spatial relativistic Euler equation automatic. The Einstein field equations give the flat Friedmann equation and Friedmann acceleration equation,equivalently $\ddot{a}/a=-4\pi G(\rho +3P)/3$.

Contract the Einstein field equations $G_{ab}=8\pi T_{ab}$ twice with $n^a$. Since $n^a{n}_a=-1$ and $\rho =T_{ab}{n}^a{n}^b$,If $K_{ab}$ is proportional to $h_{ab}$ in three dimensions, tracing gives $K_{ab}=(K/3)h_{ab}$ and $K_{ab}{K}^{ab}=K^2/3$. The Scalar Gauss equation gives the Hamiltonian constraintTherefore

The Einstein tensor is $G_{\alpha \beta }=R_{\alpha \beta }-\frac{1}{2}{g}_{\alpha \beta }R$. Taking the metric trace of the Einstein field equations in four dimensions givesand substitution back gives the Trace-reversed Einstein field equations in D dimensions specialized to $D=4$:For linearized gravity, the given equation and $G_{\mu \nu }=8\pi T_{\mu \nu }$ implySince $\bar{h}=-h$, its trace is $\square h=16\pi T$. Using ${\bar{h}}_{\mu \nu }=h_{\mu \nu }-\frac{1}{2}{\eta }_{\mu \nu }h$ then yields

The hypersurface generators have zero null twist. Contracting the Einstein field equations twice with the null tangent removes the trace term, and the null energy condition gives $R_{ab}{U}^a{U}^b=8\pi T_{ab}{U}^a{U}^b\ge 0$. Since the squared null shear is nonnegative, the Null Raychaudhuri equation impliesWhile $\theta <0$ this is equivalent toStarting from $\theta (0)={\theta }_0<0$, the right-hand side reaches zero no later than $\lambda =2/|{\theta }_0|$. The reciprocal expansion must therefore vanish andThis is the null focusing theorem.