The Spectral Geometry of Plasma Confinement: A Davis Field Equations Framework for Fusion Stability, Transport Bottlenecks, and Cross-Domain Universality

This repository contains the complete manuscript, computational code (12 experiments), figures, and supporting materials for The Spectral Geometry of Plasma Confinement, which applies the Davis Field Equations framework to magnetic and inertial fusion — the fourth physical domain (after Yang-Mills gauge theory, Navier-Stokes fluid mechanics, and Bose-Einstein condensate superfluidity) where the framework's central law, C = τ/K, is computationally validated.

The Problem

Fusion confinement has no single geometric invariant. Kink modes, tearing modes, ballooning modes, and edge-localized modes are each analyzed with separate mathematical machinery, separate stability criteria, and separate codes. The field has been solving the same problem — is this magnetic bottle good enough? — one instability at a time for seventy years.

What This Paper Does

We identify the spectral gap λ₁ of the MHD stability operator as that missing invariant and validate it through twelve computational experiments with zero free parameters:

The Cheeger Effect (Experiments 1–2). The Suydam local stability criterion is violated at 100% of the plasma radius — every single flux surface fails the local test — yet the global ballooning eigenvalue is positive. The plasma is stable. The toroidal geometry provides 15.6× more stability than a cylinder through holonomy coupling invisible to any surface-by-surface analysis. This is the spectral-geometric generalization of magnetic shear stabilization, captured by the Cheeger inequality λ₁ ≥ h²/4.

The Non-Decoupling Theorem, Quantified (Experiment 3). Sweeping the inverse aspect ratio ε = a/R₀ from 0.15 to 0.70: β_{N,crit} increases monotonically (Pearson r = 0.92, p = 3.9 × 10⁻⁴). Stabilization factor from 13× to 30×. This is the Yang-Mills action S_YM ∝ ε made measurable in a tokamak: more toroidal curvature → more holonomy → more stability. It matches the known experimental trend — spherical tokamaks (MAST, NSTX) achieve β_N > 5, far above conventional aspect ratios.

The Cheeger Cut Eats the Plasma from the Outside In (Experiment 4). The stability boundary migrates inward from x = 0.80 (the pedestal) to x = 0.27 (deep core) as β_N increases. It crosses the q = 2 surface — the most operationally dangerous surface in real tokamaks, site of the (2,1) tearing mode and disruption precursors — at β_N ≈ 1.2. The transition width captures H-mode (sharp barrier, Δx = 0) vs L-mode (diffuse, Δx = 0.35) phenomenology.

The Varadhan Redemption (Experiment 5). The first attempt at a Varadhan heat kernel test (Experiment 1v1) failed catastrophically: ρ = −0.53 on a hand-tuned graph that built its own bottleneck. Experiment 5 constructs the transport operator from actual ballooning physics — diffusivity derived from λ_min(x) with zero free parameters — and achieves ρ = −0.977 in the stable regime, approaching the Shidoku benchmark of −1.000. The Riemannian manifold structure is real, not imposed.

One Mode, Spectrally Isolated (Experiment 6). At every unstable surface, at every β_N: exactly one eigenvalue is negative. Not two. Not a cascade. One. And the gap between it and the next eigenvalue grows with β_N — the unstable mode decouples from the stable spectrum. The eigenvalue spacing ratio ⟨r⟩ = 0.121, far below Poisson (0.386): the MHD operator is more ordered than an integrable system. The Trichotomy transition is a single-eigenvalue crossing with topological protection from above.

NIF's Geometric Deficit (Experiment 9). The Branch IX exact-regime completion condition Σ codim(Cᵢ) ≥ dim(M) applied to inertial confinement: seven dangerous Rayleigh-Taylor modes, four independent beam constraints, codimension deficit = 3. Robust across every parameter scan. NIF achieved ignition not by closing this geometric deficit but by reducing dim(M) through target quality — the codimension framework explains why it took decades.

The Davis Law Applied to Ballooning: C = s/α (Experiment 12 — the central result). The SUPERFLUID framework identifies magnetic shear s as the tolerance τ and normalized pressure gradient α as the curvature K. The critical ratio C_crit = s/α_crit is:

The s-α stability boundary is a curve, not a line through the origin. At s = 0.2, C = 0.56; at s = 2.5, C = 1.66. The Davis framework predicts that this curve asymptotes to a ray at moderate shear, and the data confirm it with 2–4% precision. This asymptotic constancy is a result, not an input to the Connor-Hastie-Taylor equation. The cross-geometry hierarchy (circular < neg-tri < elongated < D-shaped) mirrors the BEC obstacle-type hierarchy (β_hard < β_soft), and the instability window (1st stability → ballooning → 2nd stability) maps exactly onto the BEC Geometric Trichotomy (subcritical → vortex nucleation → sonic stabilization). The topological correspondence π₁(S¹) = ℤ → π₁(T²) = ℤ × ℤ connects vortex winding numbers to rational-surface labels through the same homotopy classification.

What Failed (and Why It Matters)

Two documented failures are retained in the paper as methodological lessons:

• Experiment 1v1: A hand-crafted MHD energy weight with five free parameters. The resonance penalty covered 74% of q-space — we built the bottleneck we claimed to discover. The Varadhan test caught it: ρ = −0.53. Lesson: no parallel lines in curved space.

• Experiment 3v2: The |∇r|² geometric variation was included in the bending coefficient f_M. This reduced stabilization by 4.4× at θ ~ π/2 — a coordinate artifact that made shaping destabilize the plasma. Lesson: the ballooning angle θ is defined such that f = 1 + Λ² IS the bending energy; multiplying by geometric factors breaks this.

Two honest null results are also reported: Experiment 7 (helicity does not predict stability in its simple integral form — profile shape matters more than total twist) and Experiment 10 (Ollivier-Ricci curvature is 66% positive on the equilibrium graph but the Lichnerowicz bound is technically vacuous).

Testable Predictions

1. ITER C_crit Test: For ITER-like D-shaped equilibria, the marginal ballooning boundary satisfies s/α_crit = 3.93 ± 0.15 for s ≥ 0.8. Testable with EFIT + MISHKA today.

2. Disruption Precursor: λ₁ → 0 via a single-eigenvalue crossing with spectral isolation. The rate |dλ₁/dt| is a geometry-aware disruption precursor. Testable on DIII-D or ASDEX-Upgrade disruption databases.

3. Non-Decoupling Scaling: β_{N,crit} ∝ ε^(1/2) in the large-aspect-ratio limit. Testable across the existing spherical tokamak dataset (MAST-U, NSTX-U, ST40).

4. C-Surface: C_crit(κ, δ) is a smooth, monotonically increasing function of elongation. The full surface can be mapped by a systematic ballooning scan across (κ, δ) space.

Repository Contents

• plasma_confinement.tex — Complete manuscript (LaTeX source)
• plasma_confinement.pdf — Compiled paper
• experiment1_v2_newcomb.py — Cylindrical Newcomb (Exp 1)
• experiment2_ballooning.py — Toroidal s-α ballooning (Exp 2)
• experiment3_shaped.py — Shaped Miller equilibrium (Exp 3)
• experiment4_cheeger_cut.py — Cheeger cut migration (Exp 4)
• experiment5_varadhan.py — Varadhan heat kernel test (Exp 5)
• experiment6_spectral.py — Full spectral analysis (Exp 6)
• experiment7_helicity.py — Helicity-spectral correlation (Exp 7)
• experiment9_icf.py — ICF codimension counting (Exp 9)
• experiment10_ricci.py — Ollivier-Ricci curvature (Exp 10)
• experiment11_universality.py — BEC-plasma universality (Exp 11)
• experiment12_superfluid_plasma.py — C = τ/K applied to CHT (Exp 12)
• plasma_confinement_outline.md — Detailed outline with all results
• exp*_results/ — Output directories with figures and JSON data

All code runs on standard Python 3.10+ with NumPy, SciPy, and Matplotlib. No institutional MHD codes, no GPU, no cluster access. Every result in this paper was computed on a single laptop.

The Through-Line

Yang-Mills. Navier-Stokes. BEC superfluidity. Plasma confinement. Four domains. One law: C = τ/K. The spectral gap of the stability operator — the smallest eigenvalue of the energy Hessian — governs the phase transition between order and disorder in every system where geometry constrains dynamics. The tolerance-to-curvature ratio is the same object seen through four lenses, and the capacity it measures is the same quantity that determines whether a magnetic bottle holds, whether a fluid stays regular, whether a condensate flows without friction, and whether a gauge field develops a mass gap.

This paper provides twelve experiments testing that claim in the domain where it matters most: the confinement of plasma for fusion energy.

Keywords: plasma confinement, MHD stability, spectral gap, Cheeger constant, ballooning modes, tokamak, Davis Field Equations, SUPERFLUID, magnetic shear, pressure gradient, Non-Decoupling Theorem, Varadhan formula, heat kernel, inertial confinement fusion, BEC universality, geometric phase transition, spectral geometry, fusion energy