Numerically stable equations for the orbital evolution of compact object binaries
Abstract
The orbital and eccentricity evolution for compact object binaries through gravitational wave emission first derived by Peters and Mathews are used extensively throughout the gravitational wave community for calculating the orbital evolution and merger time of compact binaries. While improved calculations of the binary merger time have been the focus of several investigations since, the orbital evolution has not received the same attention. As the equations lack a closed form solution, a numerical integrator is required, but standard methods typically break when the point of merger is overstepped. We present a rewrite of Peters' equations in -space, which allows common numerical solvers to converge. This leads to a more numerically robust and computationally efficient method for evolving compact binaries due to gravitational wave emission, reducing the number of function evaluations by 60\% to 70\% in our tests.
Keywords
Cite
@article{arxiv.2603.20124,
title = {Numerically stable equations for the orbital evolution of compact object binaries},
author = {Max M. Briel and Jeff J. Andrews},
journal= {arXiv preprint arXiv:2603.20124},
year = {2026}
}