A Comprehensive Comparison of Relativistic Particle Integrators
- KU Leuven, Leuven (Belgium)
- Institut fur Theoretische Physik, Frankfurt (Germany)
- Columbia Univ., New York, NY (United States)
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. Here, 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 onto the electrodynamic fields, collisions, pair creation, or radiation. The test cases include uniform electric and magnetic fields, $${\boldsymbol{E}}\times {\boldsymbol{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.
- Research Organization:
- Columbia Univ., New York, NY (United States)
- Sponsoring Organization:
- USDOE Office of Science (SC)
- Grant/Contract Number:
- SC0016542
- OSTI ID:
- 1542033
- Journal Information:
- The Astrophysical Journal. Supplement Series (Online), Vol. 235, Issue 1; ISSN 1538-4365
- Publisher:
- American Astronomical Society/IOPCopyright Statement
- Country of Publication:
- United States
- Language:
- English
Web of Science
Similar Records
Plasmoid statistics in relativistic magnetic reconnection
Particle Acceleration during Magnetic Reconnection in a Low-beta Plasma