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 novel lattice results for the form factors of the quenched three-gluon vertex of QCD, in two special kinematic configurations that depend on a single momentum scale. We consider three form factors, two associated with a classical tensor structure and one without tree-level counterpart, exhibiting markedly different infrared behaviors. Specifically, while the former display the typical suppression driven by a negative logarithmic singularity at the origin, the latter saturates at a small negative constant. These exceptional features are analyzed within the Schwinger-Dyson framework, with the aid of special relations obtained from the Slavnov-Taylor identities of the theory. The emerging picture of the underlying dynamics is thoroughly corroborated by the lattice results, both qualitatively as well as quantitatively.

Abstract: For special kinematic configurations involving a single momentum scale, certain standard relations, originating from the Slavnov-Taylor identities of the theory, may be interpreted as ordinary differential equations for the “kinetic term” of the gluon propagator. The exact solutions of these equations exhibit poles at the origin, which are incompatible with the physical answer, known to diverge only logarithmically; their elimination hinges on the validity of two integral conditions that we denominate “asymmetric” and “symmetric” sum rules, depending on the kinematics employed in their derivation. The corresponding integrands contain components of the three-gluon vertex and the ghost-gluon kernel, whose dynamics are constrained when the sum rules are imposed. For the numerical treatment we single out the asymmetric sum rule, given that its support stems predominantly from low and intermediate energy regimes of the defining integral, which are physically more interesting. Adopting a combined approach based on Schwinger-Dyson equations and lattice simulations, we demonstrate how the sum rule clearly favors the suppression of an effective form factor entering in the definition of its kernel. The results of the present work offer an additional vantage point into the rich and complex structure of the three-point sector of QCD.

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.

Abstract: We present the detailed derivation of the longitudinal part of the three-gluon vertex from the Slavnov-Taylor identities that it satisfies, by means of a nonperturbative implementation of the Ball-Chiu construction; the latter, in its original form, involves the inverse gluon propagator, the ghost dressing function, and certain form factors of the ghost-gluon kernel. The main conceptual subtlety that renders this endeavor nontrivial is the infrared finiteness of the gluon propagator, and the resulting need to separate the vertex into two pieces, one that is intimately connected with the emergence of a gluonic mass scale, and one that satisfies the original set of Slavnov-Taylor identities, but with the inverse gluon propagator replaced by its “kinetic” term. The longitudinal form factors obtained by this construction are presented for arbitrary Euclidean momenta, as well as special kinematic configurations, parametrized by a single momentum. A particularly preeminent feature of the components comprising the tree-level vertex is their considerable suppression for momenta below 1 GeV, and the appearance of the characteristic “zero-crossing” in the vicinity of 100-200 MeV. Special combinations of the form factors derived with this method are compared with the results of recent large-volume lattice simulations, and are found to capture faithfully the rather complicated curves formed by the data. A similar comparison with results obtained from Schwinger-Dyson equations reveals a fair overall agreement, but with appreciable differences at intermediate energies. A variety of issues related to the distribution of the pole terms responsible for the gluon mass generation are discussed in detail, and their impact on the structure of the transverse parts is elucidated. In addition, a brief account of several theoretical and phenomenological possibilities involving these newly acquired results is presented.