Spectrum of scalar and pseudoscalar glueballs from functional methods

We provide results for the spectrum of scalar and pseudoscalar glueballs in pure Yang–Mills theory using a parameter-free fully self-contained truncation of Dyson–Schwinger and Bethe–Salpeter equations. The only input, the scale, is fixed by comparison with lattice calculations. We obtain ground state masses of [Formula: see text] and [Formula: see text] for the scalar and pseudoscalar glueballs, respectively, and [Formula: see text] and [Formula: see text] for the corresponding first excited states. This is in very good quantitative agreement with available lattice results. Furthermore, we predict masses for the second excited states at [Formula: see text] and [Formula: see text] . The quality of the results hinges crucially on the self-consistency of the employed input. The masses are independent of a specific choice for the infrared behavior of the ghost propagator providing further evidence that this only reflects a nonperturbative gauge completion.