{"id":"4dbf26fe-cb41-4132-b0ab-7a9137223960","arxiv_id":"2505.20605","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Improved perturbative QCD predictions for B_c -> chi_cJ (P,V) decays match LHCb for the chi_c2 and chi_c1 pion ratios, while the chi_c0 ratio is underpredicted.","lead":"This paper calculates the branching ratios for B_c mesons decaying into charmonium states plus pions, kaons, rho mesons, and K* mesons using an improved version of perturbative QCD. Two predicted ratios match recent LHCb measurements, but a third predicted ratio is about eight times smaller than the value inferred from data.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Central R_chi_c1/chi_c2 agreement relies on a twist-2/twist-3 cancellation in the chi_c1 amplitude whose model normalization is unconstrained.","rationale":"The reader flagged the model distribution amplitudes as the weakest assumption, focusing on f_chi_c0. My stress-test targets the part of that structure that actually carries the central claim, namely the chi_c1 amplitude with its self-reported twist-2/twist-3 cancellation. That the paper itself notes the cancellation is a point in its favor for transparency, but it also marks the exact place where the result becomes an accident of a model normalization. The proposed v^2 variation is a single, cheap analytical check: it leaves decay constants and all other inputs untouched and changes only the shape of the Coulomb-model DAs. The paper's error budget (Eq. (30)) does not cover this. I do not allege any impropriety; the issue is that the central claim is not yet robustly supported. Since the reader already returned CONDITIONAL, my recommendation is UNCHANGED: the paper should be accepted only with the robustness test performed or with the central claim phrased as 'compatible with data for the adopted DAs'.","tokens_in":18042,"tokens_out":5251,"duration_ms":52173,"concrete_test":"Recompute Eqs. (36)-(39) and (46) with v^2 in the model wave function C(x), Eq. (10), varied from 0.3 to 0.1 and 0.5, holding all other inputs fixed. If R_chi_c1/chi_c2 changes by more than its quoted +/- 0.01 or crosses the LHCb bound 0.49, the central 'agreement with data' is not robust to the model DA shape; if the ratio stays within quoted errors, the concern is weakened.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim (abstract; Eqs. (34), (39)) is the agreement of R_chi_c2/J/psi = 0.32 +/- 0.05 and R_chi_c1/chi_c2 = 0.11 +/- 0.01 with LHCb data. The latter rests on BR(B_c -> chi_c1 pi+) = (4.16^{+0.68}_{-0.62}) x 10^{-5} in Eq. (37). The text explicitly identifies a destructive interference between twist-2 and twist-3 contributions in the chi_c1 pi mode (Section III.A). A destructive interference means the amplitude is a difference of comparable terms, so the small BR and hence the ratio are hypersensitive to the relative normalization of the twist-3 distribution amplitude phi^t_chi_c1(x) in Eq. (17), whose coefficient 23.16 is fixed only by a model Coulomb wave function C(x) with v^2 = 0.3. The quoted 10% uncertainty on decay constants (Eq. (30)) rescales all DAs together and does not probe this relative normalization. No independent lattice or QCD-sum-rule determination of these charmonium DAs is cited to pin the twist-3 strength. A moderate change in the DA shape can therefore move R_chi_c1/chi_c2 from 0.11 toward or beyond the LHCb upper limit 0.49, undermining the claimed validation of the iPQCD framework. The chi_c2 ratio is less obviously cancellation-dominated, but it shares the same unvalidated DA-model dependence.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The authors study the two-body weak decays B_c^+ -> chi_cJ (P,V)^+ (J=0,1,2; P=pi,K; V=rho,K*) in the improved perturbative QCD (iPQCD) formalism at leading order in alpha_s. They employ model distribution amplitudes for the chi_cJ mesons taken from earlier work, compute branching ratios, the ratios R_chi_c2/Jpsi and R_chi_c1/chi_c2, longitudinal polarization fractions, and translate BESIII strong-decay branching ratios into predictions for multibody final states under the narrow-width approximation. The headline results are R_chi_c2/Jpsi = 0.32 ± 0.05 and R_chi_c1/chi_c2 = 0.11 ± 0.01, which they claim agree with recent LHCb data; they also report R_chi_c0/Jpsi = 0.17 ± 0.02, which they acknowledge deviates from the LHCb-inferred value.","tokens_in":18429,"tokens_out":6415,"duration_ms":59548,"significance":"If the iPQCD framework is validated, these calculations provide a coherent set of predictions for B_c -> chi_cJ transitions and useful guidance for LHCb searches, including multibody chains. The paper is thorough in comparing with the existing theoretical literature and in exploiting the BESIII measurements. However, the validation is partial: one headline ratio is an upper limit, the chi_c0 ratio fails by a factor of about eight, and the agreement for chi_c1 rests on an interference that is sensitive to an unconstrained twist-3 distribution amplitude. The paper's own suggestion that f_chi_c0 ~ 0.3 GeV would repair the chi_c0 discrepancy shows that the current DA models are not predictive for that channel.","major_comments":[{"comment":"The statement that R_chi_c1/chi_c2 = 0.11 ± 0.01 'agree well with the data within uncertainties' overstates the test: Eq. (2) is a 90% confidence upper limit, not a measured value. A prediction below an upper limit is consistent but cannot validate the framework at the claimed precision. Moreover, the ±0.01 error includes only variations of the inputs beta_Bc, decay constants, Gegenbauer moments, and CKM elements at fixed leading order; no renormalization-scale variation or next-to-leading-order estimate is provided. Please rephrase the comparison as 'satisfies the current upper limit' and add an estimate of the perturbative uncertainty.","section":"Sec. III.A, Eq. (39)"},{"comment":"The smallness of BR(B_c -> chi_c1 pi+) = (4.16^{+0.68}_{-0.62}) x 10^{-5} is produced by a destructive interference between twist-2 and twist-3 contributions, as the paper itself emphasizes. The relative normalization of the twist-3 amplitude phi^t_chi_c1 is fixed by the model coefficient 23.16 in Eq. (17) together with the Coulomb-potential factor C(x) with v^2 = 0.3; the paper cites no independent lattice or QCD sum-rule determination of this amplitude. Since a difference of two comparable terms is hypersensitive to their relative normalization, the prediction R_chi_c1/chi_c2 = 0.11 is not robust. The uncertainty analysis in Eqs. (32)-(38) rescales all chi_cJ decay constants by a common ±10%, so it does not test this relative normalization. Please provide a sensitivity scan over the twist-3 coefficient or v^2 and report the resulting range of R_chi_c1/chi_c2.","section":"Sec. II, Eqs. (15)-(18); Sec. III.A, Eq. (37)"},{"comment":"The paper explicitly acknowledges that R_chi_c0/Jpsi = 0.17 ± 0.02 lies below the LHCb-inferred value 1.41^{+0.50}_{-0.45} by roughly a factor of eight, and then states that f_chi_c0 around 0.3 GeV would amplify BR(B_c -> chi_c0 pi+) by a factor of ten. Because f_chi_c0 is an input, not predicted, this sensitivity means the framework currently makes no robust prediction for the chi_c0 mode; the range 0.09-0.3 GeV for f_chi_c0 changes the central result by an order of magnitude. The paper should either carry f_chi_c0 as a varying parameter with a prior from independent determinations, or present the chi_c0 discrepancy as evidence that the DA model in Eqs. (7)-(9) needs revision, rather than leaving it as a post hoc fix.","section":"Sec. III.A, Eqs. (43)-(44)"}],"minor_comments":[{"comment":"The phrase 'the relative ratios among the B_c -> (J/psi, chi_c1,c2) pi+ BRs from the iPQCD also match the LHCb measurements' is inaccurate because the chi_c1 constraint is only an upper limit; please adjust the wording.","section":"Sec. I, last paragraph"},{"comment":"The formula for C(x) is typeset in a confusing way; please clarify the power of the numerator, the domain of x, and the role of v^2 = 0.3.","section":"Eq. (10)"},{"comment":"Please specify whether m_b and m_c are pole masses or MSbar masses, since the decay constants are evolved using Eq. (31) and the hard kernels depend on the mass scheme.","section":"Eq. (28)"},{"comment":"The ratio R_chi_c1_K*/rho is quoted as (1.78 ± 0.04) x 10^{-2} with a relative uncertainty of about 2%, yet the text contrasts it with the naive expectation; the uncertainty estimate appears to exclude the dominant DA-model uncertainties. Please state which errors are included.","section":"Sec. III.B, Eq. (83)"},{"comment":"Ref. [36] provides QCD sum-rule estimates of the chi_cJ decay constants; it would be helpful to also quote lattice determinations or at least mention that no lattice results are used, given the 10% uncertainty in Eq. (30) is an assumption.","section":"Ref. [36] and Eq. (30)"},{"comment":"The notation f_cbar c is shorthand for the set of chi_cJ decay constants but is not defined; please define it at first use.","section":"Sec. III.A, Eqs. (32)-(38)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the scope of the journal. The main issue is that the paper's abstract and summary overstate the agreement with data: only R_chi_c2/Jpsi is compared with a measured value, and the chi_c0 discrepancy plus the unconstrained twist-3 normalization are more serious than the 'promising' tone suggests. I recommend major revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's the short version. The paper does something new: it applies the iPQCD framework (previous work by the same group) to B_c -> chi_cJ (P,V) decays, and it produces a large menu of testable numbers. The one genuinely successful comparison is R_chi_c2/Jpsi = 0.32 +/- 0.05 against LHCb's 0.37 +/- 0.06. That is a real achievement in a channel where predictions span three orders of magnitude.\n\nWhat it also does well: the paper is transparent about the chi_c0 failure (predicted 0.17 vs inferred 1.41), it uses external charmonium DAs rather than fitting them, and it goes the extra step of folding in BESIII strong-decay BRs to predict multibody rates. The polarization fractions near unity are useful benchmarks. Citation practice looks fine.\n\nThe soft spots are in the chi_c1 channel. The abstract claims R_chi_c1/chi_c2 = 0.11 +/- 0.01 agrees with data, but that 'data' is an upper limit (<0.49), and almost every other calculation in Table II also satisfies it. More importantly, the small chi_c1 pi+ BR is produced by destructive interference between twist-2 and twist-3 amplitudes of comparable size. That makes the prediction hypersensitive to the relative normalization of the twist-3 DA phi^t_chi_c1, which is fixed only by a Coulomb-wave-function model with v^2 = 0.3. The 10% variation on the decay constants rescales all DAs together and does not probe this relative normalization. So the chi_c1/chi_c2 ratio is not a robust validation of the framework; a moderate DA reshuffling could push it to or beyond the LHCb limit. The stress-test note is right about this.\n\nThe chi_c0 discrepancy is honestly reported, but the suggestion that a larger f_chi_c0 (0.3 GeV) would fix it is post hoc. And all uncertainties are input-parameter variations at leading order; there is no hard-scale variation or higher-order estimate, so the error bars are likely too small. No code or explicit hard-kernel formulas are provided—they refer to Ref. [13]—which makes independent checks harder.\n\nOverall, this is a competent calculation of phenomenological interest, not a breakthrough. The central validation claim needs softening, and the authors should be pushed to estimate the twist-3 DA uncertainty and scale sensitivity. I'd send it to a referee; the flaws are identifiable and the prediction set is useful enough to merit serious scrutiny.","headline":"First iPQCD application to B_c -> chi_cJ(P,V) yields one solid ratio (chi_c2/Jpsi) and a fragile chi_c1/chi_c2 prediction that rests on unvalidated twist-3 DA normalization.","tokens_in":18956,"tokens_out":3018,"would_cite":true,"duration_ms":27972,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["13.25.Hw","12.38.Bx","14.40.Nd"],"model":"deepseek-v4-flash","headline":"The improved perturbative QCD formalism, keeping charm quark mass effects in the Sudakov factor, reproduces the LHCb ratios for B_c^+ → χ_{c2}π^+ and χ_{c1}π^+ decays, and extends to predict χ_{cJ}ρ^+ modes and polarization fractions.","keywords":["B_c meson decays","χ_cJ charmonia","improved perturbative QCD","Sudakov resummation","branching ratio ratios","longitudinal polarization","LHCb","QCD factorization"],"falsifier":"High-precision LHCb measurements of $R_{\\chi_{c2}/J/\\psi}$ and $R_{\\chi_{c1}/\\chi_{c2}}$ with errors below about 20% would settle the claim: if the values move outside $0.32\\pm0.05$ and $0.11\\pm0.01$ respectively by more than the combined uncertainties, the iPQCD prediction is falsified.","tokens_in":17817,"feed_emoji":"⚛️","tokens_out":22552,"duration_ms":166388,"temperature":0.7,"pith_summary":"The paper aims to establish that the improved perturbative QCD (iPQCD) formalism—a version of perturbative QCD that keeps charm quark mass effects in the Sudakov factor—can describe the measured $B_c^+$ decays into $p$-wave charmonia plus a light meson. Its central results are the ratios $R_{\\chi_{c2}/J/\\psi} = 0.32 \\pm 0.05$ and $R_{\\chi_{c1}/\\chi_{c2}} = 0.11 \\pm 0.01$, which agree with LHCb data within uncertainties, supporting the framework's predictive ability. The paper also predicts $R_{\\chi_{c0}/J/\\psi} = 0.17 \\pm 0.02$, lower than the value $1.41^{+0.50}_{-0.45}$ inferred from present data, and attributes this tension to the poorly known $\\chi_{c0}$ decay constant. It further computes the $\\chi_{c1,c2}\\rho^+$ branching ratios, finds nearly unity longitudinal polarization fractions, and identifies a constructive (destructive) twist-2/twist-3 interference in the $\\chi_{c1}\\rho^+$ ($\\chi_{c1}\\pi^+$) mode. If correct, these predictions give concrete targets for future precision measurements and probe the internal structure of charmonium states.","feed_headline":"A QCD calculation matches LHCb on B_c to χ_c2 and χ_c1 pion ratios","feed_subtitle":"It also predicts near-unity polarization and a χ_c0 ratio open to future data.","key_machinery":"The machinery is the iPQCD factorization formula, $A \\sim \\int dx_1 dx_2 dx_3\\, db_1 db_2 db_3\\, \\mathrm{Tr}\\,[C\\, \\Phi_{B_c}\\Phi_{P,V}\\Phi_{\\chi_{cJ}} H\\, e^{-S}]$, in which the Sudakov factor $e^{-S}$ now includes charm quark mass effects, and $\\Phi_{\\chi_{cJ}}$ are model distribution amplitudes for the $p$-wave charmonia (twist-2 and twist-3) built from Coulombic wave functions with small relativistic corrections. The symmetry properties of these amplitudes under $x \\leftrightarrow 1-x$, together with normalization to the decay constants $f_{\\chi_{cJ}}$, control the relative factorizable and nonfactorizable contributions and produce the distinctive ratio pattern, notably the huge $R^{\\chi_{c1}}_{\\rho/\\pi}\\approx 34$ and the near-unity longitudinal polarization fractions.","core_discovery":"The paper's central discovery is that, at leading order in the strong coupling $\\alpha_s$, the improved perturbative QCD formalism reproduces the two LHCb ratios for $p$-wave charmonium production in $B_c^+$ decays: $R_{\\chi_{c2}/J/\\psi} = 0.32 \\pm 0.05$ against the measured $0.37 \\pm 0.06$, and $R_{\\chi_{c1}/\\chi_{c2}} = 0.11 \\pm 0.01$, which satisfies the experimental upper bound $0.49$. The framework also predicts $R_{\\chi_{c0}/J/\\psi} = 0.17 \\pm 0.02$, well below the value $1.41^{+0.50}_{-0.45}$ inferred from LHCb evidence, and the paper traces this tension to the uncertain $\\chi_{c0}$ decay constant: a value around $0.3$ GeV would boost the $\\chi_{c0}\\pi$ branching ratio by a factor of ten. For the vector modes, the paper predicts longitudinal polarization fractions of about $95\\%$ ($\\chi_{c1}\\rho^+$) and $93\\%$ ($\\chi_{c2}\\rho^+$), and identifies a twist-2/twist-3 interference that is constructive in $B_c^+\\to\\chi_{c1}\\rho^+$ and destructive in $B_c^+\\to\\chi_{c1}\\pi^+$, a pattern absent in the $\\chi_{c0}$ and $\\chi_{c2}$ channels.","pith_inferences":["A lattice or sum-rule determination of $f_{\\chi_{c0}}$ and of the $\\chi_{c0}$ twist-2/twist-3 distribution amplitudes would resolve the $\\chi_{c0}$ tension; the paper's reliance on the model amplitude is the fragile link.","If future data confirm $R^{\\chi_{c1}}_{\\rho/\\pi}\\approx34$ (far above the $J/\\psi$ value around 3.15), that would corroborate the identified constructive/destructive twist interference; a value near the $J/\\psi$ ratio would call for revised $\\chi_{c1}$ distribution amplitudes.","The same iPQCD framework could be applied to semileptonic $B_c\\to\\chi_{cJ}\\ell\\nu$ decays or other two-body $B_c$ modes as a consistency check of the Sudakov treatment of charm quark mass.","Because the ratios are computed at leading order in $\\alpha_s$, a next-to-leading-order calculation would test the stability of the agreement; the current match could in principle be a truncation accident."],"forward_implications":["The measured ratios $R_{\\chi_{c2}/J/\\psi}$ and $R_{\\chi_{c1}/\\chi_{c2}}$ become tests of iPQCD: future LHCb data with smaller uncertainties will either confirm or shift the calculation.","The predicted $R_{\\chi_{c0}/J/\\psi}=0.17\\pm0.02$ implies the present LHCb evidence for $B_c^+\\to\\chi_{c0}\\pi^+$ is in tension with the framework unless $f_{\\chi_{c0}}$ is near $0.3$ GeV, which would raise the rate tenfold.","Longitudinal polarization fractions $f_L\\approx95\\%$ ($\\chi_{c1}\\rho^+$) and $93\\%$ ($\\chi_{c2}\\rho^+$) mean angular analyses should find these modes dominated by the longitudinal component.","The predicted multibody branching ratios for $B_c^+\\to\\chi_{cJ}(\\to\\pi^+\\pi^-,\\, K^+K^-,\\, 2\\pi^+\\pi^-,\\, \\pi^+\\pi^-K^+K^-)\\pi^+/\\rho^+$ give concrete, searchable final states at LHCb."],"supporting_citations":[{"why":"The LHCb measurement of the ratio R(χ_c2/J/ψ) and the upper limit on R(χ_c1/χ_c2) that the paper's central ratios must match.","marker":"[4]"},{"why":"LHCb evidence for B_c^+ → χ_c0 π^+ that the paper converts into the experimental value of R(χ_c0/J/ψ) for comparison.","marker":"[5]"},{"why":"Supplies the model distribution amplitudes for χ_cJ and the factorizable/nonfactorizable emission amplitude expressions used in the factorization.","marker":"[13]"},{"why":"Establishes the iPQCD framework and provides the B_c^+ → J/ψ π^+ branching ratio used as the denominator for the central ratios.","marker":"[19]"},{"why":"Derives the Sudakov resummation with charm quark mass effects that distinguishes iPQCD from conventional PQCD.","marker":"[20]"},{"why":"BESIII measurement of the strong decay branching ratios χ_c0,c2 → π^+π^-/K^+K^- used for the secondary decay chain predictions.","marker":"[23]"},{"why":"Particle Data Group inputs: masses, CKM elements, B_c lifetime, and strong decay branching ratios of χ_cJ.","marker":"[24]"},{"why":"BESIII measurement of χ_cJ → 2(π^+π^-) and π^+π^-K^+K^- branching ratios used for the multibody predictions.","marker":"[25]"},{"why":"Provides the χ_cJ decay constants whose values, especially f(χ_c0), control the predicted branching ratios.","marker":"[36]"}],"fun_headline_variants":["iPQCD reproduces LHCb χc2/χc1 ratios, predicts χc0 lower","χc0 ratio tension: iPQCD predicts 0.17 vs LHCb's 1.41","Twist-2/3 interference separates χc1 from χc0,χc2 in Bc decays","Bc→χcJ decays: iPQCD matches LHCb ratios, predicts polarization","iPQCD predicts ~95% longitudinal polarization for Bc→χc1ρ"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the model distribution amplitudes for the $\\chi_{cJ}$ mesons—especially the $\\chi_{c0}$ amplitude with decay constant $f_{\\chi_{c0}}=0.093$ GeV—are accurate enough, since the paper shows that changing $f_{\\chi_{c0}}$ to about $0.3$ GeV multiplies the $\\chi_{c0}\\pi$ branching ratio by ten.","fun_headline_variants_meta":{"raw":{"variants":["iPQCD reproduces LHCb χc2/χc1 ratios, predicts χc0 lower","χc0 ratio tension: iPQCD predicts 0.17 vs LHCb's 1.41","Twist-2/3 interference separates χc1 from χc0,χc2 in Bc decays","Bc→χcJ decays: iPQCD matches LHCb ratios, predicts polarization","iPQCD predicts ~95% longitudinal polarization for Bc→χc1ρ"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001339,"raw_usage":{"total_tokens":5734,"prompt_tokens":1530,"completion_tokens":4204,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":1146,"completion_tokens_details":{"reasoning_tokens":4074}},"tokens_in":1146,"tokens_out":4204,"duration_ms":31372,"temperature":1.0,"reasoning_tokens":4074,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T13:51:05.987764+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"High-precision LHCb measurements of $R_{\\chi_{c2}/J/\\psi}$ and $R_{\\chi_{c1}/\\chi_{c2}}$ with errors below about 20% would settle the claim: if the values move outside $0.32\\pm0.05$ and $0.11\\pm0.01$ respectively by more than the combined uncertainties, the iPQCD prediction is falsified.","supporting_citations":[{"cited_title":"Aaijet al.(LHCb Collaboration), J","cited_arxiv_id":null,"evidence_quote":"The LHCb measurement of the ratio R(χ_c2/J/ψ) and the upper limit on R(χ_c1/χ_c2) that the paper's central ratios must match."},{"cited_title":"Aaijet al.(LHCb Collaboration), Phys","cited_arxiv_id":null,"evidence_quote":"LHCb evidence for B_c^+ → χ_c0 π^+ that the paper converts into the experimental value of R(χ_c0/J/ψ) for comparison."},{"cited_title":"Rui, Phys","cited_arxiv_id":null,"evidence_quote":"Supplies the model distribution amplitudes for χ_cJ and the factorizable/nonfactorizable emission amplitude expressions used in the factorization."},{"cited_title":"Liu, Phys","cited_arxiv_id":null,"evidence_quote":"Establishes the iPQCD framework and provides the B_c^+ → J/ψ π^+ branching ratio used as the denominator for the central ratios."},{"cited_title":"Liu, H.-n","cited_arxiv_id":null,"evidence_quote":"Derives the Sudakov resummation with charm quark mass effects that distinguishes iPQCD from conventional PQCD."},{"cited_title":"Ablikimet al.(BESIII Collaboration), arXiv:2502.08929 [hep-ex]","cited_arxiv_id":null,"evidence_quote":"BESIII measurement of the strong decay branching ratios χ_c0,c2 → π^+π^-/K^+K^- used for the secondary decay chain predictions."},{"cited_title":"Ablikimet al.(BESIII Collaboration), Phys","cited_arxiv_id":null,"evidence_quote":"BESIII measurement of χ_cJ → 2(π^+π^-) and π^+π^-K^+K^- branching ratios used for the multibody predictions."},{"cited_title":"Olpak, A","cited_arxiv_id":null,"evidence_quote":"Provides the χ_cJ decay constants whose values, especially f(χ_c0), control the predicted branching ratios."}],"review_version":1}