{"id":"ba1fc6ba-4b77-44f3-98a5-8d3c319b7702","arxiv_id":"2607.24432","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":5,"one_line_summary":"Full final-state kinematic expansion yields scheme-independent O(v²) corrections of about −14.6% (cc̄) and +12% (non-cc̄) for inclusive J/ψ at B factories, yet tension with Belle data remains.","lead":"A complete kinematic expansion for three-body J/ψ production flips the size and sign pattern of O(v²) relativistic corrections relative to earlier incomplete calculations. The revised rates still leave a theory–experiment gap at B factories, pointing to missing higher-order NRQCD pieces.","discovery_kind":"extension","skeptic_critique":{"model":"moonshotai/kimi-k3","headline":"The dropped O(q²) phase-space Jacobian (Eq. 18) is of the same 1/s order as the z_i^(1) amplitude terms (Eqs. 16–21) the paper keeps and credits with restoring scheme independence; neglecting it is a power-counting inconsistency that is asserted, not demonstrated, to be numerically small.","rationale":"The reader's weakest_assumption identified exactly this locus — the neglect of the O(q²) phase-space Jacobian — and set CONDITIONAL pending its quantification. My pass sharpens the same concern into an internal power-counting argument: the kept z_i-derivative terms and the dropped Jacobian are both linear in z_i^(1) ∝ 1/s, so the \"suppressed by 1/s\" justification, taken at face value, would also argue for dropping the terms the paper presents as essential. The distinction can only be numerical, and no number is provided. This does not rise to REJECT: the construction is otherwise well-motivated, the scheme-independence demonstration in the Table I footnote is a genuine and checkable advance, the fragmentation cross-check in Fig. 1 provides independent support for the retained terms, and the dropped Jacobian may well turn out to be small (it is a total-derivative-like structure that could integrate to little). The honest verdict remains CONDITIONAL with the condition being exactly the test above. I add the percentage-inconsistency observation (abstract vs. Sec. III vs. summary vs. Table I) as a secondary, easily-fixed item that should be resolved in the same revision, since the paper's headline phenomenological numbers currently differ across its own sections.","tokens_in":16210,"tokens_out":3290,"duration_ms":120844,"concrete_test":"Evaluate the dropped Jacobian analytically or numerically: compute J₁ ≡ ∂z_3^(1)/∂z_3^(0) + ∂z_4^(1)/∂z_4^(0) from Eqs. 16–17 as a function of (z_3^(0), z_4^(0)), then integrate ⟨v²⟩·|M(0)|²·J₁ over the LO phase space of Eq. 8 at √s = 10.58 GeV for both channels and add the result to the O(v²) columns of Table I. If the O(v²) entries move by more than ~10% of their stated values (i.e., ≳1.5 percentage points on the 14.55%/12.18% figures), the neglect is unjustified and the headline percentages and the Fig. 1 fragmentation comparison must be redone. Simultaneously, reconcile the four different quoted percentage pairs (abstract, Sec. III, summary, Table I) against the table entries.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim is that adding the z_i-derivative terms in Eqs. 20–21 yields the \"complete\" O(v²) correction, scheme-independent and consistent with fragmentation. Those terms enter through z_i^(1) = 8ba_i/s (Eq. 16), i.e., they are explicitly O(1/s)-suppressed kinematic-expansion effects — precisely the class of effects the paper argues must be retained for consistency (the footnote to Table I shows they contribute −25.9 to −28.0 fb, dominating the total correction). Yet in Eq. 18 the authors drop the O(q²) Jacobian of the three-body phase space, ∂z_3^(1)/∂z_3^(0) + ∂z_4^(1)/∂z_4^(0), on the grounds that it is \"suppressed by the s term in the denominators of z_i^(1).\" But that suppression is identical in power counting to the z_i^(1) terms they keep: both are linear in the same z_i^(1) coefficients. Keeping one and dropping the other is only defensible if the Jacobian piece is numerically negligible after phase-space integration, and at B-factory energies 4m_c²/s ≈ 0.08 is not parametrically tiny. No estimate, bound, or number is given for the dropped term — it is omitted between Eq. 18 and Eq. 19 without evaluation. If the Jacobian contribution is even a few fb, the quoted percentages (14.55%, 12.18%) shift at the level of their quoted precision, and the claimed quantitative agreement with the fragmentation curves in Fig. 1 (which use a two-body expansion with its own phase-space treatment) could be partly accidental. Secondary issue: the quoted percentages are internally inconsistent across the paper — abstract: −14.55% (cc̄), +12.18% (non-cc̄); Sec. III text: +14.6% (non-cc̄), −12.2% (cc̄), i.e., the two numbers are swapped between channels; summary: −14.55%, +12.54%. Table I implies roughly −14.7%/−14.9% and +11.5% depending on μ_r. This looks like editing sloppiness rather than a physics error, but it means the headline numbers cannot currently be traced unambiguously to Table I.","agreement_with_reader":"agree"},"referee_report":{"model":"moonshotai/kimi-k3","summary":"The manuscript revisits O(v²) relativistic corrections to inclusive J/ψ production in e⁺e⁻ annihilation at B-factory energies, focusing on the three-body quark-level subprocesses e⁺e⁻ → J/ψ + c c̄ and e⁺e⁻ → J/ψ + gg. The central technical claim is that previous calculations were incomplete because they expanded the amplitude only in the relative momentum q and E_q while holding final-state kinematic parameters fixed; the authors derive closed-form expansion coefficients for the dimensionless energy fractions z_i (Eqs. 12–17) and add the corresponding ∂M/∂z_i terms to the amplitude expansion (Eqs. 19–21). They demonstrate numerically (Table I footnote) that with these terms the total O(v²) correction is invariant under the choice of phase-space integrand, whereas the q, E_q-only result is not. They find the full O(v²) correction suppresses σ(J/ψ + X_cc̄) by ~14.55% and enhances σ(J/ψ + X_non-cc̄) by ~12%, show high-energy agreement with fragmentation-function calculations (Fig. 1), and conclude that theory–experiment tension persists after combining O(α_s), color-octet, two-photon, and feed-down contributions.","tokens_in":16622,"tokens_out":6504,"duration_ms":222955,"significance":"If the result holds, it resolves a genuine and long-standing technical defect in the literature: prior O(v²) calculations for three-body quarkonium production gave scheme-dependent answers, and this paper identifies the missing z_i-derivative terms as the cure. The strengths are concrete: closed-form kinematic identities valid to arbitrary order in q² (Eqs. 16–17); a direct numerical demonstration of integrand independence, with the partial corrections (−0.692 vs +1.389 fb; −25.942 vs −28.023 fb) canceling to a scheme-independent total (Table I footnote); a falsifiable cross-check against fragmentation functions in both normalization and line shape (Figs. 1–2); and channel-resolved quantitative predictions relevant to the Belle tension. The downward revision of the X_cc̄ relativistic correction from prior values and the reduction of the X_non-cc̄ enhancement from 20–30% to ~12% would materially change the phenomenological ledger for prompt J/ψ production.","major_comments":[{"comment":"Eq. 18 and the paragraph following it: the O(q²) phase-space Jacobian term (∂z₃⁽¹⁾/∂z₃⁽⁰⁾ + ∂z₄⁽¹⁾/∂z₄⁽⁰⁾)𝒒² is dropped as 'suppressed by the s term in the denominators of z_i⁽¹⁾, and ... negligible.' This justification is internally inconsistent: by Eq. (16), z_i⁽¹⁾ = 8ba_i/s, so the z_i-derivative amplitude terms retained in Eqs. (20)–(21) carry exactly the same 1/s suppression. The paper's own Table I footnote shows those equally-suppressed terms contribute −25.9 to −28.0 fb — the dominant piece of the total correction — which refutes '1/s suppressed ⇒ negligible' as a power-counting argument for this class of terms. A consistent O(v²) expansion of ∫dΦ₃(q²)|M(q²)|² contains both the measure term and the amplitude terms at the same order. No estimate, bound, or integrated value of the Jacobian contribution is given; at B-factory energies r = 4m_c²/s ≈ 0.08 this is not parametrically ti","section":"§II, Eq. (18) vs. Eqs. (16), (20)–(21)"},{"comment":"Footnote 3 cites Ref. [49] (Li, Liu, Huang, Sang), described as giving semi-analytical results for J/ψ + X_non-cc̄ through O(α_s v²) with the 'conventional expansion approach' that 'yield negligible contributions to the O(v²) and O(α_s v²) corrections.' If accurate, an independent calculation finding a negligible LO O(v²) correction for the non-cc̄ channel directly contradicts the ~12% enhancement reported here — one of the paper's two headline numbers. Burying this in a footnote is not adequate. The manuscript needs a dedicated comparison: do the authors reproduce/disagree with Ref. [49]'s LO v² SDC, and if they disagree, where (e.g., the z_i-derivative terms, the treatment of the phase-space measure)? This is a correctness-risk issue for the central quantitative claim, not a consensus issue.","section":"Footnote 3 / Ref. [49]; §III"},{"comment":"The two headline percentages are stated inconsistently. The abstract gives −14.55% (X_cc̄) and +12.18% (X_non-cc̄). The first paragraph of §III states the reverse assignment ('increases the LO cross section by 14.6%' for X_non-cc̄ and 'suppresses the LO yield by 12.2%' for X_cc̄). §IV gives 14.55% and 12.54%. Meanwhile Table I implies 0.027/0.183 ≈ 14.8% and 0.018/0.121 ≈ 14.9% for X_cc̄, and 0.038/0.329 ≈ 11.6% and 0.025/0.218 ≈ 11.5% for X_non-cc̄ — none of which equals 14.55%, 12.18%, or 12.54% exactly. For a paper whose deliverable is precision percentages, the channel assignment in §III must be corrected and every quoted number reconciled with Table I (specifying μ_r and rounding).","section":"Abstract; §III first paragraph; §IV; Table I"}],"minor_comments":[{"comment":"The statement that the final-state scattering angles 'are Lorentz scalars, which ensures manifest covariance order-by-order' is misleading: CM-frame angles are frame-defined quantities, not Lorentz invariants. The expansion procedure (fixed angles, uniformly scaled three-momenta) is standard and defensible on its own; please reword rather than claim manifest covariance.","section":"§II, paragraph preceding Eq. (8)"},{"comment":"The value ⟨v²⟩_J/ψ = 0.23 is introduced without a citation or provenance (potential-model estimate? which one?). Please cite the source, and likewise state the source for ⟨v²⟩_χcJ = 0.23 used for Eq. (25).","section":"§III, Eq. (23)"},{"comment":"Table I footnote: 'amount to −0.692 fb for the J/ψ + X_cc̄ channels, respectively' — 'respectively' with a single channel; also state which μ_r the fb-level values correspond to. In the caption of Fig. 1, 'NLOa' is defined via a superscript footnote on the figure, which is easy to miss; define it in the caption text.","section":"Table I footnote; Fig. 1 caption"},{"comment":"The scheme-independence claim rests on a numerical cancellation between two specific integrands (Eq. 8 vs. Refs. [9,13]). An analytic argument (e.g., that the z_i-derivative terms are precisely the chain-rule completion making the expanded integral a total derivative under reparametrization) would substantially strengthen §III and would also clarify whether a retained Jacobian term preserves the invariance.","section":"§III, Table I footnote"},{"comment":"Notation: m_Q appears in Eqs. (20)–(21) while m_c is used everywhere else; Π^{αβ} = −g^{αβ} + p^αp^β/p² uses p where P (the meson momentum) is presumably meant. 'form the derivative terms' → 'from the derivative terms' (Table I footnote); 'both at both LO and NLO' (§III); 'One may tune the values of m_c and the renormalization scale can bring...' (§III) needs grammatical repair.","section":"Eqs. (20)–(21); various"},{"comment":"The χ_cJ cross sections in Eq. (25) appear with minimal setup (inputs only 'same as Ref. [10]'); please state the P-wave LDMEs used and whether the suppression factors 0.661/0.542/0.697 carry any scale dependence. Also, Ref. [48] is cited as 'in preparation' while being invoked in §IV as supporting evidence for the importance of complete expansions in P-wave production — soften this claim or wait for a citable version.","section":"Eq. (25); §IV"},{"comment":"Fig. 2 would benefit from stating the μ_r choice and from indicating the fragmentation-calculation curve for the non-cc̄ channel if available, parallel to what is shown for the cc̄ channel; currently the 'twofold validation' statement in §III references fragmentation only for the cc̄ line shape.","section":"Fig. 2; §III"}],"recommendation":"major_revision","confidential_remarks":"The core idea — that the z_i-derivative terms are the missing chain-rule completion for three-body relativistic corrections — is plausible and, if it survives scrutiny of the dropped phase-space Jacobian (Eq. 18), valuable. The Jacobian issue is not pedantic: the paper's own numbers show that terms with the identical 1/s suppression are the dominant correction, so the omission must be quantified rather than asserted. The unexplained contradiction with the independent semi-analytic calculation of Ref. [49] (footnote 3) is the second thing I would want resolved before recommending acceptance; the authors are best placed to identify whose expansion is incomplete. Neither issue seems unfixable within the manuscript's scope."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The real news is the z_i expansion for three-body final states. Prior B-factory papers (He–Fan–Chao, Jia) expanded only q and E_q; this work keeps the kinematic derivatives of the dimensionless energies, shows the total O(v²) is independent of the phase-space integrand (Table I footnote), and recovers high-energy agreement with fragmentation both in rate and in the shape of the momentum corrections (Fig. 1). That is a legitimate, checkable improvement. The old ~0.5% / 20–30% numbers move to roughly −15% (cc̄) and +12% (non-cc̄) once those terms are in, and the sign flip across the spectrum is physically sensible.\n\nWhat they do well: the kinematic identities (12–17) and the amplitude expansion (19–21) are clean, the scheme-independence check is explicit, and they are honest that tension with Belle survives after stitching in published α_s, CO, two-photon, and feed-down pieces. They also flag the missing joint O(α_s v²) calculation themselves.\n\nSoft spots, in proportion. First, the dropped O(q²) Jacobian in Eq. 18. The z_i^(1) pieces they keep are O(1/s); the Jacobian they discard is built from the same coefficients and the same 1/s. At √s = 10.58 GeV, 4m_c²/s ~ 0.08 is not tiny. They assert it is negligible and never quote a size. That is a power-counting inconsistency, not a fatal hole, but a referee should demand a bound or a numerical estimate—especially since those z_i terms dominate the correction (−26 fb scale). Second, the quoted percentages wander: abstract −14.55% / +12.18%, body text swapped and rounded differently, summary +12.54%, Table I slightly different again. Looks like editing, not physics, but the headline numbers are currently untraceable. Third, the phenomenology is assembled from heterogeneous LDMEs and older NLO pieces; the “tension persists” conclusion is fair for the inputs chosen, not a new prediction.\n\nWho it is for: people who compute or fit NRQCD rates at B factories or who care about consistent O(v²) three-body phase space. Worth a serious referee. I would engage the technical core and ask for the Jacobian estimate plus cleaned numbers. Not a desk reject.","headline":"Solid technical fix to three-body O(v²) expansions—scheme independence and fragmentation matching are real—but they drop a same-order phase-space Jacobian without a number, and the headline percentages are internally inconsistent.","tokens_in":15724,"tokens_out":627,"would_cite":true,"duration_ms":21098,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"A complete expansion of three-body final-state kinematics yields self-consistent O(v²) corrections that suppress J/ψ + open-charm by ~15% and enhance the non-charm channel by ~12%.","keywords":["NRQCD","J/ψ production","relativistic corrections","B factories","three-body phase space","fragmentation functions","color-singlet","O(v²) expansion"],"falsifier":"A fixed-order calculation that keeps the phase-space Jacobian and does not force angular independence of q², or a complete O(α_s v²) evaluation of both channels whose result either closes or deepens the remaining gap with Belle’s measured cross sections.","tokens_in":15236,"feed_emoji":"⚛️","tokens_out":1030,"duration_ms":48017,"temperature":0.7,"pith_summary":"Inclusive J/ψ production at B-factory energies has long shown tension between NRQCD theory and Belle measurements. Earlier relativistic corrections for the three-body quark processes omitted the O(v²) dependence of the final-state energy and momentum variables, so the answers depended on how phase space was integrated. This paper expands every final-state kinematic parameter, obtaining cross sections that are independent of the integration variables and that match fragmentation-function results in both magnitude and momentum-line-shape at high energy. The new O(v²) pieces cut the J/ψ + X_c¯c rate by roughly 14.55% and raise the J/ψ + X_non-c¯c rate by about 12%. After all published QCD corrections, color-octet terms, two-photon channels and feed-down are added, the discrepancy with experiment remains, pointing to the need for complete O(α_s, v²) or higher-order prompt calculations.","feed_headline":"Full O(v²) expansion cuts J/ψ+charm by 15%, lifts non-charm by 12%","feed_subtitle":"Three-body kinematic terms restore scheme independence and match fragmentation, yet Belle tension remains.","key_machinery":"The first-order expansion coefficients z_i^(1) of the dimensionless final-state energy variables, together with the new term they generate in the squared amplitude; that term restores integration-scheme independence and high-energy consistency with fragmentation.","core_discovery":"For three-body quark-level subprocesses that produce J/ψ, expanding every final-state kinematic variable—not only the relative momentum q and the quark energy E_q—produces O(v²) relativistic corrections that are theoretically self-consistent and independent of the choice of phase-space integration variables. Numerically these corrections suppress the J/ψ + X_c¯c cross section by ~14.55% and enhance the J/ψ + X_non-c¯c cross section by ~12%, while both the high-energy cross-section ratios and the shapes of the momentum distributions agree with fragmentation-function calculations.","pith_inferences":["The same z_i expansion should become standard for any three-body NRQCD process at moderate energies where 1/s suppression is not tiny.","If future O(α_s v²) results still undershoot the double-charm channel, production mechanisms beyond ordinary color-singlet plus color-octet may be required.","The sign change of the O(v²) correction across the J/ψ momentum spectrum is a differential signature that improved Belle-II spectra could test directly.","Scheme independence under the full expansion can serve as a practical diagnostic for incomplete relativistic calculations in other exclusive and inclusive quarkonium channels."],"forward_implications":["Earlier O(v²) numbers that kept only q and E_q derivatives must be replaced by the full-expansion values (~–14.55% and ~+12%).","High-energy O(v²) fragmentation results for S-wave quarkonium are now corroborated by fixed-order three-body calculations.","Residual theory–experiment tension shows that incomplete O(α_s, v²) prompt sums are still insufficient.","Color-octet matrix elements fitted at hadron colliders remain inconsistent with B-factory data once the new corrections are included.","Analogous full kinematic expansions are required for other three-body quarkonium processes such as inclusive h_c production."],"fun_headline_variants":["Full kinematic O(v²) expansion cuts J/ψ+c¯c by 14.55%, lifts non-c¯c 12%","Three-body O(v²) terms suppress charm channel ~15%, enhance non-charm ~12%","Complete final-state expansion yields scheme-independent J/ψ O(v²) corrections","O(v²) full kinematics match fragmentation, cut J/ψ+Xc¯c cross section 14.55%","Expanding all three-body variables: −14.55% charm, +12% non-charm for J/ψ"],"cache_read_input_tokens":128,"weakest_assumption_plain":"Outgoing particle angles are treated as independent of the relative quark momentum squared, three-momenta are assumed to scale uniformly, and the O(v²) piece of the three-body phase-space Jacobian is dropped as negligible.","fun_headline_variants_meta":{"raw":{"variants":["Full kinematic O(v²) expansion cuts J/ψ+c¯c by 14.55%, lifts non-c¯c 12%","Three-body O(v²) terms suppress charm channel ~15%, enhance non-charm ~12%","Complete final-state expansion yields scheme-independent J/ψ O(v²) corrections","O(v²) full kinematics match fragmentation, cut J/ψ+Xc¯c cross section 14.55%","Expanding all three-body variables: −14.55% charm, +12% non-charm for J/ψ"]},"model":"grok-4.5","effort":"low","cost_usd":0.005401,"raw_usage":{"total_tokens":1538,"prompt_tokens":848,"num_sources_used":0,"completion_tokens":130,"cost_in_usd_ticks":54008000,"prompt_tokens_details":{"text_tokens":848,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":560,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":848,"tokens_out":130,"duration_ms":9326,"temperature":1.0,"reasoning_tokens":560,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-31T15:01:03.597946+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"A fixed-order calculation that keeps the phase-space Jacobian and does not force angular independence of q², or a complete O(α_s v²) evaluation of both channels whose result either closes or deepens the remaining gap with Belle’s measured cross sections.","supporting_citations":[],"review_version":1}