  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.
