English

Continuous-Time Modelling of Black Hole Binary Evolution with Neural ODEs

Astrophysics of Galaxies 2026-01-21 v1 Instrumentation and Methods for Astrophysics

Abstract

Pulsar timing arrays (PTAs) can detect the low-frequency stochastic gravitational-wave background (GWB) generated by an ensemble of supermassive black hole binaries (BHBs). Accurate determination of BHB merger timescales is essential for interpreting GWBs and constraining key astrophysical quantities such as black hole (BH) occupation fractions and galaxy coalescence rates. High-accuracy NN-body codes such as \texttt{Griffin} can resolve sub-pc BHB dynamics but are too costly to explore a wide range of initial conditions, motivating the need for surrogate models that emulate their long-term evolution at much lower computational cost. We investigate neural ordinary differential equations (NODEs) as surrogates for the secular orbital evolution of BHBs. Our primary contribution is a parameterised NODE (PNODE) trained on an ensemble of NN-body simulations of galaxy mergers spanning a two-dimensional parameter space defined by the initial orbital eccentricity and particle resolution (ei,N)(e_i, N), with the learned vector field explicitly conditioned on these parameters. A single PNODE thereby learns a simulation-parameter-conditioned dynamical model for the coupled evolution of the BH pair's orbital state across the ensemble, yielding smooth trajectories from which stable hardening and eccentricity growth rates can be extracted. The PNODE accurately reproduces the secular evolution of the specific orbital energy and angular momentum, and the corresponding Keplerian orbital elements, for held-out trajectories, with modest generalisation to a partially unseen high-resolution case. Combining PNODE predictions with semi-analytical prescriptions for stellar hardening and gravitational-wave emission yields BHB merger timescales consistent with those obtained from direct NN-body inputs within current theoretical uncertainties.

Keywords

Cite

@article{arxiv.2601.13019,
  title  = {Continuous-Time Modelling of Black Hole Binary Evolution with Neural ODEs},
  author = {Julian Chan and Alessia Gualandris and Payel Das},
  journal= {arXiv preprint arXiv:2601.13019},
  year   = {2026}
}

Comments

accepted for publication in Monthly Notices of the Royal Astronomical Society

R2 v1 2026-07-01T09:10:32.941Z