REVIEW 2 major objections 4 minor 42 references
Hybrid static potentials and gluelumps on $N_f=3+1$ ensembles
T0 review · 2 major / 4 minor · reviewed 2026-08-10 · deepseek-v4-flash
Pith's one-line read This paper reports new lattice QCD measurements of hybrid static potentials and gluelump masses on 3+1-flavor ensembles with 420 MeV pions, including string-breaking thresholds.
desk verdict New N_f=3+1 hybrid potential and gluelump data with honest caveats, but the truncated Laplace basis needs a convergence check. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the Laplace trial state: a correlator in which the spatial Wilson line of a Wilson loop is replaced by products of Laplace eigenvectors $v_i(\vec x)$ and static perambulators $\tau_{ij}=v_i^\dagger U_t v_j$, with gluonic excitations inserted as symmetric covariant derivatives $\nabla_{\vec k} v(\vec x)=\frac12[U_k(\vec x)v(\vec x+\hat k)-U_k^\dagger(\vec x-\hat k)v(\vec x-\hat k)]$ of the eigenvectors. Gaussian profile weights $\rho_k(\lambda_i)=e^{-\lambda_i^2/4\sigma_k^2}$ assign different weights to eigenmodes, and a generalized eigenvalue problem selects the optimal profiles. This machinery lets one compute correlations at arbitrary on- and off-axis separations without expensive loop re-evaluations. For gluelumps the paper uses a complementary tool: a 35-shape basis of spatial Wilson loops, projected to cubic-group irreps $A_1,A_2,E,T_1,T_2$ with parity and charge conjugation, connected by adjoint static lines, with the adjoint self-energy related to the fundamental one by the Casimir ratio $C_A/C_F=N^2/(N^2-1)$.
What would settle it
Compute the same hybrid potentials and gluelump masses on the A2 ensemble with twice the number of kept Laplace eigenvectors, or with standard Wilson-loop hybrid operators, and look for differences beyond statistical errors; agreement would validate the truncated basis, and disagreement would show it biases the results.
Extended reading notes
Core claim
The paper claims that hybrid static potentials extracted from Laplace trial states built from covariant derivatives of lattice Laplace eigenmodes reproduce the expected gluelump multiplets in the $R\to0$ limit and, at pion masses near 420 MeV, show the string-breaking thresholds where $\Sigma_g$ and $\Pi_u$ hybrids meet static-light meson states. On the $N_f=3+1$ ensemble A2 ($m_\pi\simeq 409$ MeV, $48^3$ spatial volume) this threshold lies within the lattice, so the potentials are measured across the crossing region. The gluelump spectrum is obtained from a large basis of spatial loop shapes projected onto irreducible representations of the cubic group, with masses given relative to the lowest $T_1^{+-}$ state; the $E^{PC}$ and $T_2^{PC}$ channels are consistent for each $PC$, which is the expected continuum degeneracy for $J=2$. The paper therefore presents these results as new lattice QCD inputs for effective-theory descriptions of hybrid and exotic mesons.
Load-bearing premise
The whole extraction rests on the assumption that keeping only the lowest 100 or 200 Laplace-eigenvector modes on each time slice gives a complete enough basis for the hybrid and gluelump states, a point the paper does not test by varying that number.
Editorial extensions
If this is right
- On ensemble A2 the $\Sigma_g$ and $\Pi_u$ hybrid potentials cross the static-light thresholds within the measured range, so hybrid string breaking can be studied directly on that ensemble.
- The $E^{PC}$ and $T_2^{PC}$ gluelump channels agree within errors on each ensemble, confirming that the loop-shape basis keeps cutoff effects small enough for the continuum $J=2$ degeneracy to be visible.
- Gluelump masses rise systematically as the pion mass is lowered from the quenched value through 2.2 GeV and 788 MeV to 406 MeV, so dynamical fermion effects matter for these masses.
- Because Laplace-trial-state correlators give easy access to off-axis separations, the hybrid spectrum can be resolved at many more $R$ values than with standard Wilson-loop operators.
Reading between the lines
- Because the Laplace-trial-state machinery computes correlations at arbitrary separations from precomputed perambulators, the same construction should transfer to tetraquark and multiquark static potentials more cheaply than Wilson-loop operators; the paper lists this as a plan, so the transfer is an expectation rather than a demonstrated result.
- The trend shown for gluelump masses — larger masses at smaller pion mass — suggests that Born-Oppenheimer hybrid models calibrated on quenched or heavy-pion gluelumps may misplace hybrid levels at physical pion mass; this extrapolation is ours, not the paper's.
- A decisive test of the operator replacement would be to compare these Laplace-trial-state potentials with Wilson-loop hybrid potentials on the same A2 ensemble; the paper does not show such a comparison.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This proceedings paper presents lattice QCD computations of hybrid static quark-antiquark potentials using a novel "Laplace trial state" construction, in which gluonic excitations are introduced through covariant derivatives of the lowest Laplacian eigenvectors, with Gaussian profiles and a GEVP used to improve overlaps. The same paper computes gluelump masses on an extended basis of 35 spatial loop shapes for ensembles with N_f=0, 2, and 3+1 dynamical quarks (pion masses down to 406 MeV). The hybrid potentials are compared with static-light S- and P-wave thresholds, and the lowest gluelump masses are overlaid at R->0 as the expected limit of the hybrid potentials. Statistical errors are estimated with the Gamma method via pyerrors.
Significance. If the results are correct, they provide new lattice-QCD inputs for Born-Oppenheimer effective models of hybrid mesons, including excited hybrid channels and string-breaking thresholds for N_f=3+1 QCD with physical-ish pions. The Laplace trial-state method offers practical advantages: off-axis separations can be computed easily, and perambulator-based contractions open the way to studies of string breaking. The internal cross-checks are a genuine strength: the E and T2 gluelump representations agree for the J=2 channel, and the A1/A2 volumes give consistent gluelump masses. No machine-checked proofs or released code accompany the paper, so reproducibility rests on the textual description. The central new systematic—truncation of the Laplacian eigenbasis—is not tested, which limits the strength of the claims at this stage.
major comments (2)
- [Section 2, Eqs. for C_Sigma and C_Pi, Table 1] The hybrid potentials are constructed from the lowest N_vec Laplacian eigenvectors (N_vec=100 on A1/A2 and 200 on A1h, Table 1), and the correlation functions sum only over these modes and their covariant derivatives. The paper provides no convergence test with respect to N_vec and no comparison of the resulting potentials with standard Wilson-loop hybrid potentials (e.g., Ref. [6]). Since the trial-state space is confined to these modes, an incomplete basis could bias every extracted V(R) and the R->0 comparison with gluelumps. The authors should add a controlled check, for example by repeating one ensemble with a larger N_vec, or by comparing at least the ground-state Sigma_g^+ potential with a Wilson-loop calculation on the same ensembles, and discuss the expected size of the truncation effect.
- [Section 4, Figure 2 and Eq. (3)] The overlay of gluelump masses on the hybrid potentials at R->0 relies on the approximate self-energy relation m_self = (N^2/(N^2-1)) V_self, and the text states that subtracting twice the static-light S-wave mass "does not completely remove" the self-energy contribution. The figure shows no systematic uncertainty for this residual self-energy difference. Because the claim that gluelump masses are the R->0 limit of the hybrid potentials is central, the authors should either quantify the residual mismatch (for instance by estimating V_self on the same ensembles through a Wilson-loop computation) or explicitly state that the R->0 correspondence is qualitative and carries an unquantified self-energy uncertainty.
minor comments (4)
- [Section 2, first paragraph] The symbol \bar{\tau}_{ij} is introduced in the text but is not used in the displayed equations, where only \tau_{ij} and \tau_{ji} appear; please clarify the notation for the conjugate perambulator.
- [Figure 2 and caption] The labels for excited channels such as \Sigma_g' and \Pi_u' are not defined in the caption or the text; please add a brief explanation, for example that these are the first excited states obtained from the GEVP.
- [Section 3, Table 1 and Section 4] The text refers to "A1 heavy" while the table uses "A1h"; to avoid confusion, a single consistent ensemble name should be used throughout.
- [Section 4, first paragraph] The statement that "the individual energies consistently increase with decreasing pion mass" is ambiguous, since the plotted quantities are mass differences with respect to the T_1^{+-} state; please specify whether this refers to the absolute gluelump masses or to the relative splittings.
Circularity Check
No significant circularity: the hybrid potentials and gluelump masses are independent correlation-function extractions, and their R->0 agreement is a consistency check, not a fitted constraint.
full rationale
The paper's central outputs are effective energies extracted from Euclidean correlators: C_Sigma and C_Pi built from Laplacian eigenvector trial states with covariant-derivative insertions, and C_Lambda^PC built from spatial Wilson loops with adjoint static lines. These are different operator sets; neither is defined in terms of the other's output. The Gaussian profiles and GEVP select superpositions within each correlator but do not enforce any target mass. The gluelump masses are not obtained by extrapolating the hybrid potentials to R=0; they are computed separately and overlaid on the potential plots ('We include the lowest gluelump masses corresponding to the R->0 limit... the limits look reasonable'), which is a comparison, not a derivation. The static-light thresholds m_S, m_P1/2, m_P3/2 are taken from the authors' companion paper [42]; they enter only as reference levels for plotting V(R)-2m_S and marking string-breaking thresholds, so they do not fix the potential values. Methodological self-citations [14,15] introduce the Laplace trial-state construction, but the present paper applies that construction to new hybrid and gluelump observables; the cited results are not the target result. The acknowledged limitation--that the gluelump self-energy subtraction 'does not completely remove' the mismatch--and the absence of a Wilson-loop cross-check are accuracy and convergence concerns, not evidence that any output is equivalent to an input by construction.
Assumptions & free parameters
free parameters (3)
- Gaussian profile widths sigma_k =
7 widths, values not listed
- Laplacian eigenvector truncation N_l_vec =
100 (A1, A2), 200 (A1h)
- APE and HYP2 smearing parameters =
APE 20 steps alpha=0.5; HYP2 alpha1=1, alpha2=1, alpha3=0.5
assumptions (4)
- domain assumption The covariant-derivative operators built from Laplacian eigenvectors project onto the claimed hybrid quantum numbers (Sigma_g/u, Pi_g/u).
- ad hoc to paper The truncated Laplacian eigenbasis (100 to 200 vectors) is sufficient for hybrid and gluelump correlation functions.
- domain assumption The adjoint self-energy is related to the fundamental self-energy by m_self = (N^2/(N^2-1)) V_self (Eq. 3), so gluelump masses can be plotted together with hybrid potentials without an additional shift.
- domain assumption Static-light S-wave and P-wave threshold masses from the companion distillation study [42] provide the correct string-breaking scales.
Cite this review
Pith. "Pith review of Hybrid static potentials and gluelumps on $N_f=3+1$ ensembles." pith.science (2026). https://pith.science/paper/PG7P2M55
@misc{pith2026250115670,
author = {Pith},
title = {Pith review of: Hybrid static potentials and gluelumps on $N_f=3+1$ ensembles},
year = {2026},
howpublished = {\url{https://pith.science/paper/PG7P2M55}},
note = {Machine review of arXiv:2501.15670}
}
abstract
QCD permits the existence of hybrid mesons that are made up of both quarks and gluons, including exotic states, i.e., quantum numbers prohibited for pure quark-antiquark states, with possible candidates found in experiments. We present static hybrid potentials measured via Laplace trial states together with static-light meson thresholds on $N_f=3+1$ dynamical fermion ensembles with 420 MeV pions. Furthermore, we measure corresponding gluelump masses which refer to the $R\rightarrow0$ limit of the hybrid potentials and are essential input parameters for effective models to describe hybrid mesons.
Figures
Reference graph
Works this paper leans on
- [13]
-
[6]
S. Capitani, O. Philipsen, C. Reisinger, C. Riehl and M. Wagner,Precision computation of hybrid static potentials in SU(3) lattice gauge theory,Phys. Rev. D99 (2019) 034502 [1811.11046]
arXiv 2019
-
[1]
Dudek,Phenomenology of exotic hadrons - hybrid mesons and pentaquarks, Ph.D
J.J. Dudek,Phenomenology of exotic hadrons - hybrid mesons and pentaquarks, Ph.D. thesis, Oxford U., Theor. Phys., 2004
work page 2004
-
[2]
E.S. Swanson,XYZ States: Theory Overview,AIP Conf. Proc.1735(2016) 020013 [1512.04853]
work page Pith review arXiv 2016
-
[3]
N. Brambilla, S. Eidelman, C. Hanhart, A. Nefediev, C.-P. Shen, C.E. Thomas et al.,The 𝑋𝑌𝑍 states: experimental and theoretical status and perspectives, Phys. Rept.873 (2020) 1 [1907.07583]
arXiv 2020
-
[4]
N. Brambilla, W.K. Lai, A. Mohapatra and A. Vairo,Heavy hybrid decays to quarkonia, Phys. Rev. D107(2023) 054034 [2212.09187]
arXiv 2023
-
[5]
M. Berwein, N. Brambilla, A. Mohapatra and A. Vairo,Hybrids, tetraquarks, pentaquarks, doubly heavy baryons, and quarkonia in Born-Oppenheimer effective theory, Phys. Rev. D 110 (2024) 094040 [2408.04719]
arXiv 2024
-
[7]
N.A. Campbell, I.H. Jorysz and C. Michael,The Adjoint Source Potential in SU(3) Lattice Gauge Theory,Phys. Lett. B167 (1986) 91
work page 1986
Show all 42 references
-
[8]
Jorysz and C
I.H. Jorysz and C. Michael,The Field Configurations of a Static Adjoint Source in SU(2) Lattice Gauge Theory, Nucl. Phys. B302 (1988) 448
1988
-
[9]
UKQCDcollaboration, Hadrons with a heavy color adjoint particle,Phys. Rev. D59(1999) 094509 [hep-lat/9811010]
1999 arXiv
-
[10]
Marsh and R
K. Marsh and R. Lewis,A lattice QCD study of generalized gluelumps, Phys. Rev. D89 (2014) 014502 [1309.1627]
2014 arXiv
-
[11]
J. Herr, C. Schlosser and M. Wagner,Gluelump masses and mass splittings from SU(3) lattice gauge theory,Phys. Rev. D109(2024) 034516 [2306.09902]
2024 arXiv
-
[12]
Herr,Computation of gluelump masses for hybrid static potentials in SU(3) lattice gauge theory using the multilevel algorithm, MSc thesis at Goethe University Frankfurt(2022)
J. Herr,Computation of gluelump masses for hybrid static potentials in SU(3) lattice gauge theory using the multilevel algorithm, MSc thesis at Goethe University Frankfurt(2022)
2022
-
[14]
Höllwieser, F
R. Höllwieser, F. Knechtli, T. Korzec, M. Peardon and J.A. Urrea-Niño,Constructing static quark-antiquark creation operators from Laplacian eigenmodes, Phys. Rev. D107 (2023) 034511 [2212.08485]. 8 Hybrid static potentials and gluelumps on𝑁𝑓 = 3+ 1 ensembles Roman Höllwieser
2023 arXiv
-
[15]
Höllwieser, F
R. Höllwieser, F. Knechtli and M. Peardon,The static energy of a quark-antiquark pair from Laplacian eigenmodes,PoS LATTICE2022(2023) 060 [2209.01239]
2023 arXiv
-
[16]
Alasiri, E
F. Alasiri, E. Braaten and A. Mohapatra,Born-Oppenheimer Potentials for𝑆𝑈(3) Gauge Theory, 2406.05123
-
[17]
Blossier, M
B. Blossier, M. Della Morte, G. von Hippel, T. Mendes and R. Sommer,On the generalized eigenvalue method for energies and matrix elements in lattice field theory, JHEP 04(2009) 094 [0902.1265]
2009 arXiv
-
[18]
Berg and A
B. Berg and A. Billoire,Glueball Spectroscopy in Four-Dimensional SU(3) Lattice Gauge Theory. 1., Nucl. Phys. B221 (1983) 109
1983
-
[19]
Georgi,Lie Algebras in Particle Physics: From Isospin to Unified Theories, CRC Press (2018)
H. Georgi,Lie Algebras in Particle Physics: From Isospin to Unified Theories, CRC Press (2018)
2018
-
[20]
P. R. Bunker and P. Jensen,Molecular symmetry and spectroscopy,NRC Press, Ottawa, ON, Canada, 2 edition(1998)
1998
-
[21]
Bali,QCD forces and heavy quark bound states,Phys
G.S. Bali,QCD forces and heavy quark bound states,Phys. Rept.343 (2001) 1 [hep-ph/0001312]
2001 arXiv
-
[22]
Bali and A
G.S. Bali and A. Pineda,QCD phenomenology of static sources and gluonic excitations at short distances,Phys. Rev. D69(2004) 094001 [hep-ph/0310130]
2004 arXiv
-
[23]
Bauer, G.S
C. Bauer, G.S. Bali and A. Pineda,Compelling Evidence of Renormalons in QCD from High Order Perturbative Expansions, Phys. Rev. Lett.108 (2012) 242002 [1111.3946]
2012 arXiv
-
[24]
G.S. Bali, C. Bauer, A. Pineda and C. Torrero,Perturbative expansion of the energy of static sources at large orders in four-dimensional SU(3) gauge theory,Phys. Rev. D87(2013) 094517 [1303.3279]
2013 arXiv
-
[25]
G.S. Bali, C. Bauer and A. Pineda,The static quark self-energy at𝑂(𝛼20) in perturbation theory,PoS LATTICE2013(2014) 371 [1311.0114]
2014 arXiv
-
[26]
ALPHAcollaboration, Scale setting for𝑁𝑓 = 3+ 1 QCD, Eur. Phys. J. C80(2020) 349 [2002.02866]
2020 arXiv
-
[27]
Lüscher and S
M. Lüscher and S. Schaefer,Lattice QCD with open boundary conditions and twisted-mass reweighting,Comput. Phys. Commun.184 (2013) 519 [1206.2809]
2013 arXiv
-
[28]
ALPHAcollaboration, Symanzik improvement with dynamical charm: a 3+1 scheme for Wilson quarks, JHEP 06(2018) 025 [1805.01661]
2018 arXiv
-
[29]
Sommer,A New way to set the energy scale in lattice gauge theories and its applications to the static force and𝛼𝑠 in SU(2) Yang-Mills theory,Nucl
R. Sommer,A New way to set the energy scale in lattice gauge theories and its applications to the static force and𝛼𝑠 in SU(2) Yang-Mills theory,Nucl. Phys. B411 (1994) 839 [hep-lat/9310022]. 9 Hybrid static potentials and gluelumps on𝑁𝑓 = 3+ 1 ensembles Roman Höllwieser
1994 arXiv
-
[30]
Lüscher,Properties and uses of the Wilson flow in lattice QCD, JHEP 08(2010) 071 [1006.4518]
M. Lüscher,Properties and uses of the Wilson flow in lattice QCD, JHEP 08(2010) 071 [1006.4518]
2010 arXiv
-
[31]
APE collaboration,Glueball Masses and String Tension in Lattice QCD,Phys. Lett. B192 (1987) 163
1987
-
[32]
Hasenfratz and F
A. Hasenfratz and F. Knechtli,Flavor symmetry and the static potential with hypercubic blocking,Phys. Rev. D64(2001) 034504 [hep-lat/0103029]
2001 arXiv
-
[33]
ALPHAcollaboration, Static quarks with improved statistical precision, Nucl. Phys. B Proc. Suppl. 129(2004) 346 [hep-lat/0309080]
2004 arXiv
-
[34]
Della Morte, A
M. Della Morte, A. Shindler and R. Sommer,On lattice actions for static quarks, JHEP 08 (2005) 051 [hep-lat/0506008]
2005 arXiv
-
[35]
Grimbach, D
A. Grimbach, D. Guazzini, F. Knechtli and F. Palombi,O(a) improvement of the HYP static axial and vector currents at one-loop order of perturbation theory, JHEP 03(2008) 039 [0802.0862]
2008 arXiv
-
[36]
Donnellan, F
M. Donnellan, F. Knechtli, B. Leder and R. Sommer,Determination of the Static Potential with Dynamical Fermions, Nucl. Phys. B849 (2011) 45 [1012.3037]
2011 arXiv
-
[37]
ALPHAcollaboration, Monte Carlo errors with less errors, Comput. Phys. Commun.156 (2004) 143 [hep-lat/0306017]
2004 arXiv
-
[38]
ALPHAcollaboration, Critical slowing down and error analysis in lattice QCD simulations, Nucl. Phys. B845 (2011) 93 [1009.5228]
2011 arXiv
-
[39]
Joswig, S
F. Joswig, S. Kuberski, J.T. Kuhlmann and J. Neuendorf,pyerrors: a python framework for error analysis of Monte Carlo data, 2209.14371
-
[40]
Ramos,Automatic differentiation for error analysis, PoS TOOLS2020(2021) 045 [2012.11183]
A. Ramos,Automatic differentiation for error analysis, PoS TOOLS2020(2021) 045 [2012.11183]
2021 arXiv
-
[41]
TXL, T(X)L collaboration,Static potentials and glueball masses from QCD simulations with Wilson sea quarks, Phys. Rev. D62(2000) 054503 [hep-lat/0003012]
2000 arXiv
-
[42]
Struckmeier, R
L. Struckmeier, R. Höllwieser, F. Knechtli, T. Korzec, M. Peardon and J.A. Urrea-Niño, Static-light meson spectroscopy with optimal distillation profiles,PoS LATTICE2024 (2025) 104 [2501.12863]. 10
2025 arXiv
Reviewed August 10, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.