A differentiable N-body code for transit timing and dynamical modeling. I. Algorithm and derivatives

Abstract

When fitting N-body models to astronomical data - including transit times, radial velocity, and astrometric positions at observed times - the derivatives of the model outputs with respect to the initial conditions can help with model optimization and posterior sampling. Here we describe a general-purpose symplectic integrator for arbitrary orbital architectures, including those with close encounters, which we have recast to maintain numerical stability and precision for small step sizes. We compute the derivatives of the N-body coordinates and velocities as a function of time with respect to the initial conditions and masses by propagating the Jacobian along with the N-body integration. For the first time we obtain the derivatives of the transit times with respect to the initial conditions and masses using the chain rule, which is quicker and more accurate than using finite differences or automatic differentiation. We implement this algorithm in an open source package, NbodyGradient.jl, written in the Julia language, which has been used in the optimization and error analysis of transit-timing variations in the TRAPPIST-1 system. We present tests of the accuracy and precision of the code, and show that it compares favorably in speed to other integrators which are written in C.

arXiv comment24 pages; 7 figures. Revised for submission to MNRAS. Links in equations are provided to code at the repository https://github.com/ericagol/NbodyGradient.jl Some documentation and example scripts for running the code are provided in the repository

References in corpus (37)

• DHP Framework: Digital Health Passports Using Blockchain -- Use case on international tourism during the COVID-19 pandemic
• On detecting terrestrial planets with timing of giant planet transits
• Kepler-16: A Transiting Circumbinary Planet
• WHFast: A fast and unbiased implementation of a symplectic Wisdom-Holman integrator for long term gravitational simulations
• The nature of the TRAPPIST-1 exoplanets
• Extracting Planet Mass and Eccentricity From TTV data
• 55 Cancri (Copernicus): A Multi-Planet System with a Hot Super-Earth and a Jupiter Analogue
• KOI-126: A Triply-Eclipsing Hierarchical Triple with Two Low-Mass Stars
• REBOUNDx: A Library for Adding Conservative and Dissipative Forces to Otherwise Symplectic N-body Integrations
• Kepler Planet Masses and Eccentricities from TTV Analysis
• TTVFast: An efficient and accurate code for transit timing inversion problems
• Secure TTV Mass Measurements: Ten Kepler Exoplanets between 3 and 8 Earth Masses with Diverse Densities and Incident Fluxes
• Measurement of planet masses with transit timing variations due to synodic "chopping" effects
• The GENGA Code: Gravitational Encounters in N-body simulations with GPU Acceleration
• Secular dynamics of hierarchical multiple systems composed of nested binaries, with an arbitrary number of bodies and arbitrary hierarchical structure. First applications to multiplanet and multistar systems
• Transit timing to first order in eccentricity
• The Compositional Diversity of Low-Mass Exoplanets
• TRADES: A new software to derive orbital parameters from observed transit times and radial velocities. Revisiting Kepler-11 and Kepler-9
• Numerical and Analytical Modelling of Transit Time Variations
• The Effect of Conjunctions on the Transit Timing Variations of Exoplanets
• HARPS-N radial velocities confirm the low densities of the Kepler-9 planets
• Dynamics and Transit Variations of Resonant Exoplanets
• Systemic: A Testbed for Characterizing the Detection of Extrasolar Planets. II. Numerical approaches to the Transit Timing Inverse Problem
• Second-order variational equations for N-body simulations
• Fast Inversion Method for Determination of Planetary Parameters from Transit Timing Variations
• A Fast and Accurate Universal Kepler Solver without Stumpff Series
• Transit timing variations for planets near eccentricity-type mean motion resonances
• A study of symplectic integrators for planetary system problems: error analysis and comparisons
• Kepler Object of Interest Network II. Photodynamical modelling of Kepler-9 over 8 years of transit observations
• Symplectic fourth-order maps for the collisional N-body problem
• Symplectic integration for the collisional gravitational $N$-body problem
• Fast and reliable symplectic integration for planetary system $N$-body problems
• Time-symmetric integration in astrophysics
• A Simplified Photodynamical Model for Planetary Mass Determination in Low-Eccentricity Multi-Transiting Systems
• Analysis of radial velocity variations in multiple planetary systems
• Modal Decomposition of TTV - Inferring Planet Masses and Eccentricities
• A Keplerian-based Hamiltonian Splitting for Gravitational $N$-body Simulations

Cited by in corpus (9)

• TOI-1136 is a Young, Coplanar, Aligned Planetary System in a Pristine Resonant Chain
• celmech: A Python package for celestial mechanics
• GRIT: a package for structure-preserving simulations of gravitationally interacting rigid-bodies
• A young progenitor for the most common planetary systems in the Galaxy
• A Radio Technosignature Search of TRAPPIST-1 with the Allen Telescope Array
• Switching integrators reversibly in the astrophysical $N$-body problem
• Reliable event detection for Taylor methods in astrodynamics
• Analytic auto-differentiable $\Lambda$CDM cosmography
• A differentiable N-body code for transit timing and dynamical modelling -- II. Photodynamics