Phase Stability Analysis of Volume-preserving Algorithms for Accurate Single Particle Orbit Simulations in Tokamak Plasmas

Second-order Volume-preserving algorithms (VPAs) for simulating charged particle motion in electromagnetic fields have been generalized to a rotating angle formulation by using the matrix decomposition methods. Based on this method, the phase stability of this class of VPAs has been analyzed by using the Discrete Fourier Transformations (DFT) technique. It is found that two prominent VPAs, namely the $G_h^2$ and the Boris algorithm, exhibit optimal phase precision for high-frequency (gyro motion) and low-frequency dynamics (transit/bounce motion), respectively. These findings have been empirically verified through numerical experiments. The insights gained from this study enable the selection of an appropriate VPA for practical simulations based on the characteristic frequencies of specific physics problems, which can substantially enhance numerical accuracy and improve computational efficiency for long-term simulations.