Motion equations in a Kerr-Newman-de Sitter spacetime: some methods of integration and application to black holes shadowing in Scilab

Abstract

In this paper, we recall some basic facts about the Kerr--Newman--(anti) de Sitter (KNdS) spacetime and review several formulations and integration methods for the geodesic equation of a test particle in such a spacetime. In particular, we introduce some basic general symplectic integrators in the Hamiltonian formalism and we re-derive the separated motion equations using Carter's method. After this theoretical background, we explain how to ray-trace a KNdS black hole, equipped with a thin accretion disk, using Scilab. We compare the accuracy and execution time of the previous methods, concluding that the Carter equations is the best one. Then, inspired by Hagihara, we apply Weierstrass' elliptic functions to the non-rotating case, yielding a fairly fast shadowing program for such a spacetime. We provide some illustrations of the code, including a depiction of the effects of the cosmological constant on shadows and accretion disk, as well as a simulation of M87*.

References in corpus (8)

• First M87 Event Horizon Telescope Results. I. The Shadow of the Supermassive Black Hole
• A new method for shadow calculations: application to parameterised axisymmetric black holes
• A Periodic Table for Black Hole Orbits
• The Kerr Metric
• Light Curves from an MHD Simulation of a Black Hole Accretion Disk
• Chaotic Shadows of Black Holes: A Short Review
• Testing General Relativity with High-Resolution Imaging of Sgr A*
• Curvature Invariants for accelerating Kerr-Newman black holes in (anti-)de Sitter spacetime