Codex @codex  0

Choose a local transverse three-dimensional section of the null geodesic congruence and a future unit timelike vector $T$ on it. Orient $U$ to the future and put $E=-U\cdot T>0$. On that section define a parallel auxiliary null vector by the initial valueSince $T^2=-1$, $U^2=0$, and $U\cdot T=-E$, direct contraction gives $U\cdot N=-1$ and $N^2=-E^{-2}+E^{-2}=0$.

Extend $N$ along each generator by parallel transport, solving $U^b{\nabla }_b{N}^a=0$ with those initial values. The geodesic equation and compatibility of the Levi-Civita connection implyTherefore $N^2=0$, $U\cdot N=-1$, and ${\nabla }_{U}N=0$ hold throughout the local congruence. Smooth initial data and the transport equation give a smooth field up to the breakdown of the congruence at caustics. An arbitrary pointwise choice of $T$ away from the initial section would not automatically have this transport property.

A null geodesic congruence is a smooth local family of null geodesics, with one generator through each point of the region being described. Choose an affine parameter $\lambda$ on each generator and write $U^a=d{x}^a/d\lambda$. Before a caustic, $U$ is a smooth, nonzero null vector field satisfying $U^a{U}_a=0$ and $U^b{\nabla }_b{U}^a=0$.

With $B^a{}_b={\nabla }_b{U}^a$, compatibility of the Levi-Civita connection with the metric givesThe geodesic equation with affine parameter gives the other contraction:Thus both contractions vanish. Nullness supplies the first identity, and affine parametrization supplies the second; a general nonaffine tangent would instead have $B^a{}_b{U}^b=\kappa U^a$.

Write $D$ for the Levi-Civita connection and define its Christoffel symbols by $D_{{\partial }_j}{\partial }_k={\Gamma }_{jk}^i{\partial }_i$. For the Riemannian metric matrix $(g_{jk})$ and its inverse $(g^{jk})$,A curve $x(t)$ with affine parameter is a geodesic when $D_{\dot{x}}\dot{x}=0$, equivalentlyThe connection is metric compatible and torsion free; in particular ${\Gamma }_{jk}^i={\Gamma }_{kj}^i$. The coordinate velocities are ${\dot{x}}^j$, and the displayed coefficients are smooth functions of position.

For $v\in T_{p}M$ whose initial-value geodesic ${\gamma }_v$ exists through time $1$, with ${\gamma }_v(0)=p$ and ${\dot{\gamma }}_v(0)=v$, define ${\exp }_p(v)={\gamma }_v(1)$. The domain is open and contains a neighborhood of zero. Indeed the geodesic equation is a smooth first-order ODE in $(x,\dot{x})$; the constant solution with initial velocity zero exists through time $1$, and existence on a compact time interval and smooth dependence persist for nearby initial data. Reparameterization and uniqueness give ${\exp }_p(tv)={\gamma }_v(t)$ wherever defined. Consequently the differential of the exponential map at zero isThe inverse function theorem gives a neighborhood $V$ of zero on which ${\exp }_p$ is a diffeomorphism onto a neighborhood $U$ of $p$. Shrink $V$ to a ball in $T_{p}M$ and choose an orthonormal basis there. Its linear coordinate functions composed with ${\exp }_p^{-1}$ are the geodesic normal coordinates on $U$. This establishes well-defined local coordinates without assuming geodesic completeness.

In these geodesic normal coordinates, every radial geodesic has coordinate expression $x(t)=tv$. Its equation at $t=0$ reads ${\Gamma }_{jk}^i(p)v^j{v}^k=0$ for every $v$. The coefficients are symmetric in $j,k$ by the torsion-free connection property. Evaluating at the coordinate unit vectors and their pairwise sums, or using the polarization identity, therefore givesThis is the fact that Christoffel symbols vanish at the center of normal coordinates. Also $g_{ij}(p)={\delta }_{ij}$, since $(d{\exp }_p{)}_0$ is the identity in an orthonormal basis.

For sufficiently small $r>0$, the geodesic sphere is $S_r(p)={\exp }_p\{v:|v{|}_{g_p}=r\}$. It is a smooth hypersurface because ${\exp }_p$ is a diffeomorphism on the ball. It equals the local distance sphere $\{q:d_g(p,q)=r\}$ for sufficiently small radius, as the length comparison below shows; no assertion of global smoothness for large radii is needed.

The full Gauss lemma isfor $v$ in a star-shaped domain of ${\exp }_p$ and $w\in T_{p}M$. To prove it, set $F(t,s)={\exp }_p(t(v+sw))$, $T={\partial }_{t}F$ and $S={\partial }_{s}F$, for $s$ near zero. Every $t$-curve is a geodesic and has squared speed $|v+sw{|}_{g_p}^2$. The torsion-free connection property and the commuting parameter fields give $D_{t}S=D_{s}T$; metric compatibility and $D_{t}T=0$ then giveAt $t=0$, $F(0,s)=p$, so $S=0$. Integration from $0$ to $1$ followed by setting $s=0$ proves the boxed identity, because $T(1,0)=(d{\exp }_p{)}_{v}v$ and $S(1,0)=(d{\exp }_p{)}_{v}w$. For any $v$ in the domain, a small neighborhood of its compact radial segment suffices for the variation. The zero vector is covered directly.

In particular, if $w\perp v$, the images of the radial and spherical tangent directions are orthogonal; taking $w=v$ also proves preservation of radial length. Thus in geodesic polar coordinates the Riemannian metric has radial part $d{r}^2$ and no mixed radial-angular term. Every curve in a small normal ball from its center to radial coordinate $r$ has length at least $r$, by integrating the absolute radial derivative; the radial geodesic realizes this length. To exclude shortcuts leaving the ball, choose a larger normal ball of radius $\rho$ with compact closure in $U$ and restrict to $r<\rho$: an escaping curve first reaches radial coordinate $\rho$ and already has length at least $\rho$. This proves the asserted small-radius distance interpretation and orthogonality of radial geodesics to geodesic spheres.

In the two null coordinates the inverse metric has $g^{UV}=g^{VU}=-f^{-1}$ and $g^{UU}=g^{VV}=0$. ThereforeBoth tangents are null and future-directed, since $f>0$ and both coordinate directions are future-directed. To prove the affine parameter property, use the null-gradient affine-geodesic identity. The one-form $k_a=-{\nabla }_{a}U$ is closed, so ${\nabla }_a{k}_b={\nabla }_b{k}_a$. ConsequentlyThe same argument applies to ${\ell }_a=-{\nabla }_{a}V$. Thus both vector fields satisfy the geodesic equation without a nonaffinity term, and generate affinely parametrized null geodesic congruences. The coordinate vectors alone need not be affine; the factors $f^{-1}$ are essential.

In Newtonian gauge, the coordinate-comoving observer has normalized four-velocity $u^{\mu }=a^{-1}(1-\psi ,0)$, hence $u_0=a(1+\psi )$. In the $(+---)$ convention the measured photon energy is $E=u_{\mu }{p}^{\mu }$. Thus to first order

Use the time component of the geodesic equation and divide by $p^0=d\eta /d\lambda$:Let $\mathcal{H}=a^{\prime }/a$ and take the photon path direction $\mathbf{e}$ at zeroth order. The left side divided by $ϵ/a^2$ is ${ϵ}^{\prime }/ϵ-2\mathcal{H}+2\mathcal{H}\psi -d\psi /d\eta$. The contraction on the right side, using the provided connection coefficients, isThe homogeneous expansion cancels. Therefore the scalar photon energy-redshift law isThe total derivative is $d/d\eta ={\partial }_{\eta }+\mathbf{e}\cdot \nabla$ on the unperturbed ray. Perturbing that path inside a first-order potential would contribute only at second order. The homogeneous $a^{-1}$ energy redshift is already removed by the comoving momentum variable $ϵ$.

The product of the magnetic field with the velocity in this question is the three-dimensional cross product. Fix ${ϵ}_{123}=1$ and write $K^i{}_j={ϵ}^i{}_{kj}{B}^k$, so $K\mathbf{v}=\mathbf{B}\times \mathbf{v}$. In coordinates $y^0=t$, $y^i=x^i$ on ${\mathbb{R}}^3\times \mathbb{R}$, prescribe an affine connection by the following Christoffel symbols:with all other Christoffel symbols zero. These are smooth globally in the given Cartesian coordinates. The lower-index symmetry makes this a torsion-free connection. It is the geometrization of the Lorentz force by an affine connection; no condition on the curl of $\mathbf{E}$ or the divergence of $\mathbf{B}$ is required for this construction.

An affinely parametrized geodesic, with primes denoting differentiation with respect to $\lambda$, satisfiesOn the branch $t^{\prime }=c\ne 0$, use $t$ itself as an affine parameter. Dividing the spatial equation by $c^2$ givesConversely, each solution of this equation makes $\lambda =t$, $y(t)=(t,\mathbf{x}(t))$ an affinely parametrized geodesic of the constructed affine connection. Therefore its image is also an unparametrized geodesic. A general change of parameter adds a term proportional to the tangent in the geodesic equation, leaving the curve unchanged. The branch $t^{\prime }=0$ is not a trajectory with time as parameter.

For the metric realization, assume $\mathbf{E}=-\nabla U$ and $\nabla \cdot \mathbf{B}=0$. Introduce the differential formThen $dF=(\nabla \cdot \mathbf{B})d{x}^1\wedge d{x}^2\wedge d{x}^3=0$. The global Poincare lemma on the contractible space ${\mathbb{R}}^3$ supplies a one-form $A=A_{i}d{x}^i$ with $dA=F$, equivalently a magnetic vector potential with $\nabla \times \mathbf{A}=\mathbf{B}$. An explicit radial-gauge potential for a divergence-free magnetic field isThis is the radial homotopy formula for a closed two-form and is smooth even at the origin. For a constant magnetic field it gives $\mathbf{A}=\frac{1}{2}\mathbf{B}\times \mathbf{x}$.

Consider the Eisenhart-Duval lift with the Lorentzian metricThe independent one-forms $d{x}^1,d{x}^2,d{x}^3,dt,du-2A_{i}d{x}^i-Udt$ exhibit three positive directions and a two-dimensional block with one positive and one negative direction. Thus the metric is nondegenerate, of signature $(4,1)$. Its coefficients do not depend on $u$, and $g_{uu}=0$, so $\xi ={\partial }_u$ is a null Killing vector field. Indeed $g_{Au}$ are constant, so ${\Gamma }^A{}_{Bu}=0$ for the Levi-Civita connection and $\xi$ is parallel.

The geodesic Lagrangian for an affine parameter $\lambda$ isThe cyclic coordinate $u$ gives the conserved quantity $p_u=\partial L/\partial u^{\prime }=t^{\prime }$. Work at a nonzero value of $p_u$ and rescale the affine parameter to set $t^{\prime }=1$. This is the essential step in the null Kaluza-Klein reduction of a stationary force: one fixes the momentum along the null isometry and projects its geodesics, rather than dividing by $g_{uu}$.

The spatial Euler-Lagrange equations, before setting $t^{\prime }=1$, areSince $t^{\prime \prime }=0$ and $F_{ij}{v}^j=-(\mathbf{B}\times \mathbf{v}{)}_i$, their reduction isEquivalently, after taking $t$ as the affine parameter, the term $\dot{u}$ is a total derivative and the reduced Lagrangian is $L_{\mathrm{red}}=\frac{1}{2}|\dot{\mathbf{x}}{|}^2-2\mathbf{A}\cdot \dot{\mathbf{x}}-U$. Its Euler-Lagrange equations give the same sign and factor two.

There is also an explicit converse using null geodesics. For any physical trajectory, setThis makes its five-dimensional tangent null. The conserved momentum of the cyclic coordinate $t$ becomesThe quantity in parentheses is conserved, because its derivative is $\dot{\mathbf{x}}\cdot (2\mathbf{B}\times \dot{\mathbf{x}})=0$. Hence the $t$ equation, as well as the spatial and $u$ Euler-Lagrange equations, is satisfied. Every trajectory admits a null geodesic lift with $p_u=1$, and every such lift projects to the required trajectory.

Finally, a time-independent gauge transformation $\mathbf{A}\mapsto \mathbf{A}+\nabla \chi$ is absorbed by $u\mapsto u+2\chi$. The one-form $du-2A_{i}d{x}^i-Udt$ and the five-dimensional metric are unchanged. These gauge transformations of an Eisenhart-Duval lift alter the reduced Lagrangian only by the total derivative $-2d\chi /dt$.

Regard complex matrices as a real vector space of dimension $2n^2$ and Hermitian matrices as a real vector space of dimension $n^2$. Set $F(A)=A{A}^{\ast }$. Its derivative is $D{F}_A(H)=H{A}^{\ast }+A{H}^{\ast }$. At a unitary $A$, any Hermitian $S$ is reached by $H=SA/2$. Hence $I$ is a regular value, and the regular level set theorem makes $\mathrm{U}(n)=F^{-1}(I)$ a smooth real manifold of dimension $\boxed{n^2}$.

The tangent space is the kernel of this derivative. Equivalently,At $I$, therefore, $H$ is skew-Hermitian. The explicit geodesic with initial point $I$ and velocity $H$ is $\gamma (t)=e^{tH}$. The matrix exponential is unitary because $e^{t{H}^{\ast }}{e}^{tH}=I$, and $\dot{\gamma }(0)=H$. Its acceleration is $\gamma H^2$. Since $H^2$ is Hermitian and every tangent vector at $\gamma$ is $\gamma K$ with $K$ skew-Hermitian, the induced real inner product satisfiesIndeed the trace of a Hermitian matrix times a skew-Hermitian matrix is purely imaginary. Its acceleration is therefore always normal to the manifold, which is the geodesic equation for an induced Euclidean metric. Diagonalizing $H=W\mathrm{diag}(i{\lambda }_j)W^{\ast }$ gives the equivalent explicit construction $W\mathrm{diag}(e^{it{\lambda }_j})W^{\ast }$.

A smooth vector field along a map $\gamma :[a,b]\to M$ is a smooth section of the pullback tangent bundle ${\gamma }^{\ast }TM$: it assigns $V(t)\in T_{\gamma (t)}M$, with smooth coefficients in each local frame. It need not be the restriction of one ambient vector field, since a curve may return to the same point with different values of $V$.

An affine connection, also called a Koszul connection, induces a pullback connection on this bundle. In coordinates, write $V=V^i(t){\partial }_i$ and ${\nabla }_{{\partial }_j}{\partial }_k={\Gamma }^i{}_{jk}{\partial }_i$. Its covariant derivative along a curve isThe connection transformation law, or the Leibniz rule applied to a change of frame, makes this expression independent of the frame. In particular it is defined even when $\dot{\gamma }=0$; it differentiates the section's coefficients as well as the frame.

The field is parallel precisely when $D_{t}V=0$. In a frame along a coordinate segment this is the linear ordinary differential equationThe Picard-Lindelof theorem gives a unique solution with prescribed initial value. Smooth coefficients on compact subintervals are bounded, and the usual norm estimate for a linear ordinary differential equation prevents finite-time blowup there. A finite sequence of coordinate segments covers the image of the compact interval $[a,b]$; solve successively and use uniqueness on overlaps. This gives a unique parallel field on the whole interval for every $V_a\in T_{\gamma (a)}M$.

Define parallel transport by ${\tau }_t(V_a)=V(t)$. Linearity of the equation and uniqueness give ${\tau }_t(c{V}_a+d{W}_a)=c{\tau }_t(V_a)+d{\tau }_t(W_a)$. Solving along the reversed curve supplies its inverse. Therefore each ${\tau }_t:T_{\gamma (a)}M\to T_{\gamma (t)}M$ is a linear isomorphism.

For two affine connections, put $\Delta (X,Y)={\nabla }_{X}Y-{\tilde{\nabla }}_{X}Y$. Their rules givesince the two $X(f)Y$ terms cancel. The tensoriality of this operation proves that the difference of affine connections is a tensor, of type $(1,2)$, with smooth coordinate coefficients ${\Gamma }^i{}_{jk}-{\tilde{\Gamma }}^i{}_{jk}$.

A parametrized geodesic satisfies $D_t\dot{\gamma }=0$, or the geodesic equationWith $v^i={\dot{\gamma }}^i$, this becomes the first-order system ${\dot{\gamma }}^i=v^i$, ${\dot{v}}^i=-{\Gamma }^i{}_{jk}(\gamma )v^j{v}^k$. Its right side is smooth, so the Picard-Lindelof theorem gives a unique local solution for each initial point and tangent vector. This establishes uniqueness with the parametrization fixed.

The two accelerations differ byIf $\Delta (v,v)=0$ for every tangent vector, either acceleration vanishes exactly when the other does. Conversely, start the $\nabla$-geodesic with arbitrary initial vector $v$ at an arbitrary point. If it is also a $\tilde{\nabla }$-geodesic with the same parameter, evaluating this identity initially gives $\Delta (v,v)=0$. HencePolarization makes the latter condition equivalent to $\Delta (u,v)+\Delta (v,u)=0$: parametrized geodesics determine the symmetric part of an affine connection. In particular torsion-free connections are determined by their parametrized geodesics, since their difference is also symmetric and must then vanish.

For radial motion, take an affine parameter $\lambda$ and use dots for $d/d\lambda$. The radial null geodesics of the Vaidya metric satisfy the two radial geodesic equationsThe radial squared norm is $-f{\dot{v}}^2+2\dot{v}\dot{r}$. If $v$ is constant it vanishes, and the geodesic equation reduces to $\ddot{r}=0$. Thus constant-$v$ radial curves are null geodesics, with $r=a\lambda +b$ an affine parametrization. Future ingoing motion has $a<0$.

For the other family let $dr/dv=f/2$. Its tangent $k={\partial }_v+(f/2){\partial }_r$ is a null vector. Using the Christoffel symbols of the Vaidya metric, direct differentiation givesHence these are also null geodesics, but $v$ is generally not affine. For completeness, put $p=dv/d\lambda$ and chooseThis makes $\ddot{v}=-(M/r^2){\dot{v}}^2$. Differentiating $\dot{r}=(f/2)\dot{v}$ then gives $\ddot{r}=(f_v/2){\dot{v}}^2=-M^{\prime }{\dot{v}}^2/r$; substitution verifies the radial $r$ equation as well. Thus the null condition and both radial geodesic equations hold, not just the null condition alone.

Set $f=1-2M(v)/r$, so $f_r=2M/r^2$ and $f_v=-2M^{\prime }/r$, where $M^{\prime }=dM/dv$. Dots in this solution mean derivatives with respect to proper time $\tau$. The geodesic Lagrangian becomesThe Euler-Lagrange equation for $r$ gives $\ddot{v}+(f_r/2){\dot{v}}^2-r({\dot{\theta }}^2+{\sin }^2\theta {\dot{\varphi }}^2)=0$. The Euler-Lagrange equation for $v$ givesEliminate $\ddot{v}$ between these equations. The two angular Euler-Lagrange equations follow by differentiating $-2r^2\dot{\theta }$ and $-2r^2{\sin }^2\theta \dot{\varphi }$. Together the four geodesic equations areReading off the symmetric quadratic coefficients gives all nonzero Christoffel symbols of the Vaidya metric:Every unlisted Christoffel symbol vanishes. The lower-index symmetry is that of the Levi-Civita connection. In particular, the time-dependent mass contributes $+M^{\prime }/r$ to ${\Gamma }^r{}_{vv}$; it is absent from ${\Gamma }^v{}_{vv}$. There are eleven independent nonzero entries, or fifteen when both lower-index orders are counted.

The radial null geodesics satisfyThe ingoing family has $v=\text{constant}$, with $r$ decreasing toward the future; its tangent $-{\partial }_r$ is affinely parametrized because ${\Gamma }^a{}_{rr}=0$. The outgoing family hasAway from the Schwarzschild event horizon, integration gives $v=2r_{\ast }+\text{constant}$. These are null geodesics up to reparametrization: a null direction in the two-dimensional radial geometry is automatically pregeodesic, and the angular geodesic equations are satisfied by constant angles. On the Schwarzschild event horizon the outgoing family instead has $r=2M$, with tangent proportional to ${\partial }_v$.

In a Finkelstein diagram, plot $T=v-r$ vertically and $r$ horizontally. Ingoing rays obey $T+r=\text{constant}$; outgoing rays obeyOutside the Schwarzschild event horizon, outgoing rays increase $r$; on it they remain at $r=2M$; inside it they decrease $r$ even though $v$ increases. Thus both radial future null directions point toward smaller $r$ inside the black hole. Every future timelike direction lies between these null directions, so it also moves toward smaller $r$ there.

Arrows show future propagation. The vertical red line is the Schwarzschild event horizon; the black boundary at $r=0$ is the Schwarzschild singularity.

For radial free fall, let $E=f{\dot{t}}_A>0$ be Alice's conserved Killing energy. The inward branch of the timelike geodesic has $\dot{r}=-\sqrt{E^2-f}$, henceWith $\delta =r-2M$, the near-horizon geodesic equations giveThe constants in the logarithms are understood to make their arguments dimensionless. Therefore $u_A=t_A-r_{\ast }=2t_A+O(1)$. For an outgoing signal reaching the fixed-radius Bob, $t_B=u_A+r_{\ast }(r_B)$, so $t_B=2t_A+O(1)$. The redshift consequently behaves asThis is the answer when $t$ is Bob's Schwarzschild time, as appropriate to the stated observation. It is the Schwarzschild surface gravity $a=\kappa =1/(4M)$. If $t$ instead means the emission Schwarzschild time $t_A$, the same redshift is proportional to $e^{t_A/(2M)}$ and gives $M=1/(2a)$. The two answers use different clocks, not different dynamics. If the measured exponential uses Bob's proper time ${\tau }_B$, its rate is $1/(4M\sqrt{f_B})$. In SI units the reception-time result is $M=c^3/(4Ga)$.

The reception-time exponent also holds for smooth radial infall crossing the future Schwarzschild event horizon with finite nonzero $dr/d\tau$: regular Ingoing Eddington-Finkelstein coordinates give finite $v$ there, $u=v-2r_{\ast }\sim -4M\log ({\tau }_H-\tau )$. Thus no special value of $E$ is needed.

Let $q^i=e^i-\frac{1}{2}{h}^i{}_j{e}^j$, so $p^0=ϵ/a^2$ and $p^i=(ϵ/a^2)q^i$. The orthonormal tetrad gives ${\delta }_{ij}{e}^i{e}^j=1$. To first order,Hence the time-index Christoffel symbols contract toThe $\mathcal{H}{h}_{ij}$ term cancels the perturbation in the coordinate direction norm. This cancellation is why the local tetrad direction must not be mistaken for the coordinate direction.

Along the photon null geodesic, $d/d\lambda =p^{0}d/d\eta$, andThe time component of the geodesic equation now gives the photon energy redshift from a tensor metric perturbation:The derivative on the left follows the ray; the dot on $h$ is a partial conformal time derivative. The background redshift has already been removed by using comoving energy. The calculation only needs ${\Gamma }^0{}_{00}$ and ${\Gamma }^0{}_{ij}$; the unrelated printed ${\Gamma }^i{}_{j0}$ should contain $\frac{1}{2}{\dot{h}}^i{}_j$, rather than the full derivative.

Use spatial indices $i,j=1,2,3$, set ${\varphi }_i={\partial }_i\varphi$, and assume $\varphi$ is smooth. The inverse metric has $g^{tt}=-c^{-2}{e}^{-2\varphi /c^2}$ and $g^{ij}={\delta }^{ij}$. The Christoffel symbols of the Levi-Civita connection follow from ${\Gamma }^{\alpha }{}_{\beta \gamma }=g^{\alpha \delta }({\partial }_{\beta }{g}_{\gamma \delta }+{\partial }_{\gamma }{g}_{\beta \delta }-{\partial }_{\delta }{g}_{\beta \gamma })/2$:All other components vanish. Although $g_{tt}=-c^2-2\varphi +O(c^{-2})$ diverges, its affine connection has a smooth limit on compact subsets:This is the Newtonian connection from an exponential lapse. It is torsion-free, and its geodesic equation, using $t$ as an affine parameter when $\dot{t}\ne 0$, is $d^2{x}^i/d{t}^2=-{\partial }_i\varphi$. Thus $\varphi$ has the role of a Newtonian gravitational potential. There is no finite limiting nondegenerate Lorentzian metric in these fixed coordinates; $-g(c)/c^2$ tends to $d{t}^2$, a rank-one tensor. The limiting connection preserves $dt$ and the contravariant spatial tensor with components $h^{ij}={\delta }^{ij}$, $h^{t\mu }=0$, which describe the degenerate temporal/spatial structures of this limit.

Use $x^0=t$ and the Levi-Civita symbol ${ϵ}^1{}_{23}=1$. Define a smooth affine connection byThe lower indices are symmetric, so it is a torsion-free connection. With $t$ itself as parameter, the time component of the geodesic equation is identically satisfied, while its spatial components becomeThe last term is $2(\mathbf{B}\times \dot{\mathbf{x}}{)}^i$, with the wedge in the source interpreted as the three-dimensional cross product. Thus the spacetime lifts $(t,\mathbf{x}(t))$ are actually affinely parametrized geodesics, which proves the requested unparametrized claim. Conversely a geodesic with nonzero constant $dt/ds$ can be parametrized by $t$ and obeys this particle equation; purely spatial geodesics with $dt/ds=0$ are not such trajectories. No field equations or metric compatibility are needed. This construction is the affine connection for a velocity-linear force.

Parallel transport the auxiliary null vector $N$ and a screen basis along each generator of the null geodesic congruence. This keeps the screen-space projector fixed under the corresponding transported derivative. Use the curvature convention $[{\nabla }_c,{\nabla }_b]V_a=-R^d{}_{acb}{V}_d$.

With $B_{ab}={\nabla }_b{U}_a$, differentiation of the affine geodesic equation givesHere the first term comes from commuting the derivative past $U^c$, and the second from the curvature commutator. The constraints $B_{ab}{U}^a=B_{ab}{U}^b=0$ follow respectively from the null normalization and affine geodesic equation. Projecting onto the screen therefore yields the optical matrix equationThe trace of the tidal matrix is $R_{ab}{U}^a{U}^b$; extra terms from replacing the screen trace by the spacetime trace vanish by curvature antisymmetry. Taking the trace givesDecompose the optical tensor as $\hat{B}=\frac{1}{2}\theta I+\hat{\sigma }+\hat{\omega }$. The null shear is symmetric and trace-free and the null twist is antisymmetric. Their cross traces vanish, andThe negative sign in the last term is the identity $\mathrm{tr}({\hat{\omega }}^2)=-{\sum }_{I,J}{\hat{\omega }}_{IJ}^2$ on the positive-definite screen. Substitution proves the four-dimensional Null Raychaudhuri equation:The coefficient is $1/(D-2)$ in $D$ spacetime dimensions; the factor $1/2$ here belongs to the four-dimensional congruence in this question, not the five-dimensional geometry of the next question.