High Order Explicit Lorentz Invariant Volume-preserving Algorithms for Relativistic Dynamics of Charged Particles

Abstract: Lorentz invariant structure-preserving algorithms possess reference-independent secular stability, which is vital for simulating relativistic multi-scale dynamical processes. The splitting method has been widely used to construct structure-preserving algorithms, but without exquisite considerations, it can easily break the Lorentz invariance of continuous systems. In this paper, we introduce a Lorentz invariant splitting technique to construct high order explicit Lorentz invariant volume-preserving algorithms (LIVPAs) for simulating charged particle dynamics. Using this method, long-term stable explicit LIVPAs with different orders are developed and their performances of Lorentz invariance and long-term stability are analyzed in detail.