{"id":"9ecd70a4-5208-44a0-b719-ad7327911b81","arxiv_id":"2608.10965","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"A method is given for obtaining the 3.5PN radiation-reaction force and orbit correction in any Hamiltonian coordinate system by transforming from harmonic coordinates.","lead":"This paper develops a shortcut for translating radiation-reaction forces between different coordinate systems used in gravitational-wave models. It lets researchers reuse one known solution instead of redoing hard two-body calculations for each gauge.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The harmonic→EOB map in Eqs. (A3)–(A4) carries the IW gauge parameters α,β,δ_i,ϵ_5, so the 'coordinate transformation' is not a fixed geometric map; the claimed coordinate-only recipe silently absorbs the balance-gauge freedom.","rationale":"The reader's weakest assumption is the identification in Eq. (3.21) and the uniqueness of the coordinate map obtained by acceleration matching. My concern sharpens rather than replaces that: the map in Appendix A is not independent of the IW gauge parameters, so the separation between 'coordinate transformation' and 'balance-gauge choice' is not clean. The paper is honest about the need to reshuffle Schott terms, and it does demonstrate concrete gauge choices, so the work is not invalidated. However, the advertised simplification is weaker than stated: one still has to specify how Schott terms are to be transformed and how regularity of F_r is enforced. This supports the reader's CONDITIONAL verdict rather than ACCEPT. I do not see evidence of internal inconsistency that would warrant REJECT; the algebraic results may well be correct, but the central methodological claim needs either a proof of existence of the Schott reshuffling or an explicit canonical/first-principle derivation of the dissipative part of the map. The proposed numerical test for a fixed gauge would settle whether the orbit-shortcut actually reproduces the dynamics generated by the new force components.","tokens_in":26380,"tokens_out":10590,"duration_ms":108449,"concrete_test":"Adopt a fixed EOB gauge, e.g., the BD gauge (5.5) with β=2 and δ2...δ5 as in (5.10). Numerically integrate the EOB equations of motion (4.1)–(4.4) with the corresponding F_r and F_ϕ from (5.4) and E_Schott from (5.3), for an eccentric orbit, and compare at 3.5PN with the orbit obtained by applying the inverse map (A3) and the transformation formulas (6.13)–(6.14) to the harmonic solution (6.1)–(6.5). If the two orbits differ beyond the quoted O(η^7) truncation, the inverse-map shortcut is not a pure coordinate transformation and the central claim must be qualified.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that only the harmonic→X coordinate transformation is needed, without re-solving the balance equations. But the map itself is not gauge-fixed: the EOB transformation coefficients C_1...C_6 in Eq. (A4) depend linearly on the IW gauge parameters α,β,δ_1...δ_5,ϵ_5. These parameters parametrize the arbitrary split of E^* and J^* into system plus Schott terms, not the choice of coordinate chart. A genuine coordinate transformation should be fixed once the two charts are specified; here the transformation absorbs exactly the balance ambiguity the paper claims to avoid. This is not merely cosmetic: in Section V.B the direct use of Eq. (3.21) in EOB coordinates produces unwelcome 1/p_r singularities in F_r, which are removed only by manually modifying the Schott energy contribution. Thus Eq. (3.21) is not a theorem about coordinate transformations; it is a convention supplemented by a regularity requirement. Unless one proves that for every admissible gauge parameter set a regular F_r and a consistent Schott reshuffling exist, the advertised 'remarkable simplification' remains conditional on an unstated gauge-fixing procedure.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a method for obtaining the 3.5PN radiation-reaction force and Schott terms in a Hamiltonian coordinate system (ADM or EOB) from the known harmonic-coordinate Iyer-Will balance results. The authors' central claim is that, once the harmonic-coordinate solution is known, only the coordinate transformation between harmonic coordinates and the new phase-space variables is needed, without re-solving the balance equations. They provide explicit coordinate maps in Appendix A, expressions for the force components and Schott terms in ADM and EOB coordinates with the IW gauge parameters left free in the EOB case, reductions to several published EOB gauge choices, and formulas for the radiation-reaction correction to the orbit obtained by transforming the known harmonic-coordinate solution.","tokens_in":26716,"tokens_out":6170,"duration_ms":63035,"significance":"If fully substantiated, the paper would provide a convenient shortcut for translating radiation-reaction results between coordinate systems, which is directly relevant to EOB waveform models. The manuscript is valuable for its explicitness: all force components, Schott terms, and coordinate maps are given in closed form, and the EOB expressions are checked against the BD minimal gauge and the gauge choices of Refs. [20,21]. However, the advertised simplification is conditional on a number of unproven regularity and uniqueness assumptions, as detailed in the major comments; the central claim is therefore not yet established at the level claimed in the abstract.","major_comments":[{"comment":"The identification in Eq. (3.21) is the load-bearing step of the method. The paper itself notes in Eq. (3.22) that this identification is exact only at leading order, and that at 3.5PN one may reshuffle contributions via Schott terms. In the EOB case, direct application of Eq. (3.21) produces 1/p_r singularities in F_r, which are removed only by manually modifying the Schott energy (Section V.B, text preceding Eq. (5.3)). The paper does not prove that for every admissible set of IW gauge parameters a regular F_r and a consistent Schott reshuffling exist. Without such an existence argument, the claim that only the coordinate transformation is needed is not a theorem but a procedure supplemented by an unstated regularity condition. The authors should either prove the existence of the reshuffling or restate the claim as conditional on this regularity requirement.","section":"Section V.B, Eq. (3.21)"},{"comment":"The map from harmonic to EOB coordinates is not a fixed chart transformation: the coefficients C_1 through C_6 in Eq. (A4) depend linearly on the IW gauge parameters α, β, δ_i, and ε_5, which parametrize the arbitrary split of E* and J* into system plus Schott terms rather than the choice of coordinate chart. Hence the map absorbs exactly the balance ambiguity that the abstract claims to avoid. If the two charts were genuinely specified, a coordinate transformation would be fixed once the charts are given; here the result is a multi-parameter family of maps labeled by the balance gauge. The paper should either demonstrate that the EOB coordinate chart itself fixes these parameters (for example by an independent definition of EOB coordinates) or explicitly characterize the result as a family of maps parameterized by the IW gauge freedom. The derivation is also presented only as 'we find' in Appendix A, which makes it difficult to verify how the gauge dependence enters.","section":"Appendix A, Eqs. (A3)-(A4)"},{"comment":"The ADM radiation-reaction force is not cross-checked against the independent first-principle ADM result of Refs. [5,15]. Since the ADM gauge parameters are fully determined (Table I), substituting them into Eqs. (5.1)-(5.2) should reproduce the known ADM rr force, providing a direct validation of the whole apparatus. The paper currently presents these formulas without comparison. This omission is load-bearing because the ADM case is the only one where the transformation can be tested against an independent derivation; without this comparison, the reader cannot distinguish a correct derivation from an internally consistent but incorrect mapping.","section":"Section V.A"}],"minor_comments":[{"comment":"The phrase 'Differently form the ADM case' should read 'Differently from the ADM case'.","section":"Section V.B"},{"comment":"There is a typographical duplication in 'Eqs. Eqs. (40) and (43) in Ref. [20]'.","section":"Section V.B.2"},{"comment":"The notation η = 1/c as a PN expansion parameter is nonstandard and may be confused with the symmetric mass ratio ν; a brief clarification or a different symbol (e.g., ε) would improve readability.","section":"Section II"},{"comment":"The expression for δ^{rr,G3}r_h(t) contains nested parentheses that appear unbalanced in the typeset version; the authors should verify the bracket structure.","section":"Section VI, Eq. (6.4)"},{"comment":"The coordinate transformations are presented as final results with no derivation shown; given that these maps are central to the paper, a brief outline of the matching procedure and the number of equations solved at each PN order would be helpful.","section":"Appendix A"}],"recommendation":"major_revision","confidential_remarks":"The paper is technically dense and appears to be internally consistent, but the central methodological claim is stronger than what is currently demonstrated. The main issues are the conditional validity of Eq. (3.21) and the fact that the harmonic-to-EOB map is not a fixed coordinate transformation but a gauge-dependent family. Both issues are fixable by reframing or by adding the missing existence argument, and the ADM case offers a clean independent check that should be performed. I would not reject the paper on the gauge-dependence issue alone, since the authors are transparent about leaving parameters free, but the advertised 'remarkable simplification' needs to be qualified."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThis paper is a solid and useful PN calculation. The new thing is the explicit 3.5PN radiation-reaction parts of the harmonic-to-ADM and harmonic-to-EOB coordinate maps, including the full set of IW gauge parameters, and the demonstration that you can transform the rr force and the rr-corrected orbit from harmonic coordinates to these frames without re-solving the balance equations or the orbital perturbation equations. The checks against the BD minimal gauge and the Khalil et al. gauge choices are real and give confidence in the algebra.\n\nWhere it is softer: The maps in Appendix A are presented as results with no derivation shown; 'we find' is not a substitute for a reproducible calculation. More conceptually, the claim that 'only the coordinate transformation is needed' is a bit overstated. The transformation itself is not gauge-fixed—the EOB map's coefficients depend on the IW gauge parameters, so it is partly a bookkeeping of the balance ambiguity, not a purely geometric map. And Eq. (3.21) is not a theorem: the authors admit that direct use in EOB coordinates gives 1/p_r singularities in F_r, which they repair by redefining the Schott energy. That is a convention, not a consequence.\n\nThe unresolved NLO disagreement with Ref. [22] is the most concrete loose end. The authors argue it may be due to the other paper's unspecified gauge fixing. That is plausible, but they don't prove it. A referee should ask for a resolution or a clear statement of the remaining freedom.\n\nThe citation pattern looks fine; they cite the relevant harmonic, ADM, and EOB literature. No invented results.\n\nThe paper is written for EOB waveform builders and PN practitioners. If the derivations were supplied (or at least sketched) and the Ref. [22] discrepancy addressed, it would be a solid reference. As is, it is still worthwhile and deserves a serious referee, but the framing should be toned down in revision.","headline":"A useful and mostly solid 3.5PN map between harmonic and ADM/EOB radiation-reaction forces, but the 'coordinate-only' framing oversells a gauge-dependent and partially unproven recipe.","tokens_in":27221,"tokens_out":3887,"would_cite":true,"duration_ms":36826,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["04.25.Nx","04.30.-w"],"model":"deepseek-v4-flash","headline":"At 3.5PN order, radiation-reaction dynamics can be moved between coordinate systems by coordinate transformation alone, without re-solving the balance equations.","keywords":["radiation reaction","post-Newtonian expansion","3.5PN order","gauge dependence","harmonic coordinates","ADM coordinates","effective one body","Schott terms"],"falsifier":"Take a specific EOB gauge choice, integrate the EOB equations of motion with the paper's transformed $F_r$ and $F_\\phi$ for an eccentric orbit, and compare the evolved $r(t)$ and $\\phi(t)$ against the orbit obtained by inverting the coordinate map on the known harmonic solution; any difference at $O(\\eta^7)$ beyond numerical error would show that the map or the force identification is incomplete.","tokens_in":26173,"feed_emoji":"🌌","tokens_out":8151,"duration_ms":71766,"temperature":0.7,"pith_summary":"This paper claims that at the third-and-a-half post-Newtonian order, the radiation-reaction force, the Schott terms, and the radiation-reaction correction to the orbit in any Hamiltonian coordinate system can be obtained from the known harmonic-coordinate results by coordinate transformation alone, without re-solving the balance equations. The argument identifies the harmonic-coordinate mechanical energy and angular-momentum loss rates with the combination $\\dot r F_r + \\dot\\phi F_\\phi$ and the component $F_\\phi$ in the new phase-space variables, leaving all residual gauge freedom in a set of unspecified gauge parameters. Explicit maps from modified harmonic coordinates to Arnowitt-Deser-Misner (ADM) and to Effective-One-Body (EOB) phase-space coordinates are constructed through $O(\\eta^7)$, and the transformed force components and Schott terms are displayed with gauge parameters left free. If the construction is correct, it turns gauge translation into an algebraic step, which matters because EOB waveform models rely on several different gauge choices and previously needed a fresh solution of the radiation-reacted dynamics for each.","feed_headline":"One coordinate map carries radiation reaction across gauges","feed_subtitle":"Harmonic-coordinate force and orbit results translate to ADM and EOB without fresh balance-equation solves.","key_machinery":"The load-bearing mechanism is the PN-expanded coordinate and momentum map from modified harmonic coordinates to the target Hamiltonian frame, $r_h = f_r(r, \\phi, p_r, L)$, $\\phi_h = f_\\phi(r, \\phi, p_r, L)$, together with the energy-loss identification Eq. (3.21). The map's coefficients are fixed order by order by requiring that the accelerations computed from Poisson brackets in the new variables reproduce the known modified-harmonic acceleration, including the radiation-reaction part with unspecified gauge parameters; with the map in hand, the harmonic-coordinate balance results are transported to the new gauge, and its inverse converts the harmonic orbit solution to the new coordinates.","core_discovery":"On its own terms, the paper's central discovery is that the coordinate-dependent parts of the 3.5PN radiation-reaction problem are already contained in the harmonic-coordinate solution: once the losses $\\partial E_{\\rm cons}/\\partial v^i\\, a^i_{\\rm rr}$ and $\\partial J_{\\rm cons}/\\partial v^i\\, a^i_{\\rm rr}$ are identified with $\\dot r F_r + \\dot\\phi F_\\phi$ and $F_\\phi$ in a new Hamiltonian frame, Eq. (3.21), only the coordinate mapping between harmonic variables and the new phase-space variables is needed to determine the radiation-reaction force and Schott terms. The paper derives that mapping for ADM and EOB coordinates to 3.5PN order, including radiation-reaction terms in the map itself, and uses its inverse to convert the known harmonic-coordinate radiation-reacted orbit into the corresponding orbit in the new coordinates, avoiding the coupled differential equations of the variation-of-constants method.","pith_inferences":["If the maps are complete, a direct numerical check is possible: integrate the EOB equations of motion with the transformed $F_r$ and $F_\\phi$ for a given gauge choice and compare with the transformed harmonic orbit; agreement at $O(\\eta^7)$ would confirm the reshuffling into Schott terms, and disagreement would localize missing terms in the map.","The construction suggests a modular workflow for waveform models: keep the harmonic-coordinate radiation-reaction engine and the gauge parameters as a separate layer from the EOB Hamiltonian, so that changing a gauge amounts to swapping a finite set of coefficients rather than re-deriving dynamics.","An implicit consequence is that all EOB gauge choices satisfying the balance equations and the same map must predict identical gauge-invariant phasing; differences seen in published waveform models would then be attributed to higher-order or nonlocal (tail) effects beyond 3.5PN.","For eccentric orbits, the transformation of $\\dot\\phi$ contributes at 3.5PN to the phase; the explicit EOB orbit corrections provided here make it possible to test whether the next-to-leading-order quasi-circular condition used in one of the cited gauge choices is sufficient for eccentric waveforms or whether the radial component needs additional fixing."],"forward_implications":["For ADM coordinates, where all gauge parameters are already fixed, the paper gives complete explicit expressions for the 3.5PN radiation-reaction force, Schott energy, and Schott angular momentum in Hamiltonian phase-space variables.","For EOB coordinates, the same prescription yields force and Schott terms as functions of the free gauge parameters, so recent and future gauge choices in waveform models can be compared and substituted without recomputing balance equations.","The radiation-reaction correction to the orbit in the new coordinate system is obtained by inverting the coordinate map and substituting the known harmonic solution, so no separate variational-equation solve is needed for each gauge.","Because only fluxes at infinity are gauge-invariant, the residual freedom is confined to the Schott terms and the force components, matching the fact that different EOB gauge choices leave observable phasing unchanged.","The stated generalization path to 4.5PN order and to spin corrections stays within the same scheme, since only the harmonic input and the maps would need updating."],"supporting_citations":[{"why":"Supplies the balance method, the Schott-term decomposition, and the gauge parameters that all force and Schott expressions here inherit.","marker":"[13]"},{"why":"Determines the 3.5PN harmonic-coordinate gauge-parameter values that the paper uses as input for the harmonic frame.","marker":"[6, 7]"},{"why":"Establishes the leading-order harmonic radiation-reaction force and the harmonic results the paper starts from.","marker":"[2–4]"},{"why":"Fixes the ADM-coordinate gauge parameter values used in the ADM section.","marker":"[15]"},{"why":"Provides the Hamiltonian balance framework, including the identification corresponding to Eq. (3.21), that the paper generalizes to arbitrary gauge parameters.","marker":"[18]"},{"why":"Gives the harmonic-to-EOB coordinate transformation through 3PN order, extended here to 3.5PN with radiation reaction.","marker":"[22]"},{"why":"Supplies the known harmonic-coordinate radiation-reacted orbit for hyperbolic motion that the paper transforms into ADM and EOB orbits.","marker":"[24]"},{"why":"Introduces the acceleration-matching method used to fix the coefficients of the coordinate maps.","marker":"[25]"},{"why":"Provides the conservative harmonic-to-ADM transformation through 3PN order that the ADM map extends to 3.5PN.","marker":"[26]"}],"fun_headline_variants":["Single map carries 3.5PN radiation reaction to new gauges","Radiation reaction switches gauge via one transform","3.5PN radiation reaction: one map enough for new gauges","Cross-gauge radiation reaction via single coordinate map"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The construction rests on the assumption that the harmonic-coordinate mechanical energy and angular-momentum loss rates translate exactly into $\\dot r F_r + \\dot\\phi F_\\phi$ and $F_\\phi$ in the new Hamiltonian frame, with any ambiguity absorbed into Schott terms, and that the polynomial coordinate-map ansatz contains every term needed at 3.5PN order in the new gauge.","fun_headline_variants_meta":{"raw":{"variants":["Single map carries 3.5PN radiation reaction to new gauges","Radiation reaction switches gauge via one transform","3.5PN radiation reaction: one map enough for new gauges","Cross-gauge radiation reaction via single coordinate map"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000841,"raw_usage":{"total_tokens":3700,"prompt_tokens":1014,"completion_tokens":2686,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":630,"completion_tokens_details":{"reasoning_tokens":2618}},"tokens_in":630,"tokens_out":2686,"duration_ms":122243,"temperature":1.0,"reasoning_tokens":2618,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T13:17:26.322529+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a specific EOB gauge choice, integrate the EOB equations of motion with the paper's transformed $F_r$ and $F_\\phi$ for an eccentric orbit, and compare the evolved $r(t)$ and $\\phi(t)$ against the orbit obtained by inverting the coordinate map on the known harmonic solution; any difference at $O(\\eta^7)$ beyond numerical error would show that the map or the force identification is incomplete.","supporting_citations":[{"cited_title":"Post-Newtonian gravitational radiation reaction for two-body systems: Nonspinning bodies,","cited_arxiv_id":null,"evidence_quote":"Supplies the balance method, the Schott-term decomposition, and the gauge parameters that all force and Schott expressions here inherit."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the acceleration-matching method used to fix the coefficients of the coordinate maps."}],"review_version":1}