{"id":"b8e53eb6-6851-43a8-9735-a57be45a5f9d","arxiv_id":"2505.00078","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Lattice QCD with bottom quarks suggests a stable doubly bottom tetraquark (Tbb) bound by about 116 MeV relative to the B B* threshold.","lead":"This paper reports a lattice QCD calculation of a doubly heavy tetraquark made of two bottom quarks and two light antiquarks. It finds evidence for a bound state about 116 MeV below the B-B* meson pair threshold. A stable tetraquark with this quantum number is a long-sought exotic hadron, and a precise lattice prediction would guide experimental searches.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 116 MeV binding energy rests on a scattering-length-only p cot δ0 fit with no effective-range term; the omitted r_e p^2 term is comparable to the pole term at the claimed binding momentum, so the central number is unprotected.","rationale":"The reader's conditional verdict is well aligned with my reading. The strongest claim is the binding energy, but the amplitude analysis in Section 3 uses the lowest-order effective-range expansion p cot δ0 = A[0] + A[1]·a, with no p^2 (effective-range) term. The pole momentum κ ≈ 0.79 GeV is not small compared with the expected hadronic scale, so a typical effective range of order 0.5 fm contributes a term comparable to 1/a0 and can shift the pole substantially. The paper gives no estimate of r_e and no evidence that the data constrain it. I therefore agree with the reader's weakest assumption. The elastic one-channel truncation is less problematic here because the claimed state lies below B*B*, but the missing momentum-dependent term is the same class of unquantified systematic. The qualitative signal of attraction is plausible and consistent with prior work, so the result should remain conditional rather than rejected. However, the quantitative 116 MeV binding energy should not be used as a definitive lattice determination until the effective-range sensitivity is checked. The companion-paper dependence and the explicit deferral of left-hand cut effects are further reasons that the current manuscript alone does not fully support the headline number.","tokens_in":4070,"tokens_out":5363,"duration_ms":60635,"concrete_test":"Re-fit the reported finite-volume spectra with the amplitude p cot δ0 = -1/a0 + (r_e/2)p^2 + c_1 a, using the same ensembles and energy levels, and report the resulting r_e and uncertainty. If r_e is consistent with zero and r_e κ/2 ≪ 1, the original binding energy is supported. If r_e ≈ 0.5 fm or larger, recompute the pole position; a shift of order 100 MeV or movement of the pole above threshold would demonstrate that the scattering-length-only extraction is not robust. If the data cannot constrain r_e, that itself settles the concern: the 116 MeV value is an unverified truncation rather than a lattice-QCD determination.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim is the bound-state pole at 116 MeV below the BB* threshold, obtained in Section 3 by fitting p cot δ0 = A[0] + A[1]·a and extrapolating to the physical point. At the claimed binding energy, the pole momentum is κ = sqrt(2μE_b) ≈ 0.79 GeV (with μ ≈ 2.65 GeV). In the effective-range expansion p cot δ0 = -1/a0 + (r_e/2)p^2, the fitted scattering-length term gives 1/a0 ≈ 0.79 GeV, while the omitted effective-range term is (r_e/2)κ^2. This omitted term is comparable to 1/a0 unless r_e ≪ 2/κ ≈ 2.5 GeV^{-1} ≈ 0.5 fm, which is not established. Typical hadronic effective ranges are of this order, so the extracted binding energy is highly sensitive to the truncation. The paper presents no estimate of r_e and no second independent constraint on the amplitude shape; the reported fits use ground-state levels only. The elastic one-channel assumption is somewhat safer because the claimed state lies below B*B* and BBπ thresholds, but the same unquantified momentum dependence is the load-bearing issue. The text itself flags that left-hand cuts are 'beyond the scope of this work', and the essential ensemble and fitting details are deferred to Ref. [7]. These self-admitted omissions do not invalidate the qualitative attraction signal but do leave the quantitative 116 MeV result model-dependent.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports a lattice QCD study of the I(JP)=0(1+) bb\\bar u\\bar d system, extracting the finite-volume spectrum on four MILC ensembles with two volumes and three lattice spacings, and performing a Lüscher-based amplitude analysis of BB* scattering. The observed negative energy shift relative to the BB* threshold and the extracted scattering length are interpreted as evidence for a bound tetraquark Tbb with binding energy -116(+30/-36) MeV. The paper builds on the authors' earlier work and on a companion paper Ref. [7], using overlap valence quarks, an NRQCD bottom quark, wall-source and box-sink smearing, and a GEVP analysis. The central quantitative claim is the deep binding energy obtained from a scattering-length-only parametrization of p cot δ0.","tokens_in":4419,"tokens_out":3859,"duration_ms":42008,"significance":"If the result holds, it would be an important quantitative confirmation of a stable doubly bottom tetraquark, a system that is central to current discussions of exotic hadrons and of the heavy-quark spin structure of QCD. The analysis uses standard and well-tested finite-volume techniques, multiple lattice volumes and spacings, and two chiral extrapolation forms, which gives credibility to the qualitative attraction signal. However, the quantitative binding energy rests on an unverified truncation of the effective-range expansion, and the manuscript defers essential fitting and ensemble details to Ref. [7]. The paper therefore currently supports the existence of a bound state at the qualitative level, but the 116 MeV number is not yet a robust lattice QCD determination.","major_comments":[{"comment":"The central result, a0_phys = 0.25(4?3) fm and binding energy -116(+30/-36) MeV, is obtained from a fit of p cot δ0 to f = A[0] + A[1]·a, i.e. a scattering-length-only model with no effective-range term and no momentum dependence. At the claimed bound-state pole, |p| ≈ 0.79 GeV, the omitted term (r_e/2)p^2 is comparable to 1/a0 unless r_e is far below 0.5 fm, a condition that is not established anywhere in the paper. Without an estimate of r_e from the data, for example by including a p^2 term constrained by the excited finite-volume levels, or a robust bound on its size, the 116 MeV binding energy is a direct consequence of the truncation rather than a robust lattice QCD result. This is load-bearing for the central claim.","section":"Section 3, 'Amplitude analysis'"},{"comment":"The manuscript defers all ensemble and fitting details to Ref. [7]: the operator sets, GEVP time ranges, individual energy levels, and the inputs to the continuum and chiral extrapolations are not shown or summarized. In a standalone paper claiming a quantitative binding energy, these details are necessary to assess whether the quoted errors include the dominant systematic uncertainties, particularly the continuum extrapolation and the chiral extrapolation. Please include a table of energy levels or a detailed summary of the fits and their ranges.","section":"Section 3, Fig. 1 and final paragraph"},{"comment":"The figure shows two chiral/continuum extrapolation ansätze, f1 = c0 + c1 Mps and f2 = c0 + c1 Mps^2, but the manuscript does not state which one is used for the central value nor the spread between the two results. If the central binding energy changes appreciably between the two forms, the quoted ±30/-36 MeV uncertainty underestimates the systematic error. The paper should report the central values and uncertainties for both fits and justify the choice of the central value.","section":"Section 3, final paragraph"}],"minor_comments":[{"comment":"There are several typographical errors, including 'within of these exotic states' in the Abstract and 'exitence' in the Conclusion; these should be corrected.","section":"Abstract and Conclusion"},{"comment":"The notation in Eq. (1) is unclear: the tilde on Φ† is not defined, and the meaning of the extra space in Φi(x,t ) is confusing. Please define the operators and their normalization explicitly.","section":"Section 3, Eq. (1)"},{"comment":"The reported value '0.25(4?3) fm' is ambiguous; the asymmetric error convention should be stated explicitly, for example as 0.25(+0.04/-0.03) fm or with the two uncertainties defined.","section":"Section 3, final paragraph"},{"comment":"The marker and color conventions are referred to Ref. [7] rather than described; since the figure is central to the argument, the caption should be self-contained enough for the reader to identify the ensembles and the different pseudoscalar masses.","section":"Figure 1 caption"}],"recommendation":"major_revision","confidential_remarks":"The manuscript appears to be a proceedings-style contribution that relies heavily on a companion paper (Ref. [7]) for technical details. For a standalone journal publication, the absence of those details is a significant obstacle, but the qualitative result is interesting. The main technical risk is the scattering-length-only amplitude model; if the authors can add an effective-range term or otherwise quantify the truncation error, the paper would be considerably stronger. I see no citation issues or invented entities; the methods are standard in the field."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a legitimate lattice calculation, and the claim that the b-b-bar-u-bar-d system has an attractive BB* interaction and likely a bound state is reasonable. The specific number, -116 MeV, is not as solid as the paper suggests. It comes from a scattering-length-only amplitude fit, and at the claimed pole momentum the omitted effective-range term is as large as the term that is kept unless r_e is unusually small. The paper never estimates r_e.\n\nWhat is genuinely new: multiple volumes, three lattice spacings, box-sink smearing, and a Luescher amplitude analysis, giving a scattering length a0 = 0.25 fm and a binding energy from the pole. That is an incremental but real step beyond the authors' 2019 paper and Meinel et al. The methodology is standard, the ensembles are public MILC ones, and the use of ratio correlators is sensible. I don't see circularity: the binding energy is derived from the finite-volume spectrum, not put in by hand. The reference list looks appropriate, and the reliance on companion paper [7] for ensemble details is fine if that paper actually contains them.\n\nThe soft spot is exactly the amplitude truncation. At -116 MeV, the pole momentum kappa is about 0.79 GeV, so 1/a0 is about 0.79 GeV. The omitted (r_e/2)kappa^2 term is comparable to that unless r_e is well below about 0.5 fm. Typical hadronic effective ranges are around that size. The elastic one-channel assumption is likely safe because the state sits below B*B* and BB-pi thresholds, and the paper says left-hand cuts are beyond scope, but that doesn't fix the missing momentum dependence. The quoted errors, +30/-36 MeV, are fit errors; they do not include the effective-range uncertainty. The paper is also too terse: essential details of the GEVP fits and the continuum/chiral extrapolation live in the companion paper, so a referee cannot check the analysis from this text alone.\n\nI'd send this to peer review, but with a clear request: estimate or bound the effective range and show the sensitivity of the binding energy to it, or soften the central claim accordingly. The result is important enough to deserve referee time; it just shouldn't be published as a first-principles prediction of 116 MeV without that caveat. For a reader, this is a useful status report and also a compact cautionary example of why scattering-length-only pole extraction is dangerous at large binding momentum.","headline":"The qualitative attraction is probably real, but the 116 MeV binding energy rests on an unquantified effective-range truncation and should not be taken at face value.","tokens_in":4943,"tokens_out":4995,"would_cite":true,"duration_ms":49811,"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":"Lattice QCD finds the doubly bottom tetraquark is bound","keywords":["lattice QCD","tetraquarks","doubly heavy tetraquarks","bottom quarks","exotic hadrons","BB* scattering","finite-volume spectrum","finite-volume quantization"],"falsifier":"Compute the effective range from the same finite-volume levels, or add the first inelastic channels ($B^*B^*$ and $BB\\pi$) to the quantization condition: if the pole satisfying $p\\cot\\delta_0=+\\sqrt{-p^2}$ no longer sits near 116 MeV below threshold, or if the extracted scattering length moves by more than the quoted errors, the central claim is falsified.","tokens_in":3897,"feed_emoji":"⚛️","tokens_out":16579,"duration_ms":156931,"temperature":0.7,"pith_summary":"This paper tries to establish that the $I(J^P)=0(1^+)$ $bb\\bar u\\bar d$ system, two bottom quarks together with an up-antidown pair, forms a stable tetraquark rather than an unbound $BB^*$ pair. The authors compute the finite-volume spectrum on multiple lattice ensembles and use a standard quantization condition to convert the lowest energy level into an S-wave scattering amplitude. After continuum and chiral extrapolations, the amplitude has a pole below the $BB^*$ threshold, corresponding to a binding energy of about $116^{+30}_{-36}$ MeV and a scattering length $a_0^{\\mathrm{phys}}=0.25(4 3)$ fm. A stable doubly bottom tetraquark would be a concrete four-quark state predicted by QCD and a target for future experiments.","feed_headline":"Lattice QCD finds the doubly bottom tetraquark is bound","feed_subtitle":"Lattice QCD reports a 116 MeV binding for the b b ū d system, giving experiment a concrete state to look for.","key_machinery":"The load-bearing object is the finite-volume spectrum together with its translation into an infinite-volume amplitude. Correlation matrices built from diquark-antidiquark and $BB^*$ meson-meson interpolating operators in the $T_{1g}$ irrep, the finite-volume counterpart of $J^P=1^+$, are analysed with a generalized eigenvalue problem (GEVP), and wall-source with box-sink smearing suppresses excited-state contamination. The extracted ground-state energy shifts relative to the $BB^*$ threshold are fed into the standard finite-volume quantization condition, $p\\cot\\delta_0(p) = 2Z_{00}(1;(pL/2\\pi)^2)/(L\\sqrt\\pi)$, to obtain the S-wave phase shift. The amplitude is modelled as a scattering length with a lattice-spacing dependence, and the physical-point scattering length is obtained by chiral extrapolation. This chain is what connects lattice energy levels to the claimed bound-state pole.","core_discovery":"On the paper's own terms, the discovery is that the ground-state energy of the $bb\\bar u\\bar d$ system in the $T_{1g}$ representation lies below the $BB^*$ threshold on every ensemble, and that the scattering amplitude extracted from these levels via the standard finite-volume quantization condition is attractive enough to bind. Parametrizing the near-threshold amplitude with a scattering-length term plus a lattice-spacing correction, fitting across all pion masses, and extrapolating to the physical pion mass gives a scattering length $a_0^{\\mathrm{phys}}=0.25(4 3)$ fm. The corresponding pole sits $116^{+30}_{-36}$ MeV below the $BB^*$ threshold. The authors therefore conclude that the $I(J^P)=0(1^+)$ $bb\\bar u\\bar d$ state, $T_{bb}$, exists as a stable tetraquark.","pith_inferences":["A natural next step is to include the effective-range term in the amplitude; if that term is not small at the bound-state momentum $|p|\\sim 0.79$ GeV, the quoted binding energy could move outside its stated errors.","The same machinery applied to the $bc\\bar u\\bar d$ channel could reveal whether binding persists when one bottom quark is replaced by a charm quark, connecting to the experimentally observed doubly charmed tetraquark.","One could test the diquark-antidiquark picture by computing the overlap of the bound state with the diquark operator $\\Phi_D$; this paper uses both operator types but does not report that decomposition."],"forward_implications":["A stable $T_{bb}$ with $I(J^P)=0(1^+)$ should exist roughly $116^{+30}_{-36}$ MeV below the $BB^*$ threshold, making it stable against strong decay and potentially long-lived enough to be reconstructed experimentally.","The physical scattering length $a_0^{\\mathrm{phys}}=0.25(4 3)$ fm is a quantitative prediction that future lattice calculations with independent actions or finer lattices can check directly.","The negative ground-state shift observed at every pion mass and volume used in this study supports a genuine attractive $BB^*$ interaction rather than a finite-volume artifact.","The multi-volume amplitude analysis demonstrates a route from finite-volume spectra to near-threshold bound-state parameters for heavy tetraquarks, a route that can be applied to other doubly heavy channels."],"supporting_citations":[{"why":"The earlier lattice study by the same group that this work extends; it supplied the operators and the initial bound-state signal that the present multi-volume analysis is designed to firm up.","marker":"[5]"},{"why":"Provides the dynamical gauge-field ensembles with $N_f=2+1+1$ flavors and multiple lattice spacings used for all measurements in this paper.","marker":"[6]"},{"why":"The companion paper whose Table 1 and marker conventions define the ensemble parameters and the spectrum presentation used here.","marker":"[7]"},{"why":"Assesses the left-hand cuts from the $BB\\pi$ three-body channel, supporting the paper's choice to treat $BB^*$ scattering as elastic and to ignore that channel.","marker":"[8]"},{"why":"The standard quantization condition connecting finite-volume energy levels to infinite-volume S-wave phase shifts, the basis of the amplitude extraction.","marker":"[11]"},{"why":"Extends the quantization formalism to the elastic $BB^*$ amplitude analysis used here.","marker":"[12]"},{"why":"Provides the parametric procedure, in its Appendix B, that the paper follows to turn the extracted phase shifts into near-threshold pole parameters.","marker":"[13]"}],"fun_headline_variants":["Lattice QCD pins down doubly bottom tetraquark binding","Doubly bottom tetraquark bound by 116 MeV in lattice QCD","Stable bottom tetraquark predicted by lattice QCD","Lattice QCD finds bottom-bottom tetraquark bound","Lattice QCD: doubly bottom tetraquark stable"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The analysis assumes that a scattering-length-only, single-channel $BB^*$ amplitude describes the lattice energy levels up to the bound-state momentum; if the effective-range term is not negligible there, or if coupling to other channels shifts the level, the 116 MeV binding energy could change.","fun_headline_variants_meta":{"raw":{"variants":["Lattice QCD pins down doubly bottom tetraquark binding","Doubly bottom tetraquark bound by 116 MeV in lattice QCD","Stable bottom tetraquark predicted by lattice QCD","Lattice QCD finds bottom-bottom tetraquark bound","Lattice QCD: doubly bottom tetraquark stable"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000677,"raw_usage":{"total_tokens":3033,"prompt_tokens":852,"completion_tokens":2181,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":468,"completion_tokens_details":{"reasoning_tokens":2095}},"tokens_in":468,"tokens_out":2181,"duration_ms":16774,"temperature":1.0,"reasoning_tokens":2095,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T04:52:15.814487+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the effective range from the same finite-volume levels, or add the first inelastic channels ($B^*B^*$ and $BB\\pi$) to the quantization condition: if the pole satisfying $p\\cot\\delta_0=+\\sqrt{-p^2}$ no longer sits near 116 MeV below threshold, or if the extracted scattering length moves by more than the quoted errors, the central claim is falsified.","supporting_citations":[],"review_version":1}