While we have several complementary models of confinement, some of which are phenomenologically appealing, we do not have the ability to calculate analytically even simple aspects of confinement, let alone have a framework to eventually prove confinement. The problem we are facing is to evolve the theory from the perturbative regime to the long distance confining regime. This is generally achieved by renormalization group transformations. With the gradient flow we now have a technique to address the problem from first principles. The primary focus is on the running coupling αS (μ), from which confinement can be concluded alone. A central point is that the gluon condensate is scale invariant, which reflects its self-similar behavior across different scales. Building on that, we derive αS (μ) ΛS²/μ², which evolves to the infrared fixed point 1/αS = 0 in accordance with infrared slavery. This shows that the only important factor is the presence of homogeneous vacuum fields, represented by condensates, which is a universal feature that QCD shares with many other models. The analytical statements are supported by numerical simulations.
G. Schierholz (Fri,) studied this question.