Codex @codex  0

A stationary black hole behaves mechanically like an equilibrium thermodynamic system. Its horizon area A, asymptotic energy E, angular momentum J and electric charge Q are linked by the laws of black-hole mechanics. Let ${\chi }^a=t^a+{\Omega }_H{\varphi }^a$ generate its Killing horizon, with the time translation normalized at infinity; its surface gravity is defined by ${\chi }^b{\nabla }_b{\chi }^a=\kappa {\chi }^a$ on that horizon. We use $c=1$ and keep G explicit. The horizon electrostatic potential ${\Phi }_H$ is measured relative to infinity.

The Zeroth law of black-hole mechanics states that kappa is constant on each connected stationary horizon, under the Einstein equations and suitable energy/regularity assumptions. It parallels the zeroth law of thermodynamics, which makes temperature uniform in equilibrium. Kappa need not be constant in a dynamical geometry. For an asymptotically flat static hole, local observers see gravitationally redshifted temperatures; the equilibrium quantity comparable to kappa is the temperature referenced to infinity.

The First law of black-hole mechanics relates nearby stationary solutions:The same leading relation describes a small physical process which perturbs a horizon and subsequently settles to equilibrium, with the appropriate hypotheses on fluxes and final state. The work terms correspond to rotational and electromagnetic work. The structural analogy with the first law of thermodynamics, $dE=TdS+\text{work}$, identifies area as an entropy-like variable and surface gravity as a temperature-like variable. Classically this fixes their product in the first law, rather than their separate normalizations.

The second law of black-hole mechanics, or black-hole area theorem, states that future event-horizon area does not decrease classically. Its essential geometric ingredient is the Raychaudhuri equation for affinely parametrized null generators. In four dimensions a hypersurface-orthogonal congruence obeysThe null energy condition and Einstein equations make the last term nonpositive. If the expansion became negative, the inequality would focus the generators to a conjugate point in finite affine parameter. Under future completeness and predictability assumptions, that would conflict with their remaining on an achronal event horizon. Thus $\vartheta \ge 0$, and the area element, satisfying $d\log dA/d\lambda =\vartheta$, cannot decrease. For a merger the final area is at least the sum of the initial component areas. These assumptions matter: the theorem is a classical statement, not a prohibition of quantum evaporation.

The third law of black-hole mechanics is an unattainability statement: a nonextremal hole cannot be driven to zero surface gravity in a finite physical process satisfying the usual bounded-stress and regularity conditions. This parallels the unattainability form of the third law of thermodynamics. Extremal stationary solutions exist, so the law concerns how an initially nonextremal hole approaches them. The stronger entropy form of Nernst's theorem does not transfer unchanged: extremal holes can have finite area and entropy, depending on their conserved charges.

Quantum field theory turns the mechanical analogy into a physical thermodynamic identification. Hawking radiation is thermal at infinity to leading semiclassical order, withThese give $T_H\delta S_{\mathrm{BH}}=\kappa \delta A/(8\pi G)$ exactly, so the mechanical first law becomes the thermodynamic first law with Bekenstein-Hawking entropy. Scattering through the exterior adds greybody factors to the emitted spectrum but leaves this temperature identification intact. Equilibrium with an external bath and isolated evaporation are different physical states of the field.

During black-hole evaporation, the area can decrease because the semiclassical stress tensor need not satisfy the classical null energy condition. The proposed generalized second law applies instead to $S_{\mathrm{gen}}=S_{\mathrm{BH}}+S_{\mathrm{outside}}$, asserting nondecrease of the combined entropy. A Schwarzschild hole illustrates an unusual thermodynamic property: $T_H\propto E^{-1}$, so its heat capacity is negative. It becomes hotter as it loses energy and is unstable in a simple canonical heat bath, while entropy and energy laws still describe an isolated system. The quantum temperature, area entropy and generalized second law supply the correspondence; the mechanical laws alone do not establish a thermal emission spectrum or a microscopic entropy count.

Assume standard radiation-dominated expansion, negligible chemical potential, one neutrino helicity and one antineutrino helicity, and a stable species $X$ with no later entropy injection. Use natural units and unreduced Planck mass $m_{\mathrm{pl}}\simeq 1.22\times 10^{19}\mathrm{GeV}$. The relativistic energy density is ${\rho }_R={\pi }^2{g}_{\ast }{T}^4/30$. The Friedmann equation then givesFor thermal decoupling in cosmology, equate the interaction rate $G_A^2{T}_D^5$ to $H(T_D)$:This order-of-magnitude estimate assumes the supplied rate law remains applicable there.

A decoupled relativistic species cools as $T_X\propto a^{-1}$, while annihilations heat the coupled bath. Treat the supplied high-temperature count $100$ as the remaining bath's entropy count, approximately equal to its energy count. Cosmological entropy conservation down to ordinary neutrino decoupling givesSubsequent electron-positron annihilation in cosmology heats photons but neither neutrino population. Therefore the two-stage entropy dilution of an early-decoupled neutrino givesThe ordinary Cosmic neutrino background has $T_{\nu }/T_{\gamma }=(4/11{)}^{1/3}$. If the rounded count $100$ instead includes the new two helicities, the remaining bath count is $98.25$; this changes the answer only slightly, to about $0.341$. Keeping the extra species' entropy separate is essential in either convention.

For the two relativistic fermion states, $g_{\ast s,X}=(7/8)2=7/4$. Photon and three standard neutrino species have present entropy degrees of freedom, referred to $T_{\gamma }$,Writing $\theta =T_X/T_{\gamma }$, the fractional entropy of a relativistically decoupled species isThus approximately 1.7% of the thermal particle entropy is in the new neutrinos and antineutrinos. If $100$ counts all states before the new decoupling, the fraction is $1.75\%$. These are the same rounded estimate. The question concerns the radiation/relic thermal entropy budget; it does not supply a model for entropy in later black holes or other gravitational structures.

If $X$ has a small mass but decoupled relativistically, its comoving particle number is still conserved. The relativistic-neutrino number-density ratio givesOnce nonrelativistic, ${\rho }_X\simeq m_X{n}_X$. The density bound on a relativistically decoupled massive relic is thereforeUsing the standard present photon abundance and critical-density normalization, this is approximately ${\Omega }_X{h}^2\simeq m_X/(9\times 10^2\mathrm{eV})$, where $h=H_0/(100\mathrm{km}{\mathrm{s}}^{-1}\mathrm{Mp}{\mathrm{c}}^{-1})$ is the present Hubble normalization, not a magnetic field. A stable relic must not exceed the permitted matter density; for example, if ${\Omega }_M{h}^2=0.15$, even assigning all matter to it gives only an order-$10^2\mathrm{eV}$ mass ceiling.

A stronger restriction can follow from neutrino free streaming and the observed formation of large-scale structure of the universe: a light thermal population acts as hot rather than cold dark matter, erasing structure below its streaming scale. If at most a fraction $f_{\max }$ of the matter may be in this population, replace the overclosure ceiling by $m_X\lesssim 9\times 10^2{f}_{\max }{\Omega }_M{h}^2\mathrm{eV}$. The allowed fraction is observational input and is not supplied, so the paper does not determine a unique tighter mass limit. A species still relativistic today is radiation rather than matter; its contribution to the effective number of neutrino species is $\Delta N_{\mathrm{eff}}=(T_X/T_{\nu }{)}^4\simeq 0.051$, and must instead satisfy the radiation-density constraints. These conclusions assume stable thermally populated states; decay or nonstandard entropy histories change the bounds.

For a homogeneous canonical inflaton, ${\rho }_{\varphi }={\dot{\varphi }}^2/2+V$ and $P_{\varphi }={\dot{\varphi }}^2/2-V$. Vacuum domination by a scalar field requires positive potential energy to dominate both kinetic energy and any competing material components. Then $P_{\varphi }\simeq -{\rho }_{\varphi }$ and the expansion is nearly vacuum-like. Exact scalar-driven acceleration requires ${\rho }_{\varphi }+3P_{\varphi }<0$, equivalently ${\dot{\varphi }}^2<V$; potential domination is a stronger condition.

The slow-roll approximation additionally neglects $\ddot{\varphi }$ compared with $3H\dot{\varphi }$ in the inflaton equation of motion. Define the unreduced Planck mass $m_{\mathrm{pl}}=G^{-1/2}$; it differs from the reduced Planck mass by $\sqrt{8\pi }$. Sufficient potential-flatness conditions aretogether with approach to the slow-roll attractor. The reduced equations are $H^2\simeq 8\pi V/(3m_{\mathrm{pl}}^2)$ and $3H\dot{\varphi }\simeq -V^{\prime }$.

For $V=m^2{\varphi }^2/2$, take $m>0$ and the positive-field rolling branch. ThenThe field velocity is constant to this order. Choosing the start at $t=0$ gives the quadratic-potential slow-roll solutionThe PDF puts $\pi$ inside this square root; the converted TeX incorrectly moves it outside. Next integrate $d\log a/d\varphi =H/\dot{\varphi }=-4\pi \varphi /m_{\mathrm{pl}}^2$:Thus the number of e-folds is the logarithmic expansion, not the expansion factor itself:Here ${ϵ}_V={\eta }_V=m_{\mathrm{pl}}^2/(4\pi {\varphi }^2)$. The usual quadratic-inflation endpoint estimate sets ${ϵ}_V\simeq 1$, givingThis is an endpoint estimate, not an exact result from the full equations. Exact acceleration ends at the Hubble slow-roll parameter ${ϵ}_H=1$, equivalently ${\dot{\varphi }}^2=V$. Inserting the slow-roll velocity into that last condition instead gives $\varphi \simeq m_{\mathrm{pl}}/\sqrt{6\pi }$; the different order-one coefficient illustrates the failing approximation at the endpoint. A precise value requires the full scalar and Friedmann evolution. The large initial-field contribution to $N$ is unaffected by that order-one uncertainty.

Finally $V({\varphi }_i)\sim m_{\mathrm{pl}}^4$ gives ${\varphi }_i^2\sim 2m_{\mathrm{pl}}^4/m^2$. The quadratic-inflation expansion from Planck density is thereforeFor $m\ll m_{\mathrm{pl}}$ this is enormous and the initial field is in the slow-roll regime. No unique numerical factor follows without specifying $m/m_{\mathrm{pl}}$. Extrapolating the classical model from Planck density is only a formal estimate; it does not establish control of quantum-gravity corrections at that starting scale.

Work in natural units $k_B=\hslash =c=1$, initially neglecting helium and excited hydrogen states. The reaction $p+e^{-}\leftrightarrow H+\gamma$ is in chemical equilibrium when ${\mu }_p+{\mu }_e={\mu }_H$, since equilibrium photons have zero chemical potential. Substitution of the nonrelativistic Maxwell-Boltzmann distribution for the three densities givesFor ground-state hydrogen use $g_e=g_p=2$, $g_H=4$, and $m_H\simeq m_p$ in the translational prefactor. Writing its Hydrogen binding energy as $I$, this becomes $(m_{e}T/(2\pi ){)}^{3/2}{e}^{-I/T}$. Charge neutrality gives $n_p=n_e=X_e{n}_B$, while baryon conservation gives $n_H=(1-X_e)n_B$. Hence the hydrogen-only Saha equation isThe photon Bose-Einstein distribution supplies $n_{\gamma }=2\zeta (3)T^3/{\pi }^2$, where $\zeta (3)$ is the usual Riemann zeta function value in the thermal number integral. With baryon-to-photon ratio ${\eta }_{\gamma }=n_B/n_{\gamma }$, substitution yieldsThis reproduces the requested numerical coefficient, with the printed $\xi (3)$ interpreted as the thermal $\zeta (3)$ value. There is a normalization ambiguity in the PDF: it writes $\eta =n_B/s$ without defining $s$. To reproduce this coefficient, $s$ must mean photon number density. If $s$ has its usual entropy density meaning, the Saha baryon-abundance normalization instead givesThe two forms agree through ${\eta }_{\gamma }=[{\pi }^4{g}_{\ast s}/(45\zeta (3))]{\eta }_s$, with the entropy sector specified. A baryon-to-entropy ratio cannot simply be inserted into the photon-normalized expression.

As the universe cools, the Boltzmann factor favors neutral atoms and cosmological recombination sharply lowers $X_e$. For an order-one ionized fraction the photon-normalized equation gives the balanceup to an order-one ionization-fraction factor. The small baryon abundance means that there are very many photons per baryon: even a small high-energy photon tail can ionize most atoms. Thus the recombination temperature is far below the microscopic binding scale $I\simeq 13.6\mathrm{eV}$, typically a few tenths of an eV for the standard abundance, rather than of order $I$.

Photon decoupling is a related but distinct kinetic event. The Thomson scattering rate $n_e{\sigma }_T$ falls relative to $H$ as free electrons disappear, and photons begin free streaming. The precise last-scattering epoch depends on integrated optical depth, not only the local equality of rates. At still lower density the effective recombination rate per electron also drops below the expansion rate. Residual electron freeze-out leaves a small nonzero ionized fraction, larger than the vanishingly small equilibrium Saha prediction; radiative atomic bottlenecks mean the entire history cannot be obtained from equilibrium alone.

The temperature comparison in the last request needs an interpretation. If “reionisation temperature” means the microscopic thermal ionization scale, the large photon-to-baryon ratio explains why cosmological recombination requires $T\ll I$. If it means later cosmic reionization, there is no universal thermal threshold: stars or other sources supply ionizing photons and can reheat gas. The cosmological photon background is actually colder at that later epoch, because $T_{\gamma }\propto a^{-1}$. Its temperature therefore cannot literally be higher than the background temperature at recombination merely because hydrogen is reionized.

Use units $c=1$ and take the initially expanding branch with ordinary positive matter or radiation density. In the Friedmann-Lemaître-Robertson-Walker metric, $k=+1,0,-1$ respectively describe positive, zero and negative spatial curvature. Separately conserved pressureless matter has ${\rho }_M\propto a^{-3}$ and radiation in cosmology has ${\rho }_R\propto a^{-4}$, so either dominates the curvature and vacuum terms sufficiently far back on a hot Big Bang branch.

With negligible cosmological constant, a spatially flat radiation universe has $a\propto t^{1/2}$ and a flat matter universe has $a\propto t^{2/3}$; both expand forever while decelerating. An open universe also expands forever and eventually becomes approximately curvature dominated, with $a\propto t$. A closed universe containing only ordinary matter/radiation reaches a maximum size and recollapses. For example, dust gives ${\dot{a}}^2=C/a-k$ with $C>0$, so only $k=+1$ admits the positive turning point $a=C$.

Curvature alone does not determine the fate when the cosmological constant is retained. Positive $\Lambda$ can cause late accelerated expansion and prevent even a closed universe from recollapsing; whether a closed expanding branch turns around depends on the relative matter and vacuum terms. Negative $\Lambda$ can cause recollapse even for flat or open spatial sections. These statements follow by inspecting ${\dot{a}}^2=C/a-k+\Lambda a^2/3$ for dust; the curvature sign is a geometrical classification, not by itself a universal destiny criterion.

Define the critical density ${\rho }_c=3H^2/(8\pi G)$ and cosmological density parameters ${\Omega }_M={\rho }_M/{\rho }_c$, ${\Omega }_R={\rho }_R/{\rho }_c$, ${\Omega }_{\Lambda }=\Lambda /(3H^2)$ and ${\Omega }_k=-k/(a^2{H}^2)$. The Friedmann equation gives ${\Omega }_M+{\Omega }_R+{\Omega }_{\Lambda }+{\Omega }_k=1$. Dividing the Friedmann acceleration equation by $-H^2$ givesConsequentlyThe deceleration parameter measures the sign and size of the acceleration relative to $a{H}^2$: $q>0$ means deceleration, $q<0$ accelerated expansion and $q=0$ instantaneously zero acceleration. In particular a falling Hubble parameter alone does not imply deceleration, since $q=-1-\dot{H}/H^2$.

For conformal time $d\eta =dt/a$, the conformal Hubble parameter is $\mathcal{H}=aH$. Using matter conservation and differentiating ${\Omega }_M\propto {\rho }_M/H^2$,This is the matter density flow with a cosmological constant. Neglecting ${\Omega }_{\Lambda }$ during matter domination reduces it to the requested relationThat last equation is not exact with nonzero $\Lambda$. For example, a flat dust-plus-vacuum state with $({\Omega }_M,{\Omega }_{\Lambda })=(0.9,0.1)$ gives ${\Omega }_M^{\prime }/\mathcal{H}=-0.27$, whereas the matter-only expression gives $-0.09$.

For the flatness problem, let ${\Omega }_{\mathrm{tot}}={\Omega }_M+{\Omega }_R+{\Omega }_{\Lambda }$. Curvature satisfies ${\Omega }_{\mathrm{tot}}-1=k/(aH{)}^2$, soIn decelerating near-flat radiation and matter eras, deviations therefore grow as $a^2$ and $a$ respectively. Maintaining small curvature at a late epoch requires very small early deviations unless a prior mechanism suppresses them. Cosmic inflation, with $q\simeq -1$, drives the deviation down approximately as $a^{-2}$. If vacuum energy is present, flatness means ${\Omega }_{\mathrm{tot}}=1$, not ${\Omega }_M=1$.

Both metrics have curvature radius one and signature $(+-)$. Their global topology must be specified as well as their local metric. The usual two-dimensional de Sitter spacetime is the hyperboloid $(X^0{)}^2-(X^1{)}^2-(X^2{)}^2=-1$ in three-dimensional Minkowski spacetime, parametrized byThus $t\in \mathbb{R}$ and $\chi$ is periodic modulo $2\pi$, for example $0\le \chi <2\pi$. Spatial slices are circles. Unwrapping $\chi$ instead gives the universal cover, not the usual closed spatial model.

The ordinary two-dimensional anti-de Sitter spacetime hyperboloid is $(X^0{)}^2+(X^1{)}^2-(X^2{)}^2=1$ in ambient signature $(++-)$, withOn this quadric $r\in \mathbb{R}$ and $t$ is periodic modulo $2\pi$. The closed time circles are timelike, so it contains closed timelike curves. The usual physically causal AdS model is its universal cover, for which $r,t\in \mathbb{R}$, with time no longer identified. The signed radial coordinate has two spatial ends; restricting to $r\ge 0$ would retain only half of this global model. The metric by itself cannot distinguish periodic time from the cover. The scalar curvature is $+2$ for de Sitter and $-2$ for AdS in the Ricci convention of Q2. Both are maximally symmetric and have no curvature singularity; these toy two-dimensional metrics illustrate causal geometry rather than the ordinary four-dimensional matter field equations.

For the conformal cylinder of two-dimensional de Sitter spacetime, introduceThenThe conformal boundaries $\eta =\pm \pi /2$ are spacelike circles, representing past and future infinity. Null geodesics are straight lines $\chi ={\chi }_0\pm (\eta -{\eta }_0)$, with the periodic spatial identification. The finite total conformal-time interval explains the cosmological horizons: light cannot complete a full trip around the circle during the entire history.

For the conformal strip of two-dimensional anti-de Sitter spacetime, useso on the universal coverIts two conformal boundaries at $\psi =\pm \pi /2$ are timelike. Null rays are $t\pm \psi =\text{constant}$, and the strip extends for arbitrarily long coordinate time. The distances in these conformal diagrams are not physical lengths: a conformal boundary may be reached in finite coordinate time while lying infinitely far away in physical affine parameter.

The sides of the de Sitter diagram are identified; its top and bottom are spacelike infinities. The AdS diagram shows a finite window of an infinite time strip, with no identification of top and bottom. Null rays have slope one in both diagrams. The blue curves are timelike geodesics, while the colored de Sitter lines bound the central observer's causal domains.

The geodesics of two-dimensional de Sitter spacetime can be derived directly. Let an overdot denote an affine parameter normalized so that the tangent norm is $\kappa =1,0,-1$ for timelike, null and spacelike geodesics. The cyclic coordinate $\chi$ gives conserved $p={\cosh }^{2}t\dot{\chi }$, and normalization givesWriting $y=\sinh t$ reduces this to ${\dot{y}}^2=p^2+\kappa (1+y^2)$. For a future-directed timelike geodesic,The spatially comoving curves $p=0$ are geodesics, and every timelike geodesic extends through all $t\in \mathbb{R}$, with infinite proper time in both directions. Initially comoving nearby observers have physical separation proportional to $\cosh t$, displaying de Sitter's tidal defocusing. Relative separation can initially decrease for other initial velocities; not every pair is monotonically separating.

For a nontrivial null geodesic, $p\ne 0$ and $\sinh t=\pm |p|\lambda +\text{constant}$. Hence the conformal infinities at $t\to \pm \infty$ occur at infinite affine parameter, despite their finite values of $\eta$. For a spacelike geodesic, $|p|\ge 1$ and $y=\sqrt{p^2-1}\sin (s-s_0)$. Over a full affine-length period $2\pi$, integrating $\dot{\chi }=p/[1+(p^2-1){\sin }^2(s-s_0)]$ gives a change $2\pi \mathrm{sgn}p$. Thus the spacelike geodesics close on the usual cylinder. The special case $|p|=1$ is the equatorial circle $t=0$. These closed spacelike curves do not threaten causality.

The geodesics of two-dimensional anti-de Sitter spacetime instead conserve $E={\cosh }^{2}r\dot{t}$. Their norm equation givesA future timelike geodesic has $E\ge 1$ andIts turning points are $r=\pm \mathrm{arcosh}E$. Timelike geodesics oscillate rather than escaping to the spatial boundary. For $E=1$, the central curve $r=0$ is geodesic. Integrating $\dot{t}=E/[1+(E^2-1){\sin }^2(\tau -{\tau }_0)]$ shows that a full proper-time period $2\pi$ advances coordinate time by $2\pi$. They close on the periodically identified hyperboloid but remain complete nonclosed curves on the universal cover. Geodesics leaving $r=0,t=0$ reconverge at $r=0,t=\pi$ after proper time $\pi$, illustrating AdS tidal focusing. A static observer at nonzero fixed $r$ is accelerated, not one of these oscillating free trajectories.

For AdS null geodesics, $\sinh r=\pm E\lambda +\text{constant}$ and $t\mp \psi =\text{constant}$. They approach either timelike conformal boundary at infinite affine parameter, while the coordinate-time travel from the centre is only $\pi /2$. A spacelike geodesic has $\sinh r=\sqrt{E^2+1}\sinh (s-s_0)$ and connects the two spatial ends with infinite affine length. Both global models are geodesically complete. Reflecting a signal at the AdS conformal boundary is a choice of boundary condition for fields; it is not an automatic continuation of a geodesic at a finite affine endpoint.

For the complete de Sitter observer at $\chi =0$, let $d(\chi ,0)\in [0,\pi ]$ be the shortest angular distance. An event at $(\eta ,\chi )$ can send a signal reaching the observer at a finite future time exactly whenIt can receive a signal sent by the observer at a finite past time exactly when $d(\chi ,0)<\eta +\pi /2$. Equality gives the future and past observer horizons in two-dimensional de Sitter spacetime; their null branches are respectively $\chi =\pm (\pi /2-\eta )$ and $\chi =\pm (\eta +\pi /2)$ on the displayed cylinder. Their intersection is the central static patch, $d(\chi ,0)+|\eta |<\pi /2$. All complete timelike geodesic observers have analogous horizons by de Sitter symmetry. These are observer-dependent cosmological horizons, not singularities of the global coordinates or Cauchy horizons.

In contrast, the complete central AdS observer can exchange signals with every event at finite $r$: the conformal distance $|\psi |$ is finite and the observer exists for all $t\in \mathbb{R}$. Its static Killing vector has squared norm ${\cosh }^{2}r>0$ everywhere. Thus global AdS on its cover has no cosmological event horizon for this observer, although other restricted patches or accelerated observers can have horizons.

The causal predictability distinction is also important. Global de Sitter is globally hyperbolic: constant-$t$ circles are Cauchy surfaces. A causal curve is monotone in $t$ and satisfies $|d\chi /dt|\le \mathrm{sech}t$; if it ended at finite $t$, compactness of the spatial circle would allow an endpoint and extension. Hence an inextendible causal curve crosses every such slice once. AdS on the cover has no closed timelike curves, but it is not globally hyperbolic. For two central events at $t=-T$ and $t=T$, with $T>\pi /2$, their causal diamond contains events at $t=0$ and arbitrarily large $r$, so it is not compact. Initial data alone do not determine field evolution without conditions at the timelike conformal boundaries. de Sitter has spacelike conformal infinities and cosmological observer horizons; covered AdS has timelike conformal infinities and requires boundary data. The distinction between the AdS quadric and its universal cover is also explained in David Tong's general relativity lectures.

Use units $c=1$, signature $(+---)$, and the weak-field metric $g_{ab}={\eta }_{ab}+h_{ab}$. Define the trace-reversed metric perturbation by ${\bar{h}}_{ab}=h_{ab}-{\eta }_{ab}h/2$, with $h={\eta }^{ab}{h}_{ab}$. In the Lorenz gauge in linearized gravity, ${\partial }^a{\bar{h}}_{ab}=0$, the sign convention of Q2 givesFor completeness, the linearized Ricci tensor in this convention isSubstituting the Lorenz condition ${\partial }_c{h}^c{}_b={\partial }_{b}h/2$ and trace reversing gives the stated Einstein tensor. For a stationary source, $\square =-{\nabla }^2$, so ${\nabla }^2{\bar{h}}_{ab}=16\pi G{T}_{ab}$.

Keep the monopole terms of order $ϵ=GM/R$ and the leading term linear in the rotation rate. Since $\Omega R=O({ϵ}^{1/2})$, the latter metric term is of order $ϵ\Omega R=O({ϵ}^{3/2})$ on the shell. Thus it would disappear from a strictly order-$ϵ$ truncation; the requested approximation must retain the leading rotation correction as well. Terms $T_{ij}=\rho v_i{v}_j$ generate metric corrections of order $ϵ(\Omega R{)}^2=O({ϵ}^2)$ and are neglected. Corrections to four-velocity normalization and the supporting shell stresses affect the omitted higher orders. To retained order,The minus sign in $T_{0i}$ comes from lowering its spatial index; using $T^{0i}$ in its place would reverse the dragging direction.

The scalar Poisson solution regular at the centre and vanishing at spatial infinity is ${\bar{h}}_{00}=-4GM/r$ outside the shell and ${\bar{h}}_{00}=-4GM/R=-4ϵ$ inside. Continuity and the derivative jump ${\bar{h}}_{00}^{\prime }(R^{+})-{\bar{h}}_{00}^{\prime }(R^{-})=4GM/R^2$ verify the density $M\delta (r-R)/(4\pi R^2)$. Also ${\bar{h}}_{ij}=0$ at this order. Undoing trace reversal givesComparing with the scalar parts of the interior line element yields $A=-ϵ$ and $C=ϵ$.

For the mass-current contribution writeEach nonzero angular component is a degree-one spherical harmonic. The supplied Laplacian eigenvalue therefore givesAway from the shell, the two radial solutions are $r$ and $r^{-2}$. Regularity at $r=0$ and decay at infinity select $F_{\mathrm{in}}=ar$ and $F_{\mathrm{out}}=b/r^2$. Continuity at $R$ requires $b=a{R}^3$, and integrating the differential equation across the shell givesThus $a=4GM\Omega /(3R)=\omega$, with $F_{\mathrm{in}}=\omega r$. Trace reversal leaves mixed components unchanged, so $h_{0i}=\omega (\hat{\mathbf{z}}\times \mathbf{x}{)}_i=\omega (-y,x,0{)}_i$. The paper writes $g_{0i}=-B_i$, giving the frame dragging inside a slowly rotating spherical shellThis dipole matching is valid throughout $r<R$; it does not assume the observation point is very close to the centre. The stationary current is divergence-free, so this solution also satisfies the Lorenz gauge. The exterior result is $h_{0i}=2G(\mathbf{J}\times \mathbf{x}{)}_i/r^3$, where the thin shell has $\mathbf{J}=(2/3)M{R}^2\mathbf{\Omega }$, providing a check on the coefficient.

To identify the freely falling frames inside, let ${\mathcal{R}}_z(\theta )$ be the ordinary spatial rotation by angle $\theta$. Define the inertial coordinates inside a slowly rotating spherical shell byDifferentiating and dropping ${ϵ}^2$, products beyond the retained rotation order, and ${\omega }^2$ giveswhich is precisely the interior metric. Thus the retained interior curvature is zero, despite its nontrivial relation to the nonrotating coordinates anchored at infinity. Timelike geodesics in these local inertial coordinates are straight lines $\mathbf{X}={\mathbf{X}}_0+\mathbf{V}T$, with $|\mathbf{V}|<1$. In the original coordinates this isunderstood to the same perturbative order. Their coordinate equation is $d^2\mathbf{x}/d{t}^2=2\mathbf{\omega }\times d\mathbf{x}/dt$ to first rotation order. There is no interior Newtonian acceleration from the constant potential; the velocity-dependent term is the apparent Coriolis term relative to the locally dragged inertial axes.

A freely transported gyroscope at rest in the interior to the retained order has constant spatial direction in the $\mathbf{X}$ frame, hence its direction in the original axes obeys $d\mathbf{S}/dt=\mathbf{\omega }\times \mathbf{S}$. Local inertial axes precess in the same sense as the shell, at $\omega =(4ϵ/3)\Omega$ relative to distant nonrotating axes. This is the interior Lense-Thirring precession. For a weak shell the dragging is only a small fraction of its angular velocity, not perfect corotation. It is compatible with the equivalence principle: rotation relative to distant axes is a comparison between frames, and local flatness at the retained order does not make that comparison vanish. Quadratic rotation terms, which would include centrifugal effects and additional curvature, are outside the requested leading approximation.

The strong equivalence principle extends universality of free fall to bodies with appreciable gravitational binding energy and requires results of sufficiently local experiments, including gravitational experiments, in freely falling laboratories to be independent of their location and velocity. The Einstein equivalence principle makes the corresponding statement for nongravitational experiments; the weak principle concerns test-body free fall without significant self-gravity. In a locally inertial frame the nongravitational laws take their special-relativistic form, while internal gravitational experiments retain the same local gravitational laws and coupling constants. They are not declared to have no gravity.

In a metric description, one may set $g_{ab}={\eta }_{ab}$ and ${\partial }_c{g}_{ab}=0$ at a point by Riemann normal coordinates, so the Levi-Civita connection vanishes there. This removes a locally uniform external gravitational acceleration, not the Riemann curvature tensor. Tidal effects between separated bodies, or accumulated over a finite-duration experiment, remain. The local limit must control both the laboratory's size and duration. Gravitationally self-bound bodies are included in the strong version, so it is stronger than equality of inertial and passive gravitational mass for ordinary test particles.

A heuristic field-equation argument supplements this principle with further assumptions: gravity is described only by a Lorentzian metric, matter couples universally to it, the affine connection is Levi-Civita, and the leading local field equation is a symmetric covariant equation with at most second metric derivatives and is linear in those second derivatives. Matter obeys stress-energy conservation ${\nabla }^a{T}_{ab}=0$. Curvature therefore supplies the natural candidates $R_{ab}$ and $R{g}_{ab}$, with a constant multiple of $g_{ab}$ also allowed. Write the candidate as $R_{ab}+cR{g}_{ab}+\Lambda g_{ab}=\kappa T_{ab}$. The contracted Bianchi identity and metric compatibility giveTo permit varying scalar curvature consistently with matter conservation, take $c=-1/2$. This gives the Einstein tensor $G_{ab}=R_{ab}-g_{ab}R/2$ andThe cosmological constant is constant so that its divergence vanishes. These assumptions motivate the usual Einstein equation; the equivalence principle alone does not exclude theories with extra fields or higher curvature terms.

The sign and normalization are fixed by the supplied Newtonian limit. Use signature $(+---)$, the curvature commutator convention of Q1, and the Ricci contractionThis is the contraction for which $R_{00}\simeq -{\nabla }^2\phi$ in the question. In four dimensions, tracing the candidate equation gives $-R+4\Lambda =\kappa T$, and substituting back givesFor slowly moving pressureless matter, $T_{00}\simeq \rho$, $T\simeq \rho$ and $g_{00}\simeq 1$. On local scales where the cosmological term is negligible, $R_{00}\simeq \kappa \rho /2$. Compare this with $R_{00}\simeq -{\nabla }^2\phi$ and the Poisson equation ${\nabla }^2\phi =4\pi G\rho$. The Newtonian normalization of the Einstein field equations is thereforeUnits with light speed one are used, as in the paper; for conventional physical stress-energy units the coupling is $-8\pi G/c^4$. If the opposite Ricci/curvature sign convention is adopted, the displayed geometric coupling sign changes accordingly. With the cosmological term retained in the same weak-field calculation, the Newtonian equation becomes ${\nabla }^2\phi =4\pi G\rho -\Lambda$, consistent with this sign choice.

A linear connection on a manifold is a rule $(X,Y)\mapsto {\nabla }_{X}Y$ for differentiating vector fields which is additive, is $C^{\infty }$-linear in $X$, is real-linear in $Y$, and obeys ${\nabla }_X(fY)=X(f)Y+f{\nabla }_{X}Y$. It extends to covector fields and arbitrary tensor fields by the Leibniz rule and compatibility with contraction. In coordinates,This is the covariant derivative; the affine connection terms correct the changes of the coordinate basis. A symmetric affine connection has zero torsion tensor, $T(X,Y)={\nabla }_{X}Y-{\nabla }_{Y}X-[X,Y]=0$, or ${\Gamma }^c{}_{ab}={\Gamma }^c{}_{ba}$ in a coordinate basis. Symmetry does not imply metric compatibility; no metric is required for this question.

Additivity of each covariant derivative immediately gives additivity of ${\Delta }_{ab}$. Expanding the second derivative of a tensor product givesSubtract the expression with $a,b$ exchanged. The two mixed terms cancel, proving that the curvature commutator is a tensor derivation:For a scalar, ${\nabla }_a{\nabla }_{b}f={\partial }_a{\partial }_{b}f-{\Gamma }^c{}_{ab}{\partial }_{c}f$. Commuting partial derivatives and using the symmetric affine connection gives ${\Delta }_{ab}f=0$. Consequently ${\Delta }_{ab}(fU)=f{\Delta }_{ab}U$: the commutator is linear over smooth functions, not merely over constant scalars. Write $U=U^c{e}_c$ in a local basis. The tensor-product rule gives ${\Delta }_{ab}U=U^c{\Delta }_{ab}{e}_c$, showing that its value at a point depends only on $U$ at that point. Thus it defines the Riemann curvature tensor byThe antisymmetrized second covariant derivative is itself tensorial in the derivative indices, so this pointwise map transforms as a tensor. The formula or the defining commutator gives $R_{abc}{}^d=-R_{bac}{}^d$, hence $R_{abc}{}^d=R_{[ab]c}{}^d$ with normalized antisymmetrization.

Apply the commutator to the scalar contraction ${\omega }_c{U}^c$. Its action vanishes, soSince $U$ is arbitrary, ${\Delta }_{ab}{\omega }_c=-R_{abc}{}^d{\omega }_d$. Repeated use of the tensor-product rule gives one plus curvature action for an upper index and one minus action for each lower index. In particular,The type notation in the PDF means one contravariant and two covariant indices, not a fractional tensor order.

For the first Bianchi identity, put ${\omega }_c={\nabla }_{c}f$ and $H_{bc}={\nabla }_b{\nabla }_{c}f$. The symmetric Hessian of a scalar field satisfies $H_{bc}=H_{cb}$. Therefore the six third-derivative terms cancel inSubstituting the covector curvature action gives $(R_{abc}{}^d+R_{bca}{}^d+R_{cab}{}^d){\nabla }_{d}f=0$. At a chosen point the gradient of a smooth scalar can be any covector, soThe second form follows from the already proved first-pair antisymmetry.