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: In this article we use the Schwinger-Dyson equations to compute the nonperturbative modifications caused to the infrared finite gluon propagator (in the Landau gauge) by the inclusion of a small number of quark families. Our basic operating assumption is that the main bulk of the effect stems from the "one-loop dressed'' quark loop contributing to the full gluon self-energy. This quark loop is then calculated, using as basic ingredients the full quark propagator and quark-gluon vertex; for the quark propagator we use the solution obtained from the quark-gap equation, while for the vertex we employ suitable Ansatze, which guarantee the transversality of the answer. The resulting effect is included as a correction to the quenched gluon propagator, obtained in recent lattice simulations. Our main finding is that the unquenched propagator displays a considerable suppression in the intermediate momentum region, which becomes more pronounced as we increase the number of active quark families. The influence of the quarks on the saturation point of the propagator cannot be reliably computed within the present scheme; the general tendency appears to be to decrease it, suggesting a corresponding increase in the effective gluon mass. The renormalization properties of our results, and the uncertainties induced by the unspecified transverse part of the quark-gluon vertex, are discussed. Finally, the gluon propagator is compared with the available unquenched lattice data, showing rather good agreement.

Abstract: We use recent lattice data on the gluon and ghost propagators, as well as the Kugo-Ojima function, in order to extract the non-perturbative behavior of two particular definitions of the QCD effective charge, one based on the pinch technique construction, and one obtained from the standard ghost-gluon vertex. The construction relies crucially on the definition of two dimensionful quantities, which are invariant under the renormalization group, and are built out of very particular combinations of the aforementioned Green's functions. The main non-perturbative feature of both effective charges, encoded in the infrared finiteness of the gluon propagator and ghost dressing function used in their definition, is the freezing at a common finite (non-vanishing) value, in agreement with a plethora of theoretical and phenomenological expectations. We discuss the sizable discrepancy between the freezing values obtained from the present lattice analysis and the corresponding estimates derived from several phenomenological studies, and attribute its origin to the difference in the gauges employed. A particular toy calculation suggests that the modifications induced to the non-perturbative gluon propagator by the gauge choice may indeed account for the observed deviation of the freezing values.