7. Beyond Ground States: The Schrödinger Equation from Ether Dynamics

7.1 The Problem

Section 5 proved that the ether's zero-point field maintains atomic ground states. This is insufficient. Without the Schrödinger equation, there are no excited states, no discrete spectra, no transition rates — no quantum mechanics. The challenge is to derive the Schrödinger equation from ether physics, not merely to show consistency with specific stationary states. As a consequence of the Schrödinger derivation, we also obtain the de Broglie–Bohm guidance equation (Section 7.6) — establishing that the ether framework encompasses Bohmian mechanics for single particles.

The gap is filled by a result due to Nelson [19]: a particle undergoing Brownian motion with diffusion coefficient $D = \hbar/(2m)$ in an external potential satisfies the Schrödinger equation. Nelson's derivation is proven mathematics — there is no conjecture in it. What Nelson could not do was explain why $D = \hbar/(2m)$. In his framework, this was a postulate.

The ether framework resolves this. We show that Boyer's proven result (Section 6.2) — the ZPF maintains the quantum ground state — determines $D = \hbar/(2m)$ uniquely, through Nelson's stochastic dynamics. The diffusion coefficient is not a free parameter; it is fixed by the ether's zero-point fluctuation spectrum.

The logical chain is:

$$\text{Ether ZPF spectrum (Thm 4.2)} \xrightarrow{\text{Boyer (Thm 6.1)}} \rho_0 = |\psi_0|^2 \xrightarrow{\text{Nelson dynamics}} D = \frac{\hbar}{2m} \xrightarrow{\text{Nelson (Thm 7.1)}} i\hbar\dot{\psi} = \hat{H}\psi$$

Every arrow is a derivation, not an assumption. We present each in full.

7.2 Nelson's Stochastic Mechanics

7.2.1 Kinematics

A particle of mass $m$ moves through the ether. The ZPF imparts random impulses, causing a diffusive motion described by a Markov process $\mathbf{x}(t)$ with diffusion coefficient $D$ (to be determined). Because Brownian trajectories are continuous but non-differentiable, the velocity $d\mathbf{x}/dt$ does not exist. Nelson [19] defined two mean derivatives:

$$D_+\mathbf{x}(t) = \lim_{\Delta t\to 0^+}\left\langle\frac{\mathbf{x}(t+\Delta t)-\mathbf{x}(t)}{\Delta t}\bigg|\mathbf{x}(t)\right\rangle \equiv \mathbf{b}_+ \tag{7.1}$$ $$D_-\mathbf{x}(t) = \lim_{\Delta t\to 0^+}\left\langle\frac{\mathbf{x}(t)-\mathbf{x}(t-\Delta t)}{\Delta t}\bigg|\mathbf{x}(t)\right\rangle \equiv \mathbf{b}_- \tag{7.2}$$

These are the forward and backward drifts, conditional on the current position. From them, two physical velocities are constructed:

$$\mathbf{v} = \frac{\mathbf{b}_+ + \mathbf{b}_-}{2}\qquad\text{(current velocity: transports probability)} \tag{7.3}$$ $$\mathbf{u} = \frac{\mathbf{b}_+ - \mathbf{b}_-}{2}\qquad\text{(osmotic velocity: driven by diffusion)} \tag{7.4}$$

The stochastic differential equation governing the process is:

$$d\mathbf{x} = \mathbf{b}_+\,dt + \sqrt{2D}\,d\mathbf{W}(t) \tag{7.5}$$

where $d\mathbf{W}$ is a Wiener process increment ($\langle dW_i\,dW_j\rangle = \delta_{ij}\,dt$).

7.2.2 The Continuity and Osmotic Equations

The probability density $\rho(\mathbf{x},t)$ satisfies two Fokker–Planck equations, one for each time direction:

$$\partial_t\rho = -\nabla\!\cdot\!(\mathbf{b}_+\rho) + D\nabla^2\rho \qquad\text{(forward)} \tag{7.6}$$ $$\partial_t\rho = -\nabla\!\cdot\!(\mathbf{b}_-\rho) - D\nabla^2\rho \qquad\text{(backward)} \tag{7.7}$$

Proof of (7.6). The forward Fokker–Planck equation for the Itô process (7.5) is standard [112]: the probability flux is $\mathbf{J} = \mathbf{b}_+\rho - D\nabla\rho$, so $\partial_t\rho = -\nabla\cdot\mathbf{J}$. (7.7) follows by time-reversal symmetry of the diffusion.

Adding (7.6) and (7.7):

$$\partial_t\rho = -\nabla\!\cdot\!(\mathbf{v}\rho) \tag{7.8}$$

This is the continuity equation — probability is conserved with current $\mathbf{j} = \rho\mathbf{v}$.

Subtracting (7.7) from (7.6):

$$0 = -\nabla\!\cdot\!(\mathbf{u}\rho) + D\nabla^2\rho \tag{7.9}$$

Since $\nabla^2\rho = \nabla\!\cdot\!(\nabla\rho)$, this gives:

$$\boxed{\mathbf{u} = D\,\frac{\nabla\rho}{\rho} = D\,\nabla\ln\rho} \tag{7.10}$$

The osmotic velocity is proportional to the probability gradient — particles are driven from low-density to high-density regions. This is a consequence of the diffusion process, not an assumption.

7.2.3 The Stochastic Newton's Law

For differentiable trajectories, Newton's second law is $m\ddot{\mathbf{x}} = \mathbf{F}$. For stochastic trajectories, the second derivative does not exist. Nelson defined the stochastic acceleration as the time-symmetric second derivative:

$$\mathbf{a} = \frac{1}{2}(D_+D_- + D_-D_+)\mathbf{x}(t) \tag{7.11}$$

and postulated the stochastic Newton's law:

$$m\mathbf{a} = -\nabla V \tag{7.12}$$

Computing $\mathbf{a}$ in terms of $\mathbf{v}$ and $\mathbf{u}$. The Itô calculus rules for the forward and backward derivatives acting on a function $f(\mathbf{x}(t), t)$ are [19]:

$$D_+ f = \partial_t f + \mathbf{b}_+\!\cdot\!\nabla f + D\nabla^2 f \tag{7.13}$$ $$D_- f = \partial_t f + \mathbf{b}_-\!\cdot\!\nabla f - D\nabla^2 f \tag{7.14}$$

The extra $\pm D\nabla^2 f$ terms arise from the Itô correction: for a Wiener process, $(dx_i)^2 = 2D\,dt$ to leading order, generating second-derivative contributions.

The stochastic acceleration (7.11), computed by applying $D_-$ to $\mathbf{b}_+$ and $D_+$ to $\mathbf{b}_-$, is:

$$D_-\mathbf{b}_+ = D_-(\mathbf{v}+\mathbf{u}) = \partial_t(\mathbf{v}+\mathbf{u}) + [(\mathbf{v}-\mathbf{u})\!\cdot\!\nabla](\mathbf{v}+\mathbf{u}) - D\nabla^2(\mathbf{v}+\mathbf{u}) \tag{7.15a}$$ $$D_+\mathbf{b}_- = D_+(\mathbf{v}-\mathbf{u}) = \partial_t(\mathbf{v}-\mathbf{u}) + [(\mathbf{v}+\mathbf{u})\!\cdot\!\nabla](\mathbf{v}-\mathbf{u}) + D\nabla^2(\mathbf{v}-\mathbf{u}) \tag{7.15b}$$

Taking $\mathbf{a} = \frac{1}{2}[(7.15a)+(7.15b)]$:

$$\mathbf{a} = \partial_t\mathbf{v} + \frac{1}{2}\big\{[(\mathbf{v}-\mathbf{u})\!\cdot\!\nabla](\mathbf{v}+\mathbf{u}) + [(\mathbf{v}+\mathbf{u})\!\cdot\!\nabla](\mathbf{v}-\mathbf{u})\big\} - D\nabla^2\mathbf{u} \tag{7.16}$$

Expanding the braced term: $[(\mathbf{v}-\mathbf{u})\!\cdot\!\nabla](\mathbf{v}+\mathbf{u}) = (\mathbf{v}\!\cdot\!\nabla)\mathbf{v} + (\mathbf{v}\!\cdot\!\nabla)\mathbf{u} - (\mathbf{u}\!\cdot\!\nabla)\mathbf{v} - (\mathbf{u}\!\cdot\!\nabla)\mathbf{u}$. Similarly for the second bracket. The cross terms $(\mathbf{v}\!\cdot\!\nabla)\mathbf{u}$ and $(\mathbf{u}\!\cdot\!\nabla)\mathbf{v}$ cancel pairwise between the two brackets. The $D\nabla^2\mathbf{v}$ terms also cancel ($-D\nabla^2\mathbf{v}$ from (7.15a) and $+D\nabla^2\mathbf{v}$ from (7.15b)). The surviving terms give:

$$\boxed{\mathbf{a} = \partial_t\mathbf{v} + (\mathbf{v}\!\cdot\!\nabla)\mathbf{v} - (\mathbf{u}\!\cdot\!\nabla)\mathbf{u} - D\nabla^2\mathbf{u}} \tag{7.17}$$

7.2.4 Derivation of the Schrödinger Equation

Step 1. Assume the current velocity is irrotational: $\mathbf{v} = \nabla S/m$ for a scalar field $S(\mathbf{x},t)$. (For single-particle systems with no magnetic field, this is equivalent to the single-valuedness of the wavefunction.)

Step 2. Write $\rho = R^2$ with $R > 0$. From (7.10):

$$\mathbf{u} = 2D\frac{\nabla R}{R} \tag{7.18}$$

Step 3. Evaluate $(\mathbf{u}\!\cdot\!\nabla)\mathbf{u} + D\nabla^2\mathbf{u}$. Using component notation with $u_j = 2D\partial_j R/R$:

$$\partial_i u_j = 2D\!\left(\frac{\partial_i\partial_j R}{R} - \frac{\partial_i R\,\partial_j R}{R^2}\right) \tag{7.19}$$ $$u_i\partial_i u_j = 4D^2\!\left(\frac{\partial_i R\,\partial_i\partial_j R}{R^2} - \frac{|\nabla R|^2\,\partial_j R}{R^3}\right) \tag{7.20}$$ $$\partial_i^2 u_j = 2D\!\left(\frac{\partial_j\nabla^2 R}{R} - \frac{\nabla^2 R\,\partial_j R}{R^2} - \frac{2\partial_i R\,\partial_i\partial_j R}{R^2} + \frac{2|\nabla R|^2\,\partial_j R}{R^3}\right) \tag{7.21}$$

Adding (7.20) and $D\times$(7.21):

$$u_i\partial_i u_j + D\partial_i^2 u_j = 2D^2\!\left(\frac{\partial_j\nabla^2 R}{R} - \frac{\nabla^2 R\,\partial_j R}{R^2}\right) = 2D^2\partial_j\!\left(\frac{\nabla^2 R}{R}\right) \tag{7.22}$$

(All other terms cancel in pairs.)

Step 4. Substitute into (7.17) with $m\mathbf{a} = -\nabla V$:

$$\nabla\!\left[\partial_t S + \frac{|\nabla S|^2}{2m} + V - 2mD^2\frac{\nabla^2 R}{R}\right] = 0 \tag{7.23}$$

where we used $(\mathbf{v}\!\cdot\!\nabla)\mathbf{v} = \nabla|\nabla S|^2/(2m^2)$. Since the gradient vanishes, the bracket equals a function of $t$ alone, absorbable into $S$:

$$\frac{\partial S}{\partial t} + \frac{|\nabla S|^2}{2m} + V - 2mD^2\frac{\nabla^2 R}{R} = 0 \tag{7.24}$$

Step 5. Set $\hbar = 2mD$ (derived in Section 7.3) and define $\psi = R\,e^{iS/\hbar}$.

7.3 Deriving $D = \hbar/(2m)$ from the Ether ZPF

This is the result that the ether framework uniquely provides: the derivation of the diffusion coefficient — and hence Planck's constant — from the ZPF.

7.3.1 The Argument

1. Boyer proved (Theorem 6.1) that a harmonic oscillator in the ZPF reaches a stationary state with $\rho_0(x) = |\psi_0(x)|^2$ and zero probability current ($\mathbf{v} = 0$).
2. In Nelson's framework, this stationary state has osmotic velocity $\mathbf{u} = D\nabla\ln\rho_0$ and stochastic acceleration (7.17) that must satisfy Newton's law.
3. The requirement that Newton's law holds uniquely determines $D$.

7.3.2 The Derivation

Consider the harmonic oscillator, $V = \frac{1}{2}m\omega_0^2 x^2$.

Given (from Boyer, Theorem 6.1):

$$\rho_0(x) = \sqrt{\frac{m\omega_0}{\pi\hbar}}\exp\!\left(-\frac{m\omega_0 x^2}{\hbar}\right), \qquad \mathbf{v} = 0 \tag{7.28}$$

Osmotic velocity (from 6.10):

$$u = D\frac{d\ln\rho_0}{dx} = -\frac{2Dm\omega_0}{\hbar}\,x \tag{7.29}$$

Stochastic acceleration (from 6.17 with $\mathbf{v} = 0$, $\partial_t\mathbf{v} = 0$):

$$a = -u\frac{du}{dx} - D\frac{d^2u}{dx^2} \tag{7.30}$$

Since $u$ is linear in $x$: $du/dx = -2Dm\omega_0/\hbar$ (constant) and $d^2u/dx^2 = 0$. Therefore:

$$a = -\!\left(-\frac{2Dm\omega_0}{\hbar}x\right)\!\left(-\frac{2Dm\omega_0}{\hbar}\right) = -\frac{4D^2m^2\omega_0^2}{\hbar^2}\,x \tag{7.31}$$

Newton's law requires $ma = -m\omega_0^2 x$:

$$-\frac{4D^2m^2\omega_0^2}{\hbar^2} = -\omega_0^2 \tag{7.32}$$ $$\frac{4D^2m^2}{\hbar^2} = 1 \tag{7.33}$$ $$\boxed{D = \frac{\hbar}{2m}} \tag{7.34}$$

Three critical features of this result:

(a) $\omega_0$ cancels. The oscillator frequency does not appear in $D$. The diffusion coefficient is a property of the ether, not of the system — any oscillator driven by the same ZPF has the same $D/m$ ratio.

(b) Universality. Because the ZPF is universal (it permeates all space and couples to all matter), every massive particle has $D = \hbar/(2m)$. This explains why neutral particles obey the same Schrödinger equation as charged ones: the diffusion is driven by the ether's gravitational/mechanical coupling (Section 3), not solely by electromagnetic coupling. We note, however, that the derivation above is for the harmonic oscillator. The universality claim — that $D = \hbar/(2m)$ holds for arbitrary potentials — is verified for the hydrogen atom (Section 7.3.3) and guaranteed for all systems by the Nelson bridge (Theorem 7.1), but a direct SED derivation for general potentials remains an open problem (see C1 in Section 11.2).

(c) The origin of $\hbar$. Planck's constant enters through Boyer's stationary distribution (7.28), where it was determined by the ZPF spectral density $\rho(\omega) = \hbar\omega^3/(2\pi^2c^3)$. The chain is: ZPF amplitude → stationary distribution → diffusion coefficient → Schrödinger equation. The quantum of action in the Schrödinger equation is the same $\hbar$ that sets the ZPF fluctuation amplitude.

7.3.3 Verification: Hydrogen Atom

We verify that the same $D$ produces the correct dynamics for a system with a completely different potential.

For the hydrogen ground state ($\ell=0$, $\mathbf{v}=0$): $R(r) = (\pi a_0^3)^{-1/2}e^{-r/a_0}$ where $a_0 = 4\pi\epsilon_0\hbar^2/(me^2)$.

The stationary dynamical (7.24) with $\partial_tS = -E$ and $|\nabla S|^2 = 0$ requires:

$$E = V(r) - \frac{\hbar^2}{2m}\frac{\nabla^2 R}{R} \tag{7.35}$$

Computing $\nabla^2 R/R$ in spherical coordinates:

$$\frac{\nabla^2 R}{R} = \frac{R'' + 2R'/r}{R} = \frac{1}{a_0^2} - \frac{2}{a_0 r} \tag{7.36}$$

Substituting:

$$E = -\frac{e^2}{4\pi\epsilon_0 r} - \frac{\hbar^2}{2m}\!\left(\frac{1}{a_0^2} - \frac{2}{a_0 r}\right) = -\frac{e^2}{4\pi\epsilon_0 r} + \frac{e^2}{4\pi\epsilon_0 r} - \frac{\hbar^2}{2ma_0^2} \tag{7.37}$$

using $\hbar^2/(ma_0) = e^2/(4\pi\epsilon_0)$. The Coulomb potential is exactly cancelled, leaving:

$$E = -\frac{me^4}{32\pi^2\epsilon_0^2\hbar^2} = -13.6\;\text{eV} \qquad \tag{7.38}$$

The value $D = \hbar/(2m)$, derived from the harmonic oscillator ZPF, produces the correct hydrogen ground state energy from a completely different potential — confirming universality.

7.4 Consequences

The Schrödinger (7.25) has been derived from ether dynamics. All of non-relativistic quantum mechanics follows.

7.4.1 Excited States and Discrete Spectra

The time-independent Schrödinger equation admits discrete eigenvalues $E_n$ for confining potentials. For hydrogen: $E_n = -13.6\;\text{eV}/n^2$. This resolves the principal limitation of Section 5: SED alone produces ground states but not excited states; the Nelson–SED bridge produces all energy levels.

Each eigenstate $\psi_n$ corresponds to a distinct stationary stochastic process: the particle diffuses in the ether with probability density $R_n^2$ and current velocity $\nabla S_n/m$. The discreteness of $E_n$ arises because only specific diffusion patterns are compatible with normalisability and smoothness — the ether supports resonant modes, like a vibrating membrane.

7.4.2 The Quantum Potential as Ether Diffusion Pressure

The term $Q = -(\hbar^2/2m)\nabla^2R/R$ in (7.24) is the quantum potential of de Broglie–Bohm theory [113]. In that framework, it is postulated. In the ether framework, it is derived:

$$Q = -2mD^2\frac{\nabla^2\sqrt{\rho}}{\sqrt{\rho}} \tag{7.39}$$

This is a diffusion pressure: it measures the curvature of the probability amplitude and produces an effective force $-\nabla Q$ on the particle. It is large where $\rho$ varies rapidly (near nuclei, potential barriers, nodes) and vanishes where $\rho$ is uniform. The "non-locality" of $Q$ — its dependence on the global shape of $\rho$ — reflects the spatial extent of the ether ZPF.

7.5 The Physical Content: What Is Resolved

The wavefunction. $\psi = R\,e^{iS/\hbar}$ encodes the statistical state of a particle diffusing in the ether: $R^2$ is the probability density, $\nabla S/m$ is the mean current velocity. It is neither a physical field nor mere information — it is the complete description of a stochastic process.

The gravity–quantum connection. The same ether whose mean flow is gravity (Part II) produces, through its fluctuations, quantum mechanics (Part IV):

Ether aspect	Physical effect	Mechanism
Mean flow	Gravity	Acoustic metric (Section 3)
Condensate density	Dark matter	Superfluid self-interaction (Section 4.2)
Phonon ZPF energy	Dark energy	Lorentz-invariant ZPF (Section 4.3)
EM constitutive response	Plasma dielectric	Charge-driven ether polarisation (Section 5)
EM mode fluctuations	Quantum ground states	SED (Section 6)
EM fluctuation-driven diffusion	Schrödinger equation	Nelson mechanics (Section 7)

This is not an analogy. It is a single medium with quantitative predictions across all six domains.

7.6 The de Broglie–Bohm Guidance Equation from Ether Dynamics

The Schrödinger equation (Theorem 7.1) was derived from the ether's stochastic dynamics. We now show that this derivation automatically encompasses the de Broglie–Bohm pilot wave theory [113] — the deterministic quantum mechanics in which particles follow definite trajectories guided by the wave function. The connection is exact: the current velocity of the Nelson diffusion process (Section 7.2.1, (7.3)) IS the de Broglie–Bohm guidance velocity. The quantum potential ((7.39)) IS the de Broglie–Bohm quantum potential. No additional postulates are required.

This result completes a circle that began with Louis de Broglie in 1927. De Broglie proposed that quantum particles are real objects guided by real waves propagating through a physical medium — the quantum ether. His pilot wave theory was suppressed at the 1927 Solvay Conference and abandoned for twenty-five years until David Bohm reconstructed it in 1952. John Bell, the strongest advocate of both pilot wave theory and the Lorentzian ether, asked: "Why is the pilot wave picture ignored in textbooks? Should it not be taught, not as the only way, but as an antidote to the prevailing complacency?" [126]. The ether framework answers Bell's question by deriving the pilot wave from ether physics.

7.6.1 The Madelung Equations

Equations (7.8) and (7.24), derived in Section 7.2 as consequences of Nelson's stochastic dynamics, constitute the Madelung equations — first obtained by Erwin Madelung in 1926 [114] by substituting the polar decomposition $\psi = R\,e^{iS/\hbar}$ into the Schrödinger equation and separating real and imaginary parts. We collect them here for reference:

Continuity equation ((7.8)):

$$\frac{\partial\rho}{\partial t} + \nabla\!\cdot\!\left(\rho\,\frac{\nabla S}{m}\right) = 0 \tag{7.40}$$

where $\rho = R^2 = |\psi|^2$ is the probability density.

Modified Hamilton–Jacobi equation ((7.24)):

$$\frac{\partial S}{\partial t} + \frac{|\nabla S|^2}{2m} + V + Q = 0 \tag{7.41}$$

where the quantum potential is ((7.39)):

$$Q = -\frac{\hbar^2}{2m}\,\frac{\nabla^2 R}{R} = -2mD^2\,\frac{\nabla^2\!\sqrt{\rho}}{\sqrt{\rho}} \tag{7.42}$$

These three equations are mathematically equivalent to the Schrödinger equation (Theorem 7.1). They have the form of the Euler equations for an irrotational, compressible fluid with density $\rho$, velocity $\mathbf{v} = \nabla S/m$, and an internal "quantum pressure" encoded in $Q$.

7.6.2 Identification with the de Broglie–Bohm Guidance Equation

The Madelung velocity field $\mathbf{v} = \nabla S/m$ is identical to the de Broglie–Bohm guidance velocity. We derive the standard Bohmian form.

Step 1. Start from the polar decomposition $\psi = R\,e^{iS/\hbar}$, where $R(\mathbf{x},t) \geq 0$ and $S(\mathbf{x},t)$ are real.

Step 2. Compute $\nabla\psi$ by the product rule and chain rule:

$$\nabla\psi = (\nabla R)\,e^{iS/\hbar} + R\,\nabla(e^{iS/\hbar}) = (\nabla R)\,e^{iS/\hbar} + R\cdot\frac{i}{\hbar}(\nabla S)\,e^{iS/\hbar} \tag{7.43}$$

Factoring $e^{iS/\hbar}$:

$$\nabla\psi = \left[\nabla R + \frac{i}{\hbar}\,R\,\nabla S\right]e^{iS/\hbar} \tag{7.44}$$

Step 3. Divide by $\psi = R\,e^{iS/\hbar}$. The exponential factors cancel ($e^{iS/\hbar} \neq 0$):

$$\frac{\nabla\psi}{\psi} = \frac{\nabla R + (i/\hbar)\,R\,\nabla S}{R} = \frac{\nabla R}{R} + \frac{i}{\hbar}\,\nabla S \tag{7.45}$$

Since $R$ and $S$ are real, the decomposition into real and imaginary parts is immediate:

$$\text{Re}\!\left(\frac{\nabla\psi}{\psi}\right) = \frac{\nabla R}{R}, \qquad \text{Im}\!\left(\frac{\nabla\psi}{\psi}\right) = \frac{\nabla S}{\hbar} \tag{7.46}$$

Step 4. Multiply the imaginary part by $\hbar/m$:

$$\frac{\hbar}{m}\,\text{Im}\!\left(\frac{\nabla\psi}{\psi}\right) = \frac{\hbar}{m}\cdot\frac{\nabla S}{\hbar} = \frac{\nabla S}{m} \tag{7.47}$$

The factors of $\hbar$ cancel exactly.

Step 5. Therefore:

$$\boxed{\frac{d\mathbf{x}}{dt} = \frac{1}{m}\,\nabla S(\mathbf{x}, t) = \frac{\hbar}{m}\,\text{Im}\!\left(\frac{\nabla\psi(\mathbf{x},t)}{\psi(\mathbf{x},t)}\right)} \tag{7.48}$$

This is the de Broglie–Bohm guidance equation [113]. It states that a particle at position $\mathbf{x}$ moves with velocity determined by the gradient of the wave function's phase. (7.48) is not an additional postulate — it is the current velocity $\mathbf{v}_c = \nabla S/m$ (Section 7.2.1) that was derived from the ether's stochastic dynamics in Section 7.2, now expressed in the equivalent wave-function form.

where $\mathbf{u}$ is the osmotic velocity ((7.10)) and $D = \hbar/(2m)$ (Section 7.3). Thus the complex ratio $\nabla\psi/\psi$ encodes both Nelson velocities:

$$\frac{\hbar}{m}\,\frac{\nabla\psi}{\psi} = \mathbf{u} + i\,\mathbf{v}_c \tag{7.50}$$

The current velocity (de Broglie–Bohm guidance) is the imaginary part; the osmotic velocity (Nelson diffusion) is the real part.

Dimensional check. $[\hbar/m] = [\text{m}^2/\text{s}]$. $[\nabla\psi/\psi] = [\text{m}^{-1}]$ (dimensions cancel in the ratio). Product: $[\text{m/s}]$ = velocity.

Verification (plane wave). For $\psi = A\,e^{i\mathbf{p}\cdot\mathbf{x}/\hbar}$: $R = A$ (constant), $S = \mathbf{p}\cdot\mathbf{x}$. From (7.48): $\mathbf{v} = \nabla S/m = \mathbf{p}/m$. This is the de Broglie velocity.

Verification (real wave function). If $\psi$ is purely real, then $S = 0$ (or $S = n\pi\hbar$), so $\nabla S = 0$ and $\mathbf{v} = 0$. The particle is at rest — consistent with a stationary state where probability flows but the current velocity vanishes.

7.6.3 The Quantum Force as Ether Pressure Gradient

The modified Hamilton–Jacobi (7.41) implies that particles following the guidance (7.48) obey a quantum Newton's law. We derive this step by step.

Step 1. Take the gradient of (7.41):

$$\nabla\!\left(\frac{\partial S}{\partial t}\right) + \nabla\!\left(\frac{|\nabla S|^2}{2m}\right) + \nabla V + \nabla Q = 0 \tag{7.51}$$

Since $\nabla$ and $\partial/\partial t$ commute (they act on different variables):

$$\frac{\partial(\nabla S)}{\partial t} + \nabla\!\left(\frac{|\nabla S|^2}{2m}\right) + \nabla V + \nabla Q = 0 \tag{7.52}$$

Step 2. Express each term using $\mathbf{v} = \nabla S/m$, so $\nabla S = m\mathbf{v}$:

First term: $\partial(\nabla S)/\partial t = m\,\partial\mathbf{v}/\partial t$ (since $m$ is constant).

Second term: $|\nabla S|^2/(2m) = m^2|\mathbf{v}|^2/(2m) = m|\mathbf{v}|^2/2$, so $\nabla(|\nabla S|^2/(2m)) = m\,\nabla(|\mathbf{v}|^2/2)$.

Substituting into (7.52):

$$m\,\frac{\partial\mathbf{v}}{\partial t} + m\,\nabla\!\left(\frac{|\mathbf{v}|^2}{2}\right) + \nabla V + \nabla Q = 0 \tag{7.53}$$

Step 3. For irrotational flow ($\nabla \times \mathbf{v} = 0$, which holds because $\mathbf{v} = \nabla S/m$ is a gradient), the standard vector identity gives:

$$(\mathbf{v}\cdot\nabla)\mathbf{v} = \nabla\!\left(\frac{|\mathbf{v}|^2}{2}\right) - \mathbf{v}\times(\nabla\times\mathbf{v}) = \nabla\!\left(\frac{|\mathbf{v}|^2}{2}\right) \tag{7.54}$$

The cross-product term vanishes because $\nabla\times\mathbf{v} = \nabla\times(\nabla S/m) = 0$ (curl of a gradient is identically zero).

Step 4. Substitute (7.54) into (7.53):

$$m\,\frac{\partial\mathbf{v}}{\partial t} + m\,(\mathbf{v}\cdot\nabla)\mathbf{v} + \nabla V + \nabla Q = 0 \tag{7.55}$$

Step 5. Recognise the material derivative $d\mathbf{v}/dt \equiv \partial\mathbf{v}/\partial t + (\mathbf{v}\cdot\nabla)\mathbf{v}$ — the acceleration of a fluid element following the flow:

$$\boxed{m\,\frac{d\mathbf{v}}{dt} = -\nabla V - \nabla Q} \tag{7.56}$$

This is Newton's second law with a quantum force $\mathbf{F}_Q = -\nabla Q$. A particle following the de Broglie–Bohm trajectory experiences the classical potential gradient $-\nabla V$ plus the quantum force $-\nabla Q$. The quantum potential $Q$ ((7.42)) depends on the shape of the probability amplitude ($\nabla^2 R/R$), not on its magnitude — explaining why quantum effects do not diminish with distance from the source.

Physical interpretation in the ether framework. In the Madelung hydrodynamic picture, (7.56) is the Euler equation for a fluid with an internal pressure. In a classical fluid, the Euler equation has a pressure gradient $-\nabla p/\rho$. In the ether's quantum fluid, $Q$ plays the role of this internal pressure: it is the superfluid condensate's resistance to density variations. The quantum force $-\nabla Q$ pushes probability away from regions of high curvature (nodes, barriers) and toward regions of low curvature — producing the interference and tunnelling effects that distinguish quantum from classical mechanics.

Dimensional check. $[m\,d\mathbf{v}/dt] = [\text{kg}\cdot\text{m/s}^2] = [\text{N}]$. $[\nabla V] = [\text{N}]$. $[\nabla Q] = [\text{N}]$ (since $[Q] = [\text{J}]$ and $[\nabla] = [\text{m}^{-1}]$).

Classical limit ($\hbar \to 0$). Since $Q \propto \hbar^2$, $Q \to 0$ and (7.56) reduces to $m\,d\mathbf{v}/dt = -\nabla V$ — Newton's second law of classical mechanics.

7.6.4 The Superfluid Velocity Connection

The identification of the Madelung flow with superfluid hydrodynamics provides the physical mechanism for the guidance equation.

In a Bose–Einstein condensate, the order parameter $\Phi(\mathbf{x},t) = \sqrt{n_0}\,e^{i\theta}$ has the same polar structure as the Madelung decomposition $\psi = R\,e^{iS/\hbar}$. The superfluid number current density is [76]:

$$\mathbf{j} = \frac{\hbar}{2Mi}\!\left[\Phi^*\nabla\Phi - \Phi\nabla\Phi^*\right] \tag{7.57}$$

Computation. From $\Phi = \sqrt{n_0}\,e^{i\theta}$:

$$\nabla\Phi = \left[\nabla\!\sqrt{n_0} + i\,(\nabla\theta)\sqrt{n_0}\right]e^{i\theta} \tag{7.58}$$ $$\Phi^*\nabla\Phi = \sqrt{n_0}\,e^{-i\theta}\cdot\left[\nabla\!\sqrt{n_0} + i\,(\nabla\theta)\sqrt{n_0}\right]e^{i\theta} = \sqrt{n_0}\,\nabla\!\sqrt{n_0} + i\,n_0\,\nabla\theta \tag{7.59}$$

Similarly:

$$\Phi\,\nabla\Phi^* = \sqrt{n_0}\,\nabla\!\sqrt{n_0} - i\,n_0\,\nabla\theta \tag{7.60}$$

Subtracting:

$$\Phi^*\nabla\Phi - \Phi\,\nabla\Phi^* = 2i\,n_0\,\nabla\theta \tag{7.61}$$

Substituting into (7.57):

$$\mathbf{j} = \frac{\hbar}{2Mi}\cdot 2i\,n_0\,\nabla\theta = \frac{\hbar\,n_0}{M}\,\nabla\theta \tag{7.62}$$

Since $\mathbf{j} = n_0\,\mathbf{v}_s$, the superfluid velocity is:

$$\mathbf{v}_s = \frac{\hbar}{M}\,\nabla\theta \tag{7.63}$$

Identification. Comparing with the Madelung velocity $\mathbf{v} = \nabla S/m = (\hbar/m)\,\nabla(S/\hbar)$:

$$\mathbf{v} = \frac{\hbar}{m}\,\nabla\!\left(\frac{S}{\hbar}\right) \quad\longleftrightarrow\quad \mathbf{v}_s = \frac{\hbar}{M}\,\nabla\theta \tag{7.64}$$

The identification is:

$$\frac{S}{\hbar} = \theta, \qquad m = M \tag{7.65}$$

The quantum mechanical phase $S/\hbar$ IS the condensate phase $\theta$. The Madelung flow IS the superfluid flow. The de Broglie–Bohm guidance equation IS the equation of motion for excitations carried by the superfluid ether.

Vortex quantisation. Since $\psi$ must be single-valued, $e^{iS/\hbar}$ must return to itself around any closed loop:

$$\oint \nabla\!\left(\frac{S}{\hbar}\right)\cdot d\mathbf{l} = 2\pi n, \qquad n \in \mathbb{Z} \tag{7.66}$$

Therefore:

$$\oint \mathbf{v}\cdot d\mathbf{l} = \frac{\hbar}{m}\cdot 2\pi n = n\,\frac{h}{m} \tag{7.67}$$

This is the Onsager–Feynman quantisation condition for superfluid vortices: circulation is quantised in units of $h/m$. The same condition governs quantised vortices in superfluid helium-4. The ether's quantum vortices and the superfluid's quantised vortices obey the same mathematics because they ARE the same physics.

Dimensional check. $[\hbar/M] = [\text{m}^2/\text{s}]$. $[\nabla\theta] = [\text{m}^{-1}]$ ($\theta$ is dimensionless). Product: $[\text{m/s}]$.

7.6.5 Bohmian Mechanics from Ether Dynamics

Remark (what is not encompassed). Proposition 7.3 applies to single-particle non-relativistic quantum mechanics. For $N$ particles, Bohmian mechanics requires a pilot wave in $3N$-dimensional configuration space, where the guidance equation for particle $k$ is $d\mathbf{Q}_k/dt = (\hbar/m_k)\,\text{Im}[\nabla_k\psi/\psi](\mathbf{Q}_1, \ldots, \mathbf{Q}_N, t)$. The velocity of each particle depends on the instantaneous positions of ALL other particles — this is the non-locality that Bell's theorem requires. The 3-dimensional ether does not straightforwardly support a $3N$-dimensional wave. This is the central open problem for extending the ether's pilot wave to many-body systems (Section 7.7, Section 11.2, open problem I6).

7.7 Open Problems

We state the limitations without euphemism.

The node problem. Excited-state wavefunctions have nodes ($\psi = 0$, $\rho = 0$, $\mathbf{u} \to \infty$). The stochastic process has singular drift at nodes; the particle cannot cross nodal surfaces. Configuration space decomposes into nodal domains; the full quantum distribution is a mixture over domains [102, 103]. This is mathematically consistent (all predictions reproduced) but physically unsatisfying: the particle's initial nodal domain is an additional hidden variable. Nelson himself identified this as a serious concern [101].

Bell's theorem. Bell [104] proved that no local hidden variable theory reproduces all quantum predictions. The ether ZPF is non-local — the field configuration $\{\alpha_\lambda(\mathbf{k})\}$ extends across all space — so it is not excluded by Bell's theorem. For continuous-variable systems, SED reproduces quantum entanglement correlations [107]. For spin-$1/2$ systems (the canonical Bell scenario), the complete calculation has not been performed. We do not claim resolution.

Spin. Intrinsic angular momentum is not derived from ether microphysics. This requires either a model of spin as rotational ether modes or extension to relativistic stochastic mechanics [111]. However, the multi-component structure required by Proposition 6.1 (Section 6.6.4) provides a concrete pathway:

This is not a conjecture but an application of Volovik's theorem [153]: in any multi-component condensate whose quasiparticle spectrum has point nodes in momentum space, the excitations near the nodes satisfy the Weyl equation $E(\mathbf{p}) = \pm|\mathbf{c}\cdot\mathbf{p}|$, where $\mathbf{c}$ is an effective velocity matrix. The $\pm$ branches are left and right chirality. The spin-$1/2$ character is topological: it arises from a Berry phase of $\pi$ accumulated around the node and cannot be removed by continuous deformations of the order parameter. Combining two Weyl nodes of opposite chirality yields a Dirac fermion — a massive spin-$1/2$ particle.

The key observation linking this to the ether programme: Corollary 6.2 established that the ether must have multi-component structure to support transverse modes at atomic frequencies. This same multi-component structure generically provides the momentum-space nodes from which spin-$1/2$ emerges. The EM cutoff problem and the spin problem are thus not independent — they are two manifestations of the same structural requirement. Specifying the ether's order parameter (vector, spinor, or tensor) would simultaneously resolve both.

The ether's order parameter determines not only spin but also the gauge structure of its low-energy excitations: which gauge bosons emerge as Goldstone modes of the broken symmetry, and what fermion mass spectrum results. We identify this as the deepest open question in the ether programme.

Relativistic quantum mechanics and QFT. The derivation is non-relativistic. The Dirac equation and second quantisation are beyond the present scope.