We report on the first self-consistent solution of the Dyson\char21{}Schwinger equation (DSE) for the three-gluon vertex. Based on earlier results for the propagators which match data from lattice Monte Carlo simulations, we obtain results for the three-gluon vertex that are in very good agreement with available lattice data likewise. Feeding these results back into the propagator DSEs leads to some changes especially in the gluon propagator. These changes allow us to assess previously used models for the three-gluon vertex and to systematically estimate the influence of neglected two-loop diagrams with four-gluon interactions. In the final step, a full iterative solution to the coupled DSEs of pure Landau gauge QCD without quarks is then obtained for the first time in an extended truncation which now dynamically includes the complete set of three-point vertex functions.

We present a self-consistent calculation of the four-gluon vertex of Landau gauge Yang–Mills theory from a truncated Dyson–Schwinger equation. The equation contains the leading diagrams in the ultraviolet and is solved using as the only input results for lower Green functions from previous Dyson–Schwinger calculations that are in good agreement with lattice data. All quantities are therefore fixed and no higher Green functions enter within this truncation. Our self-consistent solution resolves the full momentum dependence of the vertex but is limited to the tree-level tensor structure at the moment. Calculations of selected dressing functions for other tensor structures from this solution are used to exemplify that they are suppressed compared to the tree-level structure except for possible logarithmic enhancements in the deep infrared. Our results furthermore allow one to extract a qualitative fit for the vertex and a running coupling.

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.

We present a more detailed picture of the infrared regime of Landau-gauge Yang-Mills theory. This is done within a novel framework that allows one to take into account the influence of finite scales within an infrared power counting analysis. We find that there are two qualitatively different infrared fixed points of the full system of Dyson-Schwinger equations. The first extends the known scaling solution, where the ghost dynamics is dominant and gluon propagation is strongly suppressed. It features in addition to the strong divergences of gluonic vertex functions in the previously considered uniform scaling limit, when all external momenta tend to zero, also weaker kinematic divergences, when only some of the external momenta vanish. The second solution represents the recently proposed decoupling scenario where the gluons become massive and the ghosts remain bare. In this case we find that none of the vertex functions is enhanced, so that the infrared dynamics is entirely suppressed. Our analysis also provides a strict argument why the Landau-gauge gluon dressing function cannot be infrared divergent.