pith. machine review for the scientific record. sign in

arxiv: 2604.04803 · v1 · submitted 2026-04-06 · ✦ hep-ph · hep-lat· hep-th· nucl-th

Recognition: 2 theorem links

· Lean Theorem

Glueballs, Constituent Gluons and Instantons

Edward Shuryak , Ismail Zahed
Authors on Pith no claims yet

Pith reviewed 2026-05-10 19:15 UTC · model grok-4.3

classification ✦ hep-ph hep-lathep-thnucl-th
keywords glueballsconstituent gluonsinstantonsYang-Mills theoryquenched latticescalar glueballtensor glueballRegge trajectories
0
0 comments X

The pith

The scalar glueball is a compact two-gluon state whose radius equals the instanton size of about one-third of a fermi.

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

The paper constructs a two-gluon model for the lowest glueball states in pure Yang-Mills theory and calibrates it directly to quenched lattice results. Instantons supply a dynamical gluon mass together with central and tensor forces, while adjoint Casimir scaling sets the long-range confinement and a Coulomb term handles short distances. The resulting calculation yields an unusually small scalar state whose size tracks the instanton scale, whereas the tensor state spreads out because of the centrifugal barrier. A reader would care because the approach links a well-known nonperturbative object, the instanton, to concrete sizes and level orderings that lattice simulations already report.

Core claim

We present a constituent two-gluon description of the lowest-lying glueball states in pure Yang-Mills theory, calibrated against quenched lattice results. The framework incorporates an instanton-induced dynamical gluon mass, Casimir-scaled adjoint confinement, the short-distance adjoint Coulomb interaction, and instanton-induced central and tensor forces. The scalar 0++ glueball is found to be exceptionally compact, with a radius of order the instanton size, ρ ∼ 1/3 fm, consistent with lattice indications. By contrast, the tensor 2++ state remains spatially extended due to the centrifugal barrier. We also discuss the role of S-D mixing. A semiclassical analysis further supports Regge behav

What carries the argument

The two-gluon wave function solved with an instanton-generated dynamical mass, adjoint Casimir confinement, short-range Coulomb force, and additional instanton central and tensor potentials.

If this is right

  • The tensor glueball stays spatially extended because the centrifugal barrier prevents it from sampling the short-range instanton attraction.
  • S-D mixing modifies the wave functions of both states.
  • Excited glueballs exhibit Regge trajectories whose slopes match those seen on the lattice.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Instantons, rather than the long-range confining string, would dominate the binding energy of the ground-state glueball.
  • The same short-distance instanton forces could influence glueball-meson mixing once light quarks are restored.
  • Higher radial excitations should show a characteristic pattern of sizes that future finer lattice simulations could check.

Load-bearing premise

The two-gluon picture calibrated solely to quenched lattice data remains accurate even when higher Fock components or dynamical quarks are present.

What would settle it

A lattice computation that measures the scalar glueball radius to be substantially larger than one-third of a fermi, or that shows the size is insensitive to the instanton density, would falsify the central claim.

Figures

Figures reproduced from arXiv: 2604.04803 by Edward Shuryak, Ismail Zahed.

Figure 1
Figure 1. Figure 1: FIG. 1. Repulsive, neutral and attractive channels in [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Scalar [PITH_FULL_IMAGE:figures/full_fig_p002_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Lattice data for glueball masses, normalized to [PITH_FULL_IMAGE:figures/full_fig_p003_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: A more detailed discussion, as well as the relation of this potential shape to our Monte-Carlo calculations with a single instanton or ¯II molecule using numerically generated Wilson lines, is given in Appendix C. C. The short-distance nonperturbative interactions Let us start with simple estimates of instanton￾induced interactions. Naively representing the gauge field as semiclassical O(1/g) plus quantum … view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. The colored points show the calculated ener [PITH_FULL_IMAGE:figures/full_fig_p006_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Wave functions (unnormalized) [PITH_FULL_IMAGE:figures/full_fig_p006_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7 [PITH_FULL_IMAGE:figures/full_fig_p013_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Glueball mass spectrum organized by [PITH_FULL_IMAGE:figures/full_fig_p014_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. Static potentials [PITH_FULL_IMAGE:figures/full_fig_p016_9.png] view at source ↗
read the original abstract

We present a constituent two-gluon description of the lowest-lying glueball states in pure Yang--Mills theory, calibrated against quenched lattice results. The framework incorporates an instanton-induced dynamical gluon mass, Casimir-scaled adjoint confinement, the short-distance adjoint Coulomb interaction, and instanton-induced central and tensor forces. The scalar $0^{++}$ glueball is found to be exceptionally compact, with a radius of order the instanton size, $\rho \sim \frac 13\,\mathrm{fm}$, consistent with lattice indications. By contrast, the tensor $2^{++}$ state remains spatially extended due to the centrifugal barrier. We also discuss the role of $S$-$D$ mixing. A semiclassical analysis further supports Regge behavior for excited states, in agreement with lattice results.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 3 minor

Summary. The manuscript develops a non-relativistic constituent two-gluon model for the lowest-lying glueball states in pure Yang-Mills theory. It incorporates an instanton-induced dynamical gluon mass, Casimir-scaled adjoint linear confinement, short-distance adjoint Coulomb interaction, and instanton-induced central and tensor forces. Parameters are calibrated to quenched lattice glueball masses. The model yields a compact 0++ ground state with rms radius of order the instanton size (∼1/3 fm), while the 2++ state is spatially extended due to the centrifugal barrier; S-D mixing and semiclassical Regge trajectories for excited states are also discussed.

Significance. If the two-gluon truncation and calibration hold, the work supplies a transparent phenomenological link between instanton physics and lattice glueball structure, with the explicit radius prediction and Regge support as concrete outputs. The approach is valuable for interpreting lattice indications of compactness and for guiding phenomenology, though its predictive reach is constrained by the fitting procedure.

major comments (3)
  1. [§3] §3: The radial Schrödinger equation for the two-gluon system and the explicit functional form of the full potential (instanton-induced mass term, Casimir-scaled linear term, Coulomb term, and central/tensor instanton forces) are not written down. Without these, the reported rms radius of ∼1/3 fm for the 0++ state cannot be independently reproduced or checked for numerical stability.
  2. [§2–3] §2–3: The gluon mass, adjoint string tension, and instanton density are adjusted to reproduce the quenched lattice spectrum; the subsequent claim that the 0++ radius matches lattice indications therefore rests on the same data used for calibration, with no reported propagation of fit uncertainties or cross-validation against independent observables.
  3. [§4] §4: The two-gluon Fock-space truncation is adopted without any estimate of mixing amplitudes with three-gluon or instanton-molecule configurations. Given that the extracted radius is comparable to the instanton size, such admixtures could enlarge the spatial extent and must be quantified to support the central compactness claim.
minor comments (3)
  1. The abstract states consistency with lattice indications for the radius but cites no specific lattice reference or numerical value for comparison in the main text.
  2. [§3] Wave-function plots (presumably in §3) lack sensitivity bands showing radius variation under small changes in the fitted parameters.
  3. The definition of the tensor-force strength parameter appears without an equation label, complicating cross-reference.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the constructive and detailed report. We address each major comment below and describe the revisions planned for the manuscript.

read point-by-point responses

  1. Referee: [§3] §3: The radial Schrödinger equation for the two-gluon system and the explicit functional form of the full potential (instanton-induced mass term, Casimir-scaled linear term, Coulomb term, and central/tensor instanton forces) are not written down. Without these, the reported rms radius of ∼1/3 fm for the 0++ state cannot be independently reproduced or checked for numerical stability.

    Authors: We agree that the explicit forms are required for reproducibility. In the revised manuscript we will add the radial Schrödinger equation for the two-gluon system together with the complete analytic expression for the potential, including the instanton-induced mass term, the Casimir-scaled adjoint linear term, the adjoint Coulomb interaction, and the central and tensor forces generated by instantons. This will permit independent numerical verification of the rms radius. revision: yes

  2. Referee: [§2–3] §2–3: The gluon mass, adjoint string tension, and instanton density are adjusted to reproduce the quenched lattice spectrum; the subsequent claim that the 0++ radius matches lattice indications therefore rests on the same data used for calibration, with no reported propagation of fit uncertainties or cross-validation against independent observables.

    Authors: The parameters are fitted to the lattice masses, as stated in the manuscript. The radius is nevertheless a derived observable obtained by solving the Schrödinger equation once the potential is fixed; its small value is directly tied to the instanton size scale. We will add to the revision a brief sensitivity analysis showing how the rms radius varies when the fitted parameters are shifted within their uncertainties. Cross-validation against independent quantities such as decay widths lies beyond the present scope but will be noted as a direction for future work. revision: partial

  3. Referee: [§4] §4: The two-gluon Fock-space truncation is adopted without any estimate of mixing amplitudes with three-gluon or instanton-molecule configurations. Given that the extracted radius is comparable to the instanton size, such admixtures could enlarge the spatial extent and must be quantified to support the central compactness claim.

    Authors: We acknowledge the limitation of the two-gluon truncation. A quantitative estimate of mixing amplitudes with three-gluon or instanton-molecule states would require a substantially enlarged Fock space and additional model ingredients that are outside the scope of this work. In the revision we will insert a short discussion of this approximation, emphasizing that the leading two-gluon component already produces a compact state whose size is set by the instanton scale and is consistent with existing lattice indications. Further quantification is left for future study. revision: partial

Circularity Check

0 steps flagged

No significant circularity: radius derived from wavefunction after mass calibration, not reduced to inputs by construction

full rationale

The paper calibrates a two-gluon Schrödinger model (with instanton mass, Casimir confinement, Coulomb and tensor forces) to quenched lattice glueball masses. The 0++ radius is then obtained as an output from the resulting ground-state wavefunction and compared to separate lattice indications of size. This is a standard phenomenological derivation, not a self-definitional loop or a fitted quantity renamed as prediction. No equations reduce the radius to the input masses by construction, and no load-bearing self-citations or ansatz smuggling are quoted. The chain remains self-contained against the external lattice benchmarks.

Axiom & Free-Parameter Ledger

3 free parameters · 2 axioms · 1 invented entities

The model rests on several fitted quantities and domain assumptions drawn from lattice QCD and instanton liquid phenomenology; without the full text the exact count is uncertain, but the abstract explicitly invokes calibration to quenched lattice data and instanton-induced mass.

free parameters (3)
  • instanton-induced gluon mass
    Dynamical mass generated by instantons; value chosen to match lattice glueball masses.
  • adjoint string tension
    Casimir-scaled confinement strength fitted to lattice.
  • instanton density and size
    Parameters controlling the strength of central and tensor forces; calibrated to reproduce lattice radii and splittings.
axioms (2)
  • domain assumption Quenched lattice results provide a reliable benchmark for pure Yang-Mills glueball spectrum and sizes.
    Invoked when the model is calibrated against and compared to lattice data.
  • domain assumption Two-gluon Fock space is sufficient for the lowest glueball states.
    Underlying the constituent description; higher components neglected.
invented entities (1)
  • constituent gluon with instanton-induced mass no independent evidence
    purpose: Effective massive degree of freedom for variational bound-state calculation
    Postulated to simplify the non-perturbative dynamics; independent evidence would be a direct lattice measurement of the gluon propagator pole.

pith-pipeline@v0.9.0 · 5433 in / 1688 out tokens · 46120 ms · 2026-05-10T19:15:51.044945+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Hybrid hadrons at rest and on the light front

    hep-ph 2026-04 unverdicted novelty 6.0

    A constituent-gluon Born-Oppenheimer model is used to derive light-front wave functions and gluon PDFs for hybrid hadrons ccg and qqqg.

Reference graph

Works this paper leans on

35 extracted references · 23 canonical work pages · cited by 1 Pith paper

  1. [1]

    Normalized trial wavefunctions For simplicity, we will use variational estimates of matrix elements using simplified wave functions. TheS-wave Gaussian trial function is ψS(r) = β2 S π 3/4 e−β2 S r2/2,(E1) with radial part RS(r) = 2β3 S√π 1/2 e−β2 S r2/2.(E2) TheD-wave trial function is ψD,m(r) =N D r2e−β2 Dr2/2Y2m(ˆr),(E3) with normalization fixed by 1 =...

  2. [2]

    Kinetic and confining expectation values Using standard Gaussian integrals, ⟨p2⟩S = 3 2 β2 S,⟨p 2⟩D = 7 2 β2 D.(E7) The linear confinement expectation values are ⟨r⟩S = 2√π 1 βS ,⟨r⟩ D = 16 3√π 1 βD .(E8) The corresponding rms radii are ⟨r2⟩S = 3 2β2 S ,⟨r 2⟩D = 7 2β2 D .(E9)

  3. [3]

    Coulomb expectation values For the Coulomb potentialVC(r) =−3α eff s /r, one finds 1 r S = 2βS√π ,(E10) and for theD-wave, 1 r D = 8 3√π βD.(E11) ThustheCoulombtermisparametricallysuppressed in the tensor channel by angular momentum

  4. [4]

    Instanton central attraction The scalar instanton term yields D e−r2/ρ2E S = β2 S β2 S +ρ −2 3/2 .(E12) For theD-wave, D e−r2/ρ2E D = β2 D β2 D +ρ −2 7/2 ,(E13) showing that instanton attraction is strongly sup- pressed for higher partial waves

  5. [5]

    Spin-spin contact term With the Gaussian regulator δ(3) Λ (r) = Λ2 π 3/2 e−Λ2r2 ,(E14) the contact expectation values are ⟨δ(3) Λ ⟩S = Λ2 Λ2 +β 2 S 3/2 β2 S π 3/2 ,(E15) and ⟨δ(3) Λ ⟩D = Λ2 Λ2 +β 2 D 7/2 β2 D π 3/2 15 8 .(E16) 20

  6. [6]

    Tensor matrix elements The tensor operator enters through VT (r) = CT (η) m2g(η) 1−e −r2/ρ2 r3 S12.(E17) The reduced angular matrix elements in the {5S2, 5D2}basis are ⟨5S2|S12|5S2⟩= 0,⟨ 5D2|S12|5D2⟩=− 2 7 ,⟨ 5S2|S12|5D2⟩= r 8 7 .(E18) The radial mixing integral is ISD = Z ∞ 0 dr RS(r)RD(r) 1−e −r2/ρ2 r .(E19) which yields ISD =C 1 (β2 S +β 2 D)1/2 − 1 (β...

  7. [7]

    Explicit 5S2-5D2 tensor Hamiltonian Inthecoupledbasis{| 5S2⟩,| 5D2⟩}thetensorglue- ball Hamiltonian takes the explicit2×2form H2++(η) =  HSS(η)H SD(η) HSD(η)H DD(η)   ,(E24) with the diagonalS-wave entry HSS(η) =E (0) S (βS;η) + CSS(η) m2g(η) D δ(3) Λ E S ,(E25) where⟨δ (3) Λ ⟩S is given in Eq. (E15). The diagonalD-wave entry reads HDD(η) =E (0) D (βD...

  8. [8]

    C. J. Morningstar and M. J. Peardon, Phys. Rev. D 60, 034509 (1999), arXiv:hep-lat/9901004

    work page arXiv 1999

  9. [9]

    H. B. Meyer and M. J. Teper, Nucl. Phys. B Proc. Suppl.129, 200 (2004), arXiv:hep-lat/0308035

    work page arXiv 2004

  10. [10]

    Chenet al., Glueball spectrum and matrix elements on anisotropic lattices, Phys

    Y. Chenet al., Phys. Rev. D73, 014516 (2006), arXiv:hep-lat/0510074

    work page arXiv 2006

  11. [11]

    Athenodorou and M

    A. Athenodorou and M. Teper, JHEP11, 172 (2020), arXiv:2007.06422 [hep-lat]

    work page arXiv 2020

  12. [12]

    C. A. Meyer and Y. Van Haarlem, Phys. Rev. C82, 025208 (2010), arXiv:1004.5516 [nucl-ex]

    work page arXiv 2010

  13. [13]

    Gravitational form fac- tors of glueballs in yang-mills theory,

    R. Abbott, D. C. Hackett, D. A. Pefkou, F. Romero- López, and P. Shanahan, “Gravitational form fac- tors of glueballs in yang-mills theory,” (2025), arXiv:2410.02706 [hep-lat]

    work page arXiv 2025

  14. [14]

    Morningstar, (2025), arXiv:2502.02547 [hep-ph]

    C. Morningstar, (2025), arXiv:2502.02547 [hep-ph]

    work page arXiv 2025

  15. [15]

    Schäfer and E

    T. Schäfer and E. V. Shuryak, Phys. Rev. Lett.75, 1707 (1995), arXiv:hep-ph/9410372

    work page arXiv 1995

  16. [16]

    Instantons in QCD

    T. Schafer and E. V. Shuryak, Rev. Mod. Phys.70, 323 (1998), arXiv:hep-ph/9610451

    work page Pith review arXiv 1998

  17. [17]

    Exploring improved holographic theories for QCD: Part I

    U. Gursoy and E. Kiritsis, JHEP02, 032 (2008), arXiv:0707.1324 [hep-th]

    work page Pith review arXiv 2008

  18. [18]

    Iatrakis, A

    I. Iatrakis, A. Ramamurti, and E. Shuryak, Phys. Rev. D92, 014011 (2015), arXiv:1503.04759 [hep- ph]

    work page arXiv 2015

  19. [19]

    Jiang, W

    X. Jiang, W. Sun, F. Chen, Y. Chen, M. Gong, Z. Liu, and R. Zhang, Phys. Rev. D107, 094510 (2023), arXiv:2205.12541 [hep-lat]. 27

    work page arXiv 2023

  20. [20]

    Alkofer and L

    R. Alkofer and L. von Smekal, Phys. Rept.353, 281 (2001), arXiv:hep-ph/0007355

    work page arXiv 2001

  21. [21]

    C. S. Fischer, J. Phys. G39, 045004 (2012), arXiv:1111.0603 [hep-ph]

    work page arXiv 2012

  22. [22]

    A. Holl, A. Krassnigg, and R. Alkofer, Phys. Rev. D91, 014031 (2015), arXiv:1410.1469 [hep-ph]

    work page arXiv 2015

  23. [24]

    Isgur and J

    N. Isgur and J. E. Paton, Phys. Rev. D31, 2910 (1985)

    1985

  24. [25]

    Johnson and C

    K. Johnson and C. B. Thorn, Phys. Rev. D13, 1934 (1976)

    1934

  25. [26]

    Brunner, D

    F. Brunner, D. Parganlija, and A. Rebhan, Phys. Rev. D91, 106002 (2015), arXiv:1501.07906 [hep- ph]

    work page arXiv 2015

  26. [27]

    M. M. Musakhanov and O. Egamberdiev, Phys. Lett. B781, 303 (2018), arXiv:1706.06270 [hep-ph]

    work page arXiv 2018

  27. [29]

    Mathieu, N

    V. Mathieu, N. Kochelev, and V. Vento, Int. J. Mod. Phys. E18, 1 (2009), arXiv:0810.4453 [hep- ph]

    work page arXiv 2009

  28. [30]

    H.B.Meyer,Glueball regge trajectories,Otherthesis (2004), arXiv:hep-lat/0508002

    work page arXiv 2004

  29. [31]

    J. M. Cornwall and A. Soni, Phys. Lett. B120, 431 (1983)

    1983

  30. [32]

    W.-Y. Liu, E. Shuryak, and I. Zahed, Phys. Rev. D 110, 054005 (2024), arXiv:2404.03047 [hep-ph]

    work page arXiv 2024

  31. [33]

    Shuryak and I

    E. Shuryak and I. Zahed, Phys. Rev. D107, 034023 (2023), arXiv:2110.15927 [hep-ph]

    work page arXiv 2023

  32. [34]

    Abbott, D.C

    R. Abbott, D. C. Hackett, D. A. Pefkou, F. Romero-López, and P. E. Shanahan, (2025), arXiv:2508.21821 [hep-lat]

    work page arXiv 2025

  33. [35]

    E. V. Shuryak, Nucl. Phys. B203, 93 (1982)

    1982

  34. [36]

    Shuryak and I

    E. Shuryak and I. Zahed, Phys. Rev. D107, 034025 (2023), arXiv:2112.15586 [hep-ph]

    work page arXiv 2023

  35. [37]

    S. J. Brodsky, T. Huang, and G. P. Lepage, Conf. Proc. C810816, 143 (1981)

    1981