pith. sign in

arxiv: 2605.16734 · v1 · pith:YYHFQWNBnew · submitted 2026-05-16 · ⚛️ physics.plasm-ph

FIRM3D: Fast ion reduced models in 3D

Elizabeth Paul , Alexey Knyazev , Michael Czekanski , Alexandra Lachmann , Abdullah Hyder , Christopher Albert , Matt Landreman This is my paper

Pith reviewed 2026-05-19 20:08 UTC · model grok-4.3

classification ⚛️ physics.plasm-ph
keywords energetic particlesguiding-center orbits3D magnetic fieldsfusion plasmasstellaratorsorbit integrationMHD wavesopen-source software
0
0 comments X

The pith

FIRM3D supplies an open-source suite that integrates guiding-center orbits of energetic particles through three-dimensional magnetic fields while linking to MHD equilibrium and wave codes.

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

The paper introduces FIRM3D as a standalone Python/C++/CUDA framework for following the collisionless motion of energetic particles in fusion devices. Particles born from fusion or heating follow orbits that the code computes efficiently by solving the guiding-center equations in prescribed electromagnetic fields. The software adds parallel CPU and GPU integrators, symplectic and Runge-Kutta options, and direct interfaces to external tools that supply equilibrium fields and wave perturbations. These features let researchers run focused studies of particle transport and losses without the overhead of full-device simulation packages.

Core claim

FIRM3D performs parallelized integration of the guiding-center orbit equation for Monte Carlo samples drawn from energetic-particle distributions, supplies orbit visualization and transport diagnostics such as Poincaré maps and weighted Birkhoff averages, and includes ready interfaces to BOOZ_XFORM for MHD equilibria and to AE3D and FAR3D for wave stability data.

What carries the argument

The extended guiding-center orbit integrator, grown from SIMSOPT routines, that advances particle trajectories in 3D fields with optional symplectic or Runge-Kutta steps and feeds results into classification and averaging diagnostics.

If this is right

  • Monte Carlo ensembles of energetic-particle orbits can be integrated at scale on both CPUs and GPUs.
  • Resonant transport driven by MHD waves can be quantified by coupling the orbit solver directly to wave-stability outputs.
  • Poincaré maps and orbit-classification routines become available for rapid identification of confined versus lost particles.
  • New physics modules for particle-wave interactions or additional diagnostics can be added with few external dependencies.

Where Pith is reading between the lines

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

  • Stellarator design teams could insert FIRM3D into optimization loops to penalize configurations that produce large fast-ion losses.
  • GPU scaling would permit statistical sampling of rare loss events across thousands of particles in a single run.
  • The same orbit infrastructure might later incorporate finite-Larmor-radius corrections or self-consistent wave evolution without rewriting the core integrator.

Load-bearing premise

The guiding-center approximation remains valid for the energetic particles of interest and external codes supply accurate, compatible field and wave data without large numerical artifacts.

What would settle it

A side-by-side comparison of FIRM3D guiding-center trajectories against full-orbit Lorentz integrations for a test particle whose gyroradius becomes comparable to the magnetic-field scale length.

Figures

Figures reproduced from arXiv: 2605.16734 by Abdullah Hyder, Alexandra Lachmann, Alexey Knyazev, Christopher Albert, Elizabeth Paul, Matt Landreman, Michael Czekanski.

Figure 1
Figure 1. Figure 1: Left: Energy as a function of time for an alpha particle in the 𝛽 = 2.5% Landreman QH configuration. The Dormand-Prince algorithm exhibits a net energy drift over time, while the symplectic algorithm exhibits a stable moving time average of the energy (over 10−4 seconds). Right: Relative error in canonical momentum 𝑃𝜂 conservation for a perfectly quasisymmetric field. 10 particles are traced in the same co… view at source ↗
Figure 2
Figure 2. Figure 2: Left: Comparison of a trapped 3.5 MeV alpha particle orbit in the precise QH equilibrium. The difference in the 𝑠 coordinate between the trajectories at 10−3 seconds is 7.8× 10−3. Right: Comparison of loss fraction for the precise QA equilibrium. 5000 3.5 MeV alpha particles are sampled from a fusion birth distribution function and traced for 10−2 seconds. The two codes report identical loss fractions at t… view at source ↗
Figure 3
Figure 3. Figure 3: Scaling of tracing a fusion birth distribution in the Wistell-A equilibrium on 1 Perlmutter CPU node (128 CPU threads) and 1 NVIDIA A100, as a function of number of Monte Carlo samples [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Left: Trapped particle Poincaré map showing chaotic layers responsible for banana-drift diffusion. Right: Measures of convective and diffusive transport indicate banana-trapped orbits undergo banana diffusion [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: A kinetic Poincaré plot colored by the digit accuracy of the weighted Birkhoff average, an indicator of chaos. FIRM3D’s orbit classification, transport diagnostics, and AE-induced transport capabilities are illustrated in [PITH_FULL_IMAGE:figures/full_fig_p005_5.png] view at source ↗
read the original abstract

The dynamics of energetic particle (EP) species, born from fusion reactions or plasma heating schemes, are critical for predicting the behavior of magnetic confinement fusion experiments and future fusion reactors. Because energetic particles are largely collisionless, the orbits of Monte Carlo samples drawn from a given distribution function can be efficiently integrated in prescribed electromagnetic fields. In addition to the static magneto-hydrodynamic (MHD) equilibrium fields produced by the electromagnetic coils of a fusion device, MHD waves can be excited by -- and resonantly transport -- energetic particle populations. FIRM3D is an open-source Python/C++/CUDA software suite for modeling energetic particle dynamics in 3D magnetic fields, available at https://github.com/ColumbiaStellaratorTheory/firm3d. The core guiding-center integration routines grew out of SIMSOPT (Landreman et al., 2021), but have been extended to include additional physics and diagnostics not typically required in the stellarator optimization context. This standalone framework enables focused development of energetic particle physics capabilities with minimal dependencies, making it accessible to the broader stellarator and plasma physics community. Components of FIRM3D include interfaces with MHD equilibrium and wave stability software (BOOZ_XFORM, AE3D, FAR3D); CPU and GPU parallelized integration of the guiding center orbit equation, with symplectic and Runge-Kutta integrator options; and orbit visualization and transport diagnostics, including Poincare maps, orbit classification, and weighted Birkhoff averaging.

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

0 major / 2 minor

Summary. The manuscript describes FIRM3D, an open-source Python/C++/CUDA software suite for modeling energetic particle dynamics in 3D magnetic fields. It extends guiding-center integration routines originally from SIMSOPT with added EP physics, diagnostics, and interfaces to MHD equilibrium and wave stability codes (BOOZ_XFORM, AE3D, FAR3D). Core capabilities include CPU/GPU-parallelized integration of the guiding-center orbit equation using symplectic and Runge-Kutta options, plus orbit visualization and transport diagnostics such as Poincaré maps, orbit classification, and weighted Birkhoff averaging.

Significance. If the described implementation and interfaces function as stated, FIRM3D would provide a valuable, focused, standalone framework for energetic particle studies in stellarators and other 3D configurations. Strengths include its open-source distribution, minimal dependencies, parallelization support, and explicit linkages to established tools, which promote accessibility and reproducibility for the broader plasma physics community.

minor comments (2)
  1. Abstract: The diagnostic 'weighted Birkhoff averaging' is referenced without a short definition or citation; adding one sentence of explanation would improve accessibility for readers outside the immediate subfield.
  2. The manuscript should include a brief section or table summarizing verification tests (e.g., comparison of integrators against analytic orbits or conservation properties) to substantiate the numerical accuracy claims.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for their positive summary of the FIRM3D manuscript and for recognizing its potential value as a focused, open-source framework for energetic particle studies in 3D configurations. We appreciate the recommendation for minor revision and have prepared corresponding updates to improve clarity and accessibility.

Circularity Check

0 steps flagged

No significant circularity; software description with no derivation chain

full rationale

The paper presents FIRM3D as an open-source software framework for guiding-center orbit integration in 3D fields, extending SIMSOPT routines with added diagnostics and interfaces to BOOZ_XFORM, AE3D, and FAR3D. No equations, predictions, or first-principles results are claimed that could reduce by construction to fitted parameters, self-citations, or ansatzes. The central content is a tool description and capability summary rather than a mathematical derivation; standard modeling prerequisites such as guiding-center validity are stated without internal reduction to the paper's own outputs. This is a self-contained software distribution with no load-bearing steps that exhibit circularity.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

This is a software description paper; the central contribution is the packaged tool rather than new physical axioms or parameters. No free parameters or invented entities are introduced in the abstract.

axioms (1)
  • domain assumption Guiding-center approximation is appropriate for the energetic particle dynamics under consideration.
    Standard assumption invoked when stating that orbits can be efficiently integrated using guiding-center equations.

pith-pipeline@v0.9.0 · 5815 in / 1221 out tokens · 40541 ms · 2026-05-19T20:08:36.441431+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

24 extracted references · 24 canonical work pages · 1 internal anchor

  1. [1]

    G., Buchholz, R., Kasilov, S

    Albert, C. G., Buchholz, R., Kasilov, S. V., Kernbichler, W., & Rath, K. (2023). Alpha particle confinement metrics based on orbit classification in stellarators. Journal of Plasma Physics, 89(3), 955890301

    work page 2023

  2. [2]

    G., Kasilov, S

    Albert, C. G., Kasilov, S. V., & Kernbichler, W. (2020). Symplectic integration with non-canonical quadrature for guiding-center orbits in magnetic confinement devices. Journal of Computational Physics, 403, 109065

    work page 2020

  3. [3]

    Boost.Numeric.Odeint. (2025). Boost.odeint: Solving ordinary differential equations in C++ . https://www.boost.org/libs/numeric/odeint

    work page 2025

  4. [4]

    Boozer, A. H. (1995). Quasi-helical symmetry in stellarators. Plasma Physics and Controlled Fusion, 37(11A), A103

    work page 1995

  5. [5]

    J., & Hudson, S

    Chambliss, A., Paul, E. J., & Hudson, S. R. (2025). Fast particle trajectories and integrability in quasiaxisymmetric and quasihelical stellarators. Journal of Plasma Physics, 91(3), E74

    work page 2025

  6. [6]

    CATAPULT: A CUDA-Accelerated Timestepper for Alpha Particles Using Local Tricubics

    Czekanski, M., Knyazev, A. R., Bindel, D., & Paul, E. J. (2026). CATAPULT: A CUDA-accelerated timestepper for alpha particles using local tricubics. arXiv Preprint arXiv:2604.07617

    work page internal anchor Pith review Pith/arXiv arXiv 2026

  7. [7]

    Duignan, N., & Meiss, J. D. (2023). Distinguishing between regular and chaotic orbits of flows by the weighted birkhoff average. Physica D: Nonlinear Phenomena, 449, 133749

    work page 2023

  8. [8]

    N., Pinches, S

    Gorelenkov, N. N., Pinches, S. D., & Toi, K. (2014). Energetic particle physics in fusion research in preparation for burning plasma experiments. Nuclear Fusion, 54(12), 125001

    work page 2014

  9. [9]

    Jakob, W., Rhinelander, J., & Moldovan, D. (2017). pybind11 -- seamless operability between C++11 and Python . Zenodo. https://doi.org/10.5281/zenodo.3748217

    work page doi:10.5281/zenodo.3748217 2017

  10. [10]

    Knyazev, A., Lachmann, A., Goodman, A., Hyder, A., Czekanski, M., Spong, D., & Paul, E. (2026). On shear alfv é n wave-induced energetic ion transport in optimized stellarators. arXiv Preprint arXiv:2603.03118

    work page arXiv 2026

  11. [11]

    Landreman, M. (2021). Booz\_xform. Zenodo. https://doi.org/10.5281/zenodo.4876646

    work page doi:10.5281/zenodo.4876646 2021

  12. [12]

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

    work page 2022

  13. [13]

    Landreman, M., Medasani, B., Wechsung, F., Giuliani, A., Jorge, R., & Zhu, C. (2021). SIMSOPT: A flexible framework for stellarator optimization. Journal of Open Source Software, 6(65), 3525

    work page 2021

  14. [14]

    Landreman, M., & Paul, E. J. (2022). Magnetic fields with precise quasisymmetry for plasma confinement. Physical Review Letters, 128(3), 035001

    work page 2022

  15. [15]

    von der, Labbate, J., Paul, E., Flom, E., Dudt, D., Wu, R., Kruger, T., Swanson, C., Kumar, S., Gates, D., & others

    Linden, J. von der, Labbate, J., Paul, E., Flom, E., Dudt, D., Wu, R., Kruger, T., Swanson, C., Kumar, S., Gates, D., & others. (2026). Alpha confinement in the quasi-axisymmetric helios fusion power plant. Fusion Engineering and Design, 229, 115801

    work page 2026

  16. [16]

    J., Bhattacharjee, A., Landreman, M., Alex, D., Velasco, J

    Paul, E. J., Bhattacharjee, A., Landreman, M., Alex, D., Velasco, J. L., & Nies, R. (2022). Energetic particle loss mechanisms in reactor-scale equilibria close to quasisymmetry. Nuclear Fusion, 62(12), 126054

    work page 2022

  17. [17]

    J., Mynick, H

    Paul, E. J., Mynick, H. E., & Bhattacharjee, A. (2023). Fast-ion transport in quasisymmetric equilibria in the presence of a resonant alfv é nic perturbation. Journal of Plasma Physics, 89(5), 905890515

    work page 2023

  18. [18]

    Press, W. H. (2007). Numerical recipes 3rd edition: The art of scientific computing. Cambridge university press

    work page 2007

  19. [19]

    Rahbarnia, K., Thomsen, H., Schilling, J., Vaz Mendes, S., Endler, M., Kleiber, R., Könies, A., Borchardt, M., Slaby, C., Bluhm, T., & others. (2020). Alfv é nic fluctuations measured by in-vessel Mirnov coils at the Wendelstein 7-X stellarator . Plasma Physics and Controlled Fusion, 63(1), 015005

    work page 2020

  20. [20]

    A., D'azevedo, E., & Todo, Y

    Spong, D. A., D'azevedo, E., & Todo, Y. (2010). Clustered frequency analysis of shear alfv é n modes in stellarators. Physics of Plasmas, 17(2)

    work page 2010

  21. [21]

    Swanson, C., Kumar, S., Dudt, D., Flom, E., Kalb, W., Kruger, T., Martin, M., Olatunji, J., Pasmann, S., Tang, L., & others. (2025). Overview of the helios design: A practical planar coil stellarator fusion power plant. arXiv Preprint arXiv:2512.08027

    work page arXiv 2025

  22. [22]

    A., & Todo, Y

    Toi, K., Ogawa, K., Isobe, M., Osakabe, M., Spong, D. A., & Todo, Y. (2011). Energetic-ion-driven global instabilities in stellarator/helical plasmas and comparison with tokamak plasmas. Plasma Physics and Controlled Fusion, 53(2), 024008

    work page 2011

  23. [23]

    Varela, J., Spong, D., Garcia, L., Ghai, Y., Ortiz, J., & Collaborators, F. P. (2024). Stability optimization of energetic particle driven modes in nuclear fusion devices: The FAR3d gyro-fluid code. Frontiers in Physics, 12, 1422411

    work page 2024

  24. [24]

    Varje, J., Särkimäki, K., Kontula, J., Ollus, P., Kurki-Suonio, T., Snicker, A., Hirvijoki, E., & Äkäslompolo, S. (2019). High-performance orbit-following code ASCOT5 for Monte Carlo simulations in fusion plasmas. arXiv Preprint arXiv:1908.02482. CSLReferences document

    work page arXiv 2019