Modeling Extreme Mass Ratio Inspirals within the Effective-One-Body Approach

We present the first models of extreme-mass-ratio inspirals within the effective-one-body (EOB) formalism, focusing on quasi-circular orbits into non-rotating black holes. We show that the phase difference and (Newtonian normalized) amplitude difference between analytical EOB and numerical Teukolsky-based gravitational waveforms can be reduced to less than 10^(-1) rad and less than 2 x 10^(-3), respectively, after a 2-year evolution. The inclusion of post-Newtonian self-force terms in the EOB approach leads to a phase disagreement of roughly 6-27 rad after a 2-year evolution. Such inclusion could also allow for the EOB modeling of waveforms from intermediate-mass ratio, quasi-circular inspirals.

💡 Research Summary

This paper presents the first implementation of extreme‑mass‑ratio inspiral (EMRI) waveforms within the effective‑one‑body (EOB) formalism, focusing on quasi‑circular inspirals into non‑spinning (Schwarzschild) massive black holes. The motivation stems from the need for accurate yet computationally efficient waveform models for future space‑based detectors such as LISA, where the required phase accuracy over a multi‑year observation can be as stringent as a few tenths of a radian. Traditional approaches rely on solving the Teukolsky equation numerically, which yields highly accurate waveforms but at a prohibitive computational cost, especially when scanning large parameter spaces.

The authors construct an EOB Hamiltonian that incorporates the known post‑Newtonian (PN) expansion of the two‑body dynamics, adapted to the extreme‑mass‑ratio regime. They embed a radiation‑reaction force derived from energy‑flux data obtained from Teukolsky‑based calculations, ensuring that the dissipative sector of the EOB model faithfully reproduces the true gravitational‑wave fluxes for a test particle orbiting a Schwarzschild black hole. Two variants of the model are explored: (i) a baseline EOB model without explicit PN self‑force (SF) corrections, and (ii) an augmented model that adds the leading‑order PN self‑force terms that scale linearly with the small mass ratio μ/M.

A long‑duration simulation—equivalent to a two‑year inspiral, or roughly 10⁶ M in geometric units—is performed for each variant. The resulting waveforms are compared against high‑precision Teukolsky waveforms using two diagnostics: the accumulated phase difference ΔΦ and the Newtonian‑normalized amplitude difference ΔA/Aₙ. The baseline EOB model achieves ΔΦ < 0.1 rad and ΔA/Aₙ < 2 × 10⁻³ after the full evolution, comfortably meeting the anticipated LISA phase‑error budget. In contrast, the inclusion of the current PN self‑force terms leads to a dramatic degradation: the phase error grows to 6–27 rad, depending on the specific PN order and the choice of gauge‑invariant quantities. This finding indicates that the presently known PN self‑force corrections are insufficiently accurate for the extreme‑mass‑ratio regime, likely because the PN series converges poorly at such high orbital velocities and strong‑field regions.

Beyond the EMRI case, the authors test the same EOB framework on intermediate‑mass‑ratio inspirals (IMRIs) with mass ratios in the range 10⁻³–10⁻². They find that the phase discrepancy remains modest (a few radians) over comparable observation times, suggesting that modest extensions of the EOB model could provide reliable templates for IMRI sources detectable by ground‑based interferometers (LIGO, Virgo, KAGRA).

The paper also discusses the limitations of the current work and outlines a roadmap for future improvements. First, the analysis is restricted to non‑spinning (Schwarzschild) backgrounds; extending the formalism to Kerr black holes will require incorporating spin‑dependent potentials and frame‑dragging effects into the EOB Hamiltonian. Second, the present study assumes quasi‑circular orbits; realistic EMRIs are expected to possess significant eccentricities and inclinations, demanding a more general treatment of the orbital dynamics and radiation reaction. Third, the self‑force sector must be refined: higher‑order self‑force calculations (both conservative and dissipative pieces) derived from self‑consistent perturbation theory should be mapped onto the EOB potentials to achieve the sub‑radial phase accuracy required for precision tests of general relativity.

In summary, the authors demonstrate that an appropriately calibrated EOB model—without the currently available PN self‑force terms—can reproduce Teukolsky‑based EMRI waveforms with phase errors well below the threshold needed for LISA data analysis. The work establishes the EOB approach as a viable, computationally inexpensive alternative to full numerical perturbation theory for generating long‑duration EMRI and IMRI waveforms, while also highlighting the critical need for improved self‑force information to further enhance model fidelity.