Abstract: Continuum and lattice analyses have revealed the existence of a mass-scale in the gluon two-point Schwinger function. It has long been conjectured that this expresses the action of a Schwinger mechanism for gauge boson mass generation in quantum chromodynamics (QCD). For such to be true, it is necessary and sufficient that a dynamically-generated, massless, colour-carrying, scalar gluon+gluon correlation emerges as a feature of the dressed three-gluon vertex. Working with results on elementary Schwinger functions obtained via the numerical simulation of lattice-regularised QCD, we establish with an extremely high level of confidence that just such a feature appears; hence, confirm the conjectured origin of the gluon mass scale.

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 present a novel analysis of the gluon gap equation, where its full nonlinear structure is duly taken into account. In particular, while in previous treatments the linearization of this homogeneous integral equation introduced an indeterminacy in the scale of the corresponding mass, the current approach determines it uniquely, once the value of the gauge coupling at a given renormalization point is used as input. A crucial ingredient for this construction is the “kinetic term” of the gluon propagator, whose form is not obtained from the complicated equation governing its evolution, but is rather approximated by suitable initial Ansatze, which are subsequently improved by means of a systematic iterative procedure. The multiplicative renormalization of the central equation is carried out following an approximate method, which is extensively employed in the studies of the standard quark gap equation. This approach amounts to the effective substitution of the vertex renormalization constants by kinematically simplified form factors of the three- and four-gluon vertices. The resulting numerical interplay, exemplified by the infrared suppression of the three-gluon vertex and the mild enhancement of the four-gluon vertex, is instrumental for obtaining positive-definite and monotonically decreasing running gluon masses. The resulting gluon propagators, put together from the gluon masses and kinetic terms obtained with this method, match rather accurately the data obtained from large-volume lattice simulations.