arxiv: 2604.20713 · v2 · submitted 2026-04-22 · ⚛️ physics.plasm-ph · math-ph· math.AP· math.MP

Recognition: unknown

An analytic formula for surface currents generating prescribed plasma equilibrium fields

Wadim Gerner

Authors on Pith no claims yet

Pith reviewed 2026-05-09 22:56 UTC · model grok-4.3

classification ⚛️ physics.plasm-ph math-phmath.APmath.MP

keywords plasma equilibriumsurface currentsanalytic formulaBiot-Savart operatorcoil winding surfacenormal magnetic curvatureplasma domainmagnetic field generation

0 comments

The pith

An analytic formula gives the surface current on a coil surface that exactly reproduces any prescribed plasma equilibrium field inside the enclosed volume when added to the plasma's own contribution.

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

The paper supplies a closed-form expression that takes a plasma domain, an equilibrium magnetic field defined inside it, and an enclosing coil winding surface, then returns the surface current distribution needed on that surface. The resulting currents produce a magnetic field that, together with the field generated by the plasma current density, equals the target equilibrium field at every point inside the plasma. A reader would care because the formula replaces iterative numerical searches for coil shapes with a direct calculation and simultaneously accounts for the tendency of current density to concentrate where the plasma boundary has the greatest normal magnetic curvature.

Core claim

Given a plasma domain P, an equilibrium field B on P, and a surrounding coil surface Σ, an analytic formula exists whose output is a surface current j on Σ such that the Biot-Savart field of j plus the Biot-Savart field of curl(B) equals B everywhere inside P.

What carries the argument

The analytic formula for the surface current j on Σ, which is constructed so that its Biot-Savart integral exactly cancels the difference between B and the Biot-Savart integral of the plasma current density curl(B) throughout the interior of P.

Load-bearing premise

A suitable coil winding surface Σ surrounding the plasma domain P exists and the Biot-Savart operator applied to the surface current and to curl(B) can be combined analytically to recover B exactly inside P.

What would settle it

Evaluate the Biot-Savart integral of the formula-derived j at an interior test point and verify whether it equals B minus the Biot-Savart integral of curl(B) at that point; any measurable discrepancy falsifies the formula.

read the original abstract

Given a plasma domain $P\subset\mathbb{R}^3$, a plasma equilibrium field $B$ on $P$ and a coil winding surface $\Sigma$ surrounding $P$, we provide an analytic formula whose output is a surface current distribution $j$ on $\Sigma$ such that $\operatorname{BS}(j)+\operatorname{BS}(\operatorname{curl}(B))=B$ in $P$, i.e. the combination of the plasma current magnetic field and the surface current magnetic field exactly produce the full plasma equilibrium field. We use this new formula to provide a theoretical explanation of the empirical phenomenon that currents tend to attain their maximal strength at points which correspond to points of highest normal magnetic curvature on the plasma boundary. Some discussions regarding aspects of numerical approximations are also included.

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.

Desk Editor's Note private letter to a colleague

The paper gives an explicit analytic formula for the surface current j on a winding surface that, together with the plasma current, produces any prescribed equilibrium B inside the volume.

read the letter

The paper's main result is a formula that takes a plasma domain P, an equilibrium field B inside it, and a surrounding surface Σ, then outputs a surface current j on Σ such that the Biot-Savart field from j plus the field from curl(B) equals B throughout P. It also uses the formula to derive why currents reach maximum strength at locations of highest normal curvature on the plasma boundary, which lines up with what people see in actual coil designs. That explanation is the part that feels most immediately useful. The work does well by staying focused on the inverse problem and by including a section on numerical approximations, which shows they know the formula will still need discretization in practice. The derivation appears to rest on standard vector calculus identities applied to the given domains rather than on fitting or iteration. The citation pattern is light but appropriate for a short analytic note. The soft spot is whether the formula is a genuine closed expression or an operator inversion that only looks analytic on paper. The Biot-Savart map from surface current to interior field is a nonlocal integral, and solving for j to match a prescribed interior field is a first-kind Fredholm problem that usually requires eigenfunction expansions or regularization. If the paper's derivation relies on a specific exterior extension or gauge choice that does not hold for every smooth Σ, the result is less general than the abstract suggests. The numerical-approximations discussion hints at this issue but does not fully resolve it. This paper is for plasma physicists and fusion engineers who work on coil design for stellarators or similar devices. A reader who wants an analytic handle on current distributions rather than another optimization code will find the formula and the curvature argument worth reading. It shows clear engagement with the inverse magnetostatics problem and is coherent on its own terms, so it deserves a serious referee. I would send it for review with the request that referees check the explicit steps of the derivation against the integral-operator concern.

Referee Report

2 major / 2 minor

Summary. Given a plasma domain P, an equilibrium field B on P, and a surrounding coil winding surface Σ, the paper provides an analytic formula for a surface current density j on Σ such that BS(j) + BS(curl(B)) = B inside P. It applies the formula to explain why surface currents tend to maximize at locations of highest normal magnetic curvature on the plasma boundary and includes discussion of numerical approximation aspects.

Significance. If the central formula is a true closed-form analytic expression valid for general smooth Σ and B, the result would be significant for stellarator and tokamak coil design by enabling direct (non-optimized) computation of required surface currents and by supplying a theoretical account of observed current distributions linked to boundary curvature. The attempt to derive an explicit expression rather than relying solely on numerical optimization is a positive feature, but its validity hinges on whether the Biot-Savart inversion is handled analytically without hidden series expansions or regularizations.

major comments (2)

[Abstract / formula derivation] Abstract and derivation (wherever the formula is stated): the claim that an 'analytic formula' for j exists such that BS(j) exactly equals B - BS(curl(B)) inside P for arbitrary smooth Σ is load-bearing. The Biot-Savart operator from a surface current is a nonlocal integral operator; its inversion on a general surface is a Fredholm integral equation of the first kind whose solution is not available in elementary closed form. The manuscript must exhibit the explicit expression for j and the vector-calculus steps that purportedly yield a closed-form result without eigenfunction expansions, regularization, or numerical inversion.
[Curvature explanation section] Application to curvature (section discussing maximal current strength): the explanation that currents attain maximal strength at points of highest normal magnetic curvature must be shown to follow directly from the formula without additional assumptions on the extension of B or the gauge. If the derivation relies on a specific continuation of the exterior field, this must be stated explicitly, as it would restrict the generality of the result.

minor comments (2)

[Abstract] The abstract refers to 'discussions regarding aspects of numerical approximations' but does not indicate the section number; a dedicated subsection or paragraph with concrete error estimates or convergence tests would improve clarity.
[Notation] Notation for the Biot-Savart operator BS(·) and the domain/surface symbols should be defined at first use with a brief reminder of the normalization (e.g., units for μ₀).

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for raising these important points about the explicitness of the central formula and the assumptions underlying the curvature discussion. We address each comment below and will revise the manuscript to provide additional detail and clarity.

read point-by-point responses

Referee: [Abstract / formula derivation] Abstract and derivation (wherever the formula is stated): the claim that an 'analytic formula' for j exists such that BS(j) exactly equals B - BS(curl(B)) inside P for arbitrary smooth Σ is load-bearing. The Biot-Savart operator from a surface current is a nonlocal integral operator; its inversion on a general surface is a Fredholm integral equation of the first kind whose solution is not available in elementary closed form. The manuscript must exhibit the explicit expression for j and the vector-calculus steps that purportedly yield a closed-form result without eigenfunction expansions, regularization, or numerical inversion.

Authors: We agree that the derivation steps and the explicit form of the formula require clearer presentation. In the revised manuscript we will add a dedicated subsection that walks through the vector-calculus derivation in full. The starting point is the observation that the required coil field inside P is the difference B − BS(curl B), which is both divergence-free and curl-free inside P. This harmonic vector field is then represented by an explicit surface-current expression on Σ obtained from the Biot-Savart integral and the jump relations for the tangential field; the resulting formula for j is a direct (non-iterative) integral expression that does not invoke eigenfunction expansions, regularization, or numerical matrix inversion. We will display the final expression for j and the intermediate identities so that readers can verify the closed-form character of the result. revision: yes
Referee: [Curvature explanation section] Application to curvature (section discussing maximal current strength): the explanation that currents attain maximal strength at points of highest normal magnetic curvature must be shown to follow directly from the formula without additional assumptions on the extension of B or the gauge. If the derivation relies on a specific continuation of the exterior field, this must be stated explicitly, as it would restrict the generality of the result.

Authors: We will revise the curvature section to state explicitly the extension chosen for B outside P and the gauge adopted for the vector potential. With these choices fixed, we will derive the local expression for |j| directly from the analytic formula and show that the leading term is proportional to the normal curvature of the plasma boundary; no further assumptions are required beyond those already stated for the extension. This will make the link between maximal current strength and boundary curvature fully transparent and self-contained. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation is self-contained

full rationale

The paper states it derives an analytic formula for j on Σ from given P, B, and Σ such that BS(j) + BS(curl(B)) equals the target field inside P. No quoted step reduces the claimed output to a fitted parameter, self-definition, or load-bearing self-citation that is itself unverified. The central result is presented as a direct construction rather than a renaming or inversion that is only formally labeled analytic. This matches the default expectation that most papers are not circular; the derivation chain stands on its own inputs without the forbidden patterns.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Based solely on the abstract, no explicit free parameters, axioms, or invented entities are identifiable; the formula is presented as analytic without disclosed fitting constants or new postulated objects.

pith-pipeline@v0.9.0 · 5425 in / 1047 out tokens · 29973 ms · 2026-05-09T22:56:48.772950+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

29 extracted references · 3 canonical work pages · 1 internal anchor

[1]

Aronszajn, A

N. Aronszajn, A. Krzywicki, and J. Szarski. A unique continuation theorem for exterior differential forms on Riemannian manifolds. Ark. Mat., 4:417–453, 1962

1962
[2]

C. B¨ ar. Zero sets of solutions to semilinear elliptic systems of first order. Inventiones mathemat- icae, 138:183–202, 1999

1999
[3]

Cantarella, D

J. Cantarella, D. DeTurck, and H. Gluck. The Biot-Savart operator for application in knot theory, fluid dynamics, and plasma physics. Journal of Mathematical Physics , 42(2):876–905, 2001

2001
[4]

Cantarella, D

J. Cantarella, D. DeTurck, and H. Gluck. Vector calculus and the topology of domains in 3-space. American Mathematical Monthly, 109(5):409–442, 2002

2002
[5]

Dalla Riva, M

M. Dalla Riva, M. Lanza de Cristoforis, and P. Musolino. Singularly Perturbed Boundary Value Problems. Springer, 2021

2021
[6]

Dudt and E

D.W. Dudt and E. Kolemen. DESC: A stellarator equilibrium solver. Phys. Plasmas, 27:102513, 2020

2020
[7]

Enciso and D

A. Enciso and D. Peralta-Salas. Knots and links in steady solutions of the Euler equation. Ann. of Math., 175(1):345–367, 2012

2012
[8]

W. Gerner. Properties of the Biot-Savart operator acting on surface currents. SIAM J. Math. Anal., 56(5):6446–6482, 2024

2024
[9]

W. Gerner. Non-optimal domains for the helicity maximisation problem. arXiv e-prints , page arXiv:2510.20533, October 2025

work page arXiv 2025
[10]

W. Gerner. Influence of coil geometry and coil-plasma distance on the magnetic field approxima- tion error. Plasma Phys. Control. Fusion , 68:045011, 2026

2026
[11]

W. Gerner. Kernel and image of the Biot-Savart operator and their applications in stellarator designs. J. Math. Phys. , 67:023504, 2026

2026
[12]

W. Gerner. Properties of surface helicity and its optimisation. Discrete and Continuous Dynamical Systems, 49:106–138, 2026

2026
[13]

Gerner, A

W. Gerner, A. Nicolopoulos-Salle, and D. Pereira Botelho. Electric current dynamics in the stellarator coil winding surface model. arXiv e-prints, page arXiv:2603.04085, March 2026

work page arXiv 2026
[14]

Girault and P.-A

V. Girault and P.-A. Raviart. Finite Element Methods for Navier-Stokes Equations . Springer, 1986

1986
[15]

J. D. Hanson. The virtual-casing principle and Helmholtz’s theorem. Plasma Phys. Control. Fusion, 57:115006, 2015

2015
[16]

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

1983
[17]

Kappell, M

J. Kappell, M. Landreman, and D. Malhotra. The magnetic gradient scale length explains why certain plasmas require close external magnetic coils. Plasma Phys. Control. Fusion , 66:025018, 2024

2024
[18]

Landreman

M. Landreman. An improved current potential method for fast computation of stellarator coil shapes. Nucl. Fusion, 57:046003, 2017

2017
[19]

E.L. Lima. The Jordan-Brouwer separation theorem for smooth hypersurfaces. The American Mathematical Monthly, 95(1):39–42, 1988

1988
[20]

Littlejohn

Robert G. Littlejohn. Variational principles of guiding centre motion. Journal of Plasma Physics , 29(1):111–125, 1983. 19

1983
[21]

K. Liu, Z. Zheng, and C. Zhu. Near-surface expansion of vacuum stellarator configurations. Nucl. Fusion, 66:016050, 2026

2026
[22]

Malhotra, A

D. Malhotra, A. J. Cerfon, M. O’Neil, and E. Toler. Efficient high-order singular quadrature schemes in magnetic fusion. Plasma Phys. Control. Fusion , 62:024004, 2020

2020
[23]

P. Merkel. Solution of stellarator boundary value problems with external currents. Nucl. Fusion, 27:867–871, 1987

1987
[24]

J.W. Milnor. Topology from the differentiable view point . Princeton University Press, 1965

1965
[25]

Perrella, D

D. Perrella, D. Pfefferl´ e, and L. Stoyanov. Rectifiability of divergence-free fields along invariant 2-tori. Partial Differential Equations and Applications , 3(50), 2022

2022
[26]

Rauch and M

J. Rauch and M. Taylor. Potential and Scattering Theory on Wildly Perturbed Domains. Journal of functional analysis , 18:27–59, 1975

1975
[27]

Estimating coil features from an equilibrium

E. Rodriguez and W. Sengupta. Estimating coil features from an equilibrium. arXiv e-prints , page arXiv:2604.12339, April 2026

work page internal anchor Pith review Pith/arXiv arXiv 2026
[28]

G. Schwarz. Hodge Decomposition - A Method for Solving Boundary Value Problems. Springer Verlag, 1995

1995
[29]

L. Spitzer. The Stellarator Concept. The Physics of Fluids , 1:253–264, 1958. 20

1958