{"id":"f9cbeba4-e6ae-416d-9164-2af38ea8333b","arxiv_id":"2412.16068","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"high","formal_verification":"none","parameter_count":10,"one_line_summary":"QCD sum rule predictions for the tetraquarks bb anti-b anti-c and cc anti-c anti-b give masses near 15.70 and 9.68 GeV and widths of about 36 and 55 MeV.","lead":"This paper uses QCD sum rules to predict the masses and strong decay widths of two hypothetical fully heavy tetraquarks, one containing three b quarks and one c quark, the other three c quarks and one b quark. The predicted T_b state sits only 25 MeV above its decay threshold, making its width fragile, while T_c has two allowed decay modes and a total width near 55 MeV.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Width predictions rest on exponential form-factor extrapolation from Q^2>0 to Q^2=-m_M^2 without a stability check, and for T_b the ±95 MeV mass uncertainty straddles threshold, so the quoted errors understate both risks.","rationale":"The paper is a competent QCD sum-rule analysis: the two-point mass extraction satisfies the stated pole-dominance and OPE-convergence criteria, uncertainties in M^2 and s0 are given, and the comparison with quark-model masses is useful context. The genuinely load-bearing spot is the path from Euclidean form factors to on-shell couplings. This is not a rejection of standard practice; QCD sum-rule form factors are routinely extrapolated. The gap is that no sensitivity test is reported for the exponential ansatz in Eq. (36), even though the target point lies outside the fitted dimensionless range and, for T_b, at a negative Q^2 where the sum-rule machinery is not operative. The kinematic issue compounds this: the paper itself notes that the lower bound of the T_b mass is below the eta_b B_c^- threshold, so the strong-decay assertion is sensitive to the 95 MeV mass uncertainty, which is excluded from the width error in Eq. (41). Both problems are specific and addressable with the proposed checks, and neither by itself invalidates the mass predictions or the qualitative picture. I therefore agree with the reader's CONDITIONAL verdict and would not change it; the required revisions are to test the extrapolation and to propagate the T_b mass uncertainty into the width, or to soften the instability claim accordingly.","tokens_in":17068,"tokens_out":8073,"duration_ms":81176,"concrete_test":"Recompute G(-m_Bc^2), g1(-m_eta_c^2), and g2(-m_J/psi^2) from the same SR data using at least two alternative parametrizations, e.g., a 3-point Pade and a z-expansion, using the same Q^2 window; if any on-shell coupling moves outside the quoted central value plus its stated uncertainty, the widths are not robust. Independently, propagate the ±95 MeV uncertainty in m into Eq. (39): evaluate the T_b width at m=15603 MeV and m=15793 MeV, and preferably marginalize over a Gaussian mass distribution; if a substantial probability mass sits below threshold, the quoted ±10.4 MeV width error is incomplete.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The single most load-bearing weakness is the uncontrolled on-shell extrapolation of the strong form factors. In Sec. III, Eq. (36) fits G(Q^2) to three-point sum-rule outputs over Q^2=2-30 GeV^2 (dimensionless x=Q^2/m^2 in [0.008,0.122]) and then evaluates it at Q^2=-m_Bc^2≈-39.4 GeV^2, i.e., x≈-0.159, where the sum rule is not valid. The same procedure is used for both Tc couplings in Sec. IV. No check is reported on the choice of the two-exponent ansatz; a different but equally good fit form can move the on-shell point substantially because the fitted region does not constrain the behavior at negative x. Since all three decay widths depend directly on these couplings, the central quantitative claims (36 MeV and 54.7 MeV) inherit this uncontrolled extrapolation. A separate but reinforcing problem is kinematic: Sec. III states that the T_b mass exceeds the eta_b B_c^- threshold by only 25 MeV, while the mass uncertainty is ±95 MeV; Eq. (41) does not propagate this mass uncertainty, so the assertion that T_b is 'unstable' is not supported by the paper's own error budget.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript applies QCD sum rules to scalar tetraquarks T_b = bb\\bar b\\bar c and T_c = cc\\bar c\\bar b modeled as axial-vector diquark-antidiquark states. Two-point sum rules give masses m = (15698 ± 95) MeV and \\tilde m = (9680 ± 102) MeV with current couplings. Three-point sum rules are then used to compute strong form factors and couplings for T_b → η_b B_c^- and T_c → η_c B_c^+, J/ψ B_c^{*+}; these couplings are converted into partial widths Γ[T_b→η_b B_c^-] = 36.0 ± 10.4 MeV and Γ[T_c] = 54.7 ± 12.6 MeV, leading to the conclusion that both states are unstable and relatively narrow, with useful predictions for experiment.","tokens_in":17588,"tokens_out":10222,"duration_ms":88647,"significance":"If correct, this is the first quantitative sum-rule treatment of the widths of asymmetric fully-heavy tetraquarks, extending a series of papers by the same authors on fully-heavy exotics. The two-point analysis is standard and carefully controlled (PC ≥ 0.5, OPE convergence, clear Borel windows), and the mass predictions are compared with existing quark-model results. The three-point sum-rule calculations are also standard in structure, and the paper reports explicit Borel windows, PC values, and fit constants, which facilitate reproducibility. However, the width predictions are not yet on the same footing as the masses: the on-shell couplings are obtained by extrapolating fitted exponential form factors far outside the computed region, and for T_b the mass uncertainty straddles the decay threshold. These issues must be addressed before the width numbers can be regarded as reliable predictions.","major_comments":[{"comment":"The decay widths are controlled by the strong couplings G, g1, and g2, which are extracted by the same procedure: a fit function of the form Eq. (36) is fitted to three-point sum-rule results at Q^2 = 2–30 GeV^2 (for T_b, x = Q^2/m^2 ∈ [0.008, 0.122]) and then evaluated at the on-shell point Q^2 = -m_{B_c}^2 ≈ -39.4 GeV^2 (x ≈ -0.159). The on-shell point lies well outside the fitted region, and the paper does not test the sensitivity of the extracted coupling to the choice of fit ansatz; a different functional form that represents the sum-rule data equally well in the fitted region can produce a substantially different value at such a distant extrapolation. Because the widths in Eqs. (41), (50), and (62) depend on the squares of these couplings, this uncontrolled analytic continuation is load-bearing for the central width predictions. The same issue applies to g1 and g2 in Sec. IV, where the on-shell points Q^2 = -m_{η_c}^2 ≈ -8.9 GeV^2 and Q^2 = -m_{J/ψ}^2 ≈ -9.6 GeV^2 are outside the fitted region Q^2 = 2–20 GeV^2.","section":"Sec. III, Eq. (36); Sec. IV, Eqs. (46) and (58)"},{"comment":"The paper notes that the central T_b mass is only 25 MeV above the η_b B_c^- threshold and that the 1σ lower value 15603 MeV is below threshold, yet the error quoted for Γ[T_b] in Eq. (41) includes only the uncertainty in G and in the final-state meson masses. The ±95 MeV uncertainty in m is not propagated through the phase-space factor λ_0 in Eqs. (39)–(40); including it would produce a strongly asymmetric error and, within 1σ, a zero width. For T_c, the ±102 MeV mass uncertainty is likewise omitted from the partial widths in Eqs. (50) and (62) and from the full width in Eq. (63). The claim that these widths are useful search parameters is therefore not supported by the paper's own error budget until the mass dependence of the phase space is incorporated.","section":"Sec. III, Eq. (41); Sec. IV, Eqs. (50), (62), (63)"}],"minor_comments":[{"comment":"The text states that the form factors g1(Q^2) and g2(Q^2) are computed at 'Q^2 = 2 − 20 MeV^2'; the units should be GeV^2.","section":"Sec. IV A and IV B"},{"comment":"The sentence 'These current are introduces' should read 'These currents are introduced'.","section":"Eq. (52)"},{"comment":"The arXiv identifier is given as hep-ph/00101765; the correct number is hep-ph/0010175.","section":"Ref. [40]"},{"comment":"The fitted curves extend into the negative-Q^2 region where no sum-rule data exist; shading or a dashed line for the extrapolation region would make the extent of the extrapolation transparent.","section":"Figs. 4 and 5"},{"comment":"Section II establishes OPE convergence for the two-point correlators, but the three-point sum rules are truncated at the same dimension-4 level without an analogous numerical estimate of the neglected higher-dimensional contributions.","section":"Secs. III and IV"},{"comment":"The fit parameters for G(Q^2) in Eq. (37) and for F_1, F_2 in Sec. IV are quoted without uncertainties; reporting fit errors would at least provide a partial measure of the extrapolation uncertainty.","section":"Eq. (37) and Sec. IV"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the scope of the journal and follows the authors' established program. The two major issues are addressable in revision: a stability check of the fit ansatz and a proper propagation of the tetraquark mass uncertainty. I do not see a fundamental obstacle, but the width predictions should not be quoted as final until these points are resolved."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You asked about arXiv:2412.16068. Bottom line: this is a competent QCD sum-rule paper, and as far as I can tell from the cited literature it is the first to give quantitative fall-apart widths for the scalar b b \\bar b \\bar c and c c \\bar c \\bar b tetraquarks. The masses are in the range of existing quark-model predictions, the two-point sum rules are set up in the standard way with pole-dominance and OPE-convergence checks, and the paper is honest about the model dependence of the diquark currents. That last point matters: Section V explicitly says the physical states could be mixtures with the alternative sigma-mu-nu current, which is the right caveat.\n\nThe soft spots are where the reader's report and the stress-test put them. The three decay widths all come from couplings obtained by fitting an exponential form factor to sum-rule outputs for Q^2 between 2 and 30 GeV^2 and then evaluating it at Q^2 = -m^2 of the final meson, around -39 GeV^2 for the T_b case. That is a long way outside the fitted region, and the paper does not test whether a different but equally good fit form changes the on-shell coupling. Since the widths are the new results, this is the main thing a referee should probe.\n\nThe second issue is kinematic rather than technical. The T_b central mass is only 25 MeV above the eta_b B_c^- threshold, and the stated mass uncertainty is ±95 MeV. The paper actually acknowledges in Section III that at the lower end of that window the state would be stable against fall-apart. But the quoted width error, ±10.4 MeV, only includes the coupling uncertainty and the final-meson mass uncertainties; it does not propagate the ±95 MeV from m. So the claim that T_b is unstable is a central-value statement, not an error-budget statement. The paper says this openly, which I respect, but it does mean the headline width should be read with a grain of salt.\n\nThe citation pattern is fine. The paper builds on the authors' own earlier sum-rule work on fully heavy tetraquarks, and it cites the relevant quark-model and other predictions for b b \\bar b \\bar c and c c \\bar c \\bar b. I did not see missing references that would change the conclusions.\n\nWho gets value: people working on exotic hadron spectroscopy and on sum-rule methods. The masses are a useful data point, and the T_c width prediction is less precarious than the T_b one because the phase space is real. I would bring it to a reading group as an example of what can and cannot be claimed from form-factor extrapolations. It deserves a serious referee; the right refereeing question is whether the width calculation can be stabilized, not whether the topic is worth publishing.","headline":"Competent sum-rule paper, first quantitative fall-apart widths for these two tetraquarks, but the widths rest on an uncontrolled form-factor extrapolation and the T_b width error omits the dominant mass uncertainty.","tokens_in":18089,"tokens_out":2980,"would_cite":false,"duration_ms":26484,"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":"QCD sum rules predict two new tetraquarks, Tb and Tc, with masses 15.7 and 9.7 GeV and widths 36 and 55 MeV.","keywords":["fully heavy tetraquarks","scalar tetraquarks","QCD sum rules","diquark-antidiquark states","exotic hadrons","strong decay width","bottom-charm tetraquarks"],"falsifier":"Search the $\\eta_b B_c^-$ invariant-mass spectrum near 15.7 GeV and the $\\eta_c B_c^+$ and $J/\\psi B_c^{*+}$ spectra near 9.7 GeV: finding no peak at the predicted masses, or a peak whose width is far from $(36.0\\pm10.4)$ MeV for $T_b$ and $(54.7\\pm12.6)$ MeV for $T_c$, would refute the paper's central claim.","tokens_in":16824,"feed_emoji":"⚛️","tokens_out":12135,"duration_ms":96634,"temperature":0.7,"pith_summary":"This paper predicts two new fully heavy exotic mesons: a negatively charged scalar state $T_b$ with quark content $bb\\bar{b}\\bar{c}$ and mass $(15698\\pm95)$ MeV, and a positively charged scalar $T_c$ with content $cc\\bar{c}\\bar{b}$ and mass $(9680\\pm102)$ MeV. They are modeled as axial-vector diquark-antidiquark pairs and computed with QCD sum rules. Because both masses lie above the thresholds for decay into two ordinary heavy mesons, the paper argues they should be seen as narrow 'fall-apart' resonances: $T_b\\to\\eta_b B_c^-$ with width $(36.0\\pm10.4)$ MeV, and $T_c\\to\\eta_c B_c^+$ plus $J/\\psi B_c^{*+}$ with combined width $(54.7\\pm12.6)$ MeV. These numbers give experimental searches a concrete target.","feed_headline":"Tetraquarks predicted at 15.7 and 9.7 GeV with widths 36 and 55 MeV","feed_subtitle":"Fully heavy, charged states made of bottom and charm quarks should decay into ordinary meson pairs.","key_machinery":"The machinery is a pair of axial-vector diquark-antidiquark interpolating currents, e.g. $J(x)=b_a^T C\\gamma_\\mu b_b(x)\\,[\\bar b_a\\gamma^\\mu C\\bar c_b^T-\\bar b_b\\gamma^\\mu C\\bar c_a^T]$, with an analogous $c$-quark current for $T_c$. These currents pick out the color-antisymmetric $[\\mathbf{3}_c]\\otimes[\\mathbf{3}_c]$ diquark configuration. Two-point QCD sum rules for these currents yield the masses and current couplings, while three-point correlation functions involving the tetraquark current and two meson currents yield the strong form factors $G(q^2)$, $g_1(q^2)$, and $g_2(q^2)$. The widths follow from reading each form factor at its physical meson mass shell through an exponential fit, $G(Q^2)=G_0\\exp[c_1 Q^2/m^2+c_2(Q^2/m^2)^2]$, to the Euclidean sum-rule points.","core_discovery":"The paper's claim is that the scalar four-quark states $T_b=bb\\bar{b}\\bar{c}$ and $T_c=cc\\bar{c}\\bar{b}$ can be described as diquark-antidiquark systems built from axial-vector diquarks, and that QCD sum rules then determine their masses, current couplings, and strong decay widths. The extracted masses, $m=(15698\\pm95)$ MeV and $\\tilde{m}=(9680\\pm102)$ MeV, lie above the lowest two-meson thresholds, $\\eta_b B_c^-$ for $T_b$ and $\\eta_c B_c^+$ for $T_c$, so both states are predicted to be unstable under fall-apart decays. Computing three-point sum rules for the relevant vertices gives $\\Gamma[T_b\\to\\eta_b B_c^-]=(36.0\\pm10.4)$ MeV and a combined $\\Gamma[T_c]=(54.7\\pm12.6)$ MeV from the $\\eta_c B_c^+$ and $J/\\psi B_c^{*+}$ channels. On this picture $T_b$ is a narrow, charged, near-threshold state, while $T_c$ decays through two open channels with comparable partial widths.","pith_inferences":["The width calculation leans on a long exponential extrapolation, so an independent lattice QCD evaluation of the same strong couplings at the physical mass shells would be the most direct check of these numbers.","Because $T_b$ sits so close to threshold, a mass shift of order the quoted uncertainty could turn it from a narrow resonance into a stable particle; precision spectroscopy just above the $\\eta_b B_c^-$ threshold would test this boundary.","If future experiments resolve two nearby scalar states in the same channels, it would point to mixing between the $C\\gamma_\\mu\\otimes\\gamma^\\mu C$ and $C\\sigma_{\\mu\\nu}\\otimes\\sigma^{\\mu\\nu}C$ diquark configurations, a possibility the paper identifies but does not analyze.","The roughly 6000 MeV gap between $T_b$ and $T_c$ tracks the difference in heavy-quark masses; applying the same calculation to other flavor asymmetries might reveal a simple mass rule for all-heavy tetraquarks."],"forward_implications":["If the $T_b$ mass prediction is correct, the state sits only about 25 MeV above the $\\eta_b B_c^-$ threshold and should appear as a narrow peak that is separable from the threshold rise.","For $T_c$, both $\\eta_c B_c^+$ and $J/\\psi B_c^{*+}$ are open channels with partial widths near 27 and 28 MeV, so a search should cover both final states near 9.68 GeV.","The predicted masses fall within the spread of existing quark-model estimates, so a measurement would help decide which diquark-antidiquark interaction models are closest to QCD.","Since $T_b$ and $T_c$ are charged, negative and positive respectively, observing them would show that fully heavy tetraquarks are not restricted to neutral states.","A measured width consistent with these values would support the axial-vector diquark-antidiquark organization over alternative scalar-diquark or mixed currents, which the paper argues would give different parameters."],"supporting_citations":[{"why":"Introduces the QCD sum rule method used for the two-point mass and coupling extractions.","marker":"[32]"},{"why":"Supplies the companion sum rule formalism, including the gluon condensate contribution used in the OPE.","marker":"[33]"},{"why":"Provides the running quark masses and experimental meson masses used as numerical inputs.","marker":"[38]"},{"why":"Establishes the double dispersion relation used to derive the three-point sum rule for the strong form factor.","marker":"[39]"},{"why":"Provides the general three-point sum rule framework for form factor calculations that the decay-width part follows.","marker":"[40]"},{"why":"Demonstrates the same three-point sum rule machinery in prior applications.","marker":"[41]"},{"why":"Supplies the eta_b and B_c decay constants that enter the T_b to eta_b B_c^- width formula.","marker":"[42]"},{"why":"Gives the theoretical B_c^* mass used to set the T_c to J/psi B_c^{*+} phase space.","marker":"[44]"},{"why":"Supplies the B_c^* decay constant used in the J/psi B_c^{*+} channel.","marker":"[45]"},{"why":"Supplies the J/psi decay constant used in the same channel.","marker":"[46]"}],"fun_headline_variants":["Heavy tetraquarks: T_b and T_c decay into meson pairs","Bottom-charm tetraquarks found unstable: masses 15.7 and 9.7 GeV","Tetraquarks T_b and T_c: masses and decay widths predicted","QCD sum rules give masses and widths for heavy tetraquarks","Fully heavy tetraquarks with asymmetric quark content decay"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that each decay form factor, fitted to QCD sum rule points in the Euclidean region $Q^2=2$–$30$ GeV$^2$, can be extrapolated through an exponential function all the way to the timelike point $Q^2=-m_{B_c}^2\\approx -39.4$ GeV$^2$; if that continuation is wrong, the predicted widths change substantially.","fun_headline_variants_meta":{"raw":{"variants":["Heavy tetraquarks: T_b and T_c decay into meson pairs","Bottom-charm tetraquarks found unstable: masses 15.7 and 9.7 GeV","Tetraquarks T_b and T_c: masses and decay widths predicted","QCD sum rules give masses and widths for heavy tetraquarks","Fully heavy tetraquarks with asymmetric quark content decay"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001177,"raw_usage":{"total_tokens":4960,"prompt_tokens":1137,"completion_tokens":3823,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":753,"completion_tokens_details":{"reasoning_tokens":3722}},"tokens_in":753,"tokens_out":3823,"duration_ms":24823,"temperature":1.0,"reasoning_tokens":3722,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T10:49:45.282178+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Search the $\\eta_b B_c^-$ invariant-mass spectrum near 15.7 GeV and the $\\eta_c B_c^+$ and $J/\\psi B_c^{*+}$ spectra near 9.7 GeV: finding no peak at the predicted masses, or a peak whose width is far from $(36.0\\pm10.4)$ MeV for $T_b$ and $(54.7\\pm12.6)$ MeV for $T_c$, would refute the paper's central claim.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the companion sum rule formalism, including the gluon condensate contribution used in the OPE."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the double dispersion relation used to derive the three-point sum rule for the strong form factor."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the general three-point sum rule framework for form factor calculations that the decay-width part follows."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the theoretical B_c^* mass used to set the T_c to J/psi B_c^{*+} phase space."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the J/psi decay constant used in the same channel."}],"review_version":1}