Researchers have introduced a novel residual-minimization method for preparing guiding-center distribution functions on finite grids, a foundational task for kinetic plasma simulations. The technique, detailed in a preprint, uses Polak-Ribiere+ iterations with exact quadratic line searches to minimize a fixed conservative discrete residual. A key feature is its matrix-free implementation; the operator and its adjoint are evaluated using JAX/XLA, which avoids the computationally expensive assembly of a four-dimensional matrix. This is achieved by batching independent magnetic-moment slices and compiling the iterations into a single execution graph, streamlining the calculation of equilibrium states for energetic particles in fusion devices. Source: arXiv plasm-ph
The method's primary goal is to find a distribution function f such that its time derivative, ∂f/∂t, is minimized, effectively reaching a steady state. The algorithm minimizes the norm of this residual in the physical Jacobian metric. The paper validates the workflow through an axisymmetric test case, using a NUBEAM fast-ion distribution and a corresponding constants-of-motion representation from EPCoM. This verification step confirms the method's ability to handle known physics before application to more complex, three-dimensional magnetic fields characteristic of modern stellarators and advanced tokamaks. Source: arXiv plasm-ph
The method's primary goal is to find a distribution function f such that its time derivative, ∂f/∂t, is minimized, effectively reaching a steady state.
Following the axisymmetric verification, the method was applied to several non-axisymmetric configurations, demonstrating its utility for contemporary fusion research. Calculations were performed for a precise quasi-axisymmetric (QA) geometry, the Helically Symmetric eXperiment (HSX), and the Wendelstein 7-X (W7-X) stellarator. For these three-dimensional cases, the method successfully reduced both the discrete residual and the finite-time variation of the distribution function when sampled along independently integrated particle trajectories. This demonstrates the method's capability to produce stable, collisionless steady-state distributions in intricate magnetic confinement systems. Source: arXiv plasm-ph
An important extension of the work addresses the physical constraint of non-negativity for the distribution function. This is treated as an optional constrained optimization problem. The authors note that enforcing non-negative projection does not alter the underlying residual operator but does change the set of feasible zero-residual solutions and the final endpoint selection. The results presented are described as finite-grid, finite-time evidence of a valid steady-state solution. This new tool offers a more robust and potentially efficient pathway for initializing high-fidelity kinetic simulations, which are critical for predicting and understanding energetic particle behavior and overall plasma performance in future fusion reactors. Source: arXiv plasm-ph