A differentiable N-body code for transit timing and dynamical modelling - I. Algorithm and derivatives - INSU - Institut national des sciences de l'Univers Accéder directement au contenu
Article Dans Une Revue Monthly Notices of the Royal Astronomical Society Année : 2021

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

David M. Hernandez
  • Fonction : Auteur
Zachary Langford
  • Fonction : Auteur

Résumé

When fitting N-body models to astronomical data - such as 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 favourably in speed to other integrators that are written in C.
Fichier principal
Vignette du fichier
stab2044.pdf (2.13 Mo) Télécharger le fichier
Origine : Fichiers éditeurs autorisés sur une archive ouverte

Dates et versions

insu-03748190 , version 1 (09-08-2022)

Identifiants

Citer

Eric Agol, David M. Hernandez, Zachary Langford. A differentiable N-body code for transit timing and dynamical modelling - I. Algorithm and derivatives. Monthly Notices of the Royal Astronomical Society, 2021, 507, pp.1582-1605. ⟨10.1093/mnras/stab2044⟩. ⟨insu-03748190⟩
11 Consultations
16 Téléchargements

Altmetric

Partager

Gmail Facebook X LinkedIn More