|
|
Krupczak, R., da Silva, T. N., Domingues, T. S., Luzum, M., Denicol, G. S., Gardim, F. G., et al. (2024). Causality violations in simulations of large and small heavy-ion collisions. Phys. Rev. C, 109(3), 034908–12pp.
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.
|
|
|
|
Aguilar, A. C., Ferreira, M. N., & Papavassiliou, J. (2022). Exploring smoking-gun signals of the Schwinger mechanism in QCD. Phys. Rev. D, 105(1), 014030–26pp.
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.
|
|
|
|
Ferreira, M. N., & Papavassiliou, J. (2023). Gauge Sector Dynamics in QCD. Particles, 6(1), 312–363.
Abstract: The dynamics of the QCD gauge sector give rise to non-perturbative phenomena that are crucial for the internal consistency of the theory; most notably, they account for the generation of a gluon mass through the action of the Schwinger mechanism, the taming of the Landau pole, the ensuing stabilization of the gauge coupling, and the infrared suppression of the three-gluon vertex. In the present work, we review some key advances in the ongoing investigation of this sector within the framework of the continuum Schwinger function methods, supplemented by results obtained from lattice simulations.
|
|
|
|
Aguilar, A. C., Ambrosio, C. O., De Soto, F., Ferreira, M. N., Oliveira, B. M., Papavassiliou, J., et al. (2021). Ghost dynamics in the soft gluon limit. Phys. Rev. D, 104(5), 054028–18pp.
Abstract: We present a detailed study of the dynamics associated with the ghost sector of quenched QCD in the Landau gauge, where the relevant dynamical equations are supplemented with key inputs originating from large-volume lattice simulations. In particular, we solve the coupled system of Schwinger-Dyson equations that governs the evolution of the ghost dressing function and the ghost-gluon vertex, using as input for the gluon propagator lattice data that have been cured from volume and discretization artifacts. In addition, we explore the soft gluon limit of the same system, employing recent lattice data for the three-gluon vertex that enters in one of the diagrams defining the Schwinger-Dyson equation of the ghost-gluon vertex. The results obtained from the numerical treatment of these equations are in excellent agreement with lattice data for the ghost dressing function, once the latter have undergone the appropriate scale-setting and artifact elimination refinements. Moreover, the coincidence observed between the ghost-gluon vertex in general kinematics and in the soft gluon limit reveals an outstanding consistency of physical concepts and computational schemes.
|
|
|
|
Aguilar, A. C., Ferreira, M. N., & Papavassiliou, J. (2021). Gluon dynamics from an ordinary differential equation. Eur. Phys. J. C, 81(1), 54–20pp.
Abstract: We present a novel method for computing the nonperturbative kinetic term of the gluon propagator from an ordinary differential equation, whose derivation hinges on the central hypothesis that the regular part of the three-gluon vertex and the aforementioned kinetic term are related by a partial Slavnov-Taylor identity. The main ingredients entering in the solution are projection of the three-gluon vertex and a particular derivative of the ghost-gluon kernel, whose approximate form is derived from a Schwinger-Dyson equation. Crucially, the requirement of a pole-free answer determines the initial condition, whose value is calculated from an integral containing the same ingredients as the solution itself. This feature fixes uniquely, at least in principle, the form of the kinetic term, once the ingredients have been accurately evaluated. In practice, however, due to substantial uncertainties in the computation of the necessary inputs, certain crucial components need be adjusted by hand, in order to obtain self-consistent results. Furthermore, if the gluon propagator has been independently accessed from the lattice, the solution for the kinetic term facilitates the extraction of the momentum-dependent effective gluon mass. The practical implementation of this method is carried out in detail, and the required approximations and theoretical assumptions are duly highlighted.
|
|
|
|
Aguilar, A. C., Ferreira, M. N., Figueiredo, C. T., & Papavassiliou, J. (2019). Gluon mass scale through nonlinearities and vertex interplay. Phys. Rev. D, 100(9), 094039–19pp.
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.
|
|
|
|
Aguilar, A. C., De Soto, F., Ferreira, M. N., Papavassiliou, J., Rodriguez-Quintero, J., & Zafeiropoulos, S. (2020). Gluon propagator and three-gluon vertex with dynamical quarks. Eur. Phys. J. C, 80(2), 154–17pp.
Abstract: We present a detailed analysis of the kinetic and mass terms associated with the Landau gauge gluon propagator in the presence of dynamical quarks, and a comprehensive dynamical study of certain special kinematic limits of the three-gluon vertex. Our approach capitalizes on results from recent lattice simulations with (2+1) domain wall fermions, a novel nonlinear treatment of the gluon mass equation, and the nonperturbative reconstruction of the longitudinal three-gluon vertex from its fundamental Slavnov-Taylor identities. Particular emphasis is placed on the persistence of the suppression displayed by certain combinations of the vertex form factors at intermediate and low momenta, already known from numerous pure Yang-Mills studies. One of our central findings is that the inclusion of dynamical quarks moderates the intensity of this phenomenon only mildly, leaving the asymptotic low-momentum behavior unaltered, but displaces the characteristic “zero crossing” deeper into the infrared region. In addition, the effect of the three-gluon vertex is explored at the level of the effective gauge coupling, whose size is considerably reduced with respect to its counterpart obtained from the ghost-gluon vertex. The main upshot of the above considerations is the further confirmation of the tightly interwoven dynamics between the two- and three-point sectors of QCD.
|
|
|
|
Aguilar, A. C., De Soto, F., Ferreira, M. N., Papavassiliou, J., & Rodriguez-Quintero, J. (2021). Infrared facets of the three-gluon vertex. Phys. Lett. B, 818, 136352–7pp.
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.
|
|
|
|
Pinto-Gomez, F., De Soto, F., Ferreira, M. N., Papavassiliou, J., & Rodriguez-Quintero, J. (2023). Lattice three-gluon vertex in extended kinematics: Planar degeneracy. Phys. Lett. B, 838, 137737–8pp.
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.
|
|
|
|
Aguilar, A. C., Cardona, J. C., Ferreira, M. N., & Papavassiliou, J. (2017). Non-Abelian Ball-Chiu vertex for arbitrary Euclidean momenta. Phys. Rev. D, 96(1), 014029–29pp.
Abstract: We determine the non-Abelian version of the four nontransverse form factors of the quark-gluon vertex, using exact expressions derived from the Slavnov-Taylor identity that this vertex satisfies. In addition to the quark and ghost propagators, a key ingredient of the present approach is the quark-ghost scattering kernel, which is computed within the one-loop dressed approximation. The vertex form factors obtained from this procedure are evaluated for arbitrary Euclidean momenta, and display features not captured by the well-known Ball-Chiu vertex, deduced from the Abelian (ghost-free) Ward identity. Particularly interesting in this analysis is the so-called soft-gluon limit, which, unlike other kinematic configurations considered, is especially sensitive to the approximations employed for the vertex entering in the quark-ghost scattering kernel, and may even be affected by a subtle numerical instability. As an elementary application of the results obtained, we evaluate and compare certain renormalization-point-independent combinations, which contribute to the interaction kernels appearing in the standard quark gap and Bethe-Salpeter equations. In doing so, even though all form factors of the quark-gluon vertex, and in particular the transverse ones which are unconstrained by our procedure, enter nontrivially in the aforementioned kernels, only the contribution of a single form factor, corresponding to the classical (tree-level) tensor, will be considered.
|
|
