“…The heavy-quarkonium system is usually described by non-relativistic QCD [1], while the gluonic hadronization has been treated using soft collinear effective theory [2], gluon distribution amplitudes [3], and perturbative QCD [4,5]. 2 Glueballs are a natural consequence of QCD, and predictions of their properties have been made using different approaches, such as potential models [6,7,8], lattice QCD calculations [9,10,11,12], bag models [13,14,15,16], flux-tube models [17], the QCD sum rules [18], the Bethe-Salpeter (B-S) equation [19,20], QCD factorization formalism models [21,22], weakly-bound-state models [23], and a three-dimensional relativistic equation [24]. However, despite intense experimental searches [25,26,27,28,29,30,31], there is no conclusive experimental evidence of their direct observation, although there are strong indications that glueballs contribute to the rich light scalar [32,33,34,35,36,37,38,39,40,41] and tensor [42,…”