{"id":"abf73897-fd74-4a3f-978b-7fec47e5b1be","arxiv_id":"2607.23163","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"A PQCD calculation predicts B(B_c→χc0ππ)=3.24×10⁻³, B(B_c→χc1ππ)=4.19×10⁻³ and R_{χc1/χc0}≈1.30 for rho-mediated decays.","lead":"This paper predicts how often a particle called the B_c meson decays into charmonium plus a pion pair or a kaon-pion pair, produced through short-lived rho or K-star resonances. The predicted rates are large enough for the LHCb experiment to check, and a key ratio provides a new test of the theory.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Central ratio R≈1.30 is driven by a twist-2/twist-3 cancellation in the χc1 amplitude that is not assigned an uncertainty; without a stability test, the headline contrast with the two-body value 4.7 is not established.","rationale":"The reader's weakest_assumption focuses on the absolute normalization (N_ρ, N_K*, F⊥/F∥ ansatz) and the NWA coherent-sum comparison. Those do affect the absolute rates and the Table IV extraction, but the central ratio R — the element that supports the 'dramatically reshapes relative yields' message — is largely independent of N_ρ and F⊥/F∥. My stress-test therefore identifies a different soft spot: the unquantified twist-2/3 cancellation in the χc1 amplitude. This is consistent with the reader's general assessment of medium correctness risk, but targets the most load-bearing quantity more precisely. The concern is not that the paper is wrong; it is that a key number is presented without an uncertainty and may be sensitive to a known model-dependent cancellation. The proposed check (vary φ^t and t, recompute R) would settle whether R is stable. Since the reader already issued CONDITIONAL based on related normalization uncertainties, my recommendation is UNCHANGED: the verdict remains CONDITIONAL, with the added requirement that the robustness of R be demonstrated before the ratio is used as a precision benchmark. I found no grounds for REJECT: the PQCD calculation is internally coherent, the channels are new, and the predictions are falsifiable.","tokens_in":22358,"tokens_out":15809,"duration_ms":141996,"concrete_test":"Recompute B(B_c^+→χc1π+π0) and B(B_c^+→χc0π+π0) with (a) the twist-3 χc1 DA φ^t set to zero, (b) its coefficient 23.16 in Eq. (31) varied by ±20%, and (c) the hard scale t set to 0.9t and 1.1t, as in the tables, and form R from each pair. If R moves outside roughly [1.0, 2.5] in any of these variations, the headline ratio 1.30 is not robust and the contrast with 4.7 needs a caveat; if R stays within that band, the central claim is supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's headline claim is that resonant ρ production flips the relative P-wave charmonium yields: R^ππ_{χc1/χc0} ≈ 1.30 versus the two-body PQCD value ≈ 4.7 (Eq. (81), Table V). This ratio is essentially independent of the fitted normalizations N_ρ/N_K* (they cancel between numerator and denominator) and only weakly affected by the F⊥/F∥ ansatz, since f_L ≈ 94%. What sets R is the relative magnitude of the χc0 and χc1 amplitudes after integrating the two-meson form factor over ω². Naively (f_{χc1}/f_{χc0})² ≈ 4, reduced to ≈ 3.6 by phase space; the drop to 1.30 comes from destructive interference between twist-2 and twist-3 contributions to the χc1 longitudinal amplitude (Eqs. (60)-(61)). This cancellation is not quantified: the paper quotes R without any uncertainty, and the hard-scale variation alone shifts the χc1 rate by +9.11/−8.21 ×10⁻⁴ (≈20%) and the χc0 rate by +6.57/−5.96 ×10⁻⁴ (Tables I-II), which can move R by tens of percent unless correlated. The twist-3 DA φ^t_{χc1} (Eq. (31)) is a model input with no stated error, and the same cancellation is known from Ref. [16] to be highly pattern-dependent. If R were actually 2–3, the qualitative 'reshaping' conclusion would be substantially weakened. This is the most load-bearing point because the absolute rates at 10⁻³ and the NWA comparison in Table IV are secondary; the novel physics claim is the ratio.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper computes quasi-two-body B_c^+→χ_{c0},χ_{c1}[ρ(K^*)→]ππ(Kπ) branching ratios in the leading-order PQCD framework. The ππ and Kπ systems are treated through P-wave two-meson distribution amplitudes with Gounaris-Sakurai and relativistic Breit-Wigner time-like form factors, normalized by factors N_ρ=1.05 and N_{K^*}=1.48 imported from Ref. [42]. For the χ_{cJ} mesons, twist-2 and twist-3 LCDA models from Refs. [6,16] are used. The main numerical outputs are B(B_c^+→χ_{c0}π^+π^0)=3.24×10^{-3}, B(B_c^+→χ_{c1}π^+π^0)=4.19×10^{-3}, f_L≈94% for χ_{c1}, and the central ratio R^{ππ}_{χ_{c1}/χ_{c0}}≈1.30, which is contrasted with the two-body PQCD value ≈4.7. The Kπ channels are predicted at the 10^{-6} level. The paper also extracts two-body B_c^+→χ_{cJ}ρ^+(K^{*+}) branching fractions under the narrow-width approximation and compares them with earlier calculations.","tokens_in":22926,"tokens_out":6893,"duration_ms":66231,"significance":"If the predictions are correct, the paper provides a sharp, falsifiable test of quasi-two-body PQCD factorization in B_c decays into P-wave charmonia. The ratio R^{ππ}_{χ_{c1}/χ_{c0}}≈1.30 versus the direct two-body value ≈4.7 is a genuinely interesting qualitative claim, and the predicted large longitudinal polarization fraction f_L≈94% offers an independent experimental handle. The manuscript is also transparent in presenting the analytic amplitudes and a detailed error decomposition for the branching ratios; no target branching ratio is fitted, which is a strength. However, the headline ratio itself is not assigned an uncertainty, and the central claim rests on a twist-2/twist-3 cancellation that is not stress-tested. The absolute rates also inherit fitted normalization factors whose transferability is not quantitatively examined. These issues are fixable, but they must be addressed before the main claim can be considered established.","major_comments":[{"comment":"The headline ratio R^{ππ}_{χ_{c1}/χ_{c0}}≈1.30 is quoted without an uncertainty, although the underlying amplitudes are subject to the errors listed in Tables I and II. The χ_{c1} longitudinal amplitude is a difference of twist-2 and twist-3 contributions (Eqs. (60)–(61)); the twist-3 DA φ^t_{χ_{c1}} (Eq. (31)) is a model input with no assigned error. The separate hard-scale variations, +9.11/−8.21×10^{-4} for χ_{c1} and +6.57/−5.96×10^{-4} for χ_{c0}, are large compared with the central rates and cannot be assumed to cancel in the ratio. Since the contrast between 1.30 and the two-body value 4.7 is the paper's main claim, the ratio needs a fully propagated uncertainty and a stability test (e.g., varying the coefficients in φ^t_{χ_{c1}}, the parameter v^2 in Eq. (24), and the hard scale in a correlated way). Without this, the qualitative 'reshaping' conclusion is not established.","section":"Eq. (81); Tables I, II, V"},{"comment":"The normalization factors N_ρ=1.05 and N_{K^*}=1.48 are fitted in Ref. [42] for other PQCD quasi-two-body analyses. The paper imports them without a quantitative discussion of channel dependence. They cancel in the ratio R, but they directly set the absolute branching ratios, which are central outputs at the 10^{-3} level. Given that the quoted uncertainties on N_ρ and N_{K^*} are only 3–4%, while the theoretical inputs vary by 20–30%, the authors should either justify transferability within the same PQCD framework or conservatively inflate the normalization uncertainty. A concrete check would be to recompute Tables I and II with N_ρ=1 and N_ρ=1.10, and analogously for N_{K^*}, reporting the shifts separately.","section":"Sec. II.B, Eqs. (48), (51); Tables I, II"},{"comment":"The narrow-width extraction uses the total coherent ρ sum (32.37×10^{-4} for χ_{c0} and 41.92×10^{-4} for χ_{c1}, Tables I and II) as B(B_c^+→χ_{cJ}ρ^+). But for the ρ(770) resonance, the NWA should use the ρ(770)-only entries, 25.86×10^{-4} and 32.62×10^{-4}. Including ρ(1450) and ρ(1700) and their constructive interference inflates the χ_{c0} two-body value by about 25% relative to ρ(770) alone. The comparison with Ref. [16] in Table V is therefore not on a like-for-like basis. Please correct this or explicitly define the Table IV quantity as the coherent ρ-series NWA and compare with correspondingly defined two-body predictions.","section":"Table IV and accompanying text"},{"comment":"The prescription F^∥(ω²)/F^⊥(ω²) ≈ f_V/f_T^V fixes the normalization of all perpendicular two-meson DAs, but no uncertainty is assigned to this relation. These DAs enter the χ_{c0} amplitude through φ^s_{ππ} (Eq. (36)) and hence the χ_{c0} rate and the ratio R. Because the central claim depends on the relative χ_{c1}/χ_{c0} normalization, this assumption needs a quantitative error estimate or a sensitivity scan. A short paragraph reporting the shift in R when f_T^V/f_V is varied by, say, ±10% would suffice.","section":"Sec. II.B after Eq. (53)"}],"minor_comments":[{"comment":"The abstract states that the dominant ρ(770) channels yield B(B_c^+→χ_{c0}π^+π^0)=3.24×10^{-3}, but Table I lists the ρ(770)-only value as 25.86×10^{-4}; the quoted number is the coherent ρ(770)+ρ(1450)+ρ(1700) sum. Please disambiguate the wording.","section":"Abstract; Table I"},{"comment":"The numerical chain 0.053 → 0.055 → 0.110 → 2.0×10^{-3} is not self-explanatory. Please state explicitly which factor is included in each step; the final drop from 0.110 to 2.0×10^{-3} is particularly nontrivial.","section":"Eq. (79)"},{"comment":"The 'Other predictions' columns use footnote letters a–d with references at the end of the table, but no legend is given in the table caption. Add an explicit mapping between the letters and Refs. [72], [17], [73], [74].","section":"Table IV"},{"comment":"Some distribution amplitudes are written with a b-dependence, e.g., φ^L_{χ_{c1}}(x_3,b_3) in Eqs. (60)–(61), while the defining equations in Sec. II.B list only x. Please clarify the convention or add a remark that the b-dependence is implicit through the Sudakov factor.","section":"General notation"}],"recommendation":"major_revision","confidential_remarks":"The paper is a solid PQCD calculation with clear formulas and input transparency, but the main qualitative claim — the ratio R≈1.30 versus 4.7 — is not yet quantitatively supported because the ratio carries no uncertainty and the underlying twist-2/twist-3 cancellation is not stress-tested. The NWA comparison in Table IV also needs correction. These are fixable within the manuscript's scope, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First, this is a genuine new computation. The leading-order PQCD amplitudes for Bc→χc0,c1(ρ/K*→)ππ(Kπ) are written out, the branching ratios and polarization fractions are new, and the central numbers follow from the stated formulas. It's a solid extension of the quasi-two-body program from J/ψππ to P-wave charmonia, and the comparison with two-body iPQCD is exactly the right frame.\n\nWhere it earns credit: the paper doesn't hide its approximations. It states that full phase-space factorization is unavailable, invokes the Watson theorem for absorbing final-state interactions into the time-like form factors, and spells out the DA models in detail. Error bars are propagated from four sources for every branching ratio, which is more than this literature often does. The Kπ suppression to ~10⁻⁶ and the decrease of fL with resonance mass are plausible, testable patterns.\n\nThe soft spots. First, the headline claim—Rχc1/χc0 ≈ 1.30 versus 4.7—is quoted without any uncertainty. The central value sits on a cancellation between twist-2 and twist-3 contributions in the χc1 amplitude; the twist-3 DA φt_χc1 is a model input with no error, and the hard-scale variation alone moves each rate by ~20%. If that cancellation is pattern-dependent, as the same group found in the two-body case, R could easily come out at 2–3. The paper should present a stability scan, or at least assign a defensible error to R, before claiming that resonant ρ production dramatically reshapes the yields. Second, the absolute rates inherit the fitted normalizations Nρ = 1.05 and NK* = 1.48 from Ref. [42], yet those uncertainties never appear in Tables I–II. They may be reasonable fixing factors within one framework, but their channel dependence is untested; for absolute benchmarks they should be varied or argued universal. The F⊥/F∥ ≈ fT_V/fV point is minor, since fL ≈ 94% dominates.\n\nThird, Table IV: the narrow-width extraction uses the coherent ρ(770)+ρ(1450)+ρ(1700) sum and labels it B(Bc→χcJρ+), then compares it with two-body predictions for ρ(770). The comparison should use the ρ(770)-only row, or at least flag the mismatch. The 1.30 ratio survives this (it is ≈1.26 with ρ(770) alone), but the absolute comparison as printed is misleading.\n\nNet: this paper deserves a serious referee. The amplitude work is careful, the channels are unmeasured, and the predictions are usable for LHCb. It needs revision—add an error bar to R, propagate the N factors, and clean up the NWA table—not rethinking.","headline":"Genuinely new PQCD predictions for Bc→χcJππ(Kπ), but the headline ratio R≈1.30 has no quoted uncertainty and rests on a twist-3 cancellation that needs a stability check before the 'reshaping' claim is taken at face value.","tokens_in":23344,"tokens_out":3402,"would_cite":true,"duration_ms":34425,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"B_c decays through a resonant ρ are predicted to produce χ_{c0,c1}ππ at the 10^{-3} level, with the χ_{c1}-to-χ_{c0} ratio near 1.3, far from the two-body prediction of ≈4.7.","keywords":["B_c meson","quasi-two-body decays","P-wave charmonia","PQCD factorization","two-meson distribution amplitudes","Gounaris-Sakurai model","branching ratios","polarization fractions"],"falsifier":"Measure B(B_c^+→χ_{c0}π^+π^0) and B(B_c^+→χ_{c1}π^+π^0) at LHCb in the χ_{cJ}→J/ψγ chain. The central claim predicts their ratio R≈1.30; if the measured ratio instead approaches the two-body PQCD value ≈4.7, then the hypothesis that resonant ρ production reshapes the P-wave charmonium yields through factorizable emission is falsified.","tokens_in":22261,"feed_emoji":"⚛️","tokens_out":10757,"duration_ms":92865,"temperature":0.7,"pith_summary":"The paper applies perturbative QCD (PQCD) factorization to quasi-two-body decays B_c^+ → χ_{c0,c1}[ρ(K^*)→]ππ(Kπ), treating the meson pair via two-meson distribution amplitudes. It predicts B(B_c^+→χ_{c0}π^+π^0)=3.24×10^{-3} and B(B_c^+→χ_{c1}π^+π^0)=4.19×10^{-3}, with ρ(770) alone supplying about 80% of the rate and coherent ρ interference adding about 25%. The headline result is the ratio R_{χc1/χc0}^{ππ}≈1.30, sharply below the two-body PQCD result ≈4.7; the authors argue this shows that resonant ρ production reshapes the relative yields of P-wave charmonia by feeding both states through the same factorizable-emission mechanism. They also predict f_L≈94% for χ_{c1} and Kπ modes at 10^{-6}. These channels should be measurable at LHCb, and the ratio provides a systematic-cancelling test of the factorization framework.","feed_headline":"B_c→χ_cππ rates reach 10^{-3}; ratio near 1.3","feed_subtitle":"The χc1/χc0 yield ratio should sit at ~1.3, far below the two-body PQCD value ~4.7, a clean LHCb test.","key_machinery":"The central machinery is the quasi-two-body PQCD factorization with two-meson P-wave distribution amplitudes (DAs) for the ππ and Kπ pairs. The time-like form factors F^{∥}_{ππ}, F^{⊥}_{ππ}, F^{∥}_{Kπ}, F^{⊥}_{Kπ} are parametrized by the Gounaris-Sakurai model for ρ(770), ρ(1450), ρ(1700) and the relativistic Breit-Wigner for K*(892), with measured normalization factors N_ρ≈1.05 and N_{K*}≈1.48 and the approximation F^{⊥}/F^{∥}≈f^{V}_{T}/f^{V}. For the χ_{cJ} mesons, both twist-2 and twist-3 light-cone DAs are included; the twist-3 contributions grow with the pair invariant mass, which drives the decrease of f_L from 94% (ρ(770)) to 65% (ρ(1700)) and controls the χc1-to-χc0 ratio. The hard k","core_discovery":"On the authors' own terms: in B_c^+→χ_{c0,c1}(ρ→ππ), the intermediate ρ couples to both charmonia through factorizable emission diagrams, so the two states share one production path. The calculation gives B(B_c^+→χ_{c0}π^+π^0)=3.24×10^{-3} and B(B_c^+→χ_{c1}π^+π^0)=4.19×10^{-3}, hence R^{ππ}_{χc1/χc0}≈1.30. The two-body transitions, by contrast, have R≈4.7 because χ_{c0} production is color-singlet suppressed and χ_{c1} suffers from a vanishing leading-twist decay constant. The paper concludes that the near-unity resonant ratio signals the bypassing of those suppressions, and predicts that a measurement of R>1 at LHCb would confirm the mechanism. For Kπ, the predicted branching ratios are ~1","pith_inferences":["Inference: A measurement of R^{ππ}_{χc1/χc0} significantly above ~1.5 would argue that non-factorizable or final-state-interaction effects differ for the two charmonia, pointing to the need to go beyond the simple factorizable-emission picture.","Inference: The predicted 25% constructive interference from ρ(1450) and ρ(1700) is testable by Dalitz-plot analysis; seeing destructive interference instead would require revising the GS-model weight factors.","Inference: Because the two normalization constants N_ρ and N_{K*} largely cancel in R_{K/π} but not in absolute rates, the ratio predictions are more robust than the individual branching ratios; experiments should prioritize measuring R_{K/π} and R_{χc1/χc0} before absolute branching fractions.","Inference: The same two-meson DA machinery should predict the corresponding B_s or B_c decays involving ω and φ resonances; if the fitted normalizations are truly universal, comparable signals should appear there, offering an independent cross-check."],"forward_implications":["If the predictions hold, B_c^+→χ_{c0}π^+π^0 and B_c^+→χ_{c1}π^+π^0 have branching ratios near 10^{-3}, comparable to B_c→J/ψπ, and should be accessible in the full LHCb dataset through χ_{cJ}→J/ψγ.","The ratio R^{ππ}_{χc1/χc0}≈1.30 (versus ≈4.7 for the two-body transition) gives a clean, systematic-cancelling test of whether resonant ρ production changes the relative P-wave charmonium yields.","The predicted f_L≈94% for χ_{c1} modes, decreasing to ~65% at ρ(1700), offers an angular-analysis discriminator for the two-meson DA model.","The Kπ channels, suppressed to ~10^{-6} with R_{K/π}≈2×10^{-3} independent of the charmonium spin, provide an SU(3)-breaking check.","Under the narrow-width approximation, the two-body estimates B(B_c→χ_{c0}ρ^+)≈3.24×10^{-3} and B(B_c→χ_{c1}ρ^+)≈4.19×10^{-3} disagree with the existing two-body PQCD numbers (2.97×10^{-4} and 1.40×10^{-3}), sharpening the question of why resonant production would be so much larger."],"fun_headline_variants":["χc1/χc0 ratio 1.3 from B_c ρ decay vs 4.7 two-body","B_c→χ_c(ρ→ππ): rates ~10⁻³, ratio near 1.3","Resonant ρ sets χc1/χc0 at 1.3, far from 4.7","PQCD: B_c→χ_cππ at 10⁻³, ratio ~1.3"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The calculation rests on the quasi-two-body factorization ansatz: all non-perturbative ππ/Kπ physics, including final-state interactions, is absorbed into two-meson distribution amplitudes with channel-independent form-factor normalizations; if those normalizations (or the F⊥/F∥≈f_T/f approximation) depend on the decaying meson or the charmonium partner, the predicted branching ratios shift, and the decisive ratio could move back toward the two-body value.","fun_headline_variants_meta":{"raw":{"variants":["χc1/χc0 ratio 1.3 from B_c ρ decay vs 4.7 two-body","B_c→χ_c(ρ→ππ): rates ~10⁻³, ratio near 1.3","Resonant ρ sets χc1/χc0 at 1.3, far from 4.7","PQCD: B_c→χ_cππ at 10⁻³, ratio ~1.3"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000338,"raw_usage":{"total_tokens":1821,"prompt_tokens":979,"completion_tokens":842,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":723,"completion_tokens_details":{"reasoning_tokens":739}},"tokens_in":723,"tokens_out":842,"duration_ms":8297,"temperature":1.0,"reasoning_tokens":739,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T03:25:21.703559+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure B(B_c^+→χ_{c0}π^+π^0) and B(B_c^+→χ_{c1}π^+π^0) at LHCb in the χ_{cJ}→J/ψγ chain. The central claim predicts their ratio R≈1.30; if the measured ratio instead approaches the two-body PQCD value ≈4.7, then the hypothesis that resonant ρ production reshapes the P-wave charmonium yields through factorizable emission is falsified.","supporting_citations":[],"review_version":1}