Correlation functions of Landau gauge Yang-Mills theory

Correlation functions of Yang-Mills theory in the Landau gauge are calculated from their equations of motion. The employed setup is completely parameter free and leads, within errors, to good quantitative agreement with corresponding lattice results for the ghost and gluon propagators as well as the ghost-gluon and three-gluon vertices. Also, the four-gluon vertex is calculated. The present setup allows for the first time for a unique subtraction of quadratic divergences in the gluon propagator Dyson-Schwinger equation. Thus, there are no ambiguities which can arise due to the use of models or auxiliary workarounds. In addition, several self-tests of the results are described that allow assessing the truncation error in a self-consistent way. This enables a new perspective on how to identify limitations of the present setup and develop future improvements.