{"id":"78e94643-c0c1-4a77-baaa-58dfe60cb360","arxiv_id":"2505.09052","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"A revised boundary-to-bound dictionary, using p_inf^2 -> p_inf^2 e^(-i pi), makes fourth-order post-Minkowskian bound-orbit observables real and matches numerical relativity binding energies for mass ratios 1 and 10.","lead":"This paper proposes a fix to the standard 'dictionary' that converts gravitational scattering calculations into predictions for bound orbits of black hole binaries. The fix uses a specific analytic continuation, and the authors compare the resulting binding energies with numerical relativity simulations.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Branch rule underdetermined: p_inf^2 -> p_inf^2 e^{-i pi} does not fix whether logs/powers use arg +pi or -pi, the two branches give different complex results, and the paper never tests the alternative nor the 3PM consistency limit.","rationale":"The central claim is that the 4PM dictionary must be replaced by the analytic continuation p_inf^2 -> p_inf^2 e^{-i pi}. For this to be correct, the continuation must be unique (or otherwise physically selected), and it must reduce to the known 3PM dictionary where the original Kalin-Porto dictionary is uncontroversial. The paper shows only that its chosen prescription yields real 2PM/4PM precession angles and a sparse EOB/NR comparison; it does not compare the two possible branch limits or the 3PM truncation. Since e^{-i pi} = -1, the substitution is literally a map to the same complex value and is therefore incomplete as a statement about non-integer powers and logarithms. This is exactly the weakest assumption identified by the reader, and it is the most load-bearing because if the alternative branch (e^{+i pi}) yields different observables, the whole dictionary revision is not uniquely defined; if they yield the same observables, the paper's phrasing is merely ambiguous and a simple clarification would suffice. The SXS comparison is not decisive: it is sparse, the EOB model uses the same 4PM inputs, and Ref. [56] (needed for Eq. (6)) is listed as 'To appear 111, 111 (2024)' with no arXiv ID, so it cannot be checked. The concern does not overturn the reader's conditional verdict; it reinforces the need for the requested derivation and validation.","tokens_in":10304,"tokens_out":12274,"duration_ms":118800,"concrete_test":"Implement the two-sided limit explicitly: for a representative bound orbit (e.g., mass ratio q = 1, a chosen gamma < 1 or the ISCO point from SXS:BBH:0066), compute the 4PM periastron advance Delta Phi_4 from Eqs. (7)-(8), first using p_inf^2 = |x| e^{+i pi} (principal log: log|x| + i pi, powers with arg +pi) and then p_inf^2 = |x| e^{-i pi} (log|x| - i pi, powers with arg -pi), with all other inputs identical. Compare the real and imaginary parts. If the real parts differ or if either branch gives a nonzero imaginary part, the dictionary is not uniquely specified by the paper's rule; report which branch reproduces the known 3PM bound-state result when P4 is set to zero, and which (if any) matches the SXS binding-energy curve.","verdict_should_be":"UNCHANGED","load_bearing_attack":"At the core of the revision is Eqs. (5)-(6): P4 contains inverse powers of p_inf^2 and transcendental functions (log, Li2, arccosh) evaluated at p_inf^2, so continuation from scattering (p_inf^2 > 0) to bound states (p_inf^2 < 0) requires fixing a branch for these functions. The paper states that p_r^2 is analytic in the upper half p_inf^2 plane and claims this is accomplished by substituting p_inf^2 -> p_inf^2 e^{-i pi}. But e^{-i pi} = -1; the substitution alone does not distinguish the two sides of the negative real axis. The observable depends on whether log(p_inf^2) is evaluated as log|x| + i pi (arg = +pi, upper side) or log|x| - i pi (arg = -pi, lower side), and similarly for non-integer powers. The paper never computes both branches and shows that only its choice yields real Delta Phi_4; it also never demonstrates that the new dictionary, truncated at 3PM (P4 = 0), reproduces the established 3PM bound-state results obtained by the original Kalin-Porto dictionary, which is a mandatory consistency check since no singularity appears at that order. The only numerical evidence is the SXS comparison, which is thin and partly self-referenced: the EOB model in Refs. [50-53] is built with the same continuation, the figure caption describes the NR point as the 'final state' rather than a well-defined ISCO binding-energy value, and Ref. [56] (needed for the h_i coefficients in Eq. (6)) is listed as 'To appear 111, 111 (2024)' with no public details, so the 4PM expression cannot be independently checked.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper addresses the failure of the Kalin-Porto dictionary when mapping scattering data to bound-state observables at 4PM order. The authors identify a singular factor (p_inf^2)^{-n/2} in the Fourier transform of the scattering amplitude, argue that the original dictionary, which connects scattering and bound states at p_inf^2=0, becomes invalid at 4PM, and propose an analytic continuation p_inf^2 -> p_inf^2 e^{-i pi}, motivated by an analogy with Hawking's treatment of black hole radiation. They then show that the new dictionary restores the expected monotonic behavior of the 2PM and 4PM precession angles and compare EOB binding energies with SXS numerical-relativity data for mass ratios q=1 and q=10, reporting good agreement.","tokens_in":10696,"tokens_out":4606,"duration_ms":46526,"significance":"If the proposed continuation is correct, the paper would resolve a concrete obstruction to using 4PM scattering information for bound orbits, which is important for PM-based waveform models and for the boundary-to-bound program generally. The paper's diagnosis of a singular factor at p_inf^2=0 is concrete and reasonable, and the figures showing that the original dictionary gives complex precession angles illustrate a genuine problem. However, the central prescription is justified more by assertion than by proof: the branch choice for multivalued functions is not specified or tested, the 3PM consistency limit is not checked, and the numerical validation relies on the author's own EOB framework without an independent calibration. These gaps are load-bearing for the central claim and need to be addressed before the result can be accepted.","major_comments":[{"comment":"The continuation p_inf^2 -> p_inf^2 e^{-i pi} does not by itself fix the branches of the multivalued functions appearing in P4, namely log(p_inf^2), Li2, arccosh, and non-integer powers. Since e^{-i pi} = e^{+i pi} = -1, the substitution alone does not distinguish the two sides of the negative real axis; the result depends on whether log(-|x|) is assigned +i pi or -i pi, and similarly for the other functions. The paper asserts that p_r^2 is analytic in the upper half p_inf^2 plane and that branch cuts lie in the lower half-plane, which would suggest approaching the negative real axis from above (arg = +pi), yet the text's notation e^{-i pi} suggests the opposite. The paper never computes both branch choices and shows that only the chosen one yields a real, physical 4PM precession angle. This ambiguity must be resolved with an explicit principal-branch definition and a demonstration that the chosen continuation is unique or physically selected.","section":"Sec. III, Eqs. (5)-(6)"},{"comment":"No consistency check at 3PM is provided. At 3PM the singular factor (p_inf^2)^{-n/2} is absent, and the original Kalin-Porto dictionary is known to work; the new dictionary must reduce to the established 3PM bound-state results in the limit P4 -> 0. The paper does not show this reduction, which is a mandatory check because the branch rule only affects terms involving multivalued functions that appear at 4PM. The absence of this check leaves open the possibility that the new continuation is only one of several branches and that its success at 4PM is accidental.","section":"Sec. III, after Eq. (8)"},{"comment":"The EOB/NR validation is not independent and the comparison is not apples-to-apples. The EOB model in Refs. [50-53] is constructed within the same PM framework and, according to the text, uses the same continuation, so agreement with SXS data cannot be used as evidence for the branch choice unless the EOB Hamiltonian is independently calibrated from other input. Moreover, the figure compares the EOB innermost stable circular orbit (ISCO) binding energy to the NR 'final state of the black hole,' which is not the same quantity; no error bars, no definition of the NR binding energy, and no description of how the ISCO point is extracted are given. A robust validation would compare well-defined ISCO binding energies (or the full Eb(j) curve) against NR data with stated uncertainties and with an independently constructed EOB model.","section":"Sec. III, Fig. 6 and Refs. [50-53]"},{"comment":"The coefficients h_i appearing in T4p are referred to Ref. [56], which is listed as 'To appear 111, 111 (2024)' with no title, journal, or accessible preprint. Since Eq. (6) is the core 4PM input and the h_i are necessary to reproduce the claimed results, this reference must be replaced by a publicly available version or the expressions must be given explicitly. As it stands, the central 4PM calculation cannot be independently checked.","section":"Eq. (6) and Ref. [56]"},{"comment":"The claim that p_r^2 is analytic in the upper half p_inf^2 plane is asserted rather than proved. The Fourier transform in Eq. (3) is defined for p_inf^2 > 0, and no argument is given that the resulting coefficients P_n, which involve logarithms, dilogarithms, and arccosh functions, extend analytically to the upper half-plane without acquiring additional branch cuts. Since the entire continuation prescription rests on this analyticity statement, a proof or a precise reference establishing it is required.","section":"Sec. III, paragraph on analytic continuation"}],"minor_comments":[{"comment":"The manuscript contains numerous grammatical errors and typos, e.g., 'These work lie', 'are should exhibit', 'one shows', and 'detial'; these should be corrected in a thorough language edit.","section":"Throughout"},{"comment":"Several figures are hard to read: Fig. 5 appears to lack axis labels or definitions of the plotted quantities, and Fig. 6's caption does not specify what quantity is plotted on each axis or how the NR 'final state' point is defined. The figures should be redrawn with clear labels, units, and error bars.","section":"Figs. 5 and 6"},{"comment":"The text states that the Fourier transform includes a factor of (p_inf^2)^{-n/2}, but this factor is not displayed explicitly in Eq. (3); the authors should indicate where in the derivation of P_n this factor arises, since it is the central object of the paper.","section":"Eq. (3)"},{"comment":"The analogy with Hawking's black-hole radiation method is presented only briefly; the paper should explain more concretely how the Bogoliubov-coefficient continuation maps onto the continuation of P_n, rather than relying on a general analogy.","section":"Sec. III"}],"recommendation":"major_revision","confidential_remarks":"The paper addresses a real and timely problem, and the identification of the singular factor is a useful observation. However, the proposed dictionary revision is underdetermined as stated, and the numerical validation is partly self-referential. The branch ambiguity and the missing 3PM limit are the main technical blockers; the reliance on an inaccessible reference and the figure quality should also be fixed before resubmission."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper identifies a genuine defect and proposes a plausible fix, but the analytic continuation is not pinned down tightly enough to accept as is. The concrete claim, that the (p_inf^2)^(-n/2) factor in the Fourier transform makes the standard Kalin-Porto dictionary produce complex 4PM precession angles, is new to me and is supported by the figures. The proposed e^(-i pi) dictionary entry may well be right, but the paper does not demonstrate it rigorously.\n\nThe genuinely new elements are the singularity diagnosis and the e^(-i pi) entry, which I do not find in the cited prior literature. The paper also shows concretely that the old dictionary gives a decreasing 2PM precession angle as the binary tightens, which is unphysical. That is worth knowing.\n\nThe main soft spot is the branch choice. e^(-i pi) = -1, so the substitution alone does not say which side of the negative real p_inf^2 axis you are on. The text says p_r^2 is analytic in the upper half-plane, which would favor +i pi, not -i pi; the two branches generally give different complex results because of logs and Li2 functions. The paper never computes both branches and shows only one is real, nor does it check the 3PM limit (P4 = 0) against the known Kalin-Porto results. That consistency test should be mandatory and is easy to do.\n\nThe validation is thinner than the text suggests. The EOB model in Refs. [50-53] is the author's own, built with the same continuation, so the SXS comparison is partly a self-consistency check. Figure 6 labels the NR marker as the final state of the black hole, not an ISCO binding-energy value, and there are only two mass ratios. Ref. [56], which supplies the h_i coefficients in Eq. (6), is listed as \"To appear\" with no public details, so the 4PM expression cannot be independently checked. The Hawking analogy is evocative but does not add rigor.\n\nThe core complaint against the old dictionary is well supported, the proposed continuation is testable, and if it holds, it removes an obstruction to using 4PM scattering data for bound-orbit observables. This is aimed at the PM/EOB/waveform community. It deserves a serious referee: the issue is real and the resolution is checkable. But the paper needs major revision: pin the branch (or show both and explain), add the 3PM consistency limit, make Eq. (6) independently checkable, and either release code/data or give error bars and a clean ISCO definition. I would not desk-reject it; I would send it for review with those requests.","headline":"Real defect, plausible fix, but branch choice underjustified and validation partly self-referenced; send to review with major revision.","tokens_in":11211,"tokens_out":3880,"would_cite":false,"duration_ms":37799,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["03.65.Nk","04.25.Nx","04.20.-q","04.20.Cv"],"model":"deepseek-v4-flash","headline":"The standard dictionary connecting scattering and bound orbits breaks at 4PM; one analytic continuation repairs it.","keywords":["post-Minkowskian expansion","boundary-to-bound dictionary","scattering angle","periastron advance","analytic continuation","binding energy","effective one-body theory","numerical relativity"],"falsifier":"The central claim would be settled by computing the 4PM periastron advance and binding energy with the opposite continuation $p_\\infty^2 e^{+i\\pi}$, and with an explicit branch-cut contribution, then comparing both sets against high-precision numerical-relativity data for small-eccentricity binaries; if the alternative branch fits the data where $e^{-i\\pi}$ does not, the proposed dictionary is wrong, and if neither branch fits, the analytic-continuation picture itself needs revision.","tokens_in":10064,"feed_emoji":"🕳️","tokens_out":11264,"duration_ms":97099,"temperature":0.7,"pith_summary":"The paper claims that the standard dictionary for converting gravitational two-body scattering results into bound-orbit observables—replace $\\beta$ by $i\\beta$ and impact parameter $b$ by $\\pm i|b|$—fails at fourth post-Minkowskian (4PM) order. The reason is that the Fourier transform of the scattering amplitude carries a factor $(p_\\infty^2)^{-n/2}$, so the radial momentum is singular at $p_\\infty^2=0$, exactly the point where the dictionary tries to join scattering states to bound states. The proposed fix is to analytically continue around that singularity as $p_\\infty^2 \\to p_\\infty^2 e^{-i\\pi}$, the same style of continuation used in black-hole radiation. With this substitution the 4PM periastron advance becomes real rather than complex, the scattering and precession angles move in the same direction as the binary separation changes, and effective-one-body binding energies agree with numerical-relativity data for mass ratios $q=1$ and $q=10$. If correct, the fix keeps scattering-derived data usable for bound-orbit and waveform modeling without ad hoc omissions of problematic terms.","feed_headline":"One analytic continuation repairs the bound-orbit dictionary at 4PM","feed_subtitle":"The new continuation yields real precession angles and binding energies that match numerical-relativity data.","key_machinery":"The load-bearing object is the radial momentum squared $p_r^2$ of the binary, expressed at 4PM as $p_r^2 = P_0 r^2 - J^2/r^2 + P_1 G/r + P_2 (G/r)^2 + P_3 (G/r)^3 + P_4 (G/r)^4$, where the coefficients $P_n$ come from the Fourier transform of the scattering amplitude and carry factors $(p_\\infty^2)^{-n/2}$. The singularity at $p_\\infty^2=0$ is why any dictionary based on $\\beta\\to i\\beta$ fails at this order. The proposed replacement is to analytically continue $p_\\infty^2$ around the singularity as $p_\\infty^2\\to p_\\infty^2 e^{-i\\pi}$, the same style of continuation used in black-hole radiation; this choice places branch cuts in the lower half-plane, keeps $p_r^2$ analytic in the upper half-plane, and supplies a corrected dictionary entry $b=J/p_\\infty \\to \\pm i|b|$.","core_discovery":"The central discovery is a diagnosis and a cure. The radial momentum squared $p_r^2$ for the two-body system, built from Fourier-transformed scattering amplitudes, contains coefficients $P_n$ that scale as $(p_\\infty^2)^{-n/2}$; at $n=4$ this introduces a genuine singularity at $p_\\infty^2=0$. The original dictionary assumes that the hyperbolic parameter $\\gamma=\\cosh\\beta$ for scattering can be connected to the bound-state parameter $\\gamma=\\cos\\beta$ through $\\beta\\to i\\beta$ at $\\beta=0$, i.e. at the singular point, and this invalid connection produces complex periastron advances and even a 2PM precession angle that decreases as the separation shrinks. The paper argues that the correct connection is instead to continue $p_\\infty^2$ from positive to negative values around the lower half-plane, $p_\\infty^2 \\to p_\\infty^2 e^{-i\\pi}$, which respects analyticity in the upper half-plane. This revised dictionary yields real bound-state quantities, consistent trends for the scattering and precession angles, and binding-energy curves that track numerical-relativity data.","pith_inferences":["If the $e^{-i\\pi}$ sheet choice is right, the same continuation should also cure complex-valued spin and tidal dictionary entries whenever their Fourier coefficients first acquire a $(p_\\infty^2)^{-n/2}$ singularity; the paper does not test that extension.","The prescription could be tested independently by computing the same bound observables with the alternative continuation $p_\\infty^2 e^{+i\\pi}$ and checking whether realness and agreement with numerical relativity uniquely select $e^{-i\\pi}$; the paper argues analyticity in the upper half-plane but does not rule out extra branch-cut contributions.","A natural next check is to compare the corrected 4PM binding energy against high-precision numerical relativity for additional mass ratios and for small nonzero eccentricities, making the comparison more than a sparse spot check."],"forward_implications":["At 4PM order, the periastron advance computed with the corrected dictionary is real for all $\\beta$, removing the complex unphysical values produced by the old dictionary.","With the corrected dictionary, the 2PM and 4PM precession angles increase as the binary separation decreases, matching the behavior of the scattering angle required by classical expectations.","Effective-one-body binding-energy curves built from 4PM scattering data and the new continuation reproduce numerical-relativity data at the innermost stable circular orbit for mass ratios $q=1$ and $q=10$.","Post-Minkowskian waveform models no longer need to manually discard, modify, or substitute problematic terms such as $\\mathrm{Li}_2$ or $\\log(\\gamma^2-1)$ to avoid complex results."],"supporting_citations":[{"why":"Supplies the original dictionary ($\\beta\\to i\\beta$, $b\\to\\pm i|b|$) that the paper identifies as failing at 4PM.","marker":"[32]"},{"why":"Provide the conservative Hamiltonian and the coefficients $P_n$ whose Fourier transform introduces $(p_\\infty^2)^{-n/2}$.","marker":"[54, 55]"},{"why":"Supplies the analytic-continuation method around a singularity used to justify $p_\\infty^2\\to p_\\infty^2 e^{-i\\pi}$.","marker":"[57]"},{"why":"Defines the $h_i$ coefficients entering the 4PM term $T_{4p}$ that must be continued.","marker":"[56]"},{"why":"Effective-one-body construction used to compute binding energy from the corrected 4PM dictionary for comparison with numerical relativity.","marker":"[50-53]"},{"why":"Establish binding energy versus angular momentum as a key waveform ingredient and provide the numerical-relativity comparison baseline.","marker":"[58, 59]"},{"why":"A prior effective-one-body waveform model that had to modify or omit terms to avoid complex results, illustrating the practical problem the fix solves.","marker":"[48]"}],"fun_headline_variants":["Analytic continuation fixes 4PM bound-orbit dictionary","Revised dictionary cures complex bound-state observables","Bound-state mapping repaired by half-plane continuation","Scattering-to-bound dictionary saved at 4PM","New continuation yields real precession and binding energy"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument depends on assuming that going around the singularity in the direction $p_\\infty^2 \\to p_\\infty^2 e^{-i\\pi}$, rather than the opposite direction, gives the physical bound-state quantities; this sheet choice is asserted rather than proved.","fun_headline_variants_meta":{"raw":{"variants":["Analytic continuation fixes 4PM bound-orbit dictionary","Revised dictionary cures complex bound-state observables","Bound-state mapping repaired by half-plane continuation","Scattering-to-bound dictionary saved at 4PM","New continuation yields real precession and binding energy"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000883,"raw_usage":{"total_tokens":3867,"prompt_tokens":1052,"completion_tokens":2815,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":668,"completion_tokens_details":{"reasoning_tokens":2742}},"tokens_in":668,"tokens_out":2815,"duration_ms":18826,"temperature":1.0,"reasoning_tokens":2742,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T21:41:42.026128+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"The central claim would be settled by computing the 4PM periastron advance and binding energy with the opposite continuation $p_\\infty^2 e^{+i\\pi}$, and with an explicit branch-cut contribution, then comparing both sets against high-precision numerical-relativity data for small-eccentricity binaries; if the alternative branch fits the data where $e^{-i\\pi}$ does not, the proposed dictionary is wrong, and if neither branch fits, the analytic-continuation picture itself needs revision.","supporting_citations":[{"cited_title":"Jing, Effective one-body theory of spinless binary evolution dynamics, To appear 111, 111 (2024)","cited_arxiv_id":null,"evidence_quote":"Supplies the analytic-continuation method around a singularity used to justify $p_\\infty^2\\to p_\\infty^2 e^{-i\\pi}$."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the $h_i$ coefficients entering the 4PM term $T_{4p}$ that must be continued."}],"review_version":1}