FANTASY: User-Friendly Symplectic Geodesic Integrator for Arbitrary Metrics with Automatic Differentiation

Abstract

We present FANTASY (Finally A Numerical Trajectory Algorithm both Straightforward and sYmplectic), a user-friendly, open-source symplectic geodesic integrator written in Python. FANTASY is designed to work "out-of-the-box" and does not require anything from the user aside from the metric and the initial conditions for the geodesics. FANTASY efficiently computes derivatives up to machine precision using automatic differentiation, allowing the integration of geodesics in arbitrary space(times) without the need for the user to manually input Christoffel symbols or any other metric derivatives. Further, FANTASY utilizes a Hamiltonian integration scheme that doubles the phase space, where two copies of the particle phase space are evolved together. This technique allows for an integration scheme that is both explicit and symplectic, even when the Hamiltonian is not separable. FANTASY comes prebuilt with second and fourth order schemes, and is easily extendible to higher order schemes. FANTASY also includes an automatic Jacobian calculator that allows for coordinate transformations to be done automatically.

Minor typos corrected, matches accepted version

References in corpus (12)

• First M87 Event Horizon Telescope Results. I. The Shadow of the Supermassive Black Hole
• First M87 Event Horizon Telescope Results. VI. The Shadow and Mass of the Central Black Hole
• First M87 Event Horizon Telescope Results. IV. Imaging the Central Supermassive Black Hole
• First M87 Event Horizon Telescope Results. V. Physical Origin of the Asymmetric Ring
• First M87 Event Horizon Telescope Results. II. Array and Instrumentation
• First M87 Event Horizon Telescope Results. III. Data Processing and Calibration
• A Metric for Rapidly Spinning Black Holes Suitable for Strong-Field Tests of the No-Hair Theorem
• Dynamical Chern-Simons Modified Gravity I: Spinning Black Holes in the Slow-Rotation Approximation
• A new method for shadow calculations: application to parameterised axisymmetric black holes
• Bumpy Black Holes in Alternate Theories of Gravity
• The existence of relativistic stars in f(R) gravity
• Algorithms And Programs For Strong Gravitational Lensing In Kerr Space-time Including Polarization

Cited by in corpus (5)

• Construction of explicit symplectic integrators in general relativity. IV. Kerr black holes
• Not all spacetime coordinates for general-relativistic ray tracing are created equal
• Dynamical Lensing Tomography of Black Hole Ringdown
• Relativistic location algorithm in curved spacetime
• Construction of second-order six-dimensional Hamiltonian-conserving scheme