Comment on "Symplectic integration of magnetic systems": a proof that the Boris algorithm is not variational

Abstract

The Boris algorithm for integrating charged particle trajectories in electric and magnetic fields is popular due to its simple implementation, rapid iteration, and observed long-term numerical fidelity. The underlying cause of this long-term fidelity has become a matter of controversy, with one article claiming the method to be symplectic [S. D. Webb, J. Comput. Phys. 270 (2014) 570], and others claiming the method to be volume preserving but not symplectic [e.g. H. Qin et al., Phys. Plasmas 20 (2013) 084503]. To resolve the discrepancy, this letter leverages a discrete Helmholtz condition to demonstrate that no variational formulation of the Boris algorithm exists, indicating that the long-term fidelity should be attributed to the volume-preserving properties of the algorithm.

References in corpus (1)

• Symplectic Integration of Magnetic Systems

Cited by in corpus (14)

• A comprehensive comparison of relativistic particle integrators
• Explicit high-order non-canonical symplectic particle-in-cell algorithms for Vlasov-Maxwell systems
• Structure-Preserving Geometric Particle-in-Cell Methods for Vlasov-Maxwell Systems
• Explicit symplectic algorithms based on generating functions for charged particle dynamics
• High-order energy-conserving Line Integral Methods for charged particle dynamics
• Explicit high-order noncanonical symplectic algorithms for ideal two-fluid systems
• Explicit Structure-Preserving Geometric Particle-in-Cell Algorithm in Curvilinear Orthogonal Coordinate Systems and Its Applications to Whole-Device 6D Kinetic Simulations of Tokamak Physics
• Field theory and structure-preserving geometric particle-in-cell algorithm for drift wave instability and turbulence
• Explicit High-Order Gauge-Independent Symplectic Algorithms for Relativistic Charged Particle Dynamics
• Geometric Electrostatic Particle-In-Cell Algorithm on Unstructured Meshes
• Advanced fuel fusion, phase space engineering, and structure-preserving geometric algorithms
• Large-stepsize integrators for charged-particle dynamics over multiple time scales
• Efficient energy-preserving methods for charged-particle dynamics
• Exponential energy-preserving methods for charged-particle dynamics in a strong and constant magnetic field