A comprehensive comparison of relativistic particle integrators

Abstract

We compare relativistic particle integrators commonly used in plasma physics showing several test cases relevant for astrophysics. Three explicit particle pushers are considered, namely the Boris, Vay, and Higuera-Cary schemes. We also present a new relativistic fully implicit particle integrator that is energy conserving. Furthermore, a method based on the relativistic guiding center approximation is included. The algorithms are described such that they can be readily implemented in magnetohydrodynamics codes or Particle-in-Cell codes. Our comparison focuses on the strengths and key features of the particle integrators. We test the conservation of invariants of motion, and the accuracy of particle drift dynamics in highly relativistic, mildly relativistic, and non-relativistic settings. The methods are compared in idealized test cases, i.e., without considering feedback on the electrodynamic fields, collisions, pair creation, or radiation. The test cases include uniform electric and magnetic fields, $\mathbf{E}\times \mathbf{B}$-fields, force-free fields, and setups relevant for high-energy astrophysics, e.g., a magnetic mirror, a magnetic dipole, and a magnetic null. These tests have direct relevance for particle acceleration in shocks and in magnetic reconnection.

Matching published version in ApJS

References in corpus (14)

• The Black Hole Accretion Code
• MPI-AMRVAC for Solar and Astrophysics
• MPI-AMRVAC 2.0 for Solar and Astrophysical Applications
• Divergence-Free Reconstruction of Magnetic Fields and WENO Schemes for Magnetohydrodynamics
• Magnetohydrodynamic-Particle-in-Cell Method for Coupling Cosmic Rays with a Thermal Plasma: Application to Non-relativistic Shocks
• Structure-preserving second-order integration of relativistic charged particle trajectories in electromagnetic fields
• Nonthermally Dominated Electron Acceleration during Magnetic Reconnection in a Low-beta Plasma
• Diffusion in Pulsar Wind Nebulae: An Investigation using Magnetohydrodynamic and Particle Transport Models
• Reconnection and particle acceleration in interacting flux ropes -- II. 3D effects on test particles in magnetically dominated plasmas
• Reconnection and particle acceleration in interacting flux ropes I. Magnetohydrodynamics and test particles in 2.5D
• Comment on "Symplectic integration of magnetic systems": a proof that the Boris algorithm is not variational
• An fully implicit scheme for numerical integration of the relativistic particle equation of motion
• An Efficient Second-Order Accurate and Continuous Interpolation for Block-Adaptive Grids
• High Order Numerical Integrators for Relativistic Charged Particle Tracking

Cited by in corpus (48)

• Magnetic Reconnection and Hot Spot Formation in Black Hole Accretion Disks
• MPI-AMRVAC 3.0: updates to an open-source simulation framework
• Afivo: a framework for quadtree/octree AMR with shared-memory parallelization and geometric multigrid methods
• General relativistic resistive magnetohydrodynamics with robust primitive variable recovery for accretion disk simulations
• On the Boris solver in particle-in-cell simulation
• Novel aspects of cosmic ray diffusion in synthetic magnetic turbulence
• PIC methods in astrophysics: Simulations of relativistic jets and kinetic physics in astrophysical systems
• PIC Simulation Methods for Cosmic Radiation and Plasma Instabilities
• Generalized, energy-conserving numerical simulations of particles in general relativity. I. Time-like and null geodesics
• Generalized, energy-conserving numerical simulations of particles in general relativity. II. Test particles in electromagnetic fields and GRMHD
• Relativistic Particle Transport and Acceleration in Structured Plasma Turbulence
• Structure-preserving numerical simulations of test particle dynamics around slowly rotating neutron stars within Hartle-Thorne approach
• Fast acceleration of transrelativistic electrons in astrophysical turbulence
• Weak Alfvénic turbulence in relativistic plasmas II: Current sheets and dissipation
• A coupled guiding center-Boris particle pusher for magnetized plasmas in compact-object magnetospheres
• Optimisation of Thin Plastic Foil Targets for Production of Laser-Generated Protons in the GeV Range
• Modeling interactions of narrowband large amplitude whistler-mode waves with electrons in the solar wind inside ~.3 AU and at 1 AU using a particle tracing code
• Turbulence and particle energization in twisted flux ropes in solar-wind conditions
• A relativistic particle pusher for ultra-strong electromagnetic fields
• Towards Synthetic Magnetic Turbulence with Coherent Structures
• Cosmic rays and random magnetic traps
• A guiding center implementation for relativistic particle dynamics in the PLUTO code
• Stochastic diffusion of electrons interacting with whistler-mode waves in the solar wind
• OSIRIS-GR: General relativistic activation of the polar cap of a compact neutron star
• Magnetic Island Merging: Two-dimensional MHD Simulation and Test-Particle Modeling
• High Order Explicit Lorentz Invariant Volume-preserving Algorithms for Relativistic Dynamics of Charged Particles
• Particle acceleration in neutron star ultra-strong electromagnetic fields
• An Analytic Boris Pusher for Plasma Simulation
• Particle orbits at the magnetopause: Kelvin-Helmholtz induced trapping
• Towards Laser Ion Acceleration With Holed Targets
• Toward Realistic Solar Flare Models: An explicit Particle-In-Cell solver in the DISPATCH framework
• An Evaluation of Different Numerical Methods to Calculate the Pitch-angle Diffusion Coefficient from Full-orbit Simulations: disentangling a rope of sand
• Diffusion of relativistic charged particles and field lines in isotropic turbulence: I. Numerical simulations
• A Parallel-Kinetic-Perpendicular-Moment Model for Magnetized Plasmas
• Nonthermal electron acceleration in turbulent post-flare coronal loops
• A fundamental derivation of two Boris solvers and the Ge-Marsden theorem
• Test particle sampling and particle acceleration in a 2D coronal plasmoid-mediated reconnecting current sheet
• Synthetic Turbulence via an Instanton Gas Approximation
• Explicit volume-preserving numerical schemes for relativistic trajectories and spin dynamics
• A New Particle Pusher with Hadronic Interactions for Modeling Multimessenger Emission from Compact Objects
• A new field solver for modeling of relativistic particle-laser interactions using the particle-in-cell algorithm
• Modelling cosmic-ray transport: magnetised versus unmagnetised motion in astrophysical magnetic turbulence
• Towards Modelling AR Sco: Generalised Particle Dynamics and Strong Radiation-Reaction Regimes
• Numerical Investigation of Efficient Electron Acceleration at an Unsteady Solar Flare Loop-Top
• Fully kinetic shearing-box simulations of magnetorotational turbulence in 2D and 3D. I. Pair plasmas
• An Electromagnetic Particle-Particle Model on Solving Relativistic Binary Collision
• High Performance Implementation of Boris Particle Pusher on DPC++. A First Look at oneAPI
• High-energy Particle Transport in Three-dimensional Anisotropic Turbulent Magnetic Fields