An Independent Computational Reproduction of SU(2) Lattice Glueball Scaling

This technical note presents an independent computational reproduction of the established scalar glueball mass scaling in four-dimensional pure SU(2) Yang–Mills theory. Using an independently implemented C++17/OpenMP lattice Monte Carlo simulation with the Wilson gauge action, five lattice spacings (β = 2.3–2.7) are investigated. A continuum extrapolation yields m(0⁺⁺)/√σ = 3.96 ± 0.29, consistent within 0.6 standard deviations with the published benchmark of 3.781 ± 0.023 (Athenodorou and Teper, 2021). The analysis includes variational glueball correlators, Wilson-loop determinations of the string tension, Sommer-scale extraction, jackknife uncertainty estimation, autocorrelation diagnostics, 32 fit-window variations, and independent-seed and finite-volume checks. An additional plaquette-susceptibility scan examines the strong–weak coupling crossover on two lattice volumes. The principal contribution is a transparent and reproducible numerical implementation of established lattice gauge theory results, rather than a new physical prediction or a mathematical proof of the Yang–Mills mass gap. The accompanying research archive contains C++ simulation sources, Python analysis scripts, numerical results, validation data, and the raw Monte Carlo measurements needed to reproduce the reported calculations.