Pith. sign in

REVIEW 4 major objections 7 minor 57 references

SPECTRE finds 3D plasma equilibria with islands and chaos by minimizing real-space force imbalance on stepped-pressure interfaces.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · grok-4.5

2026-07-30 11:38 UTC pith:ROPE2FH4

load-bearing objection Solid methods upgrade to SPEC that actually works on strongly shaped stellarators; the QI–HINT island claim is softer than the abstract suggests but the verification ladder otherwise holds. the 4 major comments →

arxiv 2607.27135 v1 pith:ROPE2FH4 submitted 2026-07-29 physics.plasm-ph physics.comp-ph

SPECTRE: A robust solver for 3D equilibria with arbitrary topology

E. Balkovic , J. Loizu , E. Lanti , C. Smiet , C. Lazzati , R. Ramasamy , A. Goodman , J. Geiger
show 1 more author
J. P. Graves
This is my paper
classification physics.plasm-ph physics.comp-ph
keywords MRxMHD3D equilibriamagnetic islandsstellaratorforce minimizationfree-boundaryquasi-isodynamicstepped-pressure
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

Three-dimensional fusion plasmas do not always have nested flux surfaces; they can form islands and chaotic field regions that standard nested-surface codes cannot represent. This paper presents SPECTRE, a solver based on multi-region relaxed MHD that divides the plasma into Taylor-relaxed volumes separated by ideal interfaces and drives the jump in total pressure across those interfaces to zero. Unlike its predecessor SPEC, SPECTRE evaluates force in real space on a grid and minimizes the squared residual with a trust-region least-squares scheme, so it does not require an almost-perfect initial guess and does not discard high-order force harmonics. The authors verify the method on vacuum and finite-beta cases with known solutions in both fixed- and free-boundary modes, including a strongly shaped quasi-axisymmetric design and W7-X, then show that a modern finite-beta quasi-isodynamic equilibrium develops a core island whose phase and width match an independent HINT calculation. If the approach holds up, it supplies a practical route to equilibria that include the topology changes actually present in optimized stellarators.

Core claim

SPECTRE computes stepped-pressure MRxMHD equilibria with general magnetic topology by representing the interface force in real space and minimizing its squared residual with a stable trust-region least-squares algorithm. The resulting solver recovers known vacuum and finite-beta fixed- and free-boundary solutions and produces a strongly shaped finite-beta quasi-isodynamic equilibrium whose core island agrees with HINT.

What carries the argument

Real-space force residual minimized by trust-region least squares: the total-pressure jump is sampled on an adjustable (θ,ζ) grid on each ideal interface and the sum of squared samples is driven down without truncating the force spectrum to the interface Fourier resolution.

Load-bearing premise

A modest number of relaxed volumes whose interfaces are spaced by a simple flux rule (here quadratic in toroidal flux) is enough to capture island phase, width, and force balance without systematically biasing the equilibrium.

What would settle it

Run the same finite-beta QI case in free-boundary mode from the coils and compare the core 8/9 island phase, O-point location, and width directly against a free-boundary HINT calculation; a clear mismatch would falsify the claim that the present SPECTRE solution is already reliable.

Watch this falsifier. Get emailed when new claim-graph text bears on it.

If this is right

  • Fixed- and free-boundary MRxMHD equilibria with islands become routinely obtainable for modern optimized stellarators without an expert initial guess.
  • Coil and plasma optimization can target configurations that deliberately retain or control islands rather than assuming perfect nesting.
  • Saturation of tearing and neoclassical tearing modes can be studied as sequences of stepped-pressure equilibria in toroidal geometry.
  • Quantitative free-boundary comparisons with HINT and resistive MHD codes become feasible once interface placement is systematized.

Where Pith is reading between the lines

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

  • Automatic interface placement that avoids low-order rationals and separatrices would remove the last major manual step and make SPECTRE usable inside optimization loops.
  • Because the Beltrami solver is inherited from SPEC, existing analytic force-gradient machinery could later be hybridized with the trust-region loop for faster local refinement once a good basin is found.
  • Agreement on a single core island does not yet guarantee divertor-island fidelity; edge topology is the natural next stress test.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

4 major / 7 minor

Summary. The manuscript presents SPECTRE, a new Multi-Region Relaxed MHD (MRxMHD) equilibrium solver that replaces SPEC's Fourier-truncated force root-finding (Powell/Newton) with a real-space force representation minimized by a trust-region least-squares algorithm, plus interface-intersection checks, map2disc-based initial guesses, and Henneberg angle parametrization. The central claim is that this framework is robust where SPEC was fragile, and that it correctly recovers MRxMHD equilibria with general topology. Verification proceeds through a ladder of tests: axisymmetric axis position vs VMEC (vacuum and β=6%, Nvol and mpol scans), a strongly shaped QA vacuum configuration (Poincaré-error convergence between 1- and 7-volume solutions, axis vs VMEC, cold-start from map2disc), free-boundary W7-X vacuum against a Biot-Savart coil field (nested surfaces, 5/5 divertor island, ι profile), a finite-β free-boundary rotating ellipse against published VMEC/SPEC results (axis agreement at the 10⁻⁴ m level), and finally a fixed-boundary finite-β QI configuration exhibiting an ι=8/9 core island compared qualitatively to HINT. The code is released publicly under MIT license.

Significance. If the results hold, this is a significant practical advance: SPEC's fragility on strongly shaped optimized stellarators is a real and well-known bottleneck, and a robust MRxMHD solver that cold-starts from map2disc guesses and runs free-boundary directly from coils would be broadly useful for optimization and divertor/island studies. The manuscript ships a public MIT-licensed code, uses external verification targets (VMEC, Biot-Savart coil fields, HINT, published SPEC ellipse results) rather than fitted ones, demonstrates Poincaré-error convergence with Fourier resolution, and — commendably — initializes interfaces far from the answer so the verification does not rest on good initial guesses. The QI/HINT island comparison would be a first-of-a-kind benchmark for the community, but it is currently the least supported element of the paper.

major comments (4)
  1. [§3.5] §3.5 (Fig. 12b): the abstract's headline claim — that the core island 'is in agreement with a HINT calculation' — rests on island phase, O-point position, and width. Phase and, to a large extent, O-point radius are substantially inherited by construction, since SPECTRE's fixed boundary is itself an approximate HINT core flux surface and the resonant harmonic content is boundary-driven. The genuinely discriminating observable is the island width, yet it is reported only as 'similar,' with no number. Please quantify the width in both codes (e.g., radial extent of the separatrix at a specified toroidal plane) and report the comparison explicitly.
  2. §3.5: in MRxMHD the interfaces are ideal (B·n̂=0), so the 8/9 island must reside entirely within a single relaxed subvolume, bounding its width by that subvolume's radial extent. With Nvol=8 and quadratic spacing in Ψt, each volume spans roughly Δr ≈ a/8, plausibly comparable to the island width itself. The only sensitivity evidence offered is the assertion that 'perturbing this spacing does not impact the overall equilibrium' — no data, no definition of 'perturbing,' and it is unclear whether the island width was among the checked quantities. If the resonant surface had landed near an interface, the island would have been squeezed or absent. A convergence/sensitivity scan in Nvol (or at least several interface-spacing laws) with the island width reported for each is needed for the HINT comparison to be load-bearing. This is acknowledged as future work, but the abstract's 'in agreement's
  3. [§2.2] §2.2: the convergence criterion for the trust-region least-squares minimization is never stated. Fig. 2 shows ⟨F⟩ decreasing by three orders of magnitude, but the reader is not told what stopping condition, tolerance, or acceptance threshold defines a converged equilibrium, nor how the force-sampling grid density (θi, ζj) is chosen relative to (mpol, ntor). Since under-sampling the real-space force would reintroduce exactly the aliasing problem diagnosed for SPEC in §2.1 (Fig. 1), the paper should state the sampling rule used in each test case and demonstrate insensitivity to it for at least one case.
  4. [§3.3–3.4] §3.3 and §3.4: the free-boundary verification is presented largely visually (Poincaré overlays in Fig. 9, axis positions and ι profiles in Figs. 9c, 10b). Given that the central selling point is robustness and accuracy relative to SPEC, quantitative error metrics for the free-boundary cases (e.g., field-line displacement or B-field difference against the Biot-Savart reference, analogous to the Poincaré-error metric used in §3.2) would substantially strengthen the claim. A direct SPEC-vs-SPECTRE robustness comparison (e.g., success rate from degraded initial guesses) is also promised by the framing ('performs significantly better') but never shown; even a single head-to-head failure case from SPEC's fragility regime would substantiate it.
minor comments (7)
  1. [§2.1] §2.1: 'B2 = (∇ × A)2)' has a stray closing parenthesis.
  2. [§2.2 / Fig. 7] Units of ⟨F⟩ are given as T² in Fig. 7b, but §2.1 quotes the SPEC residual in kPa; a consistent unit (pressure units are physically natural for [p + B²/2μ0]) would ease comparison between sections.
  3. [Fig. 1] Fig. 1: the red box indicating the truncated window is not visible in the reproduced caption description; please verify the figure renders correctly and state the case parameters (resolution, interface shown) in the caption.
  4. [§3.2/3.4] §3.2: wall-clock/CPU costs are reported for QA (9–120 cpu-hours) and the ellipse (~10⁴ cpu-hours); a sentence on how these scale with Nvol and grid sampling, and how they compare to SPEC runtimes, would help readers assess practicality.
  5. [§3.5] §3.5: state the HINT resolution, relaxation parameters, and convergence diagnostics used, so the comparison is reproducible by independent groups; also clarify how the 'approximate HINT core flux surface' was extracted (which surface label, interpolation method).
  6. [§1.3] Eqs. (1.4)–(1.5) and §1.3: 'Equation 1.4' is cited in §1.3 where Eq. (1.2) (Beltrami) is meant — likely a mislabel; please check cross-references.
  7. [§1.1/§4] Consider citing the recent independent MRxMHD developments (e.g., DESC-based stepped-pressure work or other relaxed-equilibrium codes) if applicable, to situate SPECTRE relative to non-SPEC-lineage tools.

Circularity Check

0 steps flagged

No significant circularity: SPECTRE is a numerical solver paper verified against external codes, not a derivation that reduces to its inputs.

full rationale

The paper’s load-bearing claims are that a real-space force residual minimized by trust-region least squares yields robust MRxMHD equilibria, and that those equilibria match known external solutions. Verification targets are independent: VMEC nested-surface solutions (axisymmetry, QA vacuum, rotating ellipse), Biot–Savart coil fields (W7-X free-boundary vacuum), prior SPEC free-boundary ellipse results started from an independent initial guess, and HINT Poincaré structure for the QI case. Force residual ⟨F⟩, Beltrami error, Poincaré error, and axis position are not fitted to reproduce those codes; interfaces are optimized from deliberately poor guesses (map2disc / inward interpolation). Self-citations to SPEC, MRxMHD (Hole/Hudson/Taylor), and analytic force gradients define the model and the baseline solver being improved—they do not supply the success metric or force the reported agreement by construction. Methodological caveats in §3.5 (HINT-extracted fixed boundary, hand-chosen interface flux spacing, unquantified island width, missing Nvol scan) are validation-quality concerns, not circular reductions of prediction to input. No self-definitional loop, fitted-input-as-prediction, uniqueness import, or ansatz smuggling is present in the derivation chain.

Axiom & Free-Parameter Ledger

5 free parameters · 5 axioms · 1 invented entities

The central claim rests on the established MRxMHD variational model and on numerical choices (force sampling, trust-region LS, interface DOFs, Nvol and flux spacing) rather than on new physical entities. Free parameters are discretization and initialization knobs; axioms are standard MHD/MRxMHD domain assumptions plus the modeling choice that stepped-pressure interfaces approximate the target physics.

free parameters (5)
  • Nvol (number of relaxed subvolumes) = case-dependent (4–64 scanned in axisymmetry)
    Chosen per case (e.g. 8 for tokamak/QI, 7 for QA, 16 for ellipse); controls stepped-pressure resolution and island capture.
  • mpol, ntor Fourier resolution = e.g. 6–14 QA; 18 W7-X; 12 ellipse; 12×11 QI
    Truncation of interface geometry and Beltrami basis; scanned for QA and fixed for other cases; limits minimum achievable force.
  • Interface toroidal-flux spacing law = quadratic in Ψt (QI)
    Quadratic in Ψt for QI (linear in approximate minor radius); authors note placement near rationals/separatrices is delicate and only lightly perturbed.
  • Real-space force sampling grid (θ,ζ)
    Adjustable mesh that defines the discrete least-squares objective; spacing trades speed vs force accuracy.
  • Trust-region least-squares hyperparameters
    Step bounds and termination of the SciPy/trust-region scheme affect path and declared convergence of ⟨F⟩.
axioms (5)
  • domain assumption MRxMHD: energy minimized subject to helicity per subvolume yields Beltrami fields ∇×B=μB in volumes and [p+B²/2μ0]=0 on interfaces
    Foundational model (§1.2, Eqs. 1.1–1.3); inherited from Hole/Taylor/Dennis; SPECTRE solves it rather than re-derives it.
  • domain assumption In the large-Nvol limit MRxMHD approaches ideal MHD
    Cited Dennis et al. 2013; used to justify comparison to VMEC when surfaces are nested.
  • ad hoc to paper Truncated real-space least-squares force minimum is an adequate practical equilibrium even if exact continuous force zero is unreachable at finite resolution
    Explicit methodological shift from SPEC root-finding (§2.1–2.2); central to claiming robustness.
  • domain assumption Beltrami solver inherited from SPEC is sufficiently accurate that force errors are dominated by geometry/resolution not by the linear field solve
    Stated in §2.3; Beltrami residual is monitored (Fig. 7a) as the limiter on ⟨F⟩.
  • ad hoc to paper Fixed plasma boundary extracted from an approximate HINT core surface is fair enough for qualitative island comparison
    §3.5 setup; authors themselves call for future free-boundary SPECTRE vs HINT.
invented entities (1)
  • SPECTRE solver (real-space force + trust-region LS MRxMHD code) independent evidence
    purpose: Compute stepped-pressure 3D equilibria more robustly than SPEC on strongly shaped geometries
    Software/method entity, not a new physical object; independent handle is public code and cross-code benchmarks.

pith-pipeline@v1.2.0-grok45-kimik3 · 18154 in / 3615 out tokens · 62904 ms · 2026-07-30T11:38:46.556663+00:00 · methodology

0 comments
read the original abstract

We present SPECTRE, a new equilibrium code based on the Multi-Region relaxed MHD model for robustly calculating 3D equilibria with general magnetic topology, allowing for flux surfaces, magnetic islands, and chaos. The code builds on a previous MRxMHD solver, SPEC, but performs significantly better thanks to a new formulation of force, the use of a stable trust-region-based least squares minimization scheme, as well as several additional features. SPECTRE is verified through application to configurations with known equilibrium solutions, in vacuum and with finite beta, both in the fixed boundary and the free boundary mode. Notably, these include vacuum equilibria of a quasi-axisymmetric (QA) device in fixed-boundary mode, and of W7-X in the free-boundary mode, along with a classical stellarator finite-beta free-boundary case. Finally, the solver is applied to a modern optimized finite-beta quasi-isodynamic (QI) configuration, where we demonstrate calculation of a strongly-shaped equilibrium with a core island which is in agreement with a HINT calculation.

Figures

Figures reproduced from arXiv: 2607.27135 by A. Goodman, C. Lazzati, C. Smiet, E. Balkovic, E. Lanti, J. Geiger, J. Loizu, J. P. Graves, R. Ramasamy.

Figure 1
Figure 1. Figure 1: Spectrum of the force for an optimized equilibrium in SPEC where only components inside the red box are targeted using root finding (data from Baillod (2023)) 2. New framework for minimization Inspired by the strong theoretical basis and the wide applicability of the MRxMHD model, we have developed a more robust framework for computing MRxMHD equilibria. This has resulted in the development of a new code, … view at source ↗
Figure 2
Figure 2. Figure 2: Evolution of ⟨F⟩ during an optimization, performed in either cylindrical Fourier R, Z (orange) or Henneberg space (red). The optimization corresponds to the calculation of the QA vacuum equilibrium from section 3.2. change the approach and use trust-region-based methods (Virtanen et al. 2020; Vogklis & Lagaris 2019) to minimize the real-space force in order to find equilibrium states in SPECTRE. The trust-… view at source ↗
Figure 3
Figure 3. Figure 3: The cross section of the D-shaped tokamak test case, including the VMEC flux surfaces, SPECTRE interfaces (Nvol = 8), and matching magnetic axes from the two codes. SPEC Beltrami fields (Hudson et al. 2012; Qu et al. 2020) carry over to SPECTRE. The code is written in modern Python and links to a Fortran-based extension for calculating the Beltrami fields. 3. Verification calculations Next, we present a se… view at source ↗
Figure 4
Figure 4. Figure 4: Convergence of the magnetic axis position in axisymmetry as a function of Nvol for a) β = 0 and b) β = 6%. A similar scan in c), now as a function of the Fourier resolution, with the y-axis showing the difference in axis position relative to the largest resolution. 3.2. QA vacuum The following verification case is of a QA reactor-scale configuration optimized for neoclassical confinement and stability (Hen… view at source ↗
Figure 5
Figure 5. Figure 5: a) Initial and optimized SPECTRE interfaces demonstrating robustness for a poor initial geometry guess, b) Comparison of SPECTRE interfaces and Poincaré trace with VMEC flux surfaces for the vacuum QA configuration [PITH_FULL_IMAGE:figures/full_fig_p008_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: A 3D depiction of the interfaces and the plasma boundary (outermost black surface) for the optimized QA vacuum equilibrium, with the magnetic axis shown in green and the color representing magnetic field strength. which is very close to the one predicted by VMEC (≈1mm). The VMEC calculation was performed with mpol = ntor = 22, NS = 220, Ftol = 10−16. The SPECTRE force minimization required 9 cpu-hours (1.3… view at source ↗
Figure 7
Figure 7. Figure 7: Convergence trends for the QA vacuum equilibrium as a function of the Fourier resolution: a) Error in the Beltrami field (maximum of per-volume error), b) Residual ⟨F⟩, c) Poincaré error between 1-volume and 7-volume cases, d) Magnetic axis position (VMEC axis shown in dashed black). 5 6 7 8 9 10 11 12 13 R [m] −2 −1 0 1 2 Z [m] Poincare (Nvol = 1) Interfaces (Nvol = 7) Plasma boundary [PITH_FULL_IMAGE:fi… view at source ↗
Figure 8
Figure 8. Figure 8: Comparison of Poincaré trace for a single-volume Beltrami field and a 7-volume SPECTRE solution with interfaces in force balance. configuration OP 1.1 (Wolf et al. 2017). The domain consists of a single toroidal plasma volume and an annular vacuum volume, solved with a different boundary condition on the outer surface B · nˆ = Bcoils · nˆ + Bplasma · nˆ ̸= 0 that is referred to as the computational boundar… view at source ↗
Figure 9
Figure 9. Figure 9: Comparison of Poincaré trace from the Biot-Savart evaluation of the coils, and from the SPECTRE free-boundary solution at a) ϕ = 0 and b) ϕ = π/10. Iota profile at θ = ϕ = 0 of the two fields shown in c). volume requires Picard iteration (Hudson et al. 2025) since the Beltrami field is solved for a given normal field, which itself depends on the Beltrami field. Figures 9a and 9b show the Poincaré plots of … view at source ↗
Figure 10
Figure 10. Figure 10: a) The input pressure profiles normalized to beta on axis, b) the calculated iota profiles of the rotating ellipse equilibria in VMEC and SPECTRE. 9.0 9.5 10.0 10.5 11.0 R [m] −0.8 −0.6 −0.4 −0.2 0.0 0.2 0.4 0.6 0.8 Z [m] VMEC VMEC axis SPEC SPECTRE Computational boundary B~ from coils [PITH_FULL_IMAGE:figures/full_fig_p011_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: The setup of the free-boundary rotating-ellipse example. Includes VMEC flux surfaces and the interfaces from SPEC and SPECTRE (exactly coincide). We study this configuration with a given pressure profile (⟨β⟩volume = 1.5%, βaxis = 3%) and zero net toroidal current (which is characteristic of QI configurations). Due to the lack of a good reference solution such as in vacuum equilibria, the goal of this stu… view at source ↗
Figure 12
Figure 12. Figure 12: Calculated QI finite-β equilibria. a) Comparison of VMEC flux surfaces and SPECTRE interfaces b) Poincaré traces of the equilibrium field from SPECTRE (red) and HINT (blue), both depicting the 8/9 core island chain. The black line is the fixed plasma boundary used in VMEC and SPECTRE [PITH_FULL_IMAGE:figures/full_fig_p013_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Rotational transform profiles of the QI configuration found using VMEC (blue line) and SPECTRE (red points). SPECTRE iota shows flattening around the ι = 8/9 resonance where the island exists. other relevant stellarator concepts, as well as comparisons to initial-value resistive MHD codes such as JOREK (Hoelzl et al. 2021) and M3D-C1 (Ferraro et al. 2018). We hope that SPECTRE can help bring about the tru… view at source ↗

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

57 extracted references · 25 canonical work pages

  1. [1]

    Hirshman, S. P. and Whitson, J. C. , title = ". The Physics of Fluids , volume =. 1983 , month =

    1983

  2. [2]

    Muraglia and A

    M. Muraglia and A. Poyé and O. Agullo and N. Dubuit and X. Garbet , doi =. Nonlinear dynamics of NTM seeding by turbulence , volume =. Plasma Physics and Controlled Fusion , keywords =

  3. [3]

    Loizu and Y

    J. Loizu and Y. M. Huang and S. R. Hudson and A. Baillod and A. Kumar and Z. S. Qu , doi =. Direct prediction of nonlinear tearing mode saturation using a variational principle , volume =. Physics of Plasmas , month =

  4. [4]

    Journal of Plasma Physics , author=

    Nonlinear saturation of resistive tearing modes in a cylindrical tokamak with and without solving the dynamics , volume=. Journal of Plasma Physics , author=. 2023 , pages=. doi:10.1017/S0022377823000934 , number=

    work page doi:10.1017/s0022377823000934 2023

  5. [5]

    S. R. Hudson and R. L. Dewar and G. Dennis and M. J. Hole and M. McGann and G. von Nessi and S. Lazerson , doi =. Computation of multi-region relaxed magnetohydrodynamic equilibria , volume =. Physics of Plasmas , month =

  6. [6]

    Tokamak Magnetohydrodynamic Equilibrium States with Axisymmetric Boundary and a 3D Helical Core , author =. Phys. Rev. Lett. , volume =. 2010 , month =

    2010

  7. [7]

    X -Point Collapse and Saturation in the Nonlinear Tearing Mode Reconnection , author =. Phys. Rev. Lett. , volume =. 2005 , month =

    2005

  8. [8]

    2004 , issn =

    Simple and rigorous solution for the nonlinear tearing mode , journal =. 2004 , issn =. doi:https://doi.org/10.1016/j.physleta.2004.02.010 , author =

    work page doi:10.1016/j.physleta.2004.02.010 2004

  9. [9]

    2023 , series =

    Fitzpatrick, Richard , title =. 2023 , series =

    2023

  10. [10]

    2020 , url =

    Towards integrated control of tokamak plasmas: physics-based control of neoclassical tearing modes in the TCV tokamak , author =. 2020 , url =

    2020

  11. [11]

    Introduction to Plasma Physics , year =

    R.J Goldston , doi =. Introduction to Plasma Physics , year =

  12. [12]

    Rutherford, P. H. , title = ". The Physics of Fluids , volume =. 1973 , month =

    1973

  13. [13]

    Helical temperature perturbations associated with tearing modes in tokamak plasmas , volume =

    Richard Fitzpatrick , doi =. Helical temperature perturbations associated with tearing modes in tokamak plasmas , volume =. Physics of Plasmas , month =

  14. [14]

    1999 , month =

    S Günter and A Gude and M Maraschek and Q Yu and the ASDEX Upgrade Team , title =. 1999 , month =. doi:10.1088/0741-3335/41/6/306 , url =

    work page doi:10.1088/0741-3335/41/6/306 1999

  15. [15]

    Bhattacharjee , title =

    Yi-Min Huang and A. Bhattacharjee , title =. 2016 , month =. doi:10.3847/0004-637X/818/1/20 , url =

    doi:10.3847/0004-637x/818/1/20 2016

  16. [16]

    Relaxation of Toroidal Plasma and Generation of Reverse Magnetic Fields , author =. Phys. Rev. Lett. , volume =. 1974 , month =

    1974

  17. [17]

    Journal of Plasma Physics , author=

    Stepped pressure profile equilibria in cylindrical plasmas via partial Taylor relaxation , volume=. Journal of Plasma Physics , author=. 2006 , pages=. doi:10.1017/S0022377806005861 , number=

    work page doi:10.1017/s0022377806005861 2006

  18. [18]

    , year =

    Militello, Fulvio and Porcelli, F. , year =. Simple analysis of the nonlinear saturation of the tearing mode , volume =. Physics of Plasmas , doi =

  19. [19]

    , journal =

    Kumar, A. , journal =. doi:10.1088/1361-6587/ac53ee , year =

    work page doi:10.1088/1361-6587/ac53ee

  20. [20]

    Journal of Plasma Physics , author=

    Computation of multi-region, relaxed magnetohydrodynamic equilibria with prescribed toroidal current profile , volume=. Journal of Plasma Physics , author=. 2021 , pages=. doi:10.1017/S0022377821000520 , number=

    work page doi:10.1017/s0022377821000520 2021

  21. [21]

    doi:10.1088/0029-5515/55/11/113014 , year =

    work page doi:10.1088/0029-5515/55/11/113014

  22. [22]

    Watanabe and S

    K.Y. Watanabe and S. Sakakibara and Y. Narushima and H. Funaba and K. Narihara and K. Tanaka and T. Yamaguchi and K. Toi and S. Ohdachi and O. Kaneko and H. Yamada and Y. Suzuki and W.A. Cooper and S. Murakami and N. Nakajima and I. Yamada and K. Kawahata and T. Tokuzawa and A. Komori and LHD experimental group , title =. doi:10.1088/0029-5515/45/11/004 , year =

    work page doi:10.1088/0029-5515/45/11/004

  23. [23]

    Helander, Per and Sigmar, D. J. , year=. Collisional transport in magnetized plasmas , publisher=

  24. [24]

    doi:10.1088/0741-3335/41/12B/312 , year =

    Mitchell A Berger , title =. doi:10.1088/0741-3335/41/12B/312 , year =

    work page doi:10.1088/0741-3335/41/12b/312

  25. [25]

    Hoelzl and G.T.A

    M. Hoelzl and G.T.A. Huijsmans and S.J.P. Pamela and M. Bécoulet and E. Nardon and F.J. Artola and B. Nkonga and C.V. Atanasiu and V. Bandaru and A. Bhole and D. Bonfiglio and A. Cathey and O. Czarny and A. Dvornova and T. Fehér and A. Fil and E. Franck and S. Futatani and M. Gruca and H. Guillard and J.W. Haverkort and I. Holod and D. Hu and S.K. Kim and...

  26. [26]

    Teng and N

    Q. Teng and N. Ferraro and D.A. Gates and R.B. White , title =. doi:10.1088/1741-4326/aad7c9 , year =

    work page doi:10.1088/1741-4326/aad7c9

  27. [27]

    Glasser, A. H. and Greene, J. M. and Johnson, J. L. , title = ". The Physics of Fluids , volume =. 1976 , month =. doi:10.1063/1.861490 , url =

    doi:10.1063/1.861490 1976

  28. [28]

    2022 , month =

    M Kong and F Felici and O Sauter and C Galperti and T Vu and C J Ham and T C Hender and M Maraschek and M Reich and the TCV team and the ASDEX Upgrade team and the MAST team and the EUROfusion MST1 Team , title =. 2022 , month =. doi:10.1088/1361-6587/ac48be , url =

    work page doi:10.1088/1361-6587/ac48be 2022

  29. [29]

    E., Weller, A., Zarnstorff, M

    Geiger J. E., Weller, A., Zarnstorff, M. C., Nührenberg, C., Werner, A. H. F., Kolesnichenko, Y. I.,. Equilibrium and Stability of High-B Plasmas in Wendelstein 7-AS , journal =. 2004 , publisher =

    2004

  30. [30]

    and Friedrichs, K

    Courant, R. and Friedrichs, K. and Lewy, H. , title=. 1967 , volume=

    1967

  31. [31]

    Brunetti and J.P

    D. Brunetti and J.P. Graves and W.A. Cooper and D. Terranova , title =. doi:10.1088/0029-5515/54/6/064017 , year =

    work page doi:10.1088/0029-5515/54/6/064017

  32. [32]

    doi:10.1088/0741-3335/44/9/315 , year =

    O Sauter and R J Buttery and R Felton and T C Hender and D F Howell and contributors to the EFDA-JET Workprogramme , title =. doi:10.1088/0741-3335/44/9/315 , year =

    work page doi:10.1088/0741-3335/44/9/315

  33. [33]

    Relaxation and magnetic reconnection in plasmas , author =. Rev. Mod. Phys. , volume =. 1986 , month =. doi:10.1103/RevModPhys.58.741 , url =

    doi:10.1103/revmodphys.58.741 1986

  34. [34]

    The Physics of Fluids , volume =

    Grad, Harold , title =. The Physics of Fluids , volume =. 1967 , month =

    1967

  35. [35]

    doi:10.1088/0029-5515/46/11/L01 , year =

    Suzuki, Yasuhiro and Nakajima, Noriyoshi and Watanabe, Kiyomasa and Nakamura, Yuji and Hayashi, Takaya , title =. doi:10.1088/0029-5515/46/11/L01 , year =

    doi:10.1088/0029-5515/46/11/l01

  36. [36]

    and Teukolsky, Saul A

    Press, William H. and Teukolsky, Saul A. and Vetterling, William T. and Flannery, Brian P. , title =. 2007 , isbn =

    2007

  37. [37]

    2025 , month =

    Babin, Robert and Hindenlang, Florian and Maj, Omar and Köberl, Robert , title =. 2025 , month =. doi:10.1088/1361-6587/adaee6 , url =

    work page doi:10.1088/1361-6587/adaee6 2025

  38. [38]

    Coordinate parameterisation and spectral method optimisation for Beltrami field solver in stellarator geometry , volume =

    Qu, Zhisong and Pfefferlé, David and Hudson, S and Baillod, A and Kumar, A and Dewar, Robert and Hole, Matthew , year =. Coordinate parameterisation and spectral method optimisation for Beltrami field solver in stellarator geometry , volume =. Plasma Physics and Controlled Fusion , doi =

  39. [39]

    and Drevlak, M

    Henneberg, S.A. and Drevlak, M. and Nührenberg, C. and Beidler, C.D. and Turkin, Y. and Loizu, J. and Helander, P. , title =. 2019 , month =. doi:10.1088/1741-4326/aaf604 , url =

    work page doi:10.1088/1741-4326/aaf604 2019

  40. [40]

    Hudson, S. R. and Panici, D. and Zhu, C. and Woodbury Saudeau, G. S. and Baillod, A. and Cianciosa, M. and Ware, A. S. , title =. Physics of Plasmas , volume =. 2025 , month =

    2025

  41. [41]

    and Ali, A

    Wolf, R.C. and Ali, A. and Alonso, A. and Baldzuhn, J. and Beidler, C. and Beurskens, M. and Biedermann, C. and Bosch, H.-S. and Bozhenkov, S. and Brakel, R. and Dinklage, A. and Feng, Y. and Fuchert, G. and Geiger, J. and Grulke, O. and Helander, P. and Hirsch, M. and Höfel, U. and Jakubowski, M. and Knauer, J. and Kocsis, G. and König, R. and Kornejew, ...

    work page doi:10.1088/1741-4326/aa770d 2017

  42. [42]

    2024 , month =

    Balkovic, E and Loizu, J and Graves, J P and Huang, Y-M and Smiet, C , title =. 2024 , month =. doi:10.1088/1361-6587/ad97dd , url =

    work page doi:10.1088/1361-6587/ad97dd 2024

  43. [43]

    Journal of Plasma Physics , author=

    Representing the boundary of stellarator plasmas , volume=. Journal of Plasma Physics , author=. 2021 , pages=. doi:10.1017/S0022377821000891 , number=

    doi:10.1017/s0022377821000891 2021

  44. [44]

    Hirshman, S. P. and Breslau, J. , title =. Physics of Plasmas , volume =. 1998 , month =. doi:10.1063/1.872954 , url =

    work page doi:10.1063/1.872954 1998

  45. [45]

    and Haberland, Matt and Reddy, Tyler and Cournapeau, David and Burovski, Evgeni and Peterson, Pearu and Weckesser, Warren and Bright, Jonathan and

    Virtanen, Pauli and Gommers, Ralf and Oliphant, Travis E. and Haberland, Matt and Reddy, Tyler and Cournapeau, David and Burovski, Evgeni and Peterson, Pearu and Weckesser, Warren and Bright, Jonathan and. Nature Methods , year =

  46. [46]

    A Rectangular Trust-Region Approach for Unconstrained and Box-Constrained Optimization Problems , isbn =

    Vogklis, Konstantinos and Lagaris, Isaac , year =. A Rectangular Trust-Region Approach for Unconstrained and Box-Constrained Optimization Problems , isbn =

  47. [47]

    Journal of Plasma Physics , author=

    Equilibrium beta-limits dependence on bootstrap current in classical stellarators , volume=. Journal of Plasma Physics , author=. 2023 , pages=. doi:10.1017/S0022377823000910 , number=

    doi:10.1017/s0022377823000910 2023

  48. [48]

    Equilibrium pressure limits in stellarators , url=

    Baillod, Antoine , year=. Equilibrium pressure limits in stellarators , url=. doi:10.5075/epfl-thesis-10070 , school=

    work page doi:10.5075/epfl-thesis-10070

  49. [49]

    Soviet Journal of Experimental and Theoretical Physics , year = 1958, month = jan, volume =

    On Magnetohydrodynamical Equilibrium Configurations. Soviet Journal of Experimental and Theoretical Physics , year = 1958, month = jan, volume =

    1958

  50. [50]

    Dennis, G. R. and Hudson, S. R. and Dewar, R. L. and Hole, M. J. , title =. Physics of Plasmas , volume =. 2013 , month =. doi:10.1063/1.4795739 , url =

    work page doi:10.1063/1.4795739 2013

  51. [51]

    and Loizu, J

    Baillod, A. and Loizu, J. and Graves, J. P. and Landreman, M. , title =. Physics of Plasmas , volume =. 2022 , month =. doi:10.1063/5.0080809 , url =

    work page doi:10.1063/5.0080809 2022

  52. [52]

    Smiet, C. B. and Loizu, J. and Balkovic, E. and Baillod, A. , title =. Physics of Plasmas , volume =. 2025 , month =. doi:10.1063/5.0226402 , url =

    work page doi:10.1063/5.0226402 2025

  53. [53]

    Qu, Z. S. and Zhou, Y. and Kumar, A. and Doak, J. and Loizu, J. and Hole, M. J. , title =. Physics of Plasmas , volume =. 2025 , month =. doi:10.1063/5.0260347 , url =

    work page doi:10.1063/5.0260347 2025

  54. [54]

    2020 , month =

    Hudson, S R and Loizu, J and Zhu, C and Qu, Z S and Nührenberg, C and Lazerson, S and Smiet, C B and Hole, M J , title =. 2020 , month =. doi:10.1088/1361-6587/ab9a61 , url =

    work page doi:10.1088/1361-6587/ab9a61 2020

  55. [55]

    doi:10.1088/0034-4885/77/8/087001 , year =

    Helander, Per , title =. doi:10.1088/0034-4885/77/8/087001 , year =

    doi:10.1088/0034-4885/77/8/087001

  56. [56]

    Journal of Plasma Physics , author=

    Equilibrium beta-limits in classical stellarators , volume=. Journal of Plasma Physics , author=. doi:10.1017/S0022377817000861 , number=

    doi:10.1017/s0022377817000861

  57. [57]

    2018 , month =

    M3D-C1 , author =. 2018 , month =. doi:10.11578/dc.20180627.9 , url =

    work page doi:10.11578/dc.20180627.9 2018