Abstract: We study in detail the impact of dynamical quarks on the gluon mass generation mechanism, in the Landau gauge, for the case of a small number of quark families. As in earlier considerations, we assume that the main bulk of the unquenching corrections to the gluon propagator originates from the fully dressed quark-loop diagram. The nonperturbative evaluation of this diagram provides the key relation that expresses the unquenched gluon propagator as a deviation from its quenched counterpart. This relation is subsequently coupled to the integral equation that controls the momentum evolution of the effective gluon mass, which contains a single adjustable parameter; this constitutes a major improvement compared to the analysis presented in Aguilar et al. [Phys. Rev. D 86, 014032 (2012)], where the behavior of the gluon propagator in the deep infrared was estimated through numerical extrapolation. The resulting nonlinear system is then treated numerically, yielding unique solutions for the modified gluon mass and the quenched gluon propagator, which fully confirms the picture put forth recently in several continuum and lattice studies. In particular, an infrared finite gluon propagator emerges, whose saturation point is considerably suppressed, due to a corresponding increase in the value of the gluon mass. This characteristic feature becomes more pronounced as the number of active quark families increases, and can be deduced from the infrared structure of the kernel entering in the gluon mass equation.

Abstract: We present a detailed study of the subtle interplay transpiring at the level of two integral equations that are instrumental for the dynamical generation of a gluon mass in pure Yang-Mills theories. The main novelty is the joint treatment of the Schwinger-Dyson equation governing the infrared behavior of the gluon propagator and of the integral equation that controls the formation of massless bound-state excitations, whose inclusion is instrumental for obtaining massive solutions from the former equation. The self-consistency of the entire approach imposes the requirement of using a single value for the gauge coupling entering in the two key equations; its fulfilment depends crucially on the details of the three-gluon vertex, which contributes to both of them, but with different weight. In particular, the characteristic suppression of this vertex at intermediate and low energies enables the convergence of the iteration procedure to a single gauge coupling, whose value is reasonably close to that extracted from related lattice simulations.

Abstract: We present a generalized theoretical framework for dealing with the important issue of dynamical mass generation in Yang-Mills theories, and, in particular, with the infrared finiteness of the gluon propagators, observed in a multitude of recent lattice simulations. Our analysis is manifestly gauge invariant, in the sense that it preserves the transversality of the gluon self-energy, and gauge independent, given that the conclusions do not depend on the choice of the gauge-fixing parameter within the linear covariant gauges. The central construction relies crucially on the subtle interplay between the Abelian Ward identities satisfied by the nonperturbative vertices and a special integral identity that enforces a vast number of “seagull cancellations” among the one-and two-loop dressed diagrams of the gluon Schwinger-Dyson equation. The key result of these considerations is that the gluon propagator remains rigorously massless, provided that the vertices do not contain (dynamical) massless poles. When such poles are incorporated into the vertices, under the pivotal requirement of respecting the gauge symmetry of the theory, the terms comprising the Ward identities conspire in such a way as to still enforce the total annihilation of all quadratic divergences, inducing, at the same time, residual contributions that account for the saturation of gluon propagators in the deep infrared.