{"id":"2ca3ab0b-48fa-42de-87f8-866b44450044","arxiv_id":"2607.05527","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Data-constrained NRQCD form factors yield R(χ_c0)=0.185(3), R(χ_c1)=0.147(26), R(h_c)=0.068(2) and related B_c branching fractions for P-wave charmonia.","lead":"The paper extracts P-wave charmonium wave-function derivatives from radiative decay data and uses them inside NRQCD to predict B_c semileptonic form factors, branching fractions and lepton-flavour-universality ratios R(χ_c0)=0.185(3), R(χ_c1)=0.147(26), R(h_c)=0.068(2). These numbers give experimenters cleaner Standard-Model targets and background estimates for ongoing LHC B_c analyses.","discovery_kind":"extension","skeptic_critique":{"model":"grok-4.5","headline":"Spin-symmetry transfer of pole-expansion slopes from S-wave to P-wave form factors is the single untested assumption that controls the quoted R ratios.","rationale":"The Reader correctly isolates the untested spin-symmetry transfer of the pole-expansion slopes as the weakest link that underpins the precision of the R ratios. All other ingredients (data-driven extraction of ψ′(0), NRQCD LO+RC expressions, BCL re-parametrisation) are either transparent or secondary once the shapes are fixed. Because the absolute form-factor uncertainties are large, the quoted few-percent errors on R rest almost entirely on the cancellation that follows from using identical α_i, β for equal-J final states. The concrete ±20 % variation test proposed above would quantify how much of that cancellation survives and therefore whether the CONDITIONAL verdict should be tightened or relaxed. No stronger internal inconsistency is present, so the Reader’s verdict stands.","tokens_in":47739,"tokens_out":575,"duration_ms":5302,"concrete_test":"Re-fit the BCL coefficients of Tables 26–28 after independently varying each α_i and β of Eq. 3.59 by ±20 % (a conservative estimate of residual spin-symmetry breaking). Recompute the three R ratios of Table 11; if any central value moves outside its quoted 1σ band, the spin-symmetry assumption is load-bearing and the precision claims must be enlarged.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The headline R(χ_c0)=0.185(3), R(χ_c1)=0.147(26), R(h_c)=0.068(2) (abstract, Table 11) are obtained after the form-factor shapes are fixed by the pole-expansion formula (Eq. 3.59). The slope parameters α_i and the common β that appear in that formula are taken unchanged from the B_c→J/ψ (J=1) and B_c→η_c (J=0) fits of Table 7 under the explicit assumption that “an initial state decaying to different final states with the same total spin J have the same values for the shape parameters” (§3.4, paragraph after Eq. 3.59). Because the absolute normalisations at q^{2}=0 already carry 50–70 % uncertainties (Tables 5–6), the ratios R remain precise only if the q^{2} dependence itself is correctly transferred. If the orbital angular momentum of the P-wave states modifies the slopes by more than a few percent, the cancellation that produces the tiny errors on R fails and the central values shift. No lattice, LCSR or independent model check of this transfer is provided.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.5","summary":"The paper computes semileptonic and non-leptonic B_c decays to η_c and the P-wave charmonia χ_c0, χ_c1 and h_c within NRQCD. Analytic form-factor expressions (including O(v^{2}) relativistic corrections) are taken from the literature; their normalisations at q^{2}=0 are fixed by extracting the radial wave-function derivatives ψ'_R(0) of the P-wave states from a global χ^{2} fit to radiative branching fractions (χ_c0→γγ, χ_c1→ργ/ϕγ/J/ψγ, h_c→η_cγ) together with lattice decay constants. The q^{2} shapes are obtained by a pole expansion whose slope parameters α_i and β are transferred from earlier B_c→J/ψ and B_c→η_c analyses under heavy-quark spin symmetry for equal total angular momentum J, then refined by a BCL z-expansion. The resulting form factors yield the LFU ratios R(χ_c0)=0.185(3), R(χ_c1)=0.147(26), R(h_c)=0.068(2), absolute branching fractions, selected non-leptonic rates, and production cross-sections in e^{+}e^{-} annihilation and Z decays.","tokens_in":48065,"tokens_out":1453,"duration_ms":10420,"significance":"If the spin-symmetry transfer of the shape parameters is reliable, the work supplies the first largely data-driven SM predictions for the P-wave LFU ratios with percent-level precision on R(χ_c0) and R(h_c). These ratios are experimentally relevant as backgrounds to R(J/ψ) and as independent tests of lepton-flavour universality in the heavy-heavy sector. The transparent χ^{2} extraction of the wave-function derivatives, the explicit quantification of BCL truncation errors, and the systematic exploitation of form-factor correlations that suppress the R uncertainties are genuine strengths. The same non-perturbative inputs are reused for non-leptonic B_c modes and for Z and e^{+}e^{-} production, giving a coherent set of falsifiable numbers that can be confronted with future LHCb, Belle II and lattice results.","major_comments":[{"comment":"§3.4 (paragraph after Eq. 3.59) and Table 7: the pole-expansion slopes α_i and the common β that control the entire q^{2} dependence of the B_c→P-wave form factors are taken unchanged from the B_c→J/ψ (J=1) and B_c→η_c (J=0) fits under the assumption that “an initial state decaying to different final states with the same total spin J have the same values for the shape parameters.” No lattice, LCSR or independent model check of this transfer is provided. Because the absolute normalisations already carry 50–70 % uncertainties (Tables 5–6), the quoted percent-level precision on R(χ_c0) and R(h_c) (Table 11) rests almost entirely on the correctness of the transferred slopes. A quantitative estimate of the residual spin-symmetry breaking (or a variation of α_i, β within a plausible range) is required before the R uncertainties can be regarded as robust.","section":"§3.4, Eq. (3.59), Table 7"},{"comment":"§3.4 and Table 6: NLO QCD corrections in α_s are available only for the B_c→χ_c0 form factors. For the J=1 channels the authors assign an ad-hoc 20 % prior uncertainty to the LO results, motivated solely by the size of the NLO correction found for χ_c0. This choice is not derived from a power-counting argument or from a partial NLO calculation, and it is unclear whether it adequately covers the missing hard-scattering kernels for V, A0, A1 and A2. The absolute branching fractions in Table 10 (and the non-leptonic rates that depend on A0) inherit this uncontrolled theoretical error; a more systematic assessment or an explicit statement that the absolute rates remain LO+RC only is needed.","section":"§3.4, Table 6"}],"minor_comments":[{"comment":"Abstract and §1 claim that shapes are constrained “without introducing any additional model-dependent inputs.” The spin-symmetry transfer of α_i, β is itself a model assumption; the wording should be softened to “with minimal additional model dependence.”","section":"Abstract, §1"},{"comment":"Table 4: the p-value of the global fit is quoted as 0.94 for χ^{2}/dof = 1.705/4. A brief discussion of whether the large theory errors on χ_c1→ργ (∼71 %) dominate the goodness-of-fit would help the reader judge the robustness of the extracted ψ'_R(0).","section":"Table 4"},{"comment":"Eqs. (5.20)–(5.27) and Table 19: the double-charmonium cross-sections are given without a clear statement of the renormalisation-scale variation that is customarily used to estimate higher-order QCD uncertainty. Adding a short scale-variation band would make the comparison with Belle/BaBar more informative.","section":"§5.2"},{"comment":"Several numerical tables (e.g. Tables 13–15) list branching fractions with uncertainties larger than 100 %. While this is a consequence of the large form-factor errors, a one-sentence remark that only the ratios (or the relative hierarchy) are presently meaningful would prevent over-interpretation.","section":"§4.2"},{"comment":"Typographical inconsistencies appear in the notation for the radial wave-function derivatives (ψ'_R versus ψ'R) and in a few equation labels (e.g. the repeated “(5.4)” after Eq. (5.3)). A uniform notation pass would improve readability.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The central technical novelty—data-driven extraction of the P-wave wave-function derivatives—is solid and publishable. The two major comments concern the uncontrolled transfer of shape parameters and the incomplete NLO treatment of the J=1 form factors; both can be addressed by additional numerical checks or by a more cautious error budget without requiring a complete rewrite. I therefore recommend major revision rather than rejection. The paper is a natural fit for JHEP."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The new piece is a simultaneous fit of the three P-wave radial derivatives ψ'_R(0) for χ_c0, χ_c1 and h_c to the radiative data in Table 3, then feeding those numbers into existing NRQCD kernels (Zhu 2018) to get form-factor normalisations at q^{2}=0 and the R ratios in Table 11. That extraction is transparent, the χ^{2} is well-behaved, and the R ratios correctly exploit the strong correlations that cancel most of the large absolute uncertainties. Absolute BRs remain 60-75 % uncertain, as expected; the ratios are the useful deliverables for LHCb/CMS background estimates.\n\nThey also update a handful of non-leptonic modes and Z/e^{+}e^{-} production rates with the same wave functions. The BCL truncation errors are quantified and small, the unitarity bounds are imposed, and the citation trail to the NRQCD expressions is clean. No invented entities, no circular fitting of the B_c rates themselves.\n\nThe soft spot that actually matters is the shape transfer. After fixing the normalisations they adopt the pole-expansion slopes α_i and β from their earlier B_c\to J/ψ and B_c\toη_c fits (Table 7) under the assumption that equal-J final states share the same q^{2} dependence. That is stated explicitly after Eq. 3.59. Because the absolute form factors already carry 50-70 % errors, the tiny errors on R(χ_c0) and R(h_c) survive only if those slopes are right to a few percent. There is no lattice, LCSR or independent-model check of the transfer. The 20 % ad-hoc NLO uncertainty assigned to the J=1 channels is a secondary, more honest limitation. Both issues are real but do not erase the value of the numbers as first data-driven targets.\n\nThis is for people who need concrete SM numbers for B_c\to P-wave modes right now. It is not a foundational advance, but it is careful enough and useful enough that a serious referee should see it. I would send it out.","headline":"Data-driven P-wave wave-function derivatives give usable R ratios and BRs for B_c, but the quoted precision rides on an untested spin-symmetry transfer of S-wave slopes.","tokens_in":48665,"tokens_out":574,"would_cite":true,"duration_ms":6970,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["13.20.He","12.39.Hg","14.40.Nd"],"model":"grok-4.5","headline":"Data-driven NRQCD form factors for Bc decays to P-wave charmonia yield precise lepton-flavour ratios R(χc0)=0.185(3), R(χc1)=0.147(26) and R(hc)=0.068(2).","keywords":["Bc meson","P-wave charmonia","semileptonic decays","lepton flavour universality","NRQCD form factors","radiative decays","R ratios"],"falsifier":"A lattice-QCD determination of any of the Bc\toχcJ or Bc\to hc form factors at a few q^{2} points, or a first experimental measurement of one of the R(χcJ) or R(hc) ratios, that falls outside the quoted uncertainty bands.","tokens_in":48636,"feed_emoji":"⚛️","tokens_out":974,"duration_ms":7906,"temperature":0.7,"pith_summary":"The paper claims that the form factors governing Bc meson decays into the P-wave charmonia χc0, χc1 and hc can be fixed almost entirely by experiment once their analytic NRQCD expressions are adopted. Radiative-decay measurements of those charmonia are used to extract the derivatives of their radial wave functions at the origin; the resulting normalisations, combined with shapes taken from S-wave channels under heavy-quark spin symmetry, produce controlled predictions for semileptonic branching fractions and the lepton-flavour-universality ratios R. Because cascade decays involving P-wave states are important backgrounds to the cleaner J/ψ modes already under experimental scrutiny, these ratios supply a theoretically clean reference that future LHC measurements can test. The same wave-function inputs also update selected non-leptonic Bc rates and several production cross-sections of P-wave charmonia.","feed_headline":"Bc\to P-wave form factors fixed by data yield R ratios","feed_subtitle":"Radiative decays pin down wave-function derivatives; R(χc0)=0.185(3) and companions follow","key_machinery":"NRQCD factorisation of the Bc\to P-wave transition form factors, with normalisations set by the derivatives of the radial wave functions at the origin that are fitted directly to radiative-decay data and with q^{2} shapes taken from equal-J S-wave channels under heavy-quark spin symmetry.","core_discovery":"By extracting the P-wave radial-wave-function derivatives at the origin from existing radiative-decay data and inserting them into NRQCD form-factor formulae, the authors obtain data-constrained Bc\to(χc0,χc1,hc) form factors whose shapes are fixed by heavy-quark spin symmetry. These form factors yield the lepton-flavour ratios R(χc0)=0.185(3), R(χc1)=0.147(26) and R(hc)=0.068(2) together with the associated semileptonic branching fractions.","pith_inferences":["Once lattice results for any of these form factors appear, the same radiative-decay fit can be repeated with the lattice points as additional constraints, converting the present hybrid prediction into a fully data-driven one.","The hierarchy of precision among the three R ratios already indicates which P-wave channel will give the cleanest new-physics probe when experimental data arrive.","The same extracted wave-function derivatives can be reused for any other exclusive process whose NRQCD factorisation involves the same long-distance matrix elements."],"forward_implications":["Future LHC measurements of R(χc0), R(χc1) or R(hc) can be compared directly with these Standard-Model benchmarks without large model-dependent theory errors.","Cascade backgrounds to the existing R(J/ψ) analyses can be estimated more reliably once the P-wave branching fractions are known.","The same wave-function derivatives update non-leptonic Bc\to P-wave+light-meson rates and several e^{+}e^{-} and Z-decay production cross-sections for P-wave charmonia.","Absolute branching fractions remain limited by the large form-factor uncertainties, so ratio observables will be the first precision tests."],"fun_headline_variants":["Radiative data fix P-wave derivatives for Bc R ratios","Data-constrained Bc→χc form factors yield R(χc0)=0.185(3)","NRQCD FFs from χc0,χc1,hc radiative decays give LFU ratios","Wave-function derivatives update Bc→P-wave form factors and Rs","Bc to χc0,χc1,hc: data-pinned form factors predict R ratios"],"cache_read_input_tokens":32896,"weakest_assumption_plain":"The q^{2}-shape parameters that control how the form factors vary with momentum transfer are assumed identical to those of the corresponding S-wave channels simply because the final states share the same total angular momentum.","fun_headline_variants_meta":{"raw":{"variants":["Radiative data fix P-wave derivatives for Bc R ratios","Data-constrained Bc→χc form factors yield R(χc0)=0.185(3)","NRQCD FFs from χc0,χc1,hc radiative decays give LFU ratios","Wave-function derivatives update Bc→P-wave form factors and Rs","Bc to χc0,χc1,hc: data-pinned form factors predict R ratios"]},"model":"grok-4.5","effort":"low","cost_usd":0.004404,"raw_usage":{"total_tokens":1310,"prompt_tokens":816,"num_sources_used":0,"completion_tokens":115,"cost_in_usd_ticks":44040000,"prompt_tokens_details":{"text_tokens":816,"audio_tokens":0,"image_tokens":0,"cached_tokens":128},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":379,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":816,"tokens_out":115,"duration_ms":3717,"temperature":1.0,"reasoning_tokens":379,"cache_read_input_tokens":128,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-11T06:17:58.553471+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"A lattice-QCD determination of any of the Bc\toχcJ or Bc\to hc form factors at a few q^{2} points, or a first experimental measurement of one of the R(χcJ) or R(hc) ratios, that falls outside the quoted uncertainty bands.","supporting_citations":[],"review_version":1}