pith. sign in

arxiv: 2604.08878 · v1 · submitted 2026-04-10 · ⚛️ physics.plasm-ph

Statistical equilibrium model for stellarators

Maximilian Ruth , Joshua W. Burby , Wrick Sengupta , Andrew Brown This is my paper

Pith reviewed 2026-05-10 17:15 UTC · model grok-4.3

classification ⚛️ physics.plasm-ph
keywords stellarator equilibriumMHDplasma fluctuationsvariational principlesmooth solutionsstatistical averagingGrad-Shafranovresonant surfaces
0
0 comments X

The pith

A statistical model of rapid ergodic magnetic fluctuations yields smooth equilibrium solutions for stellarators by averaging the force balance.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

Standard MHD equilibria in three-dimensional toroidal domains without symmetry produce singular plasma currents on resonant surfaces that violate length-scale separation and hinder numerical convergence. This paper starts from the assumption that the magnetic field fluctuates ergodically and rapidly relative to the MHD timescale, then averages the resulting force to obtain a variational principle for the statistical mean magnetic field that depends on fluctuation variance. Asymptotics, numerical simulations, and a Grad-Shafranov-type argument show that the new principle admits smooth solutions when the fluctuation statistics are chosen to change the standard model as little as possible. Physically, the singularities are smoothed over a length scale fixed by the amplitude of the fluctuations. The result keeps the equilibrium modeling framework largely intact while removing the singular-current pathology.

Core claim

By treating the plasma magnetic field as ergodically and rapidly fluctuating relative to the MHD time scale and averaging the force, one obtains a closed variational equilibrium problem for the statistical mean magnetic field whose solutions depend on the fluctuation variance; asymptotics, simulations, and a Grad-Shafranov-type argument then establish that this principle supports smooth solutions for specific fluctuation statistics chosen to minimally modify the standard equilibrium modeling paradigm, thereby smoothing singular current sheets with a length scale set by the magnetic field fluctuations.

What carries the argument

Averaged force balance obtained from ergodic rapid fluctuations in the magnetic field, producing a variance-dependent variational principle for the mean field.

If this is right

  • Smooth solutions exist in fully three-dimensional toroidal domains without imposed symmetry.
  • Singular current sheets are replaced by continuous current distributions whose width is controlled by the fluctuation amplitude.
  • Numerical approximations of the equilibrium converge under grid refinement.
  • The modeling change is confined to the choice of fluctuation statistics, leaving the overall variational structure close to conventional MHD.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The fluctuation variance could be calibrated against turbulence or transport calculations to make the equilibrium self-consistent with observed fluctuation levels.
  • The same averaging procedure might be applied to other singular features in plasma models, such as current sheets in reconnection or edge plasmas.
  • If the required fluctuation statistics prove realizable in experiment, the model would allow stellarator design codes to avoid artificial symmetry assumptions without introducing new free parameters.

Load-bearing premise

The plasma magnetic field must fluctuate ergodically and rapidly compared with the MHD time scale so that the force can be averaged to a closed variational problem.

What would settle it

Direct numerical comparison of the averaged variational solutions against standard MHD equilibria, checking whether singularities on resonant surfaces disappear as fluctuation variance is increased while keeping the mean field fixed.

Figures

Figures reproduced from arXiv: 2604.08878 by Andrew Brown, Joshua W. Burby, Maximilian Ruth, Wrick Sengupta.

Figure 1
Figure 1. Figure 1: A schematic of the considered geometry. On the left, the fluid domain Q is bounded by r bottom and r top in dark blue and red respectively. The back-to-labels map G maps Q to the reference domain Q0, with the boundary condition that the top and bottom fluid boundaries must map to the corresponding top and bottom boundaries in the reference configuration. magnetic confinement system. The second parameter λ … view at source ↗
Figure 2
Figure 2. Figure 2: Comparison of a direct numerical solution of Eq. (20) (black line) and the asymptotic solution (25) (red dashed line). for a resonance at the surface v(rs) = 0.5 with a rotational transform ι = 1/2. The leading-order solution is found by directly minimizing Eq. (10) on the undeformed domain, where θ = ζ = 0 and v0 is represented by a degree-200 Legendre polynomial. From the leading-order solution, we find … view at source ↗
Figure 3
Figure 3. Figure 3: A plot demonstrating the convergence of the numerical method for ϵ = 10−3 . The top row shows convergence of the force-balance residual (26) as a function of λ and of the (a) radial, (b) toroidal, and (c) poloidal resolutions. The bottom row shows the self-convergence (27) of the solution of the variational equilibrium problem in the (d) radial, (e) toroidal, and (f) poloidal resolutions. 4.3 Results This … view at source ↗
Figure 4
Figure 4. Figure 4: (a) Number of L-BFGS iterations to the solution and (b) the time to solution, both as a function of the fluctuation parameter λ and perturbation magnitude ϵ [PITH_FULL_IMAGE:figures/full_fig_p014_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Current sheets plotted on the y = 0 plane. The fluctuation parameter is scanned from left to right, with λ = 0.005 (a,d), λ = 0.01 (b,e), and λ = 0.02 (c,f). Two values of the perturbation magnitude, ϵ = 0.1 (a,b,c) and ϵ = 0.025 (d,e,f), are plotted. Lines showing the predicted current sheet widths are plotted at v = v0(rs) ± λv′ 0 (rs)Lmn(rs) for the resonant values of v0(rs) ∈ {0.35, 0.65}. convergence … view at source ↗
read the original abstract

In three dimensional toroidal domains without symmetry, the standard magnetohydrodynamic (MHD) equilibrium model used for magnetic confinement fusion does not generally support smooth solutions. Instead, solutions have singular plasma currents on resonant magnetic surfaces that violate the MHD assumption of length-scale separation, further leading to the non- or slow convergence of numerical approximations under refinement. In this work, we present an improved equilibrium principle derived from a statistical model for plasma fluctuations. Instead of being static, we assume that the plasma magnetic field is ergodically and rapidly fluctuating relative to the MHD time scale. By averaging the resulting force, we derive a variational equilibrium problem for the statistical mean magnetic field which depends on fluctuation variance. Then, through asymptotics, numerical simulations, and a Grad-Shafranov type argument, we show that the variational principle supports smooth solutions for specific fluctuation statistics chosen to minimally modify the standard equilibrium modeling paradigm. Physically, this model smooths singular current sheets with a length scale determined by the magnetic field fluctuations.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 2 minor

Summary. The paper proposes a statistical equilibrium model for stellarators in non-symmetric 3D toroidal domains. Assuming the magnetic field fluctuates ergodically and rapidly relative to the MHD timescale, the Lorentz force is averaged to derive a variational principle for the statistical mean field that depends on fluctuation variance. Through asymptotics, numerical simulations, and a Grad-Shafranov-type reduction, the authors show that this principle admits smooth solutions for specific fluctuation statistics chosen to minimally modify the standard MHD paradigm, thereby smoothing singular current sheets on a scale set by the fluctuations.

Significance. If the claims are substantiated, the work provides a physically motivated approach to regularizing 3D MHD equilibria by incorporating statistical fluctuations, potentially resolving convergence issues with singular currents in stellarator modeling. Credit is due for combining asymptotic analysis, simulations, and a reduction argument to support the existence of smooth solutions. The model introduces fluctuation variance as a free parameter controlling the smoothing length, which limits predictive power unless independently constrained.

major comments (3)
  1. [Abstract and derivation of the averaged variational principle] The fluctuation variance is presented as an input parameter (Abstract) that determines the smoothing scale of singular currents. The abstract specifies that fluctuation statistics are 'chosen to minimally modify the standard equilibrium modeling paradigm,' indicating post-hoc selection. This creates a potential circularity: the variational principle may reduce to a regularization scheme whose smoothness depends on tuning the variance rather than on independent physical input.
  2. [Asymptotics and Grad-Shafranov reduction] The central claim that the variational principle supports smooth solutions rests on the asymptotic analysis and Grad-Shafranov-type argument. Without explicit error estimates for the averaging procedure or verification that the statistical assumptions are preserved under the reduction, the robustness of the smoothness result cannot be fully assessed. The weakest assumption (ergodic rapid fluctuations relative to MHD timescale) requires stronger justification for stellarator applicability.
  3. [Numerical simulations section] Numerical simulations are invoked to support smooth solutions, but the choice of fluctuation statistics in these runs appears tailored to the desired outcome. A sensitivity study showing how results change when statistics deviate from the 'minimal modification' choice would be needed to demonstrate that the smoothing is not an artifact of the specific selection.
minor comments (2)
  1. [Notation and definitions] Clarify the precise mathematical definition of fluctuation variance and its statistical properties (e.g., how it enters the averaged force) at first introduction to avoid ambiguity in later sections.
  2. [Introduction] Add references to existing literature on statistical averaging in MHD or fluctuation-based equilibrium models to better situate the novelty of the approach.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. We address each major comment point by point below, indicating where revisions will be made.

read point-by-point responses

  1. Referee: [Abstract and derivation of the averaged variational principle] The fluctuation variance is presented as an input parameter (Abstract) that determines the smoothing scale of singular currents. The abstract specifies that fluctuation statistics are 'chosen to minimally modify the standard equilibrium modeling paradigm,' indicating post-hoc selection. This creates a potential circularity: the variational principle may reduce to a regularization scheme whose smoothness depends on tuning the variance rather than on independent physical input.

    Authors: The fluctuation variance is introduced as a physically motivated parameter that quantifies the amplitude of rapid, ergodic magnetic fluctuations relative to the MHD timescale, consistent with observed plasma turbulence in stellarators. The specific choice of statistics is made to ensure that the model recovers the standard MHD equilibrium exactly in the zero-variance limit, thereby isolating the regularization effect due to fluctuations without introducing extraneous physics. This is not post-hoc tuning for smoothness but a deliberate minimal-extension approach. We will revise the abstract and add a clarifying paragraph in the introduction to better articulate the physical motivation and experimental grounding for the fluctuation model. revision: partial

  2. Referee: [Asymptotics and Grad-Shafranov reduction] The central claim that the variational principle supports smooth solutions rests on the asymptotic analysis and Grad-Shafranov-type argument. Without explicit error estimates for the averaging procedure or verification that the statistical assumptions are preserved under the reduction, the robustness of the smoothness result cannot be fully assessed. The weakest assumption (ergodic rapid fluctuations relative to MHD timescale) requires stronger justification for stellarator applicability.

    Authors: The asymptotic analysis is formal, relying on scale separation between rapid fluctuations and the MHD equilibrium. We agree that explicit error estimates would add rigor but are technically demanding and beyond the scope of the present work, which prioritizes derivation of the model and demonstration via complementary methods. We will expand the discussion of the ergodic rapid-fluctuation assumption by citing relevant literature on fluctuation timescales in stellarator plasmas. We will also explicitly confirm in the revised text that the statistical assumptions remain consistent under the Grad-Shafranov reduction. revision: partial

  3. Referee: [Numerical simulations section] Numerical simulations are invoked to support smooth solutions, but the choice of fluctuation statistics in these runs appears tailored to the desired outcome. A sensitivity study showing how results change when statistics deviate from the 'minimal modification' choice would be needed to demonstrate that the smoothing is not an artifact of the specific selection.

    Authors: We concur that robustness to the choice of fluctuation statistics should be demonstrated. In the revised manuscript we will add a sensitivity study presenting results for alternative fluctuation statistics that deviate from the minimal-modification case. These will show that the qualitative smoothing of singular currents persists, with the smoothing length scaling according to the variance, confirming that the effect is not an artifact of the specific selection. revision: yes

Circularity Check

0 steps flagged

No significant circularity identified

full rationale

The derivation begins from the explicit physical assumption of ergodic, rapid magnetic fluctuations relative to the MHD timescale. Averaging the Lorentz force then produces a variational principle for the statistical mean field that depends on fluctuation variance as an input parameter. Asymptotics, numerical simulations, and a Grad-Shafranov-type argument are subsequently applied to demonstrate existence of smooth solutions when fluctuation statistics are selected to minimally modify the standard paradigm. None of these steps reduces by construction to its own inputs: the variance remains an independent modeling choice rather than a fitted or self-defined quantity, and the existence proofs are independent analyses rather than tautological re-statements. No self-citations, uniqueness theorems, or renamings of known results are invoked in a load-bearing way. The chain is self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The model rests on the ergodic rapid-fluctuation assumption and the existence of a well-defined average force; fluctuation variance is introduced as a controlling parameter whose value is not derived from first principles.

free parameters (1)
  • fluctuation variance
    Controls the smoothing length scale of current sheets; chosen to minimally modify standard MHD while producing smooth solutions.
axioms (1)
  • domain assumption plasma magnetic field fluctuates ergodically and rapidly relative to MHD time scale
    Invoked to justify averaging the force balance and closing the equations for the mean field.

pith-pipeline@v0.9.0 · 5471 in / 1262 out tokens · 35443 ms · 2026-05-10T17:15:49.945650+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

42 extracted references · 42 canonical work pages

  1. [1]

    Imbert-G´ erard, E

    L.-M. Imbert-G´ erard, E. J. Paul, and A. M. Wright.An Introduction to Stellarators: From Magnetic Fields to Symmetries and Optimization. Other Titles in Applied Mathematics. Society for Industrial and Applied Mathematics, January 2024

    work page 2024

  2. [2]

    S. P. Hirshman and J. C. Whitson. Steepest-descent moment method for three-dimensional magnetohydrodynamic equilibria.The Physics of Fluids, 26(12):3553–3568, December 1983

    work page 1983

  3. [3]

    D. W. Dudt and E. Kolemen. DESC: A stellarator equilibrium solver.Physics of Plasmas, 27(10):102513, 2020

    work page 2020

  4. [4]

    Panici, R

    D. Panici, R. Conlin, D. W. Dudt, K. Unalmis, and E. Kolemen. The DESC stellarator code suite. Part 1. Quick and accurate equilibria computations.Journal of Plasma Physics, 89:955890303, 2023

    work page 2023

  5. [5]

    S. A. Lazerson, J. Loizu, S. Hirshman, and S. R. Hudson. Verification of the ideal magnetohydrodynamic response at rational surfaces in the VMEC code.Physics of Plasmas, 23(1):012507, January 2016

    work page 2016

  6. [6]

    Suzuki, N

    Y. Suzuki, N. Nakajima, K. Watanabe, Y. Nakamura, and T. Hayashi. Development and application of HINT2 to helical system plasmas.Nuclear Fusion, 46(11):L19–L24, November 2006. 21 REFERENCES Burbyet al

    work page 2006

  7. [7]

    S. P. Hirshman, R. Sanchez, and C. R. Cook. SIESTA: A scalable iterative equilibrium solver for toroidal applications.Physics of Plasmas, 18(6):062504, June 2011

    work page 2011

  8. [8]

    Blickhan, J

    T. Blickhan, J. Stratton, and A. A. Kaptanoglu. MRX: A differentiable 3D MHD equilibrium solver without nested flux surfaces, 2025

    work page 2025

  9. [9]

    Reiman and H

    A. Reiman and H. Greenside. Calculation of three-dimensional MHD equilibria with islands and stochastic regions.Computer Physics Communications, 43(1):157–167, December 1986

    work page 1986

  10. [10]

    Drevlak, D

    M. Drevlak, D. Monticello, and A. Reiman. PIES free boundary stellarator equilibria with improved initial conditions.Nuclear Fusion, 45(7):731–740, July 2005

    work page 2005

  11. [11]

    B. F. Kraus and S. R. Hudson. Theory and discretization of ideal magnetohydrodynamic equilibria with fractal pressure profiles.Physics of Plasmas, 24(9):092519, September 2017

    work page 2017

  12. [12]

    S. R. Hudson and N. Nakajima. Pressure, chaotic magnetic fields, and magnetohydrodynamic equilibria.Physics of Plasmas, 17(5):052511, May 2010

    work page 2010

  13. [13]

    S. R. Hudson, R. L. Dewar, G. Dennis, M. J. Hole, M. McGann, G. von Nessi, and S. Lazerson. Computation of multi-region relaxed magnetohydrodynamic equilibria.Physics of Plasmas, 19(11):112502, November 2012

    work page 2012

  14. [14]

    O. P. Bruno and P. Laurence. Existence of three-dimensional toroidal MHD equilibria with nonconstant pressure.Communications on Pure and Applied Mathematics, 49(7):717–764, 1996

    work page 1996

  15. [15]

    I. B. Bernstein, E. A. Frieman, M. D. Kruskal, and R. M. Kulsrud. An energy principle for hydromagnetic stability problems.Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences, 244(1236):17–40, February 1958

    work page 1958

  16. [16]

    H. Grad. Toroidal Containment of a Plasma.The Physics of Fluids, 10(1):137–154, January 1967

    work page 1967

  17. [17]

    Watanabe, N

    K. Watanabe, N. Nakajima, M. Okamoto, Y. Nakamura, and M. Wakatani. Three-dimensional MHD equilibrium in the presence of bootstrap current for the Large Helical Device.Nuclear Fusion, 32(9):1499, September 1992

    work page 1992

  18. [18]

    Landreman, S

    M. Landreman, S. Buller, and M. Drevlak. Optimization of quasi-symmetric stellarators with self-consistent bootstrap current and energetic particle confinement.Physics of Plasmas, 29(8):082501, August 2022

    work page 2022

  19. [19]

    D. D. Holm, J. E. Marsden, and T. S. Ratiu. The Euler–Poincar´ e Equations and Semidirect Products with Applications to Continuum Theories.Advances in Mathematics, 137(1):1–81, July 1998

    work page 1998

  20. [20]

    Grad and H

    H. Grad and H. Rubin. Hydromagnetic Equilibria and Force-Free Fields. InProceedings of the Second United Nations International Conference on the Peaceful Uses of Atomic Energy, volume 31, page 190, Geneva, 1958. United Nations

    work page 1958

  21. [21]

    Some New Variational Properties of Hydromagnetic Equilibria.The Physics of Fluids, 7(8):1283–1292, August 1964

    Harold Grad. Some New Variational Properties of Hydromagnetic Equilibria.The Physics of Fluids, 7(8):1283–1292, August 1964

    work page 1964

  22. [22]

    Rahbarnia, H

    K. Rahbarnia, H. Thomsen, J. Schilling, S. vaz Mendes, M. Endler, R. Kleiber, A. K¨ onies, M. Borchardt, C. Slaby, T. Bluhm, M. Zilker, B. B. Carvalho, and Wendelstein 7-X Team. Alfv´ enic fluctuations measured by in-vessel Mirnov coils at the Wendelstein 7-X stellarator. Plasma Physics and Controlled Fusion, 63(1):015005, November 2020

    work page 2020

  23. [23]

    S. P. Hirshman and H. K. Meier. Optimized Fourier representations for three-dimensional magnetic surfaces.The Physics of Fluids, 28(5):1387–1391, May 1985

    work page 1985

  24. [24]

    J. P. Boyd.Chebyshev and Fourier spectral methods. Dover books on mathematics. Dover Publ, Mineola, NY, second edition, 2001

    work page 2001

  25. [25]

    Weitzner

    H. Weitzner. Ideal magnetohydrodynamic equilibrium in a non-symmetric topological torus. Physics of Plasmas, 21(2):022515, February 2014. 22 Burbyet al

    work page 2014

  26. [26]

    Classics in Mathematics

    Lars H¨ ormander.The Analysis of Linear Partial Differential Operators III: Pseudo-Differential Operators. Classics in Mathematics. Springer Berlin Heidelberg, Berlin, Heidelberg, 2007

    work page 2007

  27. [27]

    Taylor.Partial Differential Equations III: Nonlinear Equations, volume 117 of Applied Mathematical Sciences

    Michael E. Taylor.Partial Differential Equations III: Nonlinear Equations, volume 117 of Applied Mathematical Sciences. Springer International Publishing, Cham, 2023

    work page 2023

  28. [28]

    Jenko, W

    F. Jenko, W. Dorland, M. Kotschenreuther, and B. N. Rogers. Electron temperature gradient driven turbulence.Physics of Plasmas, 7(5):1904–1910, May 2000

    work page 1904

  29. [29]

    G¨ orler, X

    T. G¨ orler, X. Lapillonne, S. Brunner, T. Dannert, F. Jenko, F. Merz, and D. Told. The global version of the gyrokinetic turbulence code GENE.Journal of Computational Physics, 230(18):7053–7071, August 2011

    work page 2011

  30. [30]

    H¨ ofler et al

    K. H¨ ofler et al. Milestone in predicting core plasma turbulence: successful multi-channel validation of the gyrokinetic code GENE.Nature Communications, 16(1):2558, March 2025

    work page 2025

  31. [31]

    Huang, Y

    Y.-M. Huang, Y. Zhou, J. Loizu, S. Hudson, and A. Bhattacharjee. Structure of pressure-gradient-driven current singularity in ideal magnetohydrodynamic equilibrium. Plasma Physics and Controlled Fusion, 65(3):034008, February 2023

    work page 2023

  32. [32]

    Sato and M

    N. Sato and M. Yamada. Nested invariant tori foliating a vector field and its curl: Toward MHD equilibria and steady Euler flows in toroidal domains without continuous Euclidean isometries.Journal of Mathematical Physics, 64(8):081505, August 2023

    work page 2023

  33. [33]

    Betancourt

    O. Betancourt. BETAS, a spectral code for three-dimensional magnetohydrodynamic equilibrium and nonlinear stability calculations.Communications on Pure and Applied Mathematics, 41(5):551–568, 1988

    work page 1988

  34. [34]

    Landreman and E

    M. Landreman and E. Paul. Magnetic Fields with Precise Quasisymmetry for Plasma Confinement.Physical Review Letters, 128(3):035001, January 2022

    work page 2022

  35. [35]

    D. A. Garren and A. H. Boozer. Existence of quasihelically symmetric stellarators.Physics of Fluids B: Plasma Physics, 3(10):2822–2834, October 1991

    work page 1991

  36. [36]

    Rodr´ ıguez and A

    E. Rodr´ ıguez and A. Bhattacharjee. Islands and current singularities in quasisymmetric toroidal plasmas.Physics of Plasmas, 28(9):092506, September 2021

    work page 2021

  37. [37]

    Rodr´ ıguez, W

    E. Rodr´ ıguez, W. Sengupta, and A. Bhattacharjee. Generalized Boozer coordinates: A natural coordinate system for quasisymmetry.Physics of Plasmas, 28(9):092510, September 2021

    work page 2021

  38. [38]

    D. D. Holm. Lagrangian averages, averaged Lagrangians, and the mean effects of fluctuations in fluid dynamics.Chaos: An Interdisciplinary Journal of Nonlinear Science, 12(2):518–530, June 2002

    work page 2002

  39. [39]

    arXiv preprint arXiv:2308.14127 (2023)

    R. Cao and F. Sch¨ afer. Information geometric regularization of the barotropic Euler equation, March 2026. arXiv:2308.14127 [math]

    work page arXiv 2026

  40. [40]

    Barham, B

    W. Barham, B. K. Tran, B. S. Southworth, and F. Sch¨ afer. Hamiltonian Information Geometric Regularization of the Compressible Euler Equations, December 2025. arXiv:2512.13948 [math-ph]

    work page arXiv 2025

  41. [41]

    S´ anchez-Vizuet, M

    T. S´ anchez-Vizuet, M. E. Solano, and A. J. Cerfon. Adaptive Hybridizable Discontinuous Galerkin discretization of the Grad–Shafranov equation by extension from polygonal subdomains.Computer Physics Communications, 255:107239, October 2020

    work page 2020

  42. [42]

    dead spot

    Z. Peng, Q. Tang, and X.-Z. Tang. An Adaptive Discontinuous Petrov–Galerkin Method for the Grad–Shafranov Equation.SIAM Journal on Scientific Computing, 42(5):B1227–B1249, January 2020. A Degeneracy of the ideal MHD potential energy at rational surfaces The MHD potential energy functionalW(G) presents a major difficulty in developing a mathematical theory...

    work page 2020