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.

arXiv commentMinor typos corrected, matches accepted version

References in corpus (32)

• 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
• Black Holes in Higher Dimensions
• Rotating Charged Black Hole Solution in Heterotic String Theory
• Measuring thermodynamic length
• Testing the No-Hair Theorem with Observations in the Electromagnetic Spectrum: II. Black-Hole Images
• 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
• Black Saturn
• General parametrization of axisymmetric black holes in metric theories of gravity
• A Fast New Public Code for Computing Photon Orbits in a Kerr Spacetime
• A public code for general relativistic, polarised radiative transfer around spinning black holes
• Neutron Stars in f(R) Gravity with Perturbative Constraints
• Bumpy Black Holes in Alternate Theories of Gravity
• Rotating black hole in extended Chern-Simons modified gravity
• The existence of relativistic stars in f(R) gravity
• A Ray-Tracing Algorithm for Spinning Compact Object Spacetimes with Arbitrary Quadrupole Moments. I. Quasi-Kerr Black Holes
• GRay: a Massively Parallel GPU-Based Code for Ray Tracing in Relativistic Spacetimes
• Explicit symplectic approximation of nonseparable Hamiltonians: algorithm and long time performance
• Systematic Study of Event Horizons and Pathologies of Parametrically Deformed Kerr Spacetimes
• RAPTOR I: Time-dependent radiative transfer in arbitrary spacetimes
• Explicit methods in extended phase space for inseparable Hamiltonian problems
• General relativistic polarized radiative transfer: building a dynamics-observations interface
• Generalized, energy-conserving numerical simulations of particles in general relativity. I. Time-like and null geodesics
• GRay2: A General Purpose Geodesic Integrator for Kerr Spacetimes
• Algorithms And Programs For Strong Gravitational Lensing In Kerr Space-time Including Polarization
• ynogkm: A New Public Code For Calculating time-like Geodesics In The Kerr-Newmann Spacetime

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