{"id":"58b4901b-5295-4447-af59-33e94ace5ce9","arxiv_id":"2412.16077","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"For circular equatorial orbits, EOB and MPD dynamics and gravitational-wave fluxes agree for Schwarzschild primaries, while for Kerr the difference grows with spin and is largest for high positive primary spin.","lead":"This paper compares two ways of modeling a small spinning object orbiting a rotating black hole: the effective-one-body approximation used in gravitational-wave templates and the Mathisson-Papapetrou-Dixon test-body equations. It finds they agree for a non-spinning black hole and differ most for fast-spinning black holes, which matters for future gravitational-wave detectors.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Flux comparison is inconsistent: a non-linear Teukolsky solver is fed with first-order-in-σ EOB/MPD orbits, so the reported flux differences for σ up to 0.9 are not differences of the full models; the Schwarzschild equality is only equality of linearized dynamics.","rationale":"I read the paper as comparing linearized-in-σ conservative EOB and MPD dynamics, then using those linearized energy/angular-momentum data to drive a Teukolsky flux calculation that is not itself linearized in σ. For the reported flux differences to describe the EOB and MPD models at σ up to 0.9, O(σ²) orbital corrections must be negligible at these spins. They are not: the paper's own Table II shows the full non-linearized energies deviate from the linearized value at the ~1% level already at σ=0.5 and r=4, and the authors note that quadratic-in-spin contributions to fluxes are incomplete for large σ. Consequently the central Schwarzschild 'no difference' result reflects the equality of the truncated input orbits, not equality of EOB and MPD dynamics; the Kerr flux differences likewise mix the linearized orbit difference with uncontrolled higher-order source terms. This is exactly the reader's weakest assumption, so I agree with it. The analytic 3PN angular-momentum comparison and the linear-order Schwarzschild identity are internally consistent and provide real support for the limited linear-order conclusions, but they do not validate the finite-σ flux plots. The missing u(x) formula and placeholder supplementary URL are secondary reproducibility issues. A conditional verdict remains appropriate; the concern should be resolved by redoing the flux comparison with non-linearized dynamics at a few representative spins or by restricting the presented claims to small σ.","tokens_in":20172,"tokens_out":6497,"duration_ms":60924,"concrete_test":"Compute the Schwarzschild (a=0) fluxes feeding the Teukolsky solver with the full, non-linearized EOB circular-orbit solution (solve ∂H_eff/∂u=0 exactly) and the full Tulczyjew-Dixon MPD circular-orbit solution, for σ = 0.1, 0.5, 0.9. If the EOB-MPD flux difference is nonzero at the percent level expected from Table II (where full and linearized energies differ by ~1% at r=4, σ=0.5), the reported equality is an artifact of σ-linearization; if the difference remains zero to numerical precision, the concern is retired.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central flux comparison requires first-order-in-σ orbital data to be sufficient for Teukolsky fluxes at σ up to 0.9. Section VI states the code is not linearized in the secondary spin, while the orbital inputs are linearized: EOB u(x) from Sec. IIB via Eqs. (25)-(26) and MPD u(x) in Eq. (48). Because the flux is quadratic in the perturbative source, feeding an O(σ)-accurate orbit into an unlinearized solver yields a mixture: some O(σ²) source contributions are included, but those depending on O(σ²) orbital corrections are omitted. Thus the EOB-MPD flux difference is not a clean comparison of the two dynamics. For a Schwarzschild primary the linearized EOB and MPD inputs coincide exactly, so the Teukolsky output is identical by construction; Appendix B and Table II show that the full non-linearized dynamics do not coincide (at r=4, σ=0.5, EOB energy 0.950686 vs MPD 0.946587 vs linearized 0.937500). The claimed equality is therefore an artifact of the truncation. The authors acknowledge the quadratic-in-spin incompleteness for large σ (Sec. VI, horizon flux for {a,σ}={0.9,0.9}), but the same caveat applies to the headline infinity-flux differences for all σ=0.9 results. The 3PN angular-momentum argument and linear-order Schwarzschild identity appear sound; what is unsupported is the interpretation of the finite-σ flux plots as EOB-vs-MPD differences.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper compares, at linear order in the secondary spin σ, the circular equatorial orbit dynamics of the effective-one-body (EOB) Hamiltonian of Ref. [15] and the Mathisson-Papapetrou-Dixon (MPD) formalism under the Tulczyjew-Dixon spin supplementary condition around a Kerr black hole. It derives radius-frequency relations u_EOB(x) and u_MPD(x), compares energies and angular momenta, discusses last stable orbits, and then feeds the linearized orbital data into a frequency-domain Teukolsky solver to obtain gravitational-wave energy fluxes at infinity and into the horizon. The main reported results are: the angular-momentum difference between EOB and MPD starts at 3PN, consistent with the EOB spin-orbit sector being complete through 2.5PN; for Schwarzschild primaries the EOB and MPD fluxes coincide; and for Kerr primaries EOB fluxes are larger than MPD for negative σ and smaller for positive σ.","tokens_in":20529,"tokens_out":2972,"duration_ms":29129,"significance":"If fully established, the comparison would be a useful benchmark for improving the radiation-reaction sector of EOB models for spinning test particles in Kerr, which is a timely step for EMRI waveform modeling. The paper contains several clean and verifiable analytical results: the PN expansions in Eqs. (49)-(50) explicitly show the expected 2.5PN limitation of the EOB spin-orbit sector, the NNLO centrifugal-radius modification is shown to matter, and the linear-order Schwarzschild equality is a transparent consistency check. The authors are also candid about the need for the anti-DJS spin gauge and about the absence of an EOB LSO for large positive primary spins. However, the significance of the finite-σ flux differences is weakened by a truncation inconsistency: the Teukolsky solver is not linearized in σ, while the orbital inputs are, so the plotted flux differences at σ = ±0.9 are not established as differences of the full EOB and MPD models.","major_comments":[{"comment":"The central flux comparison is affected by an inconsistent perturbative truncation. Section VI explicitly states that the Teukolsky solver is not linearized in the secondary spin, but the orbital data fed into it are linearized in σ via the u(x) and energy/angular-momentum expressions of Secs. IIB, III C and Eq. (48). Since the gravitational-wave flux is quadratic in the source, a solver that is not linearized will include some O(σ²) source contributions while omitting the O(σ²) orbital corrections. Consequently, the flux differences shown in Figs. 5 and 6 for σ up to 0.9 are not clean differences between the full EOB and full MPD dynamics. This is not merely a cosmetic issue: for the Schwarzschild case, the equality in Fig. 7 is guaranteed by construction once both inputs are linearized, because u_EOB(x) and u_MPD(x) coincide at linear order. Appendix B, Table II, shows that the full, non-linearized EOB and MPD energies do not coincide (e.g., at r=4, σ=0.5: 0.950686 vs 0.946587), while both linearized values equal 0.937500. The authors acknowledge incomplete quadratic-in-spin flux contributions for the horizon flux at {a,σ}={0.9,0.9}, but the same caveat applies to the headline asymptotic-flux differences for all large-σ cases. The paper should either restrict the flux conclusions to spin magnitudes where O(σ²) contributions are demonstrably negligible, or supply a consistent linearized Teukolsky computation, or include nonlinear-in-σ orbital data.","section":"Sec. VI and Figs. 5-7"},{"comment":"The central EOB radius-frequency relation u_EOB(x) is not shown in the text: the reader is told that 'for practical reasons we do not present the final formula here' and is referred to the supplementary material, but reference [27] contains only 'URL-will-be-inserted' and the supplementary notebook is not available in the arXiv version. Since u_EOB(x) is the input to every subsequent EOB energy, angular-momentum and flux evaluation, the paper is currently not independently reproducible at the point where its EOB results begin. The formula should be included in an appendix, or at least a machine-readable expression should be made available with the submission.","section":"Sec. IIB and Ref. [27]"},{"comment":"The claimed linear-in-σ horizon-flux behavior rests on a numerical fit: the statement 'we have found numerically that E_H^1 = ...' appears without a derivation or an error estimate. This is used to infer that the linear-in-σ part starts 3/2 PN orders higher than the nonspinning term, and it feeds the qualitative discussion in Appendix A. As a fitted auxiliary relation it does not by itself invalidate the main conclusions, but the paper should present the fitting procedure, the range and number of data points, and the residual quality, so that the reader can judge whether Eq. (60) is a PN result or an empirical interpolation.","section":"Eq. (60) and Appendix A"}],"minor_comments":[{"comment":"There are several typographical and wording issues: 'liner order' in the Introduction, 'normaliszd' in the caption of Fig. 8, 'Schwarzchild' in Sec. VII, 'if he had not chosen' in Sec. III (should be 'if we had not chosen'), and inconsistent capitalization of sigma (σ vs lowercase 'sigma') in several places.","section":"Throughout"},{"comment":"In the text describing Fig. 5 the sentence 'As for â = −0.9, we actually consider the fluxes only up to the largest LSO value, which is ∼ 0.135' appears twice in slightly different forms; this should be condensed to avoid duplication.","section":"Sec. VI and Fig. 5"},{"comment":"The conclusion that 'EOB spin-orbit interaction is stronger than the MPD one' is presented as an explanation of the flux ordering, but the paper also notes that spin-spin contributions play a role and that the ordering changes with parameters. The wording should be softened to reflect the fact that the flux ordering is not uniquely determined by the spin-orbit term alone.","section":"Sec. VII"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a reasonable step toward benchmarking EOB radiation reaction against MPD/Tulczyjew-Dixon results for spinning test particles in Kerr, and the analytical parts (radius-frequency relations, PN expansions, LSO discussion) are worth publishing after revision. The main issue is the mismatch between linearized orbital inputs and a nonlinear Teukolsky solver; the authors should either restrict the flux comparison to small σ or provide a consistent truncation analysis. I would also urge the editor to require that the supplementary notebook be made available at submission, since the central u_EOB(x) formula is only in that file."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe genuinely new thing here is the closed-form u_MPD(x) relation for a spinning particle on Kerr, Eq. (48), plus the first EOB-vs-MPD flux comparison on Kerr using a frequency-domain Teukolsky solver. The derivations are clear, the 3PN angular-momentum disagreement matches the known 2.5PN truncation of the EOB spin-orbit sector, and the LSO analysis—especially the absence of an EOB LSO for large positive a and negative sigma—is interesting and worth pursuing.\n\nThe paper is honest about its main limitation: the orbital data fed into the Teukolsky code are linearized in the secondary spin sigma, while the solver is not. That means the O(sigma^2) flux contributions are incomplete. For sigma up to 0.9 the difference between EOB and MPD fluxes is not a clean comparison of the full dynamics. The Schwarzschild \"equality\" is by construction: at linear order the two orbits coincide, so the Teukolsky output is identical. Appendix B actually shows the non-linearized energies differ. The authors do flag the quadratic-in-spin caveat for the {a,sigma}={0.9,0.9} horizon flux, so the issue is not hidden, but the abstract's unqualified \"no difference\" is misleading, and the finite-sigma flux curves should not be read as exact benchmarks.\n\nOther soft spots: the key EOB u(x) formula is only in a supplementary notebook whose URL is a placeholder, and there are no numerical error estimates on the flux differences. For small sigma (say, |sigma| <= 0.2) the comparison is probably reliable in its qualitative and even semi-quantitative behavior, and the analytic PN results stand independent of the flux computation.\n\nBottom line: this is a useful paper for people building EOB radiation reaction for EMRIs, and the new u_MPD(x) relation is worth having. But it needs a revision: either restrict the flux comparison to small sigma or feed non-linearized orbits, supply the missing formula, and add error estimates. I'd send it to a serious referee; it's not ready as is.\n\nI would bring it to the reading group and will probably cite it for Eq. (48).","headline":"The new u_MPD(x) relation for Kerr is solid and the paper is a useful benchmark, but the headline flux comparison is contaminated by feeding linearized-in-sigma orbits into an unlinearized Teukolsky solver, so the finite-sigma differences should not be taken as exact EOB-vs-MPD benchmarks.","tokens_in":21090,"tokens_out":2928,"would_cite":true,"duration_ms":26341,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["83C57","83C10","83C25"],"pacs":["04.30.-w","04.70.-s"],"model":"deepseek-v4-flash","headline":"For a spinning test particle on circular equatorial orbits, the effective-one-body and Mathisson-Papapetrou-Dixon formalisms give identical gravitational-wave fluxes for a Schwarzschild primary, and differ for Kerr primaries, with the…","keywords":["effective-one-body","Mathisson-Papapetrou-Dixon","Tulczyjew-Dixon spin supplementary condition","Kerr black hole","circular equatorial orbits","Teukolsky equation","gravitational-wave fluxes","spin-orbit coupling"],"falsifier":"Compute the same Teukolsky fluxes using fully nonlinear-in-$\\sigma$ EOB and MPD orbital data at $\\sigma=0.9$ for a Schwarzschild primary; if the EOB and MPD fluxes no longer coincide, the reported Schwarzschild equality is an artifact of the linear-in-$\\sigma$ truncation.","tokens_in":19915,"feed_emoji":"🕳️","tokens_out":8056,"duration_ms":69360,"temperature":0.7,"pith_summary":"The paper asks whether the effective-one-body (EOB) approach and the Mathisson-Papapetrou-Dixon (MPD) description of a spinning test particle predict the same gravitational-wave emission for circular equatorial orbits around a Kerr black hole. It derives the radius–frequency relation, energy, and angular momentum at linear order in the secondary spin $\\sigma$, then feeds these to a frequency-domain Teukolsky solver. The central finding is that for a Schwarzschild primary the EOB and MPD fluxes agree exactly, while for a Kerr primary EOB fluxes are larger than MPD for negative $\\sigma$ and smaller for positive $\\sigma$, with the gap growing at high primary spin. This matters for building EOB radiation-reaction models for extreme-mass-ratio inspirals with a spinning secondary.","feed_headline":"Two black-hole orbital descriptions agree on Schwarzschild, not Kerr","feed_subtitle":"Around Kerr black holes the two descriptions diverge at 3PN order, with largest flux gaps at high positive spin.","key_machinery":"The comparison runs through the linear-in-$\\sigma$ radius–frequency relation $u(x)$, where $x = \\Omega^{2/3}$ is the orbital-frequency parameter. For EOB this relation follows from imposing circular conditions on the effective Hamiltonian, whose spin-orbit sector uses the gyrogravitomagnetic functions $G_S$ and $G_{S*}$ with the complete zeroth-order self-force expression and an NNLO centrifugal radius; for MPD it follows from expanding the orbital-frequency formula of the MPD equations under the Tulczyjew-Dixon condition. These $u(x)$ relations convert both formalisms to the same gauge-invariant frequency parameter, allowing direct comparison of energy, angular momentum, and last-stable-orbit location, and providing the input for the frequency-domain Teukolsky equation solver that computes the fluxes.","core_discovery":"At linear order in the secondary spin $\\sigma$, and using the Tulczyjew-Dixon spin supplementary condition, the paper establishes that the test-mass EOB Hamiltonian and the MPD equations produce the same circular-orbit dynamics when the central black hole is Schwarzschild: the $u(x)$ relations, energies, and angular momenta coincide, so a Teukolsky-based flux computation returns identical infinity and horizon fluxes for every $x$ and $\\sigma$ considered. For a Kerr background, the formalisms differ: the EOB energy and angular momentum are larger/smaller than MPD for negative/positive $\\sigma$, the EOB asymptotic and horizon fluxes are correspondingly larger/smaller, and the angular-momentum difference begins at 3PN order, consistent with the EOB spin-orbit sector being complete only through 2.5PN. The paper also shows that an earlier reported EOB/MPD flux difference on Schwarzschild came from not linearizing the dynamics in $\\sigma$.","pith_inferences":["If the linear-order agreement on Schwarzschild persists at higher order in $\\sigma$, it would mean the EOB test-mass spin-orbit sector and the Tulczyjew-Dixon MPD description are aligned in the non-spinning-primary limit, and the Kerr discrepancy is primarily a spin-orbit truncation effect rather than a fundamental formalism mismatch.","A natural testable extension is to recompute both the dynamics and the fluxes without linearizing in $\\sigma$; the paper's own Appendix B suggests that at $\\sigma=0.5$ the nonlinear EOB and MPD energies already disagree on Schwarzschild, so the equality likely degrades as spin grows.","The horizon-flux results hint that spin-spin effects dominate horizon absorption near the last stable orbit; feeding the Teukolsky solver with fully nonlinear orbital data for $\\{\\hat a,\\sigma\\} = \\{0.9,0.9\\}$ would separate spin-orbit from spin-spin contributions.","Extending the same comparison to eccentric equatorial orbits or to second order in the secondary spin would test whether the EOB/MPD flux difference observed here is a robust feature of the spin sector or a peculiarity of circular orbits."],"forward_implications":["For a nonspinning primary, EOB and MPD (Tulczyjew-Dixon) give the same asymptotic and horizon energy fluxes at linear order in the secondary spin, providing a benchmark for EOB radiation reaction in the Schwarzschild limit.","On Kerr, the flux mismatch is systematic: EOB fluxes exceed MPD for $\\sigma<0$ and fall below for $\\sigma>0$, so the choice of formalism visibly changes predicted gravitational-wave emission from spinning extreme-mass-ratio inspirals.","The 3PN start of the angular-momentum difference pinpoints the EOB spin-orbit sector's known truncation at 2.5PN as the source of the flux gap for Kerr primaries.","The linearized EOB last-stable-orbit behavior for large positive $\\hat a$ and negative $\\sigma$ hides a genuine absence of a last stable orbit in the full non-linearized EOB dynamics, a limitation for templates in that parameter region.","The results extend the previous Schwarzschild-only comparison to Kerr and, unlike that earlier work, explain the Schwarzschild equality as a consequence of consistent linearization in $\\sigma$."],"supporting_citations":[{"why":"Supplies the EOB Hamiltonian with enhanced spin-orbit sector whose test-mass limit is used throughout.","marker":"[15]"},{"why":"Provides the MPD circular-orbit frequency expression and asymptotic-flux results that the paper expands to linear order in the secondary spin.","marker":"[4]"},{"why":"The precursor Schwarzschild comparison; Appendix B contrasts its non-linearized dynamics with the linearized approach that yields equal fluxes.","marker":"[5]"},{"why":"The frequency-domain Teukolsky equation solver used to evaluate the asymptotic and horizon fluxes.","marker":"[36]"},{"why":"Supplies the analytical post-Newtonian horizon-flux formula used in Appendix A to benchmark the numerical horizon fluxes.","marker":"[40]"},{"why":"Gives the complete zeroth-order self-force spin-orbit coupling used for the EOB gyrogravitomagnetic function $G_{S*}$.","marker":"[24]"},{"why":"Provides the NNLO centrifugal-radius correction whose inclusion is needed to reproduce the 3PN angular-momentum agreement.","marker":"[23]"},{"why":"Provides the MPD last-stable-orbit expressions used for the $x_{\\rm LSO}$ comparison.","marker":"[38]"}],"fun_headline_variants":["Black hole orbits: two theories match on Schwarzschild, not Kerr","Spinning test particles: EOB and MPD agree on non-rotating holes","Orbit models agree for Schwarzschild, diverge for Kerr at high spin","Kerr black holes expose mismatch in orbital descriptions","EOB and MPD fluxes match on Schwarzschild, differ on Kerr"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The comparison is made only to first order in the secondary particle's spin $\\sigma$, yet fluxes are evaluated for $\\sigma$ up to 0.9, where quadratic-in-spin effects are sizable and incomplete.","fun_headline_variants_meta":{"raw":{"variants":["Black hole orbits: two theories match on Schwarzschild, not Kerr","Spinning test particles: EOB and MPD agree on non-rotating holes","Orbit models agree for Schwarzschild, diverge for Kerr at high spin","Kerr black holes expose mismatch in orbital descriptions","EOB and MPD fluxes match on Schwarzschild, differ on Kerr"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000211,"raw_usage":{"total_tokens":1402,"prompt_tokens":922,"completion_tokens":480,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":538,"completion_tokens_details":{"reasoning_tokens":386}},"tokens_in":538,"tokens_out":480,"duration_ms":4490,"temperature":1.0,"reasoning_tokens":386,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T10:48:30.603643+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the same Teukolsky fluxes using fully nonlinear-in-$\\sigma$ EOB and MPD orbital data at $\\sigma=0.9$ for a Schwarzschild primary; if the EOB and MPD fluxes no longer coincide, the reported Schwarzschild equality is an artifact of the linear-in-$\\sigma$ truncation.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the MPD circular-orbit frequency expression and asymptotic-flux results that the paper expands to linear order in the secondary spin."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The frequency-domain Teukolsky equation solver used to evaluate the asymptotic and horizon fluxes."}],"review_version":1}