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.

Abstract: Dynamical gluon mass generation has been traditionally plagued with seagull divergences, and all regularization procedures proposed over the years yield finite but scheme-dependent gluon masses. In this work we show how such divergences can be eliminated completely by virtue of a characteristic identity, valid in dimensional regularization. The ability to trigger the aforementioned identity hinges crucially on the particular Ansatz employed for the three-gluon vertex entering into the Schwinger-Dyson equation governing the gluon propagator. The use of the appropriate three-gluon vertex brings about an additional advantage: one obtains two separate (but coupled) integral equations, one for the effective charge and one for the gluon mass. This system of integral equations has a unique solution, which unambiguously determines these two quantities. Most notably, the effective charge freezes in the infrared, and the gluon mass displays power-law running in the ultraviolet, in agreement with earlier considerations.

Abstract: We study chiral symmetry breaking using the standard gap equation, supplemented with the infrared-finite gluon propagator and ghost dressing function obtained from large-volume lattice simulations. One of the most important ingredients of this analysis is the non-Abelian quark-gluon vertex, which controls the way the ghost sector enters into the gap equation. Specifically, this vertex introduces a numerically crucial dependence on the ghost dressing function and the quark-ghost scattering amplitude. This latter quantity satisfies its own, previously unexplored, dynamical equation, which may be decomposed into individual integral equations for its various form factors. In particular, the scalar form factor is obtained from an approximate version of the “one-loop dressed” integral equation, and its numerical impact turns out to be rather considerable. The detailed numerical analysis of the resulting gap equation reveals that the constituent quark mass obtained is about 300 MeV, while fermions in the adjoint representation acquire a mass in the range of (750-962) MeV.

Abstract: The gauge-invariant generation of an effective gluon mass proceeds through the well-known Schwinger mechanism, whose key dynamical ingredient is the nonperturbative formation of longitudinally coupled massless bound-state excitations. These excitations introduce poles in the vertices of the theory, in such a way as to maintain the Slavnov-Taylor identities intact in the presence of massive gluon propagators. In the present work we first focus on the modifications induced to the nonperturbative three-gluon vertex by the inclusion of massless two-gluon bound states into the kernels appearing in its skeleton expansion. Certain general relations between the basic building blocks of these bound states and the gluon mass are then obtained from the Slavnov-Taylor identities and the Schwinger-Dyson equation governing the gluon propagator. The homogeneous Bethe-Salpeter equation determining the wave function of the aforementioned bound state is then derived, under certain simplifying assumptions. It is then shown, through a detailed analytical and numerical study, that this equation admits nontrivial solutions, indicating that the QCD dynamics support indeed the formation of such massless bound states. These solutions are subsequently used, in conjunction with the aforementioned relations, to determine the momentumdependence of the dynamical gluon mass. Finally, further possibilities and open questions are briefly discussed.

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.