Generalized, energy-conserving numerical simulations of particles in general relativity. II. Test particles in electromagnetic fields and GRMHD

Abstract

Direct observations of compact objects, in the form of radiation spectra, gravitational waves from VIRGO/LIGO, and forthcoming direct imaging, are currently one of the primary source of information on the physics of plasmas in extreme astrophysical environments. The modeling of such physical phenomena requires numerical methods that allow for the simulation of microscopic plasma dynamics in presence of both strong gravity and electromagnetic fields. In Bacchini et al. (2018) we presented a detailed study on 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).

References in corpus (57)

• Event-horizon-scale structure in the supermassive black hole candidate at the Galactic Centre
• HARM: A Numerical Scheme for General Relativistic Magnetohydrodynamics
• Relativistic Reconnection: an Efficient Source of Non-Thermal Particles
• Formation of Hard Power-laws in the Energetic Particle Spectra Resulting from Relativistic Magnetic Reconnection
• Three-Dimensional General Relativistic Radiation Magnetohydrodynamical Simulation of Super-Eddington Accretion, using a new code HARMRAD with M1 Closure
• Acceleration of Particles at the Termination Shock of a Relativistic Striped Wind
• The Black Hole Accretion Code
• Simulations of particle acceleration beyond the classical synchrotron burnoff limit in magnetic reconnection: An explanation of the Crab flares
• Very High Order $\PNM$ Schemes on Unstructured Meshes for the Resistive Relativistic MHD Equations
• Resolved Magnetic-Field Structure and Variability Near the Event Horizon of Sagittarius A*
• Non-thermal particle acceleration in collisionless relativistic electron-proton reconnection
• The Energy Conserving Particle-in-Cell Method
• A public code for general relativistic, polarised radiative transfer around spinning black holes
• First-Principles Plasma Simulations of Black-Hole Jet Launching
• Circular orbits and related quasi-harmonic oscillatory motion of charged particles around weakly magnetized rotating black holes
• Acceleration of charged particles due to chaotic scattering in the combined black hole gravitational field and asymptotically uniform magnetic field
• Simulations of relativistic collisionless shocks: shock structure and particle acceleration
• Motion of charged particles around a rotating black hole in a magnetic field
• BlackHoleCam: fundamental physics of the Galactic center
• IllinoisGRMHD: An Open-Source, User-Friendly GRMHD Code for Dynamical Spacetimes
• The role of electron heating physics in images and variability of the Galactic Centre black hole Sagittarius A*
• The Disc-Jet Symbiosis Emerges: Modeling the Emission of Sagittarius A* with Electron Thermodynamics
• Radiative, two-temperature simulations of low luminosity black hole accretion flows in general relativity
• A fully covariant mean-field dynamo closure for numerical 3+1 resistive GRMHD
• Exactly Energy Conserving Semi-Implicit Particle in Cell Formulation
• Coupled jet-disk model for Sgr A*: explaining the flat-spectrum radio core with GRMHD simulations of jets
• Electron and proton heating in trans-relativistic magnetic reconnection
• Nonthermal Particle Acceleration in 3D Relativistic Magnetic Reconnection in Pair Plasma
• General-relativistic resistive-magnetohydrodynamic simulations of binary neutron stars
• Transition from Regular to Chaotic Circulation in Magnetized Coronae near Compact Objects
• Synchrotron radiation of self-collimating relativistic MHD jets
• Explicit symplectic approximation of nonseparable Hamiltonians: algorithm and long time performance
• A comprehensive comparison of relativistic particle integrators
• Structure-preserving second-order integration of relativistic charged particle trajectories in electromagnetic fields
• Global three-dimensional simulation of Earth's dayside reconnection using a two-way coupled magnetohydrodynamics with embedded particle-in-cell model: initial results
• General relativistic magnetohydrodynamical $\kappa$-jet models for Sgr A*
• Chaotic motion of Charged Particles in an Electromagnetic Field Surrounding a Rotating Black Hole
• Charged particle motion in Kerr-Newmann space-times
• Ab-initio pulsar magnetosphere: the role of general relativity
• RAPTOR I: Time-dependent radiative transfer in arbitrary spacetimes
• Modeling magnetized neutron stars using resistive MHD
• Particle-in-cell simulations of pair discharges in a starved magnetosphere of a Kerr black hole
• The Radiative Efficiency and Spectra of Slowly Accreting Black Holes from Two-Temperature GRRMHD Simulations
• Inducing chaos by breaking axial symmetry in a black hole magnetosphere
• The Properties of Reconnection Current Sheets in GRMHD Simulations of Radiatively Inefficient Accretion Flows
• rHARM: Accretion and Ejection in Resistive GR-MHD
• Energy conserving schemes for the simulation of musical instrument contact dynamics
• Generalized, energy-conserving numerical simulations of particles in general relativity. I. Time-like and null geodesics
• Near-horizon structure of escape zones of electrically charged particles around weakly magnetized rotating black hole
• The Fluid-Kinetic Particle-in-Cell Solver for Plasma Simulations
• 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
• 3+1 Magnetodynamics
• Multiple scale kinetic simulations with the energy conserving semi implicit particle in cell (PIC) method
• GRay2: A General Purpose Geodesic Integrator for Kerr Spacetimes
• PolyPIC: the Polymorphic-Particle-in-Cell Method for Fluid-Kinetic Coupling
• Test particles in relativistic resistive magnetohydrodynamics

Cited by in corpus (28)

• The Event Horizon General Relativistic Magnetohydrodynamic Code Comparison Project
• Magnetic Reconnection and Hot Spot Formation in Black Hole Accretion Disks
• MPI-AMRVAC 3.0: updates to an open-source simulation framework
• Plasmoid formation in global GRMHD simulations and AGN flares
• General relativistic resistive magnetohydrodynamics with robust primitive variable recovery for accretion disk simulations
• Construction of Explicit Symplectic Integrators in General Relativity. I. Schwarzschild Black Holes
• Construction of explicit symplectic integrators in general relativity. II. Reissner-Nordstrom black holes
• Charged particle dynamics in parabolic magnetosphere around Schwarzschild black hole
• Synthetic gamma-ray lightcurves of Kerr black-hole magnetospheric activity from particle-in-cell simulations
• Construction of explicit symplectic integrators in general relativity. IV. Kerr black holes
• General relativistic MHD simulations of non-thermal flaring in Sagittarius A*
• Acceleration and Escape Processes of High-energy Particles in Turbulence inside Hot Accretion Flows
• Construction of explicit symplectic integrators in general relativity. III. Reissner-Nordstrom-(anti)-de Sitter black holes
• Radiation Transport Two-Temperature GRMHD Simulations of Warped Accretion Disks
• Spinning black holes magnetically connected to a Keplerian disk -- Magnetosphere, reconnection sheet, particle acceleration and coronal heating
• Structure-preserving numerical simulations of test particle dynamics around slowly rotating neutron stars within Hartle-Thorne approach
• A coupled guiding center-Boris particle pusher for magnetized plasmas in compact-object magnetospheres
• VHE Emission from Magnetic Reconnection in the RIAF of SgrA*
• Particle Diffusion and Acceleration in Magnetorotational Instability Turbulence
• A novel energy-conserving scheme for eight-dimensional systems
• The use of Kepler solver in numerical integrations of quasi-Keplerian orbits
• Energy-conserving integrator for conservative Hamiltonian systems with ten-dimensional phase space
• Explicit symplectic integrators with adaptive time steps in curved spacetimes
• Impact of electric charges on chaos in magnetized $Reissner-Nordstr\ddot{o}m$ spacetimes
• A New Particle Pusher with Hadronic Interactions for Modeling Multimessenger Emission from Compact Objects
• Construction of second-order six-dimensional Hamiltonian-conserving scheme
• A new correction method for quasi-Keplerian orbits
• FPIC: a new Particle-In-Cell code for stationary and axisymmetric black-hole spacetimes