An fully implicit scheme for numerical integration of the relativistic particle equation of motion

Abstract

Relativistic strongly magnetized plasmas are produced in laboratories thanks to state-of-the-art laser technology but can naturally be found around compact objects such as neutron stars and black holes. Detailed studies of the behaviour of relativistic plasmas require accurate computations able to catch the full spatial and temporal dynamics of the system. Numerical simulations of ultra-relativistic plasmas face severe restrictions due to limitations in the maximum possible Lorentz factors that current algorithm can reproduce to good accuracy. In order to circumvent this flaw and repel the limit to $\gamma \approx 10^9$, we design a new fully implicit scheme to solve the relativistic particle equation of motion in an external electromagnetic field using a three dimensional Cartesian geometry. We show some examples of numerical integrations in constant electromagnetic fields to prove the efficiency of our algorithm. The code is also able to follow the electric drift motion for high Lorentz factors. In the most general case of spatially and temporally varying electromagnetic fields, the code performs extremely well as shown by comparison with exact analytical solutions for the relativistic electrostatic Kepler problem as well as for linearly and circularly polarized plane waves.

Submitted to Journal of Plasma Physics Special Issue "The Vlasov equation: from space to laboratory plasmas"

References in corpus (2)

• High order volume-preserving algorithms for relativistic charged particles in general electromagnetic fields
• High Order Numerical Integrators for Relativistic Charged Particle Tracking

Cited by in corpus (10)

• A comprehensive comparison of relativistic particle integrators
• On the Boris solver in particle-in-cell simulation
• A novel energy-conserving scheme for eight-dimensional systems
• A relativistic particle pusher for ultra-strong electromagnetic fields
• A guiding center implementation for relativistic particle dynamics in the PLUTO code
• Particle acceleration in neutron star ultra-strong electromagnetic fields
• An Evaluation of Different Numerical Methods to Calculate the Pitch-angle Diffusion Coefficient from Full-orbit Simulations: disentangling a rope of sand
• Explicit volume-preserving numerical schemes for relativistic trajectories and spin dynamics
• Towards Modelling AR Sco: Generalised Particle Dynamics and Strong Radiation-Reaction Regimes
• A Comparative Analysis of Relativistic Particle Pushers vis-à-vis Computation Time & Accuracy