pith. sign in

arxiv: 2605.13847 · v2 · pith:CJUJY5UUnew · submitted 2026-05-13 · 🌀 gr-qc · astro-ph.HE· hep-ph· hep-th

Analytical Fluxes from Generic Schwarzschild Geodesics

Majed Khalaf , Chris Kavanagh , Ofri Telem This is my paper

Pith reviewed 2026-05-22 09:34 UTC · model grok-4.3

classification 🌀 gr-qc astro-ph.HEhep-phhep-th
keywords gravitational wave fluxesSchwarzschild geodesicseccentric orbitsanalytic methodsChebyshev expansionextreme mass ratio inspiralspost-Newtonian expansion
0
0 comments X

The pith

A Chebyshev expansion reduces Fourier coefficients of Schwarzschild radiation to sums of Keplerian terms for arbitrary eccentricity.

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 technique for calculating gravitational-wave energy and angular-momentum fluxes from bound geodesics around a Schwarzschild black hole at any eccentricity. The central step expands the Fourier coefficients of the emitted radiation in a Chebyshev basis, converting each coefficient into a sum of simpler Keplerian-like coefficients already known from the Quantum Spectral Method. No small-eccentricity approximation is required, so the same construction works for a wide range of bound eccentric orbits once a post-Newtonian input is supplied. A sympathetic reader cares because the resulting frequency-domain fluxes can be fed directly into waveform models for extreme-mass-ratio inspirals without time-domain integration.

Core claim

The authors demonstrate that the Fourier coefficients of the gravitational radiation from generic bound Schwarzschild geodesics can be systematically expanded in a Chebyshev basis and thereby reduced to finite sums of Keplerian-like Fourier coefficients previously derived in the Quantum Spectral Method, furnishing analytic expressions for the fluxes that remain valid for arbitrary eccentricity when a 15PN-expanded input is used.

What carries the argument

Chebyshev basis expansion of the Fourier coefficients of the emitted radiation, which reduces them to sums of Keplerian-like coefficients from the Quantum Spectral Method.

If this is right

  • The method applies to a broad range of bound eccentric orbits without any small-eccentricity restriction.
  • With 15PN input it reproduces the total flux for (p,e)=(12.5,0.5) to relative accuracy 10^{-5}.
  • For the stronger-field case (p,e)=(10,0.8) it yields weighted mode-by-mode errors below 10^{-6} for selected dominant modes.
  • It supplies an analytic route to frequency-domain flux calculations needed for extreme-mass-ratio inspirals.

Where Pith is reading between the lines

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

  • The same Chebyshev reduction could be applied to Kerr geodesics once an analogous set of base coefficients becomes available.
  • Accuracy for stronger-field orbits might increase by substituting exact numerical fluxes for the 15PN input inside the expansion.
  • These closed-form fluxes could lower the computational cost of generating long, accurate waveforms for space-based detectors.

Load-bearing premise

The Fourier coefficients of the radiation admit a convergent Chebyshev expansion when the input is supplied as a 15PN series for the chosen orbital parameters.

What would settle it

A high-precision numerical computation of the total flux for the orbit (p,e)=(12.5,0.5) that differs from the analytic result by more than 10^{-5} relative error would falsify the claimed accuracy.

Figures

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

Figure 1
Figure 1. Figure 1: FIG. 1. Radiation Flux per ( [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
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 extreme-mass-ratio inspirals.

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 manuscript presents an analytic method for computing gravitational-wave fluxes from bound Schwarzschild geodesics with arbitrary eccentricity. The approach expands the Fourier coefficients of the emitted radiation in a Chebyshev basis, reducing them to sums of Keplerian-like Fourier coefficients previously derived in the Quantum Spectral Method. The construction avoids small-eccentricity expansions and is illustrated by implementing the method with a 15PN-expanded input, which reproduces the total flux to 10^{-5} relative accuracy for (p,e)=(12.5,0.5) and yields weighted mode-by-mode errors below 10^{-6} for dominant modes at (10,0.8).

Significance. If the Chebyshev reduction is shown to be robust and the input quantities are demonstrated to be sufficiently accurate, the method would provide a valuable analytic framework for frequency-domain flux calculations relevant to extreme-mass-ratio inspirals. The avoidance of small-eccentricity approximations and the systematic reduction to prior Keplerian coefficients represent clear strengths that could aid waveform modeling.

major comments (2)
  1. [Illustration section (abstract)] Illustration section (as described in the abstract): The reported accuracies rely on a 15PN-expanded input. For the case (p,e)=(10,0.8), r_min ≈ 5.55M implies a post-Newtonian expansion parameter of order 0.18. At this value, the 15PN series may exhibit poor convergence, so the observed 10^{-5}–10^{-6} agreement may primarily confirm numerical stability of the Chebyshev summation on an approximate input rather than validating the analytic reduction procedure for exact Schwarzschild geodesics. A PN-order convergence test or direct comparison against exact geodesic quantities is needed to support the central claim.
  2. [Method description] Method description: The reduction of the Chebyshev-expanded Fourier coefficients to sums of Quantum Spectral Method coefficients is load-bearing for the analytic claim, yet the manuscript supplies no explicit construction details, basis truncation criteria, or error bounds on the mapping. Without these, it is difficult to assess whether the procedure remains exact (or controlled) for arbitrary eccentricity.
minor comments (2)
  1. The abstract refers to 'weighted mode errors' without defining the weighting; this definition should be supplied explicitly in the main text near the reported results.
  2. Consider including a brief table or figure showing the Chebyshev truncation order used and its effect on the flux accuracy for the two illustrated orbits.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful and constructive review. The comments highlight important points regarding the validation of our illustration and the need for greater detail in the method. We address each major comment below and will incorporate revisions to strengthen the manuscript.

read point-by-point responses

  1. Referee: [Illustration section (abstract)] Illustration section (as described in the abstract): The reported accuracies rely on a 15PN-expanded input. For the case (p,e)=(10,0.8), r_min ≈ 5.55M implies a post-Newtonian expansion parameter of order 0.18. At this value, the 15PN series may exhibit poor convergence, so the observed 10^{-5}–10^{-6} agreement may primarily confirm numerical stability of the Chebyshev summation on an approximate input rather than validating the analytic reduction procedure for exact Schwarzschild geodesics. A PN-order convergence test or direct comparison against exact geodesic quantities is needed to support the central claim.

    Authors: We agree that the 15PN-expanded input may have limited accuracy for the (p,e)=(10,0.8) case, where the post-Newtonian parameter is not small, and that the reported errors could in part reflect the numerical stability of the Chebyshev procedure applied to an approximate input. To address this, we will add a dedicated subsection to the revised manuscript that includes a PN-order convergence study for the input quantities at the reported orbital parameters. We will also compare the resulting fluxes against available numerical geodesic data for selected cases to better isolate the performance of the analytic reduction itself. revision: yes

  2. Referee: [Method description] Method description: The reduction of the Chebyshev-expanded Fourier coefficients to sums of Quantum Spectral Method coefficients is load-bearing for the analytic claim, yet the manuscript supplies no explicit construction details, basis truncation criteria, or error bounds on the mapping. Without these, it is difficult to assess whether the procedure remains exact (or controlled) for arbitrary eccentricity.

    Authors: We acknowledge that additional explicit details on the reduction procedure would improve the manuscript. In the revised version we will expand the method section to provide the explicit algebraic construction of the Chebyshev-to-QSM mapping, specify the truncation criteria (including the number of Chebyshev modes retained and the rationale for that choice), and derive analytic error bounds based on the convergence properties of Chebyshev expansions for the relevant periodic functions. These additions will clarify that the reduction is exact in the infinite-basis limit and that the truncation error is controlled independently of eccentricity. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation introduces independent Chebyshev reduction step

full rationale

The paper's core step is a new Chebyshev-basis expansion of radiation Fourier coefficients that mathematically permits reduction to sums of Keplerian-like coefficients previously obtained in the Quantum Spectral Method. This is presented as a consequence of the basis choice rather than a redefinition or fit of the target fluxes. The 15PN input is explicitly an implementation choice for numerical illustration, not a load-bearing part of the analytic derivation. No equations are shown to be equivalent by construction to the inputs, no fitted parameters are relabeled as predictions, and the self-citation to the Quantum Spectral Method is not required to be load-bearing for the new reduction technique itself. The derivation chain therefore remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The central claim rests on the ability to expand radiation Fourier coefficients in a Chebyshev basis and on the accuracy of the 15PN input used for validation; no explicit free parameters beyond the PN order are named.

free parameters (1)
  • PN expansion order
    15PN-expanded input chosen for the numerical illustration; affects the accuracy benchmark reported.
axioms (1)
  • domain assumption Fourier coefficients of emitted radiation admit a convergent Chebyshev expansion for bound eccentric Schwarzschild geodesics
    Invoked as the core construction that allows reduction to prior Keplerian coefficients.

pith-pipeline@v0.9.0 · 5688 in / 1387 out tokens · 41990 ms · 2026-05-22T09:34:54.415687+00:00 · methodology

Review history (2 revisions) →

discussion (0)

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

Lean theorems connected to this paper

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

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

38 extracted references · 38 canonical work pages · 14 internal anchors

  1. [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...

    work page

  2. [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...

    work page

  3. [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. [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. [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. [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. [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 internal anchor Pith review arXiv 2023

  8. [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. [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. [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. [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. [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. [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. [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. [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. [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. [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. [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. [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. [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. [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. [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. [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. [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. [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. [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. [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]

    work page 2007

  28. [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. [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. [30]

    Black Hole Perturbation Toolkit, (bhptoolkit.org)

    work page

  31. [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. [32]

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

    work page 1972

  33. [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)

    work page 1974

  34. [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. [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. [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. [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. [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