Generalized, Energy-conserving Numerical Simulations of Particles in General Relativity II. Test Particles in Electromagnetic Fields and GRMHD

Observations of compact objects, in the form of radiation spectra, gravitational waves by LIGO/Virgo, and direct imaging with the Event Horizon Telescope, are currently the main information sources on plasma physics in extreme gravity. Modeling such physical phenomena Requires numerical methods that allow for the simulation of microscopic plasma dynamics in the presence of both strong gravity and electromagnetic fields. In Bacchini et al. (2018), we presented a detailed study of numerical techniques for the integration of free geodesic motion. Here, we extend the study by introducing electromagnetic forces in the simulation of charged particles in curved spacetimes. We extend the Hamiltonian energy-conserving method presented in Bacchini et al. (2018) to include the Lorentz force, and we test its performance compared to that of standard explicit Runge–Kutta and implicit midpoint rule schemes against analytic solutions. Then, we show the application of the numerical schemes to the integration of test particle trajectories in general relativistic magnetohydrodynamic (GRMHD) simulations by modifying the algorithms to handle grid-based electromagnetic fields. We test this approach by simulating ensembles of charged particles in a static GRMHD configuration obtained with the black hole accretion code (BHAC).