Semirelativistic potential model for glueball states
  The masses of two-gluon glueballs are studied with a semirelativistic potential model whose interaction is a scalar linear confinement supplemented by a one-gluon exchange mechanism. The gluon is massless but the leading corrections of the dominant part of the Hamiltonian are expressed in terms of a state dependent constituent gluon mass. The Hamiltonian depends only on 3 parameters: the strong coupling constant, the string tension, and a gluon size which removes all singularities in the leading corrections of the potential. Accurate numerical calculations are performed with a Lagrange mesh method. The masses predicted are in rather good agreement with lattice results and with some experimental glueball candidates. 
F. Brau, C. Semay, B. Silvestre-Brac