Binary black hole scattering in the extreme-mass-ratio limit: time-domain waveform

Andrea Geralico

REVIEW

Binary black hole scattering in the extreme-mass-ratio limit: time-domain waveform

Andrea Geralico

General Relativity and Quantum Cosmology · Submitted 24 Sep 2026

Pith's verdict: Genuinely new 5PM/7PN time-domain 1SF scattering waveform, more solid than its ancillary-file placement suggests, but the per-mode Hadamard-scale cancellation is asserted rather than shown and a referee should demand it be demonstrated.

T0 · AI review

How this review was made: An AI referee (deepseek/deepseek-v4-flash) read the paper on September 25, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

T0 AI review An AI read the paper and answered a fixed list of questions.
T1 Audit Automatic checks that the references and files the paper relies on exist and match.
T2 Compute The paper's code or calculations re-run with the same results.
T3 Formal statement The claims written out in a formal math language a computer can read.
T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 2 major objections · 3 minor · Ask about this paper · The compiled record · Read the full review →

Challenge this review Comment (0) Read on arXiv ↗

Cite ▾

arXiv ID2609.30118 (version 2)

Pith IDK5RIBQJX, Pith's permanent record of this paper and its review

CitationAndrea Geralico. "Binary black hole scattering in the extreme-mass-ratio limit: time-domain waveform." arXiv:2609.30118 (2026). https://pith.science/paper/K5RIBQJX

Cite this reviewBibTeX for Pith's review of this paper

Claim paper

Why Pith says this

Highlight any text to comment on it.

Editor's note

read the letter

The paper is more solid than the ancillary-file placement makes it look, but the stress-test note names the one place a referee must actually dig: the Hadamard-scale cancellation in the hereditary modes is asserted, not shown, and the integrated flux checks cannot see a per-mode error.

What is actually new: the time-domain 1SF scattering waveform at 5PM/7PN for nonspinning bodies. The same author's frequency-domain result reaches 8PM/6PN, and amplitude-based time-domain computations sit at leading order (Kovacs-Thorne), so this fills a real gap. The time-domain form is genuinely useful for cross-checking independent calculations and for modeling the non-adiabatic burst near closest approach. The derivation is worked through with unusual care: Appendix A gives the distributional inverse Fourier transforms of ln^k omega in full through the parameter-free analytically continued kernel t^{-alpha-1}_+, and Appendix B presents the Hadamard machinery honestly, including the explicit sigma(y) scale in Eq. B15. Credit where due: energy and angular momentum losses derived from this waveform agree with independent amplitude-based results (Refs. 71 and 69), which is real external validation even if it enters through the author's own flux work.

The soft spot is the one the stress test names, and it is load-bearing rather than cosmetic. The polynomial-plus-log ansatz (Eq. 26) and the Hadamard finite part with arbitrary epsilon must leave each local mode W_lm(t_r) independent of epsilon. The paper asserts this but never demonstrates it. The printed hereditary quadrupole (Eq. 37) shows the physical combination ln[(y^2+1)^2/y] - ln(2 p_inf / b), but where the formal ln epsilon has been absorbed is not shown. The flux checks are time integrals over combinations built from this same waveform; a residual term that is a total derivative in retarded time, or orthogonal to those combinations, would pass both checks while leaving W_lm(t_r) wrong. That is a genuine referee question, not a nitpick. The encouraging part is that Appendix A supplies a concrete way to settle it: the Appendix B sigma-route must agree with the Appendix A route per mode, and that agreement should be exhibited for at least one representative (l, m) and one PM/PN order.

A separate, smaller problem: the headline deliverable, the full l=2..7 mode set, lives only in an unprovided ancillary file. The body shows essentially the l=2, m=2 modes, much of it to O(p_inf^2). That is a reproducibility gate that must come with the paper.

This deserves a serious referee and likely acceptance after revision. The core result looks right, the external flux checks are genuine, and the open questions are demonstrable rather than structural. Require the ancillary file, a per-mode demonstration of the epsilon cancellation (or the App A/App B agreement), and an explicit statement of where ln epsilon is absorbed. For the 1SF/PM scattering community this is a useful benchmark; I would want the regularization demonstration before leaning on the local waveform.

Referee report

Significance

If correct, this is a genuine 5PM/7PN time-domain 1SF scattering waveform, complementary to the frequency-domain MPM results and substantially beyond the current amplitude-based time-domain state of the art (tree level). The construction is built on the standard Teukolsky/MST foundation rather than fitted to the advertised result, and the integrated energy and angular-momentum losses pass independent amplitude-based checks (Refs. 71, 69), which is a real strength. The analytically continued kernel of App. A (Eqs. A8-A12) is a clean, parameter-free derivation of the needed distributions and is a useful reference result. The risk is isolated and specific: the local waveform modes inherit a regularization-scale dependence from the Hadamard prescription of App. B, and the paper never demonstrates that this scale drops out mode by mode; the ancillary-file delivery of the central claim compounds the difficulty of verification.

Major objections (2)

[Sec. IV, [b93]] The advertised deliverable, the complete l=2..7 mode set W_lm at 5PM/7PN, is contained only in an ancillary file: 'The explicit expressions of the rescaled waveform modes ... are given in an ancillary file up to O(G^4)' (Sec. IV, [b93]). The printed material in the manuscript covers only W_22^{(even)} at leading PN orders (Eqs. 35-37). The ancillary file was not available with the manuscript text I reviewed, so the central claim of the paper cannot be checked from the paper itself. The file contents must be supplied to the referees and representative modes (at least one higher-l and one odd-parity mode at full PN order) should be printed so that the claimed 5PM/7PN accuracy is substantiated in the published record.
[App. B, Eqs. B8-B15] The hereditary integrals are evaluated with a Hadamard finite part that depends on an arbitrary scale sigma(y)=epsilon f(y), with ln sigma(y) given explicitly in Eq. B15. For the local modes W_lm(t_r) to be well-defined the epsilon-dependence must cancel between W_inst and W_hered, but the paper neither states nor shows this cancellation: after B12 it notes only that epsilon is 'eventually replaced by a finite value playing the role of an arbitrary time scale', and the printed hereditary quadrupole in Eq. 37 exhibits the physical IR combination ln[(y^2+1)^2/y] - ln(2p_inf/b) with no indication of where the formal ln(epsilon) has been absorbed. Moreover, App. A supplies a parameter-free analytically continued kernel t_+^{-alpha-1} (Eqs. A8-A12) whose finite part carries no free scale; the App. B prescription agrees with it only if the sigma-dependence cancels exactly in W_inst + W_hered. Because the integrated flux checks are built from the same waveform via the author's own prior construction (Ref. [100]; [b54]), a residual epsilon-dependent term that is a total derivative in time, or orthogonal to the specific flux combinations used, would pass those checks while leaving the advertised local modes W_lm(t_r) wrong. A mode-by-mode demonstration of regularization-scale independence is required before the central claim can be accepted.

Minor objections (3)

[Ref. [15]] The reference list entry [15] contains garbled inline citation keys ('Favata:2010zu,Talbot:2018sgr,Yoo:2023spi,Elhashash:2025hqi') embedded in the text of the bibliographic entry; this formatting error should be repaired.
[Sec. I] In the sentence 'the PM approach is especially convenient because the bodies can be treated as asymptotically free... (see, e.g., Refs. [29-31] ad references therein)', 'ad' should read 'and'.
[Eqs. 36-37] The notation for the logarithm is inconsistent: Eq. 36 uses 'log(y)' throughout, while Eq. 37 mixes 'ln(...)' for the IR-scale combinations and 'log(y)' for the orbital logarithms; one convention should be used consistently, with a single definition of the time variable y (Eq. 12) restated at the point of first use in Sec. IV.

How Pith checks papers → · Full audit of this paper →

The paper

The paper claims

Gravitational bremsstrahlung from a small body scattering off a Schwarzschild black hole is written directly in the time domain, giving the complete first-order self-force waveform at fifth post-Minkowskian (three-loop) order and seventh…

The reading

This paper computes the gravitational waveform that a small body emits while scattering off a much larger Schwarzschild black hole, working to first order in the mass ratio and expressing the signal directly as a function of time. The result includes all radiative modes with l=2,...,7 and reaches fifth post-Minkowskian (three-loop) order with seventh post-Newtonian accuracy, with each mode split into an instantaneous part and a hereditary part that carries tail effects. Because a scattering event is intrinsically non-adiabatic, a time-domain representation is the natural one, and it complements the frequency-domain 1SF waveforms previously computed. The integrated energy and angular momentum losses derived from these modes already match independent amplitude-based calculations, so the waveform serves as a benchmark for other methods at the 1SF level.

What carries the argument

The argument runs through four linked objects. The first is the Weyl scalar $\psi_4$ of the Teukolsky formalism and its asymptotic amplitude $Z^\infty_{lm\omega}$, which carries all information on the emitted radiation. The second is the set of Mano–Suzuki–Takasugi (MST) homogeneous solutions to the radial Teukolsky equation used in the source factor $F_{lm\omega}(r_p)$, which — unlike pure post-Newtonian solutions — introduce powers of $\ln \omega$ alongside powers of $\omega$ into the frequency integral. The third is the structural ansatz for the integrand, $K_{lm\omega}(t') = \sum_{n,k\ge 1} [c_{lmn}(t') + c^{\ln}_{lmn}(t') \ln^k \omega]\, \omega^n$ (Eq. 26), with no other frequency dependence. The fourth is the inversion machinery: integration over $\omega^n$ produces Dirac-delta distributions and their derivatives (the instantaneous part), while the logarithmic kernels give the integrals $I^{\ln}_k(t)$ of Appendix A, whose divergent pieces at $t'=t_r$ are handled by Hadamard regularization (Appendix B). The convolution theorem separates the final mode into its instantaneous and hereditary parts.

Extended reading notes

Core claim

The central claim is that the full first-order self-force (1SF) time-domain scattering waveform can be obtained analytically by taking the frequency-domain amplitude $Z^\infty_{lm\omega}$ built from Mano–Suzuki–Takasugi (MST) solutions of the radial Teukolsky equation and inverting the Fourier transform mode by mode. After the post-Minkowskian and post-Newtonian expansion, the frequency integral's integrand reduces to a sum of terms of the form $\omega^n$ and $\omega^n \ln^k \omega$, so the inversion is carried by two elementary objects: Dirac-delta terms giving the instantaneous part of the waveform, and the tempered distributions $I^{\ln}_k(t)$ giving the hereditary part. The hereditary integrals diverge at the retarded time and are regularized with a Hadamard finite part that introduces an arbitrary infrared time scale; the construction requires that scale to cancel from the physical waveform modes. The explicit modes $W_{lm}$ for $l=2,\dots,7$ through $O(G^4)$ with 7PN accuracy are collected in an ancillary file.

Load-bearing premise

The whole result rests on the claim that, after the post-Minkowskian and post-Newtonian expansion, the integrand of the frequency integral is exactly a sum of integer powers of $\omega$ times powers of $\ln \omega$ (Eq. 26), with no other frequency dependence; if the MST radial solutions carried additional frequency structure at the attained order, the time-domain modes would inherit an error that the flux-level checks would not expose.

Editorial extensions

If this is right

The three-loop time-domain waveform at 1SF order provides a concrete target for other methods, since its integrated losses already agree with independent amplitude-based energy and angular momentum results.
It gives the first direct time-domain versus frequency-domain comparison for the same 1SF scattering observable at high post-Minkowskian order.
The mode decomposition $W_{lm}$, with its instantaneous/hereditary split and even/odd parity parts, isolates the contribution of tail and memory effects to the time-domain signal.
Because the result far exceeds the leading-order Kovacs–Thorne waveform, it constrains the current generation of amplitude-based time-domain computations.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

If the infrared scale really cancels as required, each physical mode is a clean prediction: one could numerically verify that $W_{lm}$ is independent of the regulator, a check stronger than the flux-level tests already passed.
A natural next check is a mode-by-mode comparison against the author's own 8PM frequency-domain waveform transformed numerically to the time domain, which would test the hereditary part directly rather than only its integrals.
The method's reliance on an ansatz with only integer powers of $\omega$ and powers of $\ln \omega$ suggests that entering 6PM order may require new machinery if non-integer or additional logarithmic structure appears in the radial solution.
The time-domain waveform is a natural input to effective-one-body constructions in the extreme-mass-ratio limit, where instantaneous and hereditary parts enter the accumulated phase of the signal.

Original abstract

The time-domain gravitational waveform emitted during the scattering of two nonspinning compact bodies is computed in the extreme-mass-ratio limit. It provides a time-domain description of the gravitational radiation emitted during the encounter and complements previous first-order self-force calculations, which have primarily been formulated in the frequency domain. In this way, the present results provide a complementary representation of the radiative dynamics and facilitate a direct comparison between time-domain and frequency-domain approaches to the gravitational-wave signal. The waveform is accurate to the fifth post-Minkowskian (three-loop) level and seventh post-Newtonian order, and will serve as a benchmark for future calculations by other methods, to first order in the mass ratio. The results have been already tested in a previous work by constructing the energy and angular momentum fluxes, thereby computing the corresponding radiative losses at the 5PM and 4PM level, respectively, with the same PN accuracy, which are in agreement with recent amplitude-based calculations.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13.0 via TypeSafe API, question bank abstract-p1-20260923, 2026-09-25: 1 request, 4,405 input tokens.

What the abstract itself states, as probabilities: claims first 0.07 · states limitation 0.53 · human participants 0.01 · controlled comparison 0.08 · randomized assignment 0.01 · states quantitative result 0.51 · compares with prior results 0.96 · mentions code or data availability 0.03.

The abstract, sentence by sentence

arXiv listing abstract v1. Results and the main claim are underlined; show every role.

The time-domain gravitational waveform emitted during the scattering of two nonspinning compact bodies is computed in the extreme-mass-ratio limit. It provides a time-domain description of the gravitational radiation emitted during the encounter and complements previous first-order self-force calculations, which have primarily been formulated in the frequency domain. In this way, the present results provide a complementary representation of the radiative dynamics and facilitate a direct comparison between time-domain and frequency-domain approaches to the gravitational-wave signal. The waveform is accurate to the fifth post-Minkowskian (three-loop) level and seventh post-Newtonian order, and will serve as a benchmark for future calculations by other methods, to first order in the mass ratio. The results have been already tested in a previous work by constructing the energy and angular momentum fluxes, thereby computing the corresponding radiative losses at the 5PM and 4PM level, respectively, with the same PN accuracy, which are in agreement with recent amplitude-based calculations.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#	sentence	role	the authors' own work	states the main conclusion	strength
s0	The time-domain gravitational waveform emitted during the scattering of two nonspinning compact bodies is computed in the extreme-mass-ratio limit.	contribution 0.93	0.94	0.48	asserted 1.00
s1	It provides a time-domain description of the gravitational radiation emitted during the encounter and complements previous first-order self-force calculations, which have primarily been formulated in the frequency domain.	contribution 0.60	0.75	0.17	asserted 1.00
s2	In this way, the present results provide a complementary representation of the radiative dynamics and facilitate a direct comparison between time-domain and frequency-domain approaches to the gravitational-wave signal.	interpretation 0.95	0.61	0.16	asserted 0.94
s3	The waveform is accurate to the fifth post-Minkowskian (three-loop) level and seventh post-Newtonian order, and will serve as a benchmark for future calculations by other methods, to first order in the mass ratio.	result 0.72	0.90	0.35	asserted 1.00
s4	The results have been already tested in a previous work by constructing the energy and angular momentum fluxes, thereby computing the corresponding radiative losses at the 5PM and 4PM level, respectively, with the same PN accuracy, which are in agreement with recent amplitude-based calculations.	result 0.62	0.63	0.09	asserted 1.00

What kind of paper the classifier reports this is: theoretical result 0.92 · evidence from simulation or computation 0.57.

What was examined, and what was not

sources read

arXiv listing abstract v1

not examined

full text
figures and tables
references and citing sentences
relations to other papers

levels run

P abstract

not run

S sections
T body sentences
C citations
R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

References

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

Analytical Time-Domain Waveforms for Black-Hole Scattering in the Self-Force Expansion
gr-qc 2026-10 conditional novelty 8.0 of 10

The time-domain gravitational waveform for hyperbolic black-hole scattering at first self-force order reduces at all PM and PN orders to one family of Lauricella-function integrals, computed here to 10PM-20PN.
Analytical Time-Domain Waveforms for Black-Hole Scattering in the Post-Minkowskian Expansion
gr-qc 2026-10 conditional novelty 8.0 of 10

Closed-form time-domain l=2 and l=3 waveform multipoles for one-post-Minkowskian black-hole scattering, with full velocity dependence, are derived for the first time.

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.