DOE PAGES 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

Abstract

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,more » 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.« less

Authors:
ORCiD logo [1]; ORCiD logo [1]; ORCiD logo [2];  [3]
  1. KU Leuven (Belgium)
  2. Princeton Univ., NJ (United States)
  3. Columbia Univ., New York, NY (United States)
Publication Date:
Research Org.:
Columbia Univ., New York, NY (United States)
Sponsoring Org.:
USDOE Office of Science (SC), Fusion Energy Sciences (FES)
OSTI Identifier:
1542035
Grant/Contract Number:  
SC0016542
Resource Type:
Accepted Manuscript
Journal Name:
The Astrophysical Journal. Supplement Series (Online)
Additional Journal Information:
Journal Name: The Astrophysical Journal. Supplement Series (Online); Journal Volume: 237; Journal Issue: 1; Journal ID: ISSN 1538-4365
Publisher:
American Astronomical Society/IOP
Country of Publication:
United States
Language:
English
Subject:
79 ASTRONOMY AND ASTROPHYSICS; gravitation; methods: numerical; relativistic processes; stars: black holes; stars: imaging

Citation Formats

Bacchini, F., Ripperda, B., Chen, A. Y., and Sironi, L. Generalized, Energy-conserving Numerical Simulations of Particles in General Relativity. I. Time-like and Null Geodesics. United States: N. p., 2018. Web. doi:10.3847/1538-4365/aac9ca.
Bacchini, F., Ripperda, B., Chen, A. Y., & Sironi, L. Generalized, Energy-conserving Numerical Simulations of Particles in General Relativity. I. Time-like and Null Geodesics. United States. https://doi.org/10.3847/1538-4365/aac9ca
Bacchini, F., Ripperda, B., Chen, A. Y., and Sironi, L. Tue . "Generalized, Energy-conserving Numerical Simulations of Particles in General Relativity. I. Time-like and Null Geodesics". United States. https://doi.org/10.3847/1538-4365/aac9ca. https://www.osti.gov/servlets/purl/1542035.
@article{osti_1542035,
title = {Generalized, Energy-conserving Numerical Simulations of Particles in General Relativity. I. Time-like and Null Geodesics},
author = {Bacchini, F. and Ripperda, B. and Chen, A. Y. and Sironi, L.},
abstractNote = {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.},
doi = {10.3847/1538-4365/aac9ca},
journal = {The Astrophysical Journal. Supplement Series (Online)},
number = 1,
volume = 237,
place = {United States},
year = {Tue Jul 03 00:00:00 EDT 2018},
month = {Tue Jul 03 00:00:00 EDT 2018}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

Citation Metrics:
Cited by: 29 works
Citation information provided by
Web of Science

Save / Share:

Works referenced in this record:

Three-dimensional general relativistic radiation magnetohydrodynamical simulation of super-Eddington accretion, using a new code harmrad with M1 closure
journal, May 2014

  • McKinney, Jonathan C.; Tchekhovskoy, Alexander; Sadowski, Aleksander
  • Monthly Notices of the Royal Astronomical Society, Vol. 441, Issue 4
  • DOI: 10.1093/mnras/stu762

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

  • Morris, Michael S.; Thorne, Kip S.
  • American Journal of Physics, Vol. 56, Issue 5
  • DOI: 10.1119/1.15620

3+1 Formalism in General Relativity
book, February 2012


A Comprehensive Comparison of Relativistic Particle Integrators
journal, March 2018

  • Ripperda, B.; Bacchini, F.; Teunissen, J.
  • The Astrophysical Journal Supplement Series, Vol. 235, Issue 1
  • DOI: 10.3847/1538-4365/aab114

Understanding Modern Molecular Dynamics:  Techniques and Applications
journal, January 2000

  • Tuckerman, Mark E.; Martyna, Glenn J.
  • The Journal of Physical Chemistry B, Vol. 104, Issue 2
  • DOI: 10.1021/jp992433y

Spectra of black hole accretion models of ultraluminous X-ray sources
journal, April 2017

  • Narayan, Ramesh; Sa̧dowski, Aleksander; Soria, Roberto
  • Monthly Notices of the Royal Astronomical Society, Vol. 469, Issue 3
  • DOI: 10.1093/mnras/stx1027

Simulations of relativistic collisionless shocks: shock structure and particle acceleration
conference, January 2005

  • Spitkovsky, Anatoly
  • ASTROPHYSICAL SOURCES OF HIGH ENERGY PARTICLES AND RADIATION, AIP Conference Proceedings
  • DOI: 10.1063/1.2141897

RAPTOR: I. Time-dependent radiative transfer in arbitrary spacetimes★
journal, May 2018


ipole – semi-analytic scheme for relativistic polarized radiative transport
journal, December 2017

  • Mościbrodzka, M.; Gammie, C. F.
  • Monthly Notices of the Royal Astronomical Society, Vol. 475, Issue 1
  • DOI: 10.1093/mnras/stx3162

Viewing the Shadow of the Black Hole at the Galactic Center
journal, January 2000

  • Falcke, Heino; Melia, Fulvio; Agol, Eric
  • The Astrophysical Journal, Vol. 528, Issue 1
  • DOI: 10.1086/312423

General relativistic radiative transfer code in rotating black hole space–time: ARTIST
journal, October 2016

  • Takahashi, Rohta; Umemura, Masayuki
  • Monthly Notices of the Royal Astronomical Society, Vol. 464, Issue 4
  • DOI: 10.1093/mnras/stw2479

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

  • Ripperda, B.; Porth, O.; Xia, C.
  • Monthly Notices of the Royal Astronomical Society
  • DOI: 10.1093/mnras/stx379

Higher order explicit symmetric integrators for inseparable forms of coordinates and momenta
journal, April 2016

  • Liu, Lei; Wu, Xin; Huang, Guoqing
  • Monthly Notices of the Royal Astronomical Society, Vol. 459, Issue 2
  • DOI: 10.1093/mnras/stw807

The black hole accretion code
journal, May 2017

  • Porth, Oliver; Olivares, Hector; Mizuno, Yosuke
  • Computational Astrophysics and Cosmology, Vol. 4, Issue 1
  • DOI: 10.1186/s40668-017-0020-2

An Extension of the Athena++ code Framework for Grmhd Based on Advanced Riemann Solvers and Staggered-Mesh Constrained Transport
journal, August 2016

  • White, Christopher J.; Stone, James M.; Gammie, Charles F.
  • The Astrophysical Journal Supplement Series, Vol. 225, Issue 2
  • DOI: 10.3847/0067-0049/225/2/22

Particle acceleration and energy conservation in particle in cell simulations
journal, July 2011

  • Lapenta, Giovanni; Markidis, Stefano
  • Physics of Plasmas, Vol. 18, Issue 7
  • DOI: 10.1063/1.3602216

Imaging the Schwarzschild-radius-scale Structure of M87 with the Event Horizon Telescope Using Sparse Modeling
journal, March 2017

  • Akiyama, Kazunori; Kuramochi, Kazuki; Ikeda, Shiro
  • The Astrophysical Journal, Vol. 838, Issue 1
  • DOI: 10.3847/1538-4357/aa6305

Why is Boris algorithm so good?
journal, August 2013

  • Qin, Hong; Zhang, Shuangxi; Xiao, Jianyuan
  • Physics of Plasmas, Vol. 20, Issue 8
  • DOI: 10.1063/1.4818428

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

  • Higuera, A. V.; Cary, J. R.
  • Physics of Plasmas, Vol. 24, Issue 5
  • DOI: 10.1063/1.4979989

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

  • Tuckerman, Mark E.; Martyna, Glenn J.
  • The Journal of Physical Chemistry B, Vol. 105, Issue 31
  • DOI: 10.1021/jp011879n

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


A new method for shadow calculations: application to parameterised axisymmetric black holes
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


Modelling the black hole silhouette in Sagittarius A* with ion tori
journal, July 2012


Studying null and time-like geodesics in the classroom
journal, March 2011


Observable Emission Features of Black Hole GRMHD Jets on Event Horizon Scales
journal, August 2017


A comprehensive comparison of relativistic particle integrators
text, January 2017


Works referencing / citing this record:

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

  • Bacchini, F.; Ripperda, B.; Porth, O.
  • The Astrophysical Journal Supplement Series, Vol. 240, Issue 2
  • DOI: 10.3847/1538-4365/aafcb3

General-relativistic Resistive Magnetohydrodynamics with Robust Primitive-variable Recovery for Accretion Disk Simulations
journal, September 2019

  • Ripperda, B.; Bacchini, F.; Porth, O.
  • The Astrophysical Journal Supplement Series, Vol. 244, Issue 1
  • DOI: 10.3847/1538-4365/ab3922

The Event Horizon General Relativistic Magnetohydrodynamic Code Comparison Project
journal, August 2019

  • Porth, Oliver; Chatterjee, Koushik; Narayan, Ramesh
  • The Astrophysical Journal Supplement Series, Vol. 243, Issue 2
  • DOI: 10.3847/1538-4365/ab29fd

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

  • Porth, Oliver; Chatterjee, Koushik; Narayan, Ramesh
  • Maryland Shared Open Access Repository
  • DOI: 10.13016/m2htjg-wbul

The Event Horizon General Relativistic Magnetohydrodynamic Code Comparison Project
text, January 2019


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

  • Bacchini, F.; Ripperda, B.; Porth, O.
  • The Astrophysical Journal Supplement Series, Vol. 240, Issue 2
  • DOI: 10.3847/1538-4365/aafcb3

Numerical methods for General Relativistic particles
text, January 2018