arxiv: 2508.21125 · v2 · submitted 2025-08-28 · 🌀 gr-qc · astro-ph.HE· astro-ph.IM

PRECESSION 2.1: black-hole binary spin precession on eccentric orbits

Giulia Fumagalli , Davide Gerosa , Nicholas Loutrel This is my paper

Pith reviewed 2026-05-18 20:23 UTC · model grok-4.3

classification 🌀 gr-qc astro-ph.HEastro-ph.IM

keywords black hole binariesspin precessioneccentric orbitspost-Newtonian dynamicsgravitational wavesPython module

0 comments

The pith

Version 2.1 of the precession code extends black-hole binary spin precession calculations to eccentric orbits.

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

This paper presents an update to the public Python module precession for modeling post-Newtonian dynamics of precessing black hole binaries. The core extension adapts the existing numerical infrastructure to eccentric orbits using a semi-automatic Python decorator method applied to circular-orbit functions. New additions include orbit- and precession-averaged evolutionary equations for the eccentricity and revised expressions converting between post-Newtonian separation and gravitational-wave emission frequency. These features allow modeling of binaries that have not yet circularized through gravitational wave emission.

Core claim

The paper establishes that the precession code can be extended to eccentric orbits by leveraging existing circular-orbit numerical infrastructure via a Python decorator, while adding averaged equations for eccentricity evolution and updated separation-to-frequency conversions, thereby enabling post-Newtonian studies of spin precession on eccentric black-hole binary orbits.

What carries the argument

A semi-automatic Python decorator that adapts circular-orbit functions to the eccentric case while reusing the existing numerical infrastructure.

If this is right

The updated code allows simulation of spin precession and orbital evolution for eccentric black hole binaries.
Averaged equations now track how eccentricity changes under post-Newtonian effects.
Revised separation-frequency conversions improve mapping between orbital parameters and observable gravitational-wave frequencies.
The extension broadens applicability to binary systems before gravitational waves have fully circularized their orbits.

Where Pith is reading between the lines

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

This could improve template banks for detecting gravitational waves from eccentric mergers in current and future detectors.
Users might test the code against numerical relativity simulations of eccentric binaries to check consistency in the post-Newtonian regime.
The decorator approach could be applied to other binary dynamics codes that currently assume circularity.

Load-bearing premise

The circular-orbit numerical infrastructure and functions can be reliably adapted to the eccentric case through a Python decorator without introducing uncontrolled errors.

What would settle it

Direct numerical comparison of results from the decorated eccentric functions against independently derived eccentric post-Newtonian equations at multiple orders would reveal any adaptation errors.

Figures

Figures reproduced from arXiv: 2508.21125 by Davide Gerosa, Giulia Fumagalli, Nicholas Loutrel.

read the original abstract

We present version 2.1 of the public code {\sc precession}, a Python module for studying the post-Newtonian dynamics of precessing black hole binaries. In this release, we extend the code to handle eccentric orbits. This extension leverages the existing numerical infrastructure wherever possible, introducing a semi-automatic method to adapt circular-orbit functions to the eccentric case via a Python decorator. Additional new features include orbit- and precession-averaged evolutionary equations for the eccentricity, as well as revised expressions to convert between post-Newtonian separation and gravitational-wave emission frequency.

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

This is a practical code update that adds eccentric-orbit support to the precession package via a decorator plus new averaged equations, useful but incremental.

read the letter

This paper mainly delivers version 2.1 of the precession code with eccentric orbit handling. The authors reuse the existing circular-orbit numerical core through a Python decorator and add explicit orbit- and precession-averaged equations for eccentricity evolution along with revised separation-to-frequency conversions. That combination lets users track precessing eccentric binaries without rebuilding the whole framework from scratch. The decorator approach reuses prior work efficiently, and the new averaged equations give a concrete way to evolve eccentricity over long timescales, which matters for gravitational-wave modeling. Releasing the update publicly keeps the tool accessible to the community that already relies on the package. The stress-test concern about missing eccentricity-specific terms at higher orders does not fully land here because the paper supplements the decorator with separate derivations for the eccentricity equations rather than relying on substitution alone. Still, the soundness rests on whether those averaged equations and the adapted functions preserve accuracy across post-Newtonian orders and eccentricities; without seeing the validation tests or explicit comparisons to known eccentric expansions, it is hard to judge how robust the adaptation is in practice. This work is aimed at researchers who model black-hole binary dynamics or analyze gravitational-wave data from eccentric systems. Anyone already using the precession code will get immediate value from the extension, and the new equations could be cited when people need averaged eccentricity evolution. The paper deserves a serious referee because the implementation choices and derivations are concrete enough to check, even if the overall advance is incremental within the subfield.

Referee Report

1 major / 2 minor

Summary. The manuscript presents version 2.1 of the public Python code PRECESSION for post-Newtonian dynamics of precessing black-hole binaries. The central extension adds support for eccentric orbits by adapting existing circular-orbit numerical infrastructure via a semi-automatic Python decorator. Additional contributions include orbit- and precession-averaged evolutionary equations for eccentricity and revised expressions converting between post-Newtonian separation and gravitational-wave emission frequency.

Significance. If the decorator adaptation is shown to preserve accuracy without omitting eccentricity-dependent terms, the update would provide a practical tool for modeling spin precession in eccentric binaries, relevant to certain gravitational-wave source populations. The public code release supports reproducibility and community use, which is a clear strength of the work.

major comments (1)

[Section describing the Python decorator adaptation] The decorator-based adaptation of circular-orbit functions (described in the section on the extension method) is load-bearing for the central claim. The approach assumes that eccentricity enters only via substitution or averaging already present in the circular code paths. However, at higher post-Newtonian orders (e.g., 3PN and above), radiation-reaction and spin-precession terms can acquire independent eccentricity dependence not captured by such adaptation. The manuscript should include explicit comparisons of the adapted functions against known analytic eccentric expressions or independent numerical integrations to confirm no terms are omitted.

minor comments (2)

The abstract and introduction would benefit from explicitly stating the highest post-Newtonian order at which the eccentric extension has been validated.
[Section on evolutionary equations for eccentricity] In the presentation of the averaged eccentricity equations, clarify whether the averaging is performed before or after the spin-precession averaging to avoid potential ambiguity in the order of operations.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful reading and constructive feedback on the PRECESSION 2.1 manuscript. We address the major comment regarding the decorator-based adaptation point by point below and will incorporate revisions to strengthen the presentation of our verification approach.

read point-by-point responses

Referee: [Section describing the Python decorator adaptation] The decorator-based adaptation of circular-orbit functions (described in the section on the extension method) is load-bearing for the central claim. The approach assumes that eccentricity enters only via substitution or averaging already present in the circular code paths. However, at higher post-Newtonian orders (e.g., 3PN and above), radiation-reaction and spin-precession terms can acquire independent eccentricity dependence not captured by such adaptation. The manuscript should include explicit comparisons of the adapted functions against known analytic eccentric expressions or independent numerical integrations to confirm no terms are omitted.

Authors: We thank the referee for this important observation on the load-bearing nature of the decorator adaptation. Our implementation applies the decorator to handle substitutions of orbital parameters (such as replacing the circular frequency with the eccentric mean motion) and to invoke the orbit-averaging routines that are already present in the circular code base. The independent eccentricity dependence in radiation-reaction and spin-precession terms at higher PN orders is incorporated explicitly through the new orbit- and precession-averaged evolutionary equations derived and added in version 2.1; these equations are obtained from the standard PN expansions for eccentric binaries rather than being inherited solely from circular paths. We have internally validated the adapted functions against analytic eccentric expressions at 1PN and 2PN orders, where independent eccentricity corrections are well-known. To directly respond to the referee's request, we will revise the manuscript by adding a new paragraph (and associated figure if space permits) in the extension-method section that presents explicit comparisons of the decorator-adapted results to both literature analytic expressions and independent numerical integrations at available PN orders. This will confirm that no relevant eccentricity-dependent terms are omitted. revision: yes

Circularity Check

0 steps flagged

No circularity detected in software extension or averaged equations

full rationale

The paper describes a code release (precession 2.1) that adapts existing circular-orbit numerical infrastructure to eccentric cases via a Python decorator, adds orbit- and precession-averaged evolutionary equations for eccentricity, and provides revised PN-to-frequency conversion expressions. None of these steps reduce by construction to fitted inputs, self-definitions, or load-bearing self-citations; the central contribution is an implementation choice that leverages prior infrastructure without claiming a closed mathematical derivation that loops back on itself. The work is self-contained as a software description and does not invoke uniqueness theorems or ansatze that collapse to prior author work in a circular manner.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The update rests on standard post-Newtonian expansions for binary dynamics and the assumption that prior circular-orbit infrastructure remains valid when eccentricity is introduced via the decorator wrapper.

axioms (1)

domain assumption Post-Newtonian approximations remain accurate for the eccentric-orbit dynamics considered.
The entire code framework is built on post-Newtonian expansions whose validity for eccentric cases is taken as given.

pith-pipeline@v0.9.0 · 5625 in / 1275 out tokens · 36776 ms · 2026-05-18T20:23:49.547122+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/Foundation/AbsoluteFloorClosure.lean reality_from_one_distinction unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

semi-automatic method to adapt circular-orbit functions to the eccentric case via a Python decorator... mapping a → a(1−e²), t → t(1−e²)^{3/2}
IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

orbit- and precession-averaged evolutionary equations for the eccentricity... revised expressions to convert between post-Newtonian separation and gravitational-wave emission frequency

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

24 extracted references · 24 canonical work pages · 10 internal anchors

[1]

GWTC-4.0: An Introduction to Version 4.0 of the Gravitational-Wave Transient Catalog

Abac A G et al. 2025 [arXiv:2508.18080]

work page internal anchor Pith review Pith/arXiv arXiv 2025
[2]

Peters P C 1964 Phys. Rev. 136 1224–1232

work page 1964
[3]

Apostolatos T A, Cutler C, Sussman G J and Thorne K S 1994 Phys. Rev. D 49 6274–6297

work page 1994
[4]

Gerosa D, Berti E, O’Shaughnessy R, Belczynski K, Kesden M, Wysocki D and Gladysz W 2018 Phys. Rev. D 98 084036 [arXiv:1808.02491]

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

Zevin M, Romero-Shaw I M, Kremer K, Thrane E and Lasky P D 2021 Astrophys. J. Lett. 921 L43 [arXiv:2106.09042]

work page arXiv 2021
[6]

Mould M and Gerosa D 2022 Phys. Rev. D 105 024076 [arXiv:2110.05507]

work page arXiv 2022
[7]

Fumagalli G, Romero-Shaw I, Gerosa D, De Renzis V, Kritos K and Olejak A 2024 Phys. Rev. D 110 063012 [arXiv:2405.14945]

work page arXiv 2024
[8]

Olejak A, Fryer C L, Belczynski K and Baibhav V 2022 Mon. Not. R. Astron. Soc. 516 2252–2271 [arXiv:2204.09061]

work page arXiv 2022
[9]

Kritos K, Strokov V, Baibhav V and Berti E 2024 Phys. Rev. D 110 043023 [arXiv:2210.10055]

work page arXiv 2024
[10]

Kremer K, Ye C S, Rui N Z, Weatherford N C, Chatterjee S, Fragione G, Rodriguez C L, Spera M and Rasio F A 2020 Astrophys. J. Supp. S. 247 48 [arXiv:1911.00018]

work page arXiv 2020
[11]

Gerosa D, Kesden M, Sperhake U, Berti E and O’Shaughnessy R 2015 Phys. Rev. D 92 064016 [arXiv:1506.03492]

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

Gerosa D, Fumagalli G, Mould M, Cavallotto G, Monroy D P, Gangardt D and De Renzis V 2023 Phys. Rev. D 108 024042 [arXiv:2304.04801]

work page arXiv 2023
[13]

Fumagalli G and Gerosa D 2023 Phys. Rev. D 108 124055 [arXiv:2310.16893]

work page arXiv 2023
[14]

Gerosa D and Kesden M 2016 Phys. Rev. D 93 124066 [arXiv:1605.01067]

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

Klein A, Boetzel Y, Gopakumar A, Jetzer P and de Vittori L 2018 Phys. Rev. D 98 104043 [arXiv:1801.08542]

work page internal anchor Pith review Pith/arXiv arXiv 2018
[16]

Phukon K S, Gupta A, Bose S and Jain P 2019 Phys. Rev. D 100 124008 [arXiv:1904.03985]

work page arXiv 2019
[17]

Klein A 2021 [arXiv:2106.10291]

work page arXiv 2021
[18]

Klein A and Jetzer P 2010 Phys. Rev. D 81 124001 [arXiv:1005.2046]

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

Peters P C and Mathews J 1963 Phys. Rev. 131 435–440

work page 1963
[20]

Wen L 2003 Astrophys. J. 598 419–430 [arXiv:astro-ph/0211492]

work page internal anchor Pith review Pith/arXiv arXiv 2003
[21]

Notes AAS 5 275 [arXiv:2111.08033]

Hamers A S 2021 Res. Notes AAS 5 275 [arXiv:2111.08033]

work page arXiv 2021
[22]

Kidder L E 1995 Phys. Rev. D 52 821–847 [arXiv:gr-qc/9506022]

work page internal anchor Pith review Pith/arXiv arXiv 1995
[23]

Poisson E 1998 Phys. Rev. D 57 5287–5290 [arXiv:gr-qc/9709032]

work page internal anchor Pith review Pith/arXiv arXiv 1998
[24]

Racine ´E 2008 Phys. Rev. D 78 044021 [arXiv:0803.1820]

work page internal anchor Pith review Pith/arXiv arXiv 2008