Self-force-consistent transition from inspiral to plunge for generic (inclined, eccentric) orbits around a Kerr black hole

For a compact body of mass $\mu$ inspiraling into a massive Kerr black hole of mass $M \gg \mu$, the final passage from slow adiabatic inspiral to plunge passes through a short 'transition' regime near the boundary separating bound from plunging orbits. Ori and Thorne (2000) solved this transition analytically only for circular, equatorial orbits, finding a universal Painlevé-I-type equation and predicting a duration scaling as $\eta^{-1/5}$ ($\eta = \mu/M$), and they posed the generalization to astrophysically realistic orbits — eccentric, inclined, and fully generic — as an open problem, blocked at the time by the unknown radiation-reaction evolution of the Carter constant and by neglected conservative self-force corrections. Later work has resolved the eccentric-equatorial and inclined-circular special cases, computed the Carter-constant evolution from first principles at both adiabatic (flux-balance) and first-order self-force level, and built self-force-corrected transition models for quasi-circular equatorial and near-extremal binaries. The surviving open problem is the transition from inspiral to plunge for fully generic inclined-and-eccentric orbits with first-principles evolution of all three orbit integrals (energy, axial angular momentum, Carter constant) and conservative self-force corrections, together with the associated waveforms and an updated assessment of the observability of the transition regime for realistic generic LISA extreme-mass-ratio inspirals.

Background

An extreme-mass-ratio inspiral (EMRI) is the months-to-years-long gravitational-wave-driven capture of a stellar-mass compact object (mass $\mu$) into a massive black hole (mass $M$, with mass ratio $\eta = \mu/M \sim 10^{-4}$–$10^{-7}$); these events are a key source class for the space-based Laser Interferometer Space Antenna (LISA). In the leading 'adiabatic' approximation the small body moves through a sequence of geodesics of the Kerr spacetime (the unique stationary, rotating vacuum black-hole solution, characterized by mass $M$ and spin parameter $a$). Bound Kerr geodesics are characterized by three conserved quantities: the orbital energy $E$, the axial angular momentum $L_z$, and the Carter constant $Q$ (a conserved quantity arising from a hidden symmetry of the Kerr metric, roughly the square of the orbital angular momentum component perpendicular to the spin, which measures the orbit's polar extent and hence its inclination). Radiation reaction — the backreaction of the emitted gravitational waves — is computed in black-hole perturbation theory by solving the Teukolsky equation (the master frequency-domain equation for perturbations of Kerr) for the emitted fluxes, and the orbit-averaged rates of change of $E$, $L_z$, and $Q$ follow from flux-balance laws. Beyond this averaged, dissipative description, the gravitational self-force (the local first-order-in-$\eta$ correction to the geodesic equation, with both dissipative and conservative pieces) governs post-adiabatic corrections. Bound orbits end at the separatrix, the boundary in the space of orbit parameters (semi-latus rectum $p$, eccentricity $e$, inclination) beyond which no stable bound geodesic exists; for circular equatorial orbits this boundary is the innermost stable circular orbit (ISCO).

Ori and Thorne (2000) analyzed the 'transition regime' in which, as the ISCO is approached, the restoring force of the radial effective potential weakens faster than radiation reaction, so the body can no longer track the sequence of shrinking geodesics. Expanding the effective potential to cubic order, $V(R,\xi) = (2\alpha/3) R^3 - 2\beta R\xi$, with $R = r - r_{\rm isco}$ (radii in units of $M$) and $\xi$ a rescaled measure of the deviation of angular momentum from its ISCO value, and taking $d\xi/d\tau = -\kappa\eta$, they reduced the motion to a single universal dimensionless equation of Painlevé-I type — Painlevé I being the first of the six classical nonlinear Painlevé equations, $d^2 y/dt^2 = 6 y^2 + t$, whose required solution is selected by asymptotic matching to the inspiral at early times — valid for all spins and small mass ratios, with radial offset scaling as $\eta^{2/5}$, duration as $\eta^{-1/5}$ (in units of $M$), and cycle count also proportional to $\eta^{-1/5}$. Their conclusions section posed the open issues that define this problem: the analysis was confined to circular equatorial orbits although realistic inspirals are expected to be strongly noncircular and nonequatorial; for generic orbits the radiation-reaction-induced rate of change of the Carter constant was not then known from first principles; the analysis neglected, as explicitly acknowledged in the paper's main text, the $O(\eta)$ conservative self-force corrections to the effective potential and to the ISCO values of $r$, $E$, and $L$; and their LISA detectability estimate (signal-to-noise ratio of order 1 at 1 Gpc for a $10\,M_\odot$ object into a $10^6\,M_\odot$ hole) used the straw-man LISA noise curve of the time and circular-equatorial geometry, a limitation whose removal their conclusions cite as motivation for generalizing the analysis.

Subsequent developments define the current frontier. The averaged evolution of all three constants for generic orbits is now known from first principles through flux-balance formulae and Teukolsky flux computations, and the first-order self-force on generic bound Kerr geodesics (dissipative and conservative) has been computed. Self-forced inspiral models exist for eccentric equatorial and for inclined quasi-circular ('spherical') orbits, and adiabatic waveforms for fully generic orbits have been generated, but these inspiral models end at the separatrix. The transition itself has been generalized to eccentric equatorial orbits (within the modern two-timescale framework, in which the slow evolution of the orbit constants and the fast orbital phases are separated — a framework introduced by Hinderer and Flanagan and developed by Hughes and collaborators) and to inclined quasi-circular orbits (via a generalized Ori-Thorne procedure formulated in Mino time, a time parameter that separates the radial and polar motions), and the Compère–Küchler–Pound program proved that for quasi-circular equatorial orbits the transition remains governed by Painlevé I even when all self-force effects are included, yielding self-force-based transition-to-plunge waveform models; near-extremal spins have also been treated. What has not been built is the corresponding transition-to-plunge model for the fully generic inclined-and-eccentric geometry with first-principles Carter-constant evolution and conservative self-force corrections, together with its waveforms and an updated observability assessment. The modern program has been extended to late-time waveforms, including the excitation of the remnant's ringdown modes and its dependence on the orbital anomaly angles (Lim et al. 2019; Becker and Hughes 2025).

Problem Statement

Construct the transition from adiabatic inspiral to plunge for a compact body of mass $\mu$ orbiting a Kerr black hole of mass $M \gg \mu$ ($\eta = \mu/M \ll 1$) on a fully generic — inclined and eccentric — bound orbit, replacing the circular-equatorial restriction of the original Ori-Thorne analysis, and establish its observable consequences. Concretely, the still-open problem comprises: (1) the worldline through the separatrix and transition regime in which all three orbit integrals (energy $E$, axial angular momentum $L_z$, and Carter constant $Q$) evolve under first-principles radiation reaction — the orbit-averaged flux-balance input at adiabatic order and, at post-adiabatic order, the self-force input — so that the transition dynamics, its scaling with $\eta$ and spin, and the polar-angle motion are determined without the approximated or parameterized $dQ/dt$ on which earlier generic prescriptions relied; and (2) incorporation of the $O(\eta)$ conservative self-force corrections to the effective potential and to the separatrix location and associated constants (known to shift the transition coefficients and hence the scaling-law prefactors), which have been computed only for quasi-circular equatorial transitions. Components (1) and (2) are precisely the generalization posed by Ori and Thorne (2000). (3) The gravitational waveform generated by the generic transition and subsequent plunge, including how the transition and the two orbital anomaly angles (the radial and polar phase parameters along the orbit) control the excitation of the remnant's ringdown modes (the damped quasinormal-mode oscillations of the settling post-merger black hole). (4) An updated assessment of the observability of the transition regime for realistic generic orbits, with a realistic LISA response and current noise curves, generalizing Ori and Thorne's signal-to-noise estimate, which used the straw-man LISA noise curve and circular-equatorial geometry of the time; their conclusions identify this generalization as a key motivation for the problem. These four components are jointly required for a complete answer and form one objective — a first-principles, usable description of the last orbits of astrophysically realistic EMRIs. In the appropriate limits the generic solution must be shown to reproduce the already-solved circular-equatorial (Ori-Thorne/Painlevé-I with self-force), eccentric-equatorial, and inclined-circular transition dynamics and their scaling with $\eta$ and spin, together with the resulting waveforms. For context, the solved special cases are cited in the Background section; the open target is the generic case and its self-force-consistent, observable consequences, not a re-derivation of the special cases.

References

1. Amos Ori and Kip S. Thorne, 'The Transition from Inspiral to Plunge for a Compact Body in a Circular Equatorial Orbit Around a Massive, Spinning Black Hole', Phys. Rev. D 62, 124022 (2000), DOI 10.1103/PhysRevD.62.124022, https://arxiv.org/abs/gr-qc/0003032
2. R. O'Shaughnessy, 'Transition from inspiral to plunge for eccentric equatorial Kerr orbits', Phys. Rev. D 67, 044004 (2003), DOI 10.1103/PhysRevD.67.044004, https://arxiv.org/abs/gr-qc/0211023
3. Éanna É. Flanagan and Tanja Hinderer, 'Evolution of the Carter constant for inspirals into a black hole: effect of the black hole quadrupole', Phys. Rev. D 75, 124007 (2007), DOI 10.1103/PhysRevD.75.124007, https://arxiv.org/abs/0704.0389
4. Katsuhiko Ganz, Wataru Hikida, Hiroyuki Nakano, Norichika Sago and Takahiro Tanaka, 'Adiabatic Evolution of three Constants of Motion for Greatly Inclined Orbits in Kerr spacetime', Prog. Theor. Phys. 117, 1041 (2007), DOI 10.1143/PTP.117.1041, https://arxiv.org/abs/gr-qc/0702054
5. Tanja Hinderer and Éanna É. Flanagan, 'Two-timescale analysis of extreme mass ratio inspirals in Kerr spacetime: Orbital motion', Phys. Rev. D 78, 064028 (2008), DOI 10.1103/PhysRevD.78.064028, https://arxiv.org/abs/0805.3337
6. Pranesh A. Sundararajan, 'The transition from adiabatic inspiral to geodesic plunge for a compact object around a massive Kerr black hole: Generic orbits', Phys. Rev. D 77, 124050 (2008), DOI 10.1103/PhysRevD.77.124050, https://arxiv.org/abs/0803.4482
7. Michael Kesden, 'Transition from adiabatic inspiral to plunge into a spinning black hole', Phys. Rev. D 83, 104011 (2011), DOI 10.1103/PhysRevD.83.104011, https://arxiv.org/abs/1101.3749
8. Maarten van de Meent, 'Gravitational self-force on generic bound geodesics in Kerr spacetime', Phys. Rev. D 97, 104033 (2018), DOI 10.1103/PhysRevD.97.104033, https://arxiv.org/abs/1711.09607
9. Anuj Apte and Scott A. Hughes, 'Exciting black hole modes via misaligned coalescences: I. Inspiral, transition, and plunge trajectories using a generalized Ori-Thorne procedure', Phys. Rev. D 100, 084031 (2019), DOI 10.1103/PhysRevD.100.084031, https://arxiv.org/abs/1901.05901
10. Halston Lim, Gaurav Khanna, Anuj Apte and Scott A. Hughes, 'Exciting black hole modes via misaligned coalescences: II. The mode content of late-time coalescence waveforms', Phys. Rev. D 100, 084032 (2019), DOI 10.1103/PhysRevD.100.084032, https://arxiv.org/abs/1901.05902
11. Soichiro Isoyama, Ryuichi Fujita, Hiroyuki Nakano, Norichika Sago and Takahiro Tanaka, 'Flux-balance formulae for extreme mass-ratio inspirals', Prog. Theor. Exp. Phys. 2019, 013E01 (2019), DOI 10.1093/ptep/pty136, https://arxiv.org/abs/1809.11118
12. Geoffrey Compère, Kwinten Fransen and Caroline Jonas, 'Transition from inspiral to plunge into a highly spinning black hole', Class. Quantum Grav. 37, 095013 (2020), DOI 10.1088/1361-6382/ab79d3, https://arxiv.org/abs/1909.12848
13. Ollie Burke, Jonathan R. Gair and Joan Simón, 'Transition from Inspiral to Plunge: A Complete Near-Extremal Trajectory and Associated Waveform', Phys. Rev. D 101, 064026 (2020), DOI 10.1103/PhysRevD.101.064026, https://arxiv.org/abs/1909.12846
14. Geoffrey Compère and Lorenzo Küchler, 'Self-consistent adiabatic inspiral and transition motion', Phys. Rev. Lett. 126, 241106 (2021), DOI 10.1103/PhysRevLett.126.241106, https://arxiv.org/abs/2102.12747
15. Scott A. Hughes, Niels Warburton, Gaurav Khanna, Alvin J. K. Chua and Michael L. Katz, 'Adiabatic waveforms for extreme mass-ratio inspirals via multivoice decomposition in time and frequency', Phys. Rev. D 103, 104014 (2021), DOI 10.1103/PhysRevD.103.104014, https://arxiv.org/abs/2102.02713
16. Geoffrey Compère and Lorenzo Küchler, 'Asymptotically matched quasi-circular inspiral and transition-to-plunge in the small mass ratio expansion', SciPost Phys. 13, 043 (2022), DOI 10.21468/SciPostPhys.13.2.043, https://arxiv.org/abs/2112.02114
17. Soichiro Isoyama, Ryuichi Fujita, Alvin J. K. Chua, Hiroyuki Nakano, Adam Pound and Norichika Sago, 'Adiabatic waveforms from extreme-mass-ratio inspirals: an analytical approach', Phys. Rev. Lett. 128, 231101 (2022), DOI 10.1103/PhysRevLett.128.231101, https://arxiv.org/abs/2111.05288
18. Philip Lynch, Maarten van de Meent and Niels Warburton, 'Eccentric self-forced inspirals into a rotating black hole', Class. Quantum Grav. 39, 145004 (2022), DOI 10.1088/1361-6382/ac7507, https://arxiv.org/abs/2112.05651
19. Philip Lynch, Maarten van de Meent and Niels Warburton, 'Self-forced inspirals with spin-orbit precession', Phys. Rev. D 109, 084072 (2024), DOI 10.1103/PhysRevD.109.084072, https://arxiv.org/abs/2305.10533
20. Lorenzo Küchler, Geoffrey Compère, Leanne Durkan and Adam Pound, 'Self-force framework for transition-to-plunge waveforms', SciPost Phys. 17, 056 (2024), DOI 10.21468/SciPostPhys.17.2.056, https://arxiv.org/abs/2405.00170
21. Devin R. Becker and Scott A. Hughes, 'Transition from adiabatic inspiral to plunge for eccentric binaries', Phys. Rev. D 111, 064003 (2025), DOI 10.1103/PhysRevD.111.064003, https://arxiv.org/abs/2410.09160
22. Guillaume Lhost and Geoffrey Compère, 'Approach to the separatrix with eccentric orbits', SciPost Phys. Core 8, 059 (2025), DOI 10.21468/SciPostPhysCore.8.3.059, https://arxiv.org/abs/2412.04249
23. Josh Mathews and Adam Pound, 'Post-adiabatic waveform-generation framework for asymmetric precessing binaries', Phys. Rev. D 112, 104078 (2025), DOI 10.1103/rbkb-qnxv, https://arxiv.org/abs/2501.01413
24. Devin R. Becker, Scott A. Hughes and Gaurav Khanna, 'Gravitational waves from the late inspiral, transition, and plunge of small-mass-ratio eccentric binaries', Phys. Rev. D 113, 084046 (2026), DOI 10.1103/jyx7-jtt4, https://arxiv.org/abs/2511.21897
25. Loïc Honet, Lorenzo Küchler, Adam Pound and Geoffrey Compère, 'Transition-to-plunge self-force waveforms with a spinning primary', Phys. Rev. D 113, 044051 (2026), DOI 10.1103/sq6y-qv8h, https://arxiv.org/abs/2510.13958
26. Ariadna Ribes Metidieri, Béatrice Bonga, Badri Krishnan and José Luis Jaramillo, 'Universality in the Transition from Inspiral to Plunge: High-Accuracy Analytic Solutions and Catastrophe Theory', arXiv:2606.13786, https://arxiv.org/abs/2606.13786
27. Jonathan R. Gair and Kostas Glampedakis, 'Improved approximate inspirals of test bodies into Kerr black holes', Phys. Rev. D 73, 064037 (2006), DOI 10.1103/PhysRevD.73.064037, https://arxiv.org/abs/gr-qc/0510129