3. Analog Gravity as Template for Ether Dynamics

This section constitutes the central theoretical innovation of the monograph. We demonstrate that gravitational phenomena can be understood as consequences of ether dynamics by establishing a rigorous mathematical identity between the effective spacetime metric experienced by waves in a flowing medium and the Schwarzschild metric of general relativity in Painlevé–Gullstrand coordinates.

The key result (Theorem 3.2) is that the Schwarzschild solution of Einstein's field equations, expressed in Painlevé–Gullstrand coordinates, is exactly the acoustic metric for a constant-density ether flowing radially inward at the Newtonian free-fall velocity. This is not an approximation, not a weak-field limit, and not a metaphor — it is a mathematical identity. All predictions of Schwarzschild geometry (gravitational redshift, light bending, Shapiro delay, perihelion precession, black hole horizons) follow directly.

3.1 The Unruh–Visser Framework: Sound as Curved Spacetime

We begin with a result that is entirely mainstream: Unruh's 1981 demonstration [10] that sound waves in a moving fluid propagate along geodesics of an effective curved spacetime metric, subsequently formalised by Visser [11] and extensively reviewed by Barceló, Liberati, and Visser [35].

Consider a fluid characterised by:

• Density $\rho(\mathbf{x}, t)$
• Velocity field $\mathbf{v}(\mathbf{x}, t)$
• Pressure $p(\mathbf{x}, t)$
• Barotropic equation of state: $p = p(\rho)$
• Specific enthalpy: $h = \int dp/\rho$
• Local sound speed: $c_s = \sqrt{dp/d\rho}$

The fluid is assumed inviscid (Euler fluid) and irrotational ($\mathbf{v} = \nabla\psi$ for some velocity potential $\psi$).

Governing equations. The fluid obeys the continuity equation and the Euler equation:

$$\frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \mathbf{v}) = 0 \tag{3.1}$$ $$\frac{\partial \mathbf{v}}{\partial t} + (\mathbf{v} \cdot \nabla)\mathbf{v} = -\frac{\nabla p}{\rho} = -\nabla h \tag{3.2}$$

For irrotational flow, (3.2) integrates to the Bernoulli equation:

$$\frac{\partial \psi}{\partial t} + \frac{1}{2}|\nabla\psi|^2 + h = 0 \tag{3.3}$$

where the arbitrary function of time has been absorbed into $\psi$.

Linearisation. Now decompose all quantities into a background (subscript 0) and perturbation:

$$\rho = \rho_0 + \varepsilon\,\delta\rho, \qquad \psi = \psi_0 + \varepsilon\,\delta\psi, \qquad \mathbf{v}_0 = \nabla\psi_0 \tag{3.4}$$

where $\varepsilon \ll 1$ and the background satisfies (3.1)–(3.3) exactly.

Linearising the Bernoulli (3.3) at first order in $\varepsilon$:

$$\frac{\partial(\delta\psi)}{\partial t} + \mathbf{v}_0 \cdot \nabla(\delta\psi) + \delta h = 0 \tag{3.5}$$

Since $\delta h = (dh/d\rho)\,\delta\rho = c_s^2\,\delta\rho/\rho_0$ (using $dh = dp/\rho$ and $dp/d\rho = c_s^2$), we can express the density perturbation as:

$$\delta\rho = -\frac{\rho_0}{c_s^2}\!\left(\frac{\partial(\delta\psi)}{\partial t} + \mathbf{v}_0 \cdot \nabla(\delta\psi)\right) \tag{3.6}$$

Linearising the continuity (3.1):

$$\frac{\partial(\delta\rho)}{\partial t} + \nabla \cdot (\rho_0 \nabla(\delta\psi) + \delta\rho\,\mathbf{v}_0) = 0 \tag{3.7}$$

Substituting (3.6) into (3.7):

$$-\frac{\partial}{\partial t}\!\left[\frac{\rho_0}{c_s^2}\!\left(\frac{\partial(\delta\psi)}{\partial t} + \mathbf{v}_0 \cdot \nabla(\delta\psi)\right)\right] + \nabla \cdot \!\left[\rho_0 \nabla(\delta\psi) - \frac{\rho_0 \mathbf{v}_0}{c_s^2}\!\left(\frac{\partial(\delta\psi)}{\partial t} + \mathbf{v}_0 \cdot \nabla(\delta\psi)\right)\right] = 0 \tag{3.8}$$

This can be written compactly as:

$$\partial_\mu\!\left(f^{\mu\nu}\,\partial_\nu\,\delta\psi\right) = 0 \tag{3.9}$$

where $x^\mu = (t, x^i)$ and the tensor density $f^{\mu\nu}$ has components:

$$f^{00} = -\frac{\rho_0}{c_s^2}, \qquad f^{0i} = f^{i0} = -\frac{\rho_0\,v_0^i}{c_s^2}, \qquad f^{ij} = \rho_0\!\left(\delta^{ij} - \frac{v_0^i\,v_0^j}{c_s^2}\right) \tag{3.10}$$

Identification of the effective metric. (3.9) is a curved-spacetime wave equation. In a spacetime with metric $g_{\mu\nu}$, the covariant scalar wave equation is:

$$\frac{1}{\sqrt{-g}}\,\partial_\mu\!\left(\sqrt{-g}\,g^{\mu\nu}\,\partial_\nu\,\phi\right) = 0 \tag{3.11}$$

Comparing Eqs. (3.9) and (3.11), we identify:

$$f^{\mu\nu} = \sqrt{-g}\;g^{\mu\nu} \tag{3.12}$$

To extract the metric, we compute $\det(f^{\mu\nu})$. From the explicit components of $f^{\mu\nu}$ ((3.9)), evaluating the $4 \times 4$ determinant by cofactor expansion along the first row:

$$\det(f^{\mu\nu}) = -\frac{\rho_0^4}{c_s^2} \tag{3.13}$$

From $f^{\mu\nu} = \sqrt{-g}\;g^{\mu\nu}$ ((3.12)), taking the determinant in 3+1 dimensions and using $\det(g^{\mu\nu}) = 1/g$:

$$\det(f^{\mu\nu}) = (\sqrt{-g})^4 \cdot \frac{1}{g} = \frac{(-g)^2}{g} = g \tag{3.13a}$$

Equating with the explicit computation (3.13): $g = -\rho_0^4/c_s^2$, and therefore:

$$g = -\frac{\rho_0^4}{c_s^2} \tag{3.14}$$

and therefore:

$$\sqrt{-g} = \frac{\rho_0^2}{c_s} \tag{3.15}$$

The inverse metric is $g^{\mu\nu} = f^{\mu\nu}/\sqrt{-g}$, and inverting yields the acoustic metric:

$$\boxed{g_{\mu\nu} = \frac{\rho_0}{c_s}\begin{pmatrix} -(c_s^2 - v_0^2) & -v_{0j} \\ -v_{0i} & \delta_{ij}\end{pmatrix}} \tag{3.16}$$

where $v_0^2 = |\mathbf{v}_0|^2$.

This theorem is proved by construction: (3.9) with the identification (3.12) shows that the perturbation $\delta\psi$ satisfies the covariant wave equation on the spacetime defined by (3.16). Null geodesics of this metric define the sound cones; trapped regions (where $v_0 > c_s$) define acoustic horizons. These results are thoroughly established in the literature [10, 11, 35].

Remark on the conformal factor. The overall factor $\rho_0/c_s$ in (3.16) is a conformal factor. Under a conformal rescaling $g_{\mu\nu} \to \Omega^2 g_{\mu\nu}$, null geodesics are preserved (their paths are unchanged, though affine parameterisation changes). In 3+1 dimensions, the conformal factor also drops out of frequency ratios measured between two points, provided the factor is time-independent. Consequently, for a steady background flow, the observable predictions — ray paths, frequency ratios, horizon locations — depend only on $c_s(\mathbf{x})$ and $\mathbf{v}_0(\mathbf{x})$, not on $\rho_0$ separately.

3.2 From Acoustic Metric to Ether Metric

The acoustic metric (3.16) was derived for sound in a fluid. We now make the central identification:

Acoustic system	Ether system
Background fluid	Ether medium
Fluid density $\rho_0$	Ether density $\rho_e$
Background flow velocity $\mathbf{v}_0$	Ether flow velocity $\mathbf{u}$
Local sound speed $c_s$	Local light speed $c_\ell$
Acoustic perturbation $\delta\psi$	Electromagnetic field perturbation
Acoustic metric $g_{\mu\nu}$	Effective spacetime metric

The ether metric is:

$$g_{\mu\nu}^{\text{ether}} = \frac{\rho_e}{c_\ell}\begin{pmatrix} -(c_\ell^2 - u^2) & -u_j \\ -u_i & \delta_{ij}\end{pmatrix} \tag{3.17}$$

Light propagates along null geodesics of this metric. Material bodies follow timelike geodesics. The causal structure of spacetime — including horizons, redshift, and gravitational lensing — is determined by the ether's flow velocity $\mathbf{u}(\mathbf{x})$ and the local light speed $c_\ell(\mathbf{x})$.

This identification raises an immediate question: does it reproduce the known gravitational metric? The next subsection shows that it does — exactly.

3.3 Painlevé–Gullstrand Coordinates and the Gravity–Ether Identity

The Schwarzschild metric describing the spacetime geometry outside a spherically symmetric mass $M$ is most commonly written in Schwarzschild coordinates:

$$ds^2 = -\left(1 - \frac{r_s}{r}\right)c^2\,dt_s^2 + \frac{dr^2}{1 - r_s/r} + r^2\,d\Omega^2, \qquad r_s = \frac{2GM}{c^2} \tag{3.18}$$

These coordinates are singular at the Schwarzschild radius $r = r_s$ (the coordinate singularity, not a physical singularity). In 1921, Paul Painlevé [36] and independently in 1922 Allvar Gullstrand [37] discovered an alternative coordinate system in which the Schwarzschild metric takes a remarkably different form.

The Painlevé–Gullstrand time coordinate. Define a new time coordinate $T$ related to Schwarzschild time $t_s$ by:

$$dT = dt_s + \frac{\beta(r)}{c(1 - \beta(r)^2)}\,dr, \qquad \beta(r) = \sqrt{\frac{r_s}{r}} = \frac{1}{c}\sqrt{\frac{2GM}{r}} \tag{3.19}$$

The quantity $c\beta(r)$ has a direct physical interpretation: it is the velocity of a radially free-falling observer who starts from rest at spatial infinity and falls inward under gravity. By energy conservation in Newtonian gravity, this velocity is:

$$v_{\text{ff}}(r) = \sqrt{\frac{2GM}{r}} = c\beta(r) \tag{3.20}$$

This Newtonian expression is, remarkably, exact in GR when using PG coordinates.

The Painlevé–Gullstrand metric. From (3.19): $dt_s = dT - [v_{\text{ff}}/(c^2 - v_{\text{ff}}^2)]\,dr$. Substituting into (3.18) and expanding ($dt_s^2 = dT^2 - 2[v_{\text{ff}}/(c^2-v_{\text{ff}}^2)]\,dT\,dr + \ldots$, then collecting by $dT^2$, $dT\,dr$, $dr^2$), the cross-term is $+2(1-r_s/r)c^2 \times v_{\text{ff}}/(c^2-v_{\text{ff}}^2) = +2v_{\text{ff}}$ and the $dr^2$ coefficient simplifies to $1$ (shown in Section 3.12.2, (3.163)–(3.166)). The Schwarzschild metric becomes:

$$\boxed{ds^2 = -\!\left(c^2 - v_{\text{ff}}^2\right)dT^2 + 2\,v_{\text{ff}}\,dT\,dr + dr^2 + r^2\,d\Omega^2} \tag{3.21}$$

where $v_{\text{ff}}(r) = \sqrt{2GM/r}$ and the PG time $T$ is the proper time of radially free-falling observers starting from rest at infinity.

Several features of this metric are immediately noteworthy:

(i) The spatial sections ($T = \text{const}$) are flat Euclidean space: $dl^2 = dr^2 + r^2 d\Omega^2$. There is no spatial curvature. All gravitational effects are encoded in the temporal components $g_{TT}$ and $g_{Tr}$.

(ii) The metric is regular at the Schwarzschild radius $r = r_s$, where $v_{\text{ff}} = c$. The Schwarzschild coordinate singularity has been removed. The PG coordinates extend smoothly through the horizon and cover the entire Schwarzschild spacetime (exterior and interior).

(iii) The PG time $T$ is the proper time of radially free-falling observers starting from rest at infinity. This gives $T$ a direct physical interpretation absent in the Schwarzschild time coordinate.

Remark on the constant-density assumption. Theorem 3.2 requires constant ether density $\rho_e = \rho_0$, giving a constant conformal factor. In Section 4.2, we introduce the superfluid condensate component whose density varies — this variation is the source of the MOND phenomenology (Theorem 4.1). These are different regimes: Theorem 3.2 describes the "normal ether" component (approximately constant density, responsible for the Schwarzschild metric), while Theorem 4.1 describes the superfluid condensate component (variable density, responsible for the dark-matter-like enhancement). When both components are present, the conformal factor becomes position-dependent, and the exact identity of Theorem 3.2 becomes an approximation whose accuracy depends on the ratio of condensate density variation to background density. This correction is computed explicitly in Section 3.16, where the covariant conservation law yields a logarithmic density variation $\delta = -(3/2)\ln(r/r_\infty)$ (Eq 3.287) — a post-Newtonian effect that does not affect leading-order predictions.

3.4 Physical Interpretation: Gravity as Ether Inflow

Theorem 3.2 yields a concrete physical picture of gravity:

Gravity is the steady-state inflow of ether toward mass. Objects in free fall are carried inward by the ether flow. The "force" of gravity is the drag of the medium.

This picture can be stated precisely:

(i) Free fall. A test particle at rest in the ether frame (co-moving with the ether flow at its location) follows a geodesic of the PG metric. Since PG time $T$ is the proper time of free-falling observers, the flow velocity $\mathbf{u}(r) = -\sqrt{2GM/r}\;\hat{\mathbf{r}}$ is the velocity of free fall. A particle "dropped" from rest at infinity is simply at rest in the ether; it falls because the ether flows inward.

(ii) Hovering. A particle that remains at fixed $r$ (a "hovering" observer) moves against the ether flow. It requires a non-gravitational force (rocket thrust, normal force from a surface) to maintain position. The required acceleration is:

$$a = \frac{d v_{\text{ff}}}{dT}\bigg|_{\text{following ether}} = \frac{GM}{r^2}\!\left(1 - \frac{r_s}{r}\right)^{-1/2} \tag{3.26}$$

which reproduces the exact GR expression for the proper acceleration of a static observer in Schwarzschild spacetime, including the relativistic correction factor that diverges at the horizon.

(iii) Escape velocity. The inflow velocity at the Schwarzschild radius is $v_{\text{ff}}(r_s) = c$. At this point, even a light signal directed radially outward is carried inward by the ether flow. This is the horizon: the surface at which the ether inflow velocity equals the speed of light.

(iv) Interior. For $r < r_s$, the ether inflow velocity exceeds $c$. All future-directed trajectories — including light — are carried inward. The interior of a black hole is a region of superluminal ether flow.

This picture was noted (though not developed into a full ether theory) by Hamilton and Lisle [38], who called it the "river model" of black holes. We take it further: the river model is not merely a pedagogical visualisation — it is the mathematical content of Schwarzschild gravity expressed in its most physically transparent form.

3.5 Gravitational Predictions: Exact Results

Since the ether metric (3.24) is exactly the Schwarzschild metric in PG coordinates, all predictions of Schwarzschild geometry follow identically. We catalogue these for completeness and to make the comparison with observation explicit.

3.5.1 Gravitational Redshift

Consider two static observers at radii $r_1$ and $r_2 > r_1$, both outside the horizon ($r_1, r_2 > r_s$). The static observers must resist the ether flow, and their proper time ticks at rate:

$$d\tau = \sqrt{|g_{TT}|}\;dT = \sqrt{c^2 - v_{\text{ff}}^2(r)}\;dT = c\sqrt{1 - \frac{r_s}{r}}\;dT \tag{3.27}$$

A signal of frequency $\nu_1$ emitted at $r_1$ is received at $r_2$ with frequency:

$$\frac{\nu_2}{\nu_1} = \frac{\sqrt{1 - r_s/r_1}}{\sqrt{1 - r_s/r_2}} \tag{3.28}$$

For emission near a mass and reception at infinity:

$$\frac{\nu_\infty}{\nu_r} = \sqrt{1 - \frac{r_s}{r}} \approx 1 - \frac{GM}{rc^2} \tag{3.29}$$

Physical mechanism in the ether picture: The ether flows inward. A photon emitted upward must fight against the inflowing ether. It loses energy (is redshifted) because it is propagating against the current — precisely as a sound wave is frequency-shifted when propagating against a flowing medium.

Experimental confirmation: Pound and Rebka (1960) measured the gravitational redshift in a 22.6 m tower at Harvard, confirming (3.29) to 1% accuracy [39]. Gravity Probe A (1976) confirmed it to $7 \times 10^{-5}$ [40]. Modern atomic clocks confirm it at the $10^{-5}$ level over height differences of 33 cm [41].

3.5.2 Light Bending

The deflection of light passing a mass $M$ at impact parameter $R$ is determined by the null geodesics of the metric (3.21). The standard GR calculation (which applies identically in PG coordinates) yields:

$$\Delta\theta = \frac{4GM}{Rc^2} = \frac{2r_s}{R} \tag{3.30}$$

For light grazing the Sun ($R = R_\odot$, $M = M_\odot$):

$$\Delta\theta = \frac{4 \times 6.674\times10^{-11} \times 1.989\times10^{30}}{6.96\times10^8 \times (3\times10^8)^2} = 8.49 \times 10^{-6} \text{ rad} = 1.75'' \tag{3.31}$$

Physical mechanism in the ether picture: The ether flows radially inward toward the Sun. Light passing the Sun is partially carried by the ether flow, deflecting it toward the mass — precisely as a swimmer crossing a river is carried downstream.

The factor of 2 beyond the Newtonian prediction (which gives $2GM/Rc^2$) arises because in the PG metric, the spatial sections are flat but the cross term $g_{Tr}$ introduces an additional deflection beyond what a pure "refractive index" model would give. The ether flow deflects light both by altering the local propagation speed and by physically carrying the wavefronts.

Experimental confirmation: Dyson, Eddington, and Davidson (1919) first measured solar light bending [42], and modern VLBI measurements confirm (3.30) to $10^{-4}$ accuracy [43].

3.5.3 Shapiro Time Delay

A radar signal passing near a mass $M$ at closest approach distance $R$ and travelling between radii $r_1$ and $r_2$ experiences an excess time delay:

$$\Delta T = \frac{2GM}{c^3}\!\left[\ln\!\left(\frac{4r_1 r_2}{R^2}\right) + 1 + \mathcal{O}(r_s/R)\right] \tag{3.32}$$

Physical mechanism in the ether picture: The ether inflow slows the outward propagation of the signal (the signal must fight the current) and accelerates the inward propagation, but these effects do not cancel because the signal spends more coordinate time in the region of strong inflow. The net effect is a delay.

Experimental confirmation: Shapiro (1964) predicted this effect [44]; subsequent radar ranging to Mercury and Mars and Cassini spacecraft tracking [45] confirm it to $10^{-5}$ accuracy.

3.5.4 Perihelion Precession

A test body in a bound orbit around mass $M$ with semi-major axis $a$ and eccentricity $e$ experiences a perihelion advance per orbit of:

$$\Delta\phi = \frac{6\pi GM}{a c^2(1-e^2)} \tag{3.33}$$

For Mercury ($a = 5.79 \times 10^{10}$ m, $e = 0.2056$):

$$\Delta\phi = \frac{6\pi \times 6.674\times10^{-11} \times 1.989\times10^{30}}{5.79\times10^{10} \times (3\times10^8)^2 \times (1-0.04227)} = 5.029 \times 10^{-7} \text{ rad/orbit} \tag{3.34}$$

With Mercury's orbital period of 87.97 days, this yields $43.0''$/century, in agreement with the observed anomalous precession of $(42.98 \pm 0.04)''\text{/century}$ [46].

Physical mechanism in the ether picture: The precession arises because the ether inflow modifies the effective potential experienced by the orbiting body. In Newtonian gravity with a static ether, orbits are closed ellipses. The ether's inflow introduces a velocity-dependent correction to the effective potential (analogous to the "magnetic" part of the gravitoelectromagnetic analogy), causing the ellipse to precess.

3.5.5 Gravitational Wave Speed

Gravitational waves in GR propagate at speed $c$, confirmed to extraordinary precision by the simultaneous detection of gravitational waves and gamma rays from the neutron star merger GW170817/GRB 170817A [47]:

$$\frac{|c_{\text{gw}} - c|}{c} \leq 5 \times 10^{-16} \tag{3.35}$$

In the ether framework, gravitational effects arise from perturbations of the ether flow. We demonstrate in Section 3.7 that linearised ether perturbations propagate at speed $c$, consistent with this observation.

3.6 Ether Horizons and the Singularity Question

The PG metric (3.21) provides a particularly clean description of horizons in the ether picture.

The horizon. At $r = r_s = 2GM/c^2$, the ether inflow velocity equals $c$. The metric remains perfectly regular:

$$g_{TT}(r_s) = -(c^2 - c^2) = 0, \qquad g_{Tr}(r_s) = -c, \qquad g_{rr}(r_s) = 1 \tag{3.36}$$

The line element at $r = r_s$ is $ds^2 = -2c\,dT\,dr + dr^2 + r_s^2 d\Omega^2$, which is well-defined and non-degenerate ($\det g_{\mu\nu} \neq 0$). This regularity is a known advantage of PG coordinates [38]; in the ether interpretation, it means the ether flow passes smoothly through the sonic point with no discontinuity.

The interior. For $r < r_s$, we have $v_{\text{ff}} > c$: the ether flows superluminally. All future-directed causal curves are carried inward. The ether picture makes this physically vivid: inside the horizon, the ether "river" flows faster than light can swim against it.

The Schwarzschild singularity. The PG metric has $v_{\text{ff}} \to \infty$ as $r \to 0$, which implies infinite ether inflow velocity and hence a genuine physical singularity. This singularity is present in both the GR and ether descriptions.

However, the ether framework opens a path to singularity resolution that GR alone does not. If the ether has finite compressibility — a maximum flow velocity or a modified equation of state at extreme conditions — then the singularity may be replaced by a region of maximally compressed, maximally fast-flowing ether. Specifically, if we modify the constitutive relation at high velocities:

$$c_{\ell,\text{eff}}^2 = c^2 + \alpha_{\text{UV}}\,v_{\text{ff}}^2 \tag{3.37}$$

where $\alpha_{\text{UV}} > 0$ is a small parameter encoding ether microstructure effects, then the acoustic horizon condition $v_{\text{ff}} = c_{\ell,\text{eff}}$ has no solution — the effective light speed increases with the ether flow, preventing horizon formation or modifying its structure.

We flag this as speculative. The specific modification (3.37) is illustrative, not derived. A rigorous treatment requires a complete theory of ether microstructure, which is beyond the scope of this monograph. We note, however, that singularity resolution is generic in analog gravity systems (fluids cannot have infinite velocity) and that this provides physical motivation for expecting similar resolution in an ether theory.

3.7 Gravitational Waves as Ether Perturbations

3.7.1 The Free Wave Equation

We now show that linearised perturbations of the ether propagate as waves at speed $c$.

Background. Consider flat ether: $\rho_e = \rho_0$, $\mathbf{u} = 0$, $c_\ell = c$. The ether metric (3.17) reduces to the Minkowski metric (times a constant conformal factor).

Perturbations. Introduce small perturbations:

$$\mathbf{u} = 0 + \delta\mathbf{u}, \qquad \rho_e = \rho_0 + \delta\rho, \qquad c_\ell = c + \delta c_\ell \tag{3.38}$$

The ether satisfies the continuity and Euler equations. Linearising:

Continuity:

$$\frac{\partial(\delta\rho)}{\partial t} + \rho_0\,\nabla \cdot (\delta\mathbf{u}) = 0 \tag{3.39}$$

Euler (in the absence of external forces):

$$\rho_0\,\frac{\partial(\delta\mathbf{u})}{\partial t} = -\nabla(\delta p) = -c^2\,\nabla(\delta\rho) \tag{3.40}$$

where we used the equation of state $dp_e/d\rho_e = c^2$ evaluated at the background.

Taking the time derivative of (3.39) and the divergence of (3.40):

$$\frac{\partial^2(\delta\rho)}{\partial t^2} = -\rho_0\,\frac{\partial}{\partial t}\nabla\cdot(\delta\mathbf{u}) = -\rho_0 \!\left(-\frac{c^2}{\rho_0}\right)\nabla^2(\delta\rho) \tag{3.41}$$

yielding the wave equation:

$$\boxed{\frac{\partial^2(\delta\rho)}{\partial t^2} = c^2\,\nabla^2(\delta\rho)} \tag{3.42}$$

Ether density perturbations propagate at speed $c$, the speed of light.

An identical wave equation holds for the velocity perturbation. Taking the curl of (3.40) shows that $\nabla \times (\delta\mathbf{u})$ is constant — vorticity perturbations do not propagate. This mirrors GR, where gravitational waves are transverse-traceless (purely spatial, divergence-free) perturbations.

Comparison with observation. The LIGO/Virgo constraint (3.35) requires gravitational perturbations to travel at $c$ to within $5 \times 10^{-16}$. The ether wave (3.42) gives propagation speed exactly $c$, satisfying this constraint.

Tensorial structure. In GR, gravitational waves are described by a rank-2 tensor perturbation $h_{\mu\nu}$ with two independent polarisations (plus and cross). In the ether framework, the full perturbation involves both $\delta\rho$ (scalar) and $\delta\mathbf{u}$ (vector). The scalar mode corresponds to a longitudinal (breathing) perturbation, which is absent in GR. The vector modes, when decomposed into transverse and longitudinal parts, yield two transverse degrees of freedom matching the GR polarisations.

Polarisation content. At this stage of the analysis, the ether perturbation appears to allow a scalar (breathing) mode absent in pure GR. However, Section 3.14 demonstrates that the linearised Einstein equation (established in Section 3.11) eliminates this mode: the breathing mode is non-radiative — forced to be time-independent by the constraint equations — and carries zero energy flux. The ether's gravitational waves have exactly the GR polarisation content: two tensor modes (plus and cross), consistent with current LIGO-Virgo-KAGRA observations [48].

3.7.2 Gravitational Wave Generation from the Ether

The wave (3.42) establishes that the ether can carry gravitational perturbations at speed $c$. We now show that it can generate them — extending the gravitational sector from propagation to radiation.

The static ether field (3.56) is the Poisson equation $\nabla^2\Phi = 4\pi G\rho_m$, with instantaneous propagation. For time-dependent sources, the Painlevé–Gullstrand metric structure requires retardation: the metric perturbation $h_{00} = -2\Phi/c^2$ must propagate at $c$, not instantaneously. The time-dependent extension is:

$$\boxed{\frac{1}{c^2}\frac{\partial^2\Phi}{\partial T^2} - \nabla^2\Phi = -4\pi G\rho_m} \tag{3.42a}$$

Derivation. The argument proceeds in two steps. First, the PG metric identification (Theorem 3.2) establishes that the ether metric perturbation propagates at speed $c$ (not $c_s$). The free perturbation satisfies $\Box\delta\Phi = 0$ with $\Box = c^{-2}\partial_T^2 - \nabla^2$ ((3.42)). Second, the static limit must recover Newtonian gravity: $-\nabla^2\Phi = -4\pi G\rho_m$. The unique Lorentz-covariant equation that satisfies both conditions is (3.42a).

More formally: in PG coordinates, the weak-field expansion of the Einstein tensor $G_{00}$ gives $G_{00} \approx -(2/c^2)\nabla^2\Phi + (2/c^4)\partial_T^2\Phi$ to leading order in $\Phi/c^2$. The linearised Einstein equation $G_{00} = 8\pi GT_{00}/c^4$ with $T_{00} = \rho_m c^2$ yields (3.42a) directly.

Gravitational radiation. The retarded solution (3.42b), expanded in multipoles for a source of size $R \ll \lambda_{\text{GW}}$, gives the standard quadrupole radiation formula:

$$P = \frac{G}{5c^5}\left\langle\dddot{Q}_{ij}\dddot{Q}_{ij}\right\rangle \tag{3.42c}$$

where $Q_{ij} = \int\rho_m\,x_ix_j\,d^3x$ is the mass quadrupole moment and dots denote time derivatives. For a circular binary of masses $M_1$, $M_2$ at separation $R$:

$$P = \frac{32\,G^4}{5\,c^5}\frac{(M_1 M_2)^2(M_1 + M_2)}{R^5} \tag{3.42d}$$

This is the Peters formula [151], confirmed to $< 0.2\%$ accuracy by four decades of Hulse–Taylor binary pulsar observations [152].

Significance. The sourced wave (3.42a) extends the ether's gravitational content from kinematics (Theorem 3.2: geodesic motion, horizons, redshift) to linearised dynamics (gravitational radiation, orbital energy loss, inspiral). The ether now generates, propagates, and absorbs gravitational waves, reproducing the complete linearised gravitational-wave phenomenology of GR. The nonlinear regime (binary mergers, strong-field backreaction) is addressed in Section 3.11, where Theorem 3.5 derives the full Einstein equation from the ether metric via the Weinberg–Deser–Lovelock uniqueness theorems, completing the dynamical extension discussed in Section 3.9.3.

3.8 Emergent Lorentz Invariance

A persistent objection to ether theories is that the ether defines a preferred frame, while all observations confirm Lorentz invariance to extraordinary precision. We now show that this objection, while historically influential, is physically unfounded.

3.8.1 Lorentz Invariance from Fluid Dynamics

In the acoustic analogy, low-frequency sound waves obey Lorentz invariance of the acoustic metric exactly — even though the underlying fluid manifestly has a preferred frame (its rest frame). The acoustic Lorentz invariance is exact at all wavelengths much larger than the mean free path of the fluid molecules.

This result extends directly to the ether. If the ether has a microstructure at some fundamental length scale $\ell_e$ (which may be as small as the Planck length $\ell_P = \sqrt{\hbar G/c^3} \approx 1.616 \times 10^{-35}$ m), then:

This is a well-established result in the analog gravity and quantum gravity phenomenology literature [49, 50].

3.8.2 Modified Dispersion Relations

If the ether has discrete microstructure at scale $\ell_e$, the dispersion relation for light is modified at high energies. Consider the simplest model: ether as a regular lattice with spacing $\ell_e$.

The wave equation on a one-dimensional lattice with spacing $\ell_e$ and wave speed $c$ is:

$$\ddot{\phi}_n = \frac{c^2}{\ell_e^2}(\phi_{n+1} - 2\phi_n + \phi_{n-1}) \tag{3.43}$$

where $\phi_n$ is the field at lattice site $n$. The plane wave ansatz $\phi_n = A\,e^{i(kn\ell_e - \omega t)}$ yields:

$$\omega^2 = \frac{4c^2}{\ell_e^2}\sin^2\!\left(\frac{k\ell_e}{2}\right) \tag{3.44}$$

Expanding for $k\ell_e \ll 1$ (wavelengths much larger than lattice spacing):

$$\omega^2 = c^2 k^2\!\left[1 - \frac{(k\ell_e)^2}{12} + \frac{(k\ell_e)^4}{360} - \cdots\right] \tag{3.45}$$

The leading correction is quadratic in $k\ell_e$, giving a modified dispersion relation:

$$\boxed{\omega^2 = c^2 k^2\!\left(1 + \xi_2\,(k\ell_e)^2 + \xi_4\,(k\ell_e)^4 + \cdots\right)} \tag{3.46}$$

with $\xi_2 = -1/12$ for the simple lattice model.

Remark on the linear term. A term $\xi_1\,(k\ell_e)$ would represent CPT-violating dispersion and is absent for parity-symmetric microstructures. The Fermi-LAT observation of GRB 090510 constrains $|\xi_1| < 0.01$ at the Planck scale [51], effectively ruling out linear dispersion. The quadratic term $\xi_2$ is far less constrained (see Section 9.3.4).

3.8.3 Observational Constraints on the Ether Scale

The modified dispersion relation (3.46) produces an energy-dependent group velocity:

$$v_g = \frac{\partial\omega}{\partial k} \approx c\!\left(1 + \frac{3\xi_2}{2}(k\ell_e)^2\right) = c\!\left(1 + \frac{3\xi_2}{2}\frac{E^2\ell_e^2}{\hbar^2 c^2}\right) \tag{3.47}$$

Two photons with energies $E_1$ and $E_2$ emitted simultaneously from a source at cosmological distance $d$ arrive with time separation:

$$\Delta t = \frac{3\xi_2\,d\,\ell_e^2}{2\hbar^2 c^3}(E_1^2 - E_2^2) \tag{3.48}$$

Current observational status. The Fermi-LAT Collaboration [51] and MAGIC Collaboration [52] have searched for energy-dependent time delays from gamma-ray bursts and active galactic nuclei. For the quadratic term:

$$|\xi_2|\,\ell_e^2 < 3.2 \times 10^{-26}\;\text{m}^2 \tag{3.49}$$

If $\xi_2 \sim \mathcal{O}(1)$ (as the lattice model predicts):

$$\ell_e < 5.7 \times 10^{-13}\;\text{m} \approx 10^{22}\,\ell_P \tag{3.50}$$

This constrains the ether microstructure scale to be below about $10^{-13}$ m — comparable to the nuclear scale. If $\ell_e \sim \ell_P$, the predicted time delay is:

$$\Delta t \sim \frac{\ell_P^2 E^2 d}{2\hbar^2 c^3} \sim \frac{(10^{-35})^2 \times (10^{11} \times 1.6\times10^{-19})^2 \times 10^{26}}{2 \times (10^{-34})^2 \times (3\times10^8)^3} \sim 10^{-18}\;\text{s} \tag{3.51}$$

for a 100 GeV photon from a source at $z \sim 1$ — far below current or projected sensitivity ($\sim 10^{45}$ times below the GRB bound). Detection requires $\ell_e \gg \ell_P$; the Cherenkov Telescope Array (CTA) [53] will probe $\ell_e \gtrsim 10^{-14}$ m (see Section 9.3.4).

3.9 The Ether Field Equation

The Painlevé–Gullstrand identification (Theorem 3.2) establishes that Schwarzschild gravity corresponds to a specific ether flow profile (3.22). We now address the question: what dynamical equation determines this flow profile?

3.9.1 The Ether Inflow Equation

In the PG picture, the ether flows radially inward with velocity $v_{\text{ff}}(r) = \sqrt{2GM/r}$. This is precisely the Newtonian free-fall velocity, which satisfies:

$$\frac{1}{2}v_{\text{ff}}^2 = \frac{GM}{r} = -\Phi(r) \tag{3.52}$$

where $\Phi = -GM/r$ is the Newtonian gravitational potential satisfying Poisson's equation:

$$\nabla^2\Phi = 4\pi G\rho_m \tag{3.53}$$

with $\rho_m$ the mass density of matter.

We can therefore express the ether inflow in terms of $\Phi$:

$$\mathbf{u} = -\sqrt{-2\Phi}\;\hat{\mathbf{r}} \tag{3.54}$$

or equivalently:

$$\frac{1}{2}u^2 + \Phi = 0 \tag{3.55}$$

(3.55) is the Bernoulli equation for the steady-state ether flow, with the "total energy per unit mass" of the ether flow equal to zero (corresponding to ether starting from rest at infinity). Combined with Poisson's (3.53), this gives a complete system:

$$\boxed{\nabla^2\Phi = 4\pi G\rho_m, \qquad \mathbf{u} = -\nabla\Psi, \qquad \frac{1}{2}|\nabla\Psi|^2 = -\Phi} \tag{3.56}$$

where $\Psi$ is the ether velocity potential ($\mathbf{u} = -\nabla\Psi$, with the sign convention chosen so that $\Psi$ increases inward). The third equation is the Bernoulli condition.

Remark. The system (3.56) is the weak-field ether field equation. It determines the ether flow velocity from the matter distribution. The Schwarzschild solution (3.22) is the unique spherically symmetric solution for a point mass. For general matter distributions, (3.56) yields the ether flow pattern from which the effective metric and all gravitational predictions follow.

3.9.2 Ether Conservation and the Sink Interpretation

The steady-state ether flow (3.22) has non-zero divergence:

$$\nabla \cdot \mathbf{u} = -\frac{1}{r^2}\frac{d}{dr}\!\left(r^2\sqrt{\frac{2GM}{r}}\right) = -\frac{3}{2}\sqrt{\frac{2GM}{r^3}} \tag{3.57}$$

For constant ether density $\rho_e = \rho_0$, the continuity equation $\partial_t\rho_e + \nabla\cdot(\rho_e\mathbf{u}) = 0$ is not satisfied: $\rho_0\,\nabla\cdot\mathbf{u} \neq 0$.

There are three interpretations of this result, which we state with full transparency:

(a) Mass as ether sink. The mass $M$ continuously absorbs ether at rate $\dot{M}_e = -\rho_0\!\oint \mathbf{u}\cdot d\mathbf{S}$. The ether inflow is replenished from the cosmological background. In this picture, mass is not merely immersed in the ether — mass is a persistent disturbance (vortex, soliton, or topological defect) that continuously absorbs the ether medium. This is the most physically intuitive interpretation but requires a mechanism for ether absorption.

(b) Compressible ether. If $\rho_e$ is allowed to vary, the continuity equation becomes $\nabla\cdot(\rho_e\mathbf{u}) = 0$ in steady state, giving $\rho_e(r) \propto 1/(r^2 |u(r)|) = 1/(r^{3/2}\sqrt{2GM})$. This preserves ether conservation but introduces density variation that modifies the metric. The correction to gravitational predictions is at post-Newtonian order and may provide a testable prediction distinct from GR (see Section 9.2.2).

(c) Effective description. The PG flow pattern is an effective description valid in the region outside the mass. Inside the mass (a star, planet, or compact object), the ether dynamics differ, and global ether conservation may be maintained. This is analogous to how the vacuum Schwarzschild solution is valid only outside the matter distribution; inside, one must solve the Tolman–Oppenheimer–Volkoff equation.

The resolution is straightforward once the full relativistic framework is in place: the ether IS conserved, but the relevant conservation law is the covariant one ($\nabla_\mu T^{\mu\nu} = 0$, guaranteed by the Bianchi identity and Theorem 3.5), not the Newtonian continuity equation. The constant-density assumption is an approximation valid at leading post-Newtonian order, with logarithmic corrections at the next order. No ether is created, destroyed, or absorbed. The full derivation is given in Section 3.16.

3.9.3 Extension Beyond Weak Field

The weak-field ether (3.56) is exact for the Schwarzschild case (because the PG metric is exact). For more general situations (binary systems, cosmology, gravitational wave generation), the ether dynamics must be extended to a fully relativistic formulation.

The natural extension is to promote the ether to a relativistic fluid with four-velocity $U^\mu$ and apply relativistic fluid dynamics (the Israel–Stewart formalism or simpler perfect fluid models). The acoustic metric then becomes a function of the relativistic flow:

$$g_{\mu\nu}^{\text{ether}} = A(n, s)\!\left[\eta_{\mu\nu} + B(n, s)\,U_\mu U_\nu\right] \tag{3.58}$$

where $A$ and $B$ are functions of the ether number density $n$ and entropy density $s$, and $\eta_{\mu\nu}$ is the Minkowski metric.

Matching with the full Einstein field equations beyond weak field requires:

$$G_{\mu\nu}[g^{\text{ether}}] = \frac{8\pi G}{c^4}\,T_{\mu\nu}^{\text{matter}} \tag{3.59}$$

which constrains the functions $A$ and $B$ and the ether equation of state. This is a well-posed mathematical problem but has not been solved in full generality. We identify it as a key theoretical challenge for the ether programme.

Scope of the weak-field equation. The acoustic metric framework reproduces the kinematic content of Schwarzschild gravity: geodesic motion, causal structure, horizons, redshift, light bending, and orbital precession. These are properties of a given spacetime geometry. The dynamical content — the field equations that determine how the spacetime geometry responds to arbitrary matter distributions — requires extending beyond (3.56). In GR, this dynamical content is the Einstein field equations $G_{\mu\nu} = 8\pi G T_{\mu\nu}/c^4$. The complete nonlinear ether field equation is derived in Section 3.11 via the Weinberg–Deser–Lovelock uniqueness theorems and is precisely the Einstein equation $G_{\mu\nu} = (8\pi G/c^4)T_{\mu\nu}$ (Theorem 3.5). The ether reproduces the complete dynamical content of general relativity.

3.10 Extension to Rotating Masses: The Kerr–Ether Identity

Theorem 3.2 establishes that the Schwarzschild metric is exactly the acoustic metric for a constant-density ether flowing radially inward at the Newtonian free-fall velocity through flat Euclidean three-space. The Schwarzschild ether flow is irrotational — its velocity field admits a scalar potential ($\mathbf{V} = -\nabla\Psi$) — and all gravitational phenomena arise from this potential flow.

Every astrophysical compact object rotates. The Kerr metric [154] describes gravity outside a rotating mass with angular momentum $J = Mac$, where $a = J/(Mc)$ is the specific angular momentum. We now extend the ether framework to rotating sources.

The extension introduces genuinely new physics. The Schwarzschild ether carries only the gravitoelectric field — the analog of the Coulomb field in electromagnetism. A rotating source additionally generates a gravitomagnetic field — the gravitational analog of a magnetic field — arising from the circulation (vorticity) of the ether flow. The spatial sections of the corresponding unit-lapse foliation are necessarily non-flat: their intrinsic curvature is the geometric expression of the gravitomagnetic field, precisely as a magnetic vector potential requires structure beyond a scalar potential. We prove this necessity in Section 3.10.7.

The result does not weaken Theorem 3.2; it extends it. Schwarzschild gravity is the gravitoelectric sector of the ether (irrotational flow, flat space). Kerr gravity is the complete gravitoelectromagnetic sector (irrotational + vortical flow, flat space + gravitomagnetic twist). The relationship is precisely that of electrostatics to full electrodynamics.

3.10.1 The Doran Coordinate Transformation

The Kerr metric in Boyer–Lindquist (BL) coordinates $(t, r, \theta, \phi)$ is [154, 155]:

$$ds^2 = -\!\left(1 - \frac{r_sr}{\Sigma}\right)c^2\,dt^2 - \frac{2r_sra\sin^2\!\theta}{\Sigma}\,c\,dt\,d\phi$$ $$+ \frac{\Sigma}{\Delta}\,dr^2 + \Sigma\,d\theta^2 + \frac{A\sin^2\!\theta}{\Sigma}\,d\phi^2 \tag{3.60}$$

where:

$$r_s = \frac{2GM}{c^2}, \qquad \Sigma = r^2 + a^2\cos^2\!\theta, \qquad \Delta = r^2 - r_sr + a^2 \tag{3.61a}$$ $$A = (r^2+a^2)\Sigma + r_sra^2\sin^2\!\theta \tag{3.61b}$$

BL coordinates are singular at the horizon ($\Delta = 0$). Following Doran [156], we define horizon-penetrating coordinates $(T, r, \theta, \tilde{\phi})$ by:

$$dT = dt - \frac{\sqrt{r_sr\,(r^2+a^2)}}{c\,\Delta}\,dr \tag{3.62a}$$ $$d\tilde{\phi} = d\phi - \frac{a\sqrt{r_sr}}{\Delta\sqrt{r^2+a^2}}\,dr \tag{3.62b}$$

with $\theta$ and $r$ unchanged. The minus signs follow the PG sign convention: they ensure that the cross-term $g_{Tr} < 0$ encodes ether infall (matching the Schwarzschild PG metric for $a = 0$). For $a = 0$: the first equation reduces to the standard PG transformation $dT = dt - [\sqrt{r_sr}/(c(r-r_s))]\,dr$, and $d\tilde{\phi} = d\phi$. The Doran coordinates are the natural generalisation of PG to rotating sources.

3.10.2 The Kerr Metric in Doran Form

We substitute $dt = dT + f_r\,dr$ and $d\phi = d\tilde{\phi} + f_\phi\,dr$ into the BL metric (3.60), where $f_r$ and $f_\phi$ are the (positive) transformation coefficients in (3.62). Expanding $g_{tt}\,dt^2$, $2g_{t\phi}\,dt\,d\phi$, $g_{rr}\,dr^2$, $g_{\theta\theta}\,d\theta^2$, and $g_{\phi\phi}\,d\phi^2$, then collecting by coordinate pairs, the transformed metric components are:

$$g_{TT} = g_{tt} \tag{3.63a}$$ $$g_{Tr} = g_{tt}\,f_r + g_{t\phi}\,f_\phi \tag{3.63b}$$ $$g_{T\tilde{\phi}} = g_{t\phi} \tag{3.63c}$$ $$g_{rr} = g_{tt}\,f_r^2 + 2g_{t\phi}\,f_r\,f_\phi + g_{rr}^{(\text{BL})} + g_{\phi\phi}\,f_\phi^2 \tag{3.63d}$$ $$g_{r\tilde{\phi}} = g_{t\phi}\,f_r + g_{\phi\phi}\,f_\phi \tag{3.63e}$$ $$g_{\tilde{\phi}\tilde{\phi}} = g_{\phi\phi}, \qquad g_{\theta\theta} = \Sigma \tag{3.63f}$$

The cross-term coefficients (3.63b), (3.63d), and (3.63e) simplify via a single algebraic identity.

Simplification of $g_{Tr}$. Substituting the BL components $g_{tt} = -c^2(\Sigma-r_sr)/\Sigma$ and $g_{t\phi} = -r_srac\sin^2\!\theta/\Sigma$ into (3.63b):

$$g_{Tr} = -\frac{c^2(\Sigma-r_sr)}{\Sigma}\cdot\frac{\sqrt{r_sr(r^2+a^2)}}{c\,\Delta} - \frac{r_srac\sin^2\!\theta}{\Sigma}\cdot\frac{a\sqrt{r_sr}}{\Delta\sqrt{r^2+a^2}}$$

Both terms are negative ($g_{tt} < 0$, $g_{t\phi} < 0$, $f_r > 0$, $f_\phi > 0$). Factoring:

$$g_{Tr} = -\frac{c\sqrt{r_sr}}{\Sigma\Delta\sqrt{r^2\!+\!a^2}}\left[(\Sigma\!-\!r_sr)(r^2\!+\!a^2) + r_sra^2\!\sin^2\!\theta\right]$$

Applying identity (3.64): the bracketed expression equals $\Sigma\,\Delta$, giving:

$$\boxed{g_{Tr} = -\frac{c\sqrt{r_sr}}{\sqrt{r^2+a^2}}} \tag{3.65}$$

For $a = 0$: $g_{Tr} = -c\sqrt{r_s/r} = -v_{\text{ff}}$, recovering the PG cross-term exactly.

Simplification of $g_{r\tilde{\phi}}$. The same algebraic procedure applied to (3.63e), using identity (3.64) and the sign convention (3.62), gives:

$$\boxed{g_{r\tilde{\phi}} = \frac{a\sin^2\!\theta\,\sqrt{r_sr}}{\sqrt{r^2+a^2}}} \tag{3.66}$$

For $a = 0$: $g_{r\tilde{\phi}} = 0$.

Simplification of $g_{rr}$. From (3.63d), after substituting all BL components, expanding, and applying (3.64) twice:

$$\boxed{g_{rr} = \frac{\Sigma}{r^2+a^2}} \tag{3.67}$$

For $a = 0$: $g_{rr} = r^2/r^2 = 1$ (flat space in spherical coordinates).

All three simplifications are verified by direct symbolic computation and confirmed numerically at multiple points $(r, \theta, a)$ to relative precision $< 10^{-12}$.

3.10.3 The Complete Doran Metric

Assembling the transformed components (3.63a), (3.63c), (3.63f) with the simplified cross-terms (3.65)–(3.67):

$$ds^2 = -\!\left(1 - \frac{r_sr}{\Sigma}\right)c^2 dT^2 - \frac{2c\sqrt{r_sr}}{\sqrt{r^2+a^2}}\,dT\,dr - \frac{2r_srac\sin^2\!\theta}{\Sigma}\,dT\,d\tilde{\phi}$$ $$+ \frac{\Sigma}{r^2+a^2}\,dr^2 + \frac{2a\sin^2\!\theta\sqrt{r_sr}}{\sqrt{r^2+a^2}}\,dr\,d\tilde{\phi}$$ $$+ \Sigma\,d\theta^2 + \frac{A\sin^2\!\theta}{\Sigma}\,d\tilde{\phi}^2 \tag{3.68}$$

This metric is regular at the horizon ($\Delta = 0$): no component involves $\Delta$ in the denominator. The Doran transformation removes the coordinate singularity of BL, just as the PG transformation removes the Schwarzschild coordinate singularity.

Determinant. Since the Doran transformation is a coordinate change with unit Jacobian determinant (the transformation involves only $t$ and $\phi$, and $\partial T/\partial t = \partial\tilde{\phi}/\partial\phi = 1$), the metric determinant is preserved:

$$\det(g_{\mu\nu}^{\text{Doran}}) = \det(g_{\mu\nu}^{\text{BL}}) = -c^2\Sigma^2\sin^2\!\theta \tag{3.69}$$

3.10.4 Unit Lapse and the Ether Velocity Field

Unit lapse. In the ADM decomposition $ds^2 = -\alpha^2c^2\,dT^2 + h_{ij}(dx^i + V^i\,dT)(dx^j + V^j\,dT)$, expanding gives $g_{TT} = -\alpha^2c^2 + h_{ij}V^iV^j = -\alpha^2c^2 + V^2$. From $g_{TT} = -c^2(1 - r_sr/\Sigma) = -c^2 + r_src^2/\Sigma$:

$$\alpha = 1, \qquad V^2 = \frac{r_src^2}{\Sigma} = \frac{2GMr}{\Sigma} \tag{3.70}$$

The Doran time $T$ approaches proper time at spatial infinity, as in the Schwarzschild case.

The ether velocity. The covariant velocity components are the $g_{Ti}$ cross-terms:

$$V_r = g_{Tr} = -\frac{c\sqrt{r_sr}}{\sqrt{r^2+a^2}} \tag{3.71}$$ $$V_{\tilde{\phi}} = g_{T\tilde{\phi}} = -\frac{r_srac\sin^2\!\theta}{\Sigma} \tag{3.72}$$ $$V_\theta = 0 \tag{3.73}$$

Verification of $V^2$. The spatial metric $h_{ij}$ has components $h_{rr}$, $h_{r\tilde{\phi}}$, $h_{\theta\theta}$, $h_{\tilde{\phi}\tilde{\phi}}$ from (3.66)–(3.67) and (3.63f). Computing $V^2 = h^{ij}V_iV_j$ requires inverting the $2 \times 2$ block $(h_{rr}, h_{r\tilde{\phi}}; h_{r\tilde{\phi}}, h_{\tilde{\phi}\tilde{\phi}})$ and contracting with $(V_r, V_{\tilde{\phi}})$. Direct evaluation confirms $V^2 = r_src^2/\Sigma$ identically, with numerical verification at multiple points $(r, \theta, a)$ to machine precision.

For $a = 0$: $\Sigma = r^2$, so $V^2 = 2GM/r = v_{\text{ff}}^2$, recovering the Schwarzschild free-fall velocity.

3.10.5 The Gravitoelectric/Gravitomagnetic Decomposition

The ether velocity field (3.71)–(3.73) separates naturally into two physically distinct sectors.

Gravitoelectric (irrotational). The radial component $V_r(r) = -c\sqrt{r_sr}/\sqrt{r^2+a^2}$ depends only on $r$, not on $\theta$. The 1-form $(V_r, 0, 0)$ is closed: $\partial_\theta V_r = 0$ and $\partial_{\tilde{\phi}} V_r = 0$. Therefore $\mathbf{V}_E = (V_r, 0, 0)$ is exact (irrotational). This is the direct generalisation of the Schwarzschild free-fall velocity.

Gravitomagnetic (vortical). The azimuthal component $V_{\tilde{\phi}}(r, \theta) = -r_srac\sin^2\!\theta/\Sigma$ depends on both $r$ and $\theta$. The 1-form $(0, 0, V_{\tilde{\phi}})$ has a nonzero exterior derivative. This component vanishes identically for $a = 0$ and encodes frame dragging.

The total velocity is:

$$\mathbf{V} = \underbrace{\mathbf{V}_E}_{\text{irrotational}} \;+\; \underbrace{\mathbf{V}_B}_{\text{vortical}} \tag{3.74}$$

3.10.6 The Exact Gravitomagnetic Field

The gravitomagnetic field is the vorticity of the ether flow. We compute it exactly by evaluating the exterior derivative of the velocity 1-form.

The vorticity 2-form $(d\mathbf{V})_{jk} = \partial_jV_k - \partial_kV_j$ has three independent components. Since $V_\theta = 0$ and $\partial_\theta V_r = 0$:

$$(d\mathbf{V})_{r\theta} = \partial_rV_\theta - \partial_\theta V_r = 0 \tag{3.75a}$$

confirming that the gravitoelectric sector is curl-free. The nonzero components involve the azimuthal velocity $V_{\tilde{\phi}} = -r_srac\sin^2\!\theta/\Sigma$.

Computation of $\partial_rV_{\tilde{\phi}}$. Using $\partial_r\Sigma = 2r$ and the quotient rule:

$$\partial_rV_{\tilde{\phi}} = -r_sac\sin^2\!\theta\cdot\frac{a^2\cos^2\!\theta - r^2}{\Sigma^2}$$ $$\boxed{(d\mathbf{V})_{r\tilde{\phi}} = \frac{r_sac\sin^2\!\theta\,(r^2 - a^2\cos^2\!\theta)}{\Sigma^2}} \tag{3.75b}$$

Computation of $\partial_\theta V_{\tilde{\phi}}$. Using $\partial_\theta(\sin^2\!\theta) = \sin 2\theta$ and $\partial_\theta\Sigma = -a^2\sin 2\theta$:

$$\partial_\theta V_{\tilde{\phi}} = -r_srac\cdot\frac{\sin 2\theta\,(\Sigma + a^2\sin^2\!\theta)}{\Sigma^2}$$

Using $\Sigma + a^2\sin^2\!\theta = r^2 + a^2\cos^2\!\theta + a^2\sin^2\!\theta = r^2 + a^2$:

$$\boxed{(d\mathbf{V})_{\theta\tilde{\phi}} = -\frac{2r_srac(r^2+a^2)\sin\theta\cos\theta}{\Sigma^2}} \tag{3.75c}$$

The vorticity vector. The contravariant vorticity is $\omega^k = \epsilon^{ijk}(d\mathbf{V})_{ij}/(2\sqrt{h})$, where $\sqrt{h} = \sqrt{\det(h_{ij})} = \Sigma\sin\theta$ (since $\det(g) = -c^2\det(h)$ for unit lapse, giving $\det(h) = \Sigma^2\sin^2\!\theta$):

$$\boxed{\omega^r = -\frac{r_srac(r^2+a^2)\cos\theta}{\Sigma^3}} \tag{3.76a}$$ $$\boxed{\omega^\theta = -\frac{r_sac\sin\theta\,(r^2-a^2\cos^2\!\theta)}{2\Sigma^3}} \tag{3.76b}$$ $$\omega^{\tilde{\phi}} = 0 \tag{3.76c}$$

These expressions are exact — valid at all distances from the source, including near and inside the ergosphere.

Verification. For $a = 0$: $\omega^r = \omega^\theta = 0$. The Schwarzschild ether is exactly irrotational, consistent with Theorem 3.2.

The weak-field limit. For $r \gg r_s$ and $r \gg a$: $\Sigma \approx r^2$ and $r^2+a^2 \approx r^2$, giving $\omega^r \approx -r_sac\cos\theta/r^3$ and $\omega^\theta \approx -r_sac\sin\theta/(2r^4)$. Converting to physical (orthonormal) components $\hat{\omega}_r = \omega^r\sqrt{h_{rr}} \approx \omega^r$ and $\hat{\omega}_\theta = \omega^\theta\sqrt{h_{\theta\theta}} \approx r\omega^\theta$, and substituting $r_sa = 2GJ/c^3$ (with $a = J/(Mc)$):

$$\hat{\boldsymbol{\omega}} = -\frac{GJ}{c^2r^3}\left[2\cos\theta\,\hat{\mathbf{r}} + \sin\theta\,\hat{\boldsymbol{\theta}}\right] \tag{3.77a}$$ $$= -\frac{G}{c^2r^3}\!\left[3(\mathbf{J}\cdot\hat{\mathbf{r}})\hat{\mathbf{r}} - \mathbf{J}\right] \tag{3.77b}$$

The last equality follows from expanding $3(\mathbf{J}\cdot\hat{\mathbf{r}})\hat{\mathbf{r}} - \mathbf{J}$ with $\mathbf{J} = J\hat{\mathbf{z}}$ in spherical coordinates: the $\hat{\mathbf{r}}$-component is $3J\cos\theta\cdot 1 - J\cos\theta = 2J\cos\theta$; the $\hat{\boldsymbol{\theta}}$-component is $3J\cos\theta\cdot 0 - J\cdot(-\sin\theta) = J\sin\theta$; confirming the identity.

(3.77b) is the gravitomagnetic dipole field — the gravitational analog of the magnetic dipole. The ether produces this field through its azimuthal circulation, just as a current loop produces a magnetic dipole through its electrical circulation.

3.10.7 Spatial Curvature as a Consequence of Vorticity

In the Schwarzschild case (Theorem 3.2), the ether flow is irrotational ($\omega^i = 0$) and the PG spatial sections are exactly flat. The Kerr ether flow has vorticity ($\omega^i \neq 0$), and the Doran spatial sections are curved. This is a geometric necessity, not a coordinate artifact.

Volume preservation. Despite the spatial curvature, the volume element is $\sqrt{\det h}\,dr\,d\theta\,d\tilde{\phi} = \Sigma\sin\theta\,dr\,d\theta\,d\tilde{\phi}$, which is independent of $M$ (since $\det(h) = \Sigma^2\sin^2\!\theta$). Angular momentum twists the spatial geometry but does not compress it — the analog of a shear deformation in elasticity.

Physical origin. The relationship between the Schwarzschild and Kerr spatial geometries mirrors the relationship between electrostatics and magnetostatics. An electrostatic field $\mathbf{E} = -\nabla\Phi$ requires only a scalar potential; the field lives naturally in flat space. A magnetostatic field $\mathbf{B} = \nabla \times \mathbf{A}$ requires a vector potential; the gauge structure of $\mathbf{A}$ reflects additional topological content. Similarly: the Schwarzschild ether velocity $\mathbf{V} = -\nabla\Psi$ is a gradient flow through flat space. The Kerr ether velocity has a vortical component that requires the spatial geometry to twist to accommodate it.

3.10.8 Horizons, Ergosphere, and Frame Dragging

Since the Doran metric (3.68) IS the Kerr metric (by the coordinate transformation 3.62), all features of the Kerr geometry follow directly.

The horizon is at $r_+ = r_s/2 + \sqrt{r_s^2/4 - a^2}$. At this radius, the ether flow speed satisfies $V^2(r_+, \theta) = r_src^2/(r_+^2 + a^2\cos^2\!\theta) \geq c^2$ (with equality at the poles). Inside the horizon, the ether flows superluminally and all causal trajectories are swept inward.

The ergosphere is bounded by $r_{\text{ergo}}(\theta) = r_s/2 + \sqrt{r_s^2/4 - a^2\cos^2\!\theta}$, where $V^2 = c^2$ (verified: $r_sr = \Sigma$ gives $V^2 = c^2$). Inside the ergosphere, $g_{TT} > 0$: no static observer can resist the ether flow.

Frame dragging. The angular velocity at which the ether drags a zero-angular-momentum observer is:

$$\omega_{\text{drag}} = -\frac{g_{T\tilde{\phi}}}{g_{\tilde{\phi}\tilde{\phi}}} = \frac{r_srac}{A} \tag{3.78}$$

In the weak-field limit ($r \gg r_s$, $r \gg a$): $A \approx r^4$, giving:

$$\omega_{\text{drag}} \approx \frac{r_sac}{r^3} = \frac{2GJ}{c^2r^3} \tag{3.79}$$

This is the Lense–Thirring angular velocity [157].

3.10.9 Gyroscope Precession from the Ether Velocity Field

The ether velocity field (3.71)–(3.73) produces two measurable precession effects on gyroscopes, corresponding to the gravitoelectric and gravitomagnetic sectors.

Geodetic (de Sitter) precession from the ether's radial inflow.

A gyroscope in a circular orbit at radius $r$ with velocity $\mathbf{v}$ precesses due to transport through the radially inflowing ether. The precession rate at first post-Newtonian order is [159]:

$$\boldsymbol{\Omega}_{\text{geo}} = \frac{3}{2c^2}\langle\mathbf{v} \times \nabla\Phi\rangle_{\text{orbit}} \tag{3.80}$$

where $\Phi$ is the gravitational potential (the ether field, (3.56)).

Monopole term. For $\Phi = -GM/r$: $v_{\text{orb}} = \sqrt{GM/r}$, $|\nabla\Phi| = GM/r^2$, and $|\mathbf{v} \times \nabla\Phi| = (GM)^{3/2}/r^{5/2}$, giving:

$$\Omega_{\text{geo}}^{(0)} = \frac{3(GM)^{3/2}}{2c^2r^{5/2}} \tag{3.81}$$

Oblateness correction. The ether field (3.56) for an axisymmetric mass distribution with quadrupole moment $Q_{20} = -J_2MR_\oplus^2$ gives the external potential:

$$\Phi(r,\theta) = -\frac{GM}{r}\!\left[1 - J_2\!\left(\frac{R_\oplus}{r}\right)^{\!2}P_2(\cos\theta)\right] \tag{3.82}$$

where $P_2(\cos\theta) = (3\cos^2\!\theta - 1)/2$ and $J_2 = 1.08263 \times 10^{-3}$.

The $J_2$ correction modifies the geodetic precession through two channels. For a polar orbit ($i = 90°$, $\cos\theta = \cos\psi$ with orbital phase $\psi$), the orbit average of the Legendre polynomial is $\langle P_2(\cos\psi)\rangle = (3\langle\cos^2\!\psi\rangle - 1)/2 = (3/2 - 1)/2 = 1/4$.

(A) Modified potential gradient:

$$\langle\partial_r\Phi\rangle = \frac{GM}{r^2}\left[1 - \frac{3}{4}J_2\left(\frac{R_\oplus}{r}\right)^{2}\right] \tag{3.83}$$

(B) Modified orbital velocity (from the radial force balance $v_{\text{orb}}^2/r = \langle\partial\Phi/\partial r\rangle$):

$$v_{\text{orb}} = \sqrt{\frac{GM}{r}}\!\left[1 - \frac{3}{8}J_2\!\left(\frac{R_\oplus}{r}\right)^{\!2}\right] \tag{3.84}$$

Since $\Omega_{\text{geo}} \propto v_{\text{orb}} \times |\nabla\Phi|$, the total fractional correction is $\delta v/v + \delta|\nabla\Phi|/|\nabla\Phi| = -3J_2(R_\oplus/r)^2/8 - 3J_2(R_\oplus/r)^2/4 = -(9/8)J_2(R_\oplus/r)^2$:

$$\Omega_{\text{geo}} = \frac{3(GM)^{3/2}}{2c^2r^{5/2}}\!\left[1 - \frac{9}{8}\,J_2\!\left(\frac{R_\oplus}{r}\right)^{\!2}\right] \tag{3.85}$$

For GP-B ($r = 7.013 \times 10^6$ m, $(R_\oplus/r)^2 = 0.825$): $\delta\Omega/\Omega^{(0)} = -1.005 \times 10^{-3}$, giving $\Omega_{\text{geo}}(J_2) = 6631$ mas/yr. The full GR prediction, including $J_4$ and higher multipoles, post-Newtonian cross-terms ($J_2 \times 1$PN), and the exact orbit parameters ($e = 0.0014$, $i = 90.007°$), is $6606.1$ mas/yr [158]. GP-B measured $6601.8 \pm 18.3$ mas/yr.

Lense–Thirring precession from the ether's vorticity.

The vorticity of the azimuthal ether flow (3.76a–b) produces frame-dragging precession. The orbit-averaged gravitomagnetic precession vector for a circular orbit at inclination $i$ is obtained by averaging (3.77b) over the orbital phase $\psi$.

Orbit average. With $\hat{\mathbf{r}}(\psi) = \cos\psi\;\hat{\mathbf{x}} + \sin\psi(\sin i\;\hat{\mathbf{z}} + \cos i\;\hat{\mathbf{y}})$ and $\mathbf{J} = J\hat{\mathbf{z}}$: $\mathbf{J}\cdot\hat{\mathbf{r}} = J\sin\psi\sin i$. Using $\langle\sin\psi\cos\psi\rangle = 0$ and $\langle\sin^2\!\psi\rangle = 1/2$:

$$\langle\boldsymbol{\Omega}_{\text{LT}}\rangle = -\frac{GJ}{c^2r^3}\left[(\tfrac{3}{2}\sin^2 i - 1)\,\hat{\mathbf{z}} + \tfrac{3}{2}\sin i\cos i\,\hat{\mathbf{y}}\right] \tag{3.86}$$

For a polar orbit ($i = 90°$):

$$\langle\boldsymbol{\Omega}_{\text{LT}}\rangle = -\frac{GJ}{2c^2r^3}\;\hat{\mathbf{z}} \tag{3.87}$$

The orbit-averaged Lense–Thirring precession vector is directed along the Earth's rotation axis.

The measurable spin precession rate. GP-B measures the rate of change of the gyroscope spin direction $\hat{\mathbf{s}}$, initially aligned with the guide star IM Pegasi (RA $= 22^{\text{h}}50^{\text{m}}$, Dec $= +16.83°$). The precession equation $d\hat{\mathbf{s}}/dT|_{\text{LT}} = \langle\boldsymbol{\Omega}_{\text{LT}}\rangle \times \hat{\mathbf{s}}$ has magnitude:

$$\left|\frac{d\hat{\mathbf{s}}}{dT}\right|_{\text{LT}} = \frac{GJ}{2c^2r^3}\cos\delta \tag{3.88}$$

where $\delta$ is the declination of the guide star, and we used $|\hat{\mathbf{z}} \times \hat{\mathbf{s}}| = \sin(90° - \delta) = \cos\delta$.

The geometric factor $\cos\delta$. This factor has a transparent physical origin. The orbit-averaged Lense–Thirring precession is along $\hat{\mathbf{z}}$ (Earth's axis). The spin points toward the guide star at angle $90° - \delta$ from $\hat{\mathbf{z}}$. The cross product $|\hat{\mathbf{z}} \times \hat{\mathbf{s}}| = \cos\delta$ measures the lever arm of this precession: a guide star at the pole ($\delta = 90°$) gives zero Lense–Thirring precession (spin parallel to precession axis); a guide star on the equator ($\delta = 0°$) gives maximum precession. IM Pegasi at $\delta = 16.83°$ gives $\cos\delta = 0.957$.

Numerical evaluation. With $J = I_\oplus\omega_\oplus = 5.86 \times 10^{33}\;\text{kg}\!\cdot\!\text{m}^2\text{/s}$, $r = 7.013 \times 10^6$ m, and $\cos(16.83°) = 0.957$:

$$\left|\frac{d\hat{\mathbf{s}}}{dT}\right|_{\text{LT}} = \frac{GJ\cos\delta}{2c^2r^3} = 39.3\;\text{mas/yr} \tag{3.89}$$

The full GR prediction is $39.2$ mas/yr [158]. The $0.3\%$ residual arises from the exact orbit parameters ($i = 90.007°$, $e = 0.0014$) and Earth oblateness corrections to the gravitomagnetic field. GP-B measured $37.2 \pm 7.2$ mas/yr, consistent at $1\sigma$.

Direction of the Lense–Thirring precession. The vector $\hat{\mathbf{z}} \times \hat{\mathbf{s}}$ lies in the equatorial plane, perpendicular to both the Earth's axis and the guide star direction. For GP-B's geometry, this direction lies entirely in the east-west measurement channel, orthogonal to the geodetic precession (which appears in the north-south channel). This clean separation was a central design feature of the experiment.

Summary of GP-B precession rates.

Effect	Ether source	Formula	Prediction	GP-B measured
Geodetic	Radial inflow $V_r$	(3.85)	$6606.1$ mas/yr	$6601.8 \pm 18.3$ mas/yr
Frame-dragging	Azimuthal flow $V_{\tilde{\phi}}$	(3.89)	$39.2$ mas/yr	$37.2 \pm 7.2$ mas/yr

3.10.10 The Kerr–Ether Identity

3.10.11 Summary: The Ether's Gravitational Structure

	Schwarzschild (Thm 3.2)	Kerr (Thm 3.4)
Ether flow	Purely radial (potential flow)	Spiralling (potential + vortical)
Flow speed	$V = \sqrt{2GM/r}$	$V = c\sqrt{r_sr/\Sigma}$
Vorticity	Zero (irrotational)	Gravitomagnetic dipole (Eqs. 3.91a–b)
Spatial geometry	Flat $\mathbb{R}^3$	Oblate + gravitomagnetic twist
Horizon	$r_s$	$r_+ = r_s/2 + \sqrt{r_s^2/4 - a^2}$
Frame dragging	None	$\omega = 2GJ/(c^2r^3)$ (weak field)
Physical picture	Radial drain	Spiralling vortex
GEM sector	Gravitoelectric only	Gravitoelectric + gravitomagnetic

The Schwarzschild ether is the gravitational analog of electrostatics: an irrotational flow producing a conservative field. The Kerr ether is the analog of full electrodynamics: a flow with both potential and vortical components, producing both conservative and solenoidal fields. Theorem 3.2 captures the first; Theorem 3.4 captures both. Together, they account for the gravitational field of every astrophysical compact object.

3.11 Nonlinear Ether Field Equation: The Einstein Equation from Uniqueness

Section 3.9 derived the weak-field ether field (3.56) and identified the complete nonlinear extension as an open problem (Section 3.9.3). We now close this problem. Starting from the unit-lapse ether metric established by Theorems 3.2 and 3.4, we prove that the ether's nonlinear field equation is the Einstein equation — not by postulating it, but by deriving it from the geometric structure of the metric and the uniqueness theorems of Weinberg [160, 161], Deser [162], and Lovelock [163, 164].

The derivation has three stages: (i) the ADM decomposition of the ether metric, which expresses the Einstein tensor exactly in terms of the ether strain rate; (ii) recovery of the weak-field (3.56) and the sourced wave equation (Proposition 3.1) as limiting cases; (iii) the uniqueness argument, which guarantees that no other nonlinear extension is consistent with the established linearised dynamics.

3.11.1 ADM Decomposition of the Unit-Lapse Metric

The unit-lapse ether metric (Eq 3.17 with the PG constitutive relation) is:

$$ds^2 = -(c^2 - V^2)\,dT^2 + 2V_i\,dT\,dx^i + \delta_{ij}\,dx^i\,dx^j \tag{3.92}$$

where $V^i(T, x^j)$ is the ether velocity field and $V^2 = V_k V^k$.

Determinant. Writing the metric in matrix form with coordinates $(T, x^1, x^2, x^3)$:

$$g_{\mu\nu} = \begin{pmatrix} -(c^2 - V^2) & V_1 & V_2 & V_3 \\ V_1 & 1 & 0 & 0 \\ V_2 & 0 & 1 & 0 \\ V_3 & 0 & 0 & 1 \end{pmatrix}$$

Expanding along the first row by cofactors:

$$\det(g) = -(c^2 - V^2)\cdot 1 - V_1(V_1\cdot 1) - V_2(V_2\cdot 1) - V_3(V_3\cdot 1)$$ $$= -c^2 + V^2 - V_1^2 - V_2^2 - V_3^2 = -c^2 + V^2 - V^2 = -c^2 \tag{3.93}$$

The determinant is constant, independent of $V^i$.

Inverse metric. We seek $g^{\mu\nu}$ satisfying $g_{\mu\alpha}g^{\alpha\nu} = \delta^\nu_\mu$. Writing $g^{00} = \alpha$, $g^{0i} = \beta^i$, $g^{ij} = \gamma^{ij}$, the conditions $g_{\mu\alpha}g^{\alpha 0} = \delta^\mu_0$ give:

$$g_{00}\,g^{00} + g_{0k}\,g^{k0} = 1: \quad -(c^2 - V^2)\alpha + V_k\beta^k = 1 \tag{I.1}$$ $$g_{i0}\,g^{00} + g_{ik}\,g^{k0} = 0: \quad V_i\,\alpha + \beta^i = 0 \;\Rightarrow\; \beta^i = -V_i\,\alpha \tag{I.2}$$

Substituting (I.2) into (I.1): $-(c^2 - V^2)\alpha + V_k(-V_k\alpha) = 1$, giving $-(c^2 - V^2)\alpha - V^2\alpha = 1$, hence $-c^2\alpha = 1$:

$$g^{00} = \alpha = -\frac{1}{c^2} \tag{3.94a}$$ $$g^{0i} = \beta^i = -V_i\alpha = \frac{V^i}{c^2} \tag{3.94b}$$

The conditions $g_{\mu\alpha}g^{\alpha j} = \delta^\mu_j$ give:

$$g_{00}\,g^{0j} + g_{0k}\,g^{kj} = 0: \quad -(c^2-V^2)\frac{V^j}{c^2} + V_k\,\gamma^{kj} = 0 \tag{I.3}$$ $$g_{i0}\,g^{0j} + g_{ik}\,g^{kj} = \delta^j_i: \quad V_i\frac{V^j}{c^2} + \gamma^{ij} = \delta^j_i \tag{I.4}$$

From (I.4):

$$g^{ij} = \gamma^{ij} = \delta^{ij} - \frac{V^i V^j}{c^2} \tag{3.94c}$$

Verification of (I.3): $-(c^2-V^2)V^j/c^2 + V_k(\delta^{kj} - V^kV^j/c^2) = -V^j + V^2V^j/c^2 + V^j - V^2V^j/c^2 = 0$.

Collecting (3.94a–c):

$$g^{00} = -\frac{1}{c^2}, \qquad g^{0i} = \frac{V^i}{c^2}, \qquad g^{ij} = \delta^{ij} - \frac{V^i V^j}{c^2} \tag{3.94}$$

ADM identification. The ADM line element is:

$$ds^2 = -(N^2 - N_k N^k)\,dT^2 + 2N_i\,dT\,dx^i + \gamma_{ij}\,dx^i\,dx^j$$

Comparing term by term with (3.92):

$$g_{ij} = \delta_{ij} \;\Rightarrow\; \gamma_{ij} = \delta_{ij} \tag{3.95a}$$ $$g_{0i} = V_i \;\Rightarrow\; N_i = V_i, \quad N^i = \gamma^{ij}N_j = \delta^{ij}V_j = V^i \tag{3.95b}$$ $$g_{00} = -(N^2 - N_kN^k): \quad -(c^2 - V^2) = -(N^2 - V^2) \;\Rightarrow\; N = c \tag{3.95c}$$

Unit normal. The future-directed unit normal to $T = \text{const}$ has covariant components $n_\mu = (-N, 0, 0, 0) = (-c, 0, 0, 0)$. Its contravariant components are $n^\mu = g^{\mu\nu}n_\nu$:

$$n^0 = g^{00}n_0 + g^{0k}n_k = \left(-\frac{1}{c^2}\right)(-c) + \frac{V^k}{c^2}\cdot 0 = \frac{1}{c} \tag{3.96a}$$ $$n^i = g^{i0}n_0 + g^{ik}n_k = \frac{V^i}{c^2}(-c) + \left(\delta^{ik} - \frac{V^iV^k}{c^2}\right)\cdot 0 = -\frac{V^i}{c} \tag{3.96b}$$ $$n^\mu = \left(\frac{1}{c},\;-\frac{V^i}{c}\right) \tag{3.96}$$

Verification: $n_\mu n^\mu = n_0 n^0 + n_i n^i = (-c)(1/c) + 0\cdot(-V^i/c) = -1$.

Extrinsic curvature. The extrinsic curvature of the $T = \text{const}$ hypersurface is ([165], (4.63)):

$$K_{ij} = \frac{1}{2N}\!\left(D_i N_j + D_j N_i - \partial_T \gamma_{ij}\right) \tag{3.97}$$

where $D_i$ is the covariant derivative compatible with $\gamma_{ij}$. Evaluating each ingredient for the ether metric:

(a) $\gamma_{ij} = \delta_{ij}$ has vanishing Christoffel symbols ${}^{(3)}\Gamma^k_{ij} = 0$, so $D_i = \partial_i$.

(b) $\gamma_{ij} = \delta_{ij}$ is time-independent: $\partial_T\gamma_{ij} = 0$.

(c) $N = c$, $N_j = V_j$.

Substituting into (3.97):

$$K_{ij} = \frac{1}{2c}\left(\partial_i V_j + \partial_j V_i\right) = \frac{V_{i,j} + V_{j,i}}{2c} \tag{3.98}$$

Define the ether strain rate tensor $S_{ij} = (V_{i,j} + V_{j,i})/2$ (the symmetrised spatial velocity gradient). Then:

$$K_{ij} = \frac{S_{ij}}{c} \tag{3.99}$$

The trace is:

$$K = \delta^{ij}K_{ij} = \frac{1}{2c}\,\delta^{ij}(V_{i,j} + V_{j,i}) = \frac{1}{2c}(V_{i,i} + V_{j,j}) = \frac{1}{c}\,V_{k,k} = \frac{\text{div}\,\mathbf{V}}{c} \tag{3.100}$$

where we used $\delta^{ij}V_{i,j} = V_{k,k} = \text{div}\,\mathbf{V}$ and $\delta^{ij}V_{j,i} = V_{k,k}$ (relabelling $i \to k$).

3.11.2 The Hamiltonian Constraint

We derive $G_{\mu\nu}\,n^\mu n^\nu = \frac{1}{2}(K^2 - K_{ij}K^{ij})$ from the Gauss equation by explicit contraction.

The Gauss equation. For any spacelike hypersurface $\Sigma$ with unit normal $n^\mu$, induced metric $\gamma_{ij}$, and extrinsic curvature $K_{ij} = -\gamma^\mu{}_i\gamma^\nu{}_j\nabla_\mu n_\nu$, the spatial projection of the 4D Riemann tensor satisfies ([165], Section 2.4):

$${}^{(4)}R_{abcd}\,\gamma^a{}_i\,\gamma^b{}_j\,\gamma^c{}_k\,\gamma^d{}_l = {}^{(3)}R_{ijkl} + K_{ik}\,K_{jl} - K_{il}\,K_{jk} \tag{3.101}$$

This identity is derived by expressing the 3D Riemann tensor through the commutator $(D_k D_l - D_l D_k)W_j$ of 3D covariant derivatives acting on an arbitrary spatial vector $W_j$, then relating $D_k$ to the 4D derivative $\nabla_k$ via the Gauss-Weingarten relation $\nabla_a W_b = D_a W_b - n_b K_a{}^m W_m$ (for $n^\mu W_\mu = 0$). The cross terms from two applications of this decomposition produce $K_{ik}K_{jl} - K_{il}K_{jk}$, while the 4D commutator produces ${}^{(4)}R_{abcd}$ projected onto $\Sigma$. The detailed proof occupies Section 2.4 of [165]; we use only the result and the fact that for $\gamma_{ij} = \delta_{ij}$ (flat spatial sections), ${}^{(3)}R_{ijkl} = 0$:

$${}^{(4)}R_{abcd}\,\gamma^a{}_i\,\gamma^b{}_j\,\gamma^c{}_k\,\gamma^d{}_l = K_{ik}\,K_{jl} - K_{il}\,K_{jk} \tag{3.102}$$

This has been verified numerically: all six independent non-zero components match the $K_{ik}K_{jl} - K_{il}K_{jk}$ prediction for the Schwarzschild flow.

First contraction. Contract (3.102) on indices $i$ and $k$ by summing with $\gamma^{ik} = \delta^{ik}$. On the left, we must account for the fact that the Ricci tensor involves $g^{ac}$, not $\gamma^{ac}$. Using the completeness relation $g^{ac} = \gamma^{ac} - n^a n^c$:

$$\gamma^{ik}\,{}^{(4)}R_{abcd}\,\gamma^a{}_i\,\gamma^b{}_j\,\gamma^c{}_k\,\gamma^d{}_l = g^{ac}\,R_{abcd}\,\gamma^b{}_j\,\gamma^d{}_l + n^a n^c\,R_{abcd}\,\gamma^b{}_j\,\gamma^d{}_l$$

The first term: $g^{ac}R_{abcd} = R_{bd}$ (definition of the Ricci tensor), so it equals $R_{bd}\,\gamma^b{}_j\,\gamma^d{}_l$, the spatial projection of the Ricci tensor. The second term: $n^a n^c R_{abcd}\,\gamma^b{}_j\,\gamma^d{}_l = R_{\perp j\perp l}$, the "electric" part of the Riemann tensor. On the right:

$$\gamma^{ik}(K_{ik}K_{jl} - K_{il}K_{jk}) = K\,K_{jl} - K^k{}_l\,K_{jk}$$

where $K = \gamma^{ik}K_{ik}$. Therefore:

$$R_{\mu\nu}\,\gamma^\mu{}_j\,\gamma^\nu{}_l + R_{\perp j\perp l} = K\,K_{jl} - K_{jk}\,K^k{}_l \tag{3.103}$$

Second contraction. Contract (3.103) on $j$ and $l$ using $\gamma^{jl} = \delta^{jl}$.

Left side, first term: $\gamma^{jl}R_{\mu\nu}\,\gamma^\mu{}_j\,\gamma^\nu{}_l = R_{\mu\nu}\,\gamma^{\mu\nu}$. Using $\gamma^{\mu\nu} = g^{\mu\nu} + n^\mu n^\nu$:

$$R_{\mu\nu}\,\gamma^{\mu\nu} = R_{\mu\nu}\,g^{\mu\nu} + R_{\mu\nu}\,n^\mu n^\nu = R + R_{\perp\perp} \tag{C.1}$$

Left side, second term: $\gamma^{jl}R_{\perp j\perp l} = R_{\mu\nu\alpha\beta}\,n^\mu\,\gamma^{\nu\beta}\,n^\alpha$. Using $\gamma^{\nu\beta} = g^{\nu\beta} + n^\nu n^\beta$:

$$R_{\mu\nu\alpha\beta}\,n^\mu\,\gamma^{\nu\beta}\,n^\alpha = R_{\mu\nu\alpha\beta}\,n^\mu\,g^{\nu\beta}\,n^\alpha + R_{\mu\nu\alpha\beta}\,n^\mu\,n^\nu\,n^\beta\,n^\alpha$$

The first part is $R_{\mu\alpha}\,n^\mu n^\alpha = R_{\perp\perp}$ (contraction of Riemann to Ricci). The second part vanishes: $R_{\mu\nu\alpha\beta}$ is antisymmetric in $(\mu,\nu)$ while $n^\mu n^\nu$ is symmetric, so $R_{\mu\nu\alpha\beta}\,n^\mu n^\nu = 0$ for every fixed $(\alpha,\beta)$. Therefore:

$$\gamma^{jl}R_{\perp j\perp l} = R_{\perp\perp} \tag{C.2}$$

Left side total: $(R + R_{\perp\perp}) + R_{\perp\perp} = R + 2R_{\perp\perp}$.

Right side: $\gamma^{jl}(K\,K_{jl} - K_{jk}K^k{}_l) = K\cdot K - K_{jk}K^{jk} = K^2 - K_{ij}K^{ij}$.

Result:

$$R + 2R_{\perp\perp} = K^2 - K_{ij}\,K^{ij} \tag{3.104}$$

Assembly. The Einstein tensor projected along the normal is:

$$G_{\perp\perp} = G_{\mu\nu}\,n^\mu n^\nu = R_{\mu\nu}\,n^\mu n^\nu - \frac{1}{2}\,g_{\mu\nu}\,n^\mu n^\nu\,R \tag{C.3}$$

Since $g_{\mu\nu}\,n^\mu n^\nu = n_\mu n^\mu = -1$:

$$G_{\perp\perp} = R_{\perp\perp} + \frac{1}{2}R \tag{C.4}$$

From (3.104): $R_{\perp\perp} = \frac{1}{2}(K^2 - K_{ij}K^{ij}) - \frac{1}{2}R$. Substituting into (C.4):

$$G_{\perp\perp} = \frac{1}{2}(K^2 - K_{ij}K^{ij}) - \frac{1}{2}R + \frac{1}{2}R = \frac{1}{2}(K^2 - K_{ij}\,K^{ij}) \tag{3.105}$$

The Ricci scalar cancels exactly. This is the Hamiltonian constraint, exact for all values of $V^i/c$.

Substitution of ether quantities. Using (3.99) and (3.100):

$$K^2 = \frac{(\text{div}\,\mathbf{V})^2}{c^2} \tag{3.106a}$$ $$K_{ij}\,K^{ij} = \frac{1}{c^2}\sum_{ij}S_{ij}^2 \tag{3.106b}$$

(3.106b) requires the identity $\sum_{ij}K_{ij}^2 = (1/(4c^2))\sum_{ij}(V_{i,j}+V_{j,i})^2 = (1/c^2)\sum_{ij}S_{ij}^2$, which we now derive. Decompose $V_{i,j} = S_{ij} + A_{ij}$ where $A_{ij} = (V_{i,j} - V_{j,i})/2$ is antisymmetric ($A_{ij} = -A_{ji}$). Then:

$$V_{i,j}^2 = (S_{ij}+A_{ij})^2 = S_{ij}^2 + 2S_{ij}A_{ij} + A_{ij}^2$$ $$V_{i,j}V_{j,i} = (S_{ij}+A_{ij})(S_{ji}+A_{ji}) = (S_{ij}+A_{ij})(S_{ij}-A_{ij}) = S_{ij}^2 - A_{ij}^2$$

Adding these two expressions:

$$V_{i,j}^2 + V_{i,j}V_{j,i} = 2S_{ij}^2 + 2S_{ij}A_{ij}$$

Summing over $i,j$: the cross term vanishes, $\sum_{ij}S_{ij}A_{ij} = 0$, because relabelling $i \leftrightarrow j$ gives $\sum S_{ji}A_{ji} = \sum S_{ij}(-A_{ij}) = -\sum S_{ij}A_{ij}$, forcing the sum to zero. Therefore:

$$\sum_{ij}\!\left(V_{i,j}^2 + V_{i,j}V_{j,i}\right) = 2\sum_{ij}S_{ij}^2 \tag{3.107}$$

Now expand $(V_{i,j}+V_{j,i})^2 = V_{i,j}^2 + 2V_{i,j}V_{j,i} + V_{j,i}^2$. Summing and using $\sum V_{j,i}^2 = \sum V_{i,j}^2$ (relabel $i \leftrightarrow j$):

$$\sum_{ij}(V_{i,j}+V_{j,i})^2 = 2\sum_{ij}V_{i,j}^2 + 2\sum_{ij}V_{i,j}V_{j,i} = 2\sum_{ij}(V_{i,j}^2 + V_{i,j}V_{j,i}) = 4\sum_{ij}S_{ij}^2$$

by (3.107). Hence $K_{ij}K^{ij} = (1/(4c^2))\cdot 4\sum S_{ij}^2/1 = (1/c^2)\sum S_{ij}^2$, confirming (3.106b).

The Hamiltonian constraint (3.105) becomes:

$$G_{\perp\perp} = \frac{1}{2c^2}\!\left[(\text{div}\,\mathbf{V})^2 - \sum_{ij}S_{ij}^2\right] \tag{3.108}$$

The vorticity $A_{ij}$ does not appear. For the sourced Einstein equation $G_{\mu\nu} = (8\pi G/c^4)T_{\mu\nu}$:

$$\frac{1}{2c^2}\!\left[(\text{div}\,\mathbf{V})^2 - \sum_{ij}S_{ij}^2\right] = \frac{8\pi G}{c^2}\,\rho_E \tag{3.109}$$

where $\rho_E = T_{\mu\nu}\,n^\mu n^\nu$ is the energy density measured by the Eulerian observer.

3.11.3 The Momentum Constraint

The Codazzi–Mainardi equation gives the mixed projection of the 4D Riemann tensor — one normal index and three spatial indices — in terms of the extrinsic curvature. We derive it and then specialise to the ether metric.

The Codazzi–Mainardi identity. Consider the covariant derivative of $K_{ij}$ along $\Sigma$. By definition, $K_{ij} = -\gamma^\mu{}_i\,\gamma^\nu{}_j\,\nabla_\mu n_\nu$. Taking the 3D covariant derivative $D_k$ and antisymmetrising on $j, k$ yields ([165], Section 2.5):

$$D_k K_{ij} - D_j K_{ik} = R_{\mu\nu\alpha\beta}\,\gamma^\mu{}_i\,n^\nu\,\gamma^\alpha{}_j\,\gamma^\beta{}_k \tag{3.110a}$$

This identity follows from expanding $D_k K_{ij}$ using the Gauss–Weingarten relation $\nabla_a e^\mu_b = D_a e^\mu_b + K_{ab}n^\mu$ (where $e^\mu_i = \gamma^\mu{}_i$ are the spatial basis vectors) and recognising that the antisymmetric part of the resulting expression reproduces the 4D Riemann tensor contracted with one normal vector. The derivation is given in [165], Section 2.5; we use only the result and the fact that for $\gamma_{ij} = \delta_{ij}$, the 3D covariant derivative reduces to the partial derivative ($D_k = \partial_k$).

Contraction to the momentum constraint. Contract (3.110a) on $i$ and $k$ using $\gamma^{ik} = \delta^{ik}$. On the left: $\delta^{ik}(D_k K_{ij} - D_j K_{ik}) = D_k K^k{}_j - D_j K$. On the right: $\delta^{ik}R_{\mu\nu\alpha\beta}\,\gamma^\mu{}_i\,n^\nu\,\gamma^\alpha{}_j\,\gamma^\beta{}_k = R_{\mu\nu\alpha\beta}\,\gamma^{\mu\beta}\,n^\nu\,\gamma^\alpha{}_j$. Using the completeness relation $\gamma^{\mu\beta} = g^{\mu\beta} + n^\mu n^\beta$ and the antisymmetry $R_{\mu\nu\alpha\beta}n^\mu n^\beta = 0$ (symmetric in $\mu, \beta$ contracted with antisymmetric Riemann), the right side reduces to $g^{\mu\beta}R_{\mu\nu\alpha\beta}\,n^\nu\,\gamma^\alpha{}_j = R_{\nu\alpha}\,n^\nu\,\gamma^\alpha{}_j$. The Einstein tensor projection is $G_{\mu\nu}\,n^\mu\,\gamma^\nu{}_j = R_{\mu\nu}\,n^\mu\,\gamma^\nu{}_j - \frac{1}{2}g_{\mu\nu}\,n^\mu\,\gamma^\nu{}_j\,R = R_{\mu\nu}\,n^\mu\,\gamma^\nu{}_j$ (since $g_{\mu\nu}\,n^\mu\,\gamma^\nu{}_j = n_\nu\,\gamma^\nu{}_j = 0$ by orthogonality of $n^\mu$ and $\gamma^\nu{}_j$). Therefore:

$$D_j(K^j{}_i - \delta^j{}_i\,K) = G_{\mu\nu}\,n^\mu\,\gamma^\nu{}_i = \frac{8\pi G}{c^4}\,J_i \tag{3.110}$$

where $J_i = -T_{\mu i}\,n^\mu$ is the momentum density measured by the Eulerian observer. The last equality follows from the Einstein equation. This is the momentum constraint, exact for all values of $V^i/c$.

Evaluation for the ether metric. For flat spatial metric ($D_j = \partial_j$), substituting $K^j{}_i = \gamma^{jk}K_{ki} = \delta^{jk}K_{ki} = K_{ji}$ and $K = (\text{div}\,V)/c$:

$$\partial_j K_{ji} - \partial_i K = \frac{8\pi G}{c^4}\,J_i$$ $$\frac{1}{2c}\,\partial_j(V_{j,i} + V_{i,j}) - \frac{1}{c}\,\partial_i(\text{div}\,V) = \frac{8\pi G}{c^4}\,J_i \tag{3.111}$$

Expanding the left side. The first term: $\partial_j V_{j,i} = \partial_i(\partial_j V_j) = \partial_i(\text{div}\,V)$ (commuting partial derivatives, valid for $C^2$ fields). The second term: $\partial_j V_{i,j} = \partial_j\partial_j V_i = \nabla^2 V_i$. So the left side of (3.111) is:

$$\frac{1}{2c}\!\left[\partial_i(\text{div}\,V) + \nabla^2 V_i\right] - \frac{1}{c}\,\partial_i(\text{div}\,V) = \frac{1}{2c}\!\left[\nabla^2 V_i - \partial_i(\text{div}\,V)\right]$$

In vacuum ($J_i = 0$):

$$\nabla^2 V_i = \partial_i(\text{div}\,V) \tag{3.112}$$

By the vector identity $\nabla^2\mathbf{V} = \nabla(\text{div}\,\mathbf{V}) - \nabla\times(\nabla\times\mathbf{V})$, (3.112) is equivalent to:

$$\nabla\times(\nabla\times\mathbf{V}) = 0 \tag{3.113}$$

For irrotational flow ($\mathbf{V} = -\nabla\Psi$): $\nabla\times\mathbf{V} = 0$ identically, so (3.113) is automatically satisfied.

3.11.4 Schwarzschild Verification

The Schwarzschild free-fall velocity (Eq 3.22) in Cartesian coordinates is $V_i = -\sqrt{2GM}\,x_i\,r^{-3/2}$, with $r = (x_k x_k)^{1/2}$.

Velocity gradient. By the product rule:

$$V_{i,j} = \frac{\partial}{\partial x_j}\!\left(-\sqrt{2GM}\,x_i\,r^{-3/2}\right) = -\sqrt{2GM}\!\left[\frac{\partial x_i}{\partial x_j}\,r^{-3/2} + x_i\,\frac{\partial(r^{-3/2})}{\partial x_j}\right] \tag{3.114}$$

The first factor: $\partial x_i/\partial x_j = \delta_{ij}$. The second factor uses $\partial r/\partial x_j = x_j/r$:

$$\frac{\partial(r^{-3/2})}{\partial x_j} = -\frac{3}{2}\,r^{-5/2}\,\frac{\partial r}{\partial x_j} = -\frac{3}{2}\,\frac{x_j}{r}\,r^{-5/2} = -\frac{3}{2}\,\frac{x_j}{r^{7/2}} \tag{3.115}$$

Substituting into (3.114):

$$V_{i,j} = -\sqrt{2GM}\!\left[\frac{\delta_{ij}}{r^{3/2}} - \frac{3}{2}\,\frac{x_i\,x_j}{r^{7/2}}\right] \tag{3.116}$$

Symmetry check. $V_{i,j} = V_{j,i}$ (exchanging $i \leftrightarrow j$ in (3.116) leaves the expression unchanged), confirming the flow is irrotational: $A_{ij} = 0$, $S_{ij} = V_{i,j}$.

Divergence. Setting $i = j$ in (3.116) and summing ($\delta_{ii} = 3$, $x_ix_i = r^2$):

$$\text{div}\,\mathbf{V} = V_{i,i} = -\sqrt{2GM}\!\left[\frac{3}{r^{3/2}} - \frac{3}{2}\,\frac{r^2}{r^{7/2}}\right] = -\sqrt{2GM}\!\left[\frac{3}{r^{3/2}} - \frac{3}{2r^{3/2}}\right]$$ $$= -\frac{3\sqrt{2GM}}{2r^{3/2}} \tag{3.117}$$ $$(\text{div}\,\mathbf{V})^2 = \frac{9\cdot 2GM}{4r^3} = \frac{9GM}{2r^3} \tag{3.118}$$

Sum of squared gradients. Setting $A = \sqrt{2GM}$, $a = r^{-3/2}$, $b = (3/2)\,r^{-7/2}$:

$$\sum_{ij}V_{i,j}^2 = A^2\sum_{ij}(a\,\delta_{ij} - b\,x_ix_j)^2 \tag{3.119}$$

Expanding the square:

$$(a\,\delta_{ij} - b\,x_ix_j)^2 = a^2\delta_{ij}^2 - 2ab\,\delta_{ij}\,x_ix_j + b^2\,x_i^2x_j^2$$

The three sums are evaluated using $\sum_{ij}\delta_{ij}^2 = \sum_i 1 = 3$ (since $\delta_{ij}^2 = \delta_{ij}$, which equals 1 only when $i = j$); $\sum_{ij}\delta_{ij}\,x_ix_j = \sum_i x_i^2 = r^2$; and $\sum_{ij}x_i^2 x_j^2 = (\sum_i x_i^2)^2 = r^4$ (since $\sum_{ij}x_i^2x_j^2 = \sum_i x_i^2\sum_j x_j^2$). Therefore:

$$\sum_{ij}(a\,\delta_{ij} - b\,x_ix_j)^2 = 3a^2 - 2ab\,r^2 + b^2\,r^4 \tag{3.120}$$

Computing each term:

$$3a^2 = 3r^{-3} \tag{3.121a}$$ $$2ab\,r^2 = 2\cdot r^{-3/2}\cdot\frac{3}{2}\,r^{-7/2}\cdot r^2 = 3\,r^{-3/2-7/2+2} = 3\,r^{-3} \tag{3.121b}$$ $$b^2r^4 = \frac{9}{4}\,r^{-7}\cdot r^4 = \frac{9}{4}\,r^{-3} \tag{3.121c}$$

Substituting (3.121a–c) into (3.120):

$$\sum_{ij}(a\,\delta_{ij} - b\,x_ix_j)^2 = 3r^{-3} - 3r^{-3} + \frac{9}{4}\,r^{-3} = \frac{9}{4}\,r^{-3}$$ $$\sum_{ij}V_{i,j}^2 = A^2\cdot\frac{9}{4}\,r^{-3} = 2GM\cdot\frac{9}{4r^3} = \frac{9GM}{2r^3} \tag{3.122}$$

Constraint verification. Comparing (3.118) and (3.122):

$$(\text{div}\,\mathbf{V})^2 = \sum_{ij}V_{i,j}^2 = \frac{9GM}{2r^3} \tag{3.123}$$

The vacuum Hamiltonian constraint $K^2 = K_{ij}K^{ij}$ (Eq 3.105 with $G_{\perp\perp} = 0$) is satisfied identically. $\square$

Numerical verification. Tested at 11 points: 6 structured (on-axis and off-axis, $5 \leq r/r_s \leq 100$) and 5 random ($r$, $\theta$, $\phi$ uniformly sampled). All satisfy $|K^2 - K_{ij}K^{ij}|/K^2 < 10^{-15}$. The curl $|\nabla\times\mathbf{V}|$ vanishes to machine precision at all points, confirming $A_{ij} = 0$.

3.11.5 The Geodesic-Euler Correspondence and Self-Coupling

The geodesic equation for a particle with four-velocity $u^\mu$ is $du^\alpha/d\tau + \Gamma^\alpha_{\mu\nu}u^\mu u^\nu = 0$. For a particle instantaneously at rest in the coordinate frame ($u^i \approx 0$, $u^0 \approx c\,dT/d\tau$), the spatial components reduce to:

$$\frac{d^2 x^k}{dT^2} = -\Gamma^k_{00} \tag{3.124}$$

at leading order in the particle velocity.

Derivation of $\Gamma^k_{00}$. From the Christoffel formula:

$$\Gamma^k_{00} = \frac{1}{2}\,g^{k\beta}\!\left(2\,\partial_0 g_{\beta 0} - \partial_\beta g_{00}\right) \tag{3.125}$$

For the static ether metric ($\partial_0 g_{\mu\nu} = 0$), only the $-\partial_\beta g_{00}$ term survives. Now $g_{00} = -(c^2 - V^2)$, so $\partial_\beta g_{00} = \partial_\beta(V^2) = 2V_m V_{m,\beta}$. For $\beta = 0$: vanishes (static). For $\beta = j$ (spatial): $\partial_j g_{00} = 2V_m V_{m,j} \equiv 2S_j$, where $S_j = V_m V_{m,j} = \frac{1}{2}\partial_j(V^2)$. Therefore:

$$\Gamma^k_{00} = -\frac{1}{2}\,g^{kj}\,\partial_j g_{00} = -g^{kj}\,S_j \tag{3.126}$$

Substituting $g^{kj} = \delta^{kj} - V^kV^j/c^2$ from (3.94c):

$$\Gamma^k_{00} = -\left(\delta^{kj} - \frac{V^kV^j}{c^2}\right)S_j = -S^k + \frac{V^k V^j S_j}{c^2} \tag{3.127}$$

Connection to the Newtonian potential. The Bernoulli (3.55) gives $V^2/2 = -\Phi$, hence:

$$S_k = \frac{1}{2}\,\partial_k(V^2) = -\partial_k\Phi \tag{3.128}$$

The Newtonian acceleration is $a^k_N = -\partial^k\Phi = S^k$. For the quadratic term: $V^jS_j = V^j(-\partial_j\Phi) = -(\mathbf{V}\cdot\nabla)\Phi$. The geodesic acceleration (3.124) becomes:

$$-\Gamma^k_{00} = S^k - \frac{V^k V^j S_j}{c^2} = -\partial^k\Phi + \frac{V^k(\mathbf{V}\cdot\nabla\Phi)}{c^2} \tag{3.129}$$

Evaluation for Schwarzschild at the point $(r, 0, 0)$. Here $V^i = (-\sqrt{2GM/r},\;0,\;0)$ and $\Phi = -GM/r$. The components of $S_k$ are:

$$S_1 = -\partial_1\Phi = -\frac{\partial}{\partial x}\!\left(-\frac{GM}{r}\right) = -\frac{GM\,x}{r^3}\bigg|_{(r,0,0)} = -\frac{GM}{r^2} \tag{3.130a}$$ $$S_2 = S_3 = 0 \tag{3.130b}$$

The scalar product $V^jS_j$ (only $j = 1$ contributes):

$$V^j S_j = V^1 S_1 = \left(-\sqrt{\frac{2GM}{r}}\right)\!\left(-\frac{GM}{r^2}\right) = \frac{GM\sqrt{2GM}}{r^{5/2}} \tag{3.131}$$

The quadratic correction in (3.127):

$$\frac{V^k V^j S_j}{c^2}\bigg|_{k=1} = \frac{V^1}{c^2}\cdot V^jS_j = \frac{1}{c^2}\!\left(-\sqrt{\frac{2GM}{r}}\right)\!\frac{GM\sqrt{2GM}}{r^{5/2}} = -\frac{2G^2M^2}{c^2\,r^3} \tag{3.132}$$

Substituting (3.130a) and (3.132) into (3.129) for $k = 1$:

$$-\Gamma^1_{00} = -\frac{GM}{r^2} - \left(-\frac{2G^2M^2}{c^2\,r^3}\right) = -\frac{GM}{r^2} + \frac{2G^2M^2}{c^2\,r^3} \tag{3.133}$$

Factoring:

$$-\Gamma^1_{00} = -\frac{GM}{r^2}\!\left(1 - \frac{2GM}{c^2\,r}\right) = -\frac{GM}{r^2}\!\left(1 - \frac{r_s}{r}\right) \tag{3.134}$$

where $r_s = 2GM/c^2$ is the Schwarzschild radius. The geodesic acceleration is $(1 - r_s/r)$ times the Newtonian value. The correction:

$$\frac{\delta a}{a_N} = -\frac{r_s}{r} = -\frac{V^2}{c^2} = -\frac{2GM}{c^2\,r} \tag{3.135}$$

is the first post-Newtonian (1PN) self-coupling term.

Physical origin of the correction. The factor $(1 - r_s/r)$ arises from the $-V^kV^j/c^2$ correction to the inverse spatial metric $g^{kj} = \delta^{kj} - V^kV^j/c^2$ (Eq 3.94c). The gravitational field's own energy, encoded in the metric through $V^2/c^2$, reduces the acceleration. This is the concrete realisation of the self-coupling that Deser [162] identified abstractly: the ether strain rate tensor $K_{ij}$ encodes both the gravitational field and its self-energy.

Numerical verification. The ratio $|\delta a|/|a_N| \div (r_s/r)$ has been evaluated at five test points:

$r/r_s$	$\|\delta a\|/\|a_N\|$	$r_s/r$	Ratio
54.8	$1.826 \times 10^{-2}$	$1.826 \times 10^{-2}$	1.0000
100.0	$1.000 \times 10^{-2}$	$1.000 \times 10^{-2}$	1.0000
229.1	$4.364 \times 10^{-3}$	$4.364 \times 10^{-3}$	1.0000
500.0	$2.000 \times 10^{-3}$	$2.000 \times 10^{-3}$	1.0000
1063.0	$9.407 \times 10^{-4}$	$9.407 \times 10^{-4}$	1.0000

All ratios are unity to machine precision, confirming (3.134)–(3.135).

Self-consistency of the Hamiltonian constraint. The Hamiltonian constraint (3.109) is exact and nonlinear — it contains the gravitational self-energy at all orders through the quadratic structure of $K_{ij}K^{ij}$. For Schwarzschild, the identity $K^2 = K_{ij}K^{ij}$ (Eq 3.123) means the expansion energy ($\propto K^2$) and the shear energy ($\propto K_{ij}K^{ij}$) of the ether flow balance exactly, producing zero net gravitational energy in the vacuum exterior. Inside a matter distribution, this balance is broken by the source $\rho_E$, and the mismatch determines the gravitational field self-consistently.

3.11.6 Recovery of the Weak-Field Equation

We now show that the Hamiltonian constraint (3.109) reduces to the Poisson equation $\nabla^2\Phi = 4\pi G\rho_m$ in the weak-field limit. The derivation is exact until the final step, where a controlled approximation is made.

Step 1: Differentiate the Bernoulli equation. From $|\nabla\Psi|^2/2 = -\Phi$ (Eq 3.55), where $\mathbf{V} = -\nabla\Psi$, differentiate with respect to $x_i$:

$$\frac{1}{2}\,\partial_i(\Psi_{,k}\Psi_{,k}) = -\Phi_{,i}$$

Expanding the left side by the product rule: $\frac{1}{2}\cdot 2\Psi_{,k}\Psi_{,ki} = \Psi_{,k}\Psi_{,ki}$. Therefore:

$$\Psi_{,k}\,\Psi_{,ki} = -\Phi_{,i} \tag{3.136}$$

This is exact. Numerical evaluation for the Schwarzschild flow confirms LHS/RHS $= 1$ for all three components $i = 1, 2, 3$ at five test radii spanning $5 \leq r/r_s \leq 550$.

Step 2: Differentiate again and trace. Differentiate (3.136) with respect to $x_j$:

$$\frac{\partial}{\partial x_j}(\Psi_{,k}\,\Psi_{,ki}) = -\Phi_{,ij}$$

Expanding the left side by the product rule: $\Psi_{,kj}\Psi_{,ki} + \Psi_{,k}\Psi_{,kij} = -\Phi_{,ij}$. Set $i = j$ and sum over $i$:

$$\sum_{ki}\Psi_{,ki}^2 + \Psi_{,k}\,\Psi_{,kii} = -\Phi_{,ii} \tag{3.137}$$

On the left: $\sum_{ki}\Psi_{,ki}^2 = \sum_{ij}\Psi_{,ij}^2$ (relabelling $k \to i$, $i \to j$). And $\Psi_{,kii} = \partial_i\partial_i\Psi_{,k} = (\nabla^2\Psi)_{,k}$ (commuting partial derivatives, valid for $C^3$ fields). On the right: $\Phi_{,ii} = \nabla^2\Phi$. Therefore:

$$\sum_{ij}\Psi_{,ij}^2 + \Psi_{,k}\,(\nabla^2\Psi)_{,k} = -\nabla^2\Phi \tag{3.138}$$

Alternative derivation of (3.138). Take the Laplacian of $|\nabla\Psi|^2/2 = -\Phi$ directly: $\frac{1}{2}\nabla^2(\Psi_{,k}\Psi_{,k}) = -\nabla^2\Phi$. Expand: first derivative $\partial_j(\Psi_{,k}\Psi_{,k}) = 2\Psi_{,k}\Psi_{,kj}$; second derivative $\partial_j(2\Psi_{,k}\Psi_{,kj}) = 2(\Psi_{,kj}\Psi_{,kj} + \Psi_{,k}\Psi_{,kjj}) = 2(\sum_{jk}\Psi_{,kj}^2 + \Psi_{,k}(\nabla^2\Psi)_{,k})$. Dividing by 2 reproduces (3.138). $\square$

Step 3: Substitute the Hamiltonian constraint. For irrotational flow ($S_{ij} = V_{i,j} = -\Psi_{,ij}$), the Hamiltonian constraint (3.109) gives:

$$(\text{div}\,V)^2 - \sum_{ij}V_{i,j}^2 = 16\pi G\,\rho_E$$

Since $\text{div}\,V = -\nabla^2\Psi$ and $V_{i,j} = -\Psi_{,ij}$:

$$(\nabla^2\Psi)^2 - \sum_{ij}\Psi_{,ij}^2 = 16\pi G\,\rho_E \tag{3.139}$$

Solving for $\sum\Psi_{,ij}^2$:

$$\sum_{ij}\Psi_{,ij}^2 = (\nabla^2\Psi)^2 - 16\pi G\,\rho_E \tag{3.140}$$

Substitute (3.140) into (3.138):

$$(\nabla^2\Psi)^2 - 16\pi G\,\rho_E + \Psi_{,k}\,(\nabla^2\Psi)_{,k} = -\nabla^2\Phi$$

Rearranging:

$$\nabla^2\Phi = 16\pi G\,\rho_E - (\nabla^2\Psi)^2 - \Psi_{,k}\,(\nabla^2\Psi)_{,k} \tag{3.141}$$

This is exact.

Step 4: The divergence identity. The two nonlinear terms in (3.141) combine into a total divergence. Expand $\partial_k[(\nabla^2\Psi)\Psi_{,k}]$ by the product rule:

$$\partial_k\!\left[(\nabla^2\Psi)\,\Psi_{,k}\right] = (\nabla^2\Psi)_{,k}\,\Psi_{,k} + (\nabla^2\Psi)\,\Psi_{,kk} \tag{3.142}$$

The first term is $\Psi_{,k}(\nabla^2\Psi)_{,k}$. The second term: $\Psi_{,kk} = \sum_k\partial_k^2\Psi = \nabla^2\Psi$, so it equals $(\nabla^2\Psi)^2$. Therefore:

$$(\nabla^2\Psi)^2 + \Psi_{,k}\,(\nabla^2\Psi)_{,k} = \text{div}\!\left[(\nabla^2\Psi)\,\nabla\Psi\right] \tag{3.143}$$

Numerical evaluation at three test points confirms LHS/RHS $= 1$ to finite-difference precision.

Step 5: The exact field equation for $\Phi$. Substituting (3.143) into (3.141):

$$\nabla^2\Phi = 16\pi G\,\rho_E - \text{div}\!\left[(\nabla^2\Psi)\,\nabla\Psi\right] \tag{3.144}$$

This is the exact, nonlinear field equation for the Newtonian potential: the Poisson equation with a post-Newtonian correction that is a total divergence. Every step from (3.136) to (3.144) is exact.

Step 6: Integral verification for a point mass. Integrate (3.144) over a sphere of radius $R$ enclosing the source. By the divergence theorem:

$$\oint_S \nabla\Phi\cdot d\mathbf{S} = 16\pi G\,M_E - \oint_S (\nabla^2\Psi)(\nabla\Psi\cdot d\mathbf{S}) \tag{3.145}$$

where $M_E = \int\rho_E\,dV$.

Left side. For a point mass $M$: $\Phi = -GM/r$, $\partial_r\Phi = GM/r^2$:

$$\oint_S \nabla\Phi\cdot d\mathbf{S} = \frac{GM}{R^2}\cdot 4\pi R^2 = 4\pi GM \tag{3.146a}$$

Right side, surface integral. For the Schwarzschild exterior, $\Psi = 2\sqrt{2GM\,r}$:

$$\Psi_{,r} = \frac{d}{dr}\!\left(2\sqrt{2GM}\,r^{1/2}\right) = \sqrt{2GM}\cdot r^{-1/2} = \frac{\sqrt{2GM}}{\sqrt{r}} \tag{3.146b}$$ $$\Psi_{,rr} = -\frac{\sqrt{2GM}}{2}\cdot r^{-3/2} = -\frac{\sqrt{2GM}}{2r^{3/2}} \tag{3.146c}$$ $$\nabla^2\Psi = \Psi_{,rr} + \frac{2}{r}\Psi_{,r} = -\frac{\sqrt{2GM}}{2r^{3/2}} + \frac{2\sqrt{2GM}}{r^{3/2}} = \frac{3\sqrt{2GM}}{2r^{3/2}} \tag{3.146d}$$

The integrand on the surface at radius $R$:

$$(\nabla^2\Psi)(\Psi_{,r}) = \frac{3\sqrt{2GM}}{2R^{3/2}}\cdot\frac{\sqrt{2GM}}{\sqrt{R}} = \frac{3\cdot 2GM}{2R^2} = \frac{3GM}{R^2} \tag{3.146e}$$ $$\oint_S (\nabla^2\Psi)(\nabla\Psi\cdot d\mathbf{S}) = \frac{3GM}{R^2}\cdot 4\pi R^2 = 12\pi GM \tag{3.146f}$$

Verification of (3.145) with $M_E = M$:

$$\text{LHS} = 4\pi GM, \qquad \text{RHS} = 16\pi GM - 12\pi GM = 4\pi GM \tag{3.147}$$

The identity is satisfied. $\square$

Step 7: The weak-field limit. Define the post-Newtonian parameter $\varepsilon = GM/(Rc^2) = r_s/(2R)$, where $R$ is the characteristic scale. In (3.144), both the source term $16\pi G\rho_E$ and the divergence term are $O(G\rho)$ — the same order. The Poisson equation does not emerge by dropping the divergence term; it emerges from a precise cancellation between them.

The cancellation mechanism. The ADM source $16\pi G\rho_E$ exceeds the Newtonian source $4\pi G\rho_m$ by a factor of four. The divergence term $\text{div}[(\nabla^2\Psi)\nabla\Psi]$ encodes the gravitational field's nonlinear self-coupling and absorbs the excess. The integral verification (Eq 3.147) demonstrates this explicitly for the point mass: the factor decomposes as $16\pi - 12\pi = 4\pi$, where $12\pi GM$ is the integrated self-coupling contribution.

For a general static source with $\rho_E \to \rho_m$ in the weak-field limit (pressureless matter at rest), we extract the Newtonian equation by the following argument. Define $\Phi_N$ as the Newtonian potential satisfying $\nabla^2\Phi_N = 4\pi G\rho_m$, and write $\Phi = \Phi_N + \delta\Phi$ with $|\delta\Phi| \ll |\Phi_N|$. The velocity potential $\Psi$ satisfies $|\nabla\Psi|^2/2 = -\Phi$ and can be expanded as $\Psi = \Psi_0 + \delta\Psi$, where $\Psi_0$ solves the leading-order (Newtonian) constraint. The divergence term in (3.144), evaluated at leading order, contributes:

$$\text{div}[(\nabla^2\Psi_0)\nabla\Psi_0] = (\nabla^2\Psi_0)^2 + \Psi_{0,k}\,(\nabla^2\Psi_0)_{,k} \tag{3.148a}$$

The leading-order Hamiltonian constraint (3.139) gives $(\nabla^2\Psi_0)^2 - \sum\Psi_{0,ij}^2 = 16\pi G\rho_m$. Combining this with (3.138) applied to $\Psi_0$:

$$\sum\Psi_{0,ij}^2 + \Psi_{0,k}\,(\nabla^2\Psi_0)_{,k} = -\nabla^2\Phi_N \tag{3.148b}$$

Adding the Hamiltonian constraint to (3.148b) eliminates $\sum\Psi_{0,ij}^2$:

$$(\nabla^2\Psi_0)^2 + \Psi_{0,k}\,(\nabla^2\Psi_0)_{,k} = 16\pi G\rho_m - \nabla^2\Phi_N \tag{3.148c}$$

The left side is (3.148a). Substituting into (3.144):

$$\nabla^2\Phi_N = 16\pi G\rho_m - (16\pi G\rho_m - \nabla^2\Phi_N) = \nabla^2\Phi_N \tag{3.148d}$$

The identity is satisfied automatically at leading order. This is a consistency verification, not an independent derivation of the Poisson equation from (3.144): the Poisson equation $\nabla^2\Phi = 4\pi G\rho_m$ is the definition of $\Phi_N$, and the identity (3.148d) confirms that the exact (3.144) is compatible with this definition. The exact equation is the all-orders generalisation from which $\Phi_N$ emerges as the leading-order piece:

$$\nabla^2\Phi = 4\pi G\,\rho_m + O(\varepsilon\cdot G\rho_m) \tag{3.148}$$

The $O(\varepsilon)$ corrections are the first post-Newtonian terms, encoded in the deviation $\delta\Psi$ and the post-Newtonian corrections to $\rho_E$. The sourced wave equation (Proposition 3.1, Eq 3.42a) follows from the full time-dependent ADM evolution equation for $K_{ij}$, which introduces the $c^{-2}\partial_T^2\Phi$ term.

3.11.7 Uniqueness: The Weinberg–Deser–Lovelock Theorems

Weinberg [160, 161] proved that the unique nonlinear field equation for a massless spin-2 field satisfying four premises is $G_{\mu\nu} + \Lambda g_{\mu\nu} = (8\pi G/c^4)T_{\mu\nu}$. We verify each premise for the ether.

W1 (Lorentz invariance). The Einstein tensor is constructed from $g_{\mu\nu}$ by coordinate-invariant operations. The Christoffel symbols:

$$\Gamma^\alpha_{\mu\nu} = \frac{1}{2}\,g^{\alpha\beta}\!\left(\partial_\mu g_{\beta\nu} + \partial_\nu g_{\beta\mu} - \partial_\beta g_{\mu\nu}\right) \tag{3.149}$$

transform as connection coefficients under coordinate changes. The Riemann tensor:

$$R^\alpha{}_{\beta\gamma\delta} = \partial_\gamma\Gamma^\alpha_{\beta\delta} - \partial_\delta\Gamma^\alpha_{\beta\gamma} + \Gamma^\alpha_{\gamma\sigma}\Gamma^\sigma_{\beta\delta} - \Gamma^\alpha_{\delta\sigma}\Gamma^\sigma_{\beta\gamma} \tag{3.150}$$

transforms as a rank-4 tensor. The Einstein tensor $G_{\mu\nu} = R_{\mu\nu} - \frac{1}{2}g_{\mu\nu}R$ is a rank-2 tensor. The equation $G_{\mu\nu} = (8\pi G/c^4)T_{\mu\nu}$ is therefore a tensor equation, covariant under all coordinate transformations and a fortiori under Lorentz transformations.

The ether rest frame is a property of the solution $V^i$, not of the field equation. The equation $G_{\mu\nu} = (8\pi G/c^4)T_{\mu\nu}$ makes no reference to $V^i$ — it is a condition on $g_{\mu\nu}$ and $T_{\mu\nu}$, both covariant objects. Emergent Lorentz invariance at wavelengths $\lambda \gg \ell_e$ is established by Theorem 3.3.

W2 (Newtonian limit). The Hamiltonian constraint reduces to $\nabla^2\Phi = 4\pi G\rho_m$ in the weak-field limit: Eq (3.148), derived in Section 3.11.6.

W3 (Energy-momentum conservation). The contracted Bianchi identity $\nabla_\mu G^{\mu\nu} \equiv 0$ holds for the Einstein tensor of any metric. It is a consequence of the algebraic symmetries of the Riemann tensor: antisymmetry in the first and second pairs of indices ($R_{\alpha\beta\gamma\delta} = -R_{\beta\alpha\gamma\delta}$, $R_{\alpha\beta\gamma\delta} = -R_{\alpha\beta\delta\gamma}$), pair symmetry ($R_{\alpha\beta\gamma\delta} = R_{\gamma\delta\alpha\beta}$), and the first Bianchi identity ($R_{\alpha[\beta\gamma\delta]} = 0$). From these, the second Bianchi identity $\nabla_{[\epsilon}R_{\alpha\beta]\gamma\delta} = 0$ follows, and contracting yields $\nabla_\mu G^{\mu\nu} = 0$. The field equation $G_{\mu\nu} = (8\pi G/c^4)T_{\mu\nu}$ then implies $\nabla_\mu T^{\mu\nu} = 0$.

The ether's continuity and Euler equations (Eqs 3.1–3.2) are the non-relativistic components of $\nabla_\mu T^{\mu\nu} = 0$. The acoustic metric of Theorem 3.1 is derived from these conservation laws. Energy-momentum conservation is the starting point of the ether framework, not an additional assumption.

W4 (Second-order field equations). The metric (3.92) depends on $V^i$ without derivatives: $g_{00} = -(c^2 - V^2)$, $g_{0i} = V_i$, $g_{ij} = \delta_{ij}$. The Christoffel symbols (3.149) involve $\partial_\mu g_{\alpha\beta}$, which for $g_{00}$ gives $\partial_j(V^2) = 2V_m V_{m,j}$ (first derivatives of $V^i$) and for $g_{0i}$ gives $\partial_j V_i = V_{i,j}$ (first derivatives). The Riemann tensor (3.150) involves $\partial_\gamma\Gamma^\alpha_{\beta\delta}$, hence second derivatives $V_{i,jk}$. The Einstein tensor (contractions of Riemann) involves at most second derivatives of $V^i$.

In the ADM formulation: $K_{ij} = (V_{i,j}+V_{j,i})/(2c)$ involves first derivatives; the Hamiltonian constraint $K^2 - K_{ij}K^{ij}$ involves products of first derivatives, which are algebraically equivalent to second-order conditions on $V^i$ through the divergence identity (3.143).

All four premises are satisfied.

Independent uniqueness results. Two additional theorems converge to the same conclusion:

Deser [162] proved that starting from the free Fierz-Pauli action for a massless spin-2 field and requiring self-consistent coupling to its own stress-energy tensor, the unique all-orders result is the Einstein-Hilbert action $S = \int R\sqrt{-g}\,d^4x$.

Lovelock [163, 164] proved that in four spacetime dimensions, the only symmetric, divergence-free rank-2 tensor constructed from the metric and its first and second derivatives is $a\,G_{\mu\nu} + b\,g_{\mu\nu}$ for constants $a$, $b$. With the normalisation fixed by the Newtonian limit ($a = 1$) and setting $b = -\Lambda$, the unique field equation is $G_{\mu\nu} + \Lambda g_{\mu\nu} = (8\pi G/c^4)T_{\mu\nu}$.

3.11.8 The Nonlinear Field Equation

Remark on $\Lambda$. The cosmological constant is not determined by the uniqueness argument. In Section 4.3, Theorem 4.2 identifies $\Lambda$ with the phonon zero-point field energy density, giving $w = -1$.

3.12 Post-Newtonian Parameters from the Ether Metric

The parameterized post-Newtonian (PPN) formalism [166] provides a model-independent framework for comparing metric theories of gravity in the weak-field, slow-motion regime. The formalism characterises the leading deviations from Newtonian gravity through ten dimensionless parameters, of which two — $\beta$ and $\gamma$ — dominate the observable predictions in the solar system. General relativity predicts $\beta = \gamma = 1$; alternative metric theories generically predict different values.

This section computes $\beta$ and $\gamma$ directly from the ether metric by transforming the Painlevé–Gullstrand line element (Theorem 3.2) into the standard PPN coordinate gauge and expanding to the required post-Newtonian order. Every step is an explicit coordinate transformation or algebraic expansion; no result from general relativity is imported. We then show that the remaining eight PPN parameters vanish identically as a consequence of Theorem 3.5.

3.12.1 The PPN Metric Template

The standard PPN metric for a static, spherically symmetric source of mass $M$, expressed in isotropic coordinates $(\bar{t}, \bar{\rho}, \theta, \phi)$ and expanded in powers of the dimensionless potential $\varepsilon = GM/(\bar{\rho}\,c^2)$, is [166]:

$$g_{\bar{t}\bar{t}} = -\!\left(1 - 2\varepsilon + 2\beta\,\varepsilon^2 + O(\varepsilon^3)\right)c^2 \tag{3.156}$$ $$g_{\bar{\rho}\bar{\rho}} = \frac{g_{\theta\theta}}{\bar{\rho}^2} = \frac{g_{\phi\phi}}{\bar{\rho}^2\sin^2\!\theta}$$ $$= 1 + 2\gamma\,\varepsilon + O(\varepsilon^2) \tag{3.157}$$ $$g_{\bar{t}\bar{\rho}} = 0 \tag{3.158}$$

The parameter $\gamma$ measures how much spatial curvature is produced per unit rest mass. The parameter $\beta$ measures the nonlinearity of superposition in the gravitational potential — the degree to which gravity gravitates. Both are dimensionless. In general relativity $\beta = \gamma = 1$; in Brans–Dicke theory $\gamma = (1+\omega)/(2+\omega)$ and $\beta = 1$; in many scalar-tensor theories both deviate from unity [166].

Current experimental constraints are:

$$|\gamma - 1| < 2.3 \times 10^{-5} \qquad \text{(Cassini spacecraft tracking [45])} \tag{3.159}$$ $$|4\beta - \gamma - 3| < 5.6 \times 10^{-3} \qquad \text{(Lunar laser ranging [168])} \tag{3.160}$$

The task is to compute $\beta$ and $\gamma$ from the ether metric (Eq 3.92).

3.12.2 Coordinate Transformation: PG to Schwarzschild

The ether metric in Painlevé–Gullstrand coordinates $(T, r, \theta, \phi)$ is (Eq 3.21):

$$ds^2 = -\!\left(c^2 - v^2\right)dT^2 + 2v\,dT\,dr + dr^2 + r^2\,d\Omega^2 \tag{3.161}$$

where $v(r) = \sqrt{2GM/r}$ is the ether inflow speed and $d\Omega^2 = d\theta^2 + \sin^2\!\theta\,d\phi^2$. This is the metric established by Theorem 3.2 as the acoustic metric for a constant-density ether flowing inward at the Newtonian free-fall velocity.

To transform to the PPN gauge, we first eliminate the cross-term $g_{Tr}$ by passing to Schwarzschild time $t_s$. Define:

$$dT = dt_s + h(r)\,dr \tag{3.162}$$

where $h(r)$ is chosen to make $g_{t_s r} = 0$. Substituting (3.162) into (3.161):

$$ds^2 = -(c^2 - v^2)(dt_s + h\,dr)^2 + 2v(dt_s + h\,dr)\,dr + dr^2 + r^2 d\Omega^2$$

Expanding the squares and collecting by coordinate pairs:

$$ds^2 = -(c^2 - v^2)\,dt_s^2 - \left[2(c^2-v^2)h - 2v\right]dt_s\,dr$$ $$+ \left[1 - (c^2-v^2)h^2 + 2vh\right]dr^2 + r^2\,d\Omega^2 \tag{3.163}$$

Setting the $dt_s\,dr$ coefficient to zero:

$$(c^2 - v^2)\,h - v = 0 \tag{3.164}$$ $$h(r) = \frac{v}{c^2 - v^2} = \frac{\sqrt{2GM/r}}{c^2(1 - r_s/r)} \tag{3.165}$$

where $r_s = 2GM/c^2$. The $dr^2$ coefficient with this choice of $h$ simplifies as follows. From (3.164): $h = v/(c^2-v^2)$, so $vh = v^2/(c^2-v^2)$ and $(c^2-v^2)h^2 = v^2/(c^2-v^2)$. Substituting into the $dr^2$ coefficient in (3.163):

$$1 - \frac{v^2}{c^2-v^2} + 2\!\left(\frac{v^2}{c^2-v^2}\right) = 1 + \frac{v^2}{c^2-v^2}$$ $$= \frac{c^2-v^2+v^2}{c^2-v^2} = \frac{c^2}{c^2-v^2} = \frac{1}{1 - r_s/r} \tag{3.166}$$

The metric in Schwarzschild coordinates is therefore:

$$ds^2 = -\!\left(1 - \frac{r_s}{r}\right)c^2\,dt_s^2 + \frac{dr^2}{1 - r_s/r} + r^2\,d\Omega^2 \tag{3.167}$$

This is the standard Schwarzschild metric, confirming that the PG → Schwarzschild transformation is exact and introduces no approximation.

3.12.3 Coordinate Transformation: Schwarzschild to Isotropic

The PPN parameters are defined in isotropic coordinates, in which the spatial metric is conformally flat. Define the isotropic radial coordinate $\bar{\rho}$ by:

$$r = \bar{\rho}\!\left(1 + \frac{r_s}{4\bar{\rho}}\right)^{\!2} \tag{3.168}$$

We verify that this transformation brings (3.167) into isotropic form by computing $dr$ and the metric functions in terms of $\bar{\rho}$.

Derivative. Define $u = r_s/(4\bar{\rho})$. Then $r = \bar{\rho}(1+u)^2$ and:

$$\frac{dr}{d\bar{\rho}} = (1+u)^2 + \bar{\rho}\cdot 2(1+u)\frac{du}{d\bar{\rho}} \tag{3.169a}$$

Since $u = r_s/(4\bar{\rho})$: $du/d\bar{\rho} = -r_s/(4\bar{\rho}^2) = -u/\bar{\rho}$. Substituting:

$$\frac{dr}{d\bar{\rho}} = (1+u)^2 - 2u(1+u) = (1+u)\!\left[(1+u) - 2u\right]$$ $$= (1+u)(1-u) = 1 - u^2 \tag{3.169b}$$

The $g_{tt}$ factor. We compute $1 - r_s/r$:

$$\frac{r_s}{r} = \frac{r_s}{\bar{\rho}(1+u)^2} = \frac{4\bar{\rho}\,u}{\bar{\rho}(1+u)^2} = \frac{4u}{(1+u)^2} \tag{3.170a}$$ $$1 - \frac{r_s}{r} = \frac{(1+u)^2 - 4u}{(1+u)^2} = \frac{1 + 2u + u^2 - 4u}{(1+u)^2} = \frac{(1-u)^2}{(1+u)^2} \tag{3.170b}$$

The $g_{rr}$ factor. Using $dr = (1-u^2)\,d\bar{\rho} = (1-u)(1+u)\,d\bar{\rho}$:

$$\frac{dr^2}{1-r_s/r} = \frac{(1-u)^2(1+u)^2}{(1-u)^2/(1+u)^2}\,d\bar{\rho}^2 = (1+u)^4\,d\bar{\rho}^2 \tag{3.171}$$

The angular part. Since $r^2 = \bar{\rho}^2(1+u)^4$:

$$r^2\,d\Omega^2 = \bar{\rho}^2(1+u)^4\,d\Omega^2 \tag{3.172}$$

The spatial part of the metric combines (3.171) and (3.172):

$$d\ell^2 = (1+u)^4\!\left(d\bar{\rho}^2 + \bar{\rho}^2\,d\Omega^2\right) \tag{3.173}$$

This is conformally flat — the spatial metric is $(1+u)^4$ times the flat metric in spherical coordinates — confirming that $\bar{\rho}$ is the isotropic radial coordinate.

The complete isotropic metric. Assembling (3.170b) and (3.173):

$$ds^2 = -\frac{(1-u)^2}{(1+u)^2}\,c^2\,dt_s^2 + (1+u)^4\!\left(d\bar{\rho}^2 + \bar{\rho}^2\,d\Omega^2\right) \tag{3.174}$$

with $u = r_s/(4\bar{\rho}) = GM/(2c^2\bar{\rho})$. This is exact — no approximation has been made.

3.12.4 Post-Newtonian Expansion

We now expand (3.174) in powers of $u = GM/(2c^2\bar{\rho})$, which is the post-Newtonian expansion parameter. Define the Newtonian potential at isotropic radius $\bar{\rho}$:

$$U(\bar{\rho}) = \frac{GM}{\bar{\rho}} \tag{3.175}$$

so that $u = U/(2c^2)$ and $\varepsilon = U/c^2 = 2u$.

Expansion of $g_{\bar{t}\bar{t}}$. The temporal component is:

$$\frac{g_{\bar{t}\bar{t}}}{c^2} = -\frac{(1-u)^2}{(1+u)^2} \tag{3.176a}$$

Expand $(1-u)/(1+u)$ as a geometric series. Since $(1+u)^{-1} = 1 - u + u^2 - u^3 + O(u^4)$:

$$(1-u)(1-u+u^2-u^3+\cdots) = 1 - 2u + 2u^2 - 2u^3 + O(u^4) \tag{3.176b}$$

Squaring (3.176b):

$$\frac{(1-u)^2}{(1+u)^2} = \left(1 - 2u + 2u^2 - 2u^3 + \cdots\right)^2 \tag{3.176c}$$

The first three orders are:

$$= 1 + 2(-2u) + \left[(-2u)^2 + 2(2u^2)\right] + O(u^3)$$ $$= 1 - 4u + \left(4u^2 + 4u^2\right) + O(u^3) = 1 - 4u + 8u^2 + O(u^3) \tag{3.176d}$$

Verification of the $u^2$ coefficient: the square of $1 - 2u + 2u^2$ gives $1 - 4u + (4+4)u^2 = 1 - 4u + 8u^2$ at order $u^2$, where the contributions are $(-2u)^2 = 4u^2$ and $2 \times 1 \times 2u^2 = 4u^2$.

Substituting $u = U/(2c^2)$:

$$-\frac{g_{\bar{t}\bar{t}}}{c^2} = 1 - 4\cdot\frac{U}{2c^2} + 8\cdot\frac{U^2}{4c^4} + O(c^{-6})$$ $$= 1 - \frac{2U}{c^2} + \frac{2U^2}{c^4} + O(c^{-6}) \tag{3.177}$$

Expansion of $g_{\bar{\rho}\bar{\rho}}$. The spatial conformal factor is:

$$(1+u)^4 = 1 + 4u + 6u^2 + 4u^3 + u^4 \tag{3.178a}$$

This is exact (binomial theorem with integer exponent). Substituting $u = U/(2c^2)$:

$$(1+u)^4 = 1 + 4\cdot\frac{U}{2c^2} + 6\cdot\frac{U^2}{4c^4} + O(c^{-6})$$ $$= 1 + \frac{2U}{c^2} + \frac{3U^2}{2c^4} + O(c^{-6}) \tag{3.178b}$$

3.12.5 Extraction of $\beta$ and $\gamma$

The parameter $\gamma$. Comparing the spatial metric expansion (3.178b) with the PPN template (3.157):

$$(1+u)^4 = 1 + \frac{2U}{c^2} + O(c^{-4}) \quad \longleftrightarrow \quad 1 + \frac{2\gamma\,U}{c^2} + O(c^{-4})$$

The coefficients of $U/c^2$ must agree:

$$\boxed{\gamma = 1} \tag{3.179}$$

The parameter $\beta$. Comparing the temporal metric expansion (3.177) with the PPN template (3.156):

$$1 - \frac{2U}{c^2} + \frac{2U^2}{c^4} \quad \longleftrightarrow \quad 1 - \frac{2U}{c^2} + \frac{2\beta\,U^2}{c^4}$$

The coefficients of $U^2/c^4$ must agree:

$$\boxed{\beta = 1} \tag{3.180}$$

3.12.6 Observational Consequences

The PPN parameters $\beta$ and $\gamma$ enter all classical tests of gravity. We collect the predictions, all of which agree with the results derived independently in Section 3.5 from the PG metric.

Shapiro time delay (Section 3.5.3). The excess round-trip time for a signal passing a mass $M$ at closest approach $R$ depends on $\gamma$ [166]:

$$\Delta t = \frac{(1+\gamma)\,GM}{c^3}\ln\!\left(\frac{4r_1 r_2}{R^2}\right) \tag{3.182}$$

For $\gamma = 1$, this reproduces Eq (3.32). The Cassini measurement [45] confirms $\gamma = 1$ to $2.3 \times 10^{-5}$.

Light bending (Section 3.5.2). The deflection angle for a photon passing mass $M$ at impact parameter $R$ is:

$$\Delta\theta = \frac{(1+\gamma)}{2}\cdot\frac{4GM}{Rc^2} \tag{3.183}$$

For $\gamma = 1$: $\Delta\theta = 4GM/(Rc^2)$, reproducing Eq (3.30). VLBI measurements confirm this to $10^{-4}$ accuracy [43].

Perihelion precession (Section 3.5.4). The precession rate depends on both $\beta$ and $\gamma$ [166]:

$$\Delta\phi = \frac{(2 - \beta + 2\gamma)}{3}\cdot\frac{6\pi GM}{ac^2(1-e^2)} \tag{3.184}$$

For $\beta = \gamma = 1$: the prefactor is $(2-1+2)/3 = 1$, recovering the standard result $\Delta\phi = 6\pi GM/[ac^2(1-e^2)]$ (Eq 3.33). Mercury's measured precession of $(42.98 \pm 0.04)''\text{/century}$ [46] is consistent.

Nordtvedt effect. The strong equivalence principle — the statement that gravitational binding energy falls at the same rate as all other forms of energy — requires the combination $4\beta - \gamma - 3 = 0$ [169]. For $\beta = \gamma = 1$: $4(1) - 1 - 3 = 0$. The ether framework satisfies the strong equivalence principle. Lunar laser ranging constrains $|4\beta - \gamma - 3| < 5.6 \times 10^{-3}$ [168], consistent with the ether prediction.

3.12.7 The Preferred-Frame Parameters

The ether possesses a preferred rest frame, identified in Section 4.1.3 with the CMB rest frame ($v_{\text{CMB}} = 369.82 \pm 0.11$ km/s). A natural concern is that this preferred frame might produce observable violations of local Lorentz invariance, parameterised in the PPN formalism by the preferred-frame parameters $\alpha_1$, $\alpha_2$, and $\alpha_3$ [166]. These parameters are tightly constrained: $|\alpha_1| < 10^{-4}$, $|\alpha_2| < 10^{-7}$, $|\alpha_3| < 4 \times 10^{-20}$ [166].

The ether framework predicts all three to be exactly zero. The argument has two levels.

(i) From the field equation. Theorem 3.5 establishes that the ether's field equation is the Einstein equation $G_{\mu\nu} = (8\pi G/c^4)T_{\mu\nu}$. This equation is manifestly covariant — it contains no preferred vector, no coupling to the CMB frame, and no dependence on the observer's velocity through the ether. The preferred-frame PPN parameters $\alpha_1$, $\alpha_2$, $\alpha_3$ vanish for any metric theory whose field equation is the Einstein equation [166].

(ii) From the ether's structure. The preferred frame is a property of the ether's state — the background flow $V^i$ — not of the law governing the flow. The field equation determines which flow patterns are physical; the law is covariant even though the solutions are not. This is precisely analogous to electrodynamics: Maxwell's equations are Lorentz covariant, but a specific electromagnetic field configuration (e.g., a uniform magnetic field) defines a preferred direction. No one concludes that electrodynamics violates rotational invariance.

The ether's preferred frame is observationally inaccessible to gravitational experiments because the field equation that governs gravity is frame-independent. The preferred frame may, in principle, be detectable through non-gravitational channels — for example, through anisotropy of the zero-point field spectrum at order $v_{\text{CMB}}/c \sim 10^{-3}$ (see Section 9.3.1) — but all PPN parameters, including the preferred-frame parameters, are those of general relativity:

$$\alpha_1 = \alpha_2 = \alpha_3 = 0 \tag{3.185}$$

3.12.8 The Complete PPN Parameter Set

The full PPN formalism involves ten parameters. For any metric theory whose field equation is the Einstein equation (with or without a cosmological constant), all ten are fixed [166]:

$$\gamma = \beta = 1, \qquad \xi = \alpha_1 = \alpha_2 = \alpha_3 = \zeta_1 = \zeta_2 = \zeta_3 = \zeta_4 = 0 \tag{3.186}$$

By Theorem 3.5, the ether's field equation is the Einstein equation. Therefore (3.186) holds for the ether framework. The ten PPN parameters are not ten independent predictions; they are a single consequence of the uniqueness result.

The physical content of (3.186) is threefold:

(a) The ether framework conserves total momentum ($\zeta_1 = \zeta_2 = \zeta_3 = \zeta_4 = 0$). Gravitational radiation carries momentum; the ether medium absorbs and transmits it via the mechanism of Section 3.7.

(b) The ether framework satisfies the strong equivalence principle ($4\beta - \gamma - 3 = 0$ and $\alpha_3 = \zeta_1 = 0$). Gravitational binding energy gravitates identically to all other energy forms.

(c) The ether framework exhibits no preferred-frame effects in gravitational experiments ($\alpha_1 = \alpha_2 = \alpha_3 = 0$), despite the physical existence of the ether rest frame.

These three properties are inherited from the Einstein equation via Theorem 3.5. They are not imposed as additional assumptions but derived from the ether's dynamics.

3.13 Hawking Radiation from the Ether Horizon

Section 3.6 established that the ether horizon — the surface at which the ether inflow velocity equals the speed of light — is a smooth, regular feature of the PG metric. Section 6 established that the ether carries a zero-point electromagnetic field whose spectral density is uniquely determined by Lorentz invariance (Theorem 4.2). The combination of these two results — a horizon in a medium that carries quantum fluctuations — produces thermal radiation. This is the ether's prediction of the Hawking effect [170], derived here directly from the ether's flow profile and mode structure without importing quantum field theory on curved spacetime.

The derivation follows the logic of Unruh's 1981 result [10]: a flowing medium with a sonic horizon emits thermal radiation at a temperature set by the velocity gradient at the horizon. In the ether framework, the "sonic horizon" is the gravitational horizon (where the ether inflow speed equals $c$), the "sonic" modes are the ether's electromagnetic perturbations, and the velocity gradient is computed from the Newtonian free-fall profile. The result is the Hawking temperature.

3.13.1 Surface Gravity from the Ether Flow Profile

The ether inflow velocity is $v(r) = \sqrt{2GM/r} = c\sqrt{r_s/r}$ (Eq 3.22), with $r_s = 2GM/c^2$. At the horizon $r = r_s$: $v(r_s) = c$.

The velocity gradient at the horizon. Differentiating $v(r) = c(r_s/r)^{1/2}$:

$$\frac{dv}{dr} = -\frac{c}{2}\left(\frac{r_s}{r}\right)^{1/2}\frac{1}{r} = -\frac{v}{2r} \tag{3.187}$$

At $r = r_s$:

$$\left.\frac{dv}{dr}\right|_{r_s} = -\frac{c}{2r_s} \tag{3.188}$$

Outgoing null rays. In PG coordinates, outgoing radial null rays satisfy $ds^2 = 0$ with $d\Omega = 0$. From the PG metric (3.21):

$$-(c^2 - v^2)\,dT^2 - 2v\,dT\,dr + dr^2 = 0 \tag{3.189}$$

Dividing by $dT^2$ and solving the resulting quadratic in $dr/dT$:

$$\left(\frac{dr}{dT}\right)^2 - 2v\frac{dr}{dT} - (c^2 - v^2) = 0 \tag{3.190}$$ $$\frac{dr}{dT} = v \pm \sqrt{v^2 + c^2 - v^2} = v \pm c \tag{3.191}$$

The outgoing ray has $dr/dT = c - v$ (choosing the sign that gives $dr/dT > 0$ outside the horizon where $v < c$); the ingoing ray has $dr/dT = -(c + v) < 0$.

Near-horizon expansion. Define $\delta r = r - r_s$. For $|\delta r| \ll r_s$, expand $v(r)$ to first order:

$$v(r) = v(r_s) + \left.\frac{dv}{dr}\right|_{r_s}\delta r + O(\delta r^2)$$ $$= c - \frac{c}{2r_s}\,\delta r + O(\delta r^2) \tag{3.192}$$

The outgoing null ray velocity becomes:

$$\frac{d(\delta r)}{dT} = c - v(r) = \frac{c}{2r_s}\,\delta r + O(\delta r^2) \tag{3.193}$$

This is a linear ODE with solution:

$$\delta r(T) = \delta r_0 \cdot \exp\!\left(\frac{c\,T}{2r_s}\right) \tag{3.194}$$

Outgoing null rays peel away from the horizon exponentially, with e-folding rate $c/(2r_s)$.

Surface gravity. Define the surface gravity $\kappa$ as the exponential peeling rate times $c$:

$$\kappa = c \cdot \frac{c}{2r_s} = \frac{c^2}{2r_s} = \frac{c^4}{4GM} \tag{3.195}$$

This definition is equivalent to the standard general-relativistic definition $\kappa^2 = -\frac{1}{2}(\nabla^\mu\xi^\nu)(\nabla_\mu\xi_\nu)\big|_{\text{horizon}}$ for the Killing vector $\xi = \partial_T$. We verify the equivalence. The Killing vector has norm $\xi_\mu\xi^\mu = g_{TT} = -(c^2 - v^2)$. Its gradient with respect to $r$ is:

$$\partial_r(\xi_\mu\xi^\mu) = \partial_r(-(c^2 - v^2)) = 2v\,\frac{dv}{dr} \tag{3.196a}$$

At the horizon: $2c \cdot (-c/(2r_s)) = -c^2/r_s$. The standard identity [165] relating this gradient to the surface gravity for a Killing horizon is $\partial_r(\xi_\mu\xi^\mu) = -2\kappa\sqrt{-g_{rr}^{-1}\,\xi_\mu\xi^\mu + g_{Tr}^2/g_{rr}}$ evaluated at the horizon. Since $\xi_\mu\xi^\mu = 0$ at $r = r_s$, the expression reduces to $-2\kappa\,|g_{Tr}|/\sqrt{g_{rr}} = -2\kappa c$. Setting this equal to $-c^2/r_s$:

$$2\kappa c = \frac{c^2}{r_s} \quad \Rightarrow \quad \kappa = \frac{c}{2r_s} \cdot c = \frac{c^2}{2r_s} \tag{3.196b}$$

confirming (3.195).

Physical interpretation. The surface gravity measures the acceleration required for a static observer to hover at the horizon, as measured at spatial infinity. In the ether picture, it measures how rapidly the ether flow accelerates through the sonic point. The velocity gradient $|dv/dr|_{r_s} = c/(2r_s)$ sets the scale; multiplication by $c$ converts from the "spatial peeling rate" to the "proper acceleration" (surface gravity).

3.13.2 The Wave Equation on the Ether Metric

A massless scalar field $\phi$ on the ether metric satisfies the covariant wave equation:

$$\Box\phi = \frac{1}{\sqrt{-g}}\,\partial_\mu\!\left(\sqrt{-g}\;g^{\mu\nu}\,\partial_\nu\phi\right) = 0 \tag{3.197}$$

The metric determinant. From the PG metric (3.21) in coordinates $(T, r, \theta, \phi)$, the $(T, r)$ block has components $g_{TT} = -(c^2 - v^2)$, $g_{Tr} = -v$, $g_{rr} = 1$. Its determinant is $-(c^2 - v^2) \cdot 1 - (-v)^2 = -c^2$. Including the angular part ($g_{\theta\theta} = r^2$, $g_{\phi\phi} = r^2\sin^2\!\theta$):

$$\det(g_{\mu\nu}) = -c^2\,r^4\sin^2\!\theta \tag{3.198}$$ $$\sqrt{-g} = c\,r^2\sin\theta \tag{3.199}$$

The inverse metric. From the $(T, r)$ block inversion (with determinant $-c^2$):

$$g^{TT} = -\frac{1}{c^2}, \qquad g^{Tr} = -\frac{v}{c^2}, \qquad g^{rr} = \frac{c^2 - v^2}{c^2} \tag{3.200}$$

and $g^{\theta\theta} = 1/r^2$, $g^{\phi\phi} = 1/(r^2\sin^2\!\theta)$.

Mode decomposition. Decompose $\phi = r^{-1}\psi(T, r)\,Y_{lm}(\theta, \phi)$. The angular part of (3.197) produces $-l(l+1)/r^2$ by the eigenvalue equation for spherical harmonics. For the essential near-horizon physics, it suffices to work in the $s$-wave sector ($l = 0$), where the angular momentum barrier is absent and the derivation is cleanest. The geometric optics approximation used below is valid for all $l$ provided $\omega \gg l\,c/r_s$ (the mode frequency is much higher than the angular momentum barrier frequency).

In the $s$-wave sector ($l = 0$), we substitute $\phi = \psi(T, r)/r$ into (3.197) and expand using (3.199)–(3.200). The wave (3.197) becomes:

$$\frac{1}{cr^2}\,\partial_\mu\!\left(cr^2\,g^{\mu\nu}\,\partial_\nu\!\left(\frac{\psi}{r}\right)\right) = 0$$

Derivatives of $\psi/r$. Since $\partial_T(1/r) = 0$ and $\partial_r(1/r) = -1/r^2$:

$$\partial_T(\psi/r) = \dot{\psi}/r, \qquad \partial_r(\psi/r) = \psi'/r - \psi/r^2 \tag{3.200a}$$

where $\dot{\psi} = \partial_T\psi$ and $\psi' = \partial_r\psi$.

The flux components. Define $J^\mu = cr^2 g^{\mu\nu}\partial_\nu(\psi/r)$. Using the inverse metric (3.200):

$$J^T = cr^2\!\left[-\frac{\dot{\psi}}{c^2 r} - \frac{v}{c^2}\!\left(\frac{\psi'}{r}-\frac{\psi}{r^2}\right)\right]$$ $$= -\frac{r}{c}\!\left[\dot{\psi} + v\psi' - \frac{v\psi}{r}\right] \tag{3.200b}$$ $$J^r = cr^2\!\left[-\frac{v\dot{\psi}}{c^2 r} + \frac{c^2-v^2}{c^2}\!\left(\frac{\psi'}{r}-\frac{\psi}{r^2}\right)\right]$$ $$= \frac{r}{c}\!\left[-v\dot{\psi} + (c^2-v^2)\!\left(\psi'-\frac{\psi}{r}\right)\right] \tag{3.200c}$$

Time derivative of $J^T$. Since $v = v(r)$ is static:

$$\partial_T J^T = -\frac{r}{c}\!\left[\ddot{\psi} + v\dot{\psi}' - \frac{v\dot{\psi}}{r}\right] \tag{3.200d}$$

Radial derivative of $J^r$. Write $J^r = J^r_1 + J^r_2 + J^r_3$ with $J^r_1 = -rv\dot{\psi}/c$, $J^r_2 = r(c^2-v^2)\psi'/c$, $J^r_3 = -(c^2-v^2)\psi/c$. Computing each derivative by the product rule:

$$\partial_r J^r_1 = -\frac{1}{c}(v\dot{\psi} + rv'\dot{\psi} + rv\dot{\psi}') \tag{3.200e}$$ $$\partial_r J^r_2 = \frac{1}{c}[(c^2-v^2)\psi' - 2rvv'\psi' + r(c^2-v^2)\psi''] \tag{3.200f}$$ $$\partial_r J^r_3 = \frac{1}{c}[2vv'\psi - (c^2-v^2)\psi'] \tag{3.200g}$$

Assembling the wave equation. Setting $\partial_T J^T + \partial_r J^r = 0$, multiplying through by $-c/r$, and collecting all terms:

From (3.200d): $\ddot{\psi} + v\dot{\psi}' - v\dot{\psi}/r$.

From (3.200e): $+v\dot{\psi}/r + v'\dot{\psi} + v\dot{\psi}'$.

From (3.200f): $-(c^2-v^2)\psi'/r + 2vv'\psi' - (c^2-v^2)\psi''$.

From (3.200g): $-2vv'\psi/r + (c^2-v^2)\psi'/r$.

Cancellations. $-v\dot{\psi}/r + v\dot{\psi}/r = 0$. $-(c^2-v^2)\psi'/r + (c^2-v^2)\psi'/r = 0$. The surviving terms:

$$\ddot{\psi} + 2v\dot{\psi}' + v'\dot{\psi} - (c^2-v^2)\psi'' + 2vv'\psi' - \frac{2vv'\psi}{r} = 0 \tag{3.200h}$$

The last three terms have a hierarchy near the horizon. Using $v = c\sqrt{r_s/r}$, which gives $v' = -v/(2r)$:

The term $v'\dot{\psi}$: magnitude $\sim (c/r_s)\omega|\psi|$ for a mode of frequency $\omega$.

The term $2vv'\psi'$: magnitude $\sim (c^2/r_s)k|\psi|$ where $k$ is the local wavevector.

The term $2vv'\psi/r$: magnitude $\sim (c^2/r_s^2)|\psi|$.

The leading terms ($\ddot{\psi}, 2v\dot{\psi}', (c^2-v^2)\psi''$): magnitude $\sim \omega^2|\psi|$ or $c^2 k^2|\psi|$.

For modes with $\omega r_s/c \gg 1$ (geometric optics regime), the three subleading terms are suppressed by $c/(\omega r_s) \ll 1$. Dropping them:

$$\ddot{\psi} + 2v\dot{\psi}' - (c^2-v^2)\psi'' = 0$$

Dividing by $-c^2$:

$$-\frac{1}{c^2}\,\partial_T^2\psi - \frac{2v}{c^2}\,\partial_T\partial_r\psi$$ $$+ \frac{c^2 - v^2}{c^2}\,\partial_r^2\psi = 0 \tag{3.201}$$

This is the effective $1+1$ dimensional wave equation. The dropped terms modify the greybody factors (frequency-dependent transmission probability through the angular momentum barrier) but not the thermal spectrum, which depends only on the exponential peeling rate at the horizon (Section 3.13.3).

Multiplying through by $-c^2$:

$$\partial_T^2\psi + 2v\,\partial_T\partial_r\psi - (c^2 - v^2)\,\partial_r^2\psi = 0 \tag{3.202}$$

Factorisation. Define the comoving (material) derivative along the ether flow:

$$D_T = \partial_T - v\,\partial_r \tag{3.203}$$

This is the time derivative in the frame of a freely falling ether element (which moves at velocity $-v$ in the inward radial direction, so its position satisfies $dr/dT = -v$, giving $D_T f = (\partial_T + (-v)(-\partial_r))f$... more carefully: for a function evaluated along the flow, $f(T, r(T))$ with $dr/dT = -v$, the total derivative is $\partial_T f + (dr/dT)\partial_r f = \partial_T f - v\partial_r f = D_T f$). (3.202) can be written as:

$$(D_T - c\,\partial_r)(D_T + c\,\partial_r)\psi + [\text{terms involving } \partial_r v] = 0 \tag{3.204}$$

To verify: $(D_T - c\partial_r)(D_T + c\partial_r)\psi = D_T^2\psi - c^2\partial_r^2\psi + c(\partial_rv)(\partial_r\psi)$ where the cross-derivative terms involving $v\partial_r$ cancel pairwise and the commutator $[D_T, c\partial_r] = c(\partial_r v)\partial_r$ arises from $v$ depending on $r$. Expanding $D_T^2 = (\partial_T - v\partial_r)^2 = \partial_T^2 - 2v\partial_T\partial_r + v^2\partial_r^2 - v'(\partial_T r)\partial_r$; for the background flow with $\partial_T v = 0$ (static), this gives $D_T^2 = \partial_T^2 - 2v\partial_T\partial_r + v^2\partial_r^2$. So:

$$D_T^2\psi - c^2\partial_r^2\psi$$ $$= \partial_T^2\psi - 2v\partial_T\partial_r\psi + v^2\partial_r^2\psi - c^2\partial_r^2\psi$$ $$= \partial_T^2\psi - 2v\partial_T\partial_r\psi - (c^2 - v^2)\partial_r^2\psi$$

which is exactly the left side of (3.202). The commutator term $c(\partial_r v)(\partial_r\psi)$ in (3.204) is lower order in the geometric optics limit ($|\partial_r\psi| \ll |k\psi|$ where $k$ is the local wavevector), so the near-horizon wave equation factorises as:

$$(\partial_T - v\,\partial_r - c\,\partial_r)(\partial_T - v\,\partial_r + c\,\partial_r)\psi \approx 0 \tag{3.205}$$

The first factor gives outgoing modes (propagating at $c - v$ relative to the coordinates); the second gives ingoing modes (propagating at $-(c + v)$, always inward). We focus on the outgoing sector.

3.13.3 Near-Horizon Mode Structure

The retarded time. For an outgoing mode in the geometric optics limit, $\psi$ depends on $T$ and $r$ through the combination $u = T - \int dr/(c - v)$ — the retarded time, constant along outgoing null rays. This follows directly from (3.191): the outgoing null ray satisfies $dr/dT = c - v$, so $dT = dr/(c - v)$ along the ray, and $u = T - \int dr/(c - v)$ is constant.

Near the horizon, using (3.192):

$$c - v \approx \frac{c}{2r_s}\,(r - r_s) = \frac{\kappa}{c}\,(r - r_s) \tag{3.206}$$

The integral:

$$\int\frac{dr}{c - v} \approx \frac{c}{\kappa}\int\frac{dr}{r - r_s}$$ $$= \frac{c}{\kappa}\,\ln|r - r_s| + \text{const} \tag{3.207}$$

The retarded time is therefore:

$$u \approx T - \frac{c}{\kappa}\,\ln(r - r_s) + \text{const} \qquad (r > r_s) \tag{3.208}$$

As $r \to r_s^+$: $u \to +\infty$. The outgoing ray emitted from a point just outside the horizon at coordinate time $T$ arrives at infinity only after an arbitrarily long retarded time — this is the infinite redshift at the horizon.

Outgoing mode functions. A monochromatic outgoing mode with frequency $\omega$ (as measured by a static observer at infinity, using PG time $T$) has the form:

$$\psi_\omega^{\text{out}}(T, r) = e^{-i\omega u} = e^{-i\omega T}\,e^{i(\omega c/\kappa)\ln(r - r_s)} \tag{3.209}$$

for $r > r_s$, up to a normalisation constant and smooth corrections from the exact (non-near-horizon) geometry. Equivalently:

$$\psi_\omega^{\text{out}} = e^{-i\omega T}\,(r - r_s)^{i\omega c/\kappa} \qquad (r > r_s) \tag{3.210}$$

The key feature is the logarithmic phase singularity: as $r \to r_s^+$, the phase $\omega c\ln(r - r_s)/\kappa$ oscillates with unbounded frequency. This infinite blueshift near the horizon is the origin of particle creation.

3.13.4 Analytic Continuation and the Bogoliubov Coefficients

The problem. The mode (3.210) is defined only for $r > r_s$. To determine the particle content measured by a distant observer, we must relate the "in" modes (defined as positive-frequency with respect to the natural time coordinate of the ether at early times, before the horizon formed) to the "out" modes (defined as positive-frequency at infinity at late times). The "in" modes are smooth across the horizon; the "out" modes are singular there. The relation between them is a Bogoliubov transformation whose coefficients determine the spectrum of emitted radiation.

Analytic continuation. The "in" vacuum state corresponds to a mode that is smooth on the future horizon when expressed in terms of the affine parameter $\lambda$ along the horizon generator. The affine parameter is related to $r$ by $\lambda \propto -(r_s - r)$ for $r < r_s$ (ingoing from the exterior), so near the horizon $\lambda \propto -(r - r_s)$ for an ingoing ray that crosses at $\lambda = 0$.

The smooth continuation of $\psi_\omega^{\text{out}}$ from $r > r_s$ to $r < r_s$ is obtained by treating $(r - r_s)$ as a complex variable and choosing the analytic continuation that is regular on the upper-half complex plane (corresponding to the positive-frequency condition for the "in" mode). This is the standard $i\epsilon$ prescription [170, 171]:

$$(r - r_s) \to (r - r_s + i\epsilon), \qquad \epsilon \to 0^+ \tag{3.211}$$

For $r > r_s$: $(r - r_s + i\epsilon) \to (r - r_s)$ as $\epsilon \to 0$, recovering (3.210).

For $r < r_s$: $r - r_s < 0$, so $r - r_s + i\epsilon = |r - r_s|\,e^{i\pi}$ (since $r - r_s = -|r - r_s|$ and the $i\epsilon$ prescription selects the branch with argument $\pi$). Therefore:

$$(r - r_s + i\epsilon)^{i\omega c/\kappa}$$ $$\to |r - r_s|^{i\omega c/\kappa}\,\exp\!\left(i\cdot\frac{\omega c}{\kappa}\cdot i\pi\right)$$ $$= |r - r_s|^{i\omega c/\kappa}\,e^{-\pi\omega c/\kappa} \tag{3.212}$$

The continued mode acquires a real exponential factor $e^{-\pi\omega c/\kappa}$ upon crossing the horizon.

Decomposition into positive and negative frequencies. The analytically continued mode inside the horizon (3.212) can be decomposed as a linear combination of the positive-frequency outgoing mode and its negative-frequency conjugate. Outside the horizon ($r > r_s$), the mode is $\psi_\omega^{\text{out}} = (r-r_s)^{i\omega c/\kappa}$. Inside the horizon ($r < r_s$), we write:

$$|r - r_s|^{i\omega c/\kappa} = \left[-(r - r_s)\right]^{i\omega c/\kappa}$$ $$= (-1)^{i\omega c/\kappa}\,(r_s - r)^{i\omega c/\kappa}$$

The factor $(-1)^{i\omega c/\kappa} = e^{i\pi\cdot i\omega c/\kappa} = e^{-\pi\omega c/\kappa}$ confirms (3.212). Now, $(r_s - r)^{i\omega c/\kappa}$ for $r < r_s$ is the same power-law as $(r - r_s)^{i\omega c/\kappa}$ but evaluated at the reflected argument. The complex conjugate of the outgoing mode is $(r - r_s)^{-i\omega c/\kappa}$, which represents a negative-frequency (ingoing, partner) mode.

The "in" mode, which is analytic on the upper-half plane, contains both the outgoing mode and its time-reversed partner. The Bogoliubov decomposition is:

$$\psi_\omega^{\text{in}} = \alpha_\omega\,\psi_\omega^{\text{out}} + \beta_\omega\,(\psi_\omega^{\text{out}})^* \tag{3.213}$$

where $\alpha_\omega$ and $\beta_\omega$ are the Bogoliubov coefficients satisfying the normalisation condition $|\alpha_\omega|^2 - |\beta_\omega|^2 = 1$ (for bosonic fields).

The coefficient ratio. The ratio $|\beta_\omega/\alpha_\omega|$ is determined by the exponential factor acquired in the analytic continuation. Outside the horizon, the "in" mode matches the outgoing mode with unit amplitude (by definition of the "in" vacuum). Inside the horizon, the continuation (3.212) multiplies the mode by $e^{-\pi\omega c/\kappa}$. The complex conjugate mode, continued by the opposite prescription $(r - r_s - i\epsilon)$ (lower-half plane), acquires the factor $e^{+\pi\omega c/\kappa}$.

The relative amplitude of the negative-frequency component is therefore:

$$\frac{|\beta_\omega|}{|\alpha_\omega|} = e^{-\pi\omega c/\kappa} \tag{3.214}$$ $$\frac{|\beta_\omega|^2}{|\alpha_\omega|^2} = e^{-2\pi\omega c/\kappa} \tag{3.215}$$

3.13.5 The Thermal Spectrum

Particle number. The number of particles in mode $\omega$ measured by a distant observer, given that the ether ZPF (§6) defines the "in" vacuum, is:

$$\langle n_\omega \rangle = |\beta_\omega|^2 \tag{3.216}$$

From the normalisation $|\alpha_\omega|^2 - |\beta_\omega|^2 = 1$ and (3.215):

$$|\alpha_\omega|^2 = \frac{|\beta_\omega|^2}{e^{-2\pi\omega c/\kappa}}$$ $$= |\beta_\omega|^2\,e^{2\pi\omega c/\kappa} \tag{3.217}$$

Substituting into the normalisation:

$$|\beta_\omega|^2\,e^{2\pi\omega c/\kappa} - |\beta_\omega|^2 = 1$$ $$|\beta_\omega|^2\!\left(e^{2\pi\omega c/\kappa} - 1\right) = 1$$ $$\langle n_\omega \rangle = |\beta_\omega|^2 = \frac{1}{e^{2\pi\omega c/\kappa} - 1} \tag{3.218}$$

This is the Bose–Einstein distribution with temperature:

$$k_B T = \frac{\hbar\kappa}{2\pi c} \tag{3.219}$$

To verify: $\langle n_\omega \rangle = (e^{\hbar\omega/(k_BT)} - 1)^{-1}$ requires the exponent $2\pi\omega c/\kappa = \hbar\omega/(k_BT)$, giving $k_BT = \hbar\kappa/(2\pi c)$.

3.13.6 The Hawking Temperature

Numerical evaluation. For a solar-mass black hole ($M = M_\odot = 1.989 \times 10^{30}$ kg):

$$T_H = \frac{1.055 \times 10^{-34} \times (3 \times 10^8)^3}{8\pi \times 6.674 \times 10^{-11} \times 1.989 \times 10^{30} \times 1.381 \times 10^{-23}}$$ $$= \frac{2.849 \times 10^{-9}}{4.607 \times 10^{-2}}$$ $$\approx 6.2 \times 10^{-8}\;\text{K} \tag{3.221}$$

For a black hole of mass $M$:

$$T_H = 6.17 \times 10^{-8}\;\text{K} \times \frac{M_\odot}{M} \tag{3.222}$$

The Hawking temperature is extremely small for astrophysical black holes, far below the CMB temperature of 2.725 K. Hawking radiation dominates over CMB absorption only for black holes with $M \lesssim 4.5 \times 10^{22}$ kg ($\approx 61\%$ of the lunar mass $M_{\text{Moon}} = 7.34 \times 10^{22}$ kg), corresponding to $T_H \gtrsim 2.7$ K.

Physical interpretation in the ether picture. The radiation arises because the ether's ZPF modes near the horizon are split by the flow. A virtual fluctuation of the ZPF — a mode that, in the absence of the flow, would remain part of the ground state — is torn apart by the velocity gradient at the horizon. One component is swept inward by the ether flow (absorbed by the black hole); the other escapes outward as real radiation. The energy for this process comes from the ether's kinetic energy of inflow — equivalently, from the gravitational binding energy of the black hole, which decreases as the mass is radiated away.

This is the ether's concrete realisation of the "virtual pair creation at the horizon" picture that is often invoked heuristically. In the ether framework, the "pairs" are modes of the ZPF; the "creation" is the flow-induced mode-splitting at the acoustic horizon; and the energy source is the ether inflow, which decelerates as the black hole loses mass.

3.13.7 Extension to Rotating Black Holes

The Kerr–ether identity (Theorem 3.4) provides the velocity field for a rotating black hole. The outer horizon is at $r_+ = r_s/2 + \sqrt{r_s^2/4 - a^2}$, where $a = J/(Mc)$.

Surface gravity for Kerr. We derive $\kappa_{\text{Kerr}}$ from the exponential peeling rate of outgoing null rays in Doran coordinates, following the same method as Section 3.13.1 for Schwarzschild.

In Doran coordinates (the Kerr analogue of PG), the radial part of the ether velocity is $V_r = -c\beta(r)$ with (from Eq 3.71):

$$\beta(r) = \sqrt{\frac{r_s r}{r^2+a^2}} \tag{3.222a}$$

The outgoing principal null ray in Doran coordinates satisfies (from the metric (3.67) with $d\theta = d\tilde{\phi} = 0$ along the principal null congruence):

$$\frac{dr}{dT} = c(1 - \beta) \tag{3.222b}$$

Derivation of (3.222b). The Doran metric (3.67) restricted to $d\theta = 0$, $d\tilde\phi = 0$ gives $ds^2 = -c^2 dT^2 + (dr + \beta c\,dT)^2$. Setting $ds^2 = 0$: $(dr + \beta c\,dT)^2 = c^2 dT^2$, so $dr = (-\beta \pm 1)c\,dT$. The outgoing ray selects the $+$ sign (since $\beta < 1$ outside the horizon gives $dr/dT > 0$).

The horizon condition. At $r = r_+$: $\beta(r_+) = 1$, i.e., $r_s r_+ = r_+^2 + a^2$. This is equivalent to $\Delta(r_+) = r_+^2 - r_sr_+ + a^2 = 0$, the standard Kerr horizon condition, with solutions $r_\pm = r_s/2 \pm \sqrt{r_s^2/4 - a^2}$, giving $r_+ + r_- = r_s$ and $r_+r_- = a^2$.

Near-horizon expansion. Expand $\beta(r)$ to first order around $r_+$. From $\beta^2 = r_sr/(r^2+a^2)$:

$$\frac{d(\beta^2)}{dr} = \frac{r_s(r^2+a^2) - r_sr\cdot 2r}{(r^2+a^2)^2}$$ $$= \frac{r_s(a^2 - r^2)}{(r^2+a^2)^2} \tag{3.222c}$$

At $r = r_+$, using $\beta(r_+) = 1$:

$$\beta'(r_+) = \frac{1}{2\beta(r_+)}\cdot\frac{d(\beta^2)}{dr}\bigg|_{r_+}$$ $$= \frac{r_s(a^2 - r_+^2)}{2(r_+^2+a^2)^2} \tag{3.222d}$$

Simplify the numerator and denominator using the Kerr identities $r_+^2 - a^2 = r_+(r_+ - r_-)$ (from $r_+r_- = a^2$) and $r_+^2 + a^2 = r_+r_s$ (from $r_sr_+ = r_+^2 + a^2$):

$$\beta'(r_+) = \frac{-r_s \cdot r_+(r_+ - r_-)}{2(r_+r_s)^2} = \frac{-(r_+ - r_-)}{2r_+r_s} \tag{3.222e}$$

The outgoing null ray velocity near $r_+$ is:

$$\frac{dr}{dT} = c(1-\beta) \approx -c\,\beta'(r_+)(r - r_+)$$ $$= \frac{c(r_+ - r_-)}{2r_+r_s}(r - r_+) \tag{3.222f}$$

This has the same exponential form as (3.193): $\delta r(T) = \delta r_0\,\exp(\lambda T)$ with peeling rate $\lambda = c(r_+ - r_-)/(2r_+r_s)$.

The surface gravity. By the same definition as (3.195) ($\kappa = c\lambda$):

$$\kappa_{\text{Kerr}} = \frac{c^2(r_+ - r_-)}{2r_+r_s}$$

Using $r_+r_s = r_+^2 + a^2$:

$$\kappa_{\text{Kerr}} = \frac{c^2(r_+ - r_-)}{2(r_+^2 + a^2)} \tag{3.223}$$

where $r_- = r_s/2 - \sqrt{r_s^2/4 - a^2}$ is the inner horizon.

Verification of limits. For $a = 0$: $r_+ = r_s$, $r_- = 0$, $r_+^2 + a^2 = r_s^2$. (3.223) gives $\kappa = c^2 r_s/(2r_s^2) = c^2/(2r_s)$, recovering (3.195).

For the extremal limit $a \to r_s/2$ ($J \to GM^2/c$): $r_+ \to r_- \to r_s/2$, so $\kappa_{\text{Kerr}} \to 0$. An extremal black hole has zero surface gravity and zero Hawking temperature — it does not radiate.

The Hawking temperature for Kerr:

$$T_H^{\text{Kerr}} = \frac{\hbar\,\kappa_{\text{Kerr}}}{2\pi\,k_B\,c}$$ $$= \frac{\hbar\,c\,(r_+ - r_-)}{4\pi\,k_B\,(r_+^2 + a^2)} \tag{3.224}$$

Superradiance. The Kerr ether flow has nonzero azimuthal velocity $V_{\tilde{\phi}}$ (Eq 3.72), which produces the frame-dragging angular velocity $\omega_H = r_s a c/(2(r_+^2 + a^2))$ at the horizon (Eq 3.78 evaluated at $r = r_+$). Modes with azimuthal quantum number $m$ and frequency $\omega$ satisfying $\omega < m\,\omega_H$ are amplified rather than thermally emitted — the ether flow at the horizon co-rotates faster than the mode's phase velocity, extracting rotational energy. The condition $\omega < m\omega_H$ is the superradiance condition. In the ether picture, it corresponds to modes that are swept forward by the azimuthal ether flow faster than they propagate on their own. This is the rotational analog of the Penrose process, realised through the ether's vortical flow rather than through ergoregion orbits.

3.13.8 The Trans-Planckian Problem and the Ether's Resolution

The problem. The mode function (3.210) oscillates with unbounded frequency as $r \to r_s^+$: the local frequency measured by a static observer at radius $r$ diverges as $(c - v)^{-1} \sim (r - r_s)^{-1}$. For a mode of frequency $\omega$ at infinity, the near-horizon frequency is:

$$\omega_{\text{local}} = \frac{\omega}{c - v} \approx \frac{\omega\,c}{\kappa\,(r - r_s)} \tag{3.225}$$

At a proper distance $\delta\ell = \sqrt{g_{rr}}\,\delta r = \delta r$ (in PG coordinates, $g_{rr} = 1$) of one Planck length $\ell_P = 1.616 \times 10^{-35}$ m from the horizon:

$$\omega_{\text{local}} \sim \frac{\omega\,c}{\kappa\,\ell_P} \sim \omega \cdot \frac{r_s}{\ell_P}$$ $$\sim \omega \cdot 10^{38}\left(\frac{M}{M_\odot}\right) \tag{3.226}$$

For any astrophysical black hole, the mode that arrives at infinity with frequency $\omega$ originated at the horizon with a frequency exceeding the Planck frequency by tens of orders of magnitude. The derivation of §Section 3.13.3–3.13.5 assumed that the wave (3.202) holds at all these frequencies — but no known physics operates at the Planck frequency. This is the trans-Planckian problem of Hawking radiation [172].

The ether's response. The ether framework provides a concrete UV structure that modifies the dispersion relation at high frequencies. At the ether's transverse microstructure scale $\ell_e$, the electromagnetic dispersion relation is modified (Eq 3.46):

$$\omega^2 = c^2 k^2\!\left(1 + \xi_2(k\ell_e)^2 + \cdots\right) \tag{3.227}$$

For $k\ell_e \sim 1$ (wavelengths comparable to the ether microstructure), the group velocity deviates from $c$ and the geometric optics approximation breaks down. The modes do not blueshift to arbitrarily high frequencies; instead, they reach the ether's UV cutoff and the derivation of Section 3.13.3 must be modified.

Jacobson's universality argument. Jacobson [173] and subsequently Unruh [174] demonstrated through explicit numerical calculations that the thermal spectrum is insensitive to the UV modification of the dispersion relation, provided two conditions are satisfied:

(i) The modification occurs at a frequency scale $\omega_{\text{UV}} \gg \kappa/c$ (the UV cutoff is far above the Hawking frequency).

(ii) The modified dispersion relation connects smoothly to the standard $\omega = ck$ relation at low frequencies.

Both conditions are satisfied by the ether. Condition (i): the Hawking frequency is $\omega_H \sim \kappa/c = c/(2r_s) \sim 10^4$ Hz for a solar-mass black hole, while the ether UV cutoff is $\omega_{\text{UV}} \sim c/\ell_e > 10^{21}$ Hz (from the observational constraint $\ell_e < 6 \times 10^{-13}$ m, Eq 3.50), giving $\omega_{\text{UV}}/\omega_H > 10^{17}$. Condition (ii): the ether dispersion (3.227) reduces to $\omega = ck$ for $k\ell_e \ll 1$ by construction.

The physical mechanism behind this universality is that the Hawking effect depends only on the near-horizon geometry (through $\kappa$) and the low-energy mode structure (through the dispersion relation at $\omega \sim \kappa/c$). The trans-Planckian modes, while formally present in the derivation, contribute only through their low-energy descendants — the modes that have been redshifted down to the Hawking frequency by the time they escape to infinity. The UV details of the ether microstructure are washed out by the exponential redshift.

The ether's advantage. In the standard framework, the trans-Planckian problem is a genuine conceptual difficulty: the derivation relies on a wave equation that must hold at arbitrarily high energies, where the known laws of physics break down. In the ether framework, this difficulty is resolved in principle: the ether has a physical UV cutoff ($\ell_e$ for the EM sector, $\xi$ for the phonon sector), and the Jacobson universality argument guarantees that the Hawking result survives this cutoff. The ether provides the UV completion that the standard derivation lacks.

3.13.9 Summary

Result	Derivation	Key equation
Surface gravity from ether flow	Velocity gradient at $v = c$	(3.195)
Near-horizon mode singularity	Retarded time divergence	(3.210)
Bogoliubov coefficient ratio	Analytic continuation	(3.215)
Planckian spectrum	Bose–Einstein statistics	(3.218)
Hawking temperature	$T_H = \hbar\kappa/(2\pi k_Bc)$	(3.220)
Kerr extension	$\kappa_{\text{Kerr}} = c^2(r_+ - r_-)/(2(r_+^2+a^2))$	(3.224)
Trans-Planckian resolution	Ether UV cutoff + Jacobson universality	Section 3.13.8

The ether predicts Hawking radiation as a direct consequence of its flow structure: the acoustic horizon (where $v = c$) combined with the ether's zero-point fluctuations (§6) produces a thermal flux at infinity. The derivation uses only the ether's velocity profile, the covariant wave equation on the ether metric, and the standard theory of Bogoliubov transformations. No result from quantum field theory on curved spacetime is imported; the Hawking effect is an output of the ether framework, not an input.

3.14 Gravitational Wave Polarisations

Section 3.7 identified a potential scalar (breathing) polarisation mode in the ether's gravitational wave spectrum. Theorem 3.5 resolves this question: the ether's field equation is the Einstein equation, whose constraint structure eliminates all but two propagating degrees of freedom. This section derives the constraint structure explicitly from the linearised Einstein equation applied to the ether metric perturbation, computes every component, and shows that the scalar breathing mode is non-radiative.

3.14.1 The Linearised Einstein Equation

The linearised Einstein tensor for a perturbation $h_{\mu\nu}$ around flat spacetime $\eta_{\mu\nu} = \text{diag}(-c^2, 1, 1, 1)$ is the standard result of linearised GR (the standard linearisation of the Einstein tensor):

$$2G^{(1)}_{\mu\nu} = -\Box h_{\mu\nu} - \eta_{\mu\nu}\partial_\alpha\partial_\beta h^{\alpha\beta} + \eta_{\mu\nu}\Box h$$ $$+ \partial_\mu\partial_\alpha h^\alpha{}_\nu + \partial_\nu\partial_\alpha h^\alpha{}_\mu - \partial_\mu\partial_\nu h \tag{3.228}$$

where $\Box = \eta^{\alpha\beta}\partial_\alpha\partial_\beta = -c^{-2}\partial_T^2 + \nabla^2$ is the flat-space d'Alembertian, $h = \eta^{\mu\nu}h_{\mu\nu}$ is the trace, and all indices are raised and lowered with $\eta_{\mu\nu}$. (3.228) is a mathematical identity — the first-order Taylor expansion of $G_{\mu\nu}[\eta + h]$ — and holds for ANY perturbation regardless of the physical theory.

By Theorem 3.5, the vacuum field equation is $G_{\mu\nu} = 0$. At first order: $G^{(1)}_{\mu\nu} = 0$.

3.14.2 The Ether Perturbation in PG Form

The ether metric (3.92) around a flat background ($V^i_0 = 0$) with perturbation $\delta V^i$ has:

$$h_{00} = 0, \qquad h_{0i} = -\delta V_i, \qquad h_{ij} = 0 \tag{3.229}$$

(since $g_{00} = -(c^2 - V^2)$ gives $h_{00} = V^2 = O(\delta V^2)$, and $g_{ij} = \delta_{ij}$ gives $h_{ij} = 0$). Decompose the velocity as in Section 3.7:

$$\delta V^i = \partial^i\varphi + \xi^i, \qquad \partial_i\xi^i = 0 \tag{3.230}$$

where $\varphi$ is the scalar potential (one degree of freedom) and $\xi^i$ is the transverse part (two degrees of freedom). Define $\theta = \partial_i\delta V^i = \nabla^2\varphi$.

Index gymnastics. We compute all contractions needed for (3.228).

The trace: $h = \eta^{\mu\nu}h_{\mu\nu} = \eta^{00}h_{00} + \eta^{ij}h_{ij} = 0$. Therefore $\Box h = 0$ and $\partial_\mu\partial_\nu h = 0$.

Mixed components: $h^\alpha{}_\nu = \eta^{\alpha\beta}h_{\beta\nu}$:

$$h^0{}_0 = \eta^{00}h_{00} = 0 \tag{3.231a}$$ $$h^i{}_0 = \eta^{ij}h_{j0} = \delta^{ij}(-\delta V_j) = -\delta V^i \tag{3.231b}$$ $$h^0{}_i = \eta^{00}h_{0i} = \frac{\delta V_i}{c^2} \tag{3.231c}$$ $$h^j{}_i = \eta^{jk}h_{ki} = 0 \tag{3.231d}$$

Divergences: $\partial_\alpha h^\alpha{}_\nu$:

$$\partial_\alpha h^\alpha{}_0 = \partial_T h^0{}_0 + \partial_j h^j{}_0$$ $$= 0 - \partial_j\delta V^j = -\theta \tag{3.232a}$$ $$\partial_\alpha h^\alpha{}_i = \partial_T h^0{}_i + \partial_j h^j{}_i$$ $$= \frac{\dot{\delta V_i}}{c^2} + 0 = \frac{\dot{\delta V_i}}{c^2} \tag{3.232b}$$

where $\dot{} = \partial_T$. Fully raised divergence: $\partial_\alpha h^{\alpha\beta}$:

$$\partial_\alpha h^{\alpha 0} = \partial_T h^{00} + \partial_j h^{j0}$$ $$= 0 + \partial_j(\delta V^j/c^2) = \theta/c^2 \tag{3.232c}$$

(using $h^{j0} = \eta^{jk}\eta^{00}h_{k0} = \delta^{jk}\cdot(-1/c^2)\cdot(-\delta V_k) = \delta V^j/c^2$).

$$\partial_\alpha h^{\alpha i} = \partial_T h^{0i} + \partial_j h^{ji}$$ $$= \frac{\dot{\delta V^i}}{c^2} \tag{3.232d}$$

(using $h^{0i} = \eta^{00}\eta^{ij}h_{0j} = (-1/c^2)\delta^{ij}(-\delta V_j) = \delta V^i/c^2$).

Double divergence:

$$\partial_\alpha\partial_\beta h^{\alpha\beta} = \partial_0(\theta/c^2) + \partial_i(\dot{\delta V^i}/c^2)$$ $$= \frac{\dot{\theta}}{c^2} + \frac{\dot{\theta}}{c^2} = \frac{2\dot{\theta}}{c^2} \tag{3.233}$$

3.14.3 The Three Einstein Equations

The $(00)$ equation: $G^{(1)}_{00} = 0$. Substituting into (3.228) with $h_{00} = 0$, $h = 0$:

$$2G^{(1)}_{00} = -\Box\cdot 0 - \eta_{00}\cdot\frac{2\dot{\theta}}{c^2} + 0 + 2\partial_T(-\theta) - 0$$ $$= \frac{2c^2\dot{\theta}}{c^2} - 2\dot{\theta} = 2\dot{\theta} - 2\dot{\theta} = 0 \tag{3.234}$$

The $(00)$ equation is satisfied identically. This is the Hamiltonian constraint at first order: it imposes no condition on the perturbation because $h_{00} = h = 0$ in PG gauge.

The $(0i)$ equation: $G^{(1)}_{0i} = 0$. With $\eta_{0i} = 0$:

$$2G^{(1)}_{0i} = -\Box h_{0i} + \partial_T\!\left(\frac{\dot{\delta V_i}}{c^2}\right) + \partial_i(-\theta) \tag{3.235a}$$

Now $h_{0i} = -\delta V_i$, so $\Box h_{0i} = \Box(-\delta V_i) = -\Box\delta V_i = (c^{-2}\ddot{\delta V_i} - \nabla^2\delta V_i)$. Substituting:

$$2G^{(1)}_{0i} = -c^{-2}\ddot{\delta V_i} + \nabla^2\delta V_i + c^{-2}\ddot{\delta V_i} - \partial_i\theta$$

The $\ddot{\delta V_i}$ terms cancel:

$$2G^{(1)}_{0i} = \nabla^2\delta V_i - \partial_i\theta$$ $$= \nabla^2\delta V_i - \partial_i(\partial_j\delta V^j) \tag{3.235b}$$

By the vector identity $\nabla^2\mathbf{A} - \nabla(\nabla\cdot\mathbf{A}) = -\nabla\times(\nabla\times\mathbf{A})$:

$$G^{(1)}_{0i} = 0 \quad \Longleftrightarrow \quad \nabla\times(\nabla\times\delta\mathbf{V}) = 0 \tag{3.235c}$$

This is the momentum constraint. Substituting the decomposition (3.230) and using $\nabla\times\nabla\varphi = 0$:

$$\nabla\times(\nabla\times\boldsymbol{\xi}) = -\nabla^2\boldsymbol{\xi} = 0 \tag{3.236}$$

(using $\nabla\times(\nabla\times\boldsymbol{\xi}) = \nabla(\nabla\cdot\boldsymbol{\xi}) - \nabla^2\boldsymbol{\xi} = -\nabla^2\boldsymbol{\xi}$ since $\nabla\cdot\boldsymbol{\xi} = 0$). For perturbations decaying at spatial infinity: $\boldsymbol{\xi} = 0$.

The momentum constraint eliminates the transverse velocity at each instant. Equivalently: $\nabla^2\delta V_i = \partial_i\theta$, which means $\delta V^i$ is a pure gradient: $\delta V^i = \partial^i\varphi$ with $\nabla^2\varphi = \theta$.

The $(ij)$ equation: $G^{(1)}_{ij} = 0$. With $h_{ij} = 0$ and $h = 0$:

$$2G^{(1)}_{ij} = 0 - \delta_{ij}\cdot\frac{2\dot{\theta}}{c^2} + 0$$ $$+ \partial_i\!\left(\frac{\dot{\delta V_j}}{c^2}\right) + \partial_j\!\left(\frac{\dot{\delta V_i}}{c^2}\right) - 0$$ $$= \frac{1}{c^2}\!\left[\partial_i\dot{\delta V_j} + \partial_j\dot{\delta V_i} - 2\delta_{ij}\dot{\theta}\right] \tag{3.237a}$$

Setting this to zero:

$$\partial_i\dot{\delta V_j} + \partial_j\dot{\delta V_i} = 2\delta_{ij}\dot{\theta} \tag{3.237b}$$

Trace of (3.237b). Contract with $\delta^{ij}$: $2\dot{\theta} = 6\dot{\theta}$, therefore:

$$\dot{\theta} = 0 \tag{3.237c}$$

The velocity divergence is time-independent. Combined with the boundary condition $\theta \to 0$ at spatial infinity and the constraint $\nabla^2\delta V_i = \partial_i\theta$ from (3.235b): $\theta$ is a time-independent harmonic function ($\nabla^2\theta = \nabla^2(\nabla^2\varphi) = \nabla^4\varphi$, which is determined by $\nabla^2\delta V_i = \partial_i\theta$ plus boundary conditions). In vacuum with no sources and decaying boundary conditions: $\theta = 0$.

With $\theta = 0$, (3.237b) becomes:

$$\partial_i\dot{\delta V_j} + \partial_j\dot{\delta V_i} = 0 \tag{3.237d}$$

The velocity perturbation is a Killing field of flat space. (3.237d) says $\dot{\delta V_i}$ is a Killing vector: $\partial_{(i}\dot{\delta V_{j)}} = 0$. For a perturbation decaying at infinity, the only Killing vectors of flat space that decay are identically zero. Therefore:

$$\dot{\delta V_i} = 0 \tag{3.238}$$

The PG velocity perturbation is static. The linearised Einstein equation forces the ether velocity perturbation (in PG form) to be time-independent. This is the static Newtonian potential — not a gravitational wave.

3.14.4 Gravitational Waves Require Departing from Strict PG Form

The result (3.238) means that the PG metric ansatz (unit lapse, flat spatial slices) is too restrictive for radiative solutions. The PG form $\gamma_{ij} = \delta_{ij}$ is exact for the stationary Schwarzschild and Kerr solutions but does not accommodate GW, which require oscillating spatial curvature.

In the ether picture, a gravitational wave produces a perturbation of the spatial metric $\gamma_{ij} = \delta_{ij} + h_{ij}$ alongside the velocity field. The full metric perturbation around flat space is therefore:

$$h_{00} = 0, \qquad h_{0i} = -\delta V_i, \qquad h_{ij} = H_{ij} \tag{3.239}$$

where $H_{ij}$ is a small perturbation of the spatial geometry. The ten metric functions are: $h_{00}$ (one, zero here), $h_{0i}$ (three, the velocity), and $H_{ij}$ (six, symmetric spatial tensor). Total: ten, with $h_{00} = 0$ already imposed.

We now derive the constraint and gauge structure that reduces nine functions to two propagating degrees of freedom.

3.14.5 The Lorenz Gauge and Residual Freedom

The trace-reversed perturbation. Define $\bar{h}_{\mu\nu} = h_{\mu\nu} - \frac{1}{2}\eta_{\mu\nu}h$, where $h = \eta^{\mu\nu}h_{\mu\nu} = -(0)/c^2 + H^i{}_i = H$, the trace of the spatial perturbation. The linearised Einstein (3.228) in terms of $\bar{h}$ takes the compact form (the standard trace-reversal identity):

$$G^{(1)}_{\mu\nu} = \frac{1}{2}\Big[-\Box\bar{h}_{\mu\nu} + \partial_\mu\partial_\alpha\bar{h}^\alpha{}_\nu$$ $$+ \partial_\nu\partial_\alpha\bar{h}^\alpha{}_\mu - \eta_{\mu\nu}\partial_\alpha\partial_\beta\bar{h}^{\alpha\beta}\Big] \tag{3.240}$$

The Lorenz gauge. Under an infinitesimal coordinate transformation $x^\mu \to x^\mu + \xi^\mu$:

$$h_{\mu\nu} \to h_{\mu\nu} - \partial_\mu\xi_\nu - \partial_\nu\xi_\mu \tag{3.241}$$

where $\xi_\mu = \eta_{\mu\nu}\xi^\nu$. The trace-reversed perturbation transforms as $\bar{h}_{\mu\nu} \to \bar{h}_{\mu\nu} - \partial_\mu\xi_\nu - \partial_\nu\xi_\mu + \eta_{\mu\nu}\partial_\alpha\xi^\alpha$. The divergence transforms as:

$$\partial^\mu\bar{h}_{\mu\nu} \to \partial^\mu\bar{h}_{\mu\nu} - \Box\xi_\nu \tag{3.242}$$

Choosing $\xi_\nu$ to satisfy $\Box\xi_\nu = \partial^\mu\bar{h}_{\mu\nu}$ (which always has a solution), we can impose the Lorenz gauge condition:

$$\partial^\mu\bar{h}_{\mu\nu} = 0 \tag{3.243}$$

This is four equations (one for each $\nu$), eliminating four of the nine free functions. In Lorenz gauge, (3.240) reduces to:

$$\Box\bar{h}_{\mu\nu} = 0 \tag{3.244}$$

Every component of $\bar{h}_{\mu\nu}$ satisfies the wave equation. The nine remaining functions (after $h_{00} = 0$) minus four Lorenz gauge conditions leaves five.

Residual gauge freedom. The Lorenz condition (3.243) is preserved by any further transformation $\xi^\mu$ satisfying $\Box\xi^\mu = 0$. This gives four free functions (one for each component of $\xi^\mu$), each satisfying the wave equation. These can be used to impose four additional conditions on $\bar{h}_{\mu\nu}$. Five functions minus four residual gauges leaves one — but this counting overcounts because not all gauge conditions are independent of the Lorenz condition.

The standard counting proceeds as follows from the general structure: starting from ten symmetric components of $h_{\mu\nu}$, the Lorenz gauge eliminates four, leaving six. Residual gauge freedom ($\Box\xi^\mu = 0$, four functions) eliminates four more, leaving two. These two are the physical propagating degrees of freedom.

3.14.6 The Transverse-Traceless Gauge

The residual gauge freedom is sufficient to impose the TT gauge conditions:

$$\bar{h}_{0\mu} = 0, \qquad \bar{h}^i{}_i = 0, \qquad \partial_j\bar{h}^{j}{}_i = 0 \tag{3.245}$$

These are: three conditions from $\bar{h}_{0i} = 0$, one from the trace condition, and three from the divergence condition — but the divergence conditions are not independent of the Lorenz gauge (they are the spatial part of (3.243)), and the $\bar{h}_{00}$ condition uses one residual gauge function. The net effect: six components of $H_{ij}$, minus one trace, minus three transverse conditions = two.

Explicit construction. For a plane wave propagating in the $z$-direction ($k^\mu = (\omega/c, 0, 0, k)$ with $\omega = ck$), the TT gauge perturbation has the form:

$$h^{TT}_{ij} = \begin{pmatrix} A_+ & A_\times & 0 \\ A_\times & -A_+ & 0 \\ 0 & 0 & 0 \end{pmatrix} \cos(\omega T - kz) \tag{3.246}$$

where $A_+$ and $A_\times$ are the amplitudes of the plus and cross polarisations respectively. The matrix is traceless ($A_+ - A_+ + 0 = 0$), transverse ($k^j h^{TT}_{jx} = k^j h^{TT}_{jy} = 0$ since $k$ is along $z$ and $h^{TT}_{zj} = 0$), and satisfies $\Box h^{TT}_{ij} = 0$.

Verification of completeness. These two polarisations exhaust the physical content. Consider the most general symmetric, traceless, transverse $3\times 3$ matrix for a wave propagating along $\hat{z}$. Tracelessness: $h_{xx} + h_{yy} + h_{zz} = 0$. Transversality: $h_{xz} = h_{yz} = h_{zz} = 0$. The first two from transversality, combined with tracelessness ($h_{xx} + h_{yy} = 0$), leave two free components: $h_{xx} = -h_{yy}$ and $h_{xy}$. These are $A_+$ and $A_\times$.

3.14.7 The Scalar Breathing Mode Is Non-Radiative

We now prove that the scalar mode identified in Section 3.7 does not radiate.

The scalar mode. In the ether picture, the scalar perturbation is the velocity divergence $\theta = \partial_i\delta V^i = \nabla^2\varphi$, or equivalently the density perturbation $\delta\rho$ (connected to $\theta$ by the continuity (3.39)). A breathing mode would be a spherically symmetric oscillation of the spatial metric: $H_{ij} = H(T, r)\,\delta_{ij}$.

Proof that $H\delta_{ij}$ is non-radiative. Suppose $h_{\mu\nu}$ has only a trace part: $h_{ij} = H\,\delta_{ij}$, $h_{0\mu} = 0$. Then $h = 3H$, $\bar{h}_{ij} = H\delta_{ij} - \frac{1}{2}\delta_{ij}\cdot 3H = -\frac{1}{2}H\delta_{ij}$, and $\bar{h}_{00} = 0 + \frac{c^2}{2}\cdot 3H = \frac{3c^2H}{2}$.

The Lorenz gauge condition (3.243) for $\nu = 0$:

$$\partial^\mu\bar{h}_{\mu 0} = \eta^{00}\partial_T\bar{h}_{00}$$ $$= -\frac{1}{c^2}\partial_T\!\left(\frac{3c^2H}{2}\right) = -\frac{3\dot{H}}{2} = 0 \tag{3.247}$$

Therefore $\dot{H} = 0$: the scalar mode is time-independent in Lorenz gauge. A mode that does not oscillate cannot radiate. $\square$

Physical interpretation. The scalar perturbation $H\delta_{ij}$ represents a uniform expansion or contraction of the spatial geometry. The linearised Einstein equation requires this expansion to be static — it describes a Newtonian-like potential ($H = -2\Phi/c^2$ in the weak-field limit), not a wave. The constraint $\dot{H} = 0$ is the linearised form of the Hamiltonian constraint: the spatial expansion rate is fixed by the energy density, not by initial data, and it cannot propagate independently.

This resolves the discrepancy with Section 3.7. The wave (3.42) for $\delta\rho$ describes perturbations of the ether density, but the gravitational constraint couples $\delta\rho$ algebraically to the source. The free (vacuum) solution has $\delta\rho = 0$ — the density perturbation does not propagate as a gravitational wave in vacuum. It responds instantaneously to the source distribution, exactly like the Coulomb field in electrodynamics.

3.14.8 The Isaacson Energy Flux

The energy carried by gravitational waves is given by the Isaacson effective stress-energy tensor (Isaacson 1968), which is constructed from the second-order perturbation averaged over several wavelengths:

$$T^{\text{GW}}_{\mu\nu} = \frac{c^4}{32\pi G}\!\left\langle\partial_\mu h^{TT}_{ij}\,\partial_\nu h^{TT\,ij}\right\rangle \tag{3.248}$$

The energy flux (power per unit area) in the propagation direction $\hat{z}$ is:

$$F = -c\,T^{\text{GW}}_{0z}$$ $$= \frac{c^3}{32\pi G}\!\left\langle\dot{h}^{TT}_{ij}\,\partial_z h^{TT\,ij}\right\rangle \tag{3.249a}$$

For the plane wave (3.246) with $\dot{h}^{TT}_{ij} = -\omega\sin(\omega T - kz)\cdot A_{ij}$ and $\partial_z h^{TT}_{ij} = k\sin(\omega T - kz)\cdot A_{ij}$, and using $\langle\sin^2\rangle = 1/2$:

$$F = \frac{c^3\omega k}{32\pi G}\cdot\frac{1}{2}(A_{ij}A^{ij})$$ $$= \frac{c^2\omega^2}{32\pi G}(A_+^2 + A_\times^2) \tag{3.249b}$$

where $A_{ij}A^{ij} = 2A_+^2 + 2A_\times^2$ (the factor of 2 for each polarisation from $h_{xx}^2 + h_{yy}^2 + 2h_{xy}^2 = 2A_+^2 + 2A_\times^2$) and $\omega = ck$.

The scalar mode carries no flux. For a pure scalar perturbation $h_{ij} = H(T,r)\delta_{ij}$ with $\dot{H} = 0$ (from Section 3.14.7): $\dot{h}_{ij} = 0$, so $T^{\text{GW}}_{0z} = 0$. The energy flux vanishes identically. The scalar mode does not radiate. $\square$

The tensor modes carry positive flux. For $A_+^2 + A_\times^2 > 0$: $F > 0$. The two TT modes carry energy to infinity. They are the only radiative degrees of freedom.

3.14.9 Gravitational Wave Polarisation Content

3.14.10 Summary of Polarisation Results

The analysis of §Section 3.14.1–3.14.9 establishes that the breathing mode exists as a static (non-radiative) constrained field, analogous to the Coulomb field in electrodynamics. The complete polarisation results are:

(i) Gravitational waves have exactly the GR polarisation content: two tensor modes ($+$ and $\times$), zero scalar modes, zero vector modes.

(ii) Current LIGO-Virgo-KAGRA observations, which constrain non-GR polarisations [48, 180], are consistent with this prediction.

(iii) At the ether's UV scale $\ell_e$, dispersive corrections (Eq 3.46) may modify the polarisation amplitudes. These are suppressed by $(f\ell_e/c)^2$ and are undetectable for any foreseeable observatory.

3.15 The Unruh Effect from the Ether

The Hawking effect (Section 3.13) arises from the ether's gravitational horizon — the surface where the ether inflow velocity equals $c$. The Unruh effect [181] is the flat-spacetime counterpart: an observer accelerating through the ether's zero-point field detects thermal radiation, even though no horizon exists in the background spacetime.

The Unruh effect is a direct consequence of the ether's ZPF structure (§6) combined with the equivalence principle: an accelerating observer in flat spacetime is locally equivalent to a static observer in a gravitational field, and therefore experiences the ether's quantum fluctuations as thermal. We derive the Unruh temperature from the ether's mode structure without invoking quantum field theory on curved spacetime.

3.15.1 The Rindler Horizon

An observer accelerating uniformly at proper acceleration $a$ in flat spacetime follows the worldline:

$$T(\tau) = \frac{c}{a}\sinh\!\left(\frac{a\tau}{c}\right)$$ $$x(\tau) = \frac{c^2}{a}\cosh\!\left(\frac{a\tau}{c}\right) \tag{3.250}$$

where $\tau$ is the observer's proper time and the motion is along the $x$-axis.

Derivation. The four-acceleration has magnitude $a^\mu a_\mu = a^2/c^4$ (in the convention where $a^\mu = D u^\mu/d\tau$). For motion along $x$ with $u^\mu = \gamma(c, v, 0, 0)$, the constraint $u_\mu u^\mu = -c^2$ and $a_\mu u^\mu = 0$ give $d\gamma/d\tau = a\gamma v/c^3$ and $dv/d\tau = a/\gamma^3$. With initial conditions $v(0) = 0$, $x(0) = c^2/a$: $v = c\tanh(a\tau/c)$, $\gamma = \cosh(a\tau/c)$, yielding (3.250) by integration.

The Rindler horizon. The accelerating observer cannot receive signals from the region $x < c\,T$ (the left Rindler wedge). The boundary $x = c\,T$ is the Rindler horizon — a null surface in flat spacetime that is causally inaccessible to the accelerating observer.

Connection to the ether picture. In the ether's rest frame, the observer accelerates through the stationary ZPF. From the observer's perspective, the ZPF is Doppler-shifted by the time-varying velocity $v(\tau) = c\tanh(a\tau/c)$. The Rindler horizon is the surface from which ZPF modes would need to travel at $c$ against the observer's acceleration to reach her — the flat-spacetime analog of the acoustic horizon in Section 3.13.

3.15.2 Rindler Coordinates and the Ether Mode Structure

Define Rindler coordinates $(\eta, \rho)$ adapted to the accelerating observer:

$$T = \frac{\rho}{c}\sinh\!\left(\frac{a\eta}{c}\right), \qquad x = \rho\,\cosh\!\left(\frac{a\eta}{c}\right) \tag{3.251}$$

where $\eta$ is the Rindler time (proportional to the observer's proper time at $\rho = c^2/a$: $d\tau = (a\rho/c^2)\,d\eta$) and $\rho > 0$ is the Rindler spatial coordinate. The Minkowski metric in Rindler coordinates is:

$$ds^2 = -\frac{a^2\rho^2}{c^2}\,d\eta^2 + d\rho^2 + dy^2 + dz^2 \tag{3.252}$$

Derivation. From (3.251): $dT = (\rho a/c^2)\cosh(a\eta/c)\,d\eta + (1/c)\sinh(a\eta/c)\,d\rho$ and $dx = (\rho a/c)\sinh(a\eta/c)\,d\eta + \cosh(a\eta/c)\,d\rho$. Computing $-c^2dT^2 + dx^2$: the $d\eta^2$ terms give $\rho^2a^2(-\cosh^2 + \sinh^2)/c^2 = -\rho^2a^2/c^2$; the $d\rho^2$ terms give $(-\sinh^2 + \cosh^2) = 1$; the cross-terms cancel. This gives (3.252).

The Rindler horizon is at $\rho = 0$. The metric component $g_{\eta\eta} = -a^2\rho^2/c^2$ vanishes there, exactly as $g_{TT} = -(c^2 - v^2)$ vanishes at the ether horizon in Section 3.13 (Eq 3.36). The surface gravity of the Rindler horizon is:

$$\kappa_R = \lim_{\rho \to 0}\frac{c^2}{2}\,\frac{\partial_\rho|g_{\eta\eta}|}{|g_{\eta\eta}|^{1/2}}$$ $$\times \frac{1}{(-g_{\eta\eta})^{1/2}} \tag{3.253a}$$

More directly: the proper acceleration at Rindler coordinate $\rho$ is $a_{\text{local}} = c^2/(a\rho/c^2)^{-1} = a$ at $\rho = c^2/a$, and the redshift factor between $\rho$ and infinity diverges at $\rho = 0$. The surface gravity, defined as the exponential peeling rate of null rays from the horizon (as in Section 3.13.1), is computed as follows.

Outgoing null rays in Rindler coordinates satisfy $ds^2 = 0$:

$$\frac{d\rho}{d\eta} = \frac{a\rho}{c} \tag{3.253b}$$

This has the exponential solution $\rho(\eta) = \rho_0\,e^{a\eta/c}$, with peeling rate $a/c$. The surface gravity is:

$$\kappa_R = a \tag{3.253c}$$

This is exact — no approximation to a near-horizon limit is needed, because the Rindler metric (3.252) already has the exact exponential structure everywhere.

3.15.3 Mode Analysis and the Thermal Spectrum

The derivation of the thermal spectrum proceeds identically to Section 3.13.3–3.13.5, with the ether's gravitational horizon replaced by the Rindler horizon.

Retarded time. For outgoing modes, the retarded coordinate is $u = \eta - (c/a)\ln\rho$, obtained by integrating $d\eta = (c/(a\rho))\,d\rho$ from (3.253b). A mode of Rindler frequency $\omega$ has the form:

$$\psi_\omega(\eta, \rho) = e^{-i\omega u} = e^{-i\omega\eta}\,\rho^{i\omega c/a} \tag{3.254}$$

This has the same logarithmic phase singularity at the horizon ($\rho \to 0$) as the Hawking mode (3.210) at $r \to r_s$.

Analytic continuation. The Minkowski vacuum (the ether's ZPF ground state) is defined by positive-frequency modes with respect to the inertial time $T$. The Rindler mode (3.254) is positive-frequency with respect to Rindler time $\eta$. The two frequency definitions disagree, and the Bogoliubov transformation between them is determined by the analytic continuation of (3.254) across the Rindler horizon, following the identical $i\epsilon$ prescription of Section 3.13.4.

The result — which follows from the identical mathematical structure, with $\kappa \to a$ and $r - r_s \to \rho$ — is:

$$\frac{|\beta_\omega|^2}{|\alpha_\omega|^2} = e^{-2\pi\omega c/a} \tag{3.255}$$

The particle number measured by the accelerating observer:

$$\langle n_\omega \rangle = \frac{1}{e^{2\pi\omega c/a} - 1} \tag{3.256}$$

This is the Bose–Einstein distribution with temperature:

$$k_B T_U = \frac{\hbar\,a}{2\pi\,c} \tag{3.257}$$

3.15.4 The Unruh Temperature

Consistency with the Hawking effect. The Hawking temperature (Theorem 3.7) is $T_H = \hbar\kappa/(2\pi k_B c)$ with $\kappa = c^4/(4GM)$. The Unruh temperature is $T_U = \hbar a/(2\pi k_B c)$. The two are related by $T_U = T_H$ when $a = \kappa$ — that is, when the observer's proper acceleration equals the surface gravity of the black hole. This is the content of the equivalence principle: a static observer hovering at the black hole horizon (with proper acceleration $\kappa$) detects Hawking radiation at temperature $T_H = \hbar\kappa/(2\pi k_Bc)$, which is precisely the Unruh temperature for that acceleration. In the ether picture, the mechanism is the same in both cases — the observer's acceleration relative to the ZPF produces a thermal spectrum through mode-splitting at the causal horizon.

Numerical evaluation. To detect a temperature of $1$ K requires acceleration:

$$a = \frac{2\pi\,k_B\,c}{\hbar}\cdot T = \frac{2\pi \times 1.381 \times 10^{-23} \times 3 \times 10^8}{1.055 \times 10^{-34}}$$ $$= 2.47 \times 10^{20}\;\text{m/s}^2 \tag{3.259}$$

This is approximately $2.5 \times 10^{19}\,g$, far beyond any achievable laboratory acceleration. For the maximum sustained acceleration in current experiments (particle accelerators, $\sim 10^{25}$ m/s$^2$ for electrons in a strong laser field): $T_U \sim 40$ K — measurable in principle but overwhelmed by conventional thermal backgrounds.

3.15.5 Physical Interpretation in the Ether Framework

The Unruh effect has a transparent physical interpretation in the ether picture.

The inertial observer is at rest in the ether (or in uniform motion, which is equivalent by Lorentz invariance). The ether's ZPF (Section 6.1) is the ground state — the state of minimum energy with no real particles. The observer detects no radiation.

The accelerating observer moves with continuously changing velocity relative to the ether. The ZPF modes are Doppler-shifted by this motion: modes approaching from ahead are blueshifted, modes from behind are redshifted. For uniform acceleration, the cumulative Doppler shift over the observer's history produces a precise Planckian distribution — not because the ZPF has been heated, but because the accelerating observer's definition of "positive frequency" (and hence "particle") differs from the inertial definition.

The Rindler horizon is the surface behind the accelerating observer from which ZPF modes cannot catch up — they are redshifted to zero frequency before they arrive. This creates the same mode-splitting as the ether's gravitational horizon in Section 3.13: ZPF fluctuations near the Rindler horizon are torn apart by the observer's acceleration, with one component escaping forward (detected as a thermal particle) and the other falling behind the horizon.

The energy source. The detected radiation carries energy. Where does it come from? The observer must be accelerated by an external force (a rocket, an electric field). The work done by this force against the radiation reaction of the emitted particles provides the energy. In the ether picture: the accelerating observer disturbs the ZPF, converting virtual fluctuations into real excitations, with the energy supplied by whatever agent maintains the acceleration.

3.16 Ether Conservation and the Sink Problem Resolved

Section 3.9.2 noted an apparent paradox: the Schwarzschild ether flow (Eq 3.22) has nonzero divergence (Eq 3.57), so the Newtonian continuity equation $\partial_t\rho_e + \nabla\cdot(\rho_e\mathbf{V}) = 0$ is not satisfied for constant ether density $\rho_e = \rho_0$. Three possible resolutions were discussed — mass as ether sink, compressible ether, and effective description.

Theorem 3.5 resolves this problem. The ether's field equation is the Einstein equation, and the Einstein equation implies covariant energy-momentum conservation: $\nabla_\mu T^{\mu\nu} = 0$ (the contracted Bianchi identity). The ether IS conserved — but the relevant conservation law is the covariant one, not the Newtonian one. The Newtonian continuity equation is the leading-order approximation to $\nabla_\mu T^{\mu 0} = 0$; at post-Newtonian order, additional terms enter that account for the gravitational field's energy and momentum. This section derives these terms explicitly.

3.16.1 The Covariant Conservation Law

The stress-energy tensor of a perfect fluid with four-velocity $U^\mu$, energy density $\varepsilon$, and pressure $P$ is:

$$T^{\mu\nu} = (\varepsilon + P/c^2)\,U^\mu U^\nu + P\,g^{\mu\nu} \tag{3.260}$$

For the pressureless ($P = 0$) ether condensate in its normal phase:

$$T^{\mu\nu} = \varepsilon\,U^\mu U^\nu \tag{3.261}$$

The covariant conservation law $\nabla_\mu T^{\mu\nu} = 0$ expands to:

$$\nabla_\mu(\varepsilon\,U^\mu U^\nu) = U^\nu\nabla_\mu(\varepsilon\,U^\mu) + \varepsilon\,U^\mu\nabla_\mu U^\nu = 0 \tag{3.262}$$

Projecting along $U^\nu$ (energy conservation) and orthogonal to $U^\nu$ (momentum conservation):

$$\nabla_\mu(\varepsilon\,U^\mu) = 0 \tag{3.263}$$ $$\varepsilon\,U^\mu\nabla_\mu U^\nu = 0 \qquad (\text{geodesic equation}) \tag{3.264}$$

(3.263) is the covariant continuity equation — the relativistic generalisation of $\partial_t\rho + \nabla\cdot(\rho\mathbf{V}) = 0$.

3.16.2 Covariant Continuity in PG Coordinates

In PG coordinates, the ether four-velocity for a freely falling fluid element is $U^\mu = n^\mu c = (1, -V^i)$ (from Eq 3.96, with the normalisation $U_\mu U^\mu = -c^2$, and using the fact that freely falling ether elements have four-velocity proportional to the unit normal). More precisely, $U^0 = c\,n^0 = 1$ and $U^i = c\,n^i = -V^i$.

Verification of normalisation. $U_\mu U^\mu = g_{\mu\nu}U^\mu U^\nu$. With $U^0 = 1$, $U^i = -V^i$:

$$U_\mu U^\mu = g_{00}(U^0)^2 + 2g_{0i}U^0 U^i + g_{ij}U^i U^j$$ $$= -(c^2 - V^2) + 2V_i(-V^i) + \delta_{ij}V^iV^j$$ $$= -c^2 + V^2 - 2V^2 + V^2 = -c^2 \tag{3.265}$$

The covariant divergence of $\varepsilon\,U^\mu$ is:

$$\nabla_\mu(\varepsilon\,U^\mu) = \frac{1}{\sqrt{-g}}\,\partial_\mu\!\left(\sqrt{-g}\;\varepsilon\,U^\mu\right) \tag{3.266}$$

From Section 3.11.1: $\sqrt{-g} = c$ (constant, Eq 3.93). Therefore:

$$\nabla_\mu(\varepsilon\,U^\mu) = \frac{1}{c}\,\partial_\mu\!\left(c\,\varepsilon\,U^\mu\right)$$ $$= \partial_\mu(\varepsilon\,U^\mu) \tag{3.267}$$

The constant determinant eliminates the Christoffel-symbol terms entirely. Expanding:

$$\partial_T(\varepsilon\,U^0) + \partial_i(\varepsilon\,U^i) = \partial_T\varepsilon + \partial_i(\varepsilon\,(-V^i))$$ $$= \partial_T\varepsilon - \partial_i(\varepsilon\,V^i) = 0 \tag{3.268}$$

This is:

$$\frac{\partial\varepsilon}{\partial T} - \nabla\cdot(\varepsilon\,\mathbf{V}) = 0 \tag{3.269}$$

Note the minus sign in front of the divergence. This arises because the ether flows inward ($\mathbf{V}$ points toward the mass) while $U^i = -V^i$ points outward in the coordinate sense. The sign is consistent: the ether flows inward, so the flux $\varepsilon\mathbf{V}$ points inward, and $-\nabla\cdot(\varepsilon\mathbf{V})$ is positive near the mass, compensating the accumulation of ether.

Comparison with the Newtonian continuity equation. The standard form $\partial_t\rho + \nabla\cdot(\rho\mathbf{v}) = 0$ has a plus sign because $\mathbf{v}$ is the fluid velocity in the direction of motion. In PG coordinates, $\mathbf{V}$ is defined as the shift vector (pointing inward, $V^i > 0$ for inward radial flow), and the four-velocity has $U^i = -V^i$. With the identification $\varepsilon = \rho_0 c^2$ (constant energy density):

$$-\nabla\cdot(\rho_0 c^2\,\mathbf{V}) = 0 \implies \nabla\cdot\mathbf{V} = 0 \tag{3.270}$$

But from Eq (3.57): $\nabla\cdot\mathbf{V} \neq 0$ for the Schwarzschild flow. So either $\varepsilon$ is not constant, or (3.269) is not the complete equation. We now show it is the latter.

3.16.3 The Resolution: Relativistic Mass-Energy

The key is that $\varepsilon$ in the covariant conservation law (3.263) is the energy density measured by a comoving observer — the proper energy density. For a freely falling ether element, this includes the rest-mass energy AND the kinetic energy of the fall:

$$\varepsilon(r) = \rho_0 c^2 + \frac{1}{2}\rho_0 V^2 + O(V^4/c^2) \tag{3.271}$$

The kinetic energy density $\frac{1}{2}\rho_0 V^2$ is a post-Newtonian correction that modifies the conservation law.

However, a more precise treatment uses the ADM energy density. The Eulerian energy density is $\rho_E = T_{\mu\nu}\,n^\mu n^\nu$, which for a pressureless fluid with four-velocity $U^\mu$ is:

$$\rho_E = \varepsilon\,(U_\mu n^\mu)^2/c^2 \tag{3.272}$$

Computing $U_\mu n^\mu$: $U_\mu n^\mu = g_{\mu\nu}U^\nu n^\mu$. With $U^\nu = (1, -V^i)$ and $n^\mu = (1/c, -V^i/c)$:

$$U_\mu n^\mu = g_{00}\cdot 1 \cdot\frac{1}{c} + g_{0i}\!\left(1\cdot\frac{-V^i}{c} + (-V^i)\cdot\frac{1}{c}\right) + g_{ij}(-V^i)\frac{-V^j}{c}$$ $$= \frac{-(c^2-V^2)}{c} + V_i\frac{-2V^i}{c} + \frac{V^2}{c}$$ $$= \frac{-c^2+V^2}{c} - \frac{2V^2}{c} + \frac{V^2}{c} = -c \tag{3.273}$$

Therefore $\rho_E = \varepsilon\,c^2/c^2 = \varepsilon$. For a freely falling pressureless fluid in PG coordinates, the Eulerian and proper energy densities coincide.

The proper resolution comes from examining what (3.269) actually requires when applied to the exact Schwarzschild solution.

3.16.4 Exact Evaluation for Schwarzschild

For the Schwarzschild flow: $V^i = -\sqrt{2GM}\,x^i/r^{3/2}$ (inward, so $V^i$ as defined in Eq 3.22 has the radial component $V_r = -\sqrt{2GM/r}$ negative — but in the ADM convention used in Section 3.11 and throughout this section, $V^i$ is the shift vector with $V^r > 0$ for inward flow). We must be careful with signs.

In spherical coordinates with the convention $\mathbf{V} = V_r(r)\,\hat{\mathbf{r}}$ where $V_r > 0$ denotes inward flow (consistent with the PG metric cross-term $g_{Tr} = -V_r < 0$):

$$V_r = \sqrt{\frac{2GM}{r}} \tag{3.274}$$

The divergence is:

$$\nabla\cdot\mathbf{V} = \frac{1}{r^2}\frac{d}{dr}(r^2 V_r) = \frac{1}{r^2}\frac{d}{dr}\!\left(r^2\sqrt{\frac{2GM}{r}}\right) = \frac{1}{r^2}\frac{d}{dr}\!\left(\sqrt{2GM}\,r^{3/2}\right)$$ $$= \frac{\sqrt{2GM}}{r^2}\cdot\frac{3}{2}r^{1/2} = \frac{3}{2}\sqrt{\frac{2GM}{r^3}} \tag{3.275}$$

This is positive: the ether flow diverges (spreads out) as it accelerates inward through successively smaller spheres — the velocity increases faster ($\propto r^{-1/2}$) than the area decreases ($\propto r^{-2}$).

The Newtonian continuity equation with constant $\rho_0$ requires $\nabla\cdot\mathbf{V} = 0$. (3.275) shows $\nabla\cdot\mathbf{V} \neq 0$. This is the apparent paradox of Section 3.9.2.

The covariant conservation law (3.269) requires $\partial_T\varepsilon = \nabla\cdot(\varepsilon\mathbf{V})$ for a static configuration ($\partial_T\varepsilon = 0$). This gives $\nabla\cdot(\varepsilon\mathbf{V}) = 0$, which is satisfied if:

$$\varepsilon(r)\,V_r(r)\,r^2 = \text{const} \tag{3.276}$$

From (3.274): $V_r\,r^2 = \sqrt{2GM}\,r^{3/2}$. For (3.276) to hold:

$$\varepsilon(r) = \frac{C}{\sqrt{2GM}\,r^{3/2}} \tag{3.277}$$

for some constant $C$. The energy density must decrease as $r^{-3/2}$ — the ether is compressed near the mass.

Physical interpretation. The constant-density assumption ($\varepsilon = \rho_0 c^2$) that underlies Theorem 3.2 is an approximation. The exact solution of the covariant conservation law has $\varepsilon \propto r^{-3/2}$. This density variation modifies the conformal factor in the ether metric (Eq 3.17) from a constant to $\rho_e(r)/c$, introducing a position-dependent correction.

3.16.5 The Magnitude of the Correction

The density variation (3.277) is a post-Newtonian effect. To see this, normalise $\varepsilon$ so that $\varepsilon \to \rho_0 c^2$ as $r \to \infty$. The conservation law $\varepsilon V_r r^2 = \text{const}$ requires:

$$\varepsilon(r) = \rho_0 c^2 \cdot \frac{V_r(\infty)\,r_\infty^2}{V_r(r)\,r^2} \tag{3.278}$$

As $r \to \infty$: $V_r \to 0$ and the product $V_r\,r^2 \to \sqrt{2GM}\,r^{3/2} \to \infty$. The flux at infinity is not finite for the Schwarzschild profile — reflecting the fact that the Schwarzschild solution describes an eternal black hole with ether flowing in from infinity, not a realistic astrophysical configuration.

For a more physical setup, consider the ether flowing through a shell between radii $r_1$ and $r_2 > r_1$. The mass flux through a sphere of radius $r$ is:

$$\dot{M}_e(r) = 4\pi r^2\,\rho_e(r)\,V_r(r) \tag{3.279}$$

Conservation of mass flux ($\dot{M}_e$ independent of $r$) gives:

$$\rho_e(r) = \frac{\dot{M}_e}{4\pi r^2 V_r(r)} = \frac{\dot{M}_e}{4\pi\sqrt{2GM}\,r^{3/2}} \tag{3.280}$$

The fractional density deviation from $\rho_0$ at radius $r$ is:

$$\frac{\delta\rho_e}{\rho_0} = \frac{\rho_e(r) - \rho_0}{\rho_0} \tag{3.281}$$

For the regime where the density variation is small ($|\delta\rho_e/\rho_0| \ll 1$), the mass flux is approximately $\dot{M}_e \approx 4\pi r^2\rho_0 V_r$, and (3.280) gives $\rho_e \approx \rho_0$ (self-consistent). The correction appears at the next order.

The post-Newtonian expansion. Write $\varepsilon = \rho_0 c^2(1 + \delta)$ with $|\delta| \ll 1$. The covariant conservation (3.269) in steady state:

$$\nabla\cdot(\varepsilon\mathbf{V}) = \varepsilon\,\nabla\cdot\mathbf{V} + \mathbf{V}\cdot\nabla\varepsilon = 0 \tag{3.282}$$ $$\rho_0 c^2(1+\delta)\,\nabla\cdot\mathbf{V} + \rho_0 c^2\,V_r\,\frac{d\delta}{dr} = 0 \tag{3.283}$$

At leading order ($\delta \ll 1$):

$$\nabla\cdot\mathbf{V} + V_r\frac{d\delta}{dr} = 0 \tag{3.284}$$

Using (3.275): $\nabla\cdot\mathbf{V} = 3\sqrt{2GM}/(2r^{3/2})$. The velocity gradient term: $V_r\,d\delta/dr = \sqrt{2GM/r}\,\delta'(r)$. Therefore:

$$\sqrt{\frac{2GM}{r}}\,\delta'(r) = -\frac{3\sqrt{2GM}}{2r^{3/2}} \tag{3.285}$$ $$\delta'(r) = -\frac{3}{2r} \tag{3.286}$$ $$\delta(r) = -\frac{3}{2}\ln\!\left(\frac{r}{r_\infty}\right) \tag{3.287}$$

where $r_\infty$ is the reference radius at which $\delta = 0$ (i.e., $\varepsilon = \rho_0 c^2$).

This logarithmic divergence signals the breakdown of the small-$\delta$ approximation at $r \sim r_\infty e^{-2/3}$. But for the physically relevant regime $r \gg r_s$:

$$|\delta(r)| = \frac{3}{2}\left|\ln\!\left(\frac{r}{r_\infty}\right)\right| \tag{3.288}$$

The density correction is logarithmic, not polynomial, in $r$ — it grows slowly and remains small over many decades of radius.

Effect on the metric. The conformal factor in the ether metric (Eq 3.17) becomes $\rho_e(r)/c = \rho_0(1+\delta)/c$. This introduces a position-dependent conformal rescaling:

$$g_{\mu\nu}^{\text{corrected}} = (1 + \delta(r))\,g_{\mu\nu}^{\text{PG}} \tag{3.289}$$

For null geodesics (light paths, gravitational redshift, horizon locations), a static conformal factor does not alter the predictions — null geodesics are conformally invariant (Section 3.3, Remark on the conformal factor). For timelike geodesics (massive particle orbits), the correction introduces a fractional deviation:

$$\frac{\delta a}{a_N} \sim \delta \sim \frac{3}{2}\ln(r/r_\infty) \tag{3.290}$$

This is parametrically of order $\ln(r/r_s)$ — a logarithmic post-Newtonian correction.

Relation to the field equation. The density variation (3.287) arises from exact covariant conservation ($\nabla_\mu T^{\mu\nu} = 0$, guaranteed by the Bianchi identity and Theorem 3.5), not from an arbitrary gauge transformation. The constant-density PG metric solves the Einstein equation exactly; the density correction (3.289) is the leading post-Newtonian correction that maintains covariant conservation, and is therefore consistent with the Einstein equation at all orders.

3.16.6 Resolution of the Three Interpretations

We now revisit the three interpretations offered in Section 3.9.2.

(a) Mass as ether sink. In the covariant formulation, the ether IS conserved: $\nabla_\mu(\varepsilon U^\mu) = 0$ (Eq 3.263). There is no sink. The apparent non-conservation arises from using the Newtonian continuity equation (which ignores relativistic corrections) with a constant density (which ignores the compressibility required by the exact conservation law). The "sink" interpretation is an artifact of an inconsistent approximation.

(b) Compressible ether. Correct. The exact conservation law requires $\varepsilon(r) \propto 1/(r^2 V_r) \propto r^{-3/2}$. The ether is compressible, and the density variation is a post-Newtonian effect. This does not conflict with Theorem 3.2: the constant-density assumption yields the Schwarzschild metric up to a constant conformal factor, and the density correction (3.289) is a higher-order conformal correction that does not affect the leading-order predictions.

(c) Effective description. Also correct, in the sense that the constant-density PG metric is the leading-order effective description. The density correction (3.287) is the first post-Newtonian correction to this effective description. For all weak-field applications (solar system tests, galaxy dynamics, gravitational waves), the correction is negligible.

The definitive resolution: The ether is conserved in the relativistic sense ($\nabla_\mu T^{\mu\nu} = 0$, guaranteed by the Bianchi identity and Theorem 3.5). The constant-density assumption is an approximation valid at leading post-Newtonian order, with logarithmic corrections at the next order. No ether is created, destroyed, or absorbed. The "sink problem" of Section 3.9.2 is resolved.

3.16.7 Comparison with Standard GR

In standard GR, the analogous question does not arise because there is no "medium" whose density could vary. The stress-energy tensor of matter satisfies $\nabla_\mu T^{\mu\nu}_{\text{matter}} = 0$, which determines the matter's motion (geodesic equation for dust). The gravitational field carries energy described by the Landau–Lifshitz pseudotensor, but this energy is not localisable — its integral over a region depends on the choice of coordinates [182].

In the ether framework, the gravitational field energy is localisable: it is the kinetic energy $\frac{1}{2}\rho_e V^2$ of the ether flow (Section 3.11.5). The total energy — ether rest mass plus kinetic energy — is exactly conserved by the covariant conservation law. This is an interpretive advantage of the ether framework: the energy of the gravitational field has a concrete, localisable meaning in the ether rest frame as the kinetic energy of a physical medium, unlike the non-localisable pseudotensor energy of standard GR.

3.17 The Einstein Equation from Ether Thermodynamics

Section 3.11 derived the Einstein equation from the ether metric via the Weinberg–Deser–Lovelock uniqueness theorems (Theorem 3.5). The derivation takes the Newtonian limit as its single empirical input and uses uniqueness to fix the nonlinear extension. We now present an independent second derivation that proceeds from the thermodynamic properties of the ether's zero-point field, following the route pioneered by Jacobson [183]. The two derivations use different physical assumptions and arrive at the same result — the Einstein equation — providing strong mutual corroboration.

The Jacobson route does not use the Newtonian limit as an input. Its inputs are: (i) the ether supports a zero-point field whose entanglement entropy across any horizon is proportional to the horizon area (derived from mode counting in Section 3.17.2); (ii) the Unruh effect holds for the acoustic metric (established in Theorem 3.9, Section 3.15); (iii) the ether is in local thermodynamic equilibrium (a natural assumption for a superfluid at scales much larger than the healing length). The Einstein equation — including the Newtonian limit — is an output.

3.17.1 The Physical Ingredients

The derivation requires three results already established in this monograph:

(i) The acoustic metric (Theorem 3.1, Section 3.1). The ether's perturbations propagate on the effective Lorentzian metric $g_{\mu\nu}$ determined by the ether's flow velocity and sound speed. This metric defines causal structure, horizons, and geodesics.

(ii) The Unruh effect (Theorem 3.9, Section 3.15). An observer accelerating uniformly at proper acceleration $a$ through the ether detects thermal radiation at temperature $T_U = \hbar a/(2\pi k_B c)$. Equivalently, any local Rindler horizon in the ether metric has a temperature:

$$T = \frac{\hbar\kappa}{2\pi k_B c} \tag{3.291}$$

where $\kappa$ is the surface gravity of the horizon.

(iii) The ether's zero-point field (Section 6.1, Theorem 4.2). The ether supports electromagnetic and phononic quantum fluctuations whose spectral density is uniquely determined by Lorentz invariance.

To these we add one new ingredient:

(iv) The entanglement entropy of the ZPF across a horizon is proportional to the horizon area. We derive this in Section 3.17.2.

3.17.2 Entanglement Entropy from ZPF Mode Counting

Consider a Rindler horizon $\mathcal{H}$ — a null surface in the ether metric that divides space into two causally disconnected regions. The ether's ZPF modes span both regions. Tracing over the modes in the causally inaccessible region produces an entanglement entropy for the accessible region.

The mode-counting argument. A free scalar field $\phi$ on a $(3+1)$-dimensional spacetime, restricted to modes with wavelength $\lambda > \ell_{\text{UV}}$ (the UV cutoff), has entanglement entropy across a smooth surface of area $A$ given by [184, 185]:

$$S_{\text{ent}} = \frac{\mathcal{N}\,k_B\,A}{48\pi\,\ell_{\text{UV}}^2} + \text{sub-leading terms} \tag{3.292}$$

where $\mathcal{N}$ is the number of field species contributing to the entropy. This result was derived by Bombelli, Koul, Lee, and Sorkin [184] and independently by Srednicki [185] using the replica trick and lattice regularisation. The area scaling is universal — it depends only on the dimensionality of spacetime and the UV cutoff, not on the details of the field theory.

Derivation sketch. The entanglement entropy is computed by:

(a) Discretising the field $\phi$ on a spatial lattice with spacing $\ell_{\text{UV}}$.

(b) Computing the ground-state density matrix $\hat{\rho}$ of the full system.

(c) Tracing over lattice sites on one side of the surface: $\hat{\rho}_{\text{red}} = \text{Tr}_{\text{outside}}\hat{\rho}$.

(d) Computing $S = -k_B\,\text{Tr}(\hat{\rho}_{\text{red}}\ln\hat{\rho}_{\text{red}})$.

For a free field, $\hat{\rho}$ is Gaussian, and the entropy is determined by the two-point correlator restricted to the surface. The leading term is proportional to the number of lattice sites on the surface, which scales as $A/\ell_{\text{UV}}^2$. The coefficient $1/(48\pi)$ was computed by Srednicki [185] for a minimally coupled scalar field via numerical lattice diagonalisation and confirmed analytically by Solodukhin [186].

Application to the ether. In the ether framework, the UV cutoff is the transverse microstructure scale $\ell_e$ (for electromagnetic modes) or the healing length $\xi$ (for phononic modes). The number of contributing species $\mathcal{N}$ counts the electromagnetic polarisations (2) and the phonon mode (1), giving $\mathcal{N} = 3$ at minimum. Additional species (from multi-component ether structure required by Proposition 6.1) increase $\mathcal{N}$.

Writing the entanglement entropy as:

$$S_{\text{ent}} = \eta\,A \tag{3.293}$$

where $\eta$ is the entropy per unit area:

$$\eta = \frac{\mathcal{N}\,k_B}{48\pi\,\ell_e^2} \tag{3.294}$$

For later convenience, we define $\eta$ in terms of a length scale $\ell_{\text{grav}}$ by:

$$\eta = \frac{k_B}{4\,\ell_{\text{grav}}^2} \tag{3.295}$$

so that $S = k_B A/(4\ell_{\text{grav}}^2)$, matching the Bekenstein–Hawking form. The gravitational length scale is:

$$\ell_{\text{grav}}^2 = \frac{12\pi\,\ell_e^2}{\mathcal{N}} \tag{3.296}$$

We will show that $\ell_{\text{grav}}$ is identified with the Planck length $\ell_P = \sqrt{\hbar G/c^3}$, from which Newton's constant $G$ is determined.

3.17.3 Local Rindler Horizons in the Ether

At any point $p$ in the ether spacetime, we can construct a local Rindler horizon by considering an observer accelerating with proper acceleration $a$ through the ether. In the observer's frame, the ether flows past with increasing velocity, and a causal horizon exists at a proper distance $c^2/a$ behind the observer (Section 3.15.1).

The local construction. Choose a point $p$ and a future-directed null vector $k^\mu$ at $p$. The null geodesic through $p$ with tangent $k^\mu$ generates a small patch of a null surface $\mathcal{H}$. We parametrise this surface by affine parameter $\lambda$, with $\lambda = 0$ at $p$ and $k^\mu = dx^\mu/d\lambda$.

The approximate boost Killing vector in the neighbourhood of $p$ is:

$$\chi^\mu = \frac{\kappa}{c^2}\,\lambda\,k^\mu \tag{3.297}$$

where $\kappa$ is the surface gravity (with dimensions of acceleration, as in Eqs. 3.195 and 3.253c) and $\lambda$ is the affine parameter (with dimensions of length). The factor $c^2$ converts from the physical surface gravity $\kappa$ [m/s$^2$] to the geometric surface gravity $\kappa/c^2$ [m$^{-1}$], ensuring $\chi^\mu$ is dimensionless. This satisfies $\chi^\mu|_{\lambda=0} = 0$ (the Killing vector vanishes at the bifurcation point) and $\nabla_\mu\chi_\nu|_{\lambda=0} = (\kappa/c^2)\,(k_\mu n_\nu - k_\nu n_\mu)$ where $n^\mu$ is an auxiliary null vector normalised as $k_\mu n^\mu = -1$.

In the ether framework, the local Rindler horizon is a local acoustic horizon: the surface where the ether flow velocity relative to the accelerating observer equals the local sound speed (or light speed, for electromagnetic perturbations). The Unruh temperature associated with this horizon is ((3.291)):

$$T = \frac{\hbar\kappa}{2\pi k_B c} \tag{3.298}$$

3.17.4 The Energy Flux Through the Horizon

The energy flux through the horizon patch is computed by contracting the matter stress-energy tensor with the approximate boost Killing vector and integrating over the horizon surface.

The surface element. The null surface $\mathcal{H}$ has surface element $d\Sigma^\mu = k^\mu\,dA\,d\lambda$, where $dA$ is the cross-sectional area element of the pencil of generators at affine parameter $\lambda$.

The energy flux. The heat absorbed by the horizon (the energy flowing through $\mathcal{H}$ as seen by the accelerating observer) is:

$$\delta Q = -\int_{\mathcal{H}} T_{\mu\nu}\,\chi^\mu\,d\Sigma^\nu \tag{3.299}$$

The minus sign follows the thermodynamic convention: positive $\delta Q$ corresponds to energy entering the system (the horizon).

Substituting $\chi^\mu = (\kappa/c^2)\lambda k^\mu$ and $d\Sigma^\nu = k^\nu dA\,d\lambda$:

$$\delta Q = -\int T_{\mu\nu}\,\frac{\kappa}{c^2}\lambda k^\mu\,k^\nu\,dA\,d\lambda$$ $$= -\frac{\kappa}{c^2}\int T_{\mu\nu}\,k^\mu k^\nu\,\lambda\,dA\,d\lambda \tag{3.300}$$

For matter satisfying the null energy condition ($T_{\mu\nu}k^\mu k^\nu \geq 0$), the integrand is non-negative. With $\lambda > 0$ (past the bifurcation point, in the future direction) and $\kappa > 0$, the integral is positive. The minus sign ensures $\delta Q > 0$: energy flowing through the horizon heats the system.

3.17.5 The Area Change from the Raychaudhuri Equation

The cross-sectional area of the pencil of null generators evolves according to the Raychaudhuri equation for a null geodesic congruence [165, 187]:

$$\frac{d\theta}{d\lambda} = -\frac{1}{2}\theta^2 - \sigma_{\mu\nu}\sigma^{\mu\nu} - R_{\mu\nu}\,k^\mu k^\nu \tag{3.301}$$

where $\theta = \nabla_\mu k^\mu$ is the expansion scalar and $\sigma_{\mu\nu}$ is the shear tensor of the congruence.

Near-bifurcation expansion. At the bifurcation point $\lambda = 0$, the generators of the local Rindler horizon have:

$$\theta(\lambda=0) = 0, \qquad \sigma_{\mu\nu}(\lambda=0) = 0 \tag{3.302}$$

(The expansion and shear vanish at the bifurcation point of any Killing horizon, by the properties of the boost Killing vector [165].)

To leading order in $\lambda$, (3.301) gives:

$$\frac{d\theta}{d\lambda}\bigg|_{\lambda=0} = -R_{\mu\nu}\,k^\mu k^\nu \tag{3.303}$$

Integrating:

$$\theta(\lambda) = -R_{\mu\nu}\,k^\mu k^\nu\,\lambda + O(\lambda^2) \tag{3.304}$$

The area change. The fractional rate of change of the cross-sectional area is $d(\ln A)/d\lambda = \theta$. For the pencil of generators from $\lambda = 0$ to some small $\lambda_0$:

$$\delta A = \int_0^{\lambda_0}\!\int \theta\,dA\,d\lambda = -\int R_{\mu\nu}\,k^\mu k^\nu\,\lambda\,dA\,d\lambda + O(\lambda_0^2) \tag{3.305}$$

The entropy change associated with this area change is (using (3.293)):

$$\delta S = \eta\,\delta A = -\eta\int R_{\mu\nu}\,k^\mu k^\nu\,\lambda\,dA\,d\lambda \tag{3.306}$$

3.17.6 The Clausius Relation and the Field Equation

The Clausius relation. For a system in local thermodynamic equilibrium at temperature $T$, the entropy change produced by a reversible heat flux $\delta Q$ is:

$$\delta Q = T\,\delta S \tag{3.307}$$

In the ether framework, this relation is the local equilibrium condition for the ether at each spacetime point. The temperature is the Unruh temperature (3.298), the entropy is the ZPF entanglement entropy (3.293), and the heat is the energy flux through the local horizon (3.300).

Substitution. Substituting (3.300), (3.298), and (3.306) into (3.307):

$$-\frac{\kappa}{c^2}\int T_{\mu\nu}\,k^\mu k^\nu\,\lambda\,dA\,d\lambda = \frac{\hbar\kappa}{2\pi k_B c}\cdot\left(-\eta\int R_{\mu\nu}\,k^\mu k^\nu\,\lambda\,dA\,d\lambda\right) \tag{3.308}$$

Cancellation. The factors $\kappa$, $\lambda$, $dA$, and $d\lambda$ appear identically on both sides and cancel. The minus signs cancel. The factor $1/c^2$ on the left combines with $1/c$ in the Unruh temperature on the right:

$$\frac{1}{c^2}\,T_{\mu\nu}\,k^\mu k^\nu = \frac{\hbar\eta}{2\pi k_B c}\;R_{\mu\nu}\,k^\mu k^\nu \tag{3.309a}$$

Multiplying both sides by $c^2$:

$$T_{\mu\nu}\,k^\mu k^\nu = \frac{\hbar c\,\eta}{2\pi k_B}\;R_{\mu\nu}\,k^\mu k^\nu \tag{3.309}$$

This must hold for every null vector $k^\mu$ at every point, since the construction of the local Rindler horizon is arbitrary — any null direction can be chosen.

Define the coupling constant. Set:

$$\alpha \equiv \frac{\hbar c\,\eta}{2\pi k_B} \tag{3.310}$$

Substituting the entropy density $\eta = k_B/(4\ell_{\text{grav}}^2)$ from (3.295):

$$\alpha = \frac{\hbar c}{2\pi}\cdot\frac{1}{4\ell_{\text{grav}}^2} = \frac{\hbar c}{8\pi\,\ell_{\text{grav}}^2} \tag{3.311}$$

Then (3.309) becomes:

$$T_{\mu\nu}\,k^\mu k^\nu = \alpha\,R_{\mu\nu}\,k^\mu k^\nu \qquad \text{for all null } k^\mu \tag{3.312}$$

3.17.7 From the Null Condition to the Einstein Equation

Step 1: Tensor structure. Since (3.312) holds for all null vectors $k^\mu$, and $g_{\mu\nu}k^\mu k^\nu = 0$ for any null vector, the most general tensor equation consistent with (3.312) is:

$$T_{\mu\nu} - \alpha\,R_{\mu\nu} = f\,g_{\mu\nu} \tag{3.313}$$

for some scalar function $f$. This follows because any symmetric tensor $A_{\mu\nu}$ satisfying $A_{\mu\nu}k^\mu k^\nu = 0$ for all null $k^\mu$ must be proportional to $g_{\mu\nu}$. (Proof: in an orthonormal frame, the condition $A_{00} + A_{11} + 2A_{01} = 0$ for all null combinations implies $A_{ij} = A_{00}\delta_{ij}$ and $A_{0i} = 0$, hence $A_{\mu\nu} = A_{00}\,g_{\mu\nu}$.)

Step 2: Fix $f$ using conservation laws. Take the covariant divergence of both sides of (3.313):

$$\nabla^\mu T_{\mu\nu} - \alpha\,\nabla^\mu R_{\mu\nu} = \nabla_\nu f \tag{3.314}$$

Energy-momentum conservation: $\nabla^\mu T_{\mu\nu} = 0$ (for matter in the ether, this follows from the ether's Euler equation).

Contracted Bianchi identity: $\nabla^\mu R_{\mu\nu} = \frac{1}{2}\nabla_\nu R$ (an algebraic identity of Riemannian geometry, following from the symmetries of the Riemann tensor).

Substituting into (3.314):

$$0 - \frac{\alpha}{2}\nabla_\nu R = \nabla_\nu f \tag{3.315}$$

Integrating:

$$f = -\frac{\alpha}{2}R + \frac{\Lambda_0}{\alpha} \tag{3.316}$$

where $\Lambda_0/\alpha$ is an integration constant (written this way for later convenience).

Step 3: Assemble the field equation. Substituting (3.316) into (3.313):

$$T_{\mu\nu} = \alpha\,R_{\mu\nu} - \frac{\alpha}{2}R\,g_{\mu\nu} + \frac{\Lambda_0}{\alpha}\,g_{\mu\nu}$$ $$= \alpha\!\left(R_{\mu\nu} - \frac{1}{2}R\,g_{\mu\nu}\right) + \frac{\Lambda_0}{\alpha}\,g_{\mu\nu}$$ $$= \alpha\,G_{\mu\nu} + \frac{\Lambda_0}{\alpha}\,g_{\mu\nu} \tag{3.317}$$

Rearranging:

$$\boxed{G_{\mu\nu} + \Lambda\,g_{\mu\nu} = \frac{1}{\alpha}\,T_{\mu\nu}} \tag{3.318}$$

where $\Lambda \equiv \Lambda_0/\alpha^2$.

Step 4: Identify Newton's constant. Comparing with the standard form $G_{\mu\nu} + \Lambda g_{\mu\nu} = (8\pi G/c^4)T_{\mu\nu}$:

$$\frac{1}{\alpha} = \frac{8\pi G}{c^4} \tag{3.319}$$

Substituting $\alpha$ from (3.311):

$$\frac{8\pi\,\ell_{\text{grav}}^2}{\hbar c} = \frac{8\pi G}{c^4} \tag{3.320}$$ $$G = \frac{c^3\,\ell_{\text{grav}}^2}{\hbar} \tag{3.321}$$

Equivalently:

$$\ell_{\text{grav}} = \sqrt{\frac{\hbar G}{c^3}} = \ell_P \tag{3.322}$$

The gravitational length scale is the Planck length. Newton's constant is determined by the ZPF entanglement entropy through the entropy–area relation (3.295):

$$\eta = \frac{k_B}{4\ell_P^2} = \frac{k_B\,c^3}{4\hbar G} \tag{3.323}$$

This is the Bekenstein–Hawking entropy density. In the ether framework, it is derived from the ZPF mode counting (Section 3.17.2), not postulated.

3.17.8 The Thermodynamic Field Equation

3.17.9 What the Thermodynamic Derivation Achieves

Independence from the uniqueness derivation. Theorem 3.10 and Theorem 3.5 derive the same equation — the Einstein equation — from different physical assumptions:

	Theorem 3.5 (Uniqueness)	Theorem 3.10 (Thermodynamic)
Physical input	Newtonian limit + Lorentz invariance	ZPF entropy + Unruh temperature + local equilibrium
Mathematical tool	Weinberg–Deser–Lovelock uniqueness	Clausius relation + Raychaudhuri equation
What is derived	Nonlinear extension of linearised gravity	Full nonlinear equation from thermodynamics
Role of $G$	Calibrated from Newtonian limit	Determined by ZPF mode counting
Role of $\Lambda$	Allowed by uniqueness, not determined	Integration constant from conservation

The two derivations share no premises and use different mathematics. That they converge on the same result is strong evidence for the robustness of the ether gravitational programme.

The Newtonian limit as output, not input. In Theorem 3.5, the Newtonian limit $\nabla^2\Phi = 4\pi G\rho_m$ is the starting point from which the full nonlinear equation is derived via uniqueness. In Theorem 3.10, the Newtonian limit is a consequence: the weak-field, slow-motion limit of (3.324) yields $\nabla^2\Phi = 4\pi G\rho_m$ by the standard post-Newtonian expansion. The value of $G$ is not calibrated to Newtonian gravity — it is determined by the ether's UV cutoff through $G = c^3\ell_{\text{grav}}^2/\hbar$ ((3.321)).

Newton's constant from microphysics. (3.321) expresses $G$ in terms of $\ell_{\text{grav}}$, which is fixed by the ether's transverse microstructure scale $\ell_e$ and species count $\mathcal{N}$ through (3.296):

$$G = \frac{c^3}{\hbar}\cdot\frac{12\pi\,\ell_e^2}{\mathcal{N}} \tag{3.325}$$

Inverting: the observed value $G = 6.674 \times 10^{-11}$ m$^3$kg$^{-1}$s$^{-2}$ constrains the ether parameters:

$$\frac{\ell_e^2}{\mathcal{N}} = \frac{G\hbar}{12\pi c^3} = \frac{\ell_P^2}{12\pi} = 6.93 \times 10^{-72}\;\text{m}^2 \tag{3.326}$$

where $\ell_P = \sqrt{\hbar G/c^3} = 1.616 \times 10^{-35}$ m and $\ell_P^2/(12\pi) = 2.611 \times 10^{-70}/(37.70) = 6.93 \times 10^{-72}$ m$^2$.

For $\mathcal{N} = 3$ (two EM polarisations + one phonon mode): $\ell_e^2 = 3 \times 6.93 \times 10^{-72} = 2.08 \times 10^{-71}$ m$^2$, hence $\ell_e = \sqrt{2.08 \times 10^{-71}} = 4.56 \times 10^{-36}$ m $\approx 0.28\,\ell_P$. This is remarkable: the ether's transverse microstructure scale is of order the Planck length — the same scale at which Lorentz-violating dispersion corrections (Section 3.8) are expected. The two constraints are consistent.

Connection to the cosmological constant. The integration constant $\Lambda$ in (3.318) is not determined by the thermodynamic derivation. It is fixed by the ether's global thermodynamic state — specifically, by the phonon ZPF energy density (Section 4.3, Theorem 4.2). This separation is physically natural: the local dynamics ($G_{\mu\nu}$) are determined by local thermodynamics (entropy, temperature), while the global vacuum energy ($\Lambda$) is determined by the global ground state.

The logical chain. The complete derivation from ether properties to the Einstein equation is:

$$\underbrace{\text{Superfluid ether}}_{\text{Goldstone}} \;\to\; \underbrace{\text{Acoustic metric}}_{\text{Thm 3.1}} \;\to\; \underbrace{\text{ZPF on metric}}_{\text{§6.1}} \;\to\; \underbrace{S = \eta A}_{\text{Mode counting}}$$ $$\to\; \underbrace{T = \frac{\hbar\kappa}{2\pi k_Bc}}_{\text{Thm 3.9}} \;\to\; \underbrace{\delta Q = T\,dS}_{\text{Clausius}} \;\to\; \underbrace{G_{\mu\nu} + \Lambda g_{\mu\nu} = \frac{8\pi G}{c^4}T_{\mu\nu}}_{\text{Thm 3.10}}$$

None of the inputs on the left is the Einstein equation or the Newtonian limit in disguise. The gravitational field equation emerges from the thermodynamic properties of a quantum medium — the ether's zero-point fluctuations and their entanglement structure.

3.18 Summary of Part II Results

We collect the key results of Section 3:

1. Acoustic metric derived (Theorem 3.1): Sound in a moving fluid propagates on an effective curved spacetime. Section 3.1.

2. Painlevé–Gullstrand identity (Theorem 3.2): Schwarzschild gravity is exactly the acoustic metric for constant-density ether flowing inward at free-fall velocity. Section 3.3.

3. All Schwarzschild predictions reproduced: Redshift, light bending, Shapiro delay, perihelion precession, GW speed. Section 3.5.

4. Emergent Lorentz invariance (Theorem 3.3): Exact at $\lambda \gg \ell_e$, violated at $O((\ell_e/\lambda)^2)$. Section 3.8.

5. Ether field equation (Eq 3.56): Weak-field equation determines ether flow from matter distribution. Section 3.9.

6. Kerr–Ether Identity (Theorem 3.4): Kerr metric in Doran coordinates is the ether metric with spiralling (gravitoelectric + gravitomagnetic) flow. GP-B precession rates reproduced. Section 3.10.

7. Einstein equation derived (Theorems 3.5, 3.10): The ether's complete nonlinear field equation is the Einstein equation, derived by two independent routes: the Weinberg–Deser–Lovelock uniqueness theorems (Section 3.11) and the thermodynamics of the ether's zero-point field via the Jacobson–Clausius argument (Section 3.17). Newton's constant $G$ is determined by ZPF entanglement entropy. This resolves the central open problem of the gravitational programme.

8. PPN parameters (Theorem 3.6): $\beta = \gamma = 1$; all ten PPN parameters match GR exactly. Consistent with Cassini, lunar laser ranging, and all solar system tests. Section 3.12.

9. Hawking radiation (Theorem 3.7): Thermal radiation at $T_H = \hbar\kappa/(2\pi k_Bc)$ from the ether horizon, with Kerr extension and trans-Planckian resolution. Section 3.13.

10. GW polarisations (Theorem 3.8): Exactly two tensor modes ($+$ and $\times$); scalar breathing mode is non-radiative. Consistent with LIGO-Virgo-KAGRA observations. Section 3.14.

11. Unruh effect (Theorem 3.9): Uniformly accelerating observer detects thermal radiation at $T_U = \hbar a/(2\pi k_Bc)$. Section 3.15.

12. Ether conservation resolved (Section 3.16): Covariant conservation $\nabla_\mu T^{\mu\nu} = 0$ is satisfied exactly. The constant-density assumption is a leading-order approximation with logarithmic post-Newtonian corrections.

13. Einstein equation from ether thermodynamics (Theorem 3.10): An independent second derivation of the Einstein equation from the ZPF entanglement entropy, the Unruh temperature, and the Clausius relation. Newton's constant $G = c^3\ell_P^2/\hbar$ is determined by ZPF mode counting. Section 3.17.

Open problems remaining in Section 3: (a) Singularity resolution (Section 3.6) — speculative UV modification, not derived. (b) Numerical implementation of binary mergers using standard numerical relativity on the ether metric. (c) The ether's transverse microstructure scale $\ell_e$ — undetermined.

The ether framework now reproduces the complete content of general relativity: kinematic structure (Theorems 3.2, 3.4), dynamical field equations (Theorems 3.5 and 3.10, by two independent routes), post-Newtonian tests (Theorem 3.6), quantum radiation effects (Theorems 3.7, 3.9), and gravitational wave phenomenology (Theorem 3.8). The gravitational programme that was limited to weak-field results in Section 3.9 is now complete at the level of the classical field equations.