Papers › Structure-preserving second-order integration of relativistic charged particle...
Structure-preserving second-order integration of relativistic charged particle trajectories in electromagnetic fields
Adam V. Higuera, John R. Cary
The archive published only this paper's code-link row. Authors, date and abstract are from arXiv's metadata (CC0), read from the Kaggle arXiv metadata snapshot of 2026-09-12 where its title matched the archive's; the title is the archive's.
Time-centered, hence second-order, methods for integrating the relativistic momentum of charged particles in an electromagnetic field are derived. A new method is found by averaging the momentum before use in the magnetic rotation term, and an implementation is presented that differs from the relativistic Boris Push [1] only in the method for calculating the Lorentz factor. This is shown to have the same second-order accuracy in time as that (Boris Push) [1] found by splitting the electric acceleration and magnetic rotation and that [2] found by averaging the velocity in the magnetic rotation term. All three methods are shown to conserve energy when there is no electric field. The Boris method and the current method are shown to be volume-preserving, while the method of [2] and the current method preserve the E⃗ ×B⃗ velocity. Thus, of these second-order relativistic momentum integrations, only the integrator introduced here both preserves volume and gives the correct E⃗ ×B⃗ velocity. While all methods have error that is second-order in time, they deviate from each other by terms that increase as the motion becomes relativistic. Numerical results show that [2] develops energy errors near resonant orbits of a test problem that neither volume-preserving integrator does. [1] J. Boris, Relativistic plasma simulation-optimization of a hybrid code, in: Proc. Fourth Conf. Num. Sim. Plasmas, Naval Res. Lab, Wash. DC, 1970, pp. 3-67. [2] J.-L. Vay, Simulation of beams or plasmas crossing at relativistic velocity, Physics of Plasmas (1994-present) 15 (5) (2008) 056701.
Code
Repository list and official/mentioned flags are the archive's, frozen 2025-07-28. Reachability, where shown, is from one Syntology probe window (2026-09-16 to 2026-09-18); repositories not probed show nothing. GitHub stars are not tracked.
Code Syntology ran Syntology
Not run by Syntology. Nothing on this page verifies that the listed code works.
Results from the paper archive 2025-07-28
No leaderboard rows for this paper in the archive.
Report a problem or propose a change · a person checks every report against the paper or source before anything changes; decisions are listed on /corrections