skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Generalized, Energy-conserving Numerical Simulations of Particles in General Relativity. I. Time-like and Null Geodesics

Journal Article · · The Astrophysical Journal. Supplement Series (Online)

The numerical integration of particle trajectories in curved spacetimes is fundamental for obtaining realistic models of the particle dynamics around massive compact objects such as black holes and neutron stars. Generalized algorithms capable of handling generic metrics are needed for studies of both standard (Schwarzschild and Kerr metrics) and nonstandard (e.g., Schwarzschild metric plus nonclassical perturbations or multiple black hole metrics) spacetimes. The most commonly employed explicit numerical schemes (e.g., Runge–Kutta) are incapable of producing highly accurate results at critical points, e.g., in the regions close to the event horizon where gravity causes extreme curvature of the spacetime, at an acceptable computational cost. In this work we describe a generalized algorithm for the numerical integration of time-like (massive particles) and null (photons) geodesics in any given 3 + 1 split spacetime. We introduce a new, exactly energy-conserving implicit integration scheme based on the preservation of the underlying Hamiltonian, and we compare its properties with a standard fourth-order Runge–Kutta explicit scheme and an implicit midpoint scheme. We test the numerical performance of the three schemes against analytic solutions of particle and photon orbits in Schwarzschild and Kerr spacetimes. We additionally prove the versatility of our framework in handling more exotic metrics, such as Morris–Thorne wormholes and quantum-perturbed Schwarzschild black holes. The generalized method is also discussed in the perspective of future extensions to more complex particle dynamics, e.g., the addition of the Lorentz force acting on charged particles, which allows for test particle diagnostics in GRMHD simulations.

Research Organization:
Columbia Univ., New York, NY (United States)
Sponsoring Organization:
USDOE Office of Science (SC), Fusion Energy Sciences (FES)
Grant/Contract Number:
SC0016542
OSTI ID:
1542035
Journal Information:
The Astrophysical Journal. Supplement Series (Online), Vol. 237, Issue 1; ISSN 1538-4365
Publisher:
American Astronomical Society/IOPCopyright Statement
Country of Publication:
United States
Language:
English
Citation Metrics:
Cited by: 29 works
Citation information provided by
Web of Science

References (64)

Three-dimensional general relativistic radiation magnetohydrodynamical simulation of super-Eddington accretion, using a new code harmrad with M1 closure journal May 2014
GeodesicViewer – A tool for exploring geodesics in the theory of relativity journal February 2010
Wormholes in spacetime and their use for interstellar travel: A tool for teaching general relativity journal May 1988
3+1 Formalism in General Relativity book February 2012
A Comprehensive Comparison of Relativistic Particle Integrators journal March 2018
Understanding Modern Molecular Dynamics:  Techniques and Applications journal January 2000
Spectra of black hole accretion models of ultraluminous X-ray sources journal April 2017
Simulations of relativistic collisionless shocks: shock structure and particle acceleration conference January 2005
RAPTOR: I. Time-dependent radiative transfer in arbitrary spacetimes★ journal May 2018
ipole – semi-analytic scheme for relativistic polarized radiative transport journal December 2017
Viewing the Shadow of the Black Hole at the Galactic Center journal January 2000
General relativistic radiative transfer code in rotating black hole space–time: ARTIST journal October 2016
Beyond the event horizon or altogether without it? journal March 2017
GYOTO: a new general relativistic ray-tracing code journal October 2011
The Mathematical Theory of Black Holes book January 1984
Modelling the black hole silhouette in Sagittarius A* with ion tori journal July 2012
Reconnection and particle acceleration in interacting flux ropes I. Magnetohydrodynamics and test particles in 2.5D journal February 2017
Higher order explicit symmetric integrators for inseparable forms of coordinates and momenta journal April 2016
The black hole accretion code journal May 2017
An Extension of the Athena++ code Framework for Grmhd Based on Advanced Riemann Solvers and Staggered-Mesh Constrained Transport journal August 2016
Particle acceleration and energy conservation in particle in cell simulations journal July 2011
Imaging the Schwarzschild-radius-scale Structure of M87 with the Event Horizon Telescope Using Sparse Modeling journal March 2017
Why is Boris algorithm so good? journal August 2013
Jet-Launching Structure Resolved Near the Supermassive Black Hole in M87 journal September 2012
Energy conserving schemes for the simulation of musical instrument contact dynamics journal March 2015
Studying null and time-like geodesics in the classroom journal March 2011
Structure-preserving second-order integration of relativistic charged particle trajectories in electromagnetic fields journal May 2017
Observable Emission Features of Black Hole GRMHD Jets on Event Horizon Scales journal August 2017
Escape of photons from two fixed extreme Reissner-Nordström black holes journal November 2008
A Ray-Tracing Algorithm for Spinning Compact Object Spacetimes with Arbitrary Quadrupole Moments. ii. Neutron Stars journal June 2012
An Improved Distance and mass Estimate for sgr a* from a Multistar Orbit Analysis journal October 2016
Photon Rings Around kerr and Kerr-Like Black Holes journal October 2013
Testing the No-Hair Theorem with Observations in the Electromagnetic Spectrum. i. Properties of a Quasi-Kerr Spacetime journal May 2010
Testing the No-Hair Theorem with Observations in the Electromagnetic Spectrum. iii. Quasi-Periodic Variability journal December 2010
Gravitational Field of a Spinning Mass as an Example of Algebraically Special Metrics journal September 1963
A periodic table for black hole orbits journal May 2008
Einstein rings as a tool for estimating distances and the mass of a Schwarzschild black hole journal June 2008
A Ray-Tracing Algorithm for Spinning Compact Object Spacetimes with Arbitrary Quadrupole Moments. i. Quasi-Kerr Black Holes journal December 2011
A QUANTITATIVE TEST OF THE NO-HAIR THEOREM WITH Sgr A* USING STARS, PULSARS, AND THE EVENT HORIZON TELESCOPE journal February 2016
ODYSSEY : A PUBLIC GPU-BASED CODE FOR GENERAL RELATIVISTIC RADIATIVE TRANSFER IN KERR SPACETIME journal March 2016
Stability of a Schwarzschild Singularity journal November 1957
Symmetric integrator for nonintegrable Hamiltonian relativistic systems journal December 2012
Spherical Photon Orbits Around a Kerr Black Hole journal November 2003
Astrophysical imaging of Kerr black holes with scalar hair journal October 2016
Shadows of Bonnor black dihole by chaotic lensing journal March 2018
New method for shadow calculations: Application to parametrized axisymmetric black holes journal October 2016
Understanding Modern Molecular Dynamics:  Techniques and Applications journal July 2001
Studying null and time-like geodesics in the classroom text January 2011
GYOTO: a new general relativistic ray-tracing code text January 2011
A Symmetric Integrator for non-integrable Hamiltonian Relativistic Systems text January 2012
Jet Launching Structure Resolved Near the Supermassive Black Hole in M87 text January 2012
Energy conserving schemes for the simulation of musical instrument contact dynamics text January 2015
Photon Rings around Kerr and Kerr-like Black Holes text January 2015
A Quantitative Test of the No-Hair Theorem with Sgr A* using stars, pulsars, and the Event Horizon Telescope text January 2015
An Extension of the Athena++ Code Framework for GRMHD Based on Advanced Riemann Solvers and Staggered-Mesh Constrained Transport text January 2015
Odyssey: A Public GPU-Based Code for General-Relativistic Radiative Transfer in Kerr Spacetime text January 2016
A new method for shadow calculations: application to parameterised axisymmetric black holes text January 2016
Reconnection and particle acceleration in interacting flux ropes I. Magnetohydrodynamics and test particles in 2.5D text January 2016
Observable Emission Features of Black Hole GRMHD Jets on Event Horizon Scales text January 2017
Shadows of Bonnor black dihole by chaotic lensing text January 2017
RAPTOR I: Time-dependent radiative transfer in arbitrary spacetimes text January 2018
Viewing the Shadow of the Black Hole at the Galactic Center text January 1999
Imaging the Schwarzschild-radius-scale Structure of M87 with the Event Horizon Telescope using Sparse Modeling text January 2017
A comprehensive comparison of relativistic particle integrators text January 2017

Cited By (9)

Generalized, Energy-conserving Numerical Simulations of Particles in General Relativity. II. Test Particles in Electromagnetic Fields and GRMHD journal February 2019
General-relativistic Resistive Magnetohydrodynamics with Robust Primitive-variable Recovery for Accretion Disk Simulations journal September 2019
The Event Horizon General Relativistic Magnetohydrodynamic Code Comparison Project journal August 2019
A Novel Energy-conserving Scheme for Eight-dimensional Hamiltonian Problems journal December 2019
The Event Horizon General Relativistic Magnetohydrodynamic Code Comparison Project collection January 2019
Generalized, energy-conserving numerical simulations of particles in general relativity. II. Test particles in electromagnetic fields and GRMHD text January 2018
The Event Horizon General Relativistic Magnetohydrodynamic Code Comparison Project text January 2019
Numerical methods for General Relativistic particles text January 2018
Structure preserving algorithms for simulation of linearly damped acoustic systems text January 2020