arxiv: 2605.13847 · v1 · submitted 2026-05-13 · 🌀 gr-qc · hep-th

Recognition: 2 theorem links

· Lean Theorem

Analytical Fluxes from Generic Schwarzschild Geodesics

Chris Kavanagh, Majed Khalaf, Ofri Telem

Authors on Pith no claims yet

Pith reviewed 2026-05-14 17:35 UTC · model grok-4.3

classification 🌀 gr-qc hep-th

keywords gravitational wave fluxesSchwarzschild geodesicseccentric orbitsChebyshev expansionpost-Newtonianextreme mass ratio inspiralsFourier coefficients

0 comments

The pith

A Chebyshev expansion of radiation coefficients yields analytic gravitational-wave fluxes for eccentric Schwarzschild geodesics.

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

The paper develops an analytic method to compute gravitational-wave energy and angular momentum fluxes from bound orbits of arbitrary eccentricity around a Schwarzschild black hole. It expands the Fourier coefficients of the emitted radiation in a Chebyshev basis, which reduces them to sums of simpler Keplerian-like coefficients previously obtained via the Quantum Spectral Method. This construction avoids any small-eccentricity approximation and therefore applies across a wide range of bound eccentric geodesics. When the input is taken from a 15 post-Newtonian expansion, the resulting fluxes match numerical benchmarks to relative accuracy 10^{-5} for (p,e)=(12.5,0.5) and show mode-by-mode errors below 10^{-6} for dominant modes at (p,e)=(10,0.8). The approach supplies a frequency-domain route to flux calculations needed for extreme-mass-ratio inspiral modeling.

Core claim

The authors show that expanding the Fourier coefficients of the radiation in a Chebyshev basis reduces them to sums of Keplerian-like Fourier coefficients from the Quantum Spectral Method, thereby furnishing an analytic expression for the fluxes emitted by generic bound Schwarzschild geodesics without invoking a small-eccentricity expansion.

What carries the argument

Chebyshev expansion of the radiation Fourier coefficients, which reduces them to sums of Keplerian-like coefficients.

If this is right

Fluxes can be obtained in the frequency domain for any bound eccentricity without small-e expansions.
Total flux is reproduced to relative accuracy 10^{-5} for the orbit (p,e)=(12.5,0.5) with 15PN input.
Dominant modes for the stronger-field orbit (p,e)=(10,0.8) agree to weighted errors below 10^{-6}.
The construction supplies analytic expressions suitable for EMRI waveform modeling.

Where Pith is reading between the lines

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

The same reduction could be tested on Kerr geodesics once analogous high-order inputs become available.
The method offers a fast benchmark against which fully numerical flux integrators can be validated.
Combining the frequency-domain fluxes with existing time-domain waveforms may improve computational efficiency for detector data analysis.

Load-bearing premise

The Chebyshev series for the Fourier coefficients converges rapidly when the radiation input is supplied by a 15PN expansion, even at eccentricities up to 0.8.

What would settle it

A direct high-precision numerical integration of the flux for an orbit with eccentricity significantly above 0.8 that deviates from the Chebyshev-reduced result by more than the reported tolerance would falsify the convergence claim.

Figures

Figures reproduced from arXiv: 2605.13847 by Chris Kavanagh, Majed Khalaf, Ofri Telem.

read the original abstract

We present an analytic method for computing gravitational-wave fluxes from bound Schwarzschild geodesics with arbitrary eccentricity. Our approach systematically expands the Fourier coefficients of the emitted radiation in a Chebyshev basis, allowing them to be reduced to sums of Keplerian-like Fourier coefficients previously derived in the Quantum Spectral Method. Because the construction does not rely on a small-eccentricity expansion, it applies to a broad range of bound eccentric orbits. As an illustration, we implement the method using a $15$PN-expanded input and find that it reproduces the total flux for the case $(p,e)=(12.5,0.5)$ to relative accuracy $10^{-5}$, while for the stronger-field case $(p,e)=(10,0.8)$ it yields weighted mode-by-mode errors below $10^{-6}$ for the selected dominant modes analyzed. These results provide an analytic route to frequency-domain flux calculations relevant to EMRI modeling.

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.

Desk Editor's Note private letter to a colleague

The Chebyshev reduction to Keplerian sums is a genuine technical step for arbitrary-eccentricity fluxes, but the checks stay narrow and the convergence is assumed rather than shown.

read the letter

The paper's main advance is a systematic way to expand the Fourier coefficients of radiation from bound Schwarzschild geodesics in a Chebyshev basis, then reduce them to sums of the Keplerian coefficients already derived in the Quantum Spectral Method. This construction does not invoke a small-eccentricity limit, so it targets the full range of bound eccentric orbits that matter for EMRI modeling. That reduction itself is new and cleanly executed on the terms given in the abstract. The numerical illustration with a 15PN input is a fair first test: it recovers the total flux to roughly 10^{-5} relative accuracy at (p,e)=(12.5,0.5) and keeps weighted mode-by-mode errors below 10^{-6} for the dominant modes at (10,0.8). Those results show the method can be practical in the regimes they checked. The soft spots are exactly where the stress-test note flags them. Validation is confined to two orbits, no analytic bound on Chebyshev truncation error is supplied, and there are no convergence tests for e>0.8 or p<10 where the PN input itself degrades and the geodesic Fourier content becomes richer. The rapid-convergence assumption for generic cases is therefore taken rather than demonstrated. Because the overall accuracy is also limited by the 15PN series, the claimed analytic route is still provisional outside the tested window. This work is aimed at people building frequency-domain flux tools for LISA EMRIs or looking for alternatives to purely numerical or small-eccentricity approximations. A reader already using the Quantum Spectral Method coefficients would see immediate value in the reduction step. I would send it to peer review. The technical construction is solid enough to merit referee time, provided the authors add broader convergence checks and error budgets in revision.

Referee Report

1 major / 0 minor

Summary. The paper presents an analytic method for computing gravitational-wave fluxes from bound Schwarzschild geodesics with arbitrary eccentricity. It expands the Fourier coefficients of the emitted radiation in a Chebyshev basis to reduce them to sums of Keplerian-like Fourier coefficients previously derived via the Quantum Spectral Method. The construction is illustrated by inserting a 15PN-expanded input, which reproduces the total flux for the orbit (p,e)=(12.5,0.5) to relative accuracy 10^{-5} and yields weighted mode-by-mode errors below 10^{-6} for selected dominant modes at (p,e)=(10,0.8).

Significance. If the Chebyshev reduction holds with controlled error for generic bound geodesics, the method would supply an analytic route to frequency-domain fluxes that avoids small-eccentricity expansions and directly re-uses existing Keplerian coefficients; this could be useful for EMRI waveform modeling.

major comments (1)

[Abstract] Abstract: the assertion that the method 'applies to a broad range of bound eccentric orbits' and provides 'an analytic route' for generic geodesics is load-bearing for the central claim, yet the only supporting evidence is numerical agreement at two isolated points ((12.5,0.5) and (10,0.8)) with a 15PN input; no analytic bound on Chebyshev truncation error, no convergence tests for e>0.8 or p<10, and no error budget across the parameter space are supplied.

Simulated Author's Rebuttal

1 responses · 1 unresolved

We thank the referee for their careful reading and constructive comments on the manuscript. We address the major comment point by point below and will revise the manuscript to strengthen the presentation of the method's scope and supporting evidence.

read point-by-point responses

Referee: [Abstract] Abstract: the assertion that the method 'applies to a broad range of bound eccentric orbits' and provides 'an analytic route' for generic geodesics is load-bearing for the central claim, yet the only supporting evidence is numerical agreement at two isolated points ((12.5,0.5) and (10,0.8)) with a 15PN input; no analytic bound on Chebyshev truncation error, no convergence tests for e>0.8 or p<10, and no error budget across the parameter space are supplied.

Authors: We agree that the abstract's phrasing overstates the demonstrated generality and that the numerical support is limited to the two illustrated cases. The Chebyshev expansion provides an exact reduction to sums of Keplerian coefficients in the limit of infinite order for any bound geodesic (assuming sufficient smoothness of the relevant functions), but practical truncation requires error control that is not analytically bounded in the current work. In revision we will: (i) tone down the abstract to state that the method supplies an analytic framework illustrated on bound eccentric orbits, (ii) add convergence tests and Chebyshev-coefficient decay plots for additional orbits with e=0.9 and p=8, and (iii) include an empirical error budget for the tested points. An analytic truncation-error bound for generic geodesics is not supplied and would require further analysis of geodesic smoothness properties. revision: partial

standing simulated objections not resolved

Deriving a general analytic bound on Chebyshev truncation error for arbitrary bound Schwarzschild geodesics

Circularity Check

0 steps flagged

No significant circularity: Chebyshev reduction to independent QSM coefficients plus external 15PN validation

full rationale

The derivation expands radiation Fourier coefficients in a Chebyshev basis and reduces them to sums of Keplerian-like coefficients previously obtained in the Quantum Spectral Method. This step is a change of basis that does not redefine or fit the target fluxes from themselves. The reported numerical accuracies (10^{-5} total flux at (12.5,0.5) and <10^{-6} mode errors at (10,0.8)) are obtained by feeding an external 15PN series into the method and comparing against known fluxes; they function as validation tests rather than predictions forced by construction. No load-bearing self-citation, fitted parameter renamed as prediction, or ansatz smuggled via prior work by the same authors appears in the chain. The method remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

The method relies on the existence and prior computation of Keplerian Fourier coefficients via the Quantum Spectral Method and on the validity of a 15PN expansion as input; no new free parameters, ad-hoc axioms, or invented entities are introduced in the abstract.

pith-pipeline@v0.9.0 · 5455 in / 1184 out tokens · 38113 ms · 2026-05-14T17:35:30.102768+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear
Our approach systematically expands the Fourier coefficients of the emitted radiation in a Chebyshev basis, allowing them to be reduced to sums of Keplerian-like Fourier coefficients previously derived in the Quantum Spectral Method.
IndisputableMonolith/Foundation/AlexanderDuality.lean alexander_duality_circle_linking unclear
the construction does not rely on a small-eccentricity expansion, it applies to a broad range of bound eccentric orbits

Reference graph

Works this paper leans on

38 extracted references · 32 canonical work pages · 13 internal anchors

[1]

The latter are crucial to ensure convergence over the entire intervalr p ≤r≤ ra

for a recent application of Gegenbauer polynomials in numerical self-force computations). The latter are crucial to ensure convergence over the entire intervalr p ≤r≤ ra. In this section, we assume thatf(r) is already in the form of a Laurent-like series, namely is a sum of mono- mialsr j withj∈Z. As we will see below, this requires some pre-processing in...
[2]

Light-Dark,

we showed how relativistic Keplerian Fourier ele- ments, computed using (13), can be analytically contin- ued to the unbound regime and benchmarked against post-Minkowski/post-Newtonian computations. Using the algorithm outlined in this manuscript, this continua- tion could be extended to the Schwarzschild case as well. Finally, the analytical insight gai...
[3]

B. P. Abbottet al.(LIGO Scientific, Virgo), Observation of Gravitational Waves from a Binary Black Hole Merger, Phys. Rev. Lett.116, 061102 (2016), arXiv:1602.03837 [gr-qc]

work page internal anchor Pith review Pith/arXiv arXiv 2016
[4]

B. P. Abbottet al.(LIGO Scientific, Virgo), GW170817: Observation of Gravitational Waves from a Binary Neu- tron Star Inspiral, Phys. Rev. Lett.119, 161101 (2017), arXiv:1710.05832 [gr-qc]

work page internal anchor Pith review Pith/arXiv arXiv 2017
[5]

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]

work page internal anchor Pith review Pith/arXiv arXiv 2021
[6]

Laser Interferometer Space Antenna

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

work page internal anchor Pith review Pith/arXiv arXiv 2017
[7]

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

work page arXiv 2023
[8]

LISA Definition Study Report

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

work page internal anchor Pith review Pith/arXiv arXiv 2024
[9]

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

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

work page internal anchor Pith review Pith/arXiv arXiv 2017
[10]

Self-force and radiation reaction in general relativity

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

work page internal anchor Pith review Pith/arXiv arXiv 2019
[11]

Pound and B

A. Pound and B. Wardell, Black hole perturbation the- ory and gravitational self-force, arXiv e-prints (2021), arXiv:2101.04592 [gr-qc]

work page arXiv 2021
[12]

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

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

work page internal anchor Pith review Pith/arXiv arXiv 2008
[13]

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

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

work page internal anchor Pith review Pith/arXiv arXiv 2026
[14]

Honet, A

L. Honet, A. Pound, and G. Comp` ere, Hybrid wave- form model for asymmetric spinning binaries: Self-force meets post-Newtonian theory, Phys. Rev. D113, 064035 (2026), arXiv:2510.16114 [gr-qc]

work page arXiv 2026
[15]

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 postadi- abatic terms for small-mass-ratio binaries, Phys. Rev. D 109, 124048 (2024), arXiv:2310.08927 [gr-qc]

work page arXiv 2024
[16]

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]

work page arXiv 2020
[17]

Munna, High-order post-Newtonian expansion of the generalized redshift invariant for eccentric-orbit, equato- rial extreme-mass-ratio inspirals with a spinning primary, Phys

C. Munna, High-order post-Newtonian expansion of the generalized redshift invariant for eccentric-orbit, equato- rial extreme-mass-ratio inspirals with a spinning primary, Phys. Rev. D108, 084012 (2023), arXiv:2307.11158 [gr- qc]

work page arXiv 2023
[18]

Calculation of radiation reaction effect on orbital parameters in Kerr spacetime

N. Sago and R. Fujita, Calculation of radiation reaction effect on orbital parameters in Kerr spacetime, PTEP 2015, 073E03 (2015), arXiv:1505.01600 [gr-qc]

work page internal anchor Pith review Pith/arXiv arXiv 2015
[19]

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

work page arXiv 2026
[20]

Munna and C

C. Munna and C. R. Evans, Eccentric-orbit extreme- mass-ratio-inspiral radiation II: 1PN correction to leading-logarithm and subleading-logarithm flux se- quences and the entire perturbative 4PN flux, Phys. Rev. 6 D102, 104006 (2020), arXiv:2009.01254 [gr-qc]

work page arXiv 2020
[21]

Gravitational wave snapshots of generic extreme mass ratio inspirals

S. Drasco and S. A. Hughes, Gravitational wave snap- shots 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

work page internal anchor Pith review Pith/arXiv arXiv 2006
[22]

Skoupy, G

V. Skoupy, G. Lukes-Gerakopoulos, L. V. Drummond, and S. A. Hughes, Asymptotic gravitational-wave fluxes from a spinning test body on generic orbits around a Kerr black hole, Phys. Rev. D108, 044041 (2023), arXiv:2303.16798 [gr-qc]

work page arXiv 2023
[23]

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

work page arXiv 2021
[24]

M. L. Katz, A. J. K. Chua, L. Speri, N. Warbur- ton, 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]

work page arXiv 2021
[25]

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

work page arXiv 2021
[26]

A new gravitational wave generation algorithm for particle perturbations of the Kerr spacetime

E. Harms, S. Bernuzzi, A. Nagar, and A. Zenginoglu, A new gravitational wave generation algorithm for particle perturbations of the Kerr spacetime, Class. Quant. Grav. 31, 245004 (2014), arXiv:1406.5983 [gr-qc]

work page internal anchor Pith review Pith/arXiv arXiv 2014
[27]

P. A. Sundararajan, G. Khanna, and S. A. Hughes, To- wards adiabatic waveforms for inspiral into kerr black holes. i. a new model of the source for the time do- main perturbation equation, Phys. Rev. D76(2007), arXiv:0703028 [gr-qc]

2007
[28]

Khalaf and O

M. Khalaf and O. Telem, The quantum spectral method: from atomic orbitals to classical self-force, JHEP09, 053, arXiv:2310.03798 [gr-qc]

work page arXiv
[29]

Khalaf, C.-H

M. Khalaf, C.-H. Shen, and O. Telem, Bound-Unbound Universality and the All-Order Semi-Classical Wave Function in Schwarzschild, JHEP10, arXiv:2503.23317 [gr-qc]

work page arXiv
[30]

Black Hole Perturbation Toolkit, (bhptoolkit.org)
[31]

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

work page arXiv 2025
[32]

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

1972
[33]

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)

1974
[34]

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

work page arXiv 1996
[35]

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

work page arXiv 1996
[36]

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 Cir- cular Orbits around a Schwarzschild Black Hole to the 22nd Post-Newtonian Order, Prog. Theor. Phys.128, 971 (2012), arXiv:1211.5535 [gr-qc]

work page internal anchor Pith review Pith/arXiv arXiv 2012
[37]

Kavanagh, A

C. Kavanagh, A. C. Ottewill, and B. Wardell, Analyt- ical high-order post-Newtonian expansions for spinning extreme mass ratio binaries, Phys. Rev. D93, 124038 (2016), arXiv:1601.03394 [gr-qc]

work page arXiv 2016
[38]

Whittall, L

C. Whittall, L. Barack, and O. Long, Frequency-domain self-force calculations using Gegenbauer reconstruction, Phys. Rev. D112, 124045 (2025), arXiv:2509.19439 [gr- qc]. S1 Supplemental Material for:Analytical Fluxes from Generic Schwarzschild Geodesics In this Supplemental Material, we provide further details on definitions, derivations and computations ...

work page arXiv 2025