pith. machine review for the scientific record. sign in

arxiv: 2605.09250 · v2 · submitted 2026-05-10 · 🌀 gr-qc

Recognition: no theorem link

Efficient and Stable Computation of Gravitational-Wave Fluxes from Generic Kerr Orbits via a Unified HeunC Framework

Authors on Pith no claims yet

Pith reviewed 2026-05-13 07:43 UTC · model grok-4.3

classification 🌀 gr-qc
keywords gravitational wave fluxesKerr black holesconfluent Heun functionsTeukolsky equationsextreme mass ratio inspiralsself-forcenumerical integration
0
0 comments X

The pith

Reformulating the Teukolsky equations with confluent Heun functions yields stable, high-precision gravitational-wave fluxes from generic Kerr orbits.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a method to compute the energy and angular momentum carried away by gravitational waves from particles orbiting Kerr black holes in arbitrary orbits. It rewrites the governing angular and radial equations using confluent Heun functions instead of the usual approaches. A hybrid analytic continuation technique then finds the required connection coefficients without searching for auxiliary parameters, producing solutions that converge everywhere. An adaptive quadrature rule handles the rapid oscillations in the source integrals for eccentric or inclined orbits. The resulting code reaches relative errors near 10 to the minus 11 for fluxes summed over 168 modes while running several times faster than current packages, especially on high-order modes.

Core claim

The authors reformulate both the angular and radial Teukolsky equations in terms of confluent Heun functions. A hybrid analytic continuation algorithm computes the connection coefficients, eliminating auxiliary parameter dependence and directly yielding globally convergent solutions and scattering amplitudes. For generic orbits, an adaptive bi-power mapping quadrature resolves highly oscillatory source integrands. Benchmarks show that for the total radiative flux summed over 168 low-order modes, relative errors reach order 10^{-11}, with costs reduced by factors of 3-13 compared to existing packages, and up to 60 times speedup for oscillatory high-order modes.

What carries the argument

Confluent Heun functions for the angular and radial Teukolsky equations, with hybrid analytic continuation to obtain connection coefficients and adaptive bi-power mapping quadrature to integrate oscillatory sources.

If this is right

  • Total radiative fluxes summed over 168 low-order modes achieve relative errors of order 10^{-11}.
  • Computational costs drop by factors of 3 to 13 relative to GeneralizedSasakiNakamura.jl and pybhpt packages.
  • Highly oscillatory high-order modes run up to 60 times faster than specialized oscillatory integrators.
  • The method supplies the numerical foundation for high-order self-force calculations and rapid high-precision waveform generation.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The reduced cost per orbit could support denser sampling of parameter space when building template banks for space-based detectors.
  • The same Heun reformulation may extend to computing other quantities such as waveforms or tidal responses in Kerr perturbation theory.
  • Global convergence without auxiliary tuning could simplify inclusion of higher-order self-force corrections in future models.

Load-bearing premise

The hybrid analytic continuation algorithm for connection coefficients in the confluent Heun framework eliminates auxiliary parameter dependence and yields globally convergent solutions without introducing instabilities or inaccuracies in the strong-field or high-frequency regimes for generic orbits.

What would settle it

An independent high-resolution numerical integration of the Teukolsky equations for a highly eccentric, inclined orbit that produces a total flux differing from the Heun-based result by more than one part in 10^{11}.

Figures

Figures reproduced from arXiv: 2605.09250 by Changkai Chen, Jiliang Jing, Zhoujian Cao.

Figure 1
Figure 1. Figure 1: FIG. 1. Schematic of the bidirectional analytic continuation scheme. [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. SWSH functions and eigenvalues with the mechanical precision( [PITH_FULL_IMAGE:figures/full_fig_p009_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Relative errors of SWSH functions computed with machine precision ( [PITH_FULL_IMAGE:figures/full_fig_p009_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Integrand [PITH_FULL_IMAGE:figures/full_fig_p012_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Integrand [PITH_FULL_IMAGE:figures/full_fig_p013_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Logarithmic error [PITH_FULL_IMAGE:figures/full_fig_p014_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. QNM eigenvalues and errors with [PITH_FULL_IMAGE:figures/full_fig_p016_7.png] view at source ↗
read the original abstract

Modeling extreme-mass-ratio inspirals hinges on the accurate and efficient computation of gravitational-wave fluxes from generic Kerr orbits. Conventional frequency-domain techniques are often limited by costly auxiliary parameter searches and numerical instabilities in the strong-field or high-frequency regimes. We address these challenges by reformulating both the angular and radial Teukolsky equations in terms of confluent Heun functions. Employing a hybrid analytic continuation algorithm to compute the connection coefficients eliminates the dependence on auxiliary parameters, directly yielding globally convergent solutions and scattering amplitudes. To resolve the highly oscillatory source integrands for generic orbits, we implement an adaptive bi-power mapping quadrature. Comprehensive benchmarks under standard double-precision arithmetic demonstrate that, for the total radiative flux summed over 168 low-order modes, our method achieves relative errors of order $10^{-11}$, with computational costs typically reduced by factors of 3--13 compared to the state-of-the-art GeneralizedSasakiNakamura. jl and pybhpt packages. Notably, for highly oscillatory high-order modes, our framework achieves a speedup of up to 60 times compared to specialized oscillatory integrators like GeneralizedSasakiNakamura. jl. These demonstrated gains in precision and efficiency establish the framework as a robust tool for strong-field perturbation theory, providing the numerical foundation for high-order self-force calculations and rapid, high-precision waveform generation.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper claims to provide an efficient and stable method for computing gravitational-wave fluxes from generic Kerr orbits by reformulating both the angular and radial Teukolsky equations in terms of confluent Heun functions. It introduces a hybrid analytic continuation algorithm for the connection coefficients that eliminates auxiliary parameter dependence, combined with an adaptive bi-power mapping quadrature to handle highly oscillatory source integrands for generic orbits. Comprehensive benchmarks under double precision report relative errors of order 10^{-11} for the total radiative flux summed over 168 low-order modes, with typical speedups of 3-13x (up to 60x for high-order oscillatory modes) relative to GeneralizedSasakiNakamura.jl and pybhpt.

Significance. If the stability and convergence claims hold, the work would advance strong-field perturbation theory by supplying a robust, parameter-free numerical foundation for high-order self-force calculations and rapid high-precision waveform generation for extreme-mass-ratio inspirals. The concrete error levels, direct use of established HeunC properties, and reported speedups over state-of-the-art packages constitute clear strengths in efficiency and precision.

major comments (2)
  1. [Methods section on connection coefficients] Methods section on connection coefficients: the hybrid analytic continuation algorithm is presented as eliminating auxiliary parameter dependence while yielding globally convergent solutions without instabilities. The manuscript supplies only empirical benchmarks (a<0.99, moderate frequencies) rather than an analytic bound on truncation or branch-cut errors for high-eccentricity, high-frequency generic orbits where the radial Teukolsky source is highly oscillatory; this is load-bearing for the claimed stability and 3-60x speedups.
  2. [Results/benchmarks] Results/benchmarks: the reported relative errors of order 10^{-11} for summed fluxes over 168 modes and the 60x speedup for high-order modes are quantified, but the manuscript must specify the exact orbit parameters (eccentricity, inclination, spin a, frequencies) used in those tests to substantiate coverage of the strong-field regime where conventional methods encounter instabilities.
minor comments (2)
  1. [Abstract] The abstract and methods should explicitly name the specialized oscillatory integrator within GeneralizedSasakiNakamura.jl that is used for the 60x comparison to improve clarity.
  2. Notation for the bi-power mapping quadrature could be defined more explicitly (e.g., the precise form of the mapping function and adaptive criterion) to aid reproducibility.

Simulated Author's Rebuttal

2 responses · 1 unresolved

We thank the referee for the careful and constructive report. We address each major comment below and have revised the manuscript to improve clarity, reproducibility, and coverage of the strong-field regime.

read point-by-point responses
  1. Referee: [Methods section on connection coefficients] Methods section on connection coefficients: the hybrid analytic continuation algorithm is presented as eliminating auxiliary parameter dependence while yielding globally convergent solutions without instabilities. The manuscript supplies only empirical benchmarks (a<0.99, moderate frequencies) rather than an analytic bound on truncation or branch-cut errors for high-eccentricity, high-frequency generic orbits where the radial Teukolsky source is highly oscillatory; this is load-bearing for the claimed stability and 3-60x speedups.

    Authors: We agree that the stability and speedup claims rest on numerical evidence rather than a closed-form analytic error bound. The hybrid continuation is constructed directly from the known global connection formulas of the confluent Heun equation, which are free of auxiliary parameters by construction; the adaptive bi-power quadrature is designed to control the oscillatory integrand error independently of the connection step. While a general analytic bound for arbitrary eccentricity and frequency would be desirable, its derivation is a substantial theoretical undertaking beyond the scope of the present work. In the revision we have added a new subsection (III.C) that (i) tabulates observed truncation and quadrature errors over an extended grid reaching e=0.9 and |ω|≈1.2, (ii) reports the maximum relative deviation from reference solutions, and (iii) explicitly delineates the parameter region where the method has been validated. These additions make the empirical foundation of the stability claim more transparent. revision: partial

  2. Referee: [Results/benchmarks] Results/benchmarks: the reported relative errors of order 10^{-11} for summed fluxes over 168 modes and the 60x speedup for high-order modes are quantified, but the manuscript must specify the exact orbit parameters (eccentricity, inclination, spin a, frequencies) used in those tests to substantiate coverage of the strong-field regime where conventional methods encounter instabilities.

    Authors: We fully concur that the benchmark parameters must be stated explicitly for reproducibility and to demonstrate coverage of the strong-field domain. The revised manuscript now includes Table II, which lists the complete set of orbital parameters (a, e, ι, p, Ω_r, Ω_θ, Ω_φ) for every flux computation reported in Section IV, including the high-order mode runs that produced the 60× speedup. These orbits include strong-field cases with a=0.99, e up to 0.85, and frequencies corresponding to near-horizon motion, precisely the regime where conventional methods encounter instabilities. revision: yes

standing simulated objections not resolved
  • A rigorous analytic bound on truncation and branch-cut errors of the hybrid analytic continuation algorithm for arbitrary high-eccentricity and high-frequency Kerr orbits

Circularity Check

0 steps flagged

No circularity: reformulation via established Heun functions and quadrature

full rationale

The paper reformulates the angular and radial Teukolsky equations in terms of confluent Heun functions, applies a hybrid analytic continuation algorithm for connection coefficients, and uses adaptive bi-power mapping quadrature for oscillatory integrands. These steps rely on known properties of Heun functions and standard numerical techniques rather than any self-referential definitions, fitted parameters renamed as predictions, or load-bearing self-citations that reduce the central claims to the inputs by construction. Benchmarks compare against external packages (GeneralizedSasakiNakamura.jl, pybhpt) with empirical error measures, providing independent validation. No step in the derivation chain exhibits the required reduction to prior outputs or self-citations.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on the mathematical properties of confluent Heun functions for the Teukolsky equations and the correctness of the hybrid continuation algorithm; these are treated as standard tools in the field.

axioms (2)
  • standard math Confluent Heun functions provide globally valid solutions to the angular and radial Teukolsky equations on Kerr backgrounds
    Invoked in the reformulation step of the abstract.
  • domain assumption The hybrid analytic continuation algorithm computes connection coefficients without auxiliary parameter searches
    Central to eliminating instabilities mentioned in the abstract.

pith-pipeline@v0.9.0 · 5545 in / 1398 out tokens · 43168 ms · 2026-05-13T07:43:43.205699+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

121 extracted references · 121 canonical work pages · 8 internal anchors

  1. [1]

    HeunC′(α+1, β+1, γ+1, δ+1, η+1; 1 2) + HeunC(α+1, β+1, γ+1, δ+1, η+1; 1

  2. [2]

    HeunC′(α−1, β−1, γ−1, δ−1, η−1; 1 2).(29) whereHeunC ′(· · ·;x)denotes thex-derivative. The condi- tionW θ = 0yields a transcendental equation that uniquely determines the eigenvalue ˆλand the angular separation con- stants sAℓm(aω); further details on the extraction procedure are provided in our previous work [101]. Although this fixes the eigenvalue, th...

  3. [3]

    This configuration serves as a crit- 11 TABLE II

    Circular orbits for schwarzschild geometry We first validate our framework against the simplest case: a test particle in an equatorial circular orbit around a Schwarzschild black hole. This configuration serves as a crit- 11 TABLE II. Comparison of three methods for GW fluxes of Schwarzschild BHs in circular orbits.Ndenotes the effective working precision...

  4. [4]

    the Fundamen- tal Research Funds for the Central Universities

    General orbits for Kerr geometry The most stringent test for any flux calculation method lies in generic (inclined and eccentric) orbits around a Kerr black hole. In this regime, the integrandG ∞,H ℓmω exhibits complex oscillatory behavior that poses significant challenges for nu- merical quadrature. To illustrate these features, Figures. 4 and 5 display ...

  5. [5]

    The symmetry of the parameterαin CHE (12) can be given as follows, HeunC(α, β, γ, δ, η;x) = e −αxHeunC(−α, β, γ, δ,η;x), β /∈Z.(A1) whereZis the set of integers

  6. [6]

    The symmetry of the parameterγin CHE (12) can be given as follows, HeunC(α, β, γ, δ, η;x) = (1−x) −γHeunC(α, β,−γ, δ, η;x), β /∈Z.(A2)

  7. [7]

    It can be seen that the param- eters(α, β, γ)in the general solution (22) of the ATE (3) and the general solution (36) of the RTE (5) are in- terchangeable

    The symmetry of the parameterβis only displayed in the Teukolsky equations. It can be seen that the param- eters(α, β, γ)in the general solution (22) of the ATE (3) and the general solution (36) of the RTE (5) are in- terchangeable. The symmetry properties of the general solutionS ℓm andR ℓm are as follows, Sℓm(α, β, γ;x) =S ℓm(−α, β, γ;x),(A3a) Sℓm(α, β,...

  8. [8]

    Substituting Eq

    Power-series summation The regular local solutionH C1(x)admits a power-series expansion [97]: HC1(x) = ∞X n=0 bnxn,(B2) where the coefficientsb n satisfy the three-term recurrence re- lation associated with the first parametrization: Pnbn =Q nbn−1 +R nbn−2,(B3) with Pn =n(n+β),(B4a) Qn = (n−1)(n+β+γ−α) +η+ 1 2 (β+ 1)(γ+ 1−α)−1 ,(B4b) Rn =α n−1 + β+γ 2 +δ....

  9. [9]

    To circumvent this, we adopt an alternative method based on the Sturm–Liouville structure of the CHE

    Alternative form via Wronskian derivative For higher-order modes, the power-series summation often requires precision significantly exceeding standard double- precision arithmetic, leading to substantial computational overhead. To circumvent this, we adopt an alternative method based on the Sturm–Liouville structure of the CHE. Berens et al.demonstrated t...

  10. [10]

    T. D. Abbottet al.(LIGO Scientific, Virgo), Improved analysis of GW150914 using a fully spin-precessing waveform Model, Phys. Rev. X6, 041014 (2016), arXiv:1606.01210 [gr-qc]. 20

  11. [11]

    B. P. Abbottet al.(LIGO Scientific, Virgo), Binary black hole mergers in the first advanced LIGO observing run, Phys. Rev. X6, 041015 (2016), [Erratum: Phys.Rev.X 8, 039903 (2018)], arXiv:1606.04856 [gr-qc]

  12. [12]

    B. P. Abbottet al.(LIGO Scientific, Virgo), GWTC-1: A gravitational-wave transient catalog of compact binary mergers observed by LIGO and Virgo during the first and second observing runs, Phys. Rev. X9, 031040 (2019), arXiv:1811.12907 [astro-ph.HE]

  13. [13]

    GWTC-2: Compact Binary Coalescences Observed by LIGO and Virgo During the First Half of the Third Observing Run

    R. Abbottet al.(LIGO Scientific, Virgo), GWTC-2: Compact binary coalescences observed by LIGO and Virgo during the first half of the third observing run, Phys. Rev. X11, 021053 (2021), arXiv:2010.14527 [gr-qc]

  14. [14]

    Abbott, T

    R. Abbottet al.(KAGRA, VIRGO, LIGO Scientific), Popula- tion of merging compact binaries inferred using gravitational waves through GWTC-3, Phys. Rev. X13, 011048 (2023), arXiv:2111.03634 [astro-ph.HE]

  15. [15]

    Laser Interferometer Space Antenna

    P. Amaro-Seoaneet al.(LISA), Laser Interferometer Space Antenna, (2017), arXiv:1702.00786 [astro-ph.IM]

  16. [16]

    LISA Definition Study Report

    M. Colpiet al.(LISA), LISA Definition Study Report, (2024), arXiv:2402.07571 [astro-ph.CO]

  17. [17]

    TianQin: a space-borne gravitational wave detector

    J. Luoet al.(TianQin), TianQin: a space-borne gravitational wave detector, Classical Quantum Gravity33, 035010 (2016), arXiv:1512.02076 [astro-ph.IM]

  18. [18]

    The TianQin project: current progress on science and technology,

    J. Meiet al.(TianQin), The TianQin project: current progress on science and technology, Prog. Theor. Exp. Phys.2021, 05A107 (2021), arXiv:2008.10332 [gr-qc]

  19. [19]

    Li et al., Gravitational wave astronomy with TianQin, Rept

    E.-K. Liet al., Gravitational wave astronomy with Tian- Qin, Rept. Prog. Phys.88, 056901 (2025), arXiv:2409.19665 [astro-ph.GA]

  20. [20]

    Ruan, Z.-K

    W.-H. Ruan, Z.-K. Guo, R.-G. Cai, and Y .-Z. Zhang, Taiji pro- gram: Gravitational-wave sources, Int. J. Mod. Phys. A35, 2050075 (2020), arXiv:1807.09495 [gr-qc]

  21. [21]

    Y . Gong, J. Luo, and B. Wang, Concepts and status of Chinese space gravitational wave detection projects, Nature Astron.5, 881 (2021), arXiv:2109.07442 [astro-ph.IM]

  22. [22]

    Amaro-Seoane, Relativistic dynamics and extreme mass ra- tio inspirals, Living Rev

    P. Amaro-Seoane, Relativistic dynamics and extreme mass ra- tio inspirals, Living Rev. Rel.21, 4 (2018), arXiv:1205.5240 [astro-ph.CO]

  23. [23]

    P. A. Seoaneet al.(LISA), Astrophysics with the Laser In- terferometer Space Antenna, Living Rev. Rel.26, 2 (2023), arXiv:2203.06016 [gr-qc]

  24. [24]

    Kawamura et al., PTEP2021, 05A105 (2021), arXiv:2006.13545 [gr-qc]

    S. Kawamuraet al., Current status of space gravitational wave antenna DECIGO and B-DECIGO, Prog. Theor. Exp. Phys. 2021, 05A105 (2021), arXiv:2006.13545 [gr-qc]

  25. [25]

    The Science of the Einstein Telescope

    A. Abacet al.(ET), The Science of the Einstein Telescope, (2025), arXiv:2503.12263 [gr-qc]

  26. [26]

    Cosmic Explorer: The U.S. Contribution to Gravitational-Wave Astronomy beyond LIGO

    D. Reitzeet al., Cosmic Explorer: The U.S. contribution to gravitational-wave astronomy beyond LIGO, Bull. Am. As- tron. Soc.51, 035 (2019), arXiv:1907.04833 [astro-ph.IM]

  27. [27]

    Punturoet al., The Einstein Telescope: A third-generation gravitational wave observatory, Classical Quantum Gravity 27, 194002 (2010)

    M. Punturoet al., The Einstein Telescope: A third-generation gravitational wave observatory, Classical Quantum Gravity 27, 194002 (2010)

  28. [28]

    Science with the space-based interferometer LISA. V: Extreme mass-ratio inspirals

    S. Babak, J. Gair, A. Sesana, E. Barausse, C. F. Sopuerta, C. P. L. Berry, E. Berti, P. Amaro-Seoane, A. Petiteau, and A. Klein, Science with the space-based interferometer LISA. V: Extreme mass-ratio inspirals, Phys. Rev. D95, 103012 (2017), arXiv:1703.09722 [gr-qc]

  29. [29]

    C. P. L. Berry, S. A. Hughes, C. F. Sopuerta, A. J. K. Chua, A. Heffernan, K. Holley-Bockelmann, D. P. Mihaylov, M. C. Miller, and A. Sesana, The unique potential of extreme mass- ratio inspirals for gravitational-wave astronomy, Bull. Am. Astron. Soc.51, 42 (2019), arXiv:1903.03686 [astro-ph.HE]

  30. [30]

    Barack and C

    L. Barack and C. Cutler, LISA capture sources: Approxi- mate waveforms, signal-to-noise ratios, and parameter esti- mation accuracy, Phys. Rev. D69, 082005 (2004), arXiv:gr- qc/0310125

  31. [31]

    J. R. Gair, M. Vallisneri, S. L. Larson, and J. G. Baker, Testing general relativity with low-frequency, space-based gravitational-wave detectors, Living Rev. Rel.16, 7 (2013), arXiv:1212.5575 [gr-qc]

  32. [32]

    Maselli, N

    A. Maselli, N. Franchini, L. Gualtieri, T. P. Sotiriou, S. Barsanti, and P. Pani, Detecting fundamental fields with LISA observations of gravitational waves from ex- treme mass-ratio inspirals, Nature Astron.6, 464 (2022), arXiv:2106.11325 [gr-qc]

  33. [33]

    Speri, S

    L. Speri, S. Barsanti, A. Maselli, T. P. Sotiriou, N. Warbur- ton, M. van de Meent, A. J. K. Chua, O. Burke, and J. Gair, Probing fundamental physics with extreme mass ratio inspi- rals: Full Bayesian inference for scalar charge, Phys. Rev. D 113, 023036 (2026), arXiv:2406.07607 [gr-qc]

  34. [34]

    Bonga, H

    B. Bonga, H. Yang, and S. A. Hughes, Tidal resonance in extreme mass-ratio inspirals, Phys. Rev. Lett.123, 101103 (2019), arXiv:1905.00030 [gr-qc]

  35. [35]

    Speri, A

    L. Speri, A. Antonelli, L. Sberna, S. Babak, E. Barausse, J. R. Gair, and M. L. Katz, Probing accretion physics with gravitational waves, Phys. Rev. X13, 021035 (2023), arXiv:2207.10086 [gr-qc]

  36. [36]

    Khalvati, A

    H. Khalvati, A. Santini, F. Duque, L. Speri, J. Gair, H. Yang, and R. Brito, Impact of relativistic waveforms in LISA’s sci- ence objectives with extreme-mass-ratio inspirals, Phys. Rev. D111, 082010 (2025), arXiv:2410.17310 [gr-qc]

  37. [37]

    J. R. Gair, C. Tang, and M. V olonteri, LISA extreme-mass- ratio inspiral events as probes of the black hole mass func- tion, Phys. Rev. D81, 104014 (2010), arXiv:1004.1921 [astro- ph.GA]

  38. [38]

    C. E. A. Chapman-Bird, C. P. L. Berry, and G. Woan, Rapid determination of LISA sensitivity to extreme mass ratio inspi- rals with machine learning, Mon. Not. Roy. Astron. Soc.522, 6043 (2023), arXiv:2212.06166 [astro-ph.HE]

  39. [39]

    Amaro-Seoane, J

    P. Amaro-Seoane, J. R. Gair, M. Freitag, M. Coleman Miller, I. Mandel, C. J. Cutler, and S. Babak, Astrophysics, detection and science applications of intermediate- and extreme mass- ratio inspirals, Classical Quantum Gravity24, R113 (2007), arXiv:astro-ph/0703495

  40. [40]

    C. L. MacLeod and C. J. Hogan, Precision of Hubble constant derived using black hole binary absolute distances and statis- tical redshift information, Phys. Rev. D77, 043512 (2008), arXiv:0712.0618 [astro-ph]

  41. [41]

    Laghi, N

    D. Laghi, N. Tamanini, W. Del Pozzo, A. Sesana, J. Gair, S. Babak, and D. Izquierdo-Villalba, Gravitational-wave cos- mology with extreme mass-ratio inspirals, Mon. Not. Roy. As- tron. Soc.508, 4512 (2021), arXiv:2102.01708 [astro-ph.CO]

  42. [42]

    Auclair et al

    P. Auclairet al.(LISA Cosmology Working Group), Cosmol- ogy with the Laser Interferometer Space Antenna, Living Rev. Rel.26, 5 (2023), arXiv:2204.05434 [astro-ph.CO]

  43. [43]

    Burke, G

    O. Burke, G. A. Piovano, N. Warburton, P. Lynch, L. Speri, C. Kavanagh, B. Wardell, A. Pound, L. Durkan, and J. Miller, Assessing the importance of first postadiabatic terms for small-mass-ratio binaries, Phys. Rev. D109, 124048 (2024), arXiv:2310.08927 [gr-qc]

  44. [44]

    Liang and H

    B. Liang and H. Wang, Recent advances in simulation-based inference for gravitational wave data analysis, Astronomical Techniques and Instruments3, 93 (2026)

  45. [45]

    Barack, Gravitational self force in extreme mass-ratio inspirals, Classical Quantum Gravity26, 213001 (2009), arXiv:0908.1664 [gr-qc]

    L. Barack, Gravitational self force in extreme mass-ratio inspirals, Classical Quantum Gravity26, 213001 (2009), arXiv:0908.1664 [gr-qc]

  46. [46]

    Pound, B

    A. Pound, B. Wardell, N. Warburton, and J. Miller, Second- 21 order self-force calculation of gravitational binding energy in compact binaries, Phys. Rev. Lett.124, 021101 (2020), arXiv:1908.07419 [gr-qc]

  47. [47]

    Warburton, A

    N. Warburton, A. Pound, B. Wardell, J. Miller, and L. Durkan, Gravitational-wave energy flux for compact binaries through second order in the mass ratio, Phys. Rev. Lett.127, 151102 (2021), arXiv:2107.01298 [gr-qc]

  48. [48]

    Wardell, A

    B. Wardell, A. Pound, N. Warburton, J. Miller, L. Durkan, and A. Le Tiec, Gravitational waveforms for compact bina- ries from second-order self-force theory, Phys. Rev. Lett.130, 241402 (2023), arXiv:2112.12265 [gr-qc]

  49. [49]

    The motion of point particles in curved spacetime

    E. Poisson, A. Pound, and I. Vega, The Motion of point par- ticles in curved spacetime, Living Rev. Rel.14, 7 (2011), arXiv:1102.0529 [gr-qc]

  50. [50]

    Self-force and radiation reaction in general relativity

    L. Barack and A. Pound, Self-force and radiation reaction in general relativity, Rept. Prog. Phys.82, 016904 (2019), arXiv:1805.10385 [gr-qc]

  51. [51]

    Pound and B

    A. Pound and B. Wardell, Black hole perturbation theory and gravitational self-force,Handbook of Gravitational Wave As- tronomy, , 1 (2020)

  52. [52]

    Mino, Perturbative approach to an orbital evolution around a supermassive black hole, Phys

    Y . Mino, Perturbative approach to an orbital evolution around a supermassive black hole, Phys. Rev. D67, 084027 (2003), arXiv:gr-qc/0302075

  53. [53]

    Drasco, E

    S. Drasco, E. E. Flanagan, and S. A. Hughes, Computing in- spirals in Kerr in the adiabatic regime. I. The scalar case, Clas- sical Quantum Gravity22, S801 (2005), arXiv:gr-qc/0505075

  54. [54]

    N. Sago, T. Tanaka, W. Hikida, K. Ganz, and H. Nakano, The Adiabatic evolution of orbital parameters in the Kerr spacetime, Prog. Theor. Phys.115, 873 (2006), arXiv:gr- qc/0511151

  55. [55]

    Flux-balance formulae

    S. Isoyama, R. Fujita, H. Nakano, N. Sago, and T. Tanaka, “Flux-balance formulae” for extreme mass-ratio inspi- rals, Prog. Theor. Exp. Phys.2019, 013E01 (2019), arXiv:1809.11118 [gr-qc]

  56. [56]

    A. M. Grant, Flux-balance laws for spinning bodies under the gravitational self-force, Phys. Rev. D111, 084015 (2025), arXiv:2406.10343 [gr-qc]

  57. [57]

    S. A. Teukolsky, Rotating black holes - separable wave equations for gravitational and electromagnetic perturbations, Phys. Rev. Lett.29, 1114 (1972)

  58. [58]

    S. A. Teukolsky, Perturbations of a rotating black hole. 1. Fun- damental equations for gravitational electromagnetic and neu- trino field perturbations, Astrophys. J.185, 635 (1973)

  59. [59]

    W. H. Press and S. A. Teukolsky, Perturbations of a rotating black hole. II. Dynamical stability of the Kerr metric, Astro- phys. J.185, 649 (1973)

  60. [60]

    S. A. Teukolsky and W. H. Press, Perturbations of a rotating black hole. III - Interaction of the hole with gravitational and electromagnetic radiation, Astrophys. J.193, 443 (1974)

  61. [61]

    R. M. Wald, On perturbations of a Kerr black hole, J. Math. Phys.14, 1453 (1973)

  62. [62]

    Two timescale analysis of extreme mass ratio inspirals in Kerr. I. Orbital Motion

    T. Hinderer and E. E. Flanagan, Two timescale analysis of ex- treme mass ratio inspirals in Kerr. I. Orbital motion, Phys. Rev. D78, 064028 (2008), arXiv:0805.3337 [gr-qc]

  63. [63]

    Miller and A

    J. Miller and A. Pound, Two-timescale evolution of extreme- mass-ratio inspirals: Waveform generation scheme for quasi- circular orbits in Schwarzschild spacetime, Phys. Rev. D103, 064048 (2021), arXiv:2006.11263 [gr-qc]

  64. [64]

    Post-adiabatic self-force waveforms: slowly spinning primary and precessing secondary

    J. Mathews, B. Wardell, A. Pound, and N. Warburton, Posta- diabatic self-force waveforms: Slowly spinning primary and precessing secondary, Phys. Rev. D113, 064034 (2026), arXiv:2510.16113 [gr-qc]

  65. [65]

    S. A. Hughes, The Evolution of circular, nonequatorial or- bits of Kerr black holes due to gravitational wave emission, Phys. Rev. D61, 084004 (2000), [Erratum: Phys.Rev.D 63, 049902 (2001), Erratum: Phys.Rev.D 65, 069902 (2002), Er- ratum: Phys.Rev.D 67, 089901 (2003), Erratum: Phys.Rev.D 78, 109902 (2008), Erratum: Phys.Rev.D 90, 109904 (2014)], arXiv:g...

  66. [66]

    S. A. Hughes, Evolution of circular, nonequatorial orbits of Kerr black holes due to gravitational wave emission. II. In- spiral trajectories and gravitational wave forms, Phys. Rev. D 64, 064004 (2001), [Erratum: Phys.Rev.D 88, 109902 (2013)], arXiv:gr-qc/0104041

  67. [67]

    Drasco and S

    S. Drasco and S. A. Hughes, Rotating black hole orbit func- tionals in the frequency domain, Phys. Rev. D69, 044015 (2004), arXiv:astro-ph/0308479

  68. [68]

    Gravitational wave snapshots of generic extreme mass ratio inspirals

    S. Drasco and S. A. Hughes, Gravitational wave snapshots of generic extreme mass ratio inspirals, Phys. Rev. D73, 024027 (2006), [Erratum: Phys.Rev.D 88, 109905 (2013), Erratum: Phys.Rev.D 90, 109905 (2014)], arXiv:gr-qc/0509101

  69. [69]

    S. A. Hughes, N. Warburton, G. Khanna, A. J. K. Chua, and M. L. Katz, Adiabatic waveforms for extreme mass-ratio in- spirals via multivoice decomposition in time and frequency, Phys. Rev. D103, 104014 (2021), [Erratum: Phys.Rev.D 107, 089901 (2023)], arXiv:2102.02713 [gr-qc]

  70. [70]

    S. Mano, H. Suzuki, and E. Takasugi, Analytic solutions of the Regge-Wheeler equation and the postMinkowskian expan- sion, Prog. Theor. Phys.96, 549 (1996), arXiv:gr-qc/9605057

  71. [71]

    S. Mano, H. Suzuki, and E. Takasugi, Analytic solutions of the Teukolsky equation and their low frequency expansions, Prog. Theor. Phys.95, 1079 (1996), arXiv:gr-qc/9603020

  72. [72]

    Fujita, W

    R. Fujita, W. Hikida, and H. Tagoshi, An efficient numerical method for computing gravitational waves induced by a parti- cle moving on eccentric inclined orbits around a Kerr black hole, Prog. Theor. Phys.121, 843 (2009), arXiv:0904.3810 [gr-qc]

  73. [73]

    Gravitational Waves from a Particle in Circular Orbits around a Schwarzschild Black Hole to the 22nd Post-Newtonian Order

    R. Fujita, Gravitational waves from a particle in circu- lar orbits around a Schwarzschild black hole to the 22nd post-Newtonian order, Prog. Theor. Phys.128, 971 (2012), arXiv:1211.5535 [gr-qc]

  74. [74]

    Calculation of radiation reaction effect on orbital parameters in Kerr spacetime

    N. Sago and R. Fujita, Calculation of radiation reaction ef- fect on orbital parameters in Kerr spacetime, Prog. Theor. Exp. Phys.2015, 073E03 (2015), arXiv:1505.01600 [gr-qc]

  75. [75]

    Munna, Analytic post-Newtonian expansion of the energy and angular momentum radiated to infinity by eccentric-orbit nonspinning extreme-mass-ratio inspirals to the 19th order, Phys

    C. Munna, Analytic post-Newtonian expansion of the energy and angular momentum radiated to infinity by eccentric-orbit nonspinning extreme-mass-ratio inspirals to the 19th order, Phys. Rev. D102, 124001 (2020), arXiv:2008.10622 [gr-qc]

  76. [76]

    N. Sago, R. Fujita, and H. Nakano, Post-Newtonian templates for phase evolution of spherical extreme mass ratio inspirals, Phys. Rev. D111, 064043 (2025), arXiv:2411.09147 [gr-qc]

  77. [77]

    J. C. Castillo, C. R. Evans, C. Kavanagh, J. Neef, A. Ottewill, and B. Wardell, Post-Newtonian expansion of gravitational en- ergy and angular momentum fluxes: Inclined spherical orbits about a Kerr black hole, Phys. Rev. D111, 084004 (2025), arXiv:2411.09700 [gr-qc]

  78. [78]

    N. Sago, R. Fujita, S. Isoyama, and H. Nakano, Secular evo- lution of orbital parameters for general bound orbits in Kerr spacetime, (2026), arXiv:2603.27941 [gr-qc]

  79. [79]

    A. J. K. Chua, M. L. Katz, N. Warburton, and S. A. Hughes, Rapid generation of fully relativistic extreme-mass- ratio-inspiral waveform templates for LISA data analysis, Phys. Rev. Lett.126, 051102 (2021), arXiv:2008.06071 [gr- qc]

  80. [80]

    M. L. Katz, A. J. K. Chua, L. Speri, N. Warburton, and S. A. Hughes, Fast extreme-mass-ratio-inspiral waveforms: New tools for millihertz gravitational-wave data analysis, Phys. Rev. D104, 064047 (2021), arXiv:2104.04582 [gr-qc]. 22

Showing first 80 references.