Abstract: We present novel results for the three-gluon vertex, obtained from an extensive quenched lattice simulation in the Landau gauge. The simulation evaluates the transversely projected vertex, spanned on a special tensorial basis, whose form factors are naturally parametrized in terms of individually Bosesymmetric variables. Quite interestingly, when evaluated in these kinematics, the corresponding form factors depend almost exclusively on a single kinematic variable, formed by the sum of the squares of the three incoming four-momenta, q, r, and p. Thus, all configurations lying on a given plane in the coordinate system (q2, r2, p2) share, to a high degree of accuracy, the same form factors, a property that we denominate planar degeneracy. We have confirmed the validity of this property through an exhaustive study of the set of configurations satisfying the condition q2 = r2, within the range [0, 5 GeV]. This drastic simplification allows for a remarkably compact description of the main bulk of the data, which is particularly suitable for future numerical applications. A semi-perturbative analysis reproduces the lattice findings rather accurately, once the inclusion of a gluon mass has cured all spurious divergences.

Abstract: Heavy-ion collisions, such as Pb-Pb or p-Pb, produce extreme conditions in temperature and density that make the hadronic matter transition to a new state, called quark-gluon plasma (QGP). Simulations of heavy-ion collisions provide a way to improve our understanding of the QGP's properties. These simulations are composed of a hybrid description that results in final observables in agreement with accelerators like LHC and RHIC. However, recent works pointed out that these hydrodynamic simulations can display acausal behavior during the evolution in certain regions, indicating a deviation from a faithful representation of the underlying QCD dynamics. To pursue a better understanding of this problem and its consequences, this work simulated two different collision systems, Pb-Pb and p-Pb at root sNN = 5.02 TeV. In this context, our results show that causality violation, even though always present, typically occurs on a small part of the system, quantified by the total energy fraction residing in the acausal region. In addition, the acausal behavior can be reduced with changes in the prehydrodynamic factors and the definition of the bulk-viscous relaxation time. Since these aspects are fairly arbitrary in current simulation models, without solid guidance from the underlying theory, it is reasonable to use the disturbing presence of acausal behavior in current simulations to guide improvements towards more realistic modeling. While this work does not solve the acausality problem, it sheds more light on this issue and also proposes a way to solve this problem in simulations of heavy-ion collisions.

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.

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.