Computing gravitational radiation from highly eccentric bound binaries in the small-mass-ratio limit
Abstract
Precise strong-gravity models describing the evolution of bound compact binary systems provide important information about the gravitational radiation such binaries generate and how they evolve due to the radiation's backreaction. The limit of very high eccentricity, e → 1, is challenging to model: frequency-domain techniques require many harmonics to accurately describe radiation in this limit, and time-domain techniques can be slow to converge due to the systems' long timescales. Many recent astrophysical analyses have nonetheless emphasized the importance of understanding gravitational-wave emission in this regime, highlighting the fact that many systems may evolve through e ≃ 1 en route to becoming detectable sources, and may enter the sensitive band of detectors with substantial eccentricity. In this paper, we show how to adapt frequency-domain methods for computing gravitational radiation in the small-mass-ratio limit to effectively handle high eccentricity. The key piece of this analysis is an integral over the source of the frequency-domain Teukolsky equation. We examine features of this integral that make it difficult to evaluate as e → 1, and study methods to circumvent these difficulties. We find a simple fix which allows us to compute gravitational-wave amplitudes and the associated fluxes of energy and angular momentum with ease up to eccentricities e = 0.99; with patience (and stamina), we have pushed to e = 0.999, though numerical accuracy appears to degrade for such extreme cases. Interestingly, we find that the well-known and often used classic formulas of Peters and Mathews are subject to a systematic error which persists at quite wide separation. This has implications for models of compact binary formation which evolve through e ≃ 1.
Related
- Pulsed Accretion onto Eccentric Binaries in Highly Misaligned Circumbinary Disks
- Pulsed Accretion onto Eccentric Binaries in Highly Misaligned Circumbinary Disks
- Parameter Estimation for Eccentric Supermassive Black Hole Binaries with Pulsar Timing Arrays
- A parametrized test of general relativity for inspiralling eccentric binaries in LISA
- Are there really some eccentric eclipsing binaries having orbital periods below 1 d?
BibTeX
@misc{orona2026computing,
title = {{Computing gravitational radiation from highly eccentric bound binaries in the small-mass-ratio limit}},
author = {Diego Andrés Rivera Orona and Scott A. Hughes},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.36327},
eprint = {2609.36327},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.36327},
}
RIS
TY - UNPB
TI - Computing gravitational radiation from highly eccentric bound binaries in the small-mass-ratio limit
AU - Diego Andrés Rivera Orona
AU - Scott A. Hughes
PY - 2026
DA - 2026/09/28
PB - arXiv
DO - 10.48550/arxiv.2609.36327
UR - https://arxiv.org/abs/2609.36327
AB - Precise strong-gravity models describing the evolution of bound compact binary systems provide important information about the gravitational radiation such binaries generate and how they evolve due to the radiation's backreaction. The limit of very high eccentricity, e → 1, is challenging to model: frequency-domain techniques require many harmonics to accurately describe radiation in this limit, and time-domain techniques can be slow to converge due to the systems' long timescales. Many recent astrophysical analyses have nonetheless emphasized the importance of understanding gravitational-wave emission in this regime, highlighting the fact that many systems may evolve through e ≃ 1 en route to becoming detectable sources, and may enter the sensitive band of detectors with substantial eccentricity. In this paper, we show how to adapt frequency-domain methods for computing gravitational radiation in the small-mass-ratio limit to effectively handle high eccentricity. The key piece of this analysis is an integral over the source of the frequency-domain Teukolsky equation. We examine features of this integral that make it difficult to evaluate as e → 1, and study methods to circumvent these difficulties. We find a simple fix which allows us to compute gravitational-wave amplitudes and the associated fluxes of energy and angular momentum with ease up to eccentricities e = 0.99; with patience (and stamina), we have pushed to e = 0.999, though numerical accuracy appears to degrade for such extreme cases. Interestingly, we find that the well-known and often used classic formulas of Peters and Mathews are subject to a systematic error which persists at quite wide separation. This has implications for models of compact binary formation which evolve through e ≃ 1.
ER -