The Schwinger-Dyson Bethe-Salpeter approach to the bound state problem is applied to the spin zero glueball spectrum. Although a moderately successful description of the two point functions and the glueball spectrum can be obtained, further work is required to reduce sensitivity to model truncations.

Gluon, ghost, and quark propagators are computed in the Schwinger-Dyson formalism. The full set of one-loop and two-loop contributions to the gap equations are evaluated for the first time. A new and efficient method for dealing with quadratic divergences is presented and a consistent renormalization method is proposed. We find that the two-loop contribution to the propagators is subleading, good fits to lattice data are possible, and reliable vertex models are important to obtaining desired ultraviolet, infrared, and gauge-symmetric behavior.