{"id":"b62d974b-cf40-404a-a51f-0ed434385e94","arxiv_id":"2608.07756","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A new 4PN gauge-invariant ISCO criterion for aligned-spin compact binaries is derived, recovering Kerr in the test-mass limit and matching numerical self-force results for retrograde and moderate spins.","lead":"This paper computes a fourth-order post-Newtonian formula for the smallest stable circular orbit of two spinning compact objects whose spins point along the orbital axis. The result extends earlier spinless and lower-order-spin calculations and is checked against exact Kerr physics and numerical self-force simulations, which makes it directly useful for modeling gravitational waves from black hole binaries.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The reduction to the 3×3 orbital block is justified only by a 1PN eigenvalue check; at 4PN, spin perturbations can couple and shift the first instability.","rationale":"We agree with the reader's weakest-assumption diagnosis. The 1PN eigenvalue ordering in Eq. (2.19) is a consistency check, not a proof: at 1PN there are no spin-orbit or spin-spin terms, so the spin sector is trivially decoupled. At 4PN the couplings in Appendices D/E are dense. The paper's two independent derivations—Hamiltonian and EoM—both invoke the same spectator-spin assumption, so their agreement does not test it. The test-mass Kerr limit and the MPD comparison validate the orbital-sector criterion in those limiting regimes, but neither checks the finite-mass, aligned-spin case where a spin/tilt instability could precede the orbital one. The GSF comparison is suggestive but does not isolate the question. Because the concern is localizable and checkable without new physics, a CONDITIONAL verdict is appropriate; we do not see grounds to reject or accept outright. We propose one explicit computation that would settle the issue: full 9×9 stability analysis at 4PN for a representative point. If the check passes, the conditional can be lifted to ACCEPT.","tokens_in":40965,"tokens_out":8468,"duration_ms":86711,"concrete_test":"From the ancillary file [40], assemble the full 9×9 Jacobian M0 of the linearized EFT EoM (2.15) for aligned-spin circular orbits at 4PN order, including the spin perturbation equations with precession vectors from Eqs. (2.4) and Appendix E. For an equal-mass representative case (ν=1/4, δ=0, e.g. χ_S=0.6, χ_Σ=0), compute the eigenfrequencies σ_p and find the smallest positive root of det(M0−iσ1)=0. Compare it with the smallest positive root of C(x) in Eq. (3.3). If the two roots coincide and all spin-sector modes have σ²>0 at that point, the decoupling assumption is validated; if a spin/tilt mode crosses zero first, the criterion in Eq. (3.3) misses the first instability.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim is that the smallest root of C(x) in Eq. (3.3) is the ISCO. This requires that, for aligned spins, the true linear instability is governed by the 3×3 orbital submatrix (2.20), with spin components as non-dynamical spectators. The only justification is the 1PN eigenvalue analysis near Eq. (2.19), where the three non-zero eigenfrequencies are ordered so that 1−6x is the first to vanish. At 1PN, however, spin couplings are absent: the leading spin-orbit acceleration enters only at 1.5PN (Eq. E4a), and the spin precession equations (2.13) couple δS⊥, δΣ⊥ to orbital variables. Appendices D/E contain numerous such couplings at 2.5PN, 3.5PN and 4PN. The Hamiltonian reduction in Eqs. (2.22)–(2.26) treats Sℓ, Σℓ as spectators but does not include δS⊥, δΣ⊥ at all. It is therefore not established that the six spin/tilt modes remain stable up to the root of C(x)=0. If a spin-coupled mode crosses zero at smaller x, Eq. (3.3) is not the stability boundary, and the central claim fails. This is a structural gap, not a coefficient error.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper derives a 4PN-order gauge-invariant stability criterion C(x) for circular orbits of compact binaries with spins aligned with the orbital angular momentum. The criterion is quoted in Eq. (3.3) for black holes and in Eq. (A1) for general compact objects, with the ISCO defined as the smallest root of C(x)=0. The derivation linearizes the EFT equations of motion and the EFT Hamiltonian about aligned-spin circular orbits, and is checked by recovering the spinless 4PN criterion, by reproducing the Kerr ISCO in the test-mass limit, and by comparing first-order-in-mass-ratio ISCO shifts with numerical GSF data and with MPD results for a spinning test particle. The paper also treats corotating black-hole binaries.","tokens_in":41176,"tokens_out":11779,"duration_ms":128078,"significance":"If correct, this is the first complete 4PN ISCO criterion for spinning, arbitrary-mass-ratio binaries and a natural extension of the authors' spinless 4PN work; it provides a gauge-invariant, parameter-free prediction that can be tested against GSF and EOB results. The manuscript has substantial strengths: the criterion is derived by two independent routes (Hamiltonian and EoM), the spinless limit and the 3.5PN harmonic-coordinate checks are reproduced, the test-mass limit exactly matches the Kerr expansion through 4PN, and the GSF/MPD comparisons are favorable in the regime where PN theory is expected to work (retrograde and moderately prograde spins). The ancillary file and appendices make the calculation transparent. The main reservation is a structural gap in the linearization argument, detailed in the major comments, which does not by itself invalidate the result but needs to be addressed.","major_comments":[{"comment":"The reduction of the linear stability problem to the 3×3 orbital matrix in Eq. (2.20) is load-bearing for the central claim, but the only justification offered is the 1PN eigenvalue statement around Eq. (2.19). Since the leading spin-orbit interaction starts at 1.5PN and spin-spin at 2PN, that 1PN check is silent about the six spin/tilt modes. The Hamiltonian reduction in Eqs. (2.22)-(2.26) simply declares S_ℓ and Σ_ℓ to be spectator variables and never includes δS_⊥, δΣ_⊥, while Appendices D and E contain couplings between transverse spin components and orbital variables at 2.5PN and higher. If any of the six spin modes has σ²=0 at an x smaller than the smallest root of Eq. (3.3), then Eq. (3.3) is not the first instability and the central claim fails. Please either prove the decoupling (for example, by showing that at the aligned-spin background the linearized spin block has purely imaginary eigenvalues and does not feed back into the orbital block) or compute and report the relevant part of the full 9×9 spectrum at 4PN.","section":"Section II.C, Eqs. (2.19)-(2.26)"},{"comment":"The favorable agreement with GSF data and with the MPD result tests the value of the first root of Eq. (3.3), but it does not test the assumption that this root is the first instability of the full linearized system. The numerical data are for the ISCO frequency itself; if a spin/tilt mode were to become unstable at smaller x, the comparison would not reveal it. The text should state this limitation explicitly or otherwise supply evidence from the full spectrum that the orbital-root criterion is the first zero of the complete system.","section":"Section III.C and Table I"}],"minor_comments":[{"comment":"The eigenvalue counting is unclear as stated: for the nine variables in Eq. (2.9), the text reports three zero modes and three nonzero modes, which accounts for only six eigenvalues; please clarify the counting or explain how the remaining modes are treated.","section":"Section II.B, after Eq. (2.19)"},{"comment":"The statement that the negative x^4 coefficient in C_corot - C_non-spin is \"sub-dominant\" would be more convincing if quantified at the relevant ISCO values (x ≈ 0.2) for the mass-ratio range 0 ≤ ν ≤ 1/4, rather than inferred purely from PN-order counting.","section":"Section III.E, Eqs. (3.41)-(3.42)"},{"comment":"The notation κ_+, κ_-, λ_+, λ_-, ι_+, ι_- is defined only in the appendix; since the main text's Eq. (3.3) already uses the BH values κ=λ=ι=1, a one-line glossary near Eq. (3.3) would help readers who do not go through Appendix A.","section":"Appendix A"}],"recommendation":"major_revision","confidential_remarks":"This is a technically impressive paper and the main calculation is very likely correct, but the decoupling of the six spin/tilt modes from the 3×3 orbital block is an assumption rather than a demonstrated result in the present text. The gap is checkable within the manuscript's own framework, so I do not see it as grounds for rejection; however, it should be resolved before publication. I would ask the authors to supply the missing spectral analysis or an explicit and clearly delimited assumption, and to adjust the wording of the stability claim accordingly. The GSF and MPD comparisons are honestly reported and the breakdown for near-extremal prograde spin is appropriately acknowledged."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nHere is my read. The headline: this is the first 4PN gauge-invariant ISCO criterion for aligned-spin, arbitrary-mass binaries, and it is mostly careful, sound PN work. The one real soft spot is the unverified decoupling of the spin perturbations from the orbital stability analysis.\n\nWhat is new and good: the authors extend Favata's 2.5PN spin criterion and their own spinless 4PN result to include NNLO spin-orbit, NNLO spin-spin, and LO cubic and quartic spin terms. They derive C(x) by two independent routes, Hamiltonian and equations of motion, and provide the spinless 4PN EFT Hamiltonian in Appendix C, which is a useful byproduct. The test-mass limit reproduces the Kerr ISCO to 4PN. The first-order-in-nu ISCO shifts match conservative GSF data to within 10% for retrograde and moderately prograde spins, and the MPD comparison for the test-particle spin is excellent for chi less than about 0.5. The ancillary file gives all formulas in machine-readable form. These checks are real and meaningful.\n\nThe soft spot is the reduction from the full linear system to the 3x3 orbital block. The authors justify it with a 1PN eigenvalue calculation at Eq. (2.19), but at 1PN spin couplings are absent; spin-orbit enters only at 1.5PN. From then on, the perpendicular spin components delta-S_perp and delta-Sigma_perp couple to the orbital sector through the precession equations. In the Hamiltonian derivation they treat S_l and Sigma_l as non-dynamical spectators and drop the other spin components entirely. Both derivations share this assumption, so the agreement between them does not resolve it. The external tests—Kerr limit, GSF, MPD—mostly probe the test-mass or test-particle limit, where the small body's spin sector is trivial; they do not validate the comparable-mass spin-mode decoupling. So the paper does not establish that the smallest root of C(x)=0 is the first instability. This is a structural gap, not a coefficient error.\n\nA secondary point: the 4PN spin sector is taken from Levi-Steinhoff without an independent gauge check at 4PN; the harmonic-coordinate cross-check stops at 3.5PN because the 4PN harmonic spin terms are not in the literature. That is a data limitation, fairly stated.\n\nBottom line: the paper deserves a serious referee. I would send it to peer review and ask the authors to close the spin-decoupling gap—ideally by performing the full eigenfrequency analysis of the 9x9 system or by providing a convincing argument why the spin modes cannot go unstable earlier. If that is done, this will be a standard reference.","headline":"First 4PN ISCO criterion for spinning comparable-mass binaries, carefully derived and well checked, but the spin-mode decoupling assumption needs to be justified before I'd take C(x)=0 as the stability boundary.","tokens_in":41790,"tokens_out":5021,"would_cite":true,"duration_ms":48612,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper derives a gauge-invariant stability function $C(x)$ that locates the innermost stable circular orbit of arbitrary-mass, aligned-spin compact binaries through fourth post-Newtonian order, and shows it recovers the Kerr ISCO in…","keywords":["innermost stable circular orbit","post-Newtonian approximation","spinning compact binaries","aligned spins","Kerr ISCO","gravitational self-force","spin-orbit coupling","spin-spin coupling"],"falsifier":"Compute the full $9\\times 9$ linearised matrix at 4PN for aligned spins using the paper's ancillary expressions and test whether any mode with $\\delta S_n$, $\\delta S_\\lambda$, $\\delta\\Sigma_n$, or $\\delta\\Sigma_\\lambda$ acquires $\\sigma^2<0$ at an $x$ smaller than the smallest root of $C(x)=0$; if such a mode exists, the ISCO is not at that root. A direct numerical alternative is to measure the conservative ISCO frequency of an equal-mass aligned-spin binary, for example with $\\chi_1=\\chi_2=0.6$, and compare it with the root of Eq. (3.3).","tokens_in":40711,"feed_emoji":"🌀","tokens_out":7557,"duration_ms":69626,"temperature":0.7,"pith_summary":"The paper aims to locate the innermost stable circular orbit (ISCO) of a compact binary when the two bodies have arbitrary mass and spins aligned with the orbital angular momentum, using the post-Newtonian approximation through fourth order (4PN). The authors derive a single gauge-invariant function $C(x)$ of the invariant orbital-frequency parameter $x=(Gm\\Omega/c^3)^{2/3}$ such that circular orbits are stable when $C(x)>0$ and the ISCO sits at the smallest root $C(x)=0$. In the test-mass limit the function reproduces the exact Kerr ISCO to 4PN order, and at first order in the mass ratio it reproduces numerical gravitational-self-force shifts of the ISCO frequency to within about 10 percent for retrograde and moderately prograde spins. The same criterion also gives the ISCO shift produced by a small body's own spin, matching the exact Mathisson-Papapetrou-Dixon prediction for a spinning particle in Kerr over the same spin range. If correct, this gives a parameter-free, analytic ISCO condition usable for comparable-mass binaries and for neutron stars once internal-structure coefficients are supplied.","feed_headline":"4PN ISCO criterion for spinning binaries matches self-force","feed_subtitle":"Gauge-invariant stability function reproduces the Kerr ISCO and agrees with numerical shifts for retrograde spins.","key_machinery":"The central object is the invariant stability function $C(x)$, a renormalised version of the determinant condition coming from the $3\\times 3$ orbital block of the linearised system, written in terms of $x=(Gm\\Omega/c^3)^{2/3}$ and the dimensionless aligned-spin variables $\\chi_S$ and $\\chi_\\Sigma$. It is computed by perturbing the 4PN effective-field-theory Hamiltonian and, independently, the EFT equations of motion around a circular orbit; the Hamiltonian route reduces to the condition $C=\\pi_0\\sigma_0-\\rho_0\\theta_0>0$, built from second functional derivatives of the Hamiltonian, while the equations-of-motion route uses the matrix $M_0$ with coefficients $\\alpha_0$, $\\beta_0$, $gamma_0$. The spins enter as non-dynamical spectators, an assumption explicitly checked at 1PN and then carried to 4PN. For black holes all spin-induced deformability coefficients are unity, producing the explicit polynomial (3.3); for other compact objects the general expression (A1) retains the parameters $\\kappa_a$, $\\lambda_a$, $\\iota_a$.","core_discovery":"In the paper's own terms, the central discovery is that the 4PN gauge-invariant ISCO criterion for aligned-spin compact binaries is the explicit function $C(x)$ of Eq. (3.3) for black holes and Eq. (A1) for arbitrary compact objects, with the ISCO at the smallest root of $C(x)=0$. The paper derives this function twice, once from the 4PN EFT Hamiltonian and once from the EFT equations of motion, and verifies the 3.5PN part in harmonic coordinates. In the test-mass limit it reduces to the 4PN expansion of the exact Kerr stability criterion, and the first-order-in-$\\nu$ ISCO-frequency shift agrees with numerical conservative gravitational-self-force calculations within about 10 percent for retrograde spins and moderately prograde spins ($\\chi_2\\lesssim 0.2$). For a spinning test particle the criterion reproduces the MPD prediction of the ISCO shift to within a few percent for $\\chi_2\\lesssim 0.3$. The paper also claims that corotating black-hole binaries are more stable than irrotational ones with the same masses and orbital frequency.","pith_inferences":["A concrete testable extension is to evaluate the full $9\\times 9$ linearised matrix at 4PN for aligned spins, using the paper's ancillary expressions, and check whether any transverse spin mode develops an instability before the root of $C(x)=0$; if one does, the 3-by-3 reduction would fail and the ISCO would shift.","The close agreement with self-force data for retrograde spins suggests the PN criterion could serve as a cheap surrogate for conservative self-force ISCO frequencies in waveform modeling of extreme-mass-ratio inspirals with retrograde spins.","For equal-mass neutron-star binaries the deformability parameters are not unity, so the paper's general formula gives a concrete prediction that future numerical-relativity or tidal-disruption calculations could test.","Because the criterion is equivalent to the inverse square of the periastron-precession factor for circular orbits, future higher-order precession calculations could be converted into ISCO estimates without repeating the full perturbation analysis."],"forward_implications":["The smallest root of $C(x)=0$ gives the 4PN ISCO frequency for any mass ratio with aligned spins, so the ISCO can be found without interpolation or fitting.","In the extreme-mass-ratio limit the criterion recovers the Kerr ISCO to 4PN order, extending the known Schwarzschild recovery to the spinning case.","The first-order-in-$\\nu$ ISCO shift matches conservative self-force numerics within about 10 percent for retrograde spins and moderate prograde spins, and the failure near maximal prograde spin is explained by the ISCO approaching the horizon.","For neutron stars or other compact objects, inserting the appropriate deformability coefficients $\\kappa_a$, $\\lambda_a$, $\\iota_a$ into Eq. (A1) yields the ISCO criterion without a new derivation.","Corotating black-hole binaries are predicted to be more stable than irrotational ones at fixed orbital frequency and mass."],"supporting_citations":[{"why":"Introduces the linear-perturbation analysis of circular orbits whose eigenvalues give the stability criterion.","marker":"[1]"},{"why":"Provides the spinless 4PN criterion that this paper extends to spins, including the non-local tail term.","marker":"[9]"},{"why":"Supplies the earlier spin-dependent criterion up to 2.5PN and the Kerr and test-spin comparisons this work generalizes.","marker":"[11]"},{"why":"Supplies the effective-field-theory Hamiltonian with all spin terms through 4PN order from which the criterion is derived.","marker":"[25–30]"},{"why":"Provides the numerical first-order gravitational-self-force ISCO shifts used for comparison.","marker":"[36]"},{"why":"Gives the exact Kerr stability criterion and ISCO used in the test-mass limit.","marker":"[44]"},{"why":"Gives the Mathisson-Papapetrou-Dixon equations for a spinning particle in Kerr used to derive the exact test-particle-spin ISCO shift.","marker":"[46]"}],"fun_headline_variants":["4PN spin ISCO matches self-force for retrograde spins","Exact 4PN ISCO for aligned spins: Kerr limit restored","Spinning binary ISCO at 4PN: self-force verified","ISCO for spinning binaries: 4PN beats test-mass limit","Corotating binaries more stable: new 4PN ISCO criterion"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the spin components can be treated as non-dynamical spectators when the circular orbit with aligned spins is perturbed, so that stability is decided by the 3-by-3 orbital block alone; this decoupling is explicitly checked only at 1PN and is assumed to persist at 4PN.","fun_headline_variants_meta":{"raw":{"variants":["4PN spin ISCO matches self-force for retrograde spins","Exact 4PN ISCO for aligned spins: Kerr limit restored","Spinning binary ISCO at 4PN: self-force verified","ISCO for spinning binaries: 4PN beats test-mass limit","Corotating binaries more stable: new 4PN ISCO criterion"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00044,"raw_usage":{"total_tokens":2278,"prompt_tokens":1040,"completion_tokens":1238,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":656,"completion_tokens_details":{"reasoning_tokens":1160}},"tokens_in":656,"tokens_out":1238,"duration_ms":10058,"temperature":1.0,"reasoning_tokens":1160,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T00:20:33.113743+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the full $9\\times 9$ linearised matrix at 4PN for aligned spins using the paper's ancillary expressions and test whether any mode with $\\delta S_n$, $\\delta S_\\lambda$, $\\delta\\Sigma_n$, or $\\delta\\Sigma_\\lambda$ acquires $\\sigma^2<0$ at an $x$ smaller than the smallest root of $C(x)=0$; if such a mode exists, the ISCO is not at that root. A direct numerical alternative is to measure the conservative ISCO frequency of an equal-mass aligned-spin binary, for example with $\\chi_1=\\chi_2=0.6$, and compare it with the root of Eq. (3.3).","supporting_citations":[{"cited_title":"Kidder, C","cited_arxiv_id":null,"evidence_quote":"Introduces the linear-perturbation analysis of circular orbits whose eigenvalues give the stability criterion."},{"cited_title":"Bardeen, W","cited_arxiv_id":null,"evidence_quote":"Gives the exact Kerr stability criterion and ISCO used in the test-mass limit."},{"cited_title":"Saijo, K.-i","cited_arxiv_id":null,"evidence_quote":"Gives the Mathisson-Papapetrou-Dixon equations for a spinning particle in Kerr used to derive the exact test-particle-spin ISCO shift."}],"review_version":1}