Symmetric multistep methods for charged particle dynamics
Ernst Hairer and Christian Lubich
Abstract. A class of explicit symmetric multistep methods is proposed for
integrating the equations of motion of charged particles in an electro-magnetic
field. The magnetic forces are built into these methods in a special way that
respects the Lagrangian structure of the problem. It is shown that such methods
approximately preserve energy and momentum over very long times, proportional
to a high power of the inverse stepsize. We explain this behaviour by studying
the modified equation of the methods and by analysing the remarkably stable
propagation of parasitic solution components.
Key Words. Linear multistep method, charged particle, magnetic field,
energy conservation, backward error analysis, modified differential equation,
modulated Fourier expansion.