Parts I–II established that the ether's longitudinal dynamics give rise to gravitational phenomena. The acoustic metric (Theorem 3.1) maps ether flow and density fields to an effective spacetime; the Painlevé–Gullstrand identity (Theorem 3.2) shows this spacetime is exactly Schwarzschild for radial ether inflow at the Newtonian free-fall velocity. The gravitational dielectric equation (Theorem 4.1) shows that the ether's self-interaction produces a modified Poisson equation, with the gravitational permittivity $\mu_e(|\mathbf{g}|/a_0)$ derived from the superfluid condensate's equation of state. The ether's longitudinal constitutive response — how the ether's density responds to the presence of matter — accounts for gravity, dark matter, and dark energy.

This section develops the complementary problem: the ether's transverse constitutive response — how the ether's electromagnetic mode structure responds to the presence of free charges.

The answer is plasma physics. A plasma is a quasi-neutral collection of free charges exhibiting collective electromagnetic behaviour. In the ether framework, a plasma is a region where free charges perturb the ether's transverse (electromagnetic) modes, modifying the local relationship between the fields $\mathbf{E}$, $\mathbf{B}$ and their sources. The result is a frequency-dependent electromagnetic permittivity $\varepsilon_r(\omega, \mathbf{k})$ — the electromagnetic analog of the gravitational permittivity $\mu_e$ of Theorem 4.1.

The parallel is not metaphorical. It is a structural correspondence between two constitutive responses of the same medium (we note that in the electromagnetic case, the plasma charges are absorbed into $\varepsilon_r$, leaving $\rho_{\text{free}}$ as the density of external charges not part of the plasma — analogous to how in Section 4.2.2, the ether's gravitational self-response is absorbed into $\mu_e$):

	Gravitational sector (Section 4)	Electromagnetic sector (this section)
Ether property modified	Density (longitudinal modes)	EM field structure (transverse modes)
Perturbation source	Baryonic matter ($\rho_m$)	Free charges ($n_e, q$)
Unperturbed constitutive relation	$\nabla^2\Phi = 4\pi G\rho_m$ (Poisson)	$\nabla \cdot \mathbf{E} = \rho_{\text{free}}/\varepsilon_0$ (Gauss)
Modified constitutive relation	$\nabla\!\cdot\![\mu_e(g/a_0)\,\mathbf{g}] = -4\pi G\rho_m$	$\nabla\!\cdot\![\varepsilon_r(\omega)\,\mathbf{E}] = \rho_{\text{free}}/\varepsilon_0$
Permittivity function	$\mu_e(x) = 1 - e^{-\sqrt{x}}$ ((4.58))	$\varepsilon_r(\omega) = 1 - \omega_p^2/\omega^2$ (derived in Section 5.4)
Characteristic scale	$a_0 \approx 1.2 \times 10^{-10}$ m/s$^2$	$\omega_p = \sqrt{n_e e^2/(\varepsilon_0 m_{\text{el}})}$
Physical meaning of scale	Acceleration below which ether self-interaction dominates	Frequency below which ether EM response is opaque
Derived from	Superfluid $P(X) = \frac{2\alpha_3}{3}X^{3/2}$ EOS	ZPF-driven charge oscillations (Section 5.4)

The section has a second purpose beyond completing the constitutive description. Plasma physics is the one domain of modern physics where medium-based language — rest frames, dielectric tensors, wave propagation through a material substrate, refractive indices — is not merely tolerated but required. No plasma physicist describes Alfvén waves as "oscillations of abstract field quantities in empty space." The medium is physically present, and the mathematics reflects this. We show that the mathematical formalism of plasma physics is, in its essentials, ether physics applied to charge-dense regions of the ether — and that this observation has specific, quantitative consequences for the ether programme.

5.1 The Ether's Electromagnetic Baseline

5.1.1 Constitutive Properties of the Unperturbed Ether

In the ether framework, the vacuum electromagnetic constitutive relations are properties of the medium. Maxwell's equations in free space ((2.8)–(2.11)) describe wave propagation through the ether at rest, with constitutive parameters:

$$\varepsilon_0 = 8.854 \times 10^{-12}\;\text{F/m}, \qquad \mu_0 = 4\pi \times 10^{-7}\;\text{H/m} \tag{5.1}$$

These are not merely unit-conversion factors. They are the ether's electromagnetic permittivity and permeability — the transverse analogs of the ether's compressibility and density that govern the longitudinal (phonon) modes. The ether interpretation table of Section 2.3 makes this explicit: $\varepsilon_0$ measures the ether's capacity for electric displacement; $\mu_0$ measures its capacity for magnetic flux; their product determines the propagation speed of transverse disturbances:

$$c = \frac{1}{\sqrt{\varepsilon_0\mu_0}} = 2.998 \times 10^8\;\text{m/s} \tag{5.2}$$

The vacuum impedance of the ether is:

$$Z_0 = \sqrt{\frac{\mu_0}{\varepsilon_0}} = 376.73\;\Omega \tag{5.3}$$

This is a measurable material property of the medium — the ratio of electric to magnetic field amplitudes in a propagating wave — with dimensions identical to an electrical impedance and a value that characterises the ether's transverse mechanical response.

In the unperturbed ether (no free charges), Maxwell's equations with the constitutive relations $\mathbf{D} = \varepsilon_0\mathbf{E}$ and $\mathbf{B} = \mu_0\mathbf{H}$ yield the wave (2.14) with plane-wave solutions:

$$\omega = c|\mathbf{k}| \tag{5.4}$$

This is the dispersion relation of the unperturbed ether's transverse modes: non-dispersive, isotropic, and valid at all wavelengths $\lambda \gg \ell_e$ (where $\ell_e$ is the ether's transverse microstructure scale, Section 3.8). The ether is transparent to all electromagnetic radiation above its microstructure cutoff.

5.1.2 The Ether Rest Frame

The ether possesses a rest frame — the frame in which its bulk velocity vanishes ($\mathbf{u} = 0$) and the ZPF is maximally isotropic. As established in Section 4.1.3, this rest frame is identified with the CMB rest frame:

$$v_{\text{CMB}} = (369.82 \pm 0.11)\;\text{km/s} \tag{5.5}$$

In the ether rest frame, Maxwell's equations hold exactly (Lorentz's Assumption 4, Section 2.4). In a frame moving at velocity $\mathbf{v}$ relative to the ether, the equations retain their form — this is the content of Poincaré's theorem (Section 2.5) — but the ether's internal state is Lorentz-boosted. For non-relativistic motion ($v \ll c$), the transformation effects are negligible, and Maxwell's equations apply with the same constitutive parameters $\varepsilon_0$, $\mu_0$ to high accuracy.

For the plasma physics of this section, the relevant rest frame is the plasma rest frame — the frame in which the bulk plasma velocity vanishes and the net current density is zero:

$$\mathbf{J}_{\text{bulk}} = \sum_s n_s q_s \mathbf{v}_s = 0 \tag{5.6}$$

where the sum runs over all charged species $s$ with number density $n_s$, charge $q_s$, and mean velocity $\mathbf{v}_s$. In standard plasma physics, this frame is distinguished operationally (as the frame of zero bulk current) but has no ontological priority. In the ether framework, it acquires physical significance: it is (locally) the frame in which the ether's electromagnetic response is isotropic and the transverse ZPF is unshifted.

Remark. The identification of the plasma rest frame with a locally preferred ether frame is structurally identical to the identification of the PG time coordinate with the proper time of free-falling observers (Section 3.3). In both cases, the ether provides a physical frame that the standard framework treats as merely one coordinate choice among many. In both cases, the physics is identical regardless of interpretation — but the ether interpretation is more economical in that it provides a physical criterion for selecting the natural frame, rather than requiring an arbitrary choice followed by a proof of frame-independence.

5.2 Plasma as a Perturbed Ether State: Definitions

5.2.1 Free Charges in the Ether

A plasma consists of free charges — particles with charge $q_s$ and mass $m_s$ for species $s$ — embedded in the ether. In the ether framework, these charges are coupled to the transverse (EM) ZPF as described in Section 6. Boyer's analysis (Theorem 6.1) shows that a single charged oscillator in the ZPF reaches a stationary state with quantum ground-state energy. The SED harmonic oscillator equation of motion ((6.7)):

$$m\ddot{x} = -m\omega_0^2 x + m\tau\dddot{x} + eE_{\text{ZPF}}(t) \tag{5.7}$$

describes a single charge interacting with the ether's EM fluctuations.

A plasma is the many-body generalisation of this picture. When many free charges are present, each charge responds to the total electric field — the sum of the ZPF, the fields of all other charges, and any externally applied fields. The collective response of the charges modifies the ether's electromagnetic mode structure. Modes that propagated freely through the unperturbed ether ((5.4)) now experience a medium with a frequency-dependent, and potentially anisotropic, electromagnetic response.

(i) Quasi-neutrality — the charge densities satisfy $\sum_s n_s q_s \approx 0$ on scales larger than the Debye length $\lambda_D$, where:

$$\lambda_D = \sqrt{\frac{\varepsilon_0 k_B T_e}{n_e e^2}} \tag{5.8}$$

with $T_e$ the electron temperature and $n_e$ the electron number density.

(ii) Collective response — the plasma frequency exceeds the collision frequency, $\omega_p \gg \nu_{\text{coll}}$, so that the charges respond collectively to electromagnetic perturbations via the ether's long-range field structure rather than through pairwise collisions.

(iii) Statistical validity — the number of charges within a Debye sphere is large, $N_D = \frac{4}{3}\pi n_e \lambda_D^3 \gg 1$, ensuring that the collective description is statistically well-defined.

These three conditions are standard [145, 146]. In the ether framework, they acquire an additional physical meaning: they ensure that the ether's transverse mode structure is modified collectively by the charges, producing a well-defined dielectric response $\varepsilon_r(\omega, \mathbf{k})$, rather than being disrupted incoherently.

5.2.2 The Electromagnetic Permittivity of the Perturbed Ether

The presence of free charges modifies the ether's constitutive relation. The unperturbed relation $\mathbf{D} = \varepsilon_0\mathbf{E}$ generalises to:

$$\mathbf{D} = \varepsilon_0\,\boldsymbol{\varepsilon}_r(\omega, \mathbf{k})\cdot\mathbf{E} \tag{5.9}$$

where $\boldsymbol{\varepsilon}_r$ is the dielectric tensor of the perturbed ether. In the absence of a background magnetic field and for an isotropic charge distribution, $\boldsymbol{\varepsilon}_r$ reduces to a scalar:

$$\varepsilon_r(\omega) = 1 - \sum_s\frac{\omega_{p,s}^2}{\omega^2 - \omega_s^2 + i\gamma_s\omega} \tag{5.10}$$

where $\omega_{p,s} = \sqrt{n_s q_s^2/(\varepsilon_0 m_s)}$ is the plasma frequency of species $s$, $\omega_s$ is the natural frequency of any restoring force (zero for free charges), and $\gamma_s$ is a damping rate. This general form is derived from the ether's charge dynamics in Section 5.4.3 ((5.29f)); we state it here to establish notation.

For a cold, collisionless, unmagnetised electron plasma with immobile ions ($\omega_s = 0$, $\gamma_s = 0$, dominant species: electrons), this reduces to:

$$\varepsilon_r(\omega) = 1 - \frac{\omega_p^2}{\omega^2} \tag{5.11}$$

with the electron plasma frequency:

$$\boxed{\omega_p = \sqrt{\frac{n_e e^2}{\varepsilon_0 m_{\text{el}}}}} \tag{5.12}$$

where $m_{\text{el}} = 9.109 \times 10^{-31}$ kg is the electron mass.

Notation. Throughout this section, $m_s$ denotes the mass of charged species $s$, and $m_{\text{el}}$ denotes the electron mass specifically. These are not to be confused with the ether quantum mass $m_e \sim 1$ eV/$c^2$ of Section 4, which governs the longitudinal (gravitational) sector. The ratio $m_e/m_{\text{el}} \sim 2 \times 10^{-6}$ reflects the vast difference in energy scales between the ether condensate and atomic physics.

(5.11) will be derived from the ether's ZPF-driven charge dynamics in Section 5.4. Here we state it to establish notation; the derivation from ether first principles — starting from the SED equation of motion (5.7) generalised to the many-body case — is the content of the next section. We emphasise that in standard plasma textbooks, (5.11) is derived from the linearised equation of motion for an electron in an oscillating electric field. Our derivation will start from the same equation of motion but within the ether's ZPF framework, establishing that the plasma dielectric function is not merely analogous to but identical with the ether's electromagnetic constitutive response in a charge-dense region.

Physical interpretation. The factor $\omega_p^2/\omega^2$ in (5.11) represents the ether's polarisation response to the free charges. At frequencies $\omega > \omega_p$, the charges cannot follow the field oscillations fast enough to screen them, and $\varepsilon_r < 1$ but positive — the ether is transparent, with a phase velocity exceeding $c$ (as in all dielectric media with $\varepsilon < 1$; the group velocity remains below $c$). At $\omega = \omega_p$, the permittivity vanishes — the ether is at resonance. At $\omega < \omega_p$, the permittivity is negative — the ether is opaque, and electromagnetic waves are evanescent. The plasma frequency is therefore the ether's electromagnetic response cutoff in the charge-dense region: the frequency below which the perturbed ether cannot sustain propagating transverse modes.

This is the electromagnetic analog of the MOND acceleration scale $a_0$, which marks the field strength below which the ether's gravitational response transitions from Newtonian ($\mu_e \approx 1$) to enhanced ($\mu_e \propto \sqrt{g/a_0}$). In both cases, the ether's response to perturbation has a characteristic scale set by the perturbing source's density.

5.2.3 The Magnetised Case: General Dielectric Tensor

In the presence of a background magnetic field $\mathbf{B}_0$, the ether's electromagnetic response becomes anisotropic. Charged particles gyrate around magnetic field lines at the cyclotron frequency:

$$\Omega_s = \frac{q_s B_0}{m_s} \tag{5.13}$$

which is signed: $\Omega_e < 0$ for electrons (charge $q_e = -e$), $\Omega_i > 0$ for positive ions. The magnitude $|\Omega_s| = |q_s|B_0/m_s$ is the cyclotron angular frequency.

The dielectric tensor in a magnetised plasma, expressed in coordinates with $\hat{\mathbf{z}} \parallel \mathbf{B}_0$, takes the Stix form [147]:

$$\boldsymbol{\varepsilon}_r(\omega) = \begin{pmatrix} S & -iD & 0 \\ iD & S & 0 \\ 0 & 0 & P \end{pmatrix} \tag{5.14}$$

where the Stix parameters are:

$$S = 1 - \sum_s\frac{\omega_{p,s}^2}{\omega^2 - \Omega_s^2}, \qquad D = \sum_s\frac{\omega_{p,s}^2\,\Omega_s/\omega}{\omega^2 - \Omega_s^2}, \qquad P = 1 - \sum_s\frac{\omega_{p,s}^2}{\omega^2} \tag{5.15}$$

The off-diagonal element $D$ (which produces Faraday rotation, whistler waves, and cyclotron resonance) arises because the magnetic field breaks the isotropy of the ether's charge response: the charges' gyromotion introduces a handedness.

In the ether interpretation, the Stix tensor (5.14) describes the full anisotropic constitutive response of the ether in a region containing both free charges and a magnetic field. The magnetic field, in the ether picture, is a rotational deformation of the ether (Section 2.3, Maxwell's interpretation table), and the Stix tensor encodes how the ether's response to additional electromagnetic perturbations is modified by this pre-existing deformation.

The derivation of (5.14)–(5.15) from ether dynamics is presented in Section 5.4.

5.3 Programme for This Section

Having established the framework and notation, we state what the remainder of this section derives and what it does not.

What is derived:

(D1) The cold plasma dielectric function $\varepsilon_r(\omega) = 1 - \omega_p^2/\omega^2$ ((5.11)) from the ether's ZPF-driven charge dynamics (Section 5.4, Theorem 5.1). The derivation starts from the SED equation of motion ((5.7)), generalised to the many-body case with self-consistent fields, and shows that the linearised collective response of charges in the ether yields the standard plasma constitutive relation.

(D2) The equivalence between Alfvén wave propagation in a magnetised plasma and transverse wave propagation in Young's elastic ether (Section 5.5, Theorem 5.2). The magnetic tension $B_0^2/\mu_0$ is identified with the effective shear modulus of the magnetised ether.

(D3) The structural identity between the Vlasov equation — the fundamental kinetic equation of plasma physics — and a preferred-frame transport equation for the ether's charge distribution function (Section 5.6).

(D4) The Langmuir wave–ZPF resonance: the identification of Landau damping as a kinetic wave–particle resonance in the ether, structurally analogous to the ZPF–matter energy exchange that maintains atomic ground states (Section 5.6, Theorem 6.1 of Boyer).

What is not derived (open problems):

(O1) The precise relationship between the plasma frequency $\omega_p$ and the ether's transverse microstructure scale $\ell_e$ (Section 3.8). Both are electromagnetic-sector quantities, but their connection requires a theory of the ether's transverse dynamics that has not been developed.

(O2) The full nonlinear plasma response from ether microphysics. The derivations of this section are linearised (small perturbations around a homogeneous background). Nonlinear plasma phenomena — wave breaking, parametric instabilities, magnetic reconnection — require the nonlinear ether dynamics, which is beyond the present scope.

(O3) The origin of the Stix tensor's off-diagonal elements from the ether's rotational microstructure. We derive the Stix tensor from the linearised equations of motion (standard), but the ether-specific interpretation of why a background magnetic field introduces anisotropy — i.e., the transverse ether's response to rotational deformation — is not developed at the microphysical level.

These open problems are consistent with the ether programme's broader limitations. The gravitational sector (Sections 3–4) similarly derives the weak-field results exactly while identifying the strong-field extension as an open problem (Section 3.9.3). The quantum sector (Section 6) derives ground states exactly while identifying excited states as requiring further work (Section 7.4.1). The plasma section follows the same pattern: exact results in the linear regime, with nonlinear and microphysical extensions flagged for future work.

---

5.4 Derivation of the Electromagnetic Dielectric from Ether Dynamics

We now derive the ether's electromagnetic constitutive response in charge-dense regions from first principles. The starting point is the SED equation of motion (Section 6.2), generalised to free charges in a background magnetic field. The result — the plasma dielectric tensor — emerges as the coherent linear response of the ether's charge population to electromagnetic perturbations, with the zero-point field maintaining the equilibrium but not contributing to the coherent response.

5.4.1 The SED Equation of Motion for Free Charges in the Ether

In Section 6.2, the SED equation of motion for a bound charged particle ((6.7)) was:

$$m\ddot{\mathbf{r}} = \mathbf{F}_{\text{binding}} + m\tau\dddot{\mathbf{r}} + q\mathbf{E}_{\text{ZPF}}(\mathbf{r}, t) \tag{6.7}$$

For a free charge of species $s$ (mass $m_s$, charge $q_s$) in a region containing other charges and a background magnetic field $\mathbf{B}_0$, this generalises to:

$$m_s\frac{d\mathbf{v}_s}{dt} = q_s\!\left(\mathbf{E}_{\text{tot}} + \mathbf{v}_s \times \mathbf{B}_0\right) + m_s\tau_s\frac{d^2\mathbf{v}_s}{dt^2} + q_s\mathbf{E}_{\text{ZPF}}(\mathbf{r}_s, t) \tag{5.16}$$

where:

• $\mathbf{E}_{\text{tot}}$ is the total coherent electric field (externally applied + self-consistently induced by the collective motion of all charges)
• $\mathbf{v}_s \times \mathbf{B}_0$ is the magnetic Lorentz force from the background field
• $\tau_s = q_s^2/(6\pi\varepsilon_0 m_s c^3)$ is the radiation reaction time for species $s$
• $\mathbf{E}_{\text{ZPF}}$ is the stochastic zero-point field ((6.2))
• The magnetic ZPF force $q_s\mathbf{v}_s \times \mathbf{B}_{\text{ZPF}}$ is omitted: it is suppressed relative to the electric force by $v_s/c \ll 1$ for non-relativistic particles, and does not contribute at linear order

(5.16) is the many-body SED equation: it couples each charge to the ether's ZPF, to the self-consistent field of all other charges, and to the background magnetic field. We now show that the linear collective response is determined entirely by the deterministic (non-stochastic) terms.

5.4.2 Separation of Coherent and Stochastic Responses

Decompose the velocity of each particle into its mean (coherent) response and a stochastic fluctuation:

$$\mathbf{v}_s = \langle\mathbf{v}_s\rangle + \delta\mathbf{v}_s \tag{5.17}$$

where $\langle\cdot\rangle$ denotes the ensemble average over ZPF realisations. The coherent field $\mathbf{E}_{\text{tot}}$ is the same in every realisation (it is the applied perturbation plus the field induced by the mean charge motion). The ZPF, by construction ((6.3)), has zero mean:

$$\langle\mathbf{E}_{\text{ZPF}}(\mathbf{r}, t)\rangle = 0 \tag{5.18}$$

Taking the ensemble average of (5.16):

$$m_s\frac{d\langle\mathbf{v}_s\rangle}{dt} = q_s\!\left(\mathbf{E}_{\text{tot}} + \langle\mathbf{v}_s\rangle \times \mathbf{B}_0\right) + m_s\tau_s\frac{d^2\langle\mathbf{v}_s\rangle}{dt^2} \tag{5.19}$$

The ZPF term vanishes by (5.18). The radiation reaction term survives but is negligible for all plasma frequencies of interest. For electrons:

$$\omega\tau_{\text{el}} = \omega \cdot \frac{e^2}{6\pi\varepsilon_0 m_{\text{el}} c^3} = \omega \times 6.26 \times 10^{-24}\;\text{s} \tag{5.20}$$

Even for the highest plasma frequencies encountered in nature ($\omega_p \sim 10^{15}$ rad/s in solid-density plasmas), $\omega\tau_{\text{el}} \sim 6 \times 10^{-9} \ll 1$. For ions, the radiation reaction is smaller by a further factor $(m_{\text{el}}/m_i)^2$. We therefore neglect the radiation reaction in the linear response, obtaining:

$$\boxed{m_s\frac{d\mathbf{v}_{s,1}}{dt} = q_s\!\left(\mathbf{E}_1 + \mathbf{v}_{s,1} \times \mathbf{B}_0\right)} \tag{5.21}$$

where we have dropped the angle brackets and written $\mathbf{v}_{s,1} \equiv \langle\mathbf{v}_s\rangle$, $\mathbf{E}_1 \equiv \mathbf{E}_{\text{tot}}$ for the first-order (linear) perturbation quantities.

Remark on the ZPF's role. The zero-point field does not appear in (5.21). This does not mean it is physically irrelevant — it means its role is in maintaining the equilibrium, not in determining the linear response. Without the ZPF, the charges would radiate and lose energy (the classical radiation catastrophe); the ZPF continuously replenishes this energy, maintaining the charges in the stationary state from which the linear perturbation (5.21) departs. This is precisely Boyer's mechanism (Theorem 6.1) applied to the many-body context: the ZPF establishes and maintains the equilibrium; the coherent perturbation acts on top of this equilibrium.

The ZPF re-enters at the kinetic level (Section 5.6), where it produces Landau damping through wave–particle resonance — the plasma analog of the energy balance that maintains atomic ground states.

5.4.3 The Unmagnetised Cold Plasma: Isotropic Dielectric

We first derive the dielectric function for the simplest case: an unmagnetised ($\mathbf{B}_0 = 0$), cold (thermal velocities negligible), collisionless plasma. This is the case stated in (5.11).

Step 1: Velocity response. Setting $\mathbf{B}_0 = 0$ in (5.21) and assuming harmonic time dependence $\mathbf{E}_1 = \mathbf{E}_0\,e^{-i\omega t}$, $\mathbf{v}_{s,1} = \mathbf{v}_{s,0}\,e^{-i\omega t}$:

$$-i\omega\,m_s\,\mathbf{v}_{s,0} = q_s\,\mathbf{E}_0 \tag{5.22}$$ $$\mathbf{v}_{s,0} = \frac{iq_s}{\omega m_s}\,\mathbf{E}_0 \tag{5.23}$$

Step 2: Current density. The macroscopic current density from species $s$ is $\mathbf{J}_s = n_s q_s \mathbf{v}_{s,1}$. Summing over all species:

$$\mathbf{J}_1 = \sum_s n_s q_s \mathbf{v}_{s,1} = \sum_s \frac{in_s q_s^2}{\omega m_s}\,\mathbf{E}_1 \equiv \boldsymbol{\sigma}(\omega)\cdot\mathbf{E}_1 \tag{5.24}$$

where the scalar conductivity is:

$$\sigma(\omega) = \sum_s \frac{in_s q_s^2}{\omega m_s} = \frac{i}{\omega}\sum_s \varepsilon_0\,\omega_{p,s}^2 = \frac{i\varepsilon_0}{\omega}\sum_s \omega_{p,s}^2 \tag{5.25}$$

with $\omega_{p,s}^2 = n_s q_s^2/(\varepsilon_0 m_s)$ the plasma frequency of species $s$.

Step 3: From conductivity to dielectric permittivity. Maxwell's equations for the perturbation give Ampère's law:

$$\nabla \times \mathbf{B}_1 = \mu_0\mathbf{J}_1 + \mu_0\varepsilon_0\frac{\partial\mathbf{E}_1}{\partial t} = \mu_0\left(\sigma\,\mathbf{E}_1 - i\omega\varepsilon_0\,\mathbf{E}_1\right) \tag{5.26}$$

This can be written as:

$$\nabla \times \mathbf{B}_1 = -i\omega\mu_0\varepsilon_0\!\left(1 + \frac{i\sigma}{\varepsilon_0\omega}\right)\mathbf{E}_1 = -i\omega\mu_0\varepsilon_0\,\varepsilon_r(\omega)\,\mathbf{E}_1 \tag{5.27}$$

where we identify the dielectric function:

$$\varepsilon_r(\omega) = 1 + \frac{i\sigma(\omega)}{\varepsilon_0\omega} = 1 + \frac{i}{\varepsilon_0\omega}\cdot\frac{i\varepsilon_0}{\omega}\sum_s\omega_{p,s}^2 = 1 - \sum_s\frac{\omega_{p,s}^2}{\omega^2} \tag{5.28}$$

For a plasma dominated by electrons (ions effectively immobile on the electron timescale, $\omega_{p,i}^2 \ll \omega_{p,e}^2$ for $\omega \gg \omega_{p,i}$):

$$\varepsilon_r(\omega) = 1 - \frac{\omega_p^2}{\omega^2} \tag{5.29}$$

with $\omega_p = \omega_{p,e} = \sqrt{n_e e^2/(\varepsilon_0 m_{\text{el}})}$, confirming (5.11). $\square$

Verification of limits:

(i) $n_e \to 0$ (no charges): $\omega_p \to 0$, $\varepsilon_r \to 1$. The ether's unperturbed constitutive relation is recovered.

(ii) $\omega \to \infty$ (high frequency): $\varepsilon_r \to 1$. The charges cannot respond to arbitrarily rapid oscillations; the ether appears unperturbed.

(iii) $\omega = \omega_p$: $\varepsilon_r = 0$. This is the resonance condition — the ether's polarisation response exactly cancels the applied field.

(iv) $\omega \to 0$ (static limit): $\varepsilon_r \to -\infty$. The charges completely screen any static electric field — Debye shielding in the ether picture.

Generalization to bound charges and collisional damping. The derivation above treats free charges ($\omega_s = 0$, $\gamma_s = 0$). We now derive the general dielectric function ((5.10)) for charges with a restoring force of natural frequency $\omega_s$ and a phenomenological damping rate $\gamma_s$ (representing collisional friction with the background).

The equation of motion for a bound, damped charge in the ether (generalising (5.21) with $\mathbf{B}_0 = 0$) is:

$$m_s\ddot{x}_s = -m_s\omega_s^2 x_s - m_s\gamma_s\dot{x}_s + q_s E_1 \tag{5.29a}$$

In Fourier space ($\partial/\partial t \to -i\omega$):

$$(-\omega^2 + \omega_s^2 - i\omega\gamma_s)\,m_s\hat{x}_s = q_s\hat{E}_1 \tag{5.29b}$$ $$\hat{x}_s = \frac{q_s\hat{E}_1}{m_s(\omega_s^2 - \omega^2 - i\omega\gamma_s)} \tag{5.29c}$$

The velocity response is $\hat{v}_s = -i\omega\hat{x}_s$, giving the conductivity:

$$\sigma_s = n_sq_s\hat{v}_s/\hat{E}_1 = \frac{-i\omega\,n_sq_s^2}{m_s(\omega_s^2 - \omega^2 - i\omega\gamma_s)} \tag{5.29d}$$

From $\varepsilon_r = 1 + i\sigma/(\varepsilon_0\omega)$ ((5.41)):

$$\varepsilon_r = 1 + \frac{i}{\varepsilon_0\omega}\sum_s\frac{-i\omega\,n_sq_s^2}{m_s(\omega_s^2 - \omega^2 - i\omega\gamma_s)} = 1 + \sum_s\frac{n_sq_s^2}{\varepsilon_0 m_s(\omega_s^2 - \omega^2 - i\omega\gamma_s)} \tag{5.29e}$$

Using $\omega_s^2 - \omega^2 - i\omega\gamma_s = -(\omega^2 - \omega_s^2 + i\omega\gamma_s)$ and $\omega_{p,s}^2 = n_sq_s^2/(\varepsilon_0 m_s)$:

$$\varepsilon_r(\omega) = 1 - \sum_s\frac{\omega_{p,s}^2}{\omega^2 - \omega_s^2 + i\gamma_s\omega} \tag{5.29f}$$

This confirms (5.10), now derived from the ether's SED framework. The free-particle limit ($\omega_s = 0$, $\gamma_s = 0$) recovers (5.29). For bound charges in a dielectric medium ($\omega_s > 0$, $\gamma_s > 0$), the permittivity has a resonance at $\omega = \omega_s$ with linewidth $\gamma_s$ — the standard Drude–Lorentz model, here derived from the ether's charge dynamics rather than postulated.

5.4.4 The Magnetised Ether: Derivation of the Stix Tensor

We now derive the full anisotropic dielectric tensor for a cold plasma in a uniform background magnetic field $\mathbf{B}_0 = B_0\hat{\mathbf{z}}$. The derivation follows from (5.21) by solving the coupled equations of motion for each species.

Step 1: Component equations of motion. For $\mathbf{B}_0 = B_0\hat{\mathbf{z}}$, the magnetic force on species $s$ is:

$$(\mathbf{v}_{s,1} \times \mathbf{B}_0)_x = v_{s,y}\,B_0, \qquad (\mathbf{v}_{s,1} \times \mathbf{B}_0)_y = -v_{s,x}\,B_0, \qquad (\mathbf{v}_{s,1} \times \mathbf{B}_0)_z = 0 \tag{5.30}$$

Substituting into (5.21) with harmonic time dependence ($\partial/\partial t \to -i\omega$):

$$-i\omega\,m_s\,v_{s,x} = q_s\!\left(E_x + v_{s,y}\,B_0\right) \tag{5.31a}$$ $$-i\omega\,m_s\,v_{s,y} = q_s\!\left(E_y - v_{s,x}\,B_0\right) \tag{5.31b}$$ $$-i\omega\,m_s\,v_{s,z} = q_s\,E_z \tag{5.31c}$$

Using the signed cyclotron frequency $\Omega_s = q_sB_0/m_s$ ((5.13)), these become:

$$-i\omega\,v_{s,x} = \frac{q_s}{m_s}E_x + \Omega_s\,v_{s,y} \tag{5.32a}$$ $$-i\omega\,v_{s,y} = \frac{q_s}{m_s}E_y - \Omega_s\,v_{s,x} \tag{5.32b}$$ $$v_{s,z} = \frac{iq_s}{\omega m_s}E_z \tag{5.32c}$$

Step 2: Solution of the transverse system. Equations (5.32a,b) form a coupled $2\times 2$ system. In matrix form:

$$\begin{pmatrix} -i\omega & -\Omega_s \\ \Omega_s & -i\omega \end{pmatrix}\begin{pmatrix} v_{s,x} \\ v_{s,y} \end{pmatrix} = \frac{q_s}{m_s}\begin{pmatrix} E_x \\ E_y \end{pmatrix} \tag{5.33}$$

The determinant is:

$$\Delta_s = (-i\omega)^2 - (-\Omega_s)(\Omega_s) = -\omega^2 + \Omega_s^2 = -(\omega^2 - \Omega_s^2) \tag{5.34}$$

Solving by Cramer's rule:

$$v_{s,x} = \frac{q_s}{m_s}\cdot\frac{-i\omega\,E_x + \Omega_s\,E_y}{-(\omega^2 - \Omega_s^2)} = \frac{q_s}{m_s}\cdot\frac{i\omega\,E_x - \Omega_s\,E_y}{\omega^2 - \Omega_s^2} \tag{5.35}$$ $$v_{s,y} = \frac{q_s}{m_s}\cdot\frac{-i\omega\,E_y - \Omega_s\,E_x}{-(\omega^2 - \Omega_s^2)} = \frac{q_s}{m_s}\cdot\frac{i\omega\,E_y + \Omega_s\,E_x}{\omega^2 - \Omega_s^2} \tag{5.36}$$

Verification of (5.35). From Cramer's rule applied to (5.33), $v_{s,x} = (q_s/m_s)\det\begin{pmatrix}E_x & -\Omega_s \\ E_y & -i\omega\end{pmatrix}/\Delta_s = (q_s/m_s)(-i\omega E_x + \Omega_s E_y)/\Delta_s$. With $\Delta_s = -(\omega^2 - \Omega_s^2)$: $v_{s,x} = (q_s/m_s)(i\omega E_x - \Omega_s E_y)/(\omega^2 - \Omega_s^2)$.

Verification of (5.36). $v_{s,y} = (q_s/m_s)\det\begin{pmatrix}-i\omega & E_x \\ \Omega_s & E_y\end{pmatrix}/\Delta_s = (q_s/m_s)(-i\omega E_y - \Omega_s E_x)/\Delta_s = (q_s/m_s)(i\omega E_y + \Omega_s E_x)/(\omega^2 - \Omega_s^2)$.

Step 3: Current density. The current density from species $s$ is $\mathbf{J}_s = n_s q_s \mathbf{v}_{s,1}$. Defining $\omega_{p,s}^2 = n_s q_s^2/(\varepsilon_0 m_s)$:

$$J_{s,x} = n_s q_s v_{s,x} = \varepsilon_0\omega_{p,s}^2\cdot\frac{i\omega\,E_x - \Omega_s\,E_y}{\omega^2 - \Omega_s^2} \tag{5.37}$$ $$J_{s,y} = \varepsilon_0\omega_{p,s}^2\cdot\frac{\Omega_s\,E_x + i\omega\,E_y}{\omega^2 - \Omega_s^2} \tag{5.38}$$ $$J_{s,z} = \varepsilon_0\omega_{p,s}^2\cdot\frac{i}{\omega}\,E_z \tag{5.39}$$

The total current $\mathbf{J}_1 = \sum_s \mathbf{J}_s$ defines the conductivity tensor $\boldsymbol{\sigma}$ through $\mathbf{J}_1 = \boldsymbol{\sigma}\cdot\mathbf{E}_1$:

$$\boldsymbol{\sigma} = \varepsilon_0\sum_s\frac{\omega_{p,s}^2}{\omega^2 - \Omega_s^2}\begin{pmatrix} i\omega & -\Omega_s & 0 \\ \Omega_s & i\omega & 0 \\ 0 & 0 & i(\omega^2-\Omega_s^2)/\omega \end{pmatrix} \tag{5.40}$$

Step 4: Dielectric tensor. From Ampère's law (as in Step 3 of Section 5.4.3), the dielectric tensor is:

$$\boldsymbol{\varepsilon}_r = \mathbf{I} + \frac{i\boldsymbol{\sigma}}{\varepsilon_0\omega} \tag{5.41}$$

Computing each component:

$$(\varepsilon_r)_{xx} = (\varepsilon_r)_{yy} = 1 + \frac{i}{\varepsilon_0\omega}\cdot\varepsilon_0\sum_s\frac{i\omega\,\omega_{p,s}^2}{\omega^2 - \Omega_s^2} = 1 - \sum_s\frac{\omega_{p,s}^2}{\omega^2 - \Omega_s^2} \equiv S \tag{5.42}$$ $$(\varepsilon_r)_{xy} = \frac{i}{\varepsilon_0\omega}\cdot\varepsilon_0\sum_s\frac{-\Omega_s\,\omega_{p,s}^2}{\omega^2 - \Omega_s^2} = -i\sum_s\frac{\omega_{p,s}^2\,\Omega_s/\omega}{\omega^2 - \Omega_s^2} \equiv -iD \tag{5.43}$$ $$(\varepsilon_r)_{yx} = \frac{i}{\varepsilon_0\omega}\cdot\varepsilon_0\sum_s\frac{\Omega_s\,\omega_{p,s}^2}{\omega^2 - \Omega_s^2} = +iD \tag{5.44}$$ $$(\varepsilon_r)_{zz} = 1 + \frac{i}{\varepsilon_0\omega}\cdot\varepsilon_0\sum_s\frac{i\omega_{p,s}^2}{\omega} = 1 - \sum_s\frac{\omega_{p,s}^2}{\omega^2} \equiv P \tag{5.45}$$

All off-diagonal elements involving $z$ vanish: $(\varepsilon_r)_{xz} = (\varepsilon_r)_{zx} = (\varepsilon_r)_{yz} = (\varepsilon_r)_{zy} = 0$.

The complete dielectric tensor is therefore:

$$\boxed{\boldsymbol{\varepsilon}_r(\omega) = \begin{pmatrix} S & -iD & 0 \\ iD & S & 0 \\ 0 & 0 & P \end{pmatrix}} \tag{5.46}$$

with the Stix parameters:

$$\boxed{S = 1 - \sum_s\frac{\omega_{p,s}^2}{\omega^2 - \Omega_s^2}, \qquad D = \sum_s\frac{\omega_{p,s}^2\,\Omega_s/\omega}{\omega^2 - \Omega_s^2}, \qquad P = 1 - \sum_s\frac{\omega_{p,s}^2}{\omega^2}} \tag{5.47}$$

This confirms (5.14)–(5.15), now derived from the ether's SED dynamics rather than stated.

Verification of limits:

(i) $B_0 \to 0$ ($\Omega_s \to 0$): $S \to 1 - \sum_s \omega_{p,s}^2/\omega^2 = P$, $D \to 0$, and the tensor becomes $\varepsilon_r\,\mathbf{I}$ with $\varepsilon_r = P = 1 - \sum_s \omega_{p,s}^2/\omega^2$. This recovers the isotropic result (5.29).

(ii) $n_s \to 0$ ($\omega_{p,s} \to 0$): $S \to 1$, $D \to 0$, $P \to 1$, $\boldsymbol{\varepsilon}_r \to \mathbf{I}$. Unperturbed ether recovered.

(iii) $\omega \to \infty$: all three parameters approach 1, giving $\boldsymbol{\varepsilon}_r \to \mathbf{I}$.

(iv) $\omega = |\Omega_s|$ (cyclotron resonance): $S$ diverges for species $s$ — the ether's response to circularly polarised waves resonates with the gyrating charges. This is the electromagnetic analog of the phonon ZPF's spectral peak at the natural frequency in Boyer's oscillator ((6.21)).

We now state the main result formally.

5.4.5 The Structural Parallel with the Gravitational Dielectric

We now make the parallel between Theorems 4.1 and 5.1 precise.

Theorem 4.1 states that any medium responding to gravitational fields by developing density enhancements produces a modified Poisson equation $\nabla\!\cdot\![\mu_e(g/a_0)\,\mathbf{g}] = -4\pi G\rho_m$. The gravitational permittivity $\mu_e$ was derived from the superfluid condensate equation of state: the ether's longitudinal (density) response to the presence of baryonic mass.

Theorem 5.1 states that a region of the ether containing free charges develops an electromagnetic dielectric tensor $\boldsymbol{\varepsilon}_r(\omega)$, modifying Maxwell's equations from $\nabla\!\cdot\!\mathbf{E} = \rho_{\text{free}}/\varepsilon_0$ to $\nabla\!\cdot\![\boldsymbol{\varepsilon}_r\cdot\mathbf{E}] = \rho_{\text{free,ext}}/\varepsilon_0$ (for external free charges, with the plasma charges absorbed into $\varepsilon_r$). This is the ether's transverse (electromagnetic) response to the presence of free charges.

The key structural identity is:

$$\text{Gravitational:}\quad \mu_e(g/a_0)\,\mathbf{g} \longleftrightarrow \quad\text{Electromagnetic:}\quad \boldsymbol{\varepsilon}_r(\omega)\cdot\mathbf{E} \tag{5.49}$$

Both are constitutive relations of the same physical medium — the ether — applied to different sectors of its excitation spectrum.

There is, however, an important structural difference. The gravitational permittivity $\mu_e$ is nonlinear: it depends on the field strength $|\mathbf{g}|/a_0$, and this nonlinearity is responsible for the MOND phenomenology (flat rotation curves, the RAR). The electromagnetic permittivity $\varepsilon_r$ derived here is linear in the field amplitude (the derivation assumed small perturbations $\mathbf{E}_1$). The nonlinearity of the gravitational case arose from the $X^{3/2}$ equation of state of the superfluid condensate; the linearity of the electromagnetic case arises because we treated the charges as test particles responding to the self-consistent field, without back-reaction on the ether's transverse microstructure.

Whether the ether's electromagnetic response becomes nonlinear at extreme field strengths (approaching the Schwinger critical field $E_{\text{crit}} = m_{\text{el}}^2c^3/(e\hbar) = 1.3 \times 10^{18}$ V/m) is an open question. QED predicts vacuum birefringence at these field strengths [148]; in the ether framework, this would correspond to a nonlinear correction to $\varepsilon_r$, analogous to the nonlinear gravitational permittivity. We flag this as an open problem.

5.4.6 The Role of the ZPF: Equilibrium Maintenance in the Plasma

The derivation of Theorem 5.1 showed that the ZPF does not contribute to the linear dielectric response ((5.18)). We now identify what it does contribute.

Equilibrium maintenance. In classical electrodynamics without ZPF, a plasma at rest is a valid equilibrium: free charges with zero mean velocity and a Maxwellian velocity distribution (in the warm case) or zero velocity (in the cold case). However, classical charges in the ZPF-free theory face two problems:

(a) Bound electrons (in neutral atoms at the plasma boundary) would spiral into nuclei in $\sim 10^{-11}$ s. The ZPF prevents this (Boyer, Theorem 6.1).

(b) Free electrons that momentarily form bound states via three-body recombination would immediately collapse classically. The ZPF maintains the atomic ground states of recombined atoms, ensuring that the ionisation/recombination equilibrium is well-defined.

In the ether framework, the plasma equilibrium state is a ZPF-maintained state, just as the atomic ground state is. The perturbation theory (Theorem 5.1) describes departures from this ZPF-maintained equilibrium.

Connection to the EM cutoff problem. The ZPF spectral density $\rho(\omega) = \hbar\omega^3/(2\pi^2c^3)$ ((6.1)) drives the fluctuations of charges in the plasma at all frequencies up to the EM cutoff $\omega_{\text{max}}^{\text{EM}}$. In the unperturbed ether, $\omega_{\text{max}}^{\text{EM}}$ is governed by the transverse microstructure scale $\ell_e$ (Section 3.8). In a plasma, the charges introduce an additional characteristic frequency — the plasma frequency $\omega_p$ — below which the ether is opaque.

The plasma frequency therefore modifies the ether's EM mode structure: modes with $\omega < \omega_p$ are evanescent and do not propagate. The effective ZPF spectrum in the plasma differs from the vacuum spectrum because the dispersion relation (5.62) modifies the density of states.

Derivation. In vacuum, the density of electromagnetic modes per unit volume per unit frequency (two polarisations) is $g_{\text{vac}}(\omega) = \omega^2/(\pi^2 c^3)$, obtained from $g = 2 \times 4\pi k^2/(2\pi)^3 \times dk/d\omega$ with $k = \omega/c$. In the plasma-ether, $k(\omega) = \sqrt{\omega^2 - \omega_p^2}/c$ (from (5.62)), so:

$$\frac{dk}{d\omega} = \frac{\omega}{c\sqrt{\omega^2 - \omega_p^2}}, \qquad k^2 = \frac{\omega^2 - \omega_p^2}{c^2} \tag{5.50a}$$

The modified density of states is:

$$g_{\text{plasma}}(\omega) = \frac{k^2}{\pi^2}\frac{dk}{d\omega} = \frac{\omega\sqrt{\omega^2 - \omega_p^2}}{\pi^2 c^3} \qquad \text{for } \omega > \omega_p \tag{5.50b}$$

and $g_{\text{plasma}}(\omega) = 0$ for $\omega < \omega_p$ (evanescent modes carry no ZPF energy). In the vacuum limit ($\omega_p \to 0$): $g_{\text{plasma}} \to \omega^2/(\pi^2c^3) = g_{\text{vac}}$.

The effective ZPF spectral energy density in the plasma is:

$$\boxed{\rho_{\text{eff}}(\omega) = \frac{1}{2}\hbar\omega\;g_{\text{plasma}}(\omega) = \frac{\hbar\omega^2\sqrt{\omega^2 - \omega_p^2}}{2\pi^2 c^3}} \qquad \text{for } \omega > \omega_p \tag{5.50}$$

For $\omega \gg \omega_p$: $\rho_{\text{eff}} \to \hbar\omega^3/(2\pi^2c^3)$, recovering the vacuum ZPF spectrum ((6.1)). Near the cutoff ($\omega \to \omega_p^+$): $\rho_{\text{eff}} \to 0$ as $\sqrt{\omega - \omega_p}$ — a smooth onset, not a sharp step. The plasma introduces a low-frequency cutoff on the ether's transverse ZPF, with the cutoff shape determined by the dispersion relation (5.62).

This observation partially addresses the EM cutoff open problem of Section 6.6.3. That section identified the need for a theory of the ether's transverse UV cutoff (at high frequencies). What we find here is a theory of the ether's transverse IR cutoff (at low frequencies) in charge-dense regions: the plasma frequency. The full mode structure of the ether's EM sector has cutoffs at both ends — $\omega_p$ from below (due to charges) and $c/\ell_e$ from above (due to microstructure) — with the ZPF populating the intermediate band.

---

5.5 Wave Propagation in the Perturbed Ether

The dielectric tensor (Theorem 5.1) determines how the ether's electromagnetic mode structure is modified by the presence of free charges. We now derive the wave modes that propagate in this perturbed ether and establish their relationship to the foundational ether physics of Sections 2–4.

5.5.1 The General Dispersion Relation

Maxwell's equations in the perturbed ether, with the plasma current absorbed into the dielectric tensor (Theorem 5.1), are:

$$\nabla \times \mathbf{E}_1 = -\frac{\partial\mathbf{B}_1}{\partial t} \tag{5.51}$$ $$\nabla \times \mathbf{H}_1 = \frac{\partial\mathbf{D}_1}{\partial t} \tag{5.52}$$

with $\mathbf{B}_1 = \mu_0\mathbf{H}_1$ (the plasma is non-magnetic at the frequencies of interest) and $\mathbf{D}_1 = \varepsilon_0\,\boldsymbol{\varepsilon}_r(\omega)\cdot\mathbf{E}_1$ ((5.48)). There are no external free currents — the plasma response is entirely encoded in $\boldsymbol{\varepsilon}_r$.

For plane-wave perturbations $\mathbf{E}_1 = \mathbf{E}_0\,e^{i(\mathbf{k}\cdot\mathbf{r} - \omega t)}$:

$$\mathbf{k} \times \mathbf{E}_1 = \omega\mathbf{B}_1 \tag{5.53}$$ $$\mathbf{k} \times \mathbf{B}_1 = -\frac{\omega}{c^2}\,\boldsymbol{\varepsilon}_r\cdot\mathbf{E}_1 \tag{5.54}$$

where we used $\mu_0\varepsilon_0 = 1/c^2$. Derivation of (5.54). From (5.52) in Fourier space: $i\mathbf{k} \times \mathbf{H}_1 = -i\omega\varepsilon_0\boldsymbol{\varepsilon}_r\cdot\mathbf{E}_1$. With $\mathbf{H}_1 = \mathbf{B}_1/\mu_0$: $i\mathbf{k} \times \mathbf{B}_1/\mu_0 = -i\omega\varepsilon_0\boldsymbol{\varepsilon}_r\cdot\mathbf{E}_1$, hence $\mathbf{k} \times \mathbf{B}_1 = -\omega\mu_0\varepsilon_0\boldsymbol{\varepsilon}_r\cdot\mathbf{E}_1 = -(\omega/c^2)\boldsymbol{\varepsilon}_r\cdot\mathbf{E}_1$.

Substituting (5.53) into (5.54) to eliminate $\mathbf{B}_1$:

$$\mathbf{k} \times (\mathbf{k} \times \mathbf{E}_1) = -\frac{\omega^2}{c^2}\,\boldsymbol{\varepsilon}_r\cdot\mathbf{E}_1 \tag{5.55}$$

Applying the vector identity $\mathbf{k} \times (\mathbf{k} \times \mathbf{E}_1) = \mathbf{k}(\mathbf{k}\cdot\mathbf{E}_1) - k^2\mathbf{E}_1$:

$$\mathbf{k}(\mathbf{k}\cdot\mathbf{E}_1) - k^2\mathbf{E}_1 = -\frac{\omega^2}{c^2}\,\boldsymbol{\varepsilon}_r\cdot\mathbf{E}_1 \tag{5.56}$$

Defining the refractive index vector $\mathbf{n} = c\mathbf{k}/\omega$ (with $n = |\mathbf{n}| = ck/\omega$) and dividing by $\omega^2/c^2$:

$$\mathbf{n}(\mathbf{n}\cdot\mathbf{E}_1) - n^2\mathbf{E}_1 + \boldsymbol{\varepsilon}_r\cdot\mathbf{E}_1 = 0 \tag{5.57}$$

In component form:

$$\boldsymbol{\Lambda}\cdot\mathbf{E}_1 = 0, \qquad \Lambda_{ij} = n_i n_j - n^2\delta_{ij} + (\varepsilon_r)_{ij} \tag{5.58}$$

Non-trivial solutions ($\mathbf{E}_1 \neq 0$) require:

$$\boxed{\det\boldsymbol{\Lambda} = 0} \tag{5.59}$$

This is the general dispersion relation for electromagnetic waves in the perturbed ether. All wave modes — transverse, longitudinal, ordinary, extraordinary, Alfvén, whistler, Bernstein — are solutions of (5.59) with the appropriate dielectric tensor.

5.5.2 Transverse Electromagnetic Waves and the Plasma Cutoff

The unmagnetised case. For $\mathbf{B}_0 = 0$, the dielectric tensor reduces to $\boldsymbol{\varepsilon}_r = \varepsilon_r\,\mathbf{I}$ with $\varepsilon_r = 1 - \omega_p^2/\omega^2$ ((5.29)). Choose $\mathbf{k} = k\hat{\mathbf{z}}$ without loss of generality (the unmagnetised plasma is isotropic). Then $\mathbf{n} = (0, 0, n)$ and:

$$\boldsymbol{\Lambda} = \begin{pmatrix} \varepsilon_r - n^2 & 0 & 0 \\ 0 & \varepsilon_r - n^2 & 0 \\ 0 & 0 & \varepsilon_r \end{pmatrix} \tag{5.60}$$

Derivation of (5.60). From (5.58): $\Lambda_{xx} = n_xn_x - n^2 + \varepsilon_r = 0 - n^2 + \varepsilon_r$; $\Lambda_{zz} = n_zn_z - n^2 + \varepsilon_r = n^2 - n^2 + \varepsilon_r = \varepsilon_r$; all off-diagonal elements vanish since $n_x = n_y = 0$ and $\varepsilon_r$ is scalar.

The dispersion relation $\det\boldsymbol{\Lambda} = 0$ gives:

$$(\varepsilon_r - n^2)^2\,\varepsilon_r = 0 \tag{5.61}$$

This factors into two physically distinct solutions.

Solution (a): Transverse waves ($\varepsilon_r = n^2$, doubly degenerate). The eigenvectors have $E_z = 0$ (from $\Lambda_{zz}E_z = \varepsilon_r E_z = 0$ with $\varepsilon_r \neq 0$), so $\mathbf{E}_1 \perp \mathbf{k}$: these are transverse modes. Restoring $n^2 = c^2k^2/\omega^2$:

$$\boxed{\omega^2 = \omega_p^2 + c^2 k^2} \tag{5.62}$$

This is the electromagnetic wave dispersion relation in the plasma-ether. It has the form of a massive relativistic dispersion ($E^2 = m^2c^4 + p^2c^2$) with the plasma frequency playing the role of the rest-mass energy.

The phase and group velocities are:

$$v_{\text{ph}} = \frac{\omega}{k} = \frac{c}{\sqrt{1 - \omega_p^2/\omega^2}} > c, \qquad v_g = \frac{d\omega}{dk} = c\sqrt{1 - \frac{\omega_p^2}{\omega^2}} < c \tag{5.63}$$

The product $v_{\text{ph}} \cdot v_g = c^2$, as required by Lorentz invariance. The group velocity (at which energy and information propagate) remains below $c$ for all $\omega > \omega_p$, consistent with the causal structure of the ether metric (Theorem 3.3).

The plasma frequency as ether EM cutoff. For $\omega < \omega_p$, (5.62) gives $k^2 = (\omega^2 - \omega_p^2)/c^2 < 0$, so $k$ is imaginary. The wave becomes evanescent:

$$\mathbf{E}_1 \propto e^{-\kappa z}\,e^{-i\omega t}, \qquad \kappa = \frac{\sqrt{\omega_p^2 - \omega^2}}{c} \tag{5.64}$$

with penetration depth $\delta = 1/\kappa = c/\sqrt{\omega_p^2 - \omega^2}$. At $\omega = 0$: $\delta = c/\omega_p$, the electromagnetic skin depth.

The plasma frequency $\omega_p$ is therefore the low-frequency cutoff of the ether's transverse mode spectrum in the charge-dense region. Below $\omega_p$, the ether cannot sustain propagating transverse electromagnetic waves — the charges' collective inertia screens the fields.

This complements the cutoff structure identified in Section 4.3: the healing length $\xi$ provides the high-frequency cutoff of the ether's longitudinal (phonon) mode spectrum at $\omega_{\max}^{\text{phonon}} = c_s/\xi$. The plasma-ether's transverse mode structure has cutoffs at both ends — $\omega_p$ from below (due to charges) and $c/\ell_e$ from above (due to transverse microstructure, Section 3.8) — with the ZPF populating the intermediate band.

Solution (b): Longitudinal oscillations ($\varepsilon_r = 0$). The eigenvector has $\mathbf{E}_1 \parallel \mathbf{k}$: this is an electrostatic charge-density oscillation. From $1 - \omega_p^2/\omega^2 = 0$:

$$\omega = \omega_p \tag{5.65}$$

This is a non-propagating oscillation in the cold limit (no $k$-dependence). It represents the natural resonance frequency of the ether's charge-driven polarisation response — the frequency at which all electrons oscillate collectively in phase. We return to its physical significance in Section 5.5.5.

5.5.3 Waves in the Magnetised Ether: Parallel Propagation

For a magnetised plasma with $\mathbf{B}_0 = B_0\hat{\mathbf{z}}$ and wave vector $\mathbf{k} = k_\parallel\hat{\mathbf{z}}$ (parallel to the magnetic field), the wave equation matrix (5.58) with the Stix tensor (5.46) becomes:

$$\boldsymbol{\Lambda} = \begin{pmatrix} S - n_\parallel^2 & -iD & 0 \\ iD & S - n_\parallel^2 & 0 \\ 0 & 0 & P \end{pmatrix} \tag{5.66}$$

where $n_\parallel = ck_\parallel/\omega$.

Derivation. With $\mathbf{n} = (0, 0, n_\parallel)$: $n_in_j = \text{diag}(0,0,n_\parallel^2)$, so $\Lambda_{xx} = 0 - n_\parallel^2 + S = S - n_\parallel^2$; $\Lambda_{xy} = 0 + (-iD) = -iD$; $\Lambda_{zz} = n_\parallel^2 - n_\parallel^2 + P = P$.

The dispersion relation $\det\boldsymbol{\Lambda} = 0$ gives:

$$\left[(S - n_\parallel^2)^2 - D^2\right]\cdot P = 0 \tag{5.67}$$

5.5.4 The Alfvén Wave and Theorem 5.2

We now derive the Alfvén wave as the low-frequency limit of the magnetised ether dispersion, and establish its formal equivalence with transverse wave propagation in Young's elastic ether.

The MHD limit. Consider frequencies far below the ion cyclotron frequency: $\omega \ll |\Omega_i| \ll |\Omega_e|$.

Step 1: Low-frequency Stix parameters. For $\omega^2 \ll \Omega_s^2$:

$$\frac{\omega_{p,s}^2}{\omega^2 - \Omega_s^2} \approx -\frac{\omega_{p,s}^2}{\Omega_s^2} \tag{5.70}$$

so that $S \approx 1 + \sum_s \omega_{p,s}^2/\Omega_s^2$. We evaluate the ratio for each species using $\omega_{p,s}^2 = n_sq_s^2/(\varepsilon_0 m_s)$ and $\Omega_s^2 = q_s^2B_0^2/m_s^2$:

$$\frac{\omega_{p,s}^2}{\Omega_s^2} = \frac{n_s q_s^2}{\varepsilon_0 m_s}\cdot\frac{m_s^2}{q_s^2 B_0^2} = \frac{n_s m_s}{\varepsilon_0 B_0^2} \tag{5.71}$$

Substituting $\varepsilon_0 = 1/(\mu_0 c^2)$ and summing over all species:

$$\sum_s \frac{\omega_{p,s}^2}{\Omega_s^2} = \frac{\mu_0 c^2}{B_0^2}\sum_s n_s m_s = \frac{\mu_0 c^2\,\rho}{B_0^2} = \frac{c^2}{v_A^2} \tag{5.72}$$

where $\rho = \sum_s n_s m_s$ is the total mass density and we define the Alfvén velocity:

$$\boxed{v_A = \frac{B_0}{\sqrt{\mu_0\rho}}} \tag{5.73}$$

Therefore:

$$S \approx 1 + \frac{c^2}{v_A^2} \tag{5.74}$$

Step 2: Low-frequency limit of $D$. From (5.47), in the limit $\omega^2 \ll \Omega_s^2$:

$$D \approx -\sum_s\frac{\omega_{p,s}^2}{\omega\,\Omega_s} = -\frac{1}{\omega\varepsilon_0 B_0}\sum_s n_s q_s \tag{5.75}$$

Derivation. $D = \sum_s \omega_{p,s}^2\Omega_s/[\omega(\omega^2 - \Omega_s^2)] \approx -\sum_s \omega_{p,s}^2/(\omega\Omega_s)$. Now $\omega_{p,s}^2/\Omega_s = [n_sq_s^2/(\varepsilon_0 m_s)]\cdot[m_s/(q_sB_0)] = n_sq_s/(\varepsilon_0 B_0)$.

By quasi-neutrality ($\sum_s n_sq_s = 0$, Definition 5.1(i)):

$$D \approx 0 \tag{5.76}$$

Step 3: Alfvén wave dispersion. With $D \approx 0$, both the R and L modes ((5.68)–(5.69)) give $n_\parallel^2 = S$:

$$\frac{c^2 k_\parallel^2}{\omega^2} = 1 + \frac{c^2}{v_A^2} = \frac{v_A^2 + c^2}{v_A^2} \tag{5.77}$$ $$\omega^2 = \frac{c^2 v_A^2}{c^2 + v_A^2}\,k_\parallel^2 \tag{5.78}$$

In the non-relativistic MHD limit $v_A \ll c$ (which holds for all laboratory and astrophysical plasmas except relativistic jets and magnetar magnetospheres):

$$\boxed{\omega^2 = k_\parallel^2\,v_A^2} \tag{5.79}$$

This is the shear Alfvén wave: a transverse electromagnetic perturbation propagating along the magnetic field at the Alfvén speed $v_A$.

Independent verification from ideal MHD. We derive the same result directly from the ideal MHD equations as an independent check.

The linearised ideal MHD equations for perturbations about a uniform equilibrium ($\rho_0$, $\mathbf{B}_0 = B_0\hat{\mathbf{z}}$, $\mathbf{v}_0 = 0$) are:

$$\rho_0\frac{\partial\mathbf{v}_1}{\partial t} = -\nabla p_1 + \frac{1}{\mu_0}(\nabla\times\mathbf{B}_1)\times\mathbf{B}_0 \tag{5.80}$$ $$\frac{\partial\mathbf{B}_1}{\partial t} = \nabla\times(\mathbf{v}_1\times\mathbf{B}_0) \tag{5.81}$$

Consider a transverse, incompressible perturbation: $\mathbf{v}_1 = v_{1x}(z,t)\,\hat{\mathbf{x}}$ (velocity perpendicular to $\mathbf{B}_0$, varying only along $\mathbf{B}_0$). The pressure gradient vanishes for incompressible transverse perturbations.

From the induction (5.81). Compute $\mathbf{v}_1 \times \mathbf{B}_0 = v_{1x}\,\hat{\mathbf{x}} \times B_0\hat{\mathbf{z}} = -v_{1x}B_0\,\hat{\mathbf{y}}$ (using $\hat{\mathbf{x}} \times \hat{\mathbf{z}} = -\hat{\mathbf{y}}$). Then $\nabla \times (-v_{1x}B_0\,\hat{\mathbf{y}})$ has components:

$$[\nabla \times (-v_{1x}B_0\,\hat{\mathbf{y}})]_x = -\frac{\partial(-v_{1x}B_0)}{\partial z} = B_0\frac{\partial v_{1x}}{\partial z}$$

and $y$- and $z$-components vanish (since $v_{1x}$ depends only on $z$ and $t$). Therefore:

$$\frac{\partial B_{1x}}{\partial t} = B_0\frac{\partial v_{1x}}{\partial z} \tag{5.82}$$

Verification of the curl. $[\nabla \times \mathbf{A}]_x = \partial A_z/\partial y - \partial A_y/\partial z$ with $A_y = -v_{1x}B_0$, $A_z = 0$: $[\nabla \times \mathbf{A}]_x = 0 - (-B_0\partial v_{1x}/\partial z) = B_0\partial v_{1x}/\partial z$.

From the momentum (5.80). Compute $\nabla \times \mathbf{B}_1$ for $\mathbf{B}_1 = B_{1x}(z,t)\,\hat{\mathbf{x}}$:

$$[\nabla \times \mathbf{B}_1]_y = \frac{\partial B_{1x}}{\partial z} - \frac{\partial B_{1z}}{\partial x} = \frac{\partial B_{1x}}{\partial z}$$

and $x$- and $z$-components vanish. Then:

$$(\nabla \times \mathbf{B}_1) \times \mathbf{B}_0 = \frac{\partial B_{1x}}{\partial z}\,\hat{\mathbf{y}} \times B_0\hat{\mathbf{z}} = B_0\frac{\partial B_{1x}}{\partial z}\,\hat{\mathbf{x}}$$

(using $\hat{\mathbf{y}} \times \hat{\mathbf{z}} = \hat{\mathbf{x}}$). The $x$-component of (5.80) gives:

$$\rho_0\frac{\partial v_{1x}}{\partial t} = \frac{B_0}{\mu_0}\frac{\partial B_{1x}}{\partial z} \tag{5.83}$$

Verification of the curl. $[\nabla \times \mathbf{B}_1]_y = \partial B_{1x}/\partial z - 0 = \partial B_{1x}/\partial z$. $\hat{\mathbf{y}} \times \hat{\mathbf{z}} = \hat{\mathbf{x}}$.

Combining (5.82) and (5.83). Take $\partial/\partial t$ of (5.83) and substitute (5.82):

$$\rho_0\frac{\partial^2 v_{1x}}{\partial t^2} = \frac{B_0}{\mu_0}\frac{\partial}{\partial z}\frac{\partial B_{1x}}{\partial t} = \frac{B_0}{\mu_0}\frac{\partial}{\partial z}\!\left(B_0\frac{\partial v_{1x}}{\partial z}\right) = \frac{B_0^2}{\mu_0}\frac{\partial^2 v_{1x}}{\partial z^2}$$

Therefore:

$$\boxed{\frac{\partial^2 v_{1x}}{\partial t^2} = v_A^2\frac{\partial^2 v_{1x}}{\partial z^2}} \tag{5.84}$$

This is the one-dimensional wave equation with propagation speed $v_A = B_0/\sqrt{\mu_0\rho_0}$, confirming (5.79). The two independent derivations — from the cold plasma dielectric (Route A, (5.79)) and from ideal MHD (Route B, (5.84)) — yield identical results.

The Alfvén–Ether equivalence. The Alfvén velocity (5.73) can be written as:

$$v_A = \sqrt{\frac{B_0^2/\mu_0}{\rho}} \equiv \sqrt{\frac{G_{\text{eff}}}{\rho_{\text{eff}}}} \tag{5.85}$$

where we identify:

• $G_{\text{eff}} = B_0^2/\mu_0$: the magnetic tension — the restoring force per unit area for transverse displacements of magnetic field lines [145, 147]. This quantity has units of pressure (Pa = N/m$^2$) and plays the mechanical role of a shear modulus.

• $\rho_{\text{eff}} = \rho = \sum_s n_s m_s$: the mass density of the charge-carrying medium.

Young's elastic ether (Section 2.1, (2.2)) supports transverse waves at phase velocity:

$$c_{\text{wave}} = \sqrt{\frac{G}{\rho_e}} \tag{2.2}$$

where $G$ is the shear modulus and $\rho_e$ the ether density. The formal identity with (5.85) is exact.

Remark on formal vs. physical identity. The formal identity does not imply that the Alfvén wave IS a vacuum ether wave — the Alfvén speed $v_A$ is typically $10^{-4}c$ to $10^{-2}c$, vastly below the speed of light. What Theorem 5.2 establishes is that the mathematical structure Young postulated for the vacuum (elastic medium supporting transverse waves via shear rigidity) is realised physically in magnetised plasmas, with magnetic tension providing the rigidity that Young could not identify mechanically.

Historical significance. The 19th-century ether was rejected partly because no material medium could simultaneously have negligible density (to avoid impeding planetary motion) and enormous shear rigidity ($G \sim 10^{13}$ Pa from (2.2), to support light waves at $c$). Theorem 5.2 shows that plasma physics routinely operates with exactly this mathematical structure — magnetic tension provides the "rigidity," and the plasma mass density provides the "inertia" — without controversy. The ether paradox was not resolved by abandoning the medium; it was resolved by recognising that the restoring force need not be mechanical elasticity but can be electromagnetic tension. A magnetised plasma is an ether with shear rigidity $B_0^2/\mu_0$.

The wave hierarchy. The three wave types derived in this section form a hierarchy of the ether's transverse excitations:

Wave	Speed	Restoring mechanism	Ether analog
Electromagnetic (5.62)	$c$ (above cutoff)	$\varepsilon_0$, $\mu_0$ (ether transverse rigidity)	Young's elastic ether
Alfvén (5.79)	$v_A = B_0/\sqrt{\mu_0\rho}$	Magnetic tension $B_0^2/\mu_0$	Magnetised elastic ether
Gravitational ((3.42))	$c$	Ether compressibility	Acoustic ether

All three are perturbations of the same physical medium, with propagation speeds determined by the relevant constitutive property.

5.5.5 Longitudinal Oscillations and the ZPF Resonance

The Langmuir oscillation. In the cold unmagnetised case, the longitudinal mode (5.65) oscillates at $\omega = \omega_p$ with no wave-like propagation. This is the natural resonance of the ether's charge-driven polarisation response — all electrons oscillate collectively in phase, displacing from the immobile ion background and creating a restoring electric field.

In a warm plasma (electron thermal velocity $v_{\text{th}} = \sqrt{k_BT_e/m_{\text{el}}}$ non-negligible), thermal pressure provides an additional restoring force. The dispersion relation acquires a $k$-dependence. From the kinetic theory of the Vlasov equation (derived in Section 5.6), the warm-plasma longitudinal dispersion is the Bohm–Gross relation [145, 146]:

$$\omega^2 = \omega_p^2 + 3k^2 v_{\text{th}}^2 \tag{5.87}$$

valid for $k\lambda_D \ll 1$ (long wavelengths compared to the Debye length). This relation is derived from the kinetic dielectric in Section 5.6.4 ((5.109j)). It is the electrostatic analog of the transverse dispersion (5.62), with the thermal velocity playing the role of $c$:

	Transverse EM wave	Longitudinal (Langmuir) wave
Dispersion	$\omega^2 = \omega_p^2 + c^2k^2$	$\omega^2 = \omega_p^2 + 3v_{\text{th}}^2 k^2$
Speed	$c$ (ether EM speed)	$\sqrt{3}\,v_{\text{th}}$ (thermal speed)
Cutoff	$\omega_p$	$\omega_p$
Nature	Transverse, EM	Longitudinal, electrostatic

Both branches share the same cutoff $\omega_p$, below which the ether cannot sustain the respective mode.

Landau damping as wave–ether resonance. The Bohm–Gross relation (5.87) describes only the real part of the frequency. The Vlasov equation (Section 5.6) also yields an imaginary part — Landau damping — arising from wave–particle resonance at $v = \omega/k$ (the phase velocity). For a Maxwellian velocity distribution $f_0(v) \propto \exp(-v^2/2v_{\text{th}}^2)$, the damping rate is [145, 149]:

$$\gamma_L = -\sqrt{\frac{\pi}{8}}\;\frac{\omega_p}{(k\lambda_D)^3}\;\exp\!\left(-\frac{1}{2k^2\lambda_D^2}\right) \tag{5.88}$$

Landau damping is a purely kinetic effect — invisible to the fluid (MHD) description. It arises because particles with velocity $v \approx \omega/k$ exchange energy resonantly with the wave: particles slightly slower than the phase velocity are accelerated (absorbing wave energy); particles slightly faster are decelerated (emitting wave energy). For a Maxwellian, there are more slow particles than fast ones at $v = \omega/k > v_{\text{th}}$, so the net energy flow is from wave to particles — the wave is damped.

The structural correspondence with Boyer's ZPF energy balance. Landau damping is the plasma analog of the SED energy balance that maintains atomic ground states (Boyer, Theorem 6.1). The correspondence is precise:

	Boyer (Theorem 6.1)	Landau damping
Mode	ZPF mode at frequency $\omega_0$	Langmuir wave at phase velocity $\omega/k$
Resonance condition	$\omega = \omega_0$ (frequency match)	$v = \omega/k$ (velocity match)
Energy exchange	Absorption from ZPF at resonance	Absorption from wave at resonance
Balance mechanism	Radiation reaction dissipates; ZPF replenishes	Faster particles emit; slower particles absorb
Net effect	Ground state energy $\hbar\omega_0/2$ maintained	Wave amplitude determined by $\partial f_0/\partial v$ at resonance

In Boyer's analysis, a single charged oscillator absorbs energy from the ZPF at its natural frequency $\omega_0$ and radiates it back via the Abraham–Lorentz reaction, maintaining the quantum ground state through a detailed energy balance. In Landau damping, a continuum of charged particles absorbs energy from a collective wave mode at the resonance velocity $v = \omega/k$, with the damping rate governed by the velocity-space gradient $\partial f_0/\partial v$ at resonance — the kinetic analog of the spectral density that governs Boyer's absorption rate.

Both mechanisms are instances of the same physical process: resonant energy transfer between electromagnetic modes of the ether and charged matter. In the SED framework, the ZPF is the ether's ground-state electromagnetic fluctuation spectrum; the Langmuir wave is a coherent excitation above this ground state. Landau damping describes how the ether's charge population thermalises this excitation back toward the ground state — precisely the relaxation mechanism that the SED programme requires but has not previously connected to plasma kinetic theory.

This connection is developed quantitatively in Section 5.6, where we derive the Vlasov equation within the ether framework and show that the Landau damping rate (5.88) follows from the ether's kinetic response.

Remark. The exponential factor $\exp(-1/(2k^2\lambda_D^2))$ in (5.88) means Landau damping is exponentially weak for long-wavelength perturbations ($k\lambda_D \ll 1$). In the ether picture, this is because the ZPF-maintained equilibrium is robust against long-wavelength disturbances — only perturbations with wavelength approaching the Debye scale can couple efficiently to individual particles. This is analogous to the robustness of Boyer's ground state: the harmonic oscillator's equilibrium energy $\hbar\omega_0/2$ is exactly maintained because the ZPF spectrum is matched to the radiation reaction at all frequencies (Corollary 6.1).

---

5.6 Kinetic Theory in the Ether: The Vlasov Equation

The fluid-level results of Sections 5.4–5.5 describe the ether's electromagnetic response in terms of macroscopic quantities — densities, mean velocities, bulk currents. This is analogous to the MHD description of ether gravity (Section 3): valid at long wavelengths but incomplete at the kinetic level. Just as the Nelson stochastic mechanics of Section 7 provides the kinetic foundation beneath the Schrödinger equation, the Vlasov equation provides the kinetic foundation beneath the plasma dielectric.

5.6.1 The Distribution Function in the Ether

In the ether framework, each charged particle follows a stochastic trajectory determined by the SED equation of motion (5.16). The complete microscopic state of a plasma containing $N$ particles is specified by the $6N$-dimensional phase-space point $\{(\mathbf{r}_i, \mathbf{v}_i)\}_{i=1}^N$. The kinetic description replaces this with the one-particle distribution function:

$$f_s(\mathbf{r}, \mathbf{v}, t)\;d^3r\;d^3v = \text{expected number of particles of species } s \text{ in } d^3r\,d^3v \tag{5.89}$$

normalised to the number density:

$$\int f_s(\mathbf{r}, \mathbf{v}, t)\;d^3v = n_s(\mathbf{r}, t) \tag{5.90}$$

The macroscopic fluid quantities of Sections 5.4–5.5 are velocity moments of $f_s$:

$$n_s = \int f_s\;d^3v, \qquad n_s\mathbf{v}_s = \int \mathbf{v}\,f_s\;d^3v, \qquad P_s = m_s\!\int|\mathbf{v} - \mathbf{v}_s|^2 f_s\;d^3v \tag{5.91}$$

where $\mathbf{v}_s(\mathbf{r}, t)$ is the mean flow velocity of species $s$ (as in (5.6); not to be confused with the Nelson osmotic velocity $\mathbf{u}$ of Section 7). The cold plasma limit of Section 5.4 corresponds to $f_s = n_s\,\delta^3(\mathbf{v} - \mathbf{v}_s)$ — all particles of species $s$ at position $\mathbf{r}$ share the single velocity $\mathbf{v}_s(\mathbf{r}, t)$.

The ether interpretation. In standard kinetic theory, $f_s$ is a statistical object describing an ensemble of possible states. In the ether framework, $f_s$ describes the distribution of charges embedded in the ether and driven by the ZPF. The ensemble average is over ZPF realisations ((6.3)), and the distribution evolves under the combined influence of the coherent fields and the stochastic ZPF.

5.6.2 Derivation of the Vlasov Equation from Ether Dynamics

Each particle of species $s$ obeys the SED equation of motion (5.16). In the collisionless kinetic regime, we make the same approximations as in Section 5.4.2:

(i) Radiation reaction negligible: $\omega\tau_s \ll 1$ ((5.20)).

(ii) ZPF averages to zero in the mean field: $\langle\mathbf{E}_{\text{ZPF}}\rangle = 0$ ((5.18)).

(iii) Self-consistent field approximation: $\mathbf{E}$ and $\mathbf{B}$ are the mean fields sourced by the moments of $f_s$.

Under these approximations, each particle obeys the deterministic equations of motion:

$$\frac{d\mathbf{r}}{dt} = \mathbf{v}, \qquad \frac{d\mathbf{v}}{dt} = \frac{q_s}{m_s}\!\left(\mathbf{E} + \mathbf{v} \times \mathbf{B}\right) \tag{5.92}$$

Phase-space incompressibility. The distribution function evolves by conservation of phase-space density along trajectories (5.92). This requires the phase-space flow to be divergence-free:

$$\frac{\partial\dot{r}_i}{\partial r_i} + \frac{\partial\dot{v}_i}{\partial v_i} = \nabla_\mathbf{r}\cdot\mathbf{v} + \nabla_\mathbf{v}\cdot\mathbf{a} = 0 + \frac{q_s}{m_s}\frac{\partial(\mathbf{v}\times\mathbf{B})_i}{\partial v_i} \tag{5.93}$$

Here $\mathbf{B} = \mathbf{B}(\mathbf{r}, t)$ is the mean field — a function of position and time only, not of velocity — determined self-consistently by (5.96). We compute the velocity divergence of the magnetic force. The $i$-th component of $\mathbf{v}\times\mathbf{B}$ is $\epsilon_{ijk}v_jB_k$. Taking $\partial/\partial v_i$:

$$\frac{\partial}{\partial v_i}(\epsilon_{ijk}v_jB_k) = \epsilon_{ijk}\delta_{ij}B_k = \epsilon_{iik}B_k = 0 \tag{5.94}$$

since $\epsilon_{iik} = 0$ (contraction of symmetric and antisymmetric indices). The phase-space flow is incompressible.

The total time derivative of $f_s$ along trajectories therefore vanishes:

$$\frac{df_s}{dt} = \frac{\partial f_s}{\partial t} + \mathbf{v}\cdot\nabla_\mathbf{r}f_s + \frac{q_s}{m_s}(\mathbf{E} + \mathbf{v}\times\mathbf{B})\cdot\nabla_\mathbf{v}f_s = 0$$

This is the Vlasov equation:

$$\boxed{\frac{\partial f_s}{\partial t} + \mathbf{v}\cdot\nabla f_s + \frac{q_s}{m_s}(\mathbf{E} + \mathbf{v}\times\mathbf{B})\cdot\nabla_\mathbf{v}f_s = 0} \tag{5.95}$$

derived from the ether's SED dynamics with the ZPF averaging to zero in the mean-field limit.

The Vlasov equation is closed self-consistently with Maxwell's equations sourced by the moments of $f_s$:

$$\nabla\cdot\mathbf{E} = \frac{1}{\varepsilon_0}\sum_s q_s\!\int f_s\;d^3v, \qquad \nabla\times\mathbf{B} = \mu_0\sum_s q_s\!\int\mathbf{v}\,f_s\;d^3v + \mu_0\varepsilon_0\frac{\partial\mathbf{E}}{\partial t} \tag{5.96}$$

This Vlasov–Maxwell system is the complete kinetic description of the plasma-ether in the collisionless, radiation-reaction-negligible regime.

The role of the ZPF. The ZPF does not appear in (5.95) but enters at two levels:

(a) Collisional corrections: ZPF-driven fluctuations produce an effective collision operator on the right-hand side of (5.95), of order $\omega_p/N_D$ [150]. For $N_D \gg 1$ (Definition 5.1(iii)), this correction is negligible.

(b) Equilibrium selection: The Vlasov equation admits infinitely many stationary solutions. The ZPF selects the physically realised equilibrium — the Maxwellian $f_0 \propto \exp(-m_sv^2/(2k_BT))$ for thermal systems. This parallels Boyer's analysis (Theorem 6.1): the ZPF picks the physically correct stationary state from among the many kinematically allowed ones.

5.6.3 The Preferred-Frame Structure

The Vlasov (5.95) is written in Newtonian phase space: $(\mathbf{r}, \mathbf{v}, t) \in \mathbb{R}^3 \times \mathbb{R}^3 \times \mathbb{R}$, with time $t$ as a universal parameter. This is a preferred-frame formulation — the equation is most naturally expressed in the rest frame of the ether.

The covariant generalization. A Lorentz-covariant kinetic equation exists: the relativistic Vlasov equation on the 7-dimensional mass-shell hypersurface in 8-dimensional phase space $(x^\mu, p^\mu)$:

$$p^\mu\frac{\partial \mathcal{F}_s}{\partial x^\mu} + q_sF^{\mu\nu}p_\nu\frac{\partial \mathcal{F}_s}{\partial p^\mu} = 0 \tag{5.97}$$

where $\mathcal{F}_s(x^\mu, p^\mu)$ is the covariant distribution function (denoted $\mathcal{F}$ to distinguish from the field tensor $F^{\mu\nu}$), and $F^{\mu\nu}$ is the electromagnetic field tensor, and the mass-shell constraint $g_{\mu\nu}p^\mu p^\nu = -m_s^2c^2$ must be imposed separately.

(5.97) is manifestly Lorentz covariant but substantially more complex than (5.95):

(i) The phase space is 8-dimensional $(x^\mu, p^\mu)$, reduced to 7 by the mass-shell constraint $g_{\mu\nu}p^\mu p^\nu = -m_s^2c^2$. At a given time slice, this yields the same 6 physical degrees of freedom $(\mathbf{r}, \mathbf{p})$ as the non-covariant formulation, but the covariant treatment must carry the constraint explicitly.

(ii) The mass-shell constraint couples spatial and temporal momentum components.

(iii) The field tensor $F^{\mu\nu}$ replaces the separate $\mathbf{E}$ and $\mathbf{B}$ fields.

(iv) Velocity moments require the invariant measure $d^3p/p^0$.

The empirical equivalence. Equations (5.95) and (5.97) make identical predictions for all observable quantities: (5.95) is (5.97) evaluated in the ether rest frame. This is precisely the pattern of Theorem 1.1: the ether framework prefers a specific frame while the relativistic framework maintains covariance, but both yield the same physics.

The pragmatic fact. Every plasma physics textbook [145, P2, P3, P5, 150] writes the Vlasov equation in the form (5.95), not (5.97). Practitioners universally work in the plasma rest frame. This is not merely pedagogical convenience — the physics is most transparent in the preferred frame:

(a) The equilibrium distribution $f_0(\mathbf{v})$ is isotropic only in the plasma rest frame.

(b) The Debye length, plasma frequency, and cyclotron frequency are rest-frame quantities.

(c) Landau damping (Section 5.6.5) is a resonance between wave phase velocity and particle velocity — both defined in the rest frame.

In the ether framework, this is not a computational shortcut but a physical reality: the plasma rest frame is the local ether frame, and the Vlasov equation describes the evolution of charges through the ether's electromagnetic mode structure.

Remark. This observation extends beyond plasma physics. The Boltzmann equation, the Fokker–Planck equation, and the BBGKY hierarchy — the foundational equations of statistical mechanics — are all formulated in preferred-frame phase space. Their covariant generalizations exist but are rarely used. Statistical mechanics, like plasma physics, is most naturally expressed in the language of a physical medium with a rest frame. The ether framework provides the ontological basis for this universally adopted practice.

5.6.4 The Kinetic Dielectric Function

We now derive the electromagnetic dielectric function from the linearised Vlasov equation, recovering the cold-plasma result ((5.29)) as a limiting case and obtaining the warm-plasma generalization that includes Landau damping.

Linearisation. Consider an unmagnetised plasma with equilibrium distribution $f_{0,s}(v)$ (homogeneous, stationary, zero mean flow) and a small one-dimensional electrostatic perturbation along $\hat{\mathbf{z}}$. The linearised Vlasov equation for the perturbation $f_{1,s}$ is:

$$\frac{\partial f_{1,s}}{\partial t} + v\frac{\partial f_{1,s}}{\partial z} + \frac{q_s E_1}{m_s}\frac{\partial f_{0,s}}{\partial v} = 0 \tag{5.98}$$

where $v$ is the $z$-component of velocity, $E_1(z,t)$ is the perturbation electric field, and $f_{0,s}(v)$ is the reduced one-dimensional equilibrium distribution normalised so that $\int f_{0,s}\,dv = n_{0,s}$.

Fourier analysis. Writing $f_{1,s} = \hat{f}_{1,s}(v)\,e^{i(kz - \omega t)}$ and $E_1 = \hat{E}_1\,e^{i(kz - \omega t)}$:

$$(-i\omega + ikv)\hat{f}_{1,s} + \frac{q_s\hat{E}_1}{m_s}\frac{\partial f_{0,s}}{\partial v} = 0 \tag{5.99}$$

Solving:

$$\hat{f}_{1,s}(v) = \frac{q_s\hat{E}_1}{m_s}\cdot\frac{\partial f_{0,s}/\partial v}{i(\omega - kv)} \tag{5.100}$$

Verification of (5.100). From (5.99): $\hat{f}_{1,s} = -[q_s\hat{E}_1/(m_s)] \cdot [\partial f_{0,s}/\partial v]/[i(kv - \omega)] = [q_s\hat{E}_1/(m_s)] \cdot [\partial f_{0,s}/\partial v]/[i(\omega - kv)]$.

Charge density. The perturbation charge density is:

$$\hat{\rho}_1 = \sum_s q_s\!\int_{-\infty}^{\infty}\hat{f}_{1,s}(v)\;dv = \sum_s\frac{q_s^2\hat{E}_1}{m_s}\int\frac{\partial f_{0,s}/\partial v}{i(\omega - kv)}\;dv \tag{5.101}$$

Poisson's equation. In Fourier space:

$$ik\varepsilon_0\hat{E}_1 = \hat{\rho}_1 \tag{5.102}$$

Substituting (5.101) and cancelling $\hat{E}_1$:

$$ik\varepsilon_0 = \sum_s\frac{q_s^2}{m_s}\int\frac{\partial f_{0,s}/\partial v}{i(\omega - kv)}\;dv \tag{5.103}$$

Converting the denominator using $i(\omega - kv) = -ik(v - \omega/k)$:

$$ik\varepsilon_0 = -\frac{1}{ik}\sum_s\frac{q_s^2}{m_s}\int\frac{\partial f_{0,s}/\partial v}{v - \omega/k}\;dv \tag{5.104}$$

Multiplying both sides by $ik$:

$$(ik)^2\varepsilon_0 = -\sum_s\frac{q_s^2}{m_s}\int\frac{\partial f_{0,s}/\partial v}{v - \omega/k}\;dv$$ $$k^2\varepsilon_0 = \sum_s\frac{q_s^2}{m_s}\int\frac{\partial f_{0,s}/\partial v}{v - \omega/k}\;dv \tag{5.105}$$

The kinetic dielectric function. Define $\varepsilon_L(\omega, k)$ through the condition that the self-consistent system (Vlasov + Poisson) has non-trivial solutions if and only if $\varepsilon_L = 0$. Using $\omega_{p,s}^2 = n_{0,s}q_s^2/(\varepsilon_0 m_s)$ and defining the normalised distribution $g_{0,s}(v) = f_{0,s}(v)/n_{0,s}$:

$$\boxed{\varepsilon_L(\omega, k) = 1 - \sum_s\frac{\omega_{p,s}^2}{k^2}\int_{C_L}\frac{g_{0,s}'(v)}{v - \omega/k}\;dv} \tag{5.106}$$

where $g_{0,s}' = dg_{0,s}/dv$ and $C_L$ is the Landau contour (Section 5.6.5).

Derivation of (5.106). From (5.105): $1 = (1/\varepsilon_0 k^2)\sum_s (q_s^2/m_s)\int [\partial f_{0,s}/\partial v]/(v - \omega/k)\,dv = \sum_s [\omega_{p,s}^2/(n_{0,s}k^2)]\int [\partial f_{0,s}/\partial v]/(v - \omega/k)\,dv = \sum_s (\omega_{p,s}^2/k^2)\int g_{0,s}'(v)/(v - \omega/k)\,dv$. Setting this equal to zero: $\varepsilon_L = 1 - \sum_s (\omega_{p,s}^2/k^2)\int g_{0,s}'/(v-\omega/k)\,dv = 0$.

Recovery of the cold-plasma dielectric. For a cold plasma, $g_{0,s}(v) = \delta(v)$ and $g_{0,s}'(v) = \delta'(v)$. The integral becomes:

$$\int\frac{\delta'(v)}{v - \omega/k}\;dv \tag{5.107}$$

Integrating by parts (boundary terms vanish since $\delta(v)$ has compact support):

$$= -\int\delta(v)\,\frac{d}{dv}\!\left(\frac{1}{v - \omega/k}\right)dv = \int\delta(v)\,\frac{1}{(v - \omega/k)^2}\;dv = \frac{1}{(\omega/k)^2} = \frac{k^2}{\omega^2} \tag{5.108}$$

Substituting into (5.106):

$$\varepsilon_L = 1 - \sum_s\frac{\omega_{p,s}^2}{k^2}\cdot\frac{k^2}{\omega^2} = 1 - \sum_s\frac{\omega_{p,s}^2}{\omega^2} \tag{5.109}$$

This recovers (5.29).

Warm-plasma generalization: the Bohm–Gross relation. For a Maxwellian equilibrium $g_0(v) = (\sqrt{2\pi}\,v_{\text{th}})^{-1}\exp(-v^2/(2v_{\text{th}}^2))$, we evaluate the principal-value integral in (5.106) for $\omega/k \gg v_{\text{th}}$ (long wavelengths, $k\lambda_D \ll 1$). Expand the integrand in powers of $kv/\omega$:

$$\frac{1}{v - \omega/k} = -\frac{k}{\omega}\cdot\frac{1}{1 - kv/\omega} = -\frac{k}{\omega}\sum_{n=0}^{\infty}\left(\frac{kv}{\omega}\right)^n \tag{5.109a}$$

valid for $|kv/\omega| < 1$, which holds in the bulk of the Maxwellian since $v \sim v_{\text{th}} \ll \omega/k$. The principal-value integral becomes:

$$\mathcal{P}\!\!\int g_0'(v)\cdot\frac{dv}{v - \omega/k} = -\frac{k}{\omega}\sum_{n=0}^{\infty}\left(\frac{k}{\omega}\right)^n\!\int_{-\infty}^{\infty} v^n\,g_0'(v)\;dv \tag{5.109b}$$

The integrals $I_n \equiv \int v^n g_0'(v)\,dv$ are evaluated by integration by parts. Since $g_0(v) \to 0$ as $v \to \pm\infty$ (Maxwellian), boundary terms vanish and $I_n = [v^n g_0]_{-\infty}^{+\infty} - n\int v^{n-1}g_0\,dv = -n\langle v^{n-1}\rangle_{g_0}$:

$$I_0 = [g_0]_{-\infty}^{+\infty} = 0 \tag{5.109c}$$ $$I_1 = -\int g_0\,dv = -1 \tag{5.109d}$$ $$I_2 = -2\int v\,g_0\,dv = 0 \qquad (\text{odd integrand}) \tag{5.109e}$$ $$I_3 = -3\int v^2 g_0\,dv = -3v_{\text{th}}^2 \tag{5.109f}$$

where we used $\int g_0\,dv = 1$ and $\int v^2 g_0\,dv = v_{\text{th}}^2$ (the variance of the normalised Maxwellian). Higher even moments contribute at order $(k\lambda_D)^4$ and beyond.

Substituting (5.109c–f) into (5.109b), retaining terms through $n = 3$:

$$\mathcal{P}\!\!\int \frac{g_0'(v)}{v - \omega/k}\,dv = -\frac{k}{\omega}\left[0 + \frac{k}{\omega}(-1) + 0 + \frac{k^3}{\omega^3}(-3v_{\text{th}}^2)\right] + \mathcal{O}\!\left(\frac{k^4 v_{\text{th}}^4}{\omega^4}\right) = \frac{k^2}{\omega^2} + \frac{3k^4 v_{\text{th}}^2}{\omega^4} + \cdots \tag{5.109g}$$

Substituting into the real part of $\varepsilon_L$ ((5.106)):

$$\text{Re}[\varepsilon_L] = 1 - \frac{\omega_p^2}{k^2}\left(\frac{k^2}{\omega^2} + \frac{3k^4 v_{\text{th}}^2}{\omega^4}\right) = 1 - \frac{\omega_p^2}{\omega^2} - \frac{3\omega_p^2 k^2 v_{\text{th}}^2}{\omega^4} \tag{5.109h}$$

Setting $\text{Re}[\varepsilon_L] = 0$ and solving perturbatively (writing $\omega^2 = \omega_p^2 + \delta$ with $\delta \ll \omega_p^2$, and approximating $\omega^4 \approx \omega_p^4$ in the last term):

$$\frac{\delta}{\omega_p^2} - \frac{3k^2 v_{\text{th}}^2}{\omega_p^2} = 0 \implies \delta = 3k^2 v_{\text{th}}^2 \tag{5.109i}$$

Therefore:

$$\boxed{\omega^2 = \omega_p^2 + 3k^2 v_{\text{th}}^2} \tag{5.109j}$$

This is the Bohm–Gross dispersion relation [145, 149], now derived from the kinetic dielectric ((5.106)) via explicit evaluation of the principal-value integral for a Maxwellian. The factor of 3 arises because the one-dimensional kinetic integral samples $\langle v^2\rangle = v_{\text{th}}^2$, and the third moment $I_3 = -3v_{\text{th}}^2$ carries the one-dimensional adiabatic index $\gamma_e = 3$.

Verification. In the cold limit ($v_{\text{th}} \to 0$): $\omega \to \omega_p$, recovering (5.65). In the long-wavelength limit ($k \to 0$): $\omega \to \omega_p$.

5.6.5 Landau Damping: Derivation and ZPF Connection

The integral in (5.106) has a pole at $v = \omega/k$. For real $\omega$, this pole lies on the real $v$-axis, and the integral requires a prescription for its evaluation.

The Landau contour. Landau (1946) [149] resolved this by treating the problem as an initial value problem. For a perturbation switched on at $t = 0$, the Laplace transform requires $\text{Im}(\omega) > 0$ for convergence. The pole at $v = \omega/k$ then lies above the real $v$-axis, and the velocity integral runs along the real axis below the pole. As $\text{Im}(\omega) \to 0^+$, the contour $C_L$ is obtained by deforming the real axis to pass below the pole:

$$\int_{C_L}\frac{g(v)}{v - \omega/k}\;dv = \mathcal{P}\!\!\int_{-\infty}^{\infty}\frac{g(v)}{v - \omega/k}\;dv + i\pi\,g(\omega/k) \tag{5.110}$$

where $\mathcal{P}$ denotes the Cauchy principal value.

Derivation of (5.110). For $\text{Im}(\omega) > 0$, the pole $v_0 = \omega/k$ lies in the upper half of the complex $v$-plane. The real-axis integral passes below the pole. As $\text{Im}(\omega) \to 0^+$, the contour is deformed by an infinitesimal semicircle below the pole, traversed clockwise. The semicircle contribution is $-\frac{1}{2}(2\pi i)\,g(v_0) = -i\pi\,g(v_0)$, but with the clockwise orientation reversed: $+i\pi\,g(v_0)$. Equivalently, by the Sokhotski–Plemelj theorem: $\lim_{\eta\to 0^+}1/(v - v_0 - i\eta) = \mathcal{P}/(v - v_0) + i\pi\delta(v - v_0)$, giving (5.110).

The imaginary part of the dielectric. From (5.106) with (5.110), applied to $g_{0,s}'$:

$$\text{Im}[\varepsilon_L(\omega, k)] = -\sum_s\frac{\omega_{p,s}^2}{k^2}\cdot\pi\,g_{0,s}'(\omega/k) = -\frac{\pi}{k^2}\sum_s\omega_{p,s}^2\,g_{0,s}'(\omega/k) \tag{5.111}$$

Sign analysis. For a Maxwellian equilibrium at $v = \omega/k > 0$: $g_{0,s}'(\omega/k) < 0$ (the distribution is monotonically decreasing). Therefore $-\omega_{p,s}^2 g_{0,s}'(\omega/k) > 0$, and $\text{Im}[\varepsilon_L] > 0$.

The Maxwellian equilibrium. For the electron species:

$$g_0(v) = \frac{1}{\sqrt{2\pi}\,v_{\text{th}}}\exp\!\left(-\frac{v^2}{2v_{\text{th}}^2}\right), \qquad g_0'(v) = -\frac{v}{v_{\text{th}}^2}\,g_0(v) \tag{5.112}$$

where $v_{\text{th}} = \sqrt{k_BT_e/m_{\text{el}}}$ is the electron thermal velocity. At the resonance velocity $v = \omega/k \approx \omega_p/k$ (for Langmuir waves with $\omega \approx \omega_p$), using $\lambda_D = v_{\text{th}}/\omega_p$:

$$\frac{\omega_p}{k} = \frac{v_{\text{th}}}{k\lambda_D}, \qquad \frac{\omega_p^2}{2k^2v_{\text{th}}^2} = \frac{1}{2k^2\lambda_D^2} \tag{5.113}$$

Therefore:

$$g_0'\!\left(\frac{\omega_p}{k}\right) = -\frac{1}{k\lambda_D\,v_{\text{th}}^2\,\sqrt{2\pi}}\;\exp\!\left(-\frac{1}{2k^2\lambda_D^2}\right) \tag{5.114}$$

Derivation of (5.114). $g_0'(v) = -(v/v_{\text{th}}^2) \cdot (1/\sqrt{2\pi}v_{\text{th}}) \exp(-v^2/(2v_{\text{th}}^2))$. At $v = v_{\text{th}}/(k\lambda_D)$: the prefactor is $-[v_{\text{th}}/(k\lambda_D)]/(v_{\text{th}}^2\sqrt{2\pi}v_{\text{th}}) = -1/(k\lambda_D\,v_{\text{th}}^2\sqrt{2\pi})$. The exponential: $\exp(-v_{\text{th}}^2/(2k^2\lambda_D^2 v_{\text{th}}^2)) = \exp(-1/(2k^2\lambda_D^2))$.

Substituting (5.114) into (5.111):

$$\text{Im}[\varepsilon_L] = \frac{\pi\omega_p^2}{k^2}\cdot\frac{1}{k\lambda_D\,v_{\text{th}}^2\,\sqrt{2\pi}}\;\exp\!\left(-\frac{1}{2k^2\lambda_D^2}\right) \tag{5.115}$$

Using $\omega_p^2/(k^2 v_{\text{th}}^2) = 1/(k^2\lambda_D^2)$ and $\omega_p/(kv_{\text{th}}) = 1/(k\lambda_D)$:

$$= \frac{\pi}{k^3\lambda_D^3\,\sqrt{2\pi}}\;\exp\!\left(-\frac{1}{2k^2\lambda_D^2}\right) = \frac{\sqrt{\pi/2}}{(k\lambda_D)^3}\;\exp\!\left(-\frac{1}{2k^2\lambda_D^2}\right) \tag{5.116}$$

Derivation of (5.116). $\pi\omega_p^2/(k^3\lambda_D v_{\text{th}}^2\sqrt{2\pi})$. Using $\lambda_D = v_{\text{th}}/\omega_p$: $\pi\omega_p^2/(k^3 \cdot v_{\text{th}}/\omega_p \cdot v_{\text{th}}^2\sqrt{2\pi}) = \pi\omega_p^3/(k^3v_{\text{th}}^3\sqrt{2\pi}) = \sqrt{\pi/2}\;\omega_p^3/(k^3v_{\text{th}}^3) = \sqrt{\pi/2}/(k\lambda_D)^3$.

The Landau damping rate. For weakly damped waves ($|\gamma| \ll \omega_r$), writing $\omega = \omega_r + i\gamma$ and expanding $\varepsilon_L(\omega_r + i\gamma, k) = 0$ to first order in $\gamma$:

$$\text{Im}[\varepsilon_L(\omega_r, k)] + \gamma\,\frac{\partial\text{Re}[\varepsilon_L]}{\partial\omega}\bigg|_{\omega_r} = 0 \tag{5.117}$$

Derivation. $\varepsilon_L(\omega_r + i\gamma) \approx \text{Re}[\varepsilon_L(\omega_r)] + i\text{Im}[\varepsilon_L(\omega_r)] + i\gamma\,\partial\varepsilon_L/\partial\omega|_{\omega_r}$. The real part gives $\omega_r$ from $\text{Re}[\varepsilon_L] = 0$ (the Bohm–Gross relation, derived as (5.109j); the functional form used here is (5.109h)). The imaginary part gives (5.117).

Therefore:

$$\gamma = -\frac{\text{Im}[\varepsilon_L(\omega_r, k)]}{\partial\text{Re}[\varepsilon_L]/\partial\omega\big|_{\omega_r}} \tag{5.118}$$

The real part of $\varepsilon_L$ near $\omega \approx \omega_p$ is $\text{Re}[\varepsilon_L] \approx 1 - \omega_p^2/\omega^2 - 3k^2v_{\text{th}}^2/\omega^2$ (from (5.109h), the Bohm–Gross dispersion). Differentiating:

$$\frac{\partial\text{Re}[\varepsilon_L]}{\partial\omega} = \frac{2\omega_p^2}{\omega^3} + \frac{6k^2v_{\text{th}}^2}{\omega^3} = \frac{2(\omega_p^2 + 3k^2v_{\text{th}}^2)}{\omega^3} \tag{5.119}$$

At $\omega_r \approx \omega_p$ and for $k\lambda_D \ll 1$ ($3k^2v_{\text{th}}^2 \ll \omega_p^2$):

$$\frac{\partial\text{Re}[\varepsilon_L]}{\partial\omega}\bigg|_{\omega_r\approx\omega_p} \approx \frac{2\omega_p^2}{\omega_p^3} = \frac{2}{\omega_p} \tag{5.119b}$$

Substituting (5.116) and (5.119b) into (5.118):

$$\gamma_L = -\frac{\sqrt{\pi/2}\,(k\lambda_D)^{-3}\exp(-1/(2k^2\lambda_D^2))}{2/\omega_p} \tag{5.120}$$ $$\boxed{\gamma_L = -\sqrt{\frac{\pi}{8}}\;\frac{\omega_p}{(k\lambda_D)^3}\;\exp\!\left(-\frac{1}{2k^2\lambda_D^2}\right)} \tag{5.121}$$

This confirms (5.88), now derived from the kinetic theory of the Vlasov equation. The negative sign confirms wave damping (energy flows from wave to particles).

Verification of limits:

(i) $k\lambda_D \to 0$ (long wavelength): $\gamma_L \to 0$ exponentially. The wave is undamped — consistent with the cold-plasma result (5.29), which has no damping.

(ii) $k\lambda_D \sim 1$ (Debye scale): $\gamma_L \sim \omega_p$ — the wave is strongly damped within one oscillation period. Waves with wavelength comparable to the Debye length are rapidly thermalised.

(iii) The damping rate depends on the velocity gradient $g_0'(\omega/k)$ at the resonance velocity ((5.111)). For a bump-on-tail distribution with $g_0'(\omega/k) > 0$, the sign reverses and the wave grows — this is the beam–plasma instability.

Physical mechanism. Particles with velocity $v \approx \omega/k$ exchange energy resonantly with the wave. Those slightly slower than the phase velocity are accelerated (absorbing wave energy); those slightly faster are decelerated (emitting wave energy). For a Maxwellian with $g_0'(\omega/k) < 0$ at $v = \omega/k > 0$ (more slow particles than fast at the resonance), the net transfer is from wave to particles: the wave is damped.

The ether interpretation: Landau damping as ZPF thermalisation. In the ether framework, Landau damping is the kinetic-level manifestation of the same resonant energy exchange that maintains Boyer's quantum ground states (Theorem 6.1). The correspondence is:

Boyer's mechanism (Section 6.2): a single charge oscillator absorbs energy from the ZPF at resonance frequency $\omega_0$ and re-radiates via the Abraham–Lorentz reaction. The balance maintains the ground-state energy $\hbar\omega_0/2$.

Landau damping: a collective wave mode (a coherent excitation above the ZPF ground state) transfers energy to particles at resonance velocity $v = \omega/k$. The mode is thermalised back toward the equilibrium — the state maintained by the ZPF.

The structural correspondence:

$$\underbrace{\text{ZPF mode at } \omega_0}_{\text{Boyer}} \longleftrightarrow \underbrace{\text{Langmuir wave at phase velocity } \omega/k}_{\text{Landau}} \tag{5.122}$$ $$\underbrace{\text{Frequency resonance: } \omega = \omega_0}_{\text{single particle}} \longleftrightarrow \underbrace{\text{Velocity resonance: } v = \omega/k}_{\text{kinetic}} \tag{5.123}$$

Both mechanisms are instances of resonant energy transfer between electromagnetic excitations of the ether and charged matter, with the equilibrium maintained by the ZPF fluctuation spectrum. In Boyer's case, the equilibrium is a single-particle quantum ground state. In Landau's case, the equilibrium is a many-body velocity distribution. Both are maintained by the ether.

Open problem. The quantitative connection between the Landau damping rate (5.121) and Boyer's absorption/emission rates (Section 6.2.3) has not been derived from first principles. Such a derivation would require computing the multi-particle SED equilibrium distribution $f_{0,s}(v)$ from the ZPF spectrum — extending Boyer's single-particle result to the many-body case — and showing that Landau damping returns perturbations to this ZPF-maintained distribution. This is a well-posed problem within the SED framework but requires the nonlinear multi-particle SED techniques that are the subject of ongoing research [150, 90]. We flag it as an open problem for the ether programme.

---

5.7 Summary and Assessment

5.7.1 Key Results

We collect the results of this section:

1. Electromagnetic dielectric equation derived (Theorem 5.1). The linearised collective response of free charges in the ether, governed by the SED equation of motion (5.16), modifies the ether's constitutive relation from $\mathbf{D} = \varepsilon_0\mathbf{E}$ to $\mathbf{D} = \varepsilon_0\,\boldsymbol{\varepsilon}_r(\omega)\cdot\mathbf{E}$ ((5.48)). The full anisotropic Stix tensor ((5.46)–(5.47)) is derived from the ether's ZPF-driven charge dynamics. The ZPF maintains the equilibrium but contributes zero to the linear response ((5.18)).

2. Structural parallel with the gravitational dielectric established (Section 5.4.5, (5.49)). The electromagnetic permittivity $\boldsymbol{\varepsilon}_r(\omega)$ is the transverse-sector analog of the gravitational permittivity $\mu_e(g/a_0)$ of Theorem 4.1. Both are constitutive responses of the same physical medium — the ether — to different perturbation sources (free charges and baryonic mass, respectively).

3. Electromagnetic wave dispersion in the plasma-ether derived ((5.62)). Transverse electromagnetic waves obey $\omega^2 = \omega_p^2 + c^2k^2$, with the plasma frequency $\omega_p = \sqrt{n_ee^2/(\varepsilon_0 m_{\text{el}})}$ as the ether's low-frequency electromagnetic cutoff. Below $\omega_p$, the perturbed ether is opaque; above $\omega_p$, it is transparent with group velocity below $c$.

4. Alfvén–Ether equivalence proved (Theorem 5.2). The shear Alfvén wave dispersion $\omega^2 = k_\parallel^2 v_A^2$ ((5.79)) is formally identical to transverse wave propagation in Young's elastic ether ((2.2)) with effective shear modulus $G_{\text{eff}} = B_0^2/\mu_0$ and density $\rho_{\text{eff}} = \rho$. The derivation was verified by two independent routes: the cold plasma dielectric in the MHD limit, and direct linearisation of the ideal MHD equations.

5. Vlasov equation derived from ether dynamics ((5.95)). The fundamental kinetic equation of collisionless plasma physics follows from the SED equation of motion under the mean-field approximation, with phase-space incompressibility proved from the vanishing of $\epsilon_{iik}$ ((5.94)).

6. Preferred-frame structure identified (Section 5.6.3). The Vlasov equation is most naturally formulated in Newtonian phase space with a preferred frame — the plasma/ether rest frame. The covariant generalization ((5.97)) exists but is substantially more complex. This is the pattern of Theorem 1.1 (LET/SR empirical equivalence) applied to kinetic theory: both formulations make identical predictions, but the preferred-frame version is more natural.

7. Kinetic dielectric function derived ((5.106)). The Vlasov equation yields a warm-plasma dielectric $\varepsilon_L(\omega, k)$ that generalizes the cold-plasma result ((5.29), recovered in (5.109); the warm-plasma Bohm–Gross relation is derived as (5.109j)) and includes the resonant wave–particle interaction responsible for Landau damping.

8. Landau damping derived and connected to Boyer's energy balance ((5.121), Section 5.6.5). The collisionless damping rate $\gamma_L = -\sqrt{\pi/8}\;\omega_p/(k\lambda_D)^3\exp(-1/(2k^2\lambda_D^2))$ follows from the Landau contour prescription applied to the kinetic dielectric. Landau damping is identified as the kinetic-level manifestation of the same resonant energy exchange — between electromagnetic excitations of the ether and charged matter — that maintains Boyer's quantum ground states (Theorem 6.1).

5.7.2 Open Problems

9. The $\omega_p$–$\ell_e$ relationship. The plasma frequency $\omega_p$ and the transverse microstructure scale $\ell_e$ are both electromagnetic-sector quantities, but their connection requires a theory of the ether's transverse dynamics that has not been developed. A complete theory would determine both from the same condensate microphysics.

10. Nonlinear electromagnetic response. The derivations of this section are linearised. Nonlinear plasma phenomena — wave breaking, parametric instabilities, magnetic reconnection — require the nonlinear ether dynamics. Whether the ether's electromagnetic response becomes nonlinear at extreme field strengths (Schwinger critical field) is an open question with implications for vacuum birefringence.

11. Stix tensor from ether rotational microstructure. The off-diagonal elements of the Stix tensor ((5.47)) are derived from the linearised equations of motion, but the ether-specific interpretation of why a background magnetic field introduces anisotropy — the transverse ether's response to rotational deformation — is not developed at the microphysical level.

12. Multi-particle SED equilibrium from ZPF spectrum. The quantitative connection between the Landau damping rate ((5.121)) and Boyer's absorption/emission rates (Section 6.2) has not been derived from first principles. This requires computing the many-body SED equilibrium distribution from the ZPF spectrum — extending Boyer's single-particle result.

13. The EM cutoff problem (partial resolution). Section 6.6.3 identified the relationship between the transverse EM cutoff and the longitudinal phonon cutoff as an open question. This section provides a partial resolution: the plasma frequency $\omega_p$ serves as the ether's low-frequency EM cutoff in charge-dense regions ((5.62)), while the transverse microstructure scale $\ell_e$ provides the high-frequency cutoff (Section 3.8). The intermediate structure — how the ether's transverse modes interpolate between these bounds — remains open.

5.7.3 The Structural Argument

The gravitational sector (Sections 3–4) established that the ether's longitudinal dynamics reproduce the empirical content of Schwarzschild gravity: the acoustic metric gives the PG metric exactly (Theorem 3.2), the superfluid equation of state gives MOND (Theorem 4.1), and the phonon ZPF gives dark energy with $w = -1$ (Theorem 4.2). The quantum sector (Sections 6–8) established that the ether's fluctuations reproduce quantum mechanics: the ZPF gives ground states (Theorem 6.1), Nelson dynamics gives the Schrödinger equation (Theorem 7.1), and the osmotic coupling gives Bell violation (Theorem 8.5).

This section fills the gap between these two pillars. It shows that the ether's transverse constitutive response — the electromagnetic analog of the gravitational dielectric — is not a theoretical construction that must be invented. It is already the working framework of plasma physics. Every plasma physicist who writes down the dielectric tensor $\boldsymbol{\varepsilon}_r(\omega)$, solves the Vlasov equation in the plasma rest frame, or invokes magnetic tension to explain Alfvén waves is performing ether physics — applying the constitutive relations of a physical medium with a preferred frame, material properties ($\varepsilon_0$, $\mu_0$, $Z_0$), and a wave spectrum determined by those properties.

The specific contributions to the ether programme are:

(a) Completion of the constitutive description. The gravitational sector had $\mu_e(g/a_0)$; the electromagnetic sector now has $\boldsymbol{\varepsilon}_r(\omega, \mathbf{k})$. Both are derived, not postulated. The ether's complete constitutive response — longitudinal (gravity + dark sector) and transverse (electromagnetism + plasma) — is now specified.

(b) Partial resolution of the EM cutoff problem. The plasma frequency identifies the ether's low-frequency EM cutoff in charge-dense regions, complementing the high-frequency cutoff from the transverse microstructure.

(c) Identification of medium-language physics in a mainstream discipline. Plasma physics provides a proof-of-concept that preferred-frame, medium-based mathematical language is physically coherent and operationally essential in at least one domain of modern physics. The philosophical argument against the ether — that medium-based physics was superseded by field-based physics — is contradicted by the daily practice of every plasma laboratory and astrophysical simulation in the world.

(d) Connection between SED and plasma kinetic theory. The Boyer–Landau structural correspondence ((5.122)–(5.123)) connects two previously disconnected research programmes: the SED analysis of quantum ground states and the kinetic theory of wave–particle resonance in plasmas. Both are instances of resonant energy transfer between electromagnetic modes of the ether and charged matter.

PART IV: QUANTUM ETHER