Lattice Calculation of Glueball Matrix Elements

Matrix elements of the form $ <0| Tr \; g^{2}GG |G> $ are calculated using the lattice QCD Monte Carlo method. Here, $|G>$ is a glueball state with quantum numbers $ 0^{++}$, $ 2^{++}$, $ 0^{-+}$ and $G$ is the gluon field strength operator. The matrix elements are obtained from the hybrid correlation functions of the fuzzy and plaquette operators performed on the $12^{4}$ and $14^{4}$ lattices at $\beta = 5.9 $ and $5.96$ respectively. These matrix elements are compared with those from the QCD sum rules and the tensor meson dominance model. They are the non-perturbative matrix elements needed in the calculation of the partial widths of $J/\Psi$ radiative decays into glueballs.