{"id":"ba83851b-4cc5-4bbc-95c4-788d80fab385","arxiv_id":"2608.04693","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"Predictions for B_c decays to χ_{c1}(nP) using a covariant Bethe-Salpeter amplitude suggest the semileptonic μν mode is more promising than the π mode for probing the nature of χc1(3872).","lead":"This paper calculates how often B_c mesons decay into excited charmonium states, assuming the known χc1(3872) particle is an ordinary 2P charmonium state. It argues the semileptonic decay B_c to χc1(3872) muon neutrino is the best channel for LHCb to look for, requiring about twice the current data, while the pion channel needs about twenty times more.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 20x/2x luminosity conclusions in Sec. V.C rest on a single model's branching ratio with only parametric errors, despite Table II showing a factor ~40 model spread for the same channel.","rationale":"The reader's conditional verdict is appropriate, but their stated weakest assumption (the charmonium assignment of chi_c1(3872)) is an explicitly declared premise rather than a hidden flaw. A more load-bearing problem is that the paper's headline experimental conclusions are not robust to the model spread it itself reports. The check I propose is cheap: it uses numbers already in Tables II and III to test whether the 20x/2x statements survive. I do not dispute the internal consistency of the Salpeter calculation or the value of an improved covariant amplitude; the issue is the claim of 'demonstrating' insufficient luminosity. If the sensitivity check shows the conclusions are stable, the paper stands. If not, the abstract and conclusions should be softened to report a model-dependent prediction. Verdict remains CONDITIONAL.","tokens_in":19923,"tokens_out":8809,"duration_ms":97554,"concrete_test":"Recompute R_{chi_c1(3872)/psi(2S)} and the implied required B_c statistics after (i) replacing B(B_c^+ -> chi_c1(3872) pi^+) by the values implied by Refs. [74] and [76] (width ratios 1.4/0.81 and 32.3/0.81 applied to the Table III branching ratio) and (ii) recomputing B(B_c^+ -> psi(2S) pi^+) with the new MB amplitude of Eq. (18) instead of the MA value from Ref. [108]. If the required luminosity changes by more than a factor of about two in either case, or if R crosses the LHCb bound 0.05, the 20x/2x statements should be replaced by a model-dependent range.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central experimental statements in Sec. V.C — that about 20 times the current B_c sample is needed for B_c^+ -> chi_c1(3872) pi^+ and only twice for B_c^+ -> chi_c1(3872) mu^+ nu_mu — follow directly from R = 0.0024, which uses B(B_c^+ -> chi_c1(3872) pi^+) = 6.31 x 10^-6 (Table III). The uncertainty quoted on this branching ratio is only the +/-5% parametric variation of quark masses and potential-model parameters; it does not cover the model dependence documented in the same paper. In Table II, B_c^+ -> chi_c1(3872) pi^+ widths are 0.81 (this work), 1.4 (Ref. [74]), 24 (Ref. [107]), and 32.3 (Ref. [76]) in units of 10^-17 GeV, a factor of about 40 between the smallest and largest. Adopting the upper end would move the predicted R above the LHCb upper limit, so the 'insufficient luminosity' explanation would cease to hold. Furthermore, R combines a numerator computed with the new amplitude MB (Eq. 18) with a denominator B(B_c^+ -> psi(2S) pi^+) = 0.0266% taken from the old MA amplitude (Ref. [108]); this mixes two versions of the same model. The feasibility claims therefore need a sensitivity study over existing model predictions and an MB-consistent denominator before they can be stated as definite.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper calculates semileptonic and color-favored nonleptonic B_c^+ decays to chi_c1(nP) (n=1,2,3) in the Bethe-Salpeter/Salpeter formalism, treating chi_c1(3872) as the conventional chi_c1(2P) charmonium state. The authors introduce a new covariant hadronic amplitude M_B that evaluates the final-state wave function in its rest frame, replacing the older M_A amplitude in large-recoil processes. Using this amplitude they compute decay widths and branching fractions, predict the ratio R_{chi_c1(3872)/psi(2S)} = 0.0024, and conclude that about 20 times the current B_c sample is needed to observe B_c^+ -> chi_c1(3872)pi^+, while B_c^+ -> chi_c1(3872)mu^+nu_mu could be accessible with roughly twice the existing data.","tokens_in":20300,"tokens_out":7869,"duration_ms":85883,"significance":"The methodological upgrade from M_A to M_B is a genuine step forward for large-recoil transitions within this model class, and the semileptonic channel B_c^+ -> chi_c1(3872)mu^+nu_mu as a more promising probe is a concrete, falsifiable suggestion for LHCb. The prediction R = 0.0024 is consistent with the existing upper limit and provides a useful target. However, the quantitative luminosity claims are based on a single model's branching fraction with only parametric error bars, while Table II shows a factor of about 40 spread among published model predictions for the same channel. The paper's central conclusions are therefore conditional on model assumptions that are acknowledged but not propagated into the stated uncertainties.","major_comments":[{"comment":"The central phenomenological statement of the paper, R_{chi_c1(3872)/psi(2S)} = 0.0024 and the inference that roughly 20 times the current B_c sample is needed, rests on B(B_c^+ -> chi_c1(3872)pi^+) = (6.31^{+0.97}_{-1.24}) x 10^-6, whose quoted uncertainty is generated only by the +/-5% scan of quark masses and potential parameters described in Section V.A. Table II, however, documents a factor of about 40 spread for the same channel among published calculations: 0.81 x 10^-17 GeV in this work, 1.4 x 10^-17 GeV in Ref. [74], 24 x 10^-17 GeV in Ref. [107], and 32.3 x 10^-17 GeV in Ref. [76]. If a prediction in the upper part of this range were adopted, the corresponding R would exceed the LHCb upper limit of 0.05, and the 'insufficient luminosity' explanation would cease to hold. The feasibility claim therefore needs a sensitivity study over the existing model predictions, or at least an explicit statement that it is conditional on the authors' model.","section":"Section V.C, Tables II and III"},{"comment":"The predicted ratio R combines a numerator computed with the new amplitude M_B (Eq. 18) with a denominator B(B_c^+ -> psi(2S)pi^+) = 0.0266% taken from Ref. [108], which is based on the older amplitude M_A (Eq. 15). The two amplitudes are not equivalent at the recoil points relevant here: Table I shows M_B versus M_A differences of 12.9 vs 15.2 (x 10^-16 GeV) for B_c^+ -> chi_c1(1P)e^+nu_e and 0.96 vs 1.53 for B_c^+ -> chi_c1(3872)e^+nu_e. Recomputing the psi(2S) branching ratio with M_B would change R by an amount that is not negligible relative to the precision with which R = 0.0024 is quoted. I ask the authors to supply an M_B-consistent denominator or to estimate the sensitivity of R to this choice.","section":"Section V.C, Eqs. (15) and (18), Ref. [108]"},{"comment":"The derivation of M_B proceeds by inserting positive-energy projectors for the final-state quarks and neglecting the negative-energy contributions in the third equality of Eq. (16), together with the spectator condition p'_1 = p_1. For chi_c1(2P), the paper itself reports v^2 = 0.39 (Section I), so the negative-energy components are not obviously small. No estimate of the truncation error is given. Because M_B is the basis of every numerical result in Tables I-IV, the authors should either quantify this error, for example by retaining the negative-energy terms in the residue integration or by comparing with a different spectator prescription, or state this as a systematic limitation in the error budget.","section":"Section III, Eqs. (16)-(18)"}],"minor_comments":[{"comment":"Reference [53] is labeled 'A. M. Sirunyan et al. (CDF Collaboration)'; Sirunyan is the CMS collaboration, so the collaboration label appears to be wrong. Please verify the citation.","section":"References"},{"comment":"Table III lists 'B_c^+ -> chi_c1(1P)tau^+tau_l'; this should presumably be tau^+nu_tau, since the decay is semileptonic.","section":"Table III"},{"comment":"The Introduction contains the typo 'corraboration' for 'collaboration'.","section":"Introduction"},{"comment":"The abbreviation 'epsilon*_mu epsilon PP_f = epsilon*_{mu alpha beta gamma} epsilon_alpha P_beta P_gamma^f' has an index mismatch: epsilon*_mu has one free Lorentz index while the right-hand side involves an additional polarization vector epsilon_alpha. Please re-express the K_4 term with consistent indices.","section":"Section II, Eq. (3)"},{"comment":"The ratio in Eq. (24), R_{chi_c1(nP)}, is a lepton-universality ratio for tau versus electron final states, but the same symbol R is used for the LHCb ratio R_{chi_c1(3872)/psi(2S)} in Eq. (1). Using different notation for these two quantities would reduce confusion.","section":"Section V, Eq. (24) and Table I"}],"recommendation":"major_revision","confidential_remarks":"The paper's real contribution is the amplitude upgrade for large-recoil transitions and the suggestion that the semileptonic channel is the more promising search mode. The headline luminosity ratios, however, are too dependent on one model to be stated as definite predictions; I would require the sensitivity analysis described in the major comments before publication. The reference list should also be checked for collaboration mislabeling beyond Ref. [53]."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper is a legitimate extension of the authors' Bethe-Salpeter program. The new covariant amplitude M^B_mu evaluates the final-state wave function in its own rest frame, a genuine improvement for large-recoil transitions, and it yields semileptonic widths for chi_c1(1P) that agree well with several independent quark models. The new chi_c1(3P) predictions and the explicit treatment of the 2P nodal suppression are useful. The paper is also transparent about its key assumption: chi_c1(3872) is taken as conventional charmonium, and the model mass is moved to the physical mass by tuning V0. That is a defensible strategy, and the resulting R = 0.0024 is consistent with the LHCb upper limit.\n\nThe soft spots are concentrated in the experimental feasibility claims. The paper quotes uncertainties of only a few percent on B(B_c -> chi_c1(3872) pi^+) from parametric variation, but Table II shows a factor of roughly 40 spread for this same width across published models. Adopting the largest prediction would make R exceed the LHCb bound, so the statement that the non-observation is simply a luminosity problem is not robust. The ratio R also mixes the new M^B amplitude for the numerator with an older M^A branching fraction for B_c -> psi(2S)pi from Ref. [108]; that is an internal inconsistency that should be fixed or justified. Finally, the 2x luminosity estimate for the semileptonic channel is a simple scaling argument and needs a real sensitivity study to be credible.\n\nNone of this undermines the core calculation. The amplitude derivation is plausible, the comparisons are honest, and the paper is careful to say it is probing the charmonium hypothesis rather than settling the nature of X(3872). The citation pattern is appropriate, and the self-citations are to the group's own prior work, which is relevant here. It is a paper for hadron phenomenologists and LHCb analysts, and it deserves a serious referee. My recommendation: send to peer review, and in the report push the authors to add a model-uncertainty band, use a consistent amplitude for both numerator and denominator, and temper or qualify the luminosity claims.","headline":"Solid Salpeter update for Bc to chi_c1(nP) with an honest charmonium assumption, but the 20x and 2x luminosity claims outrun the model uncertainty the paper itself documents.","tokens_in":20886,"tokens_out":2047,"would_cite":true,"duration_ms":26021,"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":"The paper predicts that $B_c^+\\to\\chi_{c1}(3872)\\mu^+\\nu_\\mu$ is the most accessible weak-decay probe of $\\chi_{c1}(3872)$'s structure.","keywords":["B_c meson weak decays","χc1(3872)","charmonium 2P state","Bethe–Salpeter equation","Salpeter wave functions","semileptonic decay","nonleptonic decay","relativistic corrections"],"falsifier":"Use twice the current $B_c$ sample to search for $B_c^+\\to\\chi_{c1}(3872)\\mu^+\\nu_\\mu$ and, if it is found, measure its branching fraction; a rate far below the predicted $\\sim7.4\\times10^{-5}$, or a clear absence at twenty times the sample for the $\\pi^+$ mode, would break the $\\chi_{c1}(2P)$ assumption that carries the calculation.","tokens_in":19725,"feed_emoji":"🔭","tokens_out":7699,"duration_ms":74841,"temperature":0.7,"pith_summary":"This paper argues that $\\chi_{c1}(3872)$ remains a viable conventional charmonium state, the $\\chi_{c1}(2P)$, and that weak decays of the $B_c$ meson provide a sharp way to test that idea. The authors compute $B_c^+\\to\\chi_{c1}(nP)\\ell^+\\nu_\\ell$ and $B_c^+\\to\\chi_{c1}(nP)X$ in the Bethe–Salpeter formalism, using an upgraded transition amplitude that evaluates the final-state wave function in its rest frame. Their central prediction is $R_{\\chi_{c1}(3872)/\\psi(2S)}=0.0024$, consistent with the experimental upper limit, which they read as a luminosity shortfall rather than evidence against the charmonium picture: roughly twenty times the current $B_c$ sample would be needed for the $\\pi^+$ mode. By contrast, they predict $B_c^+\\to\\chi_{c1}(3872)\\mu^+\\nu_\\mu$ should be observable with about twice the existing data, making it the most practical channel to probe the state's structure.","feed_headline":"Muon decay could reveal what χc1(3872) is made of","feed_subtitle":"A new calculation says the muon channel needs twice the current data; the pion channel needs twenty times.","key_machinery":"The load-bearing object is the covariant hadronic transition amplitude $M^B_\\mu$ (Eq. 18), which expresses $\\langle\\chi_{c1}|J_\\mu|B_c\\rangle$ as an overlap of the $B_c$ Salpeter wave function and the final-state Salpeter wave function evaluated in its own rest frame, with a factor $(M_f-\\tilde{\\omega}_1-\\tilde{\\omega}_2)/(E_f-\\omega'_1-\\omega'_2)$ from the quark propagators. This replaces the older amplitude $M^A_\\mu$, which parameterized the final-state momentum in the initial-state rest frame and handled large recoil poorly. The wave functions are constructed by $J^P$ rather than by $^{2S+1}L_J$, and the Cornell potential with screening is used; the radial wave functions are obtained by solving the Salpeter equation.","core_discovery":"Under the assignment $\\chi_{c1}(3872)=\\chi_{c1}(2P)$, the paper claims that all $B_c\\to\\chi_{c1}(nP)$ widths shrink sharply with $n$ because the 2P and 3P wave functions have radial nodes that cancel much of the overlap integral. The semileptonic width $B_c^+\\to\\chi_{c1}(3872)e^+\\nu_e$ is predicted as $(0.96^{+0.61}_{-0.42})\\times10^{-16}\\ \\mathrm{GeV}$, about an order of magnitude below the 1P channel, and many nonleptonic channels are even smaller. Using its own $B_c\\to\\psi(2S)\\pi^+$ branching fraction for normalization, the paper obtains $R_{\\chi_{c1}(3872)/\\psi(2S)}=0.0024$ and infers that the observed non-observation of $B_c^+\\to\\chi_{c1}(3872)\\pi^+$ is expected from statistics, not from an exotic structure. The semileptonic mode is singled out as the decisive channel, since $R(\\pi^+)=0.085$ makes the pion mode about twelve times harder to see.","pith_inferences":["Extension beyond the paper: repeating the same $M^B_\\mu$ calculation with $\\chi_{c1}(3872)$ wave functions from molecular or tetraquark models would test whether the ratio $R(X)$ is more discriminating than any single width.","Extension beyond the paper: the mass-tuning step is the fragile point, so a sensitivity study of $R_{\\chi_{c1}(3872)/\\psi(2S)}$ to $V_0$ within the model's own uncertainty would sharpen the twenty-times luminosity estimate.","Extension beyond the paper: if the semileptonic rate is confirmed near prediction, it would indirectly validate the large relativistic corrections, which enter at the 43% level for the 2P state."],"forward_implications":["The $\\pi^+$ channel is not a failure of the charmonium picture; observing it requires roughly twenty times the current $B_c$ data.","The semileptonic muon channel should be the experimental priority, needing only about twice the current data.","The ratio $R_{\\chi_{c1}(3872)/\\psi(2S)}=0.0024$ gives a concrete normalization target for future searches.","Nonleptonic channels with $\\rho^+$ and $D_s^{*+}$ have larger branching fractions but are harder to access experimentally, so the semileptonic mode remains the clean probe.","A measured value of $R(\\pi^+)$ near 0.085 would support the conventional charmonium assignment."],"supporting_citations":[{"why":"Supplies the experimental upper limit $R_{\\chi_{c1}(3872)/\\psi(2S)}<0.05$ that the paper's central ratio must be consistent with.","marker":"[55]"},{"why":"Provides the earlier Bethe–Salpeter calculation of $B_c\\to\\chi_{c1}(nP)$ widths using the old amplitude formula $M^A_\\mu$, which the new formula is designed to improve.","marker":"[57]"},{"why":"Defines the Salpeter equation, the instantaneous approximation of the Bethe–Salpeter equation used to solve the meson wave functions.","marker":"[79]"},{"why":"Gives the $J^P$-based Salpeter wave function representation for the pseudoscalar $B_c$ meson.","marker":"[88]"},{"why":"Gives the $J^P$-based Salpeter wave function representation for the $1^{++}$ state, including the relativistic $F_3$ term.","marker":"[90]"},{"why":"Justifies constructing wave functions by $J^P$ rather than $^{2S+1}L_J$, the basis for the relativistic wave functions used here.","marker":"[92]"},{"why":"Supplies the approach of evaluating the final-state wave function in its rest frame, the key upgrade leading to $M^B_\\mu$.","marker":"[95]"},{"why":"Provides the $B_c\\to\\psi(2S)\\pi^+$ branching fraction used as the normalization in the predicted ratio $R_{\\chi_{c1}(3872)/\\psi(2S)}$.","marker":"[108]"},{"why":"Motivates the use of nonleptonic-to-semileptonic ratios like $R(X)$ as relatively universal and reliably predicted quantities.","marker":"[58]"}],"fun_headline_variants":["Semileptonic B_c decay: 2x data to spot χc1(3872)","Pion limit on χc1(3872) just lacks data, muon mode shines","χc1(3872) as 2P charmonium: B_c semileptonic test","Twice data for muon, 20x for pion: χc1(3872) from B_c","Muon channel 10x more sensitive for χc1(3872) than pion"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation assumes that $\\chi_{c1}(3872)$ really is the ordinary $\\chi_{c1}(2P)$ charmonium state and that moving the model's predicted mass from 3928.7 MeV to the observed 3871.6 MeV by tuning the free potential parameter $V_0$ yields a trustworthy wave function; if the state is a tetraquark, molecule, hybrid, or mixture, the predicted branching fractions do not describe it.","fun_headline_variants_meta":{"raw":{"variants":["Semileptonic B_c decay: 2x data to spot χc1(3872)","Pion limit on χc1(3872) just lacks data, muon mode shines","χc1(3872) as 2P charmonium: B_c semileptonic test","Twice data for muon, 20x for pion: χc1(3872) from B_c","Muon channel 10x more sensitive for χc1(3872) than pion"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001143,"raw_usage":{"total_tokens":4826,"prompt_tokens":1113,"completion_tokens":3713,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":729,"completion_tokens_details":{"reasoning_tokens":3585}},"tokens_in":729,"tokens_out":3713,"duration_ms":32057,"temperature":1.0,"reasoning_tokens":3585,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T18:46:54.695606+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Use twice the current $B_c$ sample to search for $B_c^+\\to\\chi_{c1}(3872)\\mu^+\\nu_\\mu$ and, if it is found, measure its branching fraction; a rate far below the predicted $\\sim7.4\\times10^{-5}$, or a clear absence at twenty times the sample for the $\\pi^+$ mode, would break the $\\chi_{c1}(2P)$ assumption that carries the calculation.","supporting_citations":[{"cited_title":"Ablikim et al.(BESIII Collaboration), Phys","cited_arxiv_id":null,"evidence_quote":"Supplies the experimental upper limit $R_{\\chi_{c1}(3872)/\\psi(2S)}<0.05$ that the paper's central ratio must be consistent with."},{"cited_title":"Aaij et al","cited_arxiv_id":null,"evidence_quote":"Provides the earlier Bethe–Salpeter calculation of $B_c\\to\\chi_{c1}(nP)$ widths using the old amplitude formula $M^A_\\mu$, which the new formula is designed to improve."},{"cited_title":"Zhang, Z.-J","cited_arxiv_id":null,"evidence_quote":"Defines the Salpeter equation, the instantaneous approximation of the Bethe–Salpeter equation used to solve the meson wave functions."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the $J^P$-based Salpeter wave function representation for the pseudoscalar $B_c$ meson."},{"cited_title":"Loringa, K","cited_arxiv_id":null,"evidence_quote":"Gives the $J^P$-based Salpeter wave function representation for the $1^{++}$ state, including the relativistic $F_3$ term."},{"cited_title":"Wang, Phys","cited_arxiv_id":null,"evidence_quote":"Justifies constructing wave functions by $J^P$ rather than $^{2S+1}L_J$, the basis for the relativistic wave functions used here."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the approach of evaluating the final-state wave function in its rest frame, the key upgrade leading to $M^B_\\mu$."},{"cited_title":"Navaset et al","cited_arxiv_id":null,"evidence_quote":"Provides the $B_c\\to\\psi(2S)\\pi^+$ branching fraction used as the normalization in the predicted ratio $R_{\\chi_{c1}(3872)/\\psi(2S)}$."},{"cited_title":"Aaij et al","cited_arxiv_id":null,"evidence_quote":"Motivates the use of nonleptonic-to-semileptonic ratios like $R(X)$ as relatively universal and reliably predicted quantities."}],"review_version":1}