High-Post-Newtonian-Order Dynamics Induced by Tail-of-Tail Interactions: The Non-Geodesic Terms

Geraint Pratten

High-Post-Newtonian-Order Dynamics Induced by Tail-of-Tail Interactions: The Non-Geodesic Terms

T0 review · 2 major / 4 minor · reviewed 2026-08-04 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read Eccentric binary dynamics pinned to 12th order in eccentricity

desk verdict A high-order PN computation with genuinely new Q-potential and 2SF redshift coefficients; the one structural risk is an unproven radial-momentum relation (Eq. 61) that all new EOB coefficients depend on. read the letter →

arxiv 2608.01774 v2 pith:XEVFSGSI submitted 2026-08-03 gr-qc

Geraint Pratten This is my paper

classification gr-qc MSC 83C2583C35 PACS 04.25.Nx04.30.-w

keywords tail-of-tailinteractionhereditaryeffectseccentricbinariespost-Newtonianexpansioneffective-one-bodyself-forceredshiftinvariantgravitationalwaves

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper aims to settle what gravitational-wave tail-of-tail interactions do to the conservative dynamics of eccentric, non-spinning compact binaries at the next correction beyond the first hereditary term. It extends the Delaunay-averaged tail-of-tail Hamiltonian from O(e_t^2) to O(e_t^12), then transcribes that Hamiltonian into the effective-one-body description. The result is a complete set of 5.5PN and 6.5PN coefficients for the non-geodesic Q potential through the p_r^12 sector, including the full dependence on the symmetric mass ratio. The linear-in-mass-ratio terms reproduce known first-order self-force results; the quadratic terms are new second-order self-force predictions. The same Hamiltonian, fed through the first law of binary mechanics at fixed orbital frequencies, yields the complete tail-of-tail contribution to the second-order inverse redshift at those PN orders through O(e^10).

What carries the argument

The central object is the time-split tail-of-tail action: a nonlocal-in-time integral built from products of multipole moments evaluated at advanced and retarded times, with gravitational-wave backscattering encoded in principal-value kernels. The calculation wraps this action in a 1PN quasi-Keplerian parametrisation of the harmonic-coordinate orbit, Delaunay-averages it (orbital average over one radial period), and matches the averaged Hamiltonian to the effective-one-body Hamiltonian order by order in eccentricity. Because the non-geodesic Q potential starts at quartic order in radial momentum, an eccentricity expansion through O(e_t^12) unlocks Q coefficients through p_r^12. A Fourier-Bes

What would settle it

Compute the second-order inverse redshift for eccentric orbits about a Schwarzschild black hole at fixed m2-rescaled frequencies, keeping terms through O(e^10) at 5.5PN and 6.5PN; the tail-of-tail coefficients in Eqs. (183)-(184) should appear unchanged. A cheaper check: rederive the O(e^2) second-order redshift from an independent effective-one-body Hamiltonian, extending the internal consistency check already reported through O(e^2).

Watch

Extended reading notes

Core claim

The paper derives the conservative tail-of-tail Hamiltonian for eccentric, non-spinning compact binaries to relative first post-Newtonian order and through O(e_t^12) in eccentricity, then matches it to the effective-one-body description. From the Delaunay-averaged Hamiltonian it extracts the 5.5PN and 6.5PN corrections to the non-geodesic potential Q(u,p_r;nu) through the p_r^12 sector, including the complete dependence on the symmetric mass ratio through O(nu^2). The linear-in-nu coefficients reproduce established first-order self-force results; the quadratic-in-nu coefficients are new eccentric second-order self-force predictions. Applied to the first law of binary mechanics at fixed frequ

Load-bearing premise

The derivation assumes that 1PN quasi-Keplerian motion in harmonic coordinates, together with the quadrupole (to 1PN), octupole, and current-quadrupole (to Newtonian order) source moments, captures every tail-of-tail contribution at relative 1PN order through O(e^12).

Editorial extensions

If this is right

Eccentric effective-one-body models can now carry analytic tail-of-tail Q-potential coefficients through p_r^12 at 5.5PN and 6.5PN.
The second-order inverse-redshift coefficients through O(e^10) provide fixed-frequency weak-field benchmarks that direct eccentric second-order self-force computations can target.
The O(nu^2) parts of the Q coefficients are qualitatively new: they encode mass-ratio dependence of the tail-of-tail interaction at second self-force order.
The spectral-norm identity shows each averaged hereditary bilinear equals minus a positive-definite sum over Fourier amplitudes, fixing signs term by term.
The matching procedure is iterative and can be carried to higher eccentricity order, extending the Q potential to higher p_r powers with the same inputs.

Reading between the lines

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

The same EOB matching route, applied to unbound orbits by analytic continuation, should yield the tail-of-tail part of the scattering angle; the paper lists this as future work, but the ingredients appear already present.
If a direct eccentric second-order self-force code reaches O(e^10), disagreement with Eqs. (183)-(184) would most likely trace to the 1PN quasi-Keplerian truncation or to omitted multipole sectors, since the first-order and spectral checks already close those channels at lower order.
Going beyond relative 1PN order would require 2PN orbital dynamics and higher multipole moments, so the present coefficient list is the natural stopping point of this harmonic-coordinate, fixed-PN-order approach.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

Geraint's paper is a careful, high-order PN calculation with real new results: the tail-of-tail Hamiltonian through O(e_t^12), the EOB Q-potential coefficients through p_r^12 including the full O(nu^2) dependence, and 5.5PN/6.5PN 2SF inverse-redshift predictions through O(e^10). The time-domain averages are cross-checked against a Fourier–Bessel spectral calculation, the 1SF limit reproduces the Munna–Evans results through O(e^12), and the EOB consistency check works through O(e^2). That is genuinely solid work and the coefficients are new.

The soft spot is exactly the one the stress-test note flags. The new Q coefficients all rest on Eq. (61), the corrected 1PN radial-momentum relation in EOB coordinates. The paper states it follows from the 1PN EOB energy-conservation equation and fixes a typo in Ref. [87], but no derivation is shown. The 1SF redshift check bypasses EOB matching, and the internal EOB check uses the same Eq. (61), so it does not independently test that relation. If Eq. (61) is wrong, every new q_{2m,n} shifts, including the O(nu^2) parts. This is a real risk, not a manufactured one.

The paper is honest about what it does not show: the full O(e^12) algebra is in the ancillary material, and the 2SF order is limited to O(e^10) because of the derivatives. Those are limitations, not fatal flaws. The claimed typo correction is exactly the kind of thing a referee should ask to see derived.

I would send this to peer review. It deserves a serious referee. The referee should ask for the derivation of Eq. (61) and, ideally, an independent verification of the higher-order Q coefficients, perhaps by extending the spectral method or comparing the 2SF redshift predictions to a direct low-eccentricity calculation. The paper is important for eccentric EOB modeling and 2SF benchmarks; it is not a desk reject.

Referee Report

2 major / 4 minor

Summary. The paper computes the conservative tail-of-tail dynamics of eccentric, non-spinning compact binaries to relative 1PN order and through O(e_t^12). In harmonic coordinates, it derives the Delaunay-averaged tail-of-tail Hamiltonian, matches it to the effective-one-body (EOB) description, and thereby extracts the 5.5PN and 6.5PN non-geodesic Q-potential coefficients through O(p_r^12) with complete dependence on the symmetric mass ratio through O(nu^2). It also derives the tail-of-tail contribution to the 2SF inverse redshift at those PN orders through O(e^10) using the first law at fixed orbital frequencies. Three validation routes are presented: a spectral Fourier–Bessel rederivation of the averaged Hamiltonian, a 1SF redshift comparison against the independent Munna–Evans calculation, and an EOB self-consistency check through O(e^2).

Significance. If the results are correct, this is a substantial step: it provides the first eccentric 2SF predictions for the non-geodesic EOB potential at 5.5PN and 6.5PN, extends the Q-potential to p_r^12, and gives weak-field benchmarks for ongoing second-order self-force programs. The paper is largely derived rather than fitted: there are no free parameters, the spectral and time-domain routes agree, and the 1SF comparison is claimed to agree through O(e^12) with an independent calculation. However, the central novelty rests on a single asserted relation for the 1PN EOB radial momentum, Eq. (61), which is not derived and is not tested independently at the orders needed. The paper's overall value is high, but the soundness of the new Q-potential coefficients depends on closing that gap.

major comments (2)

[IV.B, Eqs. (70)-(79)] The corrected 1PN radial-momentum relation p_r^(0) = p_r^N [1 + eta^2(1/r - 1/(2 a_e_r))] is stated without derivation, the text noting only that it corrects a typo in Ref. [87]. This relation is load-bearing: every new Q-potential coefficient in Eqs. (70)-(79) and the 2SF inverse-redshift coefficients in Eqs. (183)-(184) flow from the EOB matching that uses Eq. (61). The internal cross-check in Sec. VI.C is not independent because it evaluates the same EOB Hamiltonian with the same Eq. (61), and it is truncated at O(e_t^2); the 1SF comparison in Sec. VI.A bypasses the EOB matching entirely. I request a derivation of Eq. (61) from the 1PN EOB energy-conservation equation, or an external check (e.g., a direct evaluation of the EOB Hamiltonian along the 1PN QK orbit to O(e_t^12), or comparison with an independent eccentric 2SF calculation).
[IV.B, Eqs. (70)-(79)] The matching procedure that converts the harmonic Delaunay averages (45)-(46) into the Q-potential coefficients (70)-(79) is not shown. The text states that the expansion is iterative and that Q first enters at O(e_t^4), but the explicit linear system at each eccentricity order is absent. Because these coefficients are the central new result, the reader cannot verify the uniqueness of the extraction or the absence of cross-contamination from the A- and D-potential sectors. Please provide the matching equations (or a detailed algorithm) and, ideally, a table showing which eccentricity orders fix each q_{2m,n}. The O(e_t^12) algebra may remain in ancillary material, but the structure of the matching should be in the main text.

minor comments (4)

[VI.A] The agreement with Munna–Evans is stated as exact to all orders in eccentricity at 5.5PN and coefficient-by-coefficient at 6.5PN. A small table comparing the eccentricity expansions of U^(1)_11/2 and U^(1)_13/2 would make this check more transparent.
[V] The spectral cross-validation is described as independent at the level of the hereditary integral and orbital average, but it uses the same source multipoles and QK dynamics. This should be stated prominently as a consistency check rather than an independent physical derivation, which the text already acknowledges.
[III] The complete O(e_t^12) multipole bilinears are relegated to ancillary material, while Tables I-III show only O(e_t^4). Please ensure the ancillary material is permanently archived and that a version is included with the published paper, as the final Hamiltonian and Q coefficients cannot be checked otherwise.
[Abstract] The phrase 'complete O(nu^2) dependence' should be qualified in the abstract: the Q-potential result is through O(p_r^12), and the 2SF inverse-redshift result is through O(e^10). The current wording could overstate the eccentricity range of the redshift invariants.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the Q-potential and 2SF redshifts are derived from an independently constructed harmonic tail-of-tail Hamiltonian and transcribed to EOB; nothing is fitted or reduced to its own input.

full rationale

The paper's central derivation is self-contained in the relevant sense: the Delaunay-averaged tail-of-tail Hamiltonian in Eqs. (45)-(46) is built from standard MPM multipoles and beta coefficients (Sec. III), and the EOB matching in Sec. IV is an algebraic transcription between two representations of the same conservative dynamics. No parameter is fitted to the target predictions. The Q-potential coefficients (70)-(79) are obtained by solving the matching equations order-by-order in eccentricity, with the circular and e_t^2 sectors fixing A and \bar{D} and the higher sectors determining Q; this is an invertible map rather than a definition of Q in terms of itself. The 2SF inverse-redshift results in Eqs. (183)-(184) are derived from the same Hamiltonian through the first law at fixed frequencies, a standard variational identity, and are not fitted. The O(\nu) Q coefficients are explicitly checked against independent GSF results from Refs. [86,88], and the 1SF redshift is checked against the independent calculation of Ref. [98]. The Fourier-Bessel spectral rederivation in Sec. V uses the same source multipoles but evaluates the hereditary integral by a distinct spectral norm, so it is a consistency check rather than circularity. The manuscript's reliance on Eq. (61), the corrected 1PN EOB radial-momentum relation, is a stated input and is load-bearing, but it is not defined in terms of the derived Q coefficients or redshift invariants; a possible error there would be a correctness/verification issue, not circularity. The paper's self-citations (e.g., Pratten et al. waveform papers) are contextual and not load-bearing for the tail-of-tail derivation. No circular step could be exhibited, so the circularity score is 0.

Assumptions & free parameters 0 free parameters · 8 assumptions · 0 invented entities

The paper's contribution is a derived expansion, not a model fit: there are no fitted constants and no new entities are introduced. The central calculation leans on a standard scaffold of 1PN QK dynamics, MPM multipoles, the time-split tail-of-tail action, EOB/DJS gauge, the hereditary first law, and Schwarzschild-geodesic frequency maps, all taken from the cited literature. The only free choices are truncation orders (e^12, relative 1PN), which the paper states are computational.

assumptions (8)

domain assumption The 1PN quasi-Keplerian parameterization in harmonic coordinates (Eqs. 1-3), including the mean-motion relation Eq. (7) and the (epsilon,j)-to-orbital-element map (Eqs. 9-17).
The orbital scaffold for all Delaunay averages and matching; taken from Refs. [107-110] as established PN results.
domain assumption The time-split tail-of-tail action Eq. (18) with the beta coefficients Eq. (25) and principal-value convention is the correct conservative hereditary interaction.
Adopted from Refs. [18,86,87] and [9,77,79,113-115]; the paper states conventions but does not re-derive the action from first principles.
domain assumption The harmonic MPM multipole moments I2, I3, J2 at the stated PN accuracy (Eq. 26) are complete.
Standard moments from Refs. [1,116]; completeness at relative 1PN is presupposed.
domain assumption The first law of binary mechanics retains its standard form for the localized Delaunay-averaged hereditary Hamiltonian.
Invoked in Sec. VI, attributed to Ref. [97]; this is the basis for the redshift derivations.
domain assumption The DJS-gauge EOB Hamiltonian (Eq. 50) with A, Dbar, Q expansions (Eqs. 52-54) fully represents the dynamics at the PN orders considered.
Framework from Ref. [89]; the Q extraction depends on this gauge choice.
standard math The spectral identities used in Sec. V: the principal-value relation Eq. (90), the diagonal action on Fourier modes, and the Kapteyn-series convergence of Hansen coefficients.
Background for the Fourier-Bessel cross-validation; standard results from Refs. [78,92,121-125].
domain assumption The reference Schwarzschild geodesic maps and O(q) frequency-shift relations, including the Darwin parameters (Eqs. 148-150, 165-168, 180-181).
Input from Refs. [96,133,134,139] used for fixed-frequency GSF comparisons.
domain assumption The 2PN periastron advance (Eq. 15) enters the 2SF redshift computation.
From Refs. [109,110]; needed for the O(epsilon) term in delta_q j, Eq. (166).

how reviews work

0 comments

Cite this review

Pith. "Pith review of High-Post-Newtonian-Order Dynamics Induced by Tail-of-Tail Interactions: The Non-Geodesic Terms." pith.science (2026). https://pith.science/paper/XEVFSGSI

@misc{pith2026260801774,
  author       = {Pith},
  title        = {Pith review of: High-Post-Newtonian-Order Dynamics Induced by Tail-of-Tail Interactions: The Non-Geodesic Terms},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/XEVFSGSI}},
  note         = {Machine review of arXiv:2608.01774}
}

abstract

We compute the tail-of-tail contribution to the conservative dynamics of eccentric, non-spinning compact binaries to relative 1PN order and to $\mathcal O(e_t^{12})$. Using the $1$PN quasi-Keplerian dynamics in harmonic coordinates, we derive the Delaunay-averaged Hamiltonian at $5.5$PN and $6.5$PN order and match it to the effective-one-body description, allowing us to determine the corresponding contributions to the non-geodesic EOB $Q$ potential through the $\mathcal{O}(p_r^{12})$, including the dependence on the symmetric mass ratio up to $\mathcal{O}(\nu^2)$. The terms linear in the mass ratio reproduce the available first-order self-force results, while the quadratic terms provide qualitatively new eccentric second-order self-force predictions arising from the tail-of-tail terms. We independently rederive the averaged Hamiltonian using a Fourier--Bessel decomposition of the hereditary interaction. Applying the first law at fixed orbital frequencies, we recover the known first-order self-force redshift through $\mathcal O(e^{12})$ and obtain the complete tail-of-tail contribution to the second-order inverse redshift at $5.5$PN and $6.5$PN through $\mathcal O(e^{10})$.

Discussion (0). Sign in to comment.

Pith tools

Reviewed August 4, 2026 · model on record in the stance chip above.