Lattice QCD Study of Doubly Heavy Bottom Tetraquarks

Bhabani Sankar Tripathy; M. Padmanath; Nilmani Mathur

REVIEW 3 major objections 4 minor 13 references

Lattice QCD Study of Doubly Heavy Bottom Tetraquarks

T0 review · 3 major / 4 minor · reviewed 2026-08-16 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Lattice QCD finds the doubly bottom tetraquark is bound

desk verdict 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. read the letter →

arxiv 2505.00078 v1 pith:IKBF4VJU submitted 2025-04-30 hep-lat hep-exhep-ph

Bhabani Sankar Tripathy , Nilmani Mathur , M. Padmanath This is my paper · ORCID

classification hep-lathep-exhep-ph

keywords latticeQCDtetraquarksdoublyheavybottomquarksexotichadronsBB*scatteringfinite-volumespectrumquantization

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

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.

What carries the argument

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.

What would settle it

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.

Watch

Extended reading notes

Core claim

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.

Load-bearing premise

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.

Editorial extensions

If this is right

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.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

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.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the 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.

What 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.

The 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.

I'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.

Referee Report

3 major / 4 minor

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.

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 (3)

[Section 3, 'Amplitude analysis'] 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 3, Fig. 1 and final paragraph] 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 3, final paragraph] 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.

minor comments (4)

[Abstract and Conclusion] There are several typographical errors, including 'within of these exotic states' in the Abstract and 'exitence' in the Conclusion; these should be corrected.
[Section 3, Eq. (1)] 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 3, final paragraph] 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.
[Figure 1 caption] 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.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the binding energy is extracted from lattice finite-volume spectra via a standard Lüscher-style amplitude analysis, not assumed or self-referential.

full rationale

The derivation chain is: finite-volume energy shifts are measured from lattice correlators (Eq. 1 and GEVP Eq. 3); these are converted to S-wave p cot δ0 through the Lüscher quantization condition; p cot δ0 values are fitted with the model f = A[0] + A[1]·a; after chiral and continuum extrapolation the leading term gives the scattering length a0 = 0.25(4)(3) fm; the binding energy is then obtained by solving the bound-state pole condition p cotδ0 = +√(−p²). Each step uses independent lattice data and standard relations; the 116 MeV binding energy is a computed output of the fitted scattering length, not an input to the fit. The paper's reliance on Refs. [5,7,13] is citation of data sources, ensemble conventions, and a quantization-condition appendix by overlapping authors, but none of these citations supplies the central result or forbids alternatives. The neglect of effective-range and left-hand-cut terms is an acknowledged systematic and model-dependence issue, not a circularity: those terms are not set to zero by definition of the fitted binding energy, they are simply omitted. The paper is self-contained against external benchmarks in the sense that the lattice data determine the extracted amplitude, and the quoted binding energy would change if the amplitude model changed, which is a correctness risk rather than a logical circularity.

Assumptions & free parameters 3 free parameters · 4 assumptions · 0 invented entities

The calculation is a standard lattice QCD computation with no invented entities. The fitted parameters are the extrapolation coefficients and the implicit zero of the effective range. The main model dependence is the elastic single-channel truncation and the scattering-length-only expansion, both assumptions rather than free parameters in the usual sense.

free parameters (3)

Continuum extrapolation coefficients A0 and A1 = not reported
f = A[0] + A[1]·a is fitted to the p cot δ0 values across lattice spacings; only the final scattering length a0 = 0.25 fm is quoted.
Chiral extrapolation coefficients c0 and c1 = not reported
Fits f1 = c0 + c1 Mps and f2 = c0 + c1 Mps^2 to the continuum-extrapolated p cot δ0; the choice between the two forms is not justified.
Scattering length effective range parameter = 0 (implicitly)
The effective-range term is dropped, leaving a single scattering length to encode the near-threshold amplitude; this is the key model assumption in the binding-energy extraction.

assumptions (4)

domain assumption MILC Nf=2+1+1 HISQ ensembles are a valid QCD discretization
Adopted from Ref [6]; no validation shown in this paper.
domain assumption NRQCD action describes the bottom quark at the lattice spacings used
The b mass is tuned via spin-averaged 1S bottomonium; this is standard but not independently verified here.
domain assumption The finite-volume spectrum is described by the single-channel elastic Lüscher quantization condition with only B B*
Section 3 states that B*B* and BBπ are beyond scope; this truncation is load-bearing for the amplitude analysis.
domain assumption The GEVP ground-state energy is dominated by the lowest B B* scattering state
The paper shows a plot but no explicit fits or saturation tests.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Lattice QCD Study of Doubly Heavy Bottom Tetraquarks." pith.science (2026). https://pith.science/paper/IKBF4VJU

@misc{pith2026250500078,
  author       = {Pith},
  title        = {Pith review of: Lattice QCD Study of Doubly Heavy Bottom Tetraquarks},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/IKBF4VJU}},
  note         = {Machine review of arXiv:2505.00078}
}

abstract

Hadrons, composed of quarks and gluons bound by Quantum Chromodynamics (QCD), traditionally classified as baryons (three quarks) and mesons (quark-antiquark pairs). Nothing in the theory of QCD stands against the existence of exotic hadrons with more complex quark contents. Recent discoveries by LHCb and Belle, such as X, Y, Z states and $T_{cc}(3875)$, have renewed interest in these states. Understanding the binding mechanism within of these exotic states provides insights into QCD's non-perturbative dynamics. This work presents lattice QCD studies of two-meson interactions, involving bottom quarks, on MILC ensembles, exploring heavy tetraquark channels.

Figures

Figures reproduced from arXiv: 2505.00078 by the authors.

**Figure 1.** (Left) Ground state energies in units of the BB∗ threshold energy on all ensembles and for all Mps. (Right) Continuum extrapolated p cot δ0 estimates of the BB∗ system as a function of M2 ps in units of EBB∗ . ground state energy estimates relative to the nearest two-body decay threshold BB∗ across various Mps and ensembles, following the marker and color conventions of Ref. [7]. Faded markers represent first exci… view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

13 extracted references · 1 canonical work pages

[7]

B. S. Tripathy, N. Mathur and M. Padmanath, [arXiv:2503.09760 [hep-lat]]

arXiv
[1]

R. F. Lebed, R. E. Mitchell and E. S. Swanson, Prog. Part. Nucl. Phys. 93, 143-194 (2017) doi:10.1016/j.ppnp.2016.11.003 [arXiv:1610.04528 [hep-ph]]

arXiv 2017
[2]

F. K. Guo, C. Hanhart, U. G. Meißner, Q. Wang, Q. Zhao and B. S. Zou, Rev. Mod. Phys. 90, no.1, 015004 (2018) [erratum: Rev. Mod. Phys. 94, no.2, 029901 (2022)] doi:10.1103/RevModPhys.90.015004 [arXiv:1705.00141 [hep-ph]]

arXiv 2018
[3]

Bicudo, Phys

P. Bicudo, Phys. Rept. 1039, 1-49 (2023) doi:10.1016/j.physrep.2023.10.001 [arXiv:2212.07793 [hep-lat]]

arXiv 2023
[4]

Francis, Prog

A. Francis, Prog. Part. Nucl. Phys. 140, 104143 (2025) doi:10.1016/j.ppnp.2024.104143

arXiv 2025
[5]

Junnarkar, N

P. Junnarkar, N. Mathur and M. Padmanath, Phys. Rev. D 99, no.3, 034507 (2019) doi:10.1103/PhysRevD.99.034507 [arXiv:1810.12285 [hep-lat]]

arXiv 2019
[6]

Bazavov et al

A. Bazavov et al. [MILC], Phys. Rev. D 87, no.5, 054505 (2013) doi:10.1103/PhysRevD.87.054505 [arXiv:1212.4768 [hep-lat]]

arXiv 2013
[8]

M. L. Du, A. Filin, V. Baru, X. K. Dong, E. Epelbaum, F. K. Guo, C. Hanhart, A. Nefediev, J. Nieves and Q. Wang, Phys. Rev. Lett. 131, no.13, 131903 (2023) doi:10.1103/PhysRevLett.131.131903 [arXiv:2303.09441 [hep-ph]]

arXiv 2023

Show all 13 references

[9]

Meinel, M

S. Meinel, M. Pflaumer and M. Wagner, Phys. Rev. D 106, no.3, 034507 (2022) doi:10.1103/PhysRevD.106.034507 [arXiv:2205.13982 [hep-lat]]

2022 arXiv
[10]

Michael, Nucl

C. Michael, Nucl. Phys. B 259, 58-76 (1985) doi:10.1016/0550-3213(85)90297-4

1985 doi
[11]

Luscher, Nucl

M. Luscher, Nucl. Phys. B 354, 531-578 (1991) doi:10.1016/0550-3213(91)90366-6

1991 doi
[12]

R. A. Briceno, Phys. Rev. D 89, no.7, 074507 (2014) doi:10.1103/PhysRevD.89.074507 [arXiv:1401.3312 [hep-lat]]

2014 arXiv
[13]

Padmanath and S

M. Padmanath and S. Prelovsek, Phys. Rev. Lett. 129, no.3, 032002 (2022) doi:10.1103/PhysRevLett.129.032002 [arXiv:2202.10110 [hep-lat]]

2022 arXiv

Pith tools

Reviewed August 16, 2026 · model on record in the stance chip above.

[1] [7]

B. S. Tripathy, N. Mathur and M. Padmanath, [arXiv:2503.09760 [hep-lat]]

arXiv

[2] [1]

R. F. Lebed, R. E. Mitchell and E. S. Swanson, Prog. Part. Nucl. Phys. 93, 143-194 (2017) doi:10.1016/j.ppnp.2016.11.003 [arXiv:1610.04528 [hep-ph]]

arXiv 2017

[3] [2]

F. K. Guo, C. Hanhart, U. G. Meißner, Q. Wang, Q. Zhao and B. S. Zou, Rev. Mod. Phys. 90, no.1, 015004 (2018) [erratum: Rev. Mod. Phys. 94, no.2, 029901 (2022)] doi:10.1103/RevModPhys.90.015004 [arXiv:1705.00141 [hep-ph]]

arXiv 2018

[4] [3]

Bicudo, Phys

P. Bicudo, Phys. Rept. 1039, 1-49 (2023) doi:10.1016/j.physrep.2023.10.001 [arXiv:2212.07793 [hep-lat]]

arXiv 2023

[5] [4]

Francis, Prog

A. Francis, Prog. Part. Nucl. Phys. 140, 104143 (2025) doi:10.1016/j.ppnp.2024.104143

arXiv 2025

[6] [5]

Junnarkar, N

P. Junnarkar, N. Mathur and M. Padmanath, Phys. Rev. D 99, no.3, 034507 (2019) doi:10.1103/PhysRevD.99.034507 [arXiv:1810.12285 [hep-lat]]

arXiv 2019

[7] [6]

Bazavov et al

A. Bazavov et al. [MILC], Phys. Rev. D 87, no.5, 054505 (2013) doi:10.1103/PhysRevD.87.054505 [arXiv:1212.4768 [hep-lat]]

arXiv 2013

[8] [8]

M. L. Du, A. Filin, V. Baru, X. K. Dong, E. Epelbaum, F. K. Guo, C. Hanhart, A. Nefediev, J. Nieves and Q. Wang, Phys. Rev. Lett. 131, no.13, 131903 (2023) doi:10.1103/PhysRevLett.131.131903 [arXiv:2303.09441 [hep-ph]]

arXiv 2023