Powered by OpenAIRE graph
Found an issue? Give us feedback
ZENODOarrow_drop_down
ZENODO
Preprint . 2026
License: CC BY
Data sources: Datacite
ZENODO
Preprint . 2026
License: CC BY
Data sources: Datacite
ZENODO
Preprint . 2026
License: CC BY
Data sources: Datacite
ZENODO
Preprint . 2026
License: CC BY
Data sources: Datacite
ZENODO
Preprint . 2026
License: CC BY
Data sources: Datacite
versions View all 5 versions
addClaim

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

Authors: Davis, Bee Rosa;

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

Abstract

The Troyon Limit Is 1/(2π) The Troyon beta limit — the maximum pressure a tokamak can sustain before the external kink mode destroys confinement — has governed fusion reactor design for 40 years as an empirical fit. We derive it from first principles. β_N,crit = l_i / (2π) in natural units. The conventional coefficient of 4 decomposes exactly: 4 = 8π × 1/(2π) where 8π is unit conversion and 1/(2π) is the inverse circumference of S¹ — a topological constant. The factor comes from the double cover RP¹ → S¹: the n=1 kink lives on the image of this map, which sends circumference π to 2π. The derivation uses the δW energy principle across 20 current profiles at each wall position, with zero free parameters: C × 2π → 1.000 as the wall recedes to infinity CoV = 2.4% across current profiles Independent of aspect ratio First derivation of the Troyon coefficient from topology This is not a fit. It is not an approximation. The number that has constrained every tokamak built since 1984 is the inverse circumference of a circle. Summary This result emerges from the Davis Field Equations framework applied to plasma confinement. Eighteen experiments, two proved theorems, and zero free parameters validate the framework across four levels of MHD physics (ideal through resistive), two confinement approaches, a cross-domain universality test, a direct application of the SUPERFLUID C = τ/K law to ballooning stability, the Troyon derivation above, and resistive extensions covering neoclassical tearing modes, vertical displacement events, and resistive wall modes. Core Results The Troyon Limit Is 1/(2π) — The external kink δW energy principle gives β_N,crit / l_i = 1/(2π) in natural units (CoV = 2.4%, 20 current profiles, zero free parameters). The conventional Troyon coefficient of 4 decomposes as 4 = 8π × 1/(2π) — unit conversion times the inverse circumference of S¹. C × 2π → 1.000 as the wall recedes to infinity. The kink stability constant is a topological invariant, independent of aspect ratio. Troyon Reframing — Self-consistent Grad-Shafranov equilibria remain ballooning-stable past the Troyon limit (β_N = 3.0 > 2.8), identifying Troyon as an external kink mode boundary, not a ballooning boundary. Davis Law Universality — C = s/α_crit is universal within each tokamak geometry class (CoV < 4% at s ≥ 0.8), with a four-class hierarchy mirroring the BEC obstacle-type structure. Plateau Maximum Theorem (Proved) — C(s) has a unique maximum; the tangent condition α'_crit(s_peak) = α_crit(s_peak)/s_peak is satisfied to 0.0% precision across all four geometry classes. The plateau [1.3, 3.0] coincides with the tokamak operating range. The Cheeger Effect — Global stability persisting where local criteria fail at 100% of the plasma radius. The Suydam criterion is violated everywhere, yet the global ballooning eigenvalue is positive. Toroidal Stabilization — 15.6× stabilization from toroidal topology with a quantitative Non-Decoupling prediction (Pearson r = 0.92, p = 3.9 × 10⁻⁴). Cheeger Cut Migration — The Cheeger cut migrates inward from x = 0.80 to 0.27, crossing the q = 2 surface at β_N ≈ 1.2. Varadhan Heat Kernel — ρ = −0.977 in the stable regime, redeeming a failed earlier attempt (ρ = −0.53). Spectral Isolation — Single-mode onset with sub-Poisson spacing statistics (⟨r⟩ = 0.121, cf. Poisson 0.386). ICF Codimension Deficit — Deficit of 3 (7 modes, 4 constraints), explaining NIF's historical difficulty. BEC–Plasma Universality — Spectral capacity agreement within 2.4% for shaped toroidal geometry. Rayleigh Quotient Decomposition (Proved) — C = τ/K = B(κ,δ)/A(κ,δ) is constant because the metric is quadratic in τ and the potential is linear in K. Verified to 0.4–3.6% across all four geometry classes. Neoclassical Tearing Mode — Island growth recast as Cheeger bottleneck collapse via the modified Rutherford equation with bootstrap drive and ECCD stabilization. Vertical Displacement Event — Decay-index criterion as λ₁ → 0 in the vertical eigenspace; controllability window from growth rate vs. feedback bandwidth. Resistive Wall Mode — Bondeson-Ward dispersion relation in the resistive wall window β_N,no-wall < β_N < β_N,wall; margin as Dirichlet-to-Neumann boundary transition. Cross-Geometry Hierarchy Geometry C_max s_peak Plateau Range Circular 1.682 1.78 [1.3, 2.8] Elongated 3.877 2.03 [1.5, 2.9] D-shaped 4.072 2.12 [1.6, 3.0] Negative triangularity 3.714 1.93 [1.4, 2.8] Kink Stability — Wall Position Scan r_w/a C = β_N,crit / l_i C × 2π CoV 1.05 0.041 0.256 74.8% (wall too close) 1.15 0.142 0.893 7.5% 1.50 0.159 1.000 3.2% 2.00 0.162 1.018 2.6% ∞ (no wall) 0.163 1.025 2.4% Instability Hierarchy β_N Boundary ≈ 0.65 Local ballooning (Experiment 3) ≈ 2.8 Troyon limit = l_i/(2π) in natural units (external kink, Experiment 18) ≥ 3.0 Self-consistent ideal ballooning limit (Experiment 13) Two documented failures and two null results are retained as lessons, reinforcing that spectral-geometric analysis must derive from first-principles physics. Resistive Extensions β_N regime Mode Spectral signature Any NTM (Exp 19) Cheeger bottleneck collapse at rational surface Near VDE threshold VDE (Exp 20) λ₁ → 0 in vertical eigenspace β_N,no-wall < β_N < β_N,wall RWM (Exp 21) Dirichlet-to-Neumann transition at resistive wall Experiments # Method Result Key Finding 18 δW energy principle ✓✓ β_N,crit / l_i = 1/(2π), CoV = 2.4% 13 FreeGS self-consistent ✓ Stable past Troyon; kink is the limit 12 C = τ/K ✓✓ C_crit universal per class (CoV < 4%) Thm Plateau Maximum Proved Tangent condition at 0.0%; C_max = invariant 1v1 Hand-tuned graph Fail ρ = −0.53; graph built its own bottleneck 1v2 Newcomb (Euler) ✓ Cheeger effect at 7% Suydam violation 2 CHT ballooning ✓ 15.6× stabilization; 100% Suydam-stable 3v2 |∇r|² metric Fail Coordinate artifact; shaping destabilized 3v3 Miller metric ✓ 20× over cylinder; ε scan r = 0.92 4 Cheeger cut ✓ x_cut: 0.80 → 0.27; q = 2 at β_N = 1.2 5 Varadhan kernel ✓ ρ = −0.977 (stable), −0.932 (mixed) 6 Full spectrum ✓ Single-mode onset; ⟨r⟩ = 0.121 7 Helicity Mixed r = 0.63; profile shape mediates 9 ICF codimension ✓ Deficit = 3 (7 modes, 4 constraints) 10 Ollivier-Ricci Mixed 66% concordance; Lichnerowicz vacuous 11 BEC universality Caveat Toroidal C_norm ≈ 1.5; cylinder 40× off 18b Rayleigh quotient Proved C = B/A forced by operator structure 19 NTM (Rutherford) ✓ Island growth as Cheeger bottleneck collapse 20 VDE (decay index) ✓ λ₁ → 0 in vertical eigenspace 21 RWM (Bondeson-Ward) ✓ Margin as Dirichlet-to-Neumann transition Repository Contents Paper theory/plasma_confinement.tex — Full 41-page paper (LaTeX source) theory/plasma_confinement.pdf — Compiled PDF Experiment Scripts (Python) experiment1_v2_newcomb.py — Cylindrical Newcomb / Cheeger effect experiment2_ballooning.py — CHT toroidal ballooning experiment3_shaped.py — Shaped Miller equilibrium experiment4_cheeger_cut.py — Cheeger cut migration with β_N experiment5_varadhan.py — Varadhan heat kernel test experiment6_spectral.py — Full spectral analysis experiment7_helicity.py — Helicity correlation experiment9_icf.py — ICF codimension counting experiment10_ricci.py — Ollivier-Ricci curvature experiment11_universality.py — BEC–plasma universality experiment12_superfluid_plasma.py — SUPERFLUID C = τ/K applied to CHT plateau_maximum_theorem.py — Plateau Maximum Theorem verification experiment13_freegs.py — FreeGS Grad-Shafranov / Troyon limit experiment13_troyon.py — Troyon limit analysis (built-in GS solver) experiment17_toroidal_kink.py — Toroidal coupled kink / Newcomb operator (Experiment 17) experiment18_dcon_kink.py — DCON energy principle kink stability (Experiment 18) experiment18b_wkb_rayleigh.py — WKB Rayleigh quotient proof (Experiment 18b) dcon_energy.py — δW energy principle: kink stability & Troyon derivation (Experiment 18) troyon_kink_experiment_14_16.py — Early symmetric kink derivation attempts (Experiments 14–16) Results exp3_results/ through exp13_results/ — Figures and JSON data for each experiment theorem_results/ — Plateau Maximum Theorem verification data and figures CHIHIRO v2.0 — Real-Time Plasma Diagnostic Live at chihiro.sh — CHIHIRO is the Confinement Hierarchy via Invariant Harmonic Instability Recognition Operator: a real-time MHD stability diagnostic built on the spectral-geometric framework of this paper. Single Rust binary, sub-6ms per diagnostic, native G-EQDSK input, zero free parameters. Version 2.0 adds a DCON external kink solver (Riccati/Δ' method) that computationally validates the topological Troyon constant $C_{\rm kink} \approx 1/(2\pi)$ derived in Experiment 18, a unified four-channel stability report (ballooning + kink + ELM-free + density limit), a 100 Hz real-time monitor with time-to-boundary alarm, and resistive extensions for NTM, VDE, and RWM instabilities. 665 tests, 0 failures. Designed for the tokamak control room. Dependencies Python 3.12+ NumPy, SciPy, Matplotlib FreeGS (pip install freegs) for Experiment 13 Related Publications Davis, B. R. (2025). The Field Equations of Semantic Coherence. Zenodo. DOI: 10.5281/zenodo.14784553 Davis, B. R. (2025). The Incompressibility of Topological Charge: A Spectral Proof of the Yang-Mills Mass Gap. Zenodo. DOI: 10.5281/zenodo.17846521 Davis, B. R. (2026). Holonomy-First Navier-Stokes Regularity. Zenodo. DOI: 10.5281/zenodo.18216597 Davis, B. R. (2026). The Davis-Landau Sonic Onset Law: Universality of the Critical Mach Number for Vortex Nucleation. Zenodo. DOI: 10.5281/zenodo.18369500 Davis, B. R. (2025). The Geometry of Sameness: SUPERFLUID. Amazon. Patent Notice Aspects of the CHIHIRO diagnostic system described herein are the subject of U.S. Provisional Patent Application No. 63/999,346. The scientific results, experimental data, and theoretical framework are freely available under CC BY 4.0. Keywords plasma confinement, magnetohydrodynamic stability, spectral gap, Cheeger constant, ballooning modes, external kink mode, tokamak, Davis Field Equations, Connor-Hastie-Taylor equation, s-alpha diagram, Plateau Maximum Theorem, Troyon limit, energy principle, delta-W, Grad-Shafranov, BEC universality, geometric phase transition, double cover, fusion energy, CHIHIRO License Creative Commons Attribution 4.0 International (CC BY 4.0) applies to the paper, experimental data, source code, and theoretical framework. The CHIHIRO diagnostic system is subject to U.S. Provisional Patent Application No. 63/999,346; commercial use of the patented methods requires a separate license.

  • BIP!
    Impact byBIP!
    selected citations
    These citations are derived from selected sources.
    This is an alternative to the "Influence" indicator, which also reflects the overall/total impact of an article in the research community at large, based on the underlying citation network (diachronically).
    0
    popularity
    This indicator reflects the "current" impact/attention (the "hype") of an article in the research community at large, based on the underlying citation network.
    Average
    influence
    This indicator reflects the overall/total impact of an article in the research community at large, based on the underlying citation network (diachronically).
    Average
    impulse
    This indicator reflects the initial momentum of an article directly after its publication, based on the underlying citation network.
    Average
Powered by OpenAIRE graph
Found an issue? Give us feedback
selected citations
These citations are derived from selected sources.
This is an alternative to the "Influence" indicator, which also reflects the overall/total impact of an article in the research community at large, based on the underlying citation network (diachronically).
BIP!Citations provided by BIP!
popularity
This indicator reflects the "current" impact/attention (the "hype") of an article in the research community at large, based on the underlying citation network.
BIP!Popularity provided by BIP!
influence
This indicator reflects the overall/total impact of an article in the research community at large, based on the underlying citation network (diachronically).
BIP!Influence provided by BIP!
impulse
This indicator reflects the initial momentum of an article directly after its publication, based on the underlying citation network.
BIP!Impulse provided by BIP!
0
Average
Average
Average
Upload OA version
Are you the author of this publication? Upload your Open Access version to Zenodo!
It’s fast and easy, just two clicks!