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: In the present work, we discuss certain characteristic features encoded in some of the fundamental QCD Green's functions, for which the origin can be traced back to the nonperturbative masslessness of the ghost field, in the Landau gauge. Specifically, the ghost loops that contribute to these Green's functions display infrared divergences, akin to those encountered in the perturbative treatment, in contradistinction to the gluonic loops, for which perturbative divergences are tamed by the dynamical generation of an effective gluon mass. In d = 4, the aforementioned divergences are logarithmic, thus causing a relatively mild impact, whereas in d = 3 they are linear, giving rise to enhanced effects. In the case of the gluon propagator, these effects do not interfere with its finiteness, but make its first derivative diverge at the origin, and introduce a maximum in the region of infrared momenta. The three-gluon vertex is also affected, and the induced divergent behavior is clearly exposed in certain special kinematic configurations, usually considered in lattice simulations; the sign of the corresponding divergence is unambiguously determined. The main underlying concepts are developed in the context of a simple toy model, which demonstrates clearly the interconnected nature of the various effects. The picture that emerges is subsequently corroborated by a detailed nonperturbative analysis, combining lattice results with the dynamical integral equations governing the relevant ingredients, such as the nonperturbative ghost loop and the momentumdependent gluon mass.

Abstract: In Quantum Chromodynamics, the Schwinger mechanism endows the gluons with an effective mass through the dynamical formation of massless bound-state poles that are longitudinally coupled. The presence of these poles affects profoundly the infrared properties of the interaction vertices, inducing crucial modifications to their fundamental Ward identities. Within this general framework, we present a detailed derivation of the non-Abelian Ward identity obeyed by the pole-free part of the three-gluon vertex in the softgluon limit, and determine the smoking-gun displacement that the onset of the Schwinger mechanism produces to the standard result. Quite importantly, the quantity that describes this distinctive feature coincides formally with the bound-state wave function that controls the massless pole formation. Consequently, this signal may be computed in two independent ways: by solving an approximate version of the pertinent BetheSalpeter integral equation, or by appropriately combining the elements that enter in the aforementioned Ward identity. For the implementation of both methods we employ two- and three-point correlation functions obtained from recent lattice simulations, and a partial derivative of the ghost-gluon kernel, which is computed from the corresponding Schwinger-Dyson equation. Our analysis reveals an excellent coincidence between the results obtained through either method, providing a highly nontrivial self-consistency check for the entire approach. When compared to the null hypothesis, where the Schwinger mechanism is assumed to be inactive, the statistical significance of the resulting signal is estimated to be 3 standard deviations.

Abstract: We study an approximate version of the Schwinger-Dyson equation that controls the nonperturbative behavior of the ghost-gluon vertex in the Landau gauge. In particular, we focus on the form factor that enters in the dynamical equation for the ghost dressing function, in the same gauge, and derive its integral equation, in the “one-loop dressed” approximation. We consider two special kinematic configurations, which simplify the momentum dependence of the unknown quantity; in particular, we study the soft gluon case and the well-known Taylor limit. When coupled with the Schwinger-Dyson equation of the ghost dressing function, the contribution of this form factor provides considerable support to the relevant integral kernel. As a consequence, the solution of this coupled system of integral equations furnishes a ghost dressing function that reproduces the standard lattice results rather accurately, without the need to artificially increase the value of the gauge coupling.

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.