REVIEW 4 major objections 4 minor 38 references
Adding at least one tetraquark operator to the correlation matrix is necessary before the finite-volume spectrum in the a0(980) and kappa channels can be trusted; without it, an energy level is missed in the K-eta subsystem below threshold.
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
T0 review · deepseek-v4-flash
2026-08-02 17:48 UTC pith:UDWLLMKD
load-bearing objection A useful but not yet airtight demonstration that omitting tetraquark operators can distort the extracted spectrum in the a0 and κ channels; the selection protocol and the threshold claim need scrutiny. the 4 major comments →
Investigating the role of tetraquark operators in lattice QCD studies of the a₀(980) and kappa resonances
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The paper's central claim is that tetraquark interpolating operators are not optional add-ons but essential elements of a reliable operator basis in these scalar channels. Working at m_pi ≈ 230 MeV and m_pi L ≈ 4.4, the authors compare spectra extracted with eleven or twelve meson operators (single-meson and meson-meson) to those obtained after adding a carefully selected tetraquark operator. In the kappa channel, the tetraquark operator resolves a previously hidden level at 2.139(62) m_K, below the K-eta threshold at 2.172 m_K, and changes the number and pattern of levels in the 1.8–2.3 m_K region. In the a0(980) channel, the added operator does not introduce a new level but substantially r
What carries the argument
The central objects are the tetraquark interpolating operators: local and extended four-quark operators constructed from two quarks and two antiquarks, in either of two color-singlet couplings (the symmetric and antisymmetric contractions of color indices). Built with various spin, displacement, and flavor structures, they provide an extra row in the temporal correlation matrix that couples to stationary states invisible to single-meson and two-meson operators. The paper scans hundreds of these operators on 25 gauge configurations and picks the one that most strongly changes the low-lying spectrum; that selected operator then resolves the previously missed level.
Load-bearing premise
The argument hinges on the selected tetraquark operator, chosen for its large effect on 25 of the same 412 configurations, genuinely coupling to a QCD stationary state rather than to noise, and on that state lying below the K-eta threshold.
What would settle it
Remove the 25 configurations used in the operator-selection scan and recompute the kappa-channel spectrum; if the extra level below the K-eta threshold disappears or shifts by more than its quoted uncertainty, the selection was overfit. Or add a second, independently chosen tetraquark operator from the same dozen that showed large effects; if the level is not recovered, the single-operator result is not robust.
If this is right
- Lattice QCD extractions of the a0(980) resonance parameters that used only single- and two-meson operators may have used an incomplete finite-volume spectrum, so their K-matrix results are suspect.
- Studies of the K0*(1430), which decays into K-eta, must include tetraquark operators or risk missing a low-lying state that changes the coupled-channel analysis.
- The kappa resonance itself, decaying primarily to K-pi, appears to be less affected: the K-pi energies change only slightly when the tetraquark operator is added.
- The new level below the K-eta threshold, if confirmed, could indicate a weakly bound or virtual K-eta state in this volume, a prediction that can be tested with more energy levels and moving frames.
Where Pith is reading between the lines
- The operator-selection procedure is the paper's main fragility: the chosen tetraquark operator was selected because it produced a dramatic change on the same 25 configurations later used in the final analysis. Repeating the selection on a training subset and confirming on a hold-out set would distinguish a genuine state from an overfit.
- A natural, testable extension is to include several of the dozen similarly-behaved tetraquark operators simultaneously; if the extra level persists with stable energy, the claim that it is a real QCD state is strengthened.
- The same scanning approach could identify analogous missing levels in other light scalar channels, such as the isoscalar sigma, where disconnected diagrams and tetraquark content are expected to matter.
- If the additional K-eta level is real, the paper's qualitative quantization-condition analysis predicts a bound K-eta state in this volume; a fully coupled multichannel fit with multiple total momenta would turn that suggestion into a falsifiable prediction.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies the role of tetraquark operators in finite-volume spectrum extraction for the isodoublet strange (κ) and isotriplet nonstrange (a0(980)) channels on an N_f=2+1 anisotropic clover ensemble with m_π≈230 MeV and m_π L≈4.4. Using stochastic LapH with all disconnected contributions, the authors construct single-meson, two-meson, and hundreds of tetraquark operators, select one tetraquark operator per channel from a 25-configuration low-statistics run, and then extract the lowest energies from 412-configuration correlation matrices with and without that operator. In the I=1/2, S=1 channel they find that inclusion of the tetraquark operator resolves an additional level near the Kη threshold; in the I=1, S=0 channel the extracted spectrum above the ground state changes dramatically. The paper concludes that spectrum determinations without at least one tetraquark operator are unreliable in these channels and discusses the consequences for Lüscher-style K-matrix parameterizations.
Significance. If correct, the result is an important caution for lattice studies of light scalar mesons: a missing operator type can cause a missed finite-volume level, with direct consequences for K-matrix analyses based on the quantization condition. The paper has real strengths: it includes all disconnected contributions, uses a systematic and flexible tetraquark operator construction, checks GEVP parameter stability (Fig. 4), and quantifies level overlaps with all operators. The K-matrix part is appropriately labeled qualitative, and the paper explicitly notes the limitations of zero-total-momentum data. However, the central claim rests on a single tetraquark operator selected in-sample from 25 of the 412 configurations used in the final analysis, and the additional level is only about 0.5σ below the Kη threshold. The result is therefore plausible and worth publishing only after the selection issue is addressed and the bound-state wording is made commensurate with the statistical evidence.
major comments (4)
- [Sec. II E and Tables VII-X] The central claim is not validated out of sample. The tetraquark operators are chosen because, in a 25-configuration subset, they produce 'a dramatic change in the low-lying energies'; the final 412-configuration analysis uses the same 25 configurations. Under the null hypothesis that the hundreds of tried operators couple mainly to noise, selecting the operator that maximizes the spectral change on 25/412 configurations will tend to produce a spectral change in the full sample because the full sample contains the selection configurations. No holdout split, no blinded analysis, and no test with a second, equally plausible tetraquark operator is reported. To support the conclusion that at least one tetraquark operator is necessary, the authors should either perform the selection without using the final data (e.g., on a different subset or ensemble) or demonstrate robustness by repeating t
- [Sec. III B and Sec. IV A (Table VIII, Fig. 11)] The additional level that carries the main message is not quantitatively below the Kη threshold. Using the numbers in Table VIII, level 3 is at 2.139(62) m_K while the Kη threshold is about 2.172 m_K, i.e., about 0.5σ below threshold; its separation from level 2 (1.951(69) m_K) is about 2σ. Thus the evidence for a bound Kη state is weak even before accounting for the look-elsewhere effect from the many operator trials in Sec. II E. The abstract's 'below the Kη threshold' and Sec. IV A's 'suggest the existence of a bound Kη state' are therefore stronger than the data support. These statements should be softened, or the authors should provide a global significance that includes the operator search and fit variations.
- [Sec. II E and Table VI] The exact definition of the selected tetraquark operators is not given: the text states that superposition coefficients and the 'SS2' construction details are available upon request. Since the entire result is tied to the specific operators used, this is not a presentation detail. The group-theoretical coefficients, spin contractions, displacement patterns, and the exact SS2 operator should appear in the paper or in an appendix; otherwise the result cannot be independently reproduced or tested with a closely related but different operator.
- [Sec. IV A, Fig. 11] The no-tetraquark comparison in the Kη channel is not quantitative. The paper shows that constant and linear forms of K^{-1} can be tuned to intersect the box matrix at energies consistent with the with-tetraquark red points, but no goodness-of-fit is reported for the no-tetraquark orange points, and no model comparison or likelihood ratio between the 'missed level' and 'no missed level' hypotheses is made. The statement that omitting the tetraquark operator 'would make it difficult' to find a suitable K^{-1} is an assertion; a simple χ² comparison for the two hypotheses would make the argument concrete.
minor comments (4)
- [Sec. II C, Eq. (14)] The text says q and q-bar are defined in Eq. (4) of Ref. [23], but the definitions appear in Eq. (5) of this paper; the cross-reference should be corrected.
- [Sec. III A] The η′ is approximated as stable; this is reasonable for the stated masses but should be listed explicitly as a systematic approximation alongside the exclusion of four-meson operators.
- [Figs. 6 and 9] The figure captions use physical η and φ for non-interacting energies while the operator labels use η̃ and φ̃. This distinction is described in the text but easily missed; a brief clarification in each caption would improve readability.
- [Sec. II E] The phrase 'the operator yielding the effective energy with the smallest statistical errors and least excited-state contamination' should include a concrete measure of 'least excited-state contamination', since this selection criterion is otherwise not reproducible.
Circularity Check
No significant circularity; the central result is a numerically extracted spectrum, not an identity derived from its inputs.
full rationale
The paper's central claim is that including a tetraquark operator changes or adds finite-volume energy levels in the a0(980) and kappa channels. This is a numerical spectrum extraction, not a derived identity, so it does not reduce to an input by construction. The closest concern is the operator selection in Sec. II E: hundreds of tetraquark operators are screened on 25 configurations, and the operator producing a 'dramatic change' is retained, while the final spectrum uses all 412 configurations including the same 25. That is an in-sample selection effect and a legitimate statistical worry, but it is not circularity: the paper does not fit a parameter to data and then rename it a prediction; it transparently reports the selection procedure and then extracts the spectrum with a larger ensemble. The additional level at 2.139(62) m_K is only mildly below the K-eta threshold, which is a statistical strength issue, not a circularity issue. No load-bearing self-citation chain is used, no uniqueness theorem is imported from the authors' previous work, and no ansatz is smuggled in via citation as a substitute for derivation. The Lüscher analysis uses the same extracted energies, but the paper explicitly labels it exploratory and qualitative, and the conclusion is conditional on the extracted spectrum. The lack of full operator definitions in the paper ('available from the authors upon request') hinders reproducibility but is not circularity. Overall, the derivation chain is self-contained in the sense required by the circularity test.
Axiom & Free-Parameter Ledger
free parameters (4)
- Selected tetraquark operator (discrete choice) =
κ: T[suss]SS2(−); a0: T[(uu+dd)du]SS2(+)
- GEVP diagonalization time τ_D =
12 (κ), 7 (a0)
- Two-exponential fit windows (τ_min, τ_max) =
Per level, Tables VII–X (e.g., (7,26), (8,26), ...)
- K-matrix ansatz coefficients in Sec. IV =
constant, linear, quadratic forms shown in Figs. 10–13
axioms (5)
- domain assumption The stochastic LapH estimates of the correlators are unbiased and the employed noise dilution schemes (TF/SF/LI8 and TI16/SF/LI8) are sufficient for the precision needed here.
- domain assumption The Lüscher quantization condition (Eq. 19) with ℓ≤0 truncation and negligible channel mixing is valid for these zero-total-momentum channels below three-particle thresholds.
- domain assumption The η and η′ mesons may be treated as stable on this ensemble.
- standard math Isospin symmetry with m_u=m_d and no electromagnetic/weak effects.
- domain assumption The single-pivot GEVP with the chosen τ_N, τ_0, τ_D cleanly separates the relevant finite-volume states.
read the original abstract
The role of tetraquark operators in studying the isodoublet strange $\kappa$ and isovector nonstrange $a_0(980)$ scalar mesons in lattice QCD is examined using an ensemble with $m_\pi\approx230$ MeV and spatial extent $L$ such that $m_\pi L\approx4.4$. Hermitian correlation matrices using both single-meson, meson-meson, and tetraquark interpolating operators are used to extract the spectrum of finite-volume stationary states in the appropriate symmetry channels. Hundreds of local and extended tetraquark operators are explored. Determinations of the spectrum in each channel are found to be unreliable without the inclusion of at least one tetraquark operator. For example, the inclusion of tetraquark operators with isospin 1/2 and strangeness 1 quantum numbers reveals the existence of an additional energy level in the $K\eta$ sub-system below the $K\eta$ threshold. The implications of this on parametrizing the scattering $K$-matrix through a well-known quantization condition to extract properties of the $\kappa$ and $a_0(980)$ scalar meson resonances are discussed.
Figures
Reference graph
Works this paper leans on
-
[1]
box matrix
This state is mainly created by the K(0)K(0) op- erator, so level 1 can be interpreted as a kaon-antikaon state with residual interactions. Levels 2 and 3 are pro- duced mainly by the tetraquark and ˜ϕ(0)π(0) operators, with small overlap factors from the ˜η(0)π(0) and ROT 0 operators. This leads to an identification of level 2 as predominantly a tetraqua...
-
[2]
G. Rupp and E. Van Beveren, Scalar mesons: fifty years of challenging the quark model, Acta Phys. Polon. Supp. 11, 455 (2018), arXiv:1806.00364 [hep-ph]
Pith/arXiv arXiv 2018
-
[3]
Godfrey and N
S. Godfrey and N. Isgur, Mesons in a Relativized Quark Model with Chromodynamics, Phys. Rev. D32, 189 (1985)
1985
-
[4]
An inspection of the overlap factors in Fig
The tetraquark operator is observed to impact levels 2 and 3 mostly. An inspection of the overlap factors in Fig. 5 shows thatKπoperators produce states with very little overlaps on levels 2 and 3. These facts indicate that the tetraquark operator is having an effect mainly in the Kηdecay channel and not in theKπchannel. We conclude that the tetraquark op...
2047
-
[5]
Jaffe, Exotica, Nucl
R. Jaffe, Exotica, Nucl. Phys. B Proc. Suppl.142, 343 (2005)
2005
-
[6]
Amsler and N
C. Amsler and N. Tornqvist, Mesons beyond the naive quark model, Phys. Rept.389, 61 (2004)
2004
-
[7]
F. E. Close and N. A. Tornqvist, Scalar mesons above and below 1-GeV, J. Phys. G28, R249 (2002), arXiv:hep- ph/0204205
arXiv 2002
-
[8]
L. Maiani, F. Piccinini, A. Polosa, and V. Riquer, A New look at scalar mesons, Phys. Rev. Lett.93, 212002 (2004), arXiv:hep-ph/0407017
Pith/arXiv arXiv 2004
-
[9]
A. Stump and J. R. Green, Position-space sampling for local multiquark operators in lattice QCD using distil- lation and the importance of tetraquark operators for Tcc(3875)+ (2025), arXiv:2510.26459 [hep-lat]
arXiv 2025
-
[10]
S. Prelovsek, E. Ortiz-Pacheco, S. Collins, L. Leskovec, M. Padmanath, and I. Vujmilovic, Doubly heavy tetraquarks from lattice QCD: Incorporating diquark- antidiquark operators and the left-hand cut, Phys. Rev. D112, 014507 (2025), arXiv:2504.03473 [hep-lat]
Pith/arXiv arXiv 2025
-
[11]
H. Alharazin, A. B. Raposo, J. Bulava, S. Dawid, J. R. Green, C. Morningstar, F. Romero-L´ opez, M. Salg, S. R. Sharpe, and A. Stump, Three-body study of the Tcc(3875)+ from lattice QCD, in42th International Sym- posium on Lattice Field Theory(2026) arXiv:2602.17204 [hep-lat]
arXiv 2026
-
[12]
S. Prelovsek, T. Draper, C. B. Lang, M. Limmer, K.-F. Liu, N. Mathur, and D. Mohler, Lattice study of light scalar tetraquarks with I=0,2,1/2,3/2: Areσ andκtetraquarks?, Phys. Rev. D82, 094507 (2010), arXiv:1005.0948 [hep-lat]
Pith/arXiv arXiv 2010
-
[13]
C. Alexandrou, J. O. Daldrop, M. Dalla Brida, M. Grav- ina, L. Scorzato, C. Urbach, and M. Wagner, Lattice in- vestigation of the scalar mesonsa 0(980) andκusing four- quark operators, JHEP04, 137, arXiv:1212.1418 [hep- lat]
-
[14]
C. Alexandrou, J. Berlin, M. Dalla Brida, J. Finkenrath, T. Leontiou, and M. Wagner, Lattice QCD investigation of the structure of thea 0(980) meson, Phys. Rev. D97, 034506 (2018), arXiv:1711.09815 [hep-lat]
Pith/arXiv arXiv 2018
-
[15]
F.-K. Guo, L. Liu, U.-G. Meissner, and P. Wang, Tetraquarks, hadronic molecules, meson-meson scatter- ing and disconnected contributions in lattice QCD, Phys. Rev. D88, 074506 (2013), arXiv:1308.2545 [hep-lat]
Pith/arXiv arXiv 2013
-
[16]
C. Morningstar, J. Bulava, J. Foley, K. J. Juge, D. Lenkner, M. Peardon, and C. H. Wong, Improved stochastic estimation of quark propagation with Lapla- cian Heaviside smearing in lattice QCD, Phys. Rev. D 83, 114505 (2011), arXiv:1104.3870 [hep-lat]. 22
Pith/arXiv arXiv 2011
-
[17]
Luscher, Two particle states on a torus and their re- lation to the scattering matrix, Nucl
M. Luscher, Two particle states on a torus and their re- lation to the scattering matrix, Nucl. Phys.B354, 531 (1991)
1991
-
[18]
K. Rummukainen and S. A. Gottlieb, Resonance scatter- ing phase shifts on a nonrest frame lattice, Nucl. Phys. B450, 397 (1995), arXiv:hep-lat/9503028 [hep-lat]
Pith/arXiv arXiv 1995
-
[19]
C. H. Kim, C. T. Sachrajda, and S. R. Sharpe, Finite- volume effects for two-hadron states in moving frames, Nucl. Phys.B727, 218 (2005), arXiv:hep-lat/0507006 [hep-lat]
Pith/arXiv arXiv 2005
-
[20]
R. A. Briceno, Two-particle multichannel systems in a fi- nite volume with arbitrary spin, Phys. Rev.D89, 074507 (2014), arXiv:1401.3312 [hep-lat]
Pith/arXiv arXiv 2014
-
[21]
G. Fox, R. Gupta, O. Martin, and S. Otto, Monte Carlo Estimates of the Mass Gap of the O(2) and O(3) Spin Models in (1+1)-dimensions, Nucl. Phys. B205, 188 (1982)
1982
-
[22]
Michael and I
C. Michael and I. Teasdale, Extracting Glueball Masses From Lattice QCD, Nucl. Phys.B215, 433 (1983)
1983
-
[23]
L¨ uscher and U
M. L¨ uscher and U. Wolff, How to Calculate the Elastic Scattering Matrix in Two-dimensional Quantum Field Theories by Numerical Simulation, Nucl. Phys. B339, 222 (1990)
1990
- [24]
-
[25]
C. Morningstar, J. Bulava, B. Fahy, J. Foley, Y. Jhang, K. Juge, D. Lenkner, and C. Wong, Extended hadron and two-hadron operators of definite momentum for spectrum calculations in lattice QCD, Phys. Rev. D88, 014511 (2013), arXiv:1303.6816 [hep-lat]
Pith/arXiv arXiv 2013
-
[26]
C. Morningstar and M. J. Peardon, Analytic smearing of SU(3) link variables in lattice QCD, Phys. Rev. D69, 054501 (2004), arXiv:hep-lat/0311018
Pith/arXiv arXiv 2004
-
[27]
R. G. Edwards, B. Joo, and H.-W. Lin, Tuning for Three-flavors of Anisotropic Clover Fermions with Stout-link Smearing, Phys. Rev. D78, 054501 (2008), arXiv:0803.3960 [hep-lat]
Pith/arXiv arXiv 2008
-
[28]
H.-W. Linet al.(Hadron Spectrum), First results from 2+1 dynamical quark flavors on an anisotropic lattice: Light-hadron spectroscopy and setting the strange-quark mass, Phys. Rev. D79, 034502 (2009), arXiv:0810.3588 [hep-lat]
Pith/arXiv arXiv 2009
-
[29]
R. Brett, J. Bulava, J. Fallica, A. Hanlon, B. H¨ orz, and C. Morningstar, Determination ofs- andp-waveI= 1/2 Kπscattering amplitudes inN f = 2 + 1 lattice QCD, Nucl. Phys. B932, 29 (2018), arXiv:1802.03100 [hep-lat]
Pith/arXiv arXiv 2018
-
[30]
J. Bulava, B. Fahy, B. H¨ orz, K. J. Juge, C. Morningstar, and C. H. Wong,I= 1 andI= 2π−πscattering phase shifts fromN f = 2 + 1 lattice QCD, Nucl. Phys. B910, 842 (2016), arXiv:1604.05593 [hep-lat]
Pith/arXiv arXiv 2016
-
[31]
D. J. Wilson, J. J. Dudek, R. G. Edwards, and C. E. Thomas, Resonances in coupledπK, ηKscatter- ing from lattice QCD, Phys. Rev. D91, 054008 (2015), arXiv:1411.2004 [hep-ph]
Pith/arXiv arXiv 2015
-
[32]
G. Rendon, L. Leskovec, S. Meinel, J. Negele, S. Paul, M. Petschlies, A. Pochinsky, G. Silvi, and S. Syritsyn, I= 1/2S-wave andP-waveKπscattering and theκ andK ∗ resonances from lattice QCD, Phys. Rev. D102, 114520 (2020), arXiv:2006.14035 [hep-lat]
Pith/arXiv arXiv 2020
-
[33]
J. J. Dudek, R. G. Edwards, and C. E. Thomas (Hadron Spectrum), Energy dependence of theρresonance in ππelastic scattering from lattice QCD, Phys. Rev. D 87, 034505 (2013), [Erratum: Phys.Rev.D 90, 099902 (2014)], arXiv:1212.0830 [hep-ph]
Pith/arXiv arXiv 2013
-
[34]
J. J. Dudek, R. G. Edwards, and D. J. Wilson (Hadron Spectrum), Ana 0 resonance in strongly coupledπη,K K scattering from lattice QCD, Phys. Rev. D93, 094506 (2016), arXiv:1602.05122 [hep-ph]
Pith/arXiv arXiv 2016
-
[35]
C. Morningstar, J. Bulava, B. Singha, R. Brett, J. Fal- lica, A. Hanlon, and B. H¨ orz, Estimating the two-particle K-matrix for multiple partial waves and decay chan- nels from finite-volume energies, Nucl. Phys. B924, 477 (2017), arXiv:1707.05817 [hep-lat]
Pith/arXiv arXiv 2017
-
[36]
Navaset al.(Particle Data Group), Review of particle physics, Phys
S. Navaset al.(Particle Data Group), Review of particle physics, Phys. Rev. D110, 030001 (2024)
2024
-
[37]
Stanzione, J
D. Stanzione, J. West, R. Evans, T. Minyard, O. Ghat- tas, and D. Panda, Frontera: The evolution of leadership computing at the national science foundation, inProceed- ings of Practice and Experience in Advanced Research Computing (PEARC ’20)(2020)
2020
-
[38]
R. G. Edwards and B. Joo (SciDAC), The Chroma soft- ware system for lattice QCD, Nucl. Phys. Proc. Suppl. 140, 832 (2005)
2005
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.