Glueballs from Dyson-Schwinger and Bethe-Salpeter equations

The quenched spectrum of glueballs with positive charge parity is calculated from two-body bound state equations. As input, a self-contained solution for the primitively divergent correlation functions from Dyson-Schwinger equations is used. It only has one parameter to be set which is the physical scale. An important feature of this setup is the consistent construction of the bound state kernels along the same lines as the equations from which the input was obtained. Keeping only the one-particle exchanges, already good agreement with lattice results is obtained. For the tensor glueball, we present first results including two-loop contributions, elevating its calculation to the same level of truncation as for the spin zero glueballs for which such calculations have been done previously.