$2^{++}$ Di-gluonium from LSR at higher order

We improve the determination of the mass and coupling of the $2^{++}$ tensor di-gluonium by using relativistic QCD Laplace sum rules (LSR). In so doing, we evaluate the next-to-leading order (NLO) corrections to the perturbative (PT) and $⟨{\alpha }_s{G}^2⟩$ condensate and the lowest order (LO) $⟨G^3⟩$ contributions to the $2^{++}$ di-gluonium two-point correlator. Within a vacuum saturation estimate ($k_G=1$) of the dimension-eight gluon condensates, we obtain: $M_T=3028(287)MeV$ and the renormalization group invariant (RGI) coupling ${\widehat{f}}_T=224(33)MeV$. Assuming that the factorization hypothesis can be violated, we study the effect of the violation factor $k_G$ on the results and obtain: $M_T=3188(337)MeV$ and ${\widehat{f}}_T=245(32)MeV$ for $k_G=(3\pm 2)$. Our estimation does not favour the interpretation of the observed $f_2(2010)$, $f_2(2300)$ and $f_2(2340)$ as pure glueball state.