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 carry out a systematic study of the renormalization properties of the integral equation that determines the momentum evolution of the effective gluon mass in pure Yang-Mills theory, without quark effects taken into account. A detailed, all-order analysis of the complete kernel appearing in this particular equation, derived in the Landau gauge, reveals that the renormalization procedure may be accomplished through the sole use of ingredients known from the standard perturbative treatment of the theory, with no additional assumptions. However, the subtle interplay of terms operating at the level of the exact equation gets distorted by the approximations usually employed when evaluating the aforementioned kernel. This fact is reflected in the form of the obtained solutions, for which the deviations from the correct behavior are best quantified by resorting to appropriately defined renormalization-group invariant quantities. This analysis, in turn, provides a solid guiding principle for improving the form of the kernel, and furnishes a well-defined criterion for discriminating between various possibilities. Certain renormalization-group inspired Ansatze for the kernel are then proposed, and their numerical implications are explored in detail. One of the solutions obtained fulfills the theoretical expectations to a high degree of accuracy, yielding a gluon mass that is positive definite throughout the entire range of physical momenta, and displays in the ultraviolet the so-called “power-law” running, in agreement with standard arguments based on the operator product expansion. Some of the technical difficulties thwarting a more rigorous determination of the kernel are discussed, and possible future directions are briefly mentioned.

Abstract: In this paper we use two different but complementary approaches in order to study the ghost propagator of a pure SU(3) Yang-Mills theory quantized in the linear covariant gauges, focusing on its dependence on the gauge-fixing parameter xi in the deep infrared. In particular, we first solve the Schwinger-Dyson equation that governs the dynamics of the ghost propagator, using a set of simplifying approximations, and under the crucial assumption that the gluon propagators for xi > 0 are infrared finite, as is the case in the Landau gauge (xi = 0). Then we appeal to the Nielsen identities, and express the derivative of the ghost propagator with respect to xi in terms of certain auxiliary Green's functions, which are subsequently computed under the same assumptions as before. Within both formalisms we find that for xi > 0 the ghost dressing function approaches zero in the deep infrared, in sharp contrast to what happens in the Landau gauge, where it is known to saturate at a finite (nonvanishing) value. The Nielsen identities are then extended to the case of the gluon propagator, and the xi-dependence of the corresponding gluon masses is derived using as input the results obtained in the previous steps. The result turns out to be logarithmically divergent in the deep infrared; the compatibility of this behavior with the basic assumption of a finite gluon propagator is discussed, and a specific Ansatz is put forth, which readily reconciles both features.

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.

Abstract: In this work we explore the applicability of a special gluon mass generating mechanism in the context of the linear covariant gauges. In particular, the implementation of the Schwinger mechanism in pure Yang-Mills theories hinges crucially on the inclusion of massless bound-state excitations in the fundamental nonperturbative vertices of the theory. The dynamical formation of such excitations is controlled by a homogeneous linear Bethe-Salpeter equation, whose nontrivial solutions have been studied only in the Landau gauge. Here, the form of this integral equation is derived for general values of the gauge-fixing parameter, under a number of simplifying assumptions that reduce the degree of technical complexity. The kernel of this equation consists of fully dressed gluon propagators, for which recent lattice data are used as input, and of three-gluon vertices dressed by a single form factor, which is modeled by means of certain physically motivated Ansatze. The gauge-dependent terms contributing to this kernel impose considerable restrictions on the infrared behavior of the vertex form factor; specifically, only infrared finite Ansatze are compatible with the existence of nontrivial solutions. When such Ansatze are employed, the numerical study of the integral equation reveals a continuity in the type of solutions as one varies the gauge-fixing parameter, indicating a smooth departure from the Landau gauge. Instead, the logarithmically divergent form factor displaying the characteristic “zero crossing,” while perfectly consistent in the Landau gauge, has to undergo a dramatic qualitative transformation away from it, in order to yield acceptable solutions. The possible implications of these results are briefly discussed.