Winkel, MathiasMathiasWinkelSpeck, RobertRobertSpeckRuprecht, DanielDanielRuprecht2021-10-142021-10-142015-08-15Journal of Computational Physics 295: 456-474 (2015-08-15)http://hdl.handle.net/11420/10522This work introduces the high-order Boris-SDC method for integrating the equations of motion for electrically charged particles in an electric and magnetic field. Boris-SDC relies on a combination of the Boris-integrator with spectral deferred corrections (SDC). SDC can be considered as preconditioned Picard iteration to compute the stages of a collocation method. In this interpretation, inverting the preconditioner corresponds to a sweep with a low-order method. In Boris-SDC, the Boris method, a second-order Lorentz force integrator based on velocity-Verlet, is used as a sweeper/preconditioner. The presented method provides a generic way to extend the classical Boris integrator, which is widely used in essentially all particle-based plasma physics simulations involving magnetic fields, to a high-order method. Stability, convergence order and conservation properties of the method are demonstrated for different simulation setups. Boris-SDC reproduces the expected high order of convergence for a single particle and for the center-of-mass of a particle cloud in a Penning trap and shows good long-term energy stability.en0021-99912015456474Boris integratorCollocation methodHigh-orderMagnetic fieldSpectral deferred corrections (SDC)Time integrationMathematics - Numerical AnalysisMathematics - Numerical AnalysisComputer Science - Numerical AnalysisJournal Article10.1016/j.jcp.2015.04.0221409.5677v2Other