{"id":"0ada9fe3-eeb2-4be6-a96b-65764b8d7d89","arxiv_id":"2412.10909","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Closed-form post-Newtonian phasing formulas for eccentric, spin-aligned compact binary inspirals are derived to 3PN order and to eighth order in initial eccentricity, with a resummation extending validity to e0 around 0.55.","lead":"This paper derives closed-form formulas for the gravitational-wave phase of eccentric, spinning binary inspirals up to the third post-Newtonian order, expanding to the eighth power in initial eccentricity. The formulas enable faster waveform generation and show that neglecting spin-eccentricity coupling can cause mismatches above 1% for eccentricities above about 0.15, which matters for detecting such systems.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The O(e8) phasing coefficients reside only in an unverified supplemental Mathematica file, and the paper's only numerical check is against the same truncated equations that produced them, so the central claim rests on unchecked algebra.","rationale":"The reader correctly located the weakest assumption in the inherited inputs and the self-referential numerical check. I sharpen this: the unshown O(e^8) expansion is the load-bearing piece because the paper's own text says the full expressions are in the supplement, and the internal comparison cannot distinguish an error in the perturbative solution from a correct one. I credit the authors for clearly stating this limitation and for providing the supplement, but the absence of independent verification leaves the central claim conditional, not rejected. A coefficient-by-coefficient symbolic reproduction is the decisive, low-cost test. The reader's conditional verdict should stand.","tokens_in":58157,"tokens_out":6528,"duration_ms":64420,"concrete_test":"Independently re-derive the O(e0^8) TaylorT2 phase from the published inputs: take Eqs. (3.18)–(3.20) of Ref. [73], Taylor-expand dy/dt and de^2/dt to 3PN and O(e0^8), change variable to y, solve de^2/dy perturbatively to O(e0^8), integrate dt/dy and d⟨ϕ⟩/dy, and compare every coefficient with the supplemental file [78]. If any coefficient differs, the distributed formulas are wrong; exact agreement would close the main verification gap.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The manuscript's headline results — the 3PN, O(e0^8) TaylorT2 and TaylorF2 phases — are not printed in the paper. Equations (20) and (22) display only the leading O(e0^2) eccentric-spin terms; the full O(e0^8) expressions are in a Mathematica supplement [78]. The paper's Fig. 1 validation compares the analytic O(e0^8) phase with numerical integration of the same orbit-averaged equations (4a) and (4b) from which the analytic expansion was derived. That comparison is circular at the level of the new terms: a coefficient error in the perturbative solution of de^2/dy (Eqs. 6 and 8) or in the subsequent integrations would appear in both the analytic phase and the numerical reference, so agreement cannot certify correctness. The same applies to the mismatch study: it compares TaylorF2 with and without the new spin-eccentric terms, but both sides inherit the inputs of Ref. [73] and the authors' expansion algebra. Thus the central claim — valid closed-form eccentric spinning phasing to 3PN and O(e^8) — is exactly as strong as unshown algebra in a supplement. This is not an accusation of error; it is a load-bearing verification gap, because a single incorrect term in the supplemental expressions would invalidate the new terms that the paper's conclusions depend on.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper derives post-Newtonian phasing formulas for eccentric, spin-aligned compact binary inspirals up to 3PN order and to O(e0^8) in the initial eccentricity. The authors use the energy and flux inputs of Ref. [73], construct TaylorT2 and TaylorF2 approximants, and also provide a resummed TaylorT2 phase intended to remain accurate up to e0 ~ 0.5. They quantify the importance of the new eccentric spin-orbit and spin-spin terms through cycle-count estimates and mismatch studies against TaylorF2Ecc, reporting mismatches above 1% for eccentricities ≳0.15 and moderate spins. Only the O(e0^2) eccentric spin terms are printed in the main text; the full O(e0^8) expressions are relegated to a supplemental Mathematica file.","tokens_in":58455,"tokens_out":10672,"duration_ms":102893,"significance":"If the supplemental algebraic results are correct, this is a useful and timely contribution: it would provide the first closed-form, eccentric, spinning 3PN TaylorT2 and TaylorF2 phasing expressions, enabling fast template generation for eccentric spinning binaries. The derivation follows standard PN and stationary-phase techniques, uses published inputs from Ref. [73], and the displayed O(e0^2) coefficients appear structurally consistent with the non-spinning limit of Moore et al. The paper also gives a concrete, phase-based assessment of when spin-eccentricity coupling matters for detection. The main caveats are verification-related: the headline O(e0^8) expressions are not inspectable from the manuscript, and the numerical validation uses the same underlying orbit-averaged equations from which the analytic expansion was built, so it cannot independently certify the input PN coefficients or the orbit-averaging approximation.","major_comments":[{"comment":"The central claim of the paper, namely closed-form 3PN, O(e0^8) TaylorT2 and TaylorF2 phasing, is not checkable from the manuscript because Eqs. (20) and (22) display only the leading O(e0^2) eccentric spin-orbit and spin-spin terms. The paper states that the full O(e0^8) expressions appear in a GitHub supplement, and Tables II and IV, Fig. 1, and the conclusions all depend on those unprinted coefficients. A single algebraic error in the supplement would invalidate the quantitative conclusions. The authors should make the supplement a permanent, versioned part of the manuscript and provide at least one explicit O(e0^4) or O(e0^6) coefficient in an appendix, together with a symbolic verification (for example, recomputation with an independent expansion/integration code, or a check that the χ→0 limit exactly reproduces Eq. (6.26) of Ref. [77]).","section":"Section II.B, Eqs. (20)-(23) and Supplemental Material [78]"},{"comment":"The numerical reference in Fig. 1 is obtained by integrating the same orbit-averaged PN equations (4a)-(4b) that are the starting point of the analytic expansion. This comparison can detect errors in the perturbative solution of de^2/dy and in the subsequent integrations, but it cannot certify the underlying 3PN energy and flux inputs of Ref. [73], nor the neglect of oscillatory phase contributions. The text should state this limitation explicitly and add an independent check, for example a comparison with direct numerical integration of the unexpanded equations of Ref. [73] or with an independent non-spinning eccentric phasing code in the χ→0 limit, so that Fig. 1 is not the sole evidence for the new terms.","section":"Section III.B, Fig. 1 and Eq. (29)"},{"comment":"The exponent '3' in the resummation ansatz (1-e0^2)^3 is chosen by trial and error against the same numerical solution used for validation, and the claimed extension of validity to e0 ≈ 0.55 is therefore a heuristic result rather than a derived property. The abstract and conclusions should present it as such. Moreover, Fig. 1 demonstrates the resummation for only one system (1.4 + 1.4 M⊙, χ1 = 0.7, χ2 = 0.8); the validity claims for the mass and spin configurations used in Tables II and IV should be checked explicitly rather than assumed to carry over.","section":"Section II.B.3, Eqs. (24)-(26)"}],"minor_comments":[{"comment":"There is a typo in 'the the orbital phase (ϕ) includes oscillatory contributions'; delete the duplicated 'the'.","section":"Section II.A, first paragraph"},{"comment":"The notation is inconsistent in several displayed formulas: χS and χs, κS and κs, and χA and χa are used interchangeably. Please unify the notation throughout.","section":"Equations (21d), (23b), (C10e)"},{"comment":"The construction of Eq. (25) is opaque: the explicit factor 3e0^2 in front of ϕ(resum)_SO,ecc and ϕ(resum)_SS,ecc changes the relation between these coefficients and those in Eq. (21), and it appears to introduce O(e0^2) corrections to circular spin terms. Please clarify how Eq. (25) follows from the matching ansatz Eq. (24).","section":"Eq. (25)"},{"comment":"Entries below 10^-3 are omitted without explicit zeros, leaving apparent gaps in the table; fill those entries with 0.000 or explain the omission in the caption.","section":"Table III"},{"comment":"The phrase 'trail and error method' should read 'trial and error method'.","section":"Section I.A and Sec. II.B.3"},{"comment":"The supplemental file is a GitHub URL without a version or commit hash. Since the central result of the paper resides in that file, please provide a persistent archived version, such as a DOI or a journal-hosted supplement, and include a commit identifier.","section":"Reference [78]"},{"comment":"The mismatch study does not specify which amplitude model is used in the SPA waveforms, nor whether higher harmonics are included. Since the conclusions concerning detection efficiency depend on the full waveform, please state explicitly that only the phasing is varied and the amplitude is taken from TaylorF2Ecc.","section":"Section III.C, Figs. 2-3"}],"recommendation":"major_revision","confidential_remarks":"I am not recommending rejection: the derivation follows standard and well-understood methods, and the displayed O(e0^2) terms look structurally sound. The verification gap is real but fixable. The 'circularity' objection should not be overstated: a direct numerical integration of the same differential equations would catch errors in the e-expansion algebra, so Fig. 1 is meaningful for the perturbative series. What it cannot do is validate the input 3PN expressions from Ref. [73] or the orbit-averaging approximation. The editor should require that the supplement be archived and independently checked, and that the authors add at least one independent validation of the new terms before acceptance."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper does what it says: it gives closed-form TaylorT2 and TaylorF2 phasing for eccentric, aligned-spin binaries through 3PN, expanded to e0^8, plus a resummed version that extends the range to e0 around 0.5. That specific combination is new and useful for fast template generation. The derivation is standard: it takes the energy and flux from Henry & Khalil, re-expresses in the y variable, solves the eccentricity evolution perturbatively, and integrates. The cycle estimates and mismatch studies do a good job showing that ignoring spin-eccentricity coupling can cost more than 1% in match for e0 around 0.15 and spins around 0.2. That is a real result, not a trivial repackaging.\n\nThere are three soft spots, in rough order of importance. First, the central new content — the full O(e8) coefficients — is not in the paper. The main text displays only O(e2) spin-eccentric terms; equations (20) and (22) are placeholders for the supplement. That is common in this field, but the verification is the issue. Fig. 1 compares the analytic e8-expanded phase with numerical integration of Eqs. (4a)–(4b), which are the same truncated orbit-averaged equations the expansion was derived from. This catches coding slips and some consistency errors, but not a wrong coefficient that appears on both sides. A single bad term in the supplemental Mathematica file would invalidate the paper's main claim. I don't suspect error; I just note the verification is weaker than it looks. Second, the abstract overstates novelty: \"no closed-form expression ... both spin and eccentricity\" is contradicted by existing eccentric spinning models, including Ref. [50] and others cited in the paper. The specific e8, 3PN, aligned-spin phasing is new; the blanket claim is not. Third, the resummation exponent is chosen by trial and error and matched to the paper's own PN series. That is fine as a practical tool, but it is not externally calibrated.\n\nThe inputs from Ref. [73] are external and published, so the derivation inherits their quality, which is good. Citation practice is otherwise normal and self-citation is minor. The paper deserves a serious referee. I would send it out, asking the referee to spot-check the supplemental algebra, ideally by comparing a few high-order coefficients against an independent computation or at least a different integration scheme, and to require the authors to soften the novelty claim. If the supplement checks out, this is a solid, citable contribution to waveform modeling.","headline":"Genuinely useful 3PN aligned-spin eccentric phasing, but the headline O(e8) terms live only in a supplement and the only check is against the same equations that produced them; worth refereeing with a demand for independent verification.","tokens_in":58979,"tokens_out":2127,"would_cite":true,"duration_ms":25466,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["83C25","83C35"],"pacs":["04.30.-w","04.25.Nx","04.30.Tv"],"model":"deepseek-v4-flash","headline":"This paper provides the first closed-form gravitational-wave phasing for eccentric, spin-aligned binaries, to 3PN order and eighth power in eccentricity, in both time and frequency domains.","keywords":["gravitational wave phasing","eccentric compact binaries","aligned spins","post-Newtonian approximation","Taylor approximants","stationary phase approximation","spin-orbit coupling","spin-spin coupling"],"falsifier":"Decisive check: numerically integrate the full 3PN equations of motion for an aligned-spin eccentric binary with initial eccentricity 0.5 without discarding the oscillatory phase pieces, and compare the accumulated number of gravitational-wave cycles with the TaylorT2 formula and its resummed version; if the difference exceeds one cycle inside the claimed validity range, the orbit-averaged approximation underpinning the closed forms is inadequate.","tokens_in":1757,"feed_emoji":"🌌","tokens_out":3957,"duration_ms":128961,"temperature":0.7,"pith_summary":"The paper derives analytic, closed-form expressions for the gravitational-wave phase of inspiralling compact binaries that are eccentric and have spins aligned with the orbital angular momentum. These are the first such formulas to include both effects together, reaching the third post-Newtonian order in the spin terms and the eighth power in the initial eccentricity $e_0^8$. The time-domain result (TaylorT2) and its stationary-phase frequency-domain version (TaylorF2) allow fast template evaluation, which matters because numerical evolution of eccentric orbits is slow. Comparing the new TaylorF2 terms with the existing eccentric-only model, the paper reports mismatches above 1% once the initial eccentricity exceeds about 0.15 and spins reach about 0.2, with mismatches up to roughly 15% in parts of the parameter space. The implication is that template banks neglecting the combined spin-eccentricity coupling will lose signal-to-noise ratio and detection efficiency.","feed_headline":"First closed-form phasing for eccentric spinning binaries","feed_subtitle":"Closed-form TaylorT2 and TaylorF2 inspiral phases to 3PN capture spin-eccentric couplings that cost over 1% in match.","key_machinery":"The central object is the closed-form TaylorT2 phase, a post-Newtonian series in the frequency parameter $y = (x/(1-e^2))^{1/2}$, with the eccentricity replaced by its perturbative solution $e^2(y) = e_0^2 (y_0/y)^{19/3}$ times spin-dependent corrections. The machinery has four parts: taking the 3PN energy and energy flux for aligned-spin eccentric binaries as input; re-expanding them in the $y$ parameter; solving the eccentricity evolution equation order by order in $e_0$; then inverting the $dy/dt$ evolution and integrating to get $t(y)$ and the orbital phase. Applying the stationary phase approximation to the TaylorT2 phase produces the TaylorF2 frequency-domain phase, and an empirical resummation of the form $\\phi = y^{-5}(1-e_0^2)^3$ times a polynomial in $e_0$ extends the time-domain result to higher eccentricities.","core_discovery":"The central claim is that the gravitational-wave phasing of an inspiralling compact binary with aligned spins and small eccentricity can be written as a fully analytic double series in the post-Newtonian parameter and the initial eccentricity $e_0$, with no numerical integration needed. The paper constructs the TaylorT2 phase and its TaylorF2 Fourier counterpart by perturbatively solving the coupled evolution of the frequency parameter $y$ and the eccentricity $e$, using the 3PN energy and flux of spinning eccentric binaries as input. The spin contributions split into spin-orbit and spin-spin pieces at the expected PN orders (1.5PN, 2PN, 2.5PN and 3PN), and the eccentricity expansion is carried to $O(e_0^8)$. A resummed TaylorT2, built from a $(1-e_0^2)^3$ ansatz, extends the range of validity to initial eccentricities near 0.55. The paper further reports that the newly computed eccentric-spinning terms in TaylorF2 change the match by more than 1% once the initial eccentricity is about 0.15 and spins are about 0.2, so neglecting them degrades searches and parameter estimation.","pith_inferences":["The same perturbative eccentricity-expansion strategy could be applied to waveform amplitudes rather than only the phase, yielding a fully analytic eccentric spinning waveform family; this paper stops at phasing.","Because mismatches grow with effective spin and chirp mass, third-generation and space-based detectors, which accumulate more cycles, will likely need these spin-eccentricity terms even for initial eccentricities below 0.15.","The resummation ansatz was chosen empirically, so a more systematic resummation scheme might extend validity beyond about 0.55, but that is a conjecture beyond the paper.","A direct cross-check of these formulas against independent numerical waveforms for initial eccentricities near 0.3 to 0.5 would test whether the orbit-averaged approximation remains adequate at the upper edge of the claimed validity range."],"forward_implications":["With these closed forms, an eccentric spinning inspiral template can be evaluated without numerical integration of the orbital evolution, making template-bank searches feasible.","The mismatch study implies that eccentricity and spin must be modelled jointly: for initial eccentricity about 0.15 and spins about 0.2, neglecting the new terms already exceeds a 1% mismatch, and in high-spin, low-chirp-mass regions the mismatch reaches about 15%.","The resummed TaylorT2 extends the usable eccentricity range from about 0.45 (for the $O(e_0^8)$ series) to about 0.55 (for the resummed version), covering many residual eccentricities expected from dynamical formation.","The full $O(e_0^8)$ coefficients, supplied in the supplemental material, let users truncate at any desired eccentricity order, and the oscillatory part of the phase is shown to contribute less than about a tenth of a GW cycle for $e_0 = 0.2$, so the secular formulas dominate."],"supporting_citations":[{"why":"Supplies the original Newtonian gravitational-wave phasing for eccentric binaries, the starting point this work extends.","marker":"[13]"},{"why":"Demonstrates that even small eccentricities can bias parameter estimation, motivating the need for accurate eccentric phasing.","marker":"[40]"},{"why":"Provides the 2PN generic-spin evolution equations in the y parametrization used as a consistency input for the spinning eccentric evolution.","marker":"[50]"},{"why":"Supplies the 3PN energy and energy-flux expressions for spinning eccentric binaries that the phasing derivation takes as input.","marker":"[73]"},{"why":"Provides the non-spinning eccentric Taylor approximants and the TaylorF2Ecc model that this paper extends and compares against.","marker":"[77]"},{"why":"Details the stationary phase approximation used to convert the time-domain phasing into the TaylorF2 frequency-domain phase.","marker":"[84, 85]"},{"why":"Introduces the y frequency parameter, whose use improves post-Newtonian series convergence for higher eccentricities.","marker":"[87]"},{"why":"Contains the full $O(e_0^8)$ phasing coefficients that are too long for the main text.","marker":"[78]"}],"fun_headline_variants":["First closed-form GW phasing for eccentric spinning binaries","3PN inspiral phase now analytic for spin and eccentricity","Eccentric spinning binaries get fast analytic phase formula","New phasing formula: spin + eccentricity to 3PN, closed form","Mismatch >1% if spin-eccentric terms ignored in GW searches"],"cache_read_input_tokens":61184,"weakest_assumption_plain":"The formulas rest on the 3PN energy and flux expressions for spinning eccentric binaries supplied by Ref. [73], together with the assumption that orbit-averaged evolution, with the small oscillatory phase ignored, stays accurate up to initial eccentricities of about 0.5; if either gives way, the closed-form phases inherit the error.","fun_headline_variants_meta":{"raw":{"variants":["First closed-form GW phasing for eccentric spinning binaries","3PN inspiral phase now analytic for spin and eccentricity","Eccentric spinning binaries get fast analytic phase formula","New phasing formula: spin + eccentricity to 3PN, closed form","Mismatch >1% if spin-eccentric terms ignored in GW searches"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000274,"raw_usage":{"total_tokens":1715,"prompt_tokens":1094,"completion_tokens":621,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":710,"completion_tokens_details":{"reasoning_tokens":531}},"tokens_in":710,"tokens_out":621,"duration_ms":6069,"temperature":1.0,"reasoning_tokens":531,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T15:29:20.077185+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Decisive check: numerically integrate the full 3PN equations of motion for an aligned-spin eccentric binary with initial eccentricity 0.5 without discarding the oscillatory phase pieces, and compare the accumulated number of gravitational-wave cycles with the TaylorT2 formula and its resummed version; if the difference exceeds one cycle inside the claimed validity range, the orbit-averaged approximation underpinning the closed forms is inadequate.","supporting_citations":[],"review_version":1}