Next Article in Journal
Feasibility Studies for the Measurement of Open-Charm Mesons at ALICE-3 Using Decay Channels with Neutral Mesons and Photons in the Final State
Next Article in Special Issue
Precision Storage Rings for Electric Dipole Moment Searches: A Tool En Route to Physics Beyond-the-Standard-Model
Previous Article in Journal
Long-Lived Particles Anomaly Detection with Parametrized Quantum Circuits
Previous Article in Special Issue
Generalised Parton Distributions in Continuum Schwinger Methods: Progresses, Opportunities and Challenges
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Gauge Sector Dynamics in QCD

by
Mauricio Narciso Ferreira
*,† and
Joannis Papavassiliou
*,†
Department of Theoretical Physics and IFIC, University of Valencia and CSIC, E-46100 Valencia, Spain
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Particles 2023, 6(1), 312-363; https://doi.org/10.3390/particles6010017
Submission received: 7 January 2023 / Revised: 3 February 2023 / Accepted: 9 February 2023 / Published: 15 February 2023

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.

Contents
 
1      Introduction2
 
2      Basic Concepts and General Theoretical Framework5
 
3      Schwinger Mechanism in Yang–Mills Theories10
 
4      Dynamical formation of Massless Poles13
 
5      Generation of the Gluon Mass16
5.1      Gluon Mass from the qμqν Component17
5.2      Gluon Mass from the gμν Component: Seagull Identity and Ward Identity Displacement18
 
6      Renormalization Group Invariant Interaction Strength20
 
7      Three-Gluon Vertex and Its Planar Degeneracy23
 
8      Ghost Dynamics from Schwinger–Dyson Equations25
 
9      Divergent Ghost Loops and Their Impact on the QCD Green’s Functions30
 
10      Ward Identity Displacement of the Three-Gluon Vertex34
 
11      The Ghost-Gluon Kernel Contribution to the Ward Identity35
 
12      Displacement Function from Lattice Inputs39
 
13      Conclusions40
A.      Appendix A41
B.      Appendix B43
 
References43

1. Introduction

The systematic exploration of Green’s functions (n-point correlation functions) of quantum chromodynamics (QCD) [1] by means of continuous Schwinger function methods [2,3,4,5,6,7,8,9], such as Schwinger–Dyson equations (SDEs) [10,11,12,13,14,15,16,17,18,19,20,21] and the functional renormalization group [22,23,24,25,26,27,28,29,30,31], together with a plethora of gauge-fixed lattice simulations [32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79,80,81,82,83,84,85,86,87,88], has afforded ample access to the dynamical mechanisms responsible for the non-perturbative properties of this remarkable theory. Particularly prominent in this quest is the notion of the emergent hadron mass (EHM) [3,8,9,89,90,91,92,93], together with its three supporting pillars: first, the generation of a gluon mass [18,32,93,94,95,96,97,98,99,100,101,102,103,104,105,106,107,108,109,110,111,112,113,114,115,116,117,118,119,120,121,122,123,124,125,126] through the action of the Schwinger mechanism [127,128]; second, the construction of the process-independent effective charge [3,16,20,79,96,129,130,131], which arises as the QCD analog of the Gell-Mann–Low charge is known from quantum electrodynamics (QED) [132,133], and is associated with a renormalization-group invariant (RGI) scale of about half of the proton mass [20,79]; and third, the dynamical breaking of chiral symmetry and the generation of constituent quark masses [10,17,134,135,136,137,138,139,140,141,142,143,144,145,146,147,148,149,150,151,152,153,154,155,156,157,158].
The dynamics of the gauge sector of QCD, which encompasses both gluonic and ghost interactions, is instrumental in the physical picture of the EHM outlined above. In fact, the basic concepts and pivotal mechanisms sustaining the first two pillars of the EHM have their original inception and most genuine realization in the realm of pure Yang–Mills theories [18,93,94,96,109,112,117,159,160,161]. Therefore, in the present review, we focus precisely on the rich dynamical content of the gauge sector, especially in relation to the generation of a gluon mass scale out of the intricate gluon self-interactions.
The formulation of the non-perturbative QCD physics in terms of Green’s functions of the fundamental degrees of freedom, such as gluon and ghost propagators and vertices, provides an intuitive framework for unraveling a wide array of subtle mechanisms; in fact, certain distinctive features of these functions have been inextricably connected with key phenomena such as gluon mass generation, violation of reflection positivity, and confinement, to name a few. Thus, the saturation of the gluon propagator in the deep infrared [37,45,46,47,48,49,52,55,56,57,58,59,63,65,66,67,77,81] has been interpreted as an unequivocal signal of a gluon mass [32,96,97,98,99,100,103,105,107,108,109,112,160,161,162,163,164,165,166]; and the existence of an inflection point in the same function has been argued to lead to a non-positive gluon spectral density [8], and the ensuing loss of reflection positivity [8,11,13,16,167,168,169,170,171] for the dressed gluons. Similarly, the masslessness of the ghost induces [172] a maximum in the gluon propagator, and a zero crossing in the form factors of the three-gluon vertex [28,50,68,69,71,72,81,84,172,173,174,175,176,177,178,179,180], followed by an infrared divergence for vanishing momenta. The dynamic origin of these special traits will be the focal point of the analysis presented in the main body of this article.
The integral equations that govern the full momentum evolution of Green’s functions, known as SDEs, constitute the indispensable formal and practical instrument for unraveling the special characteristics mentioned above. In their primordial form, the SDEs are rigorously derived from the generating functional of the theory [133,181], and encode all dynamical information on the correlation functions, within the entire range of physical momenta. In practice, due to the enormous complexity of these equations, truncation approximations need to be implemented; but, unlike perturbation theory, no expansion parameter is available in the strongly coupled regime of the theory for carrying out such a task. Despite this intrinsic shortcoming, in recent years, the SDE predictions have become particularly robust, in part due to various theoretical advances, and in part thanks to the intense synergy with gauge-fixed lattice simulations, as will be evidenced in subsequent sections.
Typically, Green’s QCD function is defined within the quantization scheme obtained by implementing the linear covariant ( R ξ ) gauges [182]. The corresponding SDEs are derived and solved within this same quantization scheme, particularly in the Landau gauge ( ξ = 0 ), where lattice simulations are almost exclusively performed; for studies away from the Landau gauge, see e.g., [55,58,66,74,75,110,114,120,183,184,185,186,187,188,189,190,191]. A great deal may be learned, however, by considering Green’s functions and corresponding SDEs formulated within the “PT-BFM” scheme [109,192], namely the framework that arises from the fusion of the pinch technique (PT) [14,96,100,193,194,195] with the background field method (BFM) [196,197,198,199,200,201,202,203,204,205,206]. The main advantage of the PT-BFM originates from the fact that certain appropriately chosen Green’s functions satisfy Abelian Slavnov–Taylor identities (STIs), whose tree-level form does not get modified by quantum corrections. This situation is to be contrasted to the standard STIs [207,208] obtained in the conventional framework of the linear covariant gauges, which are deformed by non-trivial contributions stemming from the gauge sector of the theory. In the present work, we will carry out computations and develop arguments within both frameworks ( R ξ and PT-BFM), and will elaborate on their connection by means of the so-called background-quantum identities (BQIs) [14,209,210,211].
The article is organized as follows:
  • In Section 2, we introduce some basic notations and review certain prominent features of Green’s functions within the linear gauges and the PT-BFM formalism [109,192]. We stress, in particular, the properties of the auxiliary function G ( q ) [16,131,212,213], which relates the gluon propagators with quantum and background gluons, and is intimately connected with the definition of the process-independent and RGI interaction strength [16], to be discussed in detail in Section 6. In addition, we elucidate (with a concrete example) the important property of “block-wise” transversality, displayed by the background gluon self-energy [18,109,112].
  • In Section 3, we review the general principles associated with the Schwinger mechanism [127,128] that endows gauge bosons with an effective mass, focusing on the details associated with its realization in the context of Yang–Mills theories. We place particular emphasis on the pivotal requirement that must be satisfied by the fundamental vertices of the theory, namely the appearance of massless poles in their form factors [18,93,109,111,112,113,117,159,214].
  • In Section 4, we examine the dynamical formation of colored composite excitations (bound states) of vanishing masses, which provide the required structures in the vertices in order for the Schwinger mechanism to be activated [18,117,159,214]. The formation of these states out of a pair of gluons or a ghost–anti-ghost pair is controlled by a set of coupled Bethe–Salpeter equations (BSEs) [18,117,124,214,215], which are found to have nontrivial solutions for the corresponding Bethe–Salpeter (BS) amplitudes, to be denoted by C ( r ) and C ( r ) , respectively.
  • In Section 5, we explain in detail how the presence of the massless poles in the dressed vertices that enter the SDE of the gluon propagator give rise to a gluon mass. The demonstration is carried out separately for the g μ ν and q μ q ν / q 2 components of the gluon self-energy. The former case requires the evasion of the so-called “seagull identity” [113,166]; this becomes possible by virtue of the crucial Ward identity (WI) displacement, to be further considered in Section 10.
  • In Section 6, we go over the basic notions underpinning the PT [14,96,100,193,194], and show how their application leads naturally to the definition of a dimensionful process-independent RGI interaction strength [3,16,20,79,96,129,130,131], denoted by d ^ ( q ) . The genuine process independence of this quantity is concretely exemplified by demonstrating its appearance in two processes involving fundamentally different external fields. Next, d ^ ( q ) is computed by combining lattice data for the gluon propagator and SDE results for the function G ( q ) . Finally, the dimensionless quantity is derived that constitutes the physical definition of the one-gluon exchange interaction appearing in standard bound-state computations [15,16,17,216,217,218,219,220,221,222].
  • In Section 7, we focus on the structure of the “transversely projected” three-gluon vertex [126,174,175,223], and discuss briefly the property of planar degeneracy [86], satisfied, at a high level of accuracy [86,87,88,174,175,223], by the vertex form factors. This special property induces a striking simplification to the structure of this vertex, captured by a particularly compact expression [86], which will be extensively used in some of the following sections.
  • In Section 8, we take a close look at the ghost sector of the theory, and solve the coupled system of SDEs governing the ghost propagator and ghost–gluon vertex [85,224,225,226,227,228]; as is well-known, the ghost remains massless, but its dressing function saturates at the origin [21,42,47,49,51,56,62,63,73,79,85,112,178,225,227,228,229,230,231,232,233], because the infrared-finite gluon propagator used in the ghost SDE provides an effective infrared cutoff. In the SDE of the ghost–gluon vertex, we employ as central input the compact expression for the three-gluon vertex presented in the previous section. The results are in excellent agreement with the available lattice data for the ghost dressing function [73,85] and the form factor of the ghost–gluon vertex evaluated in the soft-gluon limit [42,43].
  • In Section 9, we discuss two important consequences of the masslessness of the ghost propagator, which manifest themselves at the level of both the gluon propagator and the three-gluon vertex. Specifically, the diagrams comprised by a ghost loop induce “unprotected” logarithms, i.e., of the type ln q 2 ; instead, gluonic loops give rise to “protected” logarithms, of the type ln ( q 2 + m 2 ) , where m is the effective gluon mass [172,234]. As q 2 0 , the unprotected contributions diverge, driving the appearance of a maximum in the gluon propagator and a divergence in its first derivative, as well as a zero-crossing and a corresponding divergence in the form factors of the three-gluon vertex. As we comment in this section, of particular phenomenological importance [234,235,236,237,238,239,240] is the relative suppression that the above features induce to the dominant vertex form factors in the intermediate range of momenta.
  • In Section 10, we discuss an outstanding feature of the WI satisfied by the pole-free part of the three-gluon vertex, namely the displacement induced by the presence of the aforementioned massless poles [93,124]. In this context, we introduce the key quantity denominated “displacement function”, whose appearance serves as a smoking gun signal of the action of the Schwinger mechanism in QCD; quite interestingly, it coincides [93,124] with the BS amplitude C ( r ) for the formation of a massless scalar out of a pair of gluons, introduced in Section 4. In addition, we derive a crucial relation, which ultimately permits the indirect determination of C ( r ) from lattice QCD [93,124,126]; an important ingredient in this relation is a partial derivative [124,241], denoted by W ( r ) , of the ghost–gluon kernel [228], to be determined in the next section.
  • In Section 11, we set up and solve the SDE that governs the evolution of W ( r ) [124,126,241,242]; the main component of this SDE is a special projection of the three-gluon vertex, which is computed by appealing to formulas established in Section 7, and allows for the accurate determination of W ( r ) in the entire range of relevant momenta [126].
  • In Section 12, we substitute into the central relation derived in Section 10 the solution for W ( r ) found in the previous section, together with the lattice data [84,85] for the gluon propagator, the ghost dressing function, and the form factor of the three-gluon vertex associated with the soft-gluon limit, in order to obtain the form of the displacement function C ( r ) [124,126]. As we discuss, the results exclude—with near-absolute certainty—the null hypothesis (absence of Schwinger mechanism, C ( r ) = 0 ), and corroborate the action of the Schwinger mechanism in QCD [126]. In addition, we show that the form of C ( r ) found is statistically completely compatible with that obtained from the BSE-based analysis presented in Section 4.
  • In Section 13, we present our conclusions.
  • Finally, in Appendix A, we derive the BQIs related to the displacement functions of the conventional and background vertices, while in Appendix B, we provide details about the renormalization scheme employed in our computations.

2. Basic Concepts and General Theoretical Framework

We start by considering the Lagrangian density of an SU(N) Yang–Mills theory, comprised of the classical part, L cl , the contribution from the ghosts, L gh , and the covariant gauge-fixing term, L gf , namely
L YM = L cl + L gh + L gf ,
where
L cl = 1 4 F μ ν a F a μ ν , L gh = c ¯ a μ D μ a b c b , L gf = 1 2 ξ ( μ A μ a ) 2 .
In the above formula, A μ a ( x ) denotes the gauge field, while c a ( x ) and c ¯ a ( x ) represent the ghost and anti-ghost fields, respectively, with a = 1 , , N 2 1 .
In addition,
F μ ν a = μ A ν a ν A μ a + g f a b c A μ b A ν c ,
is the antisymmetric field tensor, where f a b c stands for the totally antisymmetric structure constants of the SU(N) gauge group, and g is the gauge coupling, while
D μ a b = μ δ a c + g f a m b A μ m ,
denotes the covariant derivative in the adjoint representation. Finally, ξ represents the gauge-fixing parameter; ξ = 0 corresponds to the Landau gauge, while ξ = 1 specifies the Feynman–´t Hooft gauge.
The transition from the pure Yang–Mills theory of Equation (1) to QCD is implemented by supplementing the corresponding kinetic and interaction terms for the quark fields. However, since throughout this work we do not consider effects due to dynamical quarks, the aforementioned terms will be omitted entirely.
The most fundamental correlation function is the gluon propagator, whose non-perturbative features are inextricably connected with key dynamical properties of the theory. In the Landau gauge that we will employ throughout, the gluon propagator, Δ μ ν a b ( q ) = i δ a b Δ μ ν ( q ) , is completely transverse, i.e.,
Δ μ ν ( q ) = Δ ( q ) P μ ν ( q ) , P μ ν ( q ) : = g μ ν q μ q ν / q 2 .
In the continuum, the dynamical properties of the gluon propagator are encoded in the corresponding SDE, given by
Δ 1 ( q ) P μ ν ( q ) = q 2 P μ ν ( q ) + i Π μ ν ( q ) ,
where Π μ ν ( q ) is the gluon self-energy, shown diagrammatically in the first row of Figure 1. The fully-dressed vertices entering the diagrams are determined by their own SDEs, obtaining finally a tower of coupled integral equations, which, for practical purposes, must be truncated or treated approximately.
Given that, by virtue of the fundamental STI satisfied by the two-point function, the self-energy Π μ ν ( q ) is transverse,
q μ Π μ ν ( q ) = 0 ,
we have that
Π μ ν ( q ) = Π ( q ) P μ ν ( q ) ,
and from Equation (6) follows that
Δ 1 ( q ) = q 2 + i Π ( q ) .
Of particular importance is the exact way that Equation (7) is enforced at the level of the SDE given in Figure 1 which governs the gluon evolution. In particular, if we were to contract the corresponding diagrams by q μ , the entire set of diagrams must be considered in order for Equation (7) to emerge from the SDE. This pattern manifests itself already at the one-loop level, where it is known that the ghost loop must be included in order to guarantee the transversality of the self-energy. The main practical drawback stemming from this observation is that truncations, in the form of the omission of certain subsets of graphs, are likely to distort this fundamental property.
Quite interestingly, within the PT-BFM framework the transversality property of Equation (7) is enforced in a very special way, which permits physically meaningful truncations. In what follows we will predominantly employ the language of the BFM; for the basic principles of the PT and its connection with the BFM, the reader is referred to the extended literature on the subject [14,96,100,193,194,211,243], as well as to Section 6 of the present work.
The BFM is a powerful quantization procedure, where the gauge-fixing is implemented without compromising explicit gauge invariance. Within this framework, gauge field A appearing in the classical action is decomposed as A = B + Q , where B and Q are the background and quantum (fluctuating) fields, respectively. Note that the variable of integration in the generating functional Z ( J ) is the quantum field, which couples to the external sources, as J · Q . The background field does not appear in loops. Instead, it couples externally to the Feynman diagrams, connecting them with the asymptotic states to form elements of the S-matrix. Then, if the gauge-fixing term
L ^ gf = 1 2 ξ Q ( D ^ μ a b Q b μ ) 2 , D ^ μ a b = μ δ a b + g f a m b B μ m ,
is used, the resulting gauge-fixed action retains its invariance under gauge transformations of the background field. As a result of this invariance, when Green’s functions are contracted by the momentum carried by a background gluon, they satisfy Abelian (ghost-free) STIs, akin to the Takahashi identities known from QED. In particular, the STIs of the BFM retain their tree-level forms in all orders, in contradistinction to the STIs of the R ξ gauges, whose forms are modified by contributions stemming from the ghost sector.
Within the BFM, one may consider three kinds of propagators, by choosing the types of incoming and outgoing gluons [244]. In particular, we have:
(i)
The propagator 0 | T [ Q μ a ( q ) Q ν b ( q ) ] | 0 that connects two quantum gluons. Notice that this propagator coincides with the conventional gluon propagator of the covariant gauges, defined in Equation (5), under the assumption that the corresponding gauge-fixing parameters, ξ and ξ Q , are identified, i.e., ξ = ξ Q .
(ii)
The propagator 0 | T [ Q μ a ( q ) B ν b ( q ) ] | 0 that connects a Q μ a ( q ) with a B ν b ( q ) , to be denoted by Δ ˜ μ ν a b ( q ) = i δ a b Δ ˜ μ ν ( q ) .
(iii)
The propagator 0 | T [ B μ a ( q ) B ν b ( q ) ] | 0 that connects a B μ a ( q ) with a B ν b ( q ) , to be denoted by Δ ^ μ ν a b ( q ) = i δ a b Δ ^ μ ν ( q ) . Note that its full definition requires an additional gauge-fixing term, with the associated “classical” gauge-fixing parameter, ξ C [14,202,206].
Given that the relations captured by Equations (5) and (6) apply also in the cases of Δ ˜ μ ν ( q ) and Δ ^ μ ν ( q ) , one may define the corresponding self-energies Π ˜ μ ν ( q ) and Π ^ μ ν ( q ) , as well as the functions Δ ˜ ( q ) and Δ ^ ( q ) .
Quite interestingly, the three propagators defined in (i)-(iii) are related by a set of exact identities, known as BQIs [14,209,210,211]. In particular, we have that (see also Table 1)
Δ ( q ) = [ 1 + G ( q ) ] Δ ˜ ( q ) = [ 1 + G ( q ) ] 2 Δ ^ ( q ) ,
where the function G ( q ) is the g μ ν component of a particular two-point ghost function, Λ μ ν ( q ) , given by [209,211,213,245]
Λ μ ν ( q ) : = i g 2 C A k Δ μ ρ ( k ) D ( k + q ) H ν ρ ( q , k + q , k ) = g μ ν G ( q ) + q μ q ν q 2 L ( q ) ,
where C A is the Casimir eigenvalue of the adjoint representation [N for SU ( N ) ], D a b ( q ) = i δ a b D ( q ) is the ghost propagator, and H ν μ ( r , p , q ) denotes the ghost–gluon kernel defined in Figure 2.
In the Landau gauge, a special identity relates the form factors of Λ μ ν ( q ) to the ghost dressing function, F ( q ) , defined as F ( q ) = q 2 D ( q ) , namely [16,131,213]
F 1 ( q ) = 1 + G ( q ) + L ( q ) ,
which is valid before renormalization. In fact, in this particular gauge, G ( q ) coincides with the so-called Kugo–Ojima function [212,245,246,247].
To determine the renormalized form of Equation (13), we introduce the renormalization constants of the conventional Green’s functions
Δ R ( q ) = Z A 1 Δ ( q ) , F R ( q ) = Z c 1 F ( q ) , I Γ μ R ( r , p , q ) = Z 1 I Γ μ ( r , p , q ) , I Γ α μ ν R ( q , r , p ) = Z 3 I Γ α μ ν ( q , r , p ) , g R = Z g 1 g , g μ ν + Λ μ ν R ( q ) = Z Λ g μ ν + Λ μ ν ( q ) , Z g 1 = Z 1 1 Z A 1 / 2 Z c = Z 3 1 Z A 3 / 2 ,
where we denote by I Γ μ a b c ( r , p , q ) = g f a b c I Γ μ ( r , p , q ) and I Γ α μ ν a b c ( q , r , p ) = g f a b c I Γ α μ ν ( q , r , p ) the conventional ghost–gluon [ Q μ a ( q ) c c ( p ) c ¯ b ( r ) ] and three-gluon [ Q α a ( q ) Q μ b ( r ) Q ν c ( p ) ] vertices, respectively. Note that, by virtue of Taylor’s theorem [207], Z 1 is finite in the Landau gauge; its precise value depends on the renormalization scheme adopted, see Section 8. Moreover, denoting by Z ^ A the (wave-function) renormalization constant of Δ ^ ( q ) , the Abelian STIs of the BFM impose the validity of the pivotal relation [14,202,206]
Z g = Z ^ A 1 / 2 ,
which is the non-Abelian analog of the textbook relation Z e = Z A 1 / 2 [133], relating the renormalization constants of the electric charge and the photon propagator in QED.
Then, since the BQIs of Equation (11) are direct consequences of the Becchi–Rouet–Stora–Tyutin (BRST) symmetry [248,249,250] of the theory [209,211,213,245], the form is preserved by renormalization. Hence, by combining Equations (11), (14) and (15), we obtain
Z Λ = Z 1 1 Z c ,
which yields (note that in the original and widely used [3,8,16,20,79,131] version of Equation (17) the renormalization is performed in the so-called Taylor scheme, where Z 1 = 1 .)
Z 1 1 F 1 ( q ) = 1 + G ( q ) + L ( q ) .
As has been shown in [131], the dynamical equation governing L ( q ) yields L ( 0 ) = 0 , provided that the gluon propagator entering it is finite at the origin. Thus, one obtains from Equation (17) the useful identity [212]
Z 1 1 F 1 ( 0 ) = 1 + G ( 0 ) .
According to numerous lattice simulations and studies in the continuum (see e.g., [21,42,47,49,51,56,62,63,73,79,85,112,178,225,227,228,229,230,231,232,233]), the ghost dressing function reaches a finite (nonvanishing) value at the origin, which, due to Equation (18), furnishes also the value of G ( 0 ) .
The final upshot of the above considerations is that one may use the BQIs in Equation (11) to express the SDE given in Equation (6) in terms of the Π ˜ μ ν ( q ) or Π ^ μ ν ( q ) , at the modest cost of introducing in the dynamics the quantities 1 + G ( q ) or [ 1 + G ( q ) ] 2 . Focusing on the former possibility, Equation (11) becomes
Δ 1 ( q ) P μ ν ( q ) = q 2 P μ ν ( q ) + i Π ˜ μ ν ( q ) 1 + G ( q ) ,
where the diagrammatic representation of the self-energy Π ˜ μ ν ( q ) is shown in the lower panel of Figure 1.
The principal advantage of this formulation is that the self-energy Π ˜ μ ν ( q ) contains fully-dressed vertices with a background gluon of momentum q exiting from them, which satisfy Abelian STIs. In particular, denoting by I Γ ˜ μ α β ( q , r , p ) , I Γ ˜ μ ( r , p , q ) , and I Γ ˜ μ α β γ m n r s ( q , r , p , t ) the BQQ, Bcc, and BQQQ vertices, respectively, we have that [14,100,109]
q μ I Γ ˜ μ α β ( q , r , p ) = Δ α β 1 ( r ) Δ α β 1 ( p ) ,
q μ I Γ ˜ μ ( r , p , q ) = D 1 ( p ) D 1 ( r ) ,
q μ I Γ ˜ μ α β γ m n r s ( q , r , p , t ) = f m s e f e r n I Γ α β γ ( r , p , q + t ) + f m n e f e s r I Γ β γ α ( p , t , q + r ) + f m r e f e n s I Γ γ α β ( t , r , q + p ) .
In contrast, the conventional three-gluon and ghost–gluon vertices, I Γ α μ ν ( q , r , p ) and I Γ α ( r , p , q ) , respectively, satisfy the STIs [1,251,252,253,254,255]
q α I Γ α μ ν ( q , r , p ) = F ( q ) Δ 1 ( p ) P ν σ ( p ) H σ μ ( p , q , r ) Δ 1 ( r ) P μ σ ( r ) H σ ν ( r , q , p ) ,
q μ F 1 ( q ) I Γ μ ( r , p , q ) + p μ F 1 ( p ) I Γ μ ( r , q , p ) = r 2 F 1 ( r ) U ( r , q , p ) ,
where U ( r , q , p ) is an interaction kernel containing only ghost fields; its tree-level value is U 0 ( r , q , p ) = 1 . The STI for the conventional four-gluon vertex is given in Equation (C.24) of [14].
The special STIs listed in Equations (20)–(22) are responsible for the remarkable property of “block-wise” transversality [109,192,244], displayed by Π ˜ μ ν ( q ) . To appreciate this point, notice that the diagrams comprising Π ˜ μ ν ( q ) in Figure 1 were separated into three different subsets (blocks), consisting of (i) one-loop dressed diagrams containing only gluons, (ii) one-loop dressed diagrams containing a ghost loop, and (iii) two-loop dressed diagrams containing only gluons. The corresponding contributions of each block to Π ˜ μ ν ( q ) are denoted by Π ˜ μ ν ( i ) ( q ) , with i = 1 , 2 , 3 .
The block-wise transversality is a stronger version of the standard transversality relation q μ Π ˜ μ ν ( q ) = 0 ; it states that each block of diagrams mentioned above is individually transverse, namely
q μ Π ˜ μ ν ( i ) ( q ) = 0 , i = 1 , 2 , 3 .
In order to appreciate in detail the reason why the STIs in Equations (20)–(22) are instrumental for the block-wise transversality, we will consider the case of Π ˜ μ ν ( 2 ) ( q ) ; the relevant diagrams are enclosed in the blue box of Figure 1.
The diagrams ( a 3 ) and ( a 4 ) are given by
( a 3 ) μ ν ( q ) = g 2 C A k ( k + q ) μ D ( k + q ) D ( k ) I Γ ˜ I Γ ν ( k , k + q , q ) ,
( a 4 ) μ ν ( q ) = g 2 C A g μ ν k D ( k ) ,
where a color factor δ a b is suppressed in both expressions. In addition, for the formal manipulations of integrals, we employ dimensional regularization [256]; to that end, we introduce the short-hand notation
k : = μ 0 ϵ ( 2 π ) d + d d k ,
where d = 4 ϵ is the dimension of the space-time, and μ 0 denotes the ’t Hooft mass.
The contraction of graph ( a 3 ) μ ν ( q ) by q ν triggers the STI satisfied by Γ ˜ Γ ν ( k , k + q , q ) [given by Equation (21)], and we obtain
q ν ( a 3 ) μ ν ( q ) = g 2 C A k ( k + q ) μ D ( k + q ) D ( k ) D 1 ( k ) D 1 ( k + q ) = g 2 C A k ( k + q ) μ D ( k + q ) D ( k ) = g 2 C A q μ k D ( k ) ,
which is precisely the negative of the contraction q ν ( a 4 ) μ ν ( q ) . Hence,
q ν ( a 3 ) μ ν ( q ) + ( a 4 ) μ ν ( q ) = 0 .

3. Schwinger Mechanism in Yang–Mills Theories

The BRST symmetry of the Yang–Mills Lagrangian given in Equation (1) prohibits the inclusion of a mass term of the form m 2 A μ 2 . Moreover, a symmetry-preserving regularization scheme, such as dimensional regularization, prevents the generation of a mass term at any finite order in perturbation theory. Nonetheless, as affirmed four decades ago [94,95,96,97,98,99], the non-perturbative Yang–Mills dynamics endow the gluons with an effective mass, which sets the scale for all dimensionful quantities, and tames the instabilities originating from the infrared divergences of the perturbative expansion ( e.g., Landau pole). In addition, the presence of this mass causes the effective decoupling (screening) of the gluonic modes beyond a “maximum gluon wavelength” [257], and leads to the dynamical suppression of the Gribov copies, see e.g., [16,258,259] and references therein.
The generation of a gluon mass proceeds through the non-perturbative realization of the Schwinger mechanism [127,128]. Even though the technical details associated with the implementation of this mechanism in a four-dimensional non-Abelian setting are particularly elaborate, the general underlying idea is relatively easy to convey.
To that end, consider the dimensionless vacuum polarization Π ( q ) , defined through Π ( q ) = q 2 Π ( q ) , such that
Δ 1 ( q ) = q 2 [ 1 + i Π ( q ) ] .
The Schwinger mechanism is based on the fundamental observation that, if Π ( q ) develops a pole at q 2 = 0 (to be referred to as “massless pole”) then the vector meson (gluon) picks up a mass, regardless of any “prohibition” imposed by the gauge symmetry at the level of the original Lagrangian. Thus, in Euclidean space, the above sequence of ideas leads to
lim q 0 Π ( q ) = m 2 / q 2 lim q 0 Δ 1 ( q ) = lim q 0 ( q 2 + m 2 ) Δ 1 ( 0 ) = m 2 ,
and the gauge boson propagator saturates to a non-zero value at the origin. This effect of infrared saturation of the propagator signifies the generation of a mass, which is identified with the positive residue of the pole.
At this descriptive level, Schwinger’s argument is completely general, making no particular reference to the specific dynamics that would lead to the appearance of the required massless pole inside Π ( q ) . In fact, depending on the particular theory, the field-theoretic circumstances that trigger the crucial sequence captured by Equation (32) may be very distinct, see e.g., [260,261]. In the case of Yang–Mills theories, the origin of the massless poles is purely non-perturbative [159]: the strong dynamics produce scalar composite excitations, which carry color and have vanishing masses. These poles are carried by the fully-dressed vertices of the theory; and since these vertices enter the gluon SDE shown in Figure 1 (upper (lower) panel for the QQ (QB) propagator), the massless poles find their way into the gluon self-energy (or, equivalently, the gluon vacuum polarization). The detailed implementation of this idea has been presented in a series of works [18,93,96,112,116,117,118,159,160,161,166,189,262], and will be summarized in the rest of this section.
Let us focus for now on the conventional three-gluon and ghost–gluon vertices, I Γ α μ ν ( q , r , p ) and I Γ α ( r , p , q ) , respectively, introduced above Equation (23). When the formation of massless poles is triggered, these vertices assume the general form (see Figure 3)
I Γ α μ ν ( q , r , p ) = Γ α μ ν ( q , r , p ) + V α μ ν ( q , r , p ) , I Γ α ( r , p , q ) = Γ α ( r , p , q ) + V α ( r , p , q ) ,
where Γ α μ ν ( q , r , p ) and Γ α ( r , p , q ) are their pole-free components, while V α μ ν ( q , r , p ) and V α ( q , r , p ) contain longitudinally coupled poles, whose special tensorial structure is given by
V α μ ν ( q , r , p ) = q α q 2 C μ ν ( q , r , p ) + r μ r 2 A α ν ( q , r , p ) + p ν p 2 B α μ ( q , r , p ) , V α ( r , p , q ) = q α q 2 C ( r , p , q ) ,
such that
P α α ( q ) P μ μ ( r ) P ν ν ( p ) V α μ ν ( q , r , p ) = 0 , P α α ( q ) V α ( r , p , q ) = 0 .
We emphasize that the reason why V α μ ν ( q , r , p ) and V α ( q , r , p ) are longitudinally coupled may be directly inferred from their special decomposition, shown in Figure 3. In particular, let us denote by I α ( q ) the transition amplitude that connects a gluon with a massless composite scalar, depicted as a gray circle in Figure 3. Since I α ( q ) depends solely on the momentum q, and carries a single Lorentz index, α , its general form is given by I α ( q ) = q α I ( q ) , where I ( q ) is a scalar form factor [117,214]. This observation accounts directly for the form of V α ( q , r , p ) given in Equation (34); to deduce the form of V α μ ν ( q , r , p ) , one must, in addition, appeal to Bose symmetry, which imposes the structures r μ / r 2 and p ν / p 2 in the remaining two channels.
Returning to the SDE of Figure 1, the component V α μ ν ( q , r , p ) will enter in it through graphs ( d 1 ) and ( d 4 ) , while the component V α ( q , r , p ) through graph ( d 3 ). Since V α μ ν ( q , r , p ) has poles for each one of its three momenta, let us point out that only the pole associated with the q-channel, i.e., the channel that carries the momentum entering the gluon propagator is relevant for the Schwinger mechanism that will generate mass for Δ ( q ) . In fact, in the Landau gauge that we employ, the gluon propagators inside the diagrams ( d 1 ) and ( d 4 ) are transverse, leading to a considerable reduction in the number of the form factors of V α μ ν ( q , r , p ) that participate actively, since
P μ μ ( r ) P ν ν ( p ) V α μ ν ( q , r , p ) = q α q 2 P μ μ ( r ) P ν ν ( p ) C μ ν ( q , r , p ) .
Consequently, for the ensuing analysis, one requires only the tensorial decomposition of the component C μ ν ( q , r , p ) in Equation (34), which is given by
C μ ν ( q , r , p ) = C 1 g μ ν + C 2 r μ r ν + C 3 p μ p ν + C 4 r μ p ν + C 5 p μ r ν ,
where C j : = C j ( q , r , p ) . Then, the substitution of Equation (37) into Equation (36), and use of the relation q + p + r = 0 , reveals that only two form factors survive inside ( d 1 ) and ( d 4 ) , namely
P μ μ ( r ) P ν ν ( p ) V α μ ν ( q , r , p ) = q α q 2 P μ μ ( r ) P ν ν ( p ) C 1 g μ ν + C 5 q μ q ν .
Since the main function of the Schwinger mechanism is to make the gluon propagator saturate at the origin, it is important to explore the properties of the structures appearing in Equation (38) near q = 0 . To that end, we expand the r.h.s. of Equation (38), keeping terms at most linear in q. After noticing that the term proportional to C 5 in Equation (38) is of order O ( q 2 ) , we end up with a single relevant form factor associated with V α μ ν ( q , r , p ) , namely C 1 ( q , r , p ) , which survives the q 0 limit of graphs ( d 1 ) and ( d 4 ) . As for V α ( r , p , q ) , its unique component, C ( q , r , p ) , enters directly in ( d 3 ).
The continuation of this analysis entails the Taylor expansion of C 1 ( q , r , p ) and C ( r , p , q ) around q = 0 . In carrying out this expansion, one employs the following two key relations,
C 1 ( 0 , r , r ) = 0 , C ( r , r , 0 ) = 0 .
The first one follows directly from the Bose symmetry of the three-gluon vertex, which implies that C 1 ( q , r , p ) = C 1 ( q , p , r ) ; as we will see in Section 10, it may also be derived in a completely independent way from the fundamental STIs satisfied by the three-gluon vertex. The justification of the second relation in Equation (39) is less straightforward; its derivation, presented in Appendix A, relies on the BQI [14,211] linking the conventional ghost–gluon vertex, I Γ α ( r , p , q ) , with its background counterpart, I Γ ˜ I Γ α ( r , p , q ) .
Thus, after taking Equation (39) into account, the Taylor expansion of C 1 ( q , r , p ) and C ( r , p , q ) around q = 0 yields
lim q 0 C 1 ( q , r , p ) = 2 ( q · r ) C ( r ) + , lim q 0 C ( r , p , q ) = 2 ( q · r ) C ( r ) + ,
with
C ( r ) : = C 1 ( q , r , p ) p 2 q = 0 , C ( r ) : = C ( r , p , q ) p 2 q = 0 .
The functions C ( r ) and C ( r ) are of central importance for the rest of this review. In particular, there are three key points related to them that will be elucidated in detail in what follows:
1.
C ( r ) and C ( r ) are the BS amplitudes describing the formation of gluon–gluon and ghost–anti-ghost colored composite bound states, respectively, see Section 4.
2.
The gluon mass is determined by certain integrals that involve C ( r ) and C ( r ) , given explicitly in Section 5.
3.
C ( r ) and C ( r ) lead to smoking-gun displacements of the WIs. In fact, the displacement induced by C ( r ) , has been confirmed by lattice QCD, by combining judiciously the results of several lattice simulations, see Section 5.2.
We emphasize that the BFM vertices develop poles in exactly the same way as their conventional counterparts. In particular, the main relations Equations (33), (34), (39) and (41) remain valid, with the only modification that all quantities carry hats or tildes; these BFM vertices will be used extensively in Section 5. Note that the conventional and background vertices, including their pole content, are related through appropriate BQIs, see e.g., Equations (A3) and (A6).
We end this section by commenting briefly on the implementation of the Schwinger mechanism away from the Landau gauge, i.e., when the gluon propagator is given by
Δ μ ν ( q , ξ ) = Δ ( q , ξ ) P μ ν ( q ) + ξ q μ q ν / q 4 , ξ 0 ;
for further details, the reader is referred to [189].
(ı) The massless poles remain longitudinally coupled for every value of ξ , i.e., the form of V α μ ν ( q , r , p ) and V α ( r , p , q ) given in Equation (34) persists, with the only difference that the form factors comprising C μ ν ( q , r , p ) , A α ν ( q , r , p ) , B α μ ( q , r , p ) , and C ( r , p , q ) depend in general on ξ . Indeed, as explained right below Equation (35), the longitudinal nature of the poles is dictated solely by Lorentz invariance, which forces the transition amplitude I α ( q , ξ ) to assume the form I α ( q ) = q α I ( q , ξ ) ; clearly, this fundamental argument holds for every ξ .
(ii) Since the gluon propagators entering the graphs ( d 1 ) and ( d 4 ) of Figure 1 are now given by Equation (42), the l.h.s. of Equation (36) becomes Δ μ μ ( r , ξ ) Δ ν ν ( p , ξ ) V α μ ν ( q , r , p ) , and, as a result, the terms in Equation (34) proportional to p ν / p 2 and r μ / r 2 are not fully annihilated. Note, however, that the presence of poles in p 2 k 2 and r 2 ( k + q ) 2 poses no problem, given that one integrates over the loop momentum k. Similar observations hold for the BSE discussed in the next section, which acquires a more complicated form, involving not only the C ( r ) and C ( r ) , but also additional form factors [189].
(iii) A general property of the massless excitations that trigger the Schwinger mechanism is that they do not induce divergences to physical amplitudes; their contributions are completely vanishing, or, at most, finite [260,261]. As was shown recently in [93], in Landau gauge QCD this property hinges on the validity of Equations (35) and (39). Away from the Landau gauge, Equation (39) persists, because its validity relies on Bose symmetry [189]. However, in Equation (35) the substitution P μ ν ( q ) Δ μ ν ( q , ξ ) must be carried out for all projectors; as a result, the r.h.s. no longer vanishes, but includes ξ -dependent longitudinal contributions. Even though this issue has not been addressed in the literature, the longitudinal nature of the additional terms heralds their cancellation through the same general mechanism that renders physical amplitudes ξ -independent.

4. Dynamical formation of Massless Poles

One crucial aspect of the implementation of the Schwinger mechanism in a Yang–Mills context is that the poles that comprise the components V α μ ν ( q , r , p ) and V α ( q , r , p ) in Equation (34) are not introduced by hand; rather, they are generated dynamically, as massless composite excitations that carry color. In fact, this subtle process is controlled by a system of coupled linear BSEs for the functions C ( r ) and C ( r ) , which play the role of the BS amplitudes for generating composite massless scalars out of two gluons and a ghost–anti-ghost pair, respectively.
The starting points for the derivations of the aforementioned BSEs are the SDEs for I Γ α μ ν ( q , r , p ) and I Γ α ( r , p , q ) , shown diagrammatically in Figure 4, and given by [124]
I Γ α μ ν = Γ 0 α μ ν λ k I Γ α β γ Δ β ρ Δ γ σ K 11 μ ν σ ρ + 2 λ k I Γ α D D K 12 μ ν , I Γ α = Γ 0 α λ k I Γ α β γ Δ β ρ Δ γ σ K 21 σ ρ λ k I Γ α D D K 22 ,
where
λ : = i g 2 C A / 2 ,
and the tree-level expressions for the vertices I Γ α μ ν and I Γ α are given by
Γ 0 α μ ν ( q , r , p ) = ( q r ) ν g α μ + ( r p ) α g μ ν + ( p q ) μ g ν α , Γ 0 α ( r , p , q ) = r α .
Note that, for compactness, all momentum arguments have been suppressed; they may be easily restored by appealing to Figure 4.
The following steps are subsequently implemented:
  • Substitute into both sides of Equation (43) the expressions for the fully-dressed vertices given in Equation (33).
  • In order to exploit Equation (38), multiply the first equation by the factor P μ μ ( r ) P ν μ ( p ) .
  • Take the limit of the system as q 0 : this activates Equation (40) and introduces the functions C ( r ) and C ( r ) .
  • Isolate the tensor structures proportional to q α , and match the terms on both sides.
  • Employ the “one-particle exchange” approximation for the kernels K i j , to be denoted by K i j 0 , shown in Figure 5.
Thus, we arrive at a system of homogeneous equations involving C ( r ) and C ( r ) ,
C ( r ) = λ 3 k C ( k ) Δ 2 ( k ) P ρ σ ( k ) P μ ν ( r ) K ˜ 11 μ ν σ ρ + 2 λ 3 k C ( k ) D 2 ( k ) P μ ν ( r ) K ˜ 12 μ ν , C ( r ) = λ k C ( k ) Δ 2 ( k ) P σ ρ ( k ) K ˜ 21 σ ρ λ k C ( k ) D 2 ( k ) K ˜ 22 ,
where K ˜ i j : = ( r · k / r 2 ) K i j 0 ( r , r , k , k ) ; the system is diagrammatically depicted in Figure 6.
Before turning to the numerical analysis, the BSE system must be passed to the Euclidean space, following standard conversion rules. In doing so, we note that the integral measure is modified according to d 4 k i d 4 k E ; this extra factor of i combines with the λ defined in Equation (44) to give real expressions.
As announced, the system of coupled equations given in Equation (46) represents the BSEs that govern the formation of massless colored bound states out of two gluons and a ghost–anti-ghost pair. The functions C ( r ) and C ( r ) are the corresponding BS amplitudes; finding nontrivial solutions for them, i.e., something other than C ( r ) = C ( r ) = 0 identically, is crucial for the implementation of the Schwinger mechanism.
The equations in Equation (46) are linear and homogeneous in the unknown functions. There are two main consequences arising from this fact. First, the numerical solution of the system will be reduced to an eigenvalue problem. Second, the overall scale of the solutions is undetermined, since the multiplication of a given solution by an arbitrary real constant produces another solution (The ambiguity originates from considering only leading terms in the expansion around q = 0 , and may be resolved if further orders in q are kept, see e.g., [219,263,264]).
It turns out that the condition for obtaining nontrivial solutions, when expressed in terms of the strong coupling, α s : = g 2 / 4 π , states that they exist for α s = 0.63 when the renormalization point μ = 4.3 GeV. The solutions obtained when α s acquires this special value are shown in Figure 7; they have undergone scale fixing (The scale was fixed by requiring the best possible matching with the result obtained for C ( r ) from the WI displacement, see Section 12), and are denoted by C 🟉 ( r ) and C 🟉 ( r ) . Observe that C 🟉 ( r ) is significantly larger in magnitude than C 🟉 ( r ) , implying that the three-gluon vertex accounts for the bulk of the gluon mass, as originally claimed in [215].
It is important to compare the value of α s = 0.63 , imposed by the BSE eigenvalue, with the expected value for α s for the renormalization scheme employed: within the asymmetric momentum subtraction (MOM) scheme (see Appendix B), we have that α s = 0.27 [71]. This numerical discrepancy in the values of α s is clearly an artifact of the truncation employed, and concretely of the approximation of the kernels K i j by their one-particle exchange diagrams, K i j 0 . A preliminary analysis reveals that mild modifications of the kernels K i j lead to considerable variations in the value of α s , but leave the form of the solutions for C 🟉 ( r ) and C 🟉 ( r ) practically unaltered. This observation suggests that, while a more complete knowledge of the BSE kernels is required in order to bring α s closer to its MOM value, the solutions obtained with the present approximations should be considered as particularly stable.

5. Generation of the Gluon Mass

We next demonstrate in detail how the presence of the massless poles in the vertices that enter the SDE of the gluon propagator generate a gluon mass.
Since the fundamental STIs of the theory remain intact under the action of the Schwinger mechanism, Equations (7) and (8) remain valid, and the mass term m 2 = Δ 1 ( 0 ) will appear in the transverse combination Δ 1 ( 0 ) P μ ν ( q ) . However, the determination of the mass proportional to g μ ν exposes an entirely different array of principles compared to the corresponding computation for the q μ q ν / q 2 component.
The calculation with respect to the q μ q ν / q 2 component is rather direct; since the massless poles in the vertices are themselves longitudinally coupled, their contribution to the q μ q ν / q 2 component of Π μ ν ( q ) is easily worked out, as will be illustrated in Section 5.1. In contrast, the emergence of a mass proportional to g μ ν is intimately connected with a powerful relation, known as seagull identity [113,166], which in the absence of the Schwinger mechanism would enforce the masslessness of the propagator, as will be discussed in Section 5.2. In fact, one main conceptual difference between the two approaches is that in the g μ ν case, the use of the PT-BFM-based version of the SDE given in Equation (19) is crucial for the emergence of the correct result.
In order to simplify the technical aspects of the calculation without compromising its conceptual content, we will determine the contribution to the gluon mass due to the pole in the ghost–gluon vertex, namely V α ( r , p , q ) in the case of I Γ α ( r , p , q ) , and V ˜ α ( r , p , q ) in the case of I Γ ˜ α ( r , p , q ) . To that end, we will focus on the subset of self-energy graphs containing only ghost loops, i.e., graph ( d 3 ) in the case of Π μ ν ( q ) , and graphs ( a 3 ) and ( a 4 ) in the case of Π ˜ μ ν ( q ) , shown in the upper and lower row of Figure 1, respectively.

5.1. Gluon Mass from the q μ q ν Component

Let us calculate the contribution to the gluon mass stemming from the ghost loop, i.e., the diagram ( d 3 ) of Figure 1, which, for general values of q, reads
( d 3 ) μ ν ( q ) = g 2 C A k ( k + q ) μ D ( k + q ) D ( k ) I Γ ν ( k , k + q , q ) .
To isolate the q μ q ν / q 2 component of Equation (47) at the origin, we first decompose the full vertex I Γ ν ( k , k + q , q ) as in Equations (33) and (34), and drop directly the pole-free part since it does not contribute at q = 0 . Then, denoting by ( d 3 V ) μ ν ( q ) the contribution of V ν ( k , k + q , q ) to ( d 3 ) μ ν ( q ) , we obtain
( d 3 V ) μ ν ( q ) = g 2 C A q ν q 2 k ( k + q ) μ D ( k + q ) D ( k ) C ( k , k + q , q ) .
Next, a Taylor expansion around q = 0 , using Equations (39) and (40), yields
( d 3 V ) μ ν ( q ) = 2 g 2 C A q ν q ρ q 2 k k μ k ρ D 2 ( k ) C ( k ) .
Evidently, the integral above can only be proportional to g μ ρ , such that
( d 3 V ) μ ν ( q ) = 2 g 2 C A d q μ q ν q 2 k k 2 D 2 ( k ) C ( k ) ,
where the tensor structure q μ q ν / q 2 is already isolated.
Then, let us denote by Δ gh 1 ( 0 ) the contribution to the mass originating in the q μ q ν / q 2 of the ghost loop. Noting that the contribution of ( d 3 V ) μ ν ( q ) to the propagator is i times the negative of its q μ q ν / q 2 form factor, we obtain that
Δ gh 1 ( 0 ) = 4 λ d k k 2 D 2 ( k ) C ( k ) .
At this point, we set d = 4 and renormalize Equation (51). This leads to the appearance of the finite renormalization constant of the ghost–gluon vertex, Z 1 .
Next, we express the result in terms of the ghost dressing function F, pass to Euclidean space, and employ hyperspherical coordinates, to obtain the final expression
Δ gh 1 ( 0 ) = λ ^ Z 1 0 d y F 2 ( y ) C ( y ) ,
where λ ^ : = C A α s / 8 π .
The derivation of the contributions from the diagrams ( d 1 ) and ( d 4 ) proceeds in a completely analogous way, but is algebraically more involved, see [166] for details.
It is instructive to consider how the result of Equation (52) emerges in the context of Equation (19). To this end, we consider the ghost block Π ˜ μ ν ( 2 ) ( q ) of Figure 1, whose diagrams have the expressions given in Equation (27); clearly, only diagram ( a 3 ) μ ν ( q ) can contribute to the q μ q ν component of Π ˜ μ ν ( 2 ) ( q ) .
Then, we decompose I Γ ˜ α ( r , p , q ) in complete analogy with Equations (33) and (34), i.e.,
I Γ ˜ I Γ α ( r , p , q ) = Γ ˜ α ( r , p , q ) + q α q 2 C ˜ ( r , p , q ) ,
and expand the ( a 3 ) μ ν ( q ) of Equation (27) around q = 0 , isolating its q μ q ν / q 2 component. These steps eventually lead to
Δ ˜ gh 1 ( 0 ) = 4 λ d k k 2 D 2 ( k ) C ˜ ( k ) ,
where C ˜ ( q ) is defined in the exact same way as C ( q ) , namely through Equation (41) but with tildes over all relevant quantities. It is now easy to establish that Equation (54) is completely equivalent to Equation (51), simply by multiplying both of its sides by Z 1 F ( 0 ) , and then using Equation (A4) on the r.h.s. and Equations (19) and (18) on the l.h.s.
Hence, when the mass is computed through the q μ q ν / q 2 component of the self-energy, the contributions originating from the ghost diagrams of either the BQ or the QQ propagator furnish the same result. The same is not true for the calculation through the g μ ν component, since the ghost diagram ( d 3 ) μ ν of the QQ propagator is not by itself transverse, and a meaningful analysis is preferably carried out within the BFM.

5.2. Gluon Mass from the g μ ν Component: Seagull Identity and Ward Identity Displacement

The fact that the activation of the Schwinger mechanism is crucial for the self-consistent generation of a gluon mass may be best appreciated in conjunction with the so-called seagull identity [113,166]. The content of this identity is that
k k 2 f ( k ) k 2 + d 2 k f ( k ) = 0 ,
for functions f ( k ) that satisfy Wilson’s criterion [265]; the cases of physical interest are f ( k ) = Δ ( k ) , D ( k ) . The general demonstration of the validity of Equation (55) has been given in [166]; for a detailed discussion of how Equation (55) prevents the photon from acquiring a mass in scalar electrodynamics, see [18].
What is so special about Equation (55) is that, within the PT-BFM formalism, the l.h.s. of Equation (55) coincides with the contributions of loop diagrams to the g μ ν component of the gluon mass. Therefore, Equation (55) enforces the non-perturbative masslessness of the gluon in the absence of the Schwinger mechanism: even if a massive gluon propagator (made “massive” through a procedure other than the Schwinger mechanism) were to be substituted inside Equation (55), one would obtain zero as a contribution to the gluon mass! For example, the simple choice f = ( k 2 m 2 ) 1 , reduces the l.h.s of Equation (55) to (dimensionally regularized) textbook integrals, which add up to give precisely zero [18].
In order to appreciate in some detail how the seagull identity prevents the g μ ν component of the propagator from acquiring a mass in the absence of the Schwinger mechanism, let us consider once again the ghost block Π ˜ μ ν ( 2 ) ( q ) of Figure 1; now both graphs, ( a 3 ) and ( a 4 ), contribute to the g μ ν component.
Let us assume that the Schwinger mechanism is turned off; at the level of the Bcc vertex this means that V ^ α ( r , p , q ) vanishes identically, and I Γ ˜ I Γ α ( r , p , q ) = Γ ˜ Γ α ( r , p , q ) . Consequently, Γ ˜ Γ α ( r , p , q ) saturates the STI of Equation (21),
q α Γ ˜ Γ α ( r , p , q ) = D 1 ( p ) D 1 ( r ) .
Since the form factors of the vertex Γ ˜ Γ α ( r , p , q ) do not contain any poles, the derivation from Equation (56) of the corresponding WI proceeds in the standard textbook way: both sides of Equation (56) undergo a Taylor expansion around q = 0 , and terms at most linear in q are retained. Thus, one arrives at the simple QED-like WI
Γ ˜ α ( r , r , 0 ) = D 1 ( r ) r α D 2 ( r ) Γ ˜ α ( r , r , 0 ) = 2 r α D ( r ) r 2 .
We now compute the g μ ν component of Π ˜ μ ν ( 2 ) ( q ) at q = 0 , or, equivalently, Δ ˜ gh 1 ( 0 ) . From Equation (27), we see that ( a 4 ) μ ν is proportional to g μ ν in its entirety. On the other hand, ( a 3 ) μ ν ( q ) contains both g μ ν and q μ q ν components; however, the latter vanishes in the limit q 0 if the vertex is pole-free. Then, it is straightforward to show that, as q 0 ,
Δ ˜ gh 1 ( 0 ) = 2 λ d k k μ D 2 ( k ) Γ ˜ μ ( k , k , 0 ) + d k D ( k ) .
At this point, employing the WI of Equation (57) (with r k ), we get
Δ ˜ gh 1 ( 0 ) = 4 λ d k k 2 D 1 ( k ) k 2 + d 2 k D ( k ) seagull identity = 0 .
Hence, the WI satisfied by the vertex in the absence of the Schwinger mechanism triggers the seagull identity, which, in turn, enforces the masslessness of the propagator.
When the Schwinger mechanism is activated, the STIs that are satisfied by the vertices of the theory retain their original forms but are resolved through the nontrivial participation of the terms containing the massless poles [96,112,159,160,161,166,262,266]. In particular, the full vertex I Γ ˜ I Γ α ( r , p , q ) precisely satisfies Equation (21), namely
q α I Γ ˜ I Γ α ( r , p , q ) = q α Γ ˜ α ( r , p , q ) + C ˜ ( r , p , q ) = D 1 ( p ) D 1 ( r ) .
Notice in particular that the contraction of I Γ ˜ I Γ α ( r , p , q ) by q α cancels the massless pole in q 2 , leading to a completely pole-free result. Therefore, the WI obeyed by Γ ˜ α ( r , p , q ) may be derived as before, through a standard Taylor expansion, leading to
q α Γ ˜ α ( r , r , 0 ) = C ˜ ( r , r , 0 ) + q α D 1 ( r ) r α C ˜ ( r , p , q ) q α q = 0 .
Evidently, the unique zeroth-order contribution appearing in Equation (61), namely C ˜ ( r , r , 0 ) , must vanish,
C ˜ ( r , r , 0 ) = 0 .
Note that this particular property may be independently derived from the antisymmetry of C ˜ ( r , p , q ) under r p , C ˜ ( r , p , q ) = C ˜ ( p , r , q ) , which is a consequence imposed by the ghost–anti-ghost symmetry of the B ( q ) c ¯ ( r ) c ( p ) vertex. The above result, together with Equation (A3), is used to prove Equation (39) in Appendix A.
Thus, Equation (61) becomes
q α Γ ˜ α ( r , r , 0 ) = q α D 1 ( r ) r α 2 r α C ˜ ( r ) , C ˜ ( r ) : = C ˜ ( r , p , q ) p 2 q = 0 ,
and the matching of the terms linear in q yields the WI
Γ ˜ α ( r , r , 0 ) = D 1 ( r ) r α 2 r α C ˜ ( r ) WI displacement .
Comparing Equations (57) and (64), it becomes clear that the Schwinger mechanism induces a characteristic displacement to the WIs that are satisfied by the pole-free parts of the vertices [166].
Returning to Equation (58), but now substituting in it the displaced version of Equation (57), namely
D 2 ( k ) Γ ˜ μ ( k , k , 0 ) = 2 k μ D ( k ) k 2 + D 2 ( k ) C ˜ ( k ) .
When Equation (65) is substituted into Equation (58), the first term of its r.h.s. triggers the seagull identity and vanishes, exactly as before; however, the second term survives, precisely furnishing the result given in Equation (54).
Completely analogous procedures may be applied to the remaining two blocks, Π ˜ μ ν ( 1 ) ( q ) and Π ˜ μ ν ( 3 ) ( q ) , by exploiting the Abelian STIs of Equations (20) and (22), respectively [161].

6. Renormalization Group Invariant Interaction Strength

The PT-BFM formalism provides the natural framework for the construction of the RGI version of the naive one-gluon exchange interaction.
To fix the ideas, recall that in QED, the one-photon exchange interaction, defined as α Δ A ( q ) , where α : = e 2 / 4 π is the hyper-fine structure constant and Δ A ( q ) the photon propagator, is an RGI combination, by virtue of the relation Z e = Z A 1 / 2 ; see comments following Equation (15). Moreover, this particular combination is universal (process-independent) because it may be identified within any two-to-two scattering process, regardless of the nature of the initial and final states (electrons, muons, taus, etc). Instead, in QCD, the corresponding combination α s Δ ( q ) is (trivially) universal but not RGI. When the vertices that connect the gluon to the external particles are “dressed” ( Γ 0 Γ ), the combination Γ α s Δ Γ becomes RGI; however, it is no longer process-independent, because the vertices Γ contain information on the characteristics of the external particles, e.g., the Γ is not the same if the external particles are quarks or gluons. This apparent conundrum may be resolved by resorting to the PT, which reconciles harmoniously the notions of RGI and process independence.
Within the PT framework, the starting point of the construction involves “on-shell” processes [14,96,100,193,194], such as those depicted in Figure 8. The fundamental observation is that the dressed vertices appearing there contain propagator-like contributions, which may be unambiguously identified by means of a well-defined diagrammatic procedure. After discarding terms that vanish on the shell, the contributions extracted from a vertex have a two-fold effect: (i) the genuine vertex contributions left behind form a new vertex, Γ ˜ , which satisfies Abelian STIs, and (ii) when the propagator-like pieces from both vertices are allotted to the conventional propagator, Δ μ ν ( q ) , the resulting effective propagator, Δ ^ μ ν ( q ) , captures all RG logarithms associated with the running of the coupling; for example, at one loop and for large q 2 , one has
Δ ^ 1 ( q ) q 2 1 + b g 2 ln ( q 2 / μ 2 ) ,
where b = 11 C A / 48 π 2 is the first coefficient of the Yang–Mills β function. We emphasize that the PT construction goes through all orders in perturbation theory, as well as non-perturbatively, and all key properties of the PT Green’s function persist unaltered [194,195].
The correspondence between the PT and the BFM may be summarized by stating that the PT rearrangement outlined above amounts effectively to replacing the Q-type gluon that is being exchanged (carrying momentum q) by a B-type gluon [193,267,268,269]; external (on-shell) fields are always of the Q-type. Thus, the notation used above for the PT effective Green’s functions (“tildes” and “hats”) corresponds precisely to the BFM notation introduced in Section 2. Note that the formal expression of all PT rearrangements implemented diagrammatically are the BQIs that relate conventional Green’s functions to their BFM counterparts [14]. For example, in the case of the quark–gluon vertex, we have that the vertices Γ μ ( q , k 1 , k 2 ) [with external fields Q μ a ( q ) q b ( k 1 ) q ¯ c ( k 2 ) ] and Γ ˜ μ ( q , k 1 , k 2 ) [ B μ a ( q ) q b ( k 1 ) q ¯ c ( k 2 ) ] are related by the BQI [270]
Γ ˜ μ ( q , k 1 , k 2 ) = [ 1 + G ( q ) ] Γ μ ( q , k 1 , k 2 ) + ,
where the ellipsis denotes terms that vanish on the shell. Similarly, the BQI of Equation (A5), when evaluated on-shell, yields a completely analogous result, to wit,
I Γ ˜ I Γ μ α ρ ( q , k 1 , k 2 ) = [ 1 + G ( q ) ] I Γ μ α ρ ( q , k 1 , k 2 ) + .
It is now clear how the PT gives rise to a process-independent propagator-like component: regardless of the process ( i.e., the type of vertex connecting the internal gluon to the external states), each vertex contributes to the conventional Δ ( q ) a factor of [ 1 + G ( q ) ] 1 , finally leading to the BQI of Equation (11) [16].
The culmination of the above sequence of ideas is reached by noting that, by virtue of Equation (15), the combination
d ^ ( q ) : = α s Δ ^ ( q ) = α s Δ ( q ) [ 1 + G ( q ) ] 2 ,
is RGI: it retains exactly the same form before and after renormalization, and, consequently, does not depend on the renormalization point μ [96]. The quantity d ^ ( q ) has a mass dimension of 2 , and is known in the literature as the “RGI running interaction strength” [16].
The steps leading to the natural appearance of d ^ ( q ) within any given process may be summarized in the case of quark–antiquark, or gluon–gluon scattering.
Consider the S-matrix elements T q q ¯ q q ¯ , for the scattering of a quark and an antiquark, and T g g g g , for the scattering of two gluons. The quark–antiquark scattering is depicted in the left panel of Figure 8. Using the BQI of Equation (11) we obtain
T q q ¯ q q ¯ = g Γ μ ( q , k 1 , k 2 ) Δ ( q ) P μ ν ( q ) g Γ ν ( q , k 3 , k 4 ) = PT g [ 1 + G ( q ) ] 1 Γ ˜ μ ( q , k 1 , k 2 ) Δ ( q ) P μ ν ( q ) g [ 1 + G ( q ) ] 1 Γ ˜ ν ( q , k 3 , k 4 ) = PT Γ ˜ μ ( q , k 1 , k 2 ) g 2 [ 1 + G ( q ) ] 2 Δ ( q ) P μ ν ( q ) Γ ˜ ν ( q , k 3 , k 4 ) = PT Γ ˜ μ ( q , k 1 , k 2 ) g 2 Δ ^ ( q ) 4 π d ^ ( q ) P μ ν ( q ) Γ ˜ ν ( q , k 3 , k 4 ) ,
where we omit color structures.
Similarly, the scattering of two gluons depicted in the right panel of Figure 8, yields
T g g g g = g Γ α μ ρ ( k 1 , q , k 2 ) Δ ( q ) P μ ν ( q ) g Γ β ν σ ( k 3 , q , k 4 ) = PT g [ 1 + G ( q ) ] 1 Γ ˜ α μ ρ ( k 1 , q , k 2 ) Δ ( q ) P μ ν ( q ) g [ 1 + G ( q ) ] 1 Γ ˜ β ν σ ( k 3 , q , k 4 ) = PT Γ ˜ α μ ρ ( k 1 , q , k 2 ) g 2 [ 1 + G ( q ) ] 2 Δ ( q ) P μ ν ( q ) Γ ˜ β ν σ ( k 3 , q , k 4 ) = PT Γ ˜ α μ ρ ( k 1 , q , k 2 ) g 2 Δ ^ ( q ) 4 π d ^ ( k ) P μ ν ( q ) Γ ˜ β ν σ ( k 3 , q , k 4 ) .
Evidently, the same d ^ ( q ) , defined in Equation (69), appears naturally in both Equations (70) and (71): it is, in that sense, a process-independent RGI interaction capturing faithfully the one-gluon exchange dynamics [3,16,20,79,96,129,130,131].
The actual determination of d ^ ( q ) proceeds by means of the second equality in Equation (69), i.e., by combining the standard gluon propagator, Δ ( q ) , together with the function 1 + G ( q ) . In the top left panel of Figure 9 we show lattice data for the conventional gluon propagator from [85] (points) and a physically motivated fit (blue continuous), given by Equation (C11) of [124]. In the top right panel of the same figure, we show the 1 + G ( q ) auxiliary function, which can be computed by contracting Equation (12) with P μ ν ( q ) / 3 (see e.g., [131]), using the results of [228] for the ghost–gluon kernel, H ν μ ( r , p , q ) . Then, in the bottom left panel of Figure 9 we show the d ^ ( q ) that results from combining the fit for Δ ( q ) and the 1 + G ( q ) shown in the top panels of the same figure and using α s = 0.27 [71] and Z 1 = 0.9333 [see Section 8].
From the d ^ ( q ) of Equation (69) one may define the dimensionless RGI interaction [16], I ( q ) ,
I ( q ) : = q 2 d ^ ( q ) .
As explained in [16], this quantity provides the strength required in order to describe ground-state hadron observables using SDEs in the matter sector of the theory. In that sense, I ( q ) bridges a longstanding gap that has existed between non-perturbative continuum QCD and ab initio predictions of basic hadron properties.

7. Three-Gluon Vertex and Its Planar Degeneracy

The three-gluon vertex, I Γ α μ ν ( q , r , p ) , plays a pivotal role in the dynamics of QCD [234], manifesting its non-Abelian nature through the gluon self-interaction. In fact, the most celebrated perturbative feature of QCD, namely asymptotic freedom, hinges on the properties of this particular interaction vertex. Its importance in the non-perturbative domain has led to an intense effort for unveiling its elaborate features [21,28,33,34,35,36,41,50,68,69,71,78,81,86,87,122,172,173,174,175,176,177,178,179,180,271]. Indeed, as we have seen in Section 3 and Section 4, the pole structure of the three-gluon vertex is crucial for the onset of the Schwinger mechanism and the dynamical generation of a gluon mass. Moreover, its pole-free part provides highly nontrivial contributions to the SDEs of several Green’s functions, most notably the gluon propagator (cf. Figure 1), as well as in the Bethe–Salpeter and Faddeev equations that determine the properties of glue balls [235,236,238,239,240] and hybrid mesons [237], respectively.
For general momenta, I Γ α μ ν ( q , r , p ) is a particularly complicated function, comprised by 14 tensor structures and their associated form factors [251]. Fortunately, in the Landau gauge, considerable simplifications take place, making the treatment of the three-gluon vertex less cumbersome. Indeed, in the latter gauge, quantities of interest require only the knowledge of the transversely projected three-gluon vertex [126,174,175,223], Γ ¯ Γ α μ ν ( q , r , p ) , defined as
Γ ¯ Γ α μ ν ( q , r , p ) = I Γ α μ ν ( q , r , p ) P α α ( q ) P μ μ ( r ) P ν ν ( p ) = Γ α μ ν ( q , r , p ) P α α ( q ) P μ μ ( r ) P ν ν ( p ) .
Note that Γ ¯ Γ α μ ν ( q , r , p ) does not contain massless poles, by virtue of Equation (35). Furthermore, Γ ¯ Γ α μ ν ( q , r , p ) can be parameterized in terms of only 4 independent tensor structures, i.e.,
Γ ¯ Γ α μ ν ( q , r , p ) = i = 1 4 Γ ˜ i ( q 2 , r 2 , p 2 ) λ ˜ i α μ ν ( q , r , p ) .
Due to the Bose symmetry of Γ ¯ Γ α μ ν ( q , r , p ) , the λ ˜ i α μ ν ( q , r , p ) can be chosen to be individually Bose symmetric, such that its form factors Γ ˜ i ( q 2 , r 2 , p 2 ) are symmetric under the exchange of any two arguments [86]. In fact, they can only depend on three totally symmetric combinations of momenta.
Quite remarkably, lattice [86,87,88] and continuum [174,175,223] studies alike, have demonstrated that, to a very good level of accuracy, the Γ ˜ i depend exclusively on a single judiciously chosen variable. Specifically, the Γ ˜ i computed on the lattice in [86,87,88] can be parameterized in terms of the special Bose symmetry combination
s 2 = 1 2 q 2 + r 2 + p 2 .
Thus, the Γ ˜ i are the same for any combination of q 2 , r 2 , and p 2 that fulfils Equation (75) for a given value of s 2 . This property has been denominated planar degeneracy, because Equation (75) with fixed s defines a plane, normal to the vector ( 1 , 1 , 1 ) , in the first octant of the coordinate system ( q 2 , r 2 , p 2 ) .
In particular, the form factor Γ ˜ 1 ( q 2 , r 2 , p 2 ) of the classical tensor structure is rather accurately approximated by
Γ ˜ 1 ( q 2 , r 2 , p 2 ) Γ ˜ 1 ( s 2 , s 2 , 0 ) L s g ( s ) .
In the above equation, L s g is the single transverse form factor of the three-gluon vertex in the soft gluon limit [124], and is obtained in lattice simulations as the q = 0 limit of the following totally transverse projection [84]
L s g ( r ) = Γ 0 α μ ν ( q , r , p ) P α α ( q ) P μ μ ( r ) P ν ν ( p ) I Γ α μ ν ( q , r , p ) Γ 0 α μ ν ( q , r , p ) P α α ( q ) P μ μ ( r ) P ν ν ( p ) Γ 0 α μ ν ( q , r , p ) q 0 .
A particular realization of the planar degeneracy property is shown in Figure 10, where we show the classical form factor Γ ˜ 1 ( q 2 , r 2 , p 2 ) , obtained from the lattice simulation of [86]; we consider three different kinematic configurations, characterized by a single momentum. Specifically, the orange stars correspond to the soft-gluon limit, q = 0 , which implies p 2 = r 2 ; the green diamonds denote the symmetric limit, where all of the momenta have the same magnitude, q 2 = p 2 = r 2 ; and the purple circles represent points with p 2 = r 2 and q 2 = 2 r 2 . When plotted against the momentum r, the three configurations of Γ ˜ 1 ( q 2 , r 2 , p 2 ) produce three clearly distinct curves; however, when plotted in terms of the Bose symmetry variable s of Equation (75), they become statistically indistinguishable, manifesting the validity of Equation (76).
In addition to the planar degeneracy property, lattice [84,86,87,88] and continuum [174,175,179,223] results show a clear dominance of the classical form factor Γ ˜ 1 over the remaining ones. Based on these considerations, the special approximation
Γ ¯ Γ α μ ν ( q , r , p ) L sg s Γ ¯ 0 α μ ν ( q , r , p ) ,
has been put forth, where Γ ¯ 0 α μ ν ( q , r , p ) is the tree-level value of Γ ¯ Γ α μ ν ( q , r , p ) , i.e., Equation (73) with Γ α μ ν ( q , r , p ) Γ 0 α μ ν ( q , r , p ) , and the form factor L sg s has been defined in Equation (77). We emphasize that the shape of L s g ( r ) has been very precisely determined through dedicated lattice studies with large-volume simulations [68,71,84,85]. The outcome of this exploration is shown in Figure 11, where we plot the lattice data of [84] for L s g ( r ) , together with a physically motivated fit given by Equation (C12) of [124] (blue continuous curve). The corresponding fitting formula is rather complicated and will not be reported here; note, however, that the simple expression given in Equation (102) captures rather well the qualitative behavior of L sg s .
Equation (78) provides an accurate and exceptionally compact approximation for Γ ¯ Γ α μ ν ( q , r , p ) in general kinematics. This approximation, with the fit for L s g shown in Figure 11, will be used explicitly in Section 8 and Section 11, where the Γ ¯ Γ α μ ν ( q , r , p ) in general kinematics will be needed as input for the determination of other physically important quantities.

8. Ghost Dynamics from Schwinger–Dyson Equations

We next turn our attention to the ghost sector of the theory, whose scrutiny is important for several reasons. First, it has been connected to particular scenarios of color confinement [272,273]. Second, Green’s functions associated with the ghost sector appear as ingredients in the SDEs of several key functions, such as the gluon propagator and the three-gluon vertex [41,50,68,69,71,81,122,172,173,174,175,176,177,178,179,274], affecting their non-perturbative behavior in nontrivial ways, as will be discussed in Section 9. Third, the SDEs governing the ghost sector are simpler than their gluonic counterparts because they are comprised by fewer diagrams; in fact, the SDE of the ghost propagator contains a single diagram, see Figure 12. Fourth, in the Landau gauge, the validity of Taylor’s theorem [207] facilitates considerably the task of renormalization.
Consequently, the SDEs of the ghost sector are an excellent testing ground for ( a ) probing the impact of the gluonic Green’s functions that contribute to them [85]; ( b ) assessing the reliability of truncation schemes [275,276]; and ( c ) testing the agreement between lattice and continuum approaches.
One of the central results of numerous studies in the continuum [21,62,85,112,178,225,227,228,229,230,231,232,233] as well as a variety of lattice simulations [42,47,49,51,56,63,73,79] may be summarized by stating that the ghost propagator, D ( q ) , remains massless, while the corresponding dressing function, F ( q ) , saturates at the origin. As we will discuss in Section 9, the non-perturbative masslessness of the ghost has important implications for the infrared behavior of the gluon propagator and the three-gluon vertex.
In what follows we provide a concrete example of the state-of-the-art SDE analysis of the ghost sector, by solving the coupled system of equations that governs the ghost-dressing function and the ghost–gluon vertex. In order to obtain a closed system of equations, we use lattice results for the gluon propagator, the three-gluon vertex, and the value of the coupling constant in the particular renormalization scheme employed.
The main points of this analysis may be summarized as follows.
(i) We begin by considering the coupled system of SDEs given in Figure 12, which determines the ghost propagator and ghost–gluon vertex. The treatment will be simplified by neglecting the diagram ( d 3 ν ) of Figure 12, thus eliminating the dependence on the ghost–ghost–gluon–gluon vertex, Γ μ σ . This is a particularly robust truncation, because the impact of the neglected diagram on the ghost–gluon vertex has been shown to be less than 2 % [275].
(ii) Note that due to the fully transverse nature of the gluon propagators in the Landau gauge, in conjunction with the fact that various projections need to be implemented during the treatment of this system, the pole parts V of all fully dressed vertices appearing in Figure 12 will be annihilated; thus, we will have throughout the replacement I Γ Γ .
(iii) We proceed by decomposing the pole-free part, Γ ν ( r , q , p ) , of the ghost–gluon vertex into its most general Lorentz structure, namely
Γ ν ( r , q , p ) = r ν B 1 ( r , q , p ) + p ν B 2 ( r , q , p ) ,
whose scalar form factors reduce to B 1 0 = 1 and B 2 0 = 0 at the tree level. Evidently, due to the transversality of the gluon propagator, only the classical tensor r ν , accompanied by the form factor B 1 , will survive in all SDE diagrams of Figure 12.
(iv) The SDE of Figure 12 is given by
F 1 ( r ) = 1 + 2 λ k f ( k , r ) B 1 ( r , k + r , k ) Δ ( k ) D ( k + r ) ,
where λ is given by Equation (44), and we define
f ( k , r ) : = 1 ( r · k ) 2 r 2 k 2 .
(v) Next, we note that the form factor B 1 ( r , q , p ) can be extracted from Γ ν ( r , q , p ) through the projection
B 1 ( r , q , p ) = ε ν Γ ν ( r , q , p ) , ε ν : = p 2 r ν ( r · p ) p ν r 2 p 2 ( r · p ) 2 .
Hence, acting with ε ν on the diagrams in the second line of Figure 12, we obtain
B 1 ( r , q , p ) = 1 λ a ( r , q , p ) b ( r , q , p ) ,
where
a ( r , q , p ) = q α r μ ε ν k D ( k ) D ( k p ) Δ ( k + r ) B 1 ( p k , q , k + r ) B 1 ( k , k p , p ) P α μ ( k + r ) k ν , b ( r , q , p ) = q α r μ ε ν k Δ ( k ) Δ ( k p ) D ( k + r ) B 1 ( k + r , q , p k ) Γ ¯ Γ ν μ α ( p , k , k p ) .
(vi) At this point, we invoke the property of the planar degeneracy of Γ ¯ Γ α μ ν ( q , r , p ) , discussed in Section 7. Employing Equation (78) into the SDE for B 1 , the term b ( r , q , p ) of Equation (84) becomes
b ( r , q , p ) = q α r μ ε ν k Δ ( k ) Δ ( k p ) D ( k + r ) B 1 ( k + r , q , p k ) Γ ¯ Γ ν μ α 0 ( p , k , k p ) L s g ( s ¯ ) ,
with s ¯ 2 = p 2 + k 2 2 ( k · p ) .
We emphasize that although Equation (78) constitutes in general an approximation, there is one particular kinematic limit in which the expression for b ( r , q , p ) given in Equation (85) becomes exact. Specifically, in the soft gluon limit ( p = 0 ), it can be shown exactly that [85]
P μ μ ( k ) P ν ν ( k ) Γ α μ ν ( 0 , k , k ) = 2 L s g ( k ) k α P μ ν ( k ) .
Then, starting from either the general expression for b ( r , q , p ) of Equation (84) and using Equation (86), or the approximate version given by Equation (85), it can easily be shown that the p = 0 limit is the same. As such, the use of Equation (78) yields not only an excellent approximation in general kinematics, but also the exact soft gluon limit for the contribution of the three-gluon vertex to the form factor B 1 .
(vii) Now we consider the renormalization of the coupled system of equations. Since the ghost–gluon vertex is finite in the Landau gauge [207], most SDE treatments [85,224,225,226,227,228] of the ghost sector employ the so-called Taylor renormalization scheme (see Appendix B), defined in such a way that the finite renormalization constant of the ghost–gluon vertex has the exact value Z 1 = 1 [54,60,80,85,207].
However, in order to employ Equation (78) most expeditiously, it is more convenient to renormalize in the so-called asymmetric MOM scheme, defined in Appendix B, because this is precisely the scheme employed in the lattice calculations of L s g [68,71,84,85]. Past this point, we denote by Z ˜ 1 the finite value of the ghost–gluon renormalization constant in the asymmetric MOM scheme. Evidently, Equations (15) and (79) imply that B 1 R = Z ˜ 1 B 1 .
The renormalization of Equations (80) and (83) proceeds by substitution of the unrenormalized quantities by their renormalized counterparts, following Equation (15), and imposing Equation (A8) for F ( μ 2 ) .
Note that, in principle, Z ˜ 1 , may be determined from the relation Z ˜ 1 = Z 3 Z c Z A 1 , imposed by the corresponding STI [277]; however, these renormalization constants are not available to us, given that Green’s functions have been obtained from the lattice. Therefore, Z ˜ 1 is treated as an adjustable parameter, whose value is determined by requiring that the solution of the SDE for F ( q ) reproduces the corresponding lattice data of [73,85] as well as possible.
(viii) Finally, we transform Equations (80) and (83) from Minkowski to the Euclidean space, using standard conversion rules. Note that, once in Euclidean space, we will express the functional dependence of B 1 ( r , q , p ) in terms of the squared momenta of the anti-ghost and gluon legs, r 2 and p 2 , and the angle, θ , between them, i.e., B 1 ( r , q , p ) B 1 ( r 2 , p 2 , θ ) .
The result of these manipulations is that Equations (80) and (83) become
F 1 ( r ) = 1 α s C A Z ˜ 1 2 π 2 0 d k 2 k 2 Δ ( k ) 0 π d ϕ s ϕ 4 × B 1 ( r 2 , k 2 , ϕ ) F ( z ) z B 1 ( μ 2 , k 2 , ϕ ) F ( u ) u ,
and
B 1 ( r 2 , p 2 , θ ) = Z ˜ 1 α s C A Z ˜ 1 8 π 2 a ¯ + 2 b ¯ ,
respectively, with
a ¯ = 1 s θ 0 d k 2 k 2 F ( k ) 0 π d ϕ s ϕ 3 Δ ( z ) z 0 π d ω s ω F ( v ) v B 1 ( k 2 , p 2 , α ) B 1 ( v , z , β ) K a , b ¯ = 1 s θ 0 d k 2 k 2 Δ ( k ) 0 π d ϕ s ϕ 3 F ( z ) z 0 π d ω s ω Δ ( v ) v B 1 ( z , v , β ) L s g ( s ¯ ) K b .
In the above equations, we employ the notation c x : = cos x and s x : = sin x , and define the following variables
r · k : = r k c ϕ , p · k : = p k ( c θ c ϕ + s θ s ϕ c ω ) , z : = r 2 + k 2 + 2 r k c ϕ , u : = μ 2 + k 2 + 2 μ k c ϕ , s ¯ 2 : = ( p 2 + k 2 + v ) / 2 , v : = p 2 + k 2 2 p k ( c θ c ϕ + s θ s ϕ c ω ) , α : = π cos 1 c θ c ϕ + s θ s ϕ c ω , β : = cos 1 k ( p c θ c ϕ + p s θ s ϕ c ω r c ϕ ) + p r c θ k 2 v z .
Finally, the kernels K a and K b are given by
K a = ( c θ c ω s ϕ c ϕ s θ ) k s ϕ ( p c θ + r ) p c θ c ω ( k c ϕ + r ) , K b = c ω k 2 p c ϕ c θ p s θ 2 ( s ϕ 2 s ω 2 4 s ϕ 2 + 1 ) + s ϕ 2 + r s ϕ 2 s θ 2 ( 2 s ϕ 2 + 1 ) k 3 s ϕ 2 r c θ 2 p s θ 2 + p + p s θ 2 + k p 2 s ϕ 2 2 s θ 2 ( p r c θ ) r c θ p + s θ 2 ( r c θ p ) c ϕ p 3 r s θ 2 + s θ s ϕ c θ p r p 2 k 2 ( s ω 2 + s ϕ 2 s ω 2 2 s ϕ 2 ) c ϕ k ( s ω 2 2 ) ( k 2 + p 2 ) + k c ϕ k 2 r c ϕ p 2 r s θ 2 ( s ω 2 2 ) + s ω 2 + k p 2 3 s θ 2 s ϕ 2 s ω 2 2 s θ 2 s ω 2 4 s θ 2 s ϕ 2 + 3 s θ 2 + ( 3 2 s ω 2 ) s ϕ 2 2 .
We are now in a position to solve Equations (87) and (88) numerically. We choose the renormalization point at μ = 4.3 GeV and employ for Δ ( q ) and L s g ( q ) the fits to the lattice data shown in Figure 9 and Figure 11, respectively. Note that for large momenta these fits recover the behaviors dictated by the corresponding anomalous dimensions [124]. For the strong coupling, we use the value α s ( 4.3 GeV ) = 0.27 , determined from the lattice simulations of [71].
Below, we discuss the main results of this analysis:
The value of Z ˜ 1 was obtained by solving the SDE system for various values of this constant until the χ 2 of the comparison between the solution for F ( q ) and the lattice data of [73,85] was minimized. This procedure yields Z ˜ 1 = 0.9333 ± 0.0075 .
In the left panel of Figure 13, we show as a blue continuous line the SDE result for F ( q ) , with the above value of Z ˜ 1 . The result is compared to the lattice data of [73,85], which have been cured from discretization artifacts. As it turns out, the SDE and lattice results for F agree within 1 % .
We next consider the form factor B 1 . In the right panel of Figure 13 we show B 1 ( r 2 , p 2 , θ ) as a surface, for arbitrary values of the magnitudes of the momenta r and p, and for the angle θ formed between them at θ = 2 π / 3 . In the same panel, we highlight as a red dot-dashed curve the soft gluon limit B 1 ( r 2 , 0 , 2 π / 3 ) of the general kinematics B 1 ( r 2 , p 2 , 2 π / 3 ) (The soft gluon limit is approached by taking p 0 in B 1 ( r 2 , p 2 , θ ) ; in the non-perturbative case, this limit is independent of the value of θ ).
The only available SU(3) lattice data for B 1 were obtained in the soft gluon limit [42,43], and have sizable error bars. Furthermore, they have been computed within the Taylor scheme, while in the present work, we used the asymmetric MOM scheme. Nevertheless, we can meaningfully compare our SDE results with those of the lattice, and perform a statistical analysis to assess their agreement.
Specifically, denoting by B 1 T the Taylor scheme value of the form factor B 1 , it can easily be shown that
B 1 ( r 2 , p 2 , θ ) = Z ˜ 1 B 1 T ( r 2 , p 2 , θ ) ,
which allows us to carry out the desired comparison.
Then, we use Equation (90) to compute B 1 T ( r 2 , 0 , θ ) from the B 1 ( r 2 , 0 , 2 π / 3 ) slice (red dot-dashed curve) in the right panel of Figure 13, and compare the result to the lattice data of [42,43] (points) in Figure 14. Evidently, the SDE determination agrees with the lattice results.
In order to quantify this agreement, we next conduct a χ 2 analysis. To this end, we consider only the 22 lattice points r i in the interval r i [ 0.3 , 2.5 ] GeV, where the signal is most pronounced. Then, we compute the χ 2 of the data through
χ j 2 = i [ B 1 lat ( r i 2 , 0 , θ ) g j ( r i ) ] 2 ϵ B 1 ( r i 2 , 0 , θ ) ,
where B 1 lat ( r i 2 , 0 , θ ) are the lattice points shown in Figure 14, ϵ B 1 ( r i 2 , 0 , θ ) are their respective errors, and g j ( r i ) are the three hypotheses that we will compare to the lattice data. Specifically, for the g j we consider the three cases
g j ( r i ) = 1 if j = 1 , SDE with Γ ¯ Γ α μ ν = Γ ¯ 0 α μ ν L s g ( s ) if j = 2 , SDE with Γ ¯ Γ α μ ν = Γ ¯ 0 α μ ν if j = 3 ,
i.e., g 1 is the tree-level value of B 1 , g 2 is the solution of the SDE using Equation (78) for dressing the three-gluon vertex, corresponding to the red dot-dashed curve of Figure 14, and g 3 is the solution of the SDE obtained by setting the three-gluon vertex to the tree level, which amounts to the substitution L s g 1 in Equation (88), and is represented by a green dashed curve in Figure 14.
Then, for each χ j 2 we compute the probability P j that normally distributed errors would yield a χ 2 at least as large as χ j 2 , through
P j = χ j 2 χ PDF 2 ( 22 , x ) d x = Γ ( n r / 2 , χ 2 / 2 ) Γ ( n r / 2 ) n r = 22 χ 2 = χ j 2 .
In the above equation, χ PDF 2 ( n , x ) = x n / 2 1 e x / 2 / [ 2 n / 2 Γ ( n / 2 ) ] denotes the χ 2 probability distribution function with n degrees of freedom, while Γ ( z , x ) is the incomplete Γ function.
The results of the above analyses are collected in Table 2. We note that the case g 1 , i.e., the tree-level value of B 1 , is discarded at the 5.1 σ confidence level. As for case g 3 , it is discarded at the 3.4 σ level. On the other hand, the SDE result with dressed three-gluon vertex, g 2 , is statistically indistinguishable from the lattice data.
Lastly, we point out that for both F and B 1 , we find a good qualitative agreement with various related studies [21,29,178,179,224,226,227,228,278,279], including kinematics other than the soft gluon limit considered in Figure 14.

9. Divergent Ghost Loops and Their Impact on the QCD Green’s Functions

The masslessness of the ghost propagator, discussed in Section 8, has important implications for the infrared behavior of Green’s functions. Specifically, while the saturation of the gluon propagator renders gluon loops infrared finite, ghost loops furnish infrared divergent contributions [172], akin to those encountered in perturbation theory. In this section, we highlight (with two characteristic examples) how the effects of ghost loops manifest themselves at the level of the two- and three-point functions. Specifically, the ghost loops induce the appearance of a moderate maximum in the gluon propagator and are responsible for the zero-crossing and the logarithmic divergence at the origin displayed by the dominant form factors of the three-gluon vertex.
The basic observation at the level of the gluon SDE shown in Figure 1 is that, the ghost loop of ( d 3 ) , due to the masslessness of its ingredients, furnishes “unprotected” logarithms, i.e., terms of the type ln q 2 , which diverge as q 2 0 . Instead, gluonic loops contain infrared finite gluon propagators and, therefore, give rise to contributions that remain finite as q 2 0 , i.e., they may be described in terms of “protected” logarithms of the type ln ( q 2 + m 2 ) .
The circumstances described above may be modeled by
Δ 1 ( q ) = q 2 + m 2 + c 1 q 2 ln q 2 + ρ m 2 Λ 2 f ( q ) + c 2 q 2 ln q 2 Λ 2 ,
where m is the gluon mass, Λ the mass scale of QCD, and c 1 , c 2 , and ρ are constants; note that Δ 1 ( 0 ) = f ( 0 ) = m 2
Differentiating Equation (94) with respect to q 2 , we obtain
d Δ 1 ( q ) d q 2 = d f ( q ) d q 2 + c 2 1 + ln q 2 Λ 2 .
The second term on the r.h.s. of Equation (95) is infrared divergent, and necessarily dominates the behavior of the derivative of the propagator for sufficiently small q. Moreover, the value of the coefficient c 2 can be computed explicitly by expanding the ghost block Π ˜ μ ν ( 2 ) ( q ) of Figure 1 around q = 0 and using Equation (19), which yields
c 2 = α s C A Z ˜ 1 2 F 2 ( 0 ) 48 π .
Therefore, d Δ 1 ( q ) / d q 2 has the asymptotic behavior
lim q 0 d Δ 1 ( q ) d q 2 = α s C A Z ˜ 1 2 F 2 ( 0 ) 48 π ln q 2 Λ 2 ,
which diverges to as q 0 . Now, since the gluon propagator is a decreasing function in the ultraviolet, we have that d Δ 1 ( q ) / d q 2 is positive for large momenta. Therefore, there must exist a special momentum, denoted by q 🟉 , such that [ d Δ ( q ) / d q 2 ] q = q 🟉 = 0 , which corresponds to a maximum of Δ ( q ) (Note that d Δ 1 ( q ) / d q 2 is an increasing function since it is negative in the infrared and positive in the ultraviolet, i.e., d 2 Δ 1 ( q ) / d ( q 2 ) 2 > 0 . Therefore, assuming that d Δ 1 ( q ) / d q 2 only crosses zero once, q = q 🟉 must be a maximum of Δ ( q ) ).
The maximum of Δ ( q ) , predicted by means of the simple arguments presented above, is observed in lattice simulations of the gluon propagator [49,56,85]. In particular, it is clearly visible in Figure 15, where the data from the two largest volume lattice setups of [49] are shown. The red dashed lines represent smooth functions, fitted to each of the data sets, in the window q [ 0 , 0.5 ] GeV. For each of the volumes considered, V = 72 4 (left panel) and V = 80 4 (right panel), the estimate obtained for q 🟉 is q 🟉 = 140 MeV.
It is interesting to observe in passing that the existence of a maximum of Δ ( q ) has an interesting implication on the form of the spectral function of the gluon propagator [280,281,282,283,284,285]. In particular, the standard Källén-Lehmann representation [286,287] states that
Δ ( q ) = 0 d λ 2 ρ ( λ 2 ) q 2 + λ 2 ,
where ρ ( λ 2 ) is the gluon spectral function (with a factor 1 / π absorbed in it). Thus, the differentiation of both sides of Equation (98) with respect to q 2 yields
d Δ ( q ) d q 2 = 0 d λ 2 ρ ( λ 2 ) ( q 2 + λ 2 ) 2 .
Then, from Equation (99) follows that the existence of a maximum for Δ ( q ) at q = q 🟉 leads necessarily to the violation of reflection positivity [11,167,168,171], because the condition
0 d λ 2 ρ ( λ 2 ) ( q 🟉 2 + λ 2 ) 2 = 0 ,
may be fulfilled only if ρ ( λ 2 ) reverses its sign. Note that an analogous argument based on the existence of an inflection point has been presented recently in [8].
Turning to the three-gluon vertex, it is well-known that the corresponding ghost loops induce characteristic features to the form factors associated with its classical (tree-level) tensors. There are two complementary continuum descriptions of the dynamics that determine the behavior of these form factors: ( i ) the SDE of the three-gluon vertex [174,175,176,226], depicted diagrammatically in Figure 16, and ( ii ) the STI of Equation (23) [172], which, in the limit of vanishing gluon momentum, and when the displacement function and the ghost sector are neglected, yields the approximate WI
I Γ α μ ν ( 0 , r , r ) Δ μ ν 1 ( r ) r α ,
which transmits the properties of the propagator derivative to the vertex form factors, as shown schematically in Figure 17.
In the simplified kinematic circumstances where only a single representative momentum is considered, to be denoted by r, the conclusions drawn by either method may be qualitatively described in terms of a simple model, namely
L ( r ) = b 0 + b gl ln r 2 + m 2 Λ 2 + b gh ln r 2 Λ 2 ,
where L ( r ) denotes the particular combination of form factors, such that, at tree level, L 0 ( r ) = 1 , and b 0 , b gl , and b gh are positive constants. The model in Equation (102) encompasses two important cases studied on the lattice [68,69,71,81], namely ( i ) the soft gluon limit, L ( r ) L s g ( r ) , corresponding to the kinematic choice q 0 , p = r , θ : = p r ^ = π , defined in Equation (77), and ( ii )  totally symmetric limit, L ( r ) L s y m ( r ) , corresponding to q 2 = p 2 = r 2 , θ : = q r ^ = q p ^ = r p ^ = 2 π / 3 .
Upon inspection of Equation (102) we note that, as r 0 , the term with the unprotected logarithm will eventually dominate, forcing L ( r ) to reverse its sign (zero crossing), and finally display a logarithmic divergence, L ( 0 ) . Given that, in practice, b gl is considerably larger than b gh , the unprotected logarithm overtakes the protected one rather deep in the infrared: the location of the zero-crossing is at about 160 MeV [71]. Consequently, in the intermediate region of momenta, which is considered relevant for the onset of non-perturbative dynamics, we have L ( r ) < 1 ; this effect is known in the literature as the infrared suppression of the three-gluon vertex.
Most importantly, the special features of infrared suppression, zero-crossing, and logarithmic divergence at the origin have been corroborated through a variety of lattice results [50,68,69,71,72,81,84], as shown, e.g., in Figure 11. The central curve of this figure is presented as the blue line in Figure 18, where the aforementioned characteristics have been explicitly marked for the benefit of the reader. Note the close proximity of the blue curve to the d Δ 1 ( r ) / d r 2 (red dashed line), especially below 1 GeV.
We end this section by pointing out that, in the case of Yang–Mills in d = 3 [28,172,223,288,289,290,291,292,293,294,295,296,297,298,299,300,301,302], the situation is qualitatively similar to the one described above, but the divergences induced due to the masslessness of the ghost are stronger. Specifically, as may be already established at the level of a simple one-loop calculation [302], the first derivative of the gluon propagator diverges at the origin as 1 / q rather than ln q 2 . Consequently, the corresponding effects are significantly enhanced; in particular, the maximum of the gluon propagator is considerably more pronounced, becoming plainly visible on the lattice [53]. Similarly, an abrupt negative divergence is observed in the corresponding vertex form factors [41,82].

10. Ward Identity Displacement of the Three-Gluon Vertex

In complete analogy to the case of the ghost–gluon vertex discussed in Section 5.2, the WI satisfied by the pole-free part of the three-gluon vertex is also displaced in the presence of longitudinally coupled massless poles. Quite importantly, the associated displacement function, C ( r ) , coincides with the BS amplitude that controls the formation of a (colored) scalar bound state with vanishing mass out of a gluon pair. The displacement of the WI circumvents the seagull cancellation involving the gluon propagator [ i.e., f = Δ in Equation (55)], furnishing to the g μ ν component the mass originating from graphs ( d 1 ) and ( d 4 ) in Figure 1. In addition, it permits the indirect determination of the displacement function C ( r ) from the lattice; this is particularly important, given that, by virtue of Equation (35), the lattice “observables” do not perceive directly the presence of the massless poles.
The starting point of the analysis is the STI satisfied by the three-gluon vertex, I Γ α μ ν ( q , r , p ) , given by Equation (23). In order to eliminate the poles in r and p, thus isolating the displacement of the WI originating from the pole in the channel q, we contract that equation with P μ μ ( r ) P ν ν ( p ) . Note that this procedure also eliminates any longitudinal pole terms in the ghost kernels H σ μ ( p , q , r ) and H σ ν ( r , q , p ) ; for the diagrammatic definition of the ghost–gluon kernel, see Figure 2.
Then, we decompose I Γ α μ ν ( q , r , p ) into pole-free and longitudinally coupled massless pole parts, as in Equation (33), and use Equation (38), to obtain
P μ μ ( r ) P ν ν ( p ) q α Γ α μ ν ( q , r , p ) + g μ ν C 1 ( q , r , p ) + q μ q ν C 5 ( q , r , p ) = P μ μ ( r ) P ν ν ( p ) R ν μ ( p , q , r ) ,
where
R ν μ ( p , q , r ) : = F ( q ) Δ 1 ( p ) P ν σ ( p ) H σ μ ( p , q , r ) Δ 1 ( r ) P μ σ ( r ) H σ ν ( r , q , p ) .
At this point, we expand Equation (103) around q = 0 and match coefficients of equal orders. At the zeroth order in this expansion, we immediately obtain that
C 1 ( 0 , r , r ) = 0 ,
in exact analogy to Equation (62). Note that we have arrived once again at the result of Equation (39), but through an entirely different path: while Equation (39) is enforced by the Bose symmetry of the three-gluon vertex, Equation (105) is a direct consequence of the STI that this vertex satisfies.
We next gather the terms in the expansion of Equation (103) that are of first order in q. Evidently, the term C 5 does not contribute to this order. Then, the expansion leads to the appearance of derivatives of the gluon propagator, in analogy to Equation (64), but also of the ghost–gluon kernel. Specifically, we obtain for the WI of the three-gluon vertex and its displacement the expression
L s g ( r ) = F ( 0 ) Z ˜ 1 d Δ 1 ( r ) d r 2 + W ( r ) r 2 Δ 1 ( r ) C ( r ) .
In the above equation, L s g ( r ) is the form factor of the three-gluon vertex defined in Equation (77) and with lattice results shown in Figure 11, while W ( r ) is a particular derivative of the ghost–gluon kernel, namely [124,241]
W ( r ) = 1 3 r 2 P μ ν ( r ) H ν μ ( p , q , r ) q α q = 0 .
For the detailed derivation of Equation (106), we refer to [93,124].
In the following section, we will use Equation (106) to determine the displacement amplitude C ( r ) from lattice inputs. To this end, we must first pass to Euclidean space, where we note that
C E ( r E 2 ) = C ( r ) r 2 = r E 2 ,
with the extra sign originating from the fact that C is a derivative [see Equation (41)]. Then, suppressing the indices “ E ” and solving for C ( r 2 ) , we obtain the central relation
C ( r ) = L s g ( r ) F ( 0 ) W ( r ) r 2 Δ 1 ( r ) + Z ˜ 1 d Δ 1 ( r ) d r 2 .
For the determination of C ( r ) , we use lattice inputs for all the quantities that appear on the r.h.s. of Equation (109), with the exception of the function W ( r ) , which will be computed from the SDE satisfied by the ghost–gluon kernel derived and analyzed in the next section.

11. The Ghost-Gluon Kernel Contribution to the Ward Identity

In this section, we derive the SDE that determines the key function W ( r ) ; the resulting SDE will be solved using lattice inputs for the various quantities entering it. In addition, the infrared behavior of W ( r ) will be analyzed in detail, following an analytic procedure.
Our discussion starts with the SDE of the ghost–gluon kernel, H μ ν ( r , q , p ) , shown diagrammatically in Figure 19, from which W ( r ) can be obtained using Equation (107).
Note that the similarity between the diagrams shown in Figure 19 and those in the bottom panel of Figure 12, depicting the SDE of the ghost–gluon vertex, is a simple reflection of the fundamental STI relating the ghost–gluon kernel with the ghost–gluon vertex,
Γ ν ( r , q , p ) = r μ H μ ν ( r , q , p ) .
Specifically, Equation (110) is preserved by the SDEs of Γ ν ( r , q , p ) and H μ ν ( r , q , p ) ; indeed, contraction of each diagram ( h i μ ν ) of Figure 19 by r μ yields the corresponding diagram ( g i ν ) of Figure 12 (up to a shift of k k r for i = 1 , introduced to simplify certain expressions). Note that, in Figure 19, the diagram corresponding to the ( g 3 ) of Figure 12 has been omitted, for the reason explained in item (i) of Section 8.
It is well-known that, in the Landau gauge, the momentum q of the ghost field in H μ ν ( r , q , p ) factors out of its quantum corrections [1], allowing us to write [124,228,241]
H μ ν ( r , q , p ) = g μ ν + q α K μ ν α ( r , q , p ) ,
where no particular assumptions are made about the structure of the function K μ ν α ( r , q , p ) . Following Equation (107), we differentiate the r.h.s. of Equation (111) with respect to q and subsequently set q = 0 , to obtain
W ( r ) = 1 3 r α P μ ν ( r ) K μ ν α ( r , 0 , r ) .
Lastly, the finite renormalization of W in the asymmetric MOM scheme proceeds through the use of Equation (15), which leads to the appearance of an overall factor of Z ˜ 1 in the equations.
The outcome of the above steps is that W ( r ) can be written as
W ( r ) = W 1 ( r ) + W 2 ( r ) ,
where the W i ( r ) are the contributions originating from the diagrams ( h i μ ν ) in Figure 19, respectively, and read
W 1 ( r ) = λ Z ˜ 1 3 k Δ ( k ) D ( k ) D ( k + r ) ( r · k ) f ( k , r ) B 1 ( k + r , k , r ) B 1 ( k , 0 , k ) , W 2 ( r ) = λ Z ˜ 1 3 k Δ ( k ) Δ ( k + r ) D ( k + r ) B 1 ( k + r , 0 , k r ) I W ( r , k , k + r ) ,
where f ( k , r ) is given by Equation (81), and we define the specific contribution of the three-gluon vertex to the kernel of W ( r 2 ) as
I W ( q , r , p ) : = 1 2 ( q r ) ν Γ ¯ α ν α ( q , r , p ) .
Note that, from Equation (115) and the Bose symmetry of the Γ ¯ Γ α μ ν ( q , r , p ) under the exchange { q , α } { r , μ } , it follows that
I W ( q , r , p ) = I W ( r , q , p ) .
At this point, by capitalizing on the planar degeneracy of Γ ¯ Γ α μ ν ( q , r , p ) discussed in Section 7, we obtain a compact, and yet accurate, approximation for I W . Specifically, using Equation (78), we find
I W ( q , r , p ) I W 0 ( q , r , p ) L s g ( s ) ,
where I W 0 ( q , r , p ) is the tree-level value of I W , given by
I W 0 ( q , r , p ) : = 2 f ( q , r ) p 2 2 q 2 r 2 ( q 2 + r 2 ) ( q · r ) ( q · r ) 2 .
We remark that the approximation given by Equation (117) becomes exact in the limit p = 0 .
Using the above approximation for I W , the contribution W 2 ( r ) reads
W 2 ( r ) = 2 λ Z ˜ 1 3 k Δ ( k ) Δ ( k + r ) D ( k + r ) ( k + r ) 2 B 1 ( k + r , 0 , k r ) f ( k , r ) × 2 r 2 k 2 ( r 2 + k 2 ) ( r · k ) ( r · k ) 2 L s g ( s ^ ) ,
where we now have s ^ 2 = r 2 + k 2 + ( r · k ) .
Lastly, we transform W 1 of Equation (114) and W 2 of Equation (119) to Euclidean space to obtain the final expression to be used for the numerical determination of W ,
W 1 ( r ) = r α s C A Z ˜ 1 12 π 2 0 d k 2 k Δ ( k ) F ( k ) B 1 ( k 2 , k 2 , π ) 0 π d ϕ s ϕ 4 c ϕ F ( z ) z B 1 ( z , r 2 , χ ) , W 2 ( r ) = r α s C A Z ˜ 1 6 π 2 0 d k 2 k 3 Δ ( k ) 0 π d ϕ s ϕ 4 Δ ( z ) B 1 ( z , z , π ) F ( z ) z 2 k r ( 2 + c ϕ 2 ) z c ϕ × L s g r 2 + k 2 + r k c ϕ ,
where z has been defined below Equation (89) and
χ : = cos 1 ( r + k c ϕ ) z .
We emphasize that for the SDEs of both B 1 and W , given by Equations (88) and (120), respectively, we used the same approximation for the three-gluon vertex, namely Equation (78). Therefore, our analyses of B 1 and W are self-consistent, in the sense that the STI in Equation (110) is strictly preserved.
Before embarking on the numerical determination of W ( r ) for the entire range of Euclidean momenta, we discuss the infrared behavior of this function and demonstrate an important self-consistency proof involving C ( r ) .
Specifically, as discussed in Section 9, the L s g ( r ) and d Δ 1 ( r ) / d r 2 that appear in Equation (109) are infrared divergent, due to the massless ghost loops present in their SDEs. Nevertheless, the BSE solutions for the amplitude C ( r ) are all found to be finite at r = 0 , (cf. Figure 7) [117,121,124,215]. Therefore, in order for the WI displacement of Equation (109) to be consistent with the finite C ( 0 ) obtained from BSE solutions, the infrared divergences of the ingredients appearing in Equation (109) must cancel against each other.
Indeed, a careful analysis of diagram ( e 2 ) of Figure 16 yields
lim r 0 L s g ( r ) = α s C A Z ˜ 1 3 F 3 ( 0 ) 96 π ln r 2 μ 2 ,
up to infrared finite terms (We note that the results identical to Equations (97) and (122) for the infrared divergences of d Δ 1 ( r ) / d r 2 and L s g ( r ) , respectively, have been previously derived within the Curci–Ferrari model [180]). Then, assuming that only L s g ( r ) and d Δ 1 / d r 2 diverge, and using the asymptotic form of d Δ ( r ) / d r 2 given in Equation (97) to Equation (109), we find that the divergences do not fully cancel. Therefore, the finiteness of C ( 0 ) demands that the term W ( r ) / r 2 appearing in the WI must be infrared divergent.
Now, it is evident from Equation (120) that W ( r ) vanishes as r 0 . Nevertheless, the ratio W ( r ) / r 2 is found to diverge at the origin. Specifically, expanding Equation (120) around r = 0 , it can be shown that W ( r ) / r 2 has the asymptotic behavior
lim r 0 W ( r ) r 2 = α s C A Z ˜ 1 3 Δ ( 0 ) F 2 ( 0 ) 96 π ln r 2 μ 2 .
Then, combining Equations (97), (122) and (123) we find that the infrared divergences in Equation (109) cancel out exactly, leaving a finite C ( 0 ) , in full agreement with the BSE results.
We finish the discussion of the infrared finiteness of C ( 0 ) with a remark. In the absence of the Schwinger mechanism, i.e., for an identically zero C ( r ) , the infrared divergences of L s g ( r ) , W ( r ) / r 2 , and d Δ 1 ( r ) / d r 2 must also cancel in Equation (109). For instance, this cancellation can be explicitly verified at the one-loop level (in the perturbative realization of Equation (109) F ( 0 ) also diverges, participating in the overall cancellation of infrared divergences), where, evidently, C ( r ) = 0 . In that case, however, the gluon propagator is also massless, causing the gluonic loops that contribute to the functions that enter Equation (109) to also diverge, such that the cancellation occurs among all radiative diagrams. In contrast, in the presence of a gluon mass, the cancellation of the remaining infrared divergences takes place at the level of the ghost loops only, as illustrated diagrammatically in Figure 20.
We now return to the numerical determination of W ( r ) from Equation (120). To this end, we employ the fits to the lattice data of [84] for Δ ( q ) and L s g ( q ) , shown in Figure 9 and Figure 11, respectively, and the SDE solution for F ( q ) is shown in the left panel of Figure 13. All of the fits employed are constructed so as to reproduce the correct ultraviolet behavior of Green’s functions. For the value of the coupling in the asymmetric MOM scheme, defined in Appendix B, we employ g 2 = 4 π α s , with α s ( 4.3 GeV ) = 0.27 , as determined in the lattice study of [71]. Lastly, for B 1 we use the SDE result of Section 8, shown in the right panel of Figure 13, which reproduces accurately the available lattice data for the ghost–gluon vertex.
Using the above ingredients in Equation (120) we obtain the W ( r ) shown as the blue solid curve in Figure 21. The blue band in Figure 21 represents the error estimate on our results; the procedures followed to obtain it are described in detail in [126].

12. Displacement Function from Lattice Inputs

In this section, we determine C ( r ) from the main relation given in Equation (109).
For W ( r ) we use the result shown in Figure 21, together with the curves for L s g ( r ) from Figure 11, Δ ( r ) and d Δ 1 ( r ) / d r 2 from Figure 9 and Figure 18, respectively, and the F ( r ) of Figure 13. The C ( r ) obtained is shown as a black solid curve in the left panel of Figure 22. In the same panel, we show as points the estimates of C ( r ) obtained by using into Equation (109) the lattice data points of [84] directly, rather than a fit. Note that these data points, as well as those used for the propagators [85], have been carefully extrapolated to the continuum, through the methods explained in [73,85]. These methods exploit the H 4 symmetry of the hypercubic lattice, and are quite effective at minimizing discretization artifacts [54,60,62,63,70,71,73,76,80,85]. As a result of this treatment, the systematic errors are expected to be small. To estimate the uncertainty in the resulting C ( r ) , we combine the error estimate of W ( r ) , represented by the blue band in Figure 21, with the statistical error of the lattice data points for L s g ( r ) of [84], and neglect the error in the gluon propagator, which is much smaller than the errors in L s g and W . Then, a conservative error propagation analysis was carried out in [126], which takes into account an observed correlation between the errors in W ( r ) and L s g ( r ) ; the results of the analysis are the error bars shown in Figure 22.
At this point, we quantify the significance of the C ( r ) obtained above, in comparison to the null hypothesis result; evidently, in the absence of the Schwinger mechanism, this latter quantity, to be denoted by C 0 in what follows, vanishes identically, namely C 0 = 0 . To this end, we first compute the χ 2 of our points as
χ 2 = i C ( r i ) C 0 ( r i ) 2 ϵ C ( r i ) 2 ,
i.e., the null hypothesis is taken as the estimator for our data. The sum runs over the n r = 515 indices i such that r i [ 0.3 , 4.3 ] GeV, the interval of momenta for which the systematic error in our calculation of W ( r ) is best known, and ϵ C ( r i ) denotes the error estimate of C ( r i ) . Then we obtain χ 2 = 2 630 , corresponding to χ d . o . f . 2 = 5.11 . The probability P C 0 that our result for C is consistent with the null hypothesis s vanishingly small, given by the formula
P C 0 = χ 2 = 2 630 χ PDF 2 ( 515 , x ) d x = Γ ( n r / 2 , χ 2 / 2 ) Γ ( n r / 2 ) n r = 515 χ 2 = 2 630 = 7.3 × 10 280 .
Naturally, further correlations in the input data, as well as residual systematic errors, may have escaped the analysis leading to the error estimates shown in Figure 7 for C ( r ) . Since P C 0 changes rapidly with χ 2 , these unknown errors can substantially alter its value. As such, Equation (125) is to be understood as meaning that in the absence of additional uncertainties, the null hypothesis C 0 is excluded. Moreover, it is apparent in Figure 22 that even if the errors had been significantly underestimated, the null hypothesis C 0 would still be unlikely. In fact, even if the errors in all data points for C ( r ) were 95 % larger, i.e., nearly doubled, we could still discard C 0 at the 5 σ confidence level.
In the right panel of Figure 22, we compare C ( r ) to the BSE prediction, C 🟉 ( r ) , of Figure 7, shown as a purple dot-dashed curve and corresponding error band. In that panel, we observe an excellent qualitative agreement between the two results. The most noticeable quantitative difference is in the position of the minimum. Specifically, C reaches the minimum value of 0.36 ± 0.11 at r = 1.93 + 0.09 0.06 GeV, while the minimum of C 🟉 is 0.341 ± 0.003 at r = 1.5 ± 0.1 .
Nevertheless, it is clear in the right panel of Figure 22 that the BSE prediction lies within the error estimate of the lattice-derived C ( r ) . In fact, defining a χ 2 measure for the discrepancy between C and C 🟉 as
χ 🟉 2 = i C ( r i ) C 🟉 ( r i ) 2 ϵ C ( r i ) 2 ,
we obtain χ 🟉 2 = 258.5 , which is smaller than the number of degrees of freedom. Then, this value of χ 🟉 2 amounts to a probability of
P C 🟉 = Γ ( n r / 2 , χ 🟉 2 / 2 ) Γ ( n r / 2 ) n r = 515 χ 🟉 2 = 258.5 = 1 2.0 × 10 23 ,
showing that C 🟉 is statistically compatible with the lattice-derived C , with probability extremely near the unit.

13. Conclusions

The gauge sector of QCD is host to a wide array of subtle mechanisms that are of vital importance for the self-consistency and infrared stability of the theory. In the present work, we offered a comprehensive review of the intricate dynamics that account for some of the most prominent infrared phenomena, such as the generation of a gluon mass through the action of the Schwinger mechanism, the non-perturbative masslessness of the ghost, and the characteristic features induced by this particular mass pattern to the form factors of the three-gluon vertex.
The SDEs, supplemented by the judicious use of certain key results from lattice QCD, provide a robust continuum framework for carrying out such a demanding investigation. In fact, the results obtained from the SDEs are increasingly reliable, passing successfully all sorts of tests imposed on them. A particularly impressive, and certainly not isolated, case of such a success has been outlined in detail in Section 6.
Symmetry and dynamics are tightly interwoven; therefore, the information encoded in the STIs and WIs of the theory is particularly decisive for unraveling basic dynamical patterns. A striking manifestation of the profound connection between symmetry and dynamics is provided by the dual role played by the function C ( r ) , i.e., the BS amplitudes of the massless states composed by a pair of gluons, and the quantity that embodies the displacement induced to the WIs by the presence of these states.
In our opinion, the determination of C ( r ) described in Section 12 represents a major success of the entire set of concepts and techniques surrounding the generation of a gluon mass through the action of the Schwinger mechanism. Thus, fifty years after the genesis of QCD, we seem to be closing in on the mechanism that the theory uses for curing the infrared instabilities endemic to perturbation theory. We hope to be able to report further progress in this direction in the near future.

Funding

The authors are supported by the Spanish MICINN grant PID2020-113334GB-I00. M.N.F. acknowledges financial support from Generalitat Valenciana through contract CIAPOS/2021/74. J.P. also acknowledges funding from the regional Prometeo/2019/087 from Generalitat Valenciana.

Data Availability Statement

Not applicable.

Acknowledgments

The authors thank A.C. Aguilar, D. Binosi, D. Ibá nez, J. Pawlowski, C.D. Roberts, and J. Rodríguez-Quintero for the collaborations.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this work:
BFMbackground field method
BQIbackground-quantum identity
BRSTBecchi–Rouet–Stora–Tyutin
BSBethe–Salpeter
BSEBethe–Salpeter equation
EHMemergent hadron mass
MOMmomentum subtraction (renormalization scheme)
PTpinch technique
QCDquantum chromodynamics
QEDquantum electrodynamics
RGIrenormalization group invariant
SDESchwinger–Dyson equation
STISlavnov–Taylor identity
WIWard identity

Appendix A. BQIs for the BSE Amplitudes

In this appendix, we use two BQIs in order to relate the displacement functions ( C and C ) with their BFM counterparts, i.e., C ˜ and C ˜ , respectively.
The ghost–gluon vertices I Γ μ ( r , p , q ) and I Γ ˜ I Γ μ ( r , p , q ) are related via a BQI [14], which reads
I Γ ˜ I Γ μ ( r , p , q ) = 1 + G ( q ) g μ ν + L ( q ) q μ q ν q 2 I Γ ν ( r , p , q ) + F 1 ( p ) p ν K μ ν ( r , q , p ) + r 2 F 1 ( r ) K μ ( r , q , p ) ,
where K μ and K μ ν are two auxiliary functions, shown diagrammatically in Figure A1, while G ( q ) and L ( q ) are the form factors of Λ μ ν ( q ) , defined in Equation (12).
Figure A1. The auxiliary functions K μ ( q , r , p ) and K μ ν ( q , r , p ) in the BQI of Equation (A1).
Figure A1. The auxiliary functions K μ ( q , r , p ) and K μ ν ( q , r , p ) in the BQI of Equation (A1).
Particles 06 00017 g0a1
Next, we decompose the I Γ ˜ I Γ μ ( r , p , q ) and I Γ μ ( r , p , q ) in Equation (A1) into their regular and pole parts, using Equations (33) and (53), respectively. Note that the second and third terms in Equation (A1) do not contain poles in q 2 ; this is so because K μ ν ( r , q , p ) can contain (longitudinally coupled) poles only in the p ν channel, whereas K μ ( r , q , p ) has no external gluon legs (and, hence, no poles).
Then, multiplying Equation (A1) by q 2 we obtain
q μ C ˜ ( r , p , q ) = q μ 1 + G ( q ) + L ( q ) C ( r , p , q ) + O ( q 2 ) .
Setting q = 0 in Equation (A2) and using Equation (18), we find
C ( r , r , 0 ) = Z 1 F ( 0 ) C ˜ ( r , r , 0 ) .
Hence, using Equation (62), we obtain the result in Equation (39).
Then, expanding Equation (A2) to first order in q, using Equation (41) for C ( r , p , q ) and Equation (63) for C ˜ ( r , p , q ) , entails
C ( r ) = Z 1 F ( 0 ) C ˜ ( r ) ,
which is one of the main results of this appendix.
A relation identical to Equation (A4) can be obtained for C ( r ) and its BFM analog, C ˜ ( r ) . The starting point of the derivation is the BQI [14]
I Γ ˜ I Γ α μ ν ( q , r , p ) = 1 + G ( q ) g α ρ + L ( q ) q α q ρ q 2 I Γ ρ μ ν ( q , r , p ) + K ρ ν α ( r , q , p ) P μ ρ ( r ) Δ 1 ( r ) K ρ μ α ( p , q , r ) P ν ρ ( p ) Δ 1 ( p ) ,
where K μ ν α ( r , q , p ) is the function defined in Equation (111).
Then, we note that the only longitudinal poles at q = 0 present in Equation (A5) are those contained in the I Γ α μ ν ( q , r , p ) and I Γ ˜ I Γ α μ ν ( q , r , p ) vertices, with the auxiliary functions K α ν ρ ( q , p , r ) having poles only in the r μ and p ν channels. As such, isolating the q α g μ ν / q 2 pole and expanding around q = 0 , one eventually finds
C ˜ 1 ( 0 , r , r ) = Z 1 1 F 1 ( 0 ) C 1 ( 0 , r , r ) = 0 ,
and
C ( r ) = Z 1 F ( 0 ) C ˜ ( r ) ,
where C ˜ 1 ( q , r , p ) and C ˜ ( r 2 ) are defined in analogy to Equations (37) and (41), and we used Equation (39).

Appendix B. The Asymmetric MOM Scheme

In this appendix, we provide a brief overview of the asymmetric MOM scheme [68,71,84,85,241] that we employ throughout this work.
The set of boundary conditions imposed on the renormalized quantities defines the renormalization scheme employed. Within the MOM schemes [277], propagators assume their tree-level values at the subtraction point μ , namely
Δ R 1 ( μ ) = μ 2 , F R ( μ ) = 1 .
Past this point, the various MOM schemes are differentiated according to the way the three-point functions are renormalized.
In Landau gauge, a common choice of renormalization prescription is the so-called “Taylor scheme” [54,60,80,85,207]. This scheme capitalizes on the Taylor theorem [207], i.e., the observation that the unrenormalized ghost–gluon vertex in the Landau gauge reduces to its tree-level form in the soft-ghost configuration,
I Γ ν ( r , 0 , r ) = r ν .
The Taylor scheme is defined by requiring Equation (A9) to hold after renormalization [54,60,80,85,207]. Using Equation (15), this requirement yields Z 1 = 1 .
Alternatively, in lattice simulations of the three-gluon vertex, it is convenient to impose a renormalization prescription for its classical tensor structure. For example, one may choose the classical form factor to reduce to tree-level in the symmetric point, q 2 = r 2 = p 2 = μ 2 . This condition defines the “symmetric scheme” [68,71,84].
In the present work, the classical form factor of the three-gluon vertex in the soft-gluon limit, which is denoted by L s g ( r ) and defined in Equation (77), plays a key role. Indeed, it is the central ingredient in the approximation of the three-gluon vertex given by Equation (78), which is used in the SDE analysis of the ghost–gluon vertex and kernel in Section 8 and Section 11, respectively. Moreover, L s g ( r ) is one of the inputs necessary for the determination of the displacement amplitude C ( r ) in Section 12, which signals the activation of the Schwinger mechanism. As such, it is convenient to employ throughout the scheme where L s g ( r ) is most readily renormalized in lattice simulations, which is the so-called “asymmetric scheme” [68,71,84,85,241].
The asymmetric MOM scheme is defined by imposing that L s g ( r ) reduces to the tree-level at q 2 = μ 2 , i.e.,
L s g ( μ ) = 1 .
Note that in this scheme the finite renormalization constant of the ghost–gluon vertex is no longer equal to 1 [85,241]. Instead, the special value of Z 1 in the asymmetric scheme is denoted by Z ˜ 1 , and is determined to be Z ˜ 1 = 0.9333 ± 0.0075 [126], at μ = 4.3 GeV, through the SDE analysis discussed in Section 8.

References

  1. Marciano, W.J.; Pagels, H. Quantum Chromodynamics: A Review. Phys. Rep. 1978, 36, 137. [Google Scholar] [CrossRef]
  2. Qin, S.X.; Roberts, C.D. Impressions of the Continuum Bound State Problem in QCD. Chin. Phys. Lett. 2020, 37, 121201. [Google Scholar] [CrossRef]
  3. Roberts, C.D. Empirical Consequences of Emergent Mass. Symmetry 2020, 12, 1468. [Google Scholar] [CrossRef]
  4. Cui, Z.F.; Ding, M.; Gao, F.; Raya, K.; Binosi, D.; Chang, L.; Roberts, C.D.; Rodríguez-Quintero, J.; Schmidt, S.M. Kaon and pion parton distributions. Eur. Phys. J. C 2020, 80, 1064. [Google Scholar] [CrossRef]
  5. Chang, L.; Roberts, C.D. Regarding the Distribution of Glue in the Pion. Chin. Phys. Lett. 2021, 38, 081101. [Google Scholar] [CrossRef]
  6. Cui, Z.F.; Ding, M.; Morgado, J.M.; Raya, K.; Binosi, D.; Chang, L.; De Soto, F.; Roberts, C.D.; Rodríguez-Quintero, J.; Schmidt, S.M. Emergence of pion parton distributions. Phys. Rev. D 2022, 105, L091502. [Google Scholar]
  7. Lu, Y.; Chang, L.; Raya, K.; Roberts, C.D.; Rodríguez-Quintero, J. Proton and pion distribution functions in counterpoint. Phys. Lett. B 2022, 830, 137130. [Google Scholar] [CrossRef]
  8. Ding, M.; Roberts, C.D.; Schmidt, S.M. Emergence of Hadron Mass and Structure. Particles 2023, 6, 57–120. [Google Scholar] [CrossRef]
  9. Roberts, C.D. Origin of the Proton Mass. arXiv 2022, arXiv:2211.09905. [Google Scholar]
  10. Roberts, C.D.; Williams, A.G. Dyson-Schwinger equations and their application to hadronic physics. Prog. Part. Nucl. Phys. 1994, 33, 477–575. [Google Scholar] [CrossRef] [Green Version]
  11. Alkofer, R.; von Smekal, L. The Infrared behavior of QCD Green’s functions: Confinement dynamical symmetry breaking, and hadrons as relativistic bound states. Phys. Rep. 2001, 353, 281. [Google Scholar] [CrossRef]
  12. Fischer, C.S. Infrared properties of QCD from Dyson-Schwinger equations. J. Phys. G 2006, 32, R253–R291. [Google Scholar] [CrossRef]
  13. Roberts, C.D. Hadron Properties and Dyson-Schwinger Equations. Prog. Part. Nucl. Phys. 2008, 61, 50–65. [Google Scholar] [CrossRef] [Green Version]
  14. Binosi, D.; Papavassiliou, J. Pinch Technique: Theory and Applications. Phys. Rep. 2009, 479, 1–152. [Google Scholar] [CrossRef] [Green Version]
  15. Bashir, A.; Chang, L.; Cloet, I.C.; El-Bennich, B.; Liu, Y.X.; Roberts, C.D.; Tandy, P.C. Collective perspective on advances in Dyson-Schwinger Equation QCD. Commun. Theor. Phys. 2012, 58, 79–134. [Google Scholar] [CrossRef]
  16. Binosi, D.; Chang, L.; Papavassiliou, J.; Roberts, C.D. Bridging a gap between continuum-QCD and ab initio predictions of hadron observables. Phys. Lett. 2015, B742, 183–188. [Google Scholar] [CrossRef] [Green Version]
  17. Cloet, I.C.; Roberts, C.D. Explanation and Prediction of Observables using Continuum Strong QCD. Prog. Part. Nucl. Phys. 2014, 77, 1–69. [Google Scholar] [CrossRef] [Green Version]
  18. Aguilar, A.C.; Binosi, D.; Papavassiliou, J. The Gluon Mass Generation Mechanism: A Concise Primer. Front. Phys. (Beijing) 2016, 11, 111203. [Google Scholar] [CrossRef] [Green Version]
  19. Binosi, D.; Chang, L.; Papavassiliou, J.; Qin, S.X.; Roberts, C.D. Symmetry preserving truncations of the gap and Bethe–Salpeter equations. Phys. Rev. 2016, D93, 096010. [Google Scholar] [CrossRef] [Green Version]
  20. Binosi, D.; Mezrag, C.; Papavassiliou, J.; Roberts, C.D.; Rodriguez-Quintero, J. Process-independent strong running coupling. Phys. Rev. 2017, D96, 054026. [Google Scholar] [CrossRef] [Green Version]
  21. Huber, M.Q. Nonperturbative properties of Yang-Mills theories. Phys. Rep. 2020, 879, 1–92. [Google Scholar] [CrossRef]
  22. Pawlowski, J.M.; Litim, D.F.; Nedelko, S.; von Smekal, L. Infrared behavior and fixed points in Landau gauge QCD. Phys. Rev. Lett. 2004, 93, 152002. [Google Scholar] [CrossRef] [PubMed]
  23. Pawlowski, J.M. Aspects of the functional renormalisation group. Ann. Phys. 2007, 322, 2831–2915. [Google Scholar] [CrossRef] [Green Version]
  24. Fischer, C.S.; Maas, A.; Pawlowski, J.M. On the infrared behavior of Landau gauge Yang-Mills theory. Ann. Phys. 2009, 324, 2408–2437. [Google Scholar] [CrossRef] [Green Version]
  25. Carrington, M.E. Renormalization group flow equations connected to the n-particle-irreducible effective action. Phys. Rev. 2013, D87, 045011. [Google Scholar] [CrossRef] [Green Version]
  26. Carrington, M.E.; Fu, W.J.; Pickering, D.; Pulver, J.W. Renormalization group methods and the 2PI effective action. Phys. Rev. D 2015, 91, 025003. [Google Scholar] [CrossRef] [Green Version]
  27. Cyrol, A.K.; Mitter, M.; Pawlowski, J.M.; Strodthoff, N. Nonperturbative quark, gluon, and meson correlators of unquenched QCD. Phys. Rev. 2018, D97, 054006. [Google Scholar] [CrossRef] [Green Version]
  28. Corell, L.; Cyrol, A.K.; Mitter, M.; Pawlowski, J.M.; Strodthoff, N. Correlation functions of three-dimensional Yang-Mills theory from the FRG. SciPost Phys. 2018, 5, 066. [Google Scholar] [CrossRef] [Green Version]
  29. Huber, M.Q. Correlation functions of Landau gauge Yang-Mills theory. Phys. Rev. D 2020, 101, 114009. [Google Scholar] [CrossRef]
  30. Dupuis, N.; Canet, L.; Eichhorn, A.; Metzner, W.; Pawlowski, J.M.; Tissier, M.; Wschebor, N. The nonperturbative functional renormalization group and its applications. Phys. Rep. 2021, 910, 1–114. [Google Scholar] [CrossRef]
  31. Blaizot, J.P.; Pawlowski, J.M.; Reinosa, U. Functional renormalization group and 2PI effective action formalism. Ann. Phys. 2021, 431, 168549. [Google Scholar] [CrossRef]
  32. Mandula, J.; Ogilvie, M. The Gluon Is Massive: A Lattice Calculation of the Gluon Propagator in the Landau Gauge. Phys. Lett. B 1987, 185, 127–132. [Google Scholar] [CrossRef]
  33. Parrinello, C. Exploratory study of the three gluon vertex on the lattice. Phys. Rev. 1994, D50, R4247–R4251. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Alles, B.; Henty, D.; Panagopoulos, H.; Parrinello, C.; Pittori, C.; Richards, D.G. αs from the nonperturbatively renormalised lattice three gluon vertex. Nucl. Phys. 1997, B502, 325–342. [Google Scholar] [CrossRef] [Green Version]
  35. Parrinello, C.; Richards, D.; Alles, B.; Panagopoulos, H.; Pittori, C. Status of alpha-s determinations from the nonperturbatively renormalized three gluon vertex. Nucl. Phys. B Proc. Suppl. 1998, 63, 245–247. [Google Scholar] [CrossRef] [Green Version]
  36. Boucaud, P.; Leroy, J.P.; Micheli, J.; Pene, O.; Roiesnel, C. Lattice calculation of alpha(s) in momentum scheme. J. High Energy Phys. 1998, 10, 017. [Google Scholar] [CrossRef] [Green Version]
  37. Alexandrou, C.; de Forcrand, P.; Follana, E. The gluon propagator without lattice Gribov copies on a finer lattice. Phys. Rev. 2002, D65, 114508. [Google Scholar] [CrossRef] [Green Version]
  38. Bowman, P.O.; Heller, U.M.; Williams, A.G. Lattice quark propagator with staggered quarks in Landau and Laplacian gauges. Phys. Rev. D 2002, 66, 014505. [Google Scholar] [CrossRef] [Green Version]
  39. Skullerud, J.I.; Bowman, P.O.; Kizilersu, A.; Leinweber, D.B.; Williams, A.G. Nonperturbative structure of the quark gluon vertex. J. High Energy Phys. 2003, 4, 047. [Google Scholar] [CrossRef] [Green Version]
  40. Bowman, P.O.; Heller, U.M.; Leinweber, D.B.; Parappilly, M.B.; Williams, A.G. Unquenched gluon propagator in Landau gauge. Phys. Rev. D 2004, 70, 034509. [Google Scholar] [CrossRef] [Green Version]
  41. Cucchieri, A.; Maas, A.; Mendes, T. Exploratory study of three-point Green’s functions in Landau-gauge Yang-Mills theory. Phys. Rev. 2006, D74, 014503. [Google Scholar] [CrossRef] [Green Version]
  42. Ilgenfritz, E.M.; Muller-Preussker, M.; Sternbeck, A.; Schiller, A.; Bogolubsky, I. Landau gauge gluon and ghost propagators from lattice QCD. Braz. J. Phys. 2007, 37, 193–200. [Google Scholar] [CrossRef]
  43. Sternbeck, A. The Infrared Behavior of Lattice QCD Green’s Functions. Ph.D. Thesis, Humboldt-University Berlin, Berlin, Germany, 2006. [Google Scholar]
  44. Furui, S.; Nakajima, H. Unquenched Kogut-Susskind quark propagator in lattice Landau gauge QCD. Phys. Rev. D 2006, 73, 074503. [Google Scholar] [CrossRef] [Green Version]
  45. Bowman, P.O.; Heller, U.M.; Leinweber, D.B.; Parappilly, M.B.; Sternbeck, A.; von Smekal, L.; Williams, A.G.; Zhang, J.b. Scaling behavior and positivity violation of the gluon propagator in full QCD. Phys. Rev. D 2007, 76, 094505. [Google Scholar] [CrossRef] [Green Version]
  46. Kamleh, W.; Bowman, P.O.; Leinweber, D.B.; Williams, A.G.; Zhang, J. Unquenching effects in the quark and gluon propagator. Phys. Rev. 2007, D76, 094501. [Google Scholar] [CrossRef] [Green Version]
  47. Cucchieri, A.; Mendes, T. What’s up with IR gluon and ghost propagators in Landau gauge? A puzzling answer from huge lattices. PoS 2007, LATTICE2007, 297. [Google Scholar] [CrossRef] [Green Version]
  48. Cucchieri, A.; Mendes, T. Constraints on the IR behavior of the gluon propagator in Yang-Mills theories. Phys. Rev. Lett. 2008, 100, 241601. [Google Scholar] [CrossRef] [Green Version]
  49. Bogolubsky, I.; Ilgenfritz, E.; Muller-Preussker, M.; Sternbeck, A. The Landau gauge gluon and ghost propagators in 4D SU(3) gluodynamics in large lattice volumes. PoS 2007, LATTICE2007, 290. [Google Scholar] [CrossRef] [Green Version]
  50. Cucchieri, A.; Maas, A.; Mendes, T. Three-point vertices in Landau-gauge Yang-Mills theory. Phys. Rev. 2008, D77, 094510. [Google Scholar] [CrossRef] [Green Version]
  51. Cucchieri, A.; Mendes, T. Constraints on the IR behavior of the ghost propagator in Yang-Mills theories. Phys. Rev. D 2008, 78, 094503. [Google Scholar] [CrossRef] [Green Version]
  52. Cucchieri, A.; Mendes, T. Landau-gauge propagators in Yang-Mills theories at beta = 0: Massive solution versus conformal scaling. Phys. Rev. 2010, D81, 016005. [Google Scholar] [CrossRef] [Green Version]
  53. Cucchieri, A.; Mendes, T. Numerical test of the Gribov-Zwanziger scenario in Landau gauge. PoS 2009, QCD-TNT09, 026. [Google Scholar] [CrossRef]
  54. Boucaud, P.; De Soto, F.; Leroy, J.P.; Le Yaouanc, A.; Micheli, J.; Pene, O.; Rodriguez-Quintero, J. Ghost-gluon running coupling, power corrections and the determination of Lambda(MS-bar). Phys. Rev. 2009, D79, 014508. [Google Scholar]
  55. Cucchieri, A.; Mendes, T.; Santos, E.M.S. Covariant gauge on the lattice: A New implementation. Phys. Rev. Lett. 2009, 103, 141602. [Google Scholar] [CrossRef] [Green Version]
  56. Bogolubsky, I.; Ilgenfritz, E.; Muller-Preussker, M.; Sternbeck, A. Lattice gluodynamics computation of Landau gauge Green’s functions in the deep infrared. Phys. Lett. 2009, B676, 69–73. [Google Scholar] [CrossRef]
  57. Oliveira, O.; Silva, P. The Lattice infrared Landau gauge gluon propagator: The Infinite volume limit. PoS 2009, LAT2009, 226. [Google Scholar] [CrossRef] [Green Version]
  58. Cucchieri, A.; Mendes, T.; Nakamura, G.M.; Santos, E.M.S. Gluon Propagators in Linear Covariant Gauge. PoS 2010, FACESQCD, 026. [Google Scholar] [CrossRef] [Green Version]
  59. Oliveira, O.; Bicudo, P. Running Gluon Mass from Landau Gauge Lattice QCD Propagator. J. Phys. G 2011, G38, 045003. [Google Scholar] [CrossRef] [Green Version]
  60. Blossier, B.; Boucaud, P.; De soto, F.; Morenas, V.; Gravina, M.; Pene, O.; Rodriguez-Quintero, J. Ghost-gluon coupling, power corrections and Λ MS ¯ from twisted-mass lattice QCD at Nf = 2. Phys. Rev. D 2010, 82, 034510. [Google Scholar] [CrossRef] [Green Version]
  61. Maas, A. Describing gauge bosons at zero and finite temperature. Phys. Rep. 2013, 524, 203–300. [Google Scholar] [CrossRef] [Green Version]
  62. Boucaud, P.; Leroy, J.P.; Yaouanc, A.L.; Micheli, J.; Pene, O.; Rodriguez-Quintero, J. The Infrared Behaviour of the Pure Yang-Mills Green Functions. Few Body Syst. 2012, 53, 387–436. [Google Scholar] [CrossRef] [Green Version]
  63. Ayala, A.; Bashir, A.; Binosi, D.; Cristoforetti, M.; Rodriguez-Quintero, J. Quark flavour effects on gluon and ghost propagators. Phys. Rev. 2012, D86, 074512. [Google Scholar] [CrossRef] [Green Version]
  64. Oliveira, O.; Silva, P.J. The lattice Landau gauge gluon propagator: Lattice spacing and volume dependence. Phys. Rev. 2012, D86, 114513. [Google Scholar] [CrossRef]
  65. Sternbeck, A.; Müller-Preussker, M. Lattice evidence for the family of decoupling solutions of Landau gauge Yang-Mills theory. Phys. Lett. B 2013, 726, 396–403. [Google Scholar] [CrossRef] [Green Version]
  66. Bicudo, P.; Binosi, D.; Cardoso, N.; Oliveira, O.; Silva, P.J. Lattice gluon propagator in renormalizable ξ gauges. Phys. Rev. 2015, D92, 114514. [Google Scholar] [CrossRef] [Green Version]
  67. Duarte, A.G.; Oliveira, O.; Silva, P.J. Lattice Gluon and Ghost Propagators, and the Strong Coupling in Pure SU(3) Yang-Mills Theory: Finite Lattice Spacing and Volume Effects. Phys. Rev. D 2016, 94, 014502. [Google Scholar] [CrossRef] [Green Version]
  68. Athenodorou, A.; Binosi, D.; Boucaud, P.; De Soto, F.; Papavassiliou, J.; Rodriguez-Quintero, J.; Zafeiropoulos, S. On the zero crossing of the three-gluon vertex. Phys. Lett. 2016, B761, 444–449. [Google Scholar] [CrossRef]
  69. Duarte, A.G.; Oliveira, O.; Silva, P.J. Further Evidence For Zero Crossing On The Three Gluon Vertex. Phys. Rev. 2016, D94, 074502. [Google Scholar] [CrossRef] [Green Version]
  70. Oliveira, O.; Kizilersu, A.; Silva, P.J.; Skullerud, J.I.; Sternbeck, A.; Williams, A.G. Lattice Landau gauge quark propagator and the quark-gluon vertex. Acta Phys. Pol. Suppl. 2016, 9, 363–368. [Google Scholar] [CrossRef] [Green Version]
  71. Boucaud, P.; De Soto, F.; Rodríguez-Quintero, J.; Zafeiropoulos, S. Refining the detection of the zero crossing for the three-gluon vertex in symmetric and asymmetric momentum subtraction schemes. Phys. Rev. 2017, D95, 114503. [Google Scholar] [CrossRef] [Green Version]
  72. Sternbeck, A.; Balduf, P.H.; Kizilersu, A.; Oliveira, O.; Silva, P.J.; Skullerud, J.I.; Williams, A.G. Triple-gluon and quark-gluon vertex from lattice QCD in Landau gauge. PoS 2017, LATTICE2016, 349. [Google Scholar] [CrossRef] [Green Version]
  73. Boucaud, P.; De Soto, F.; Raya, K.; Rodríguez-Quintero, J.; Zafeiropoulos, S. Discretization effects on renormalized gauge-field Green’s functions, scale setting, and the gluon mass. Phys. Rev. 2018, D98, 114515. [Google Scholar] [CrossRef] [Green Version]
  74. Cucchieri, A.; Dudal, D.; Mendes, T.; Oliveira, O.; Roelfs, M.; Silva, P.J. Lattice Computation of the Ghost Propagator in Linear Covariant Gauges. PoS 2018, LATTICE2018, 252. [Google Scholar] [CrossRef]
  75. Cucchieri, A.; Dudal, D.; Mendes, T.; Oliveira, O.; Roelfs, M.; Silva, P.J. Faddeev-Popov Matrix in Linear Covariant Gauge: First Results. Phys. Rev. D 2018, 98, 091504. [Google Scholar] [CrossRef] [Green Version]
  76. Oliveira, O.; Silva, P.J.; Skullerud, J.I.; Sternbeck, A. Quark propagator with two flavors of O(a)-improved Wilson fermions. Phys. Rev. D 2019, 99, 094506. [Google Scholar] [CrossRef] [Green Version]
  77. Dudal, D.; Oliveira, O.; Silva, P.J. High precision statistical Landau gauge lattice gluon propagator computation vs. the Gribov–Zwanziger approach. Ann. Phys. 2018, 397, 351–364. [Google Scholar] [CrossRef] [Green Version]
  78. Vujinovic, M.; Mendes, T. Probing the tensor structure of lattice three-gluon vertex in Landau gauge. Phys. Rev. 2019, D99, 034501. [Google Scholar] [CrossRef] [Green Version]
  79. Cui, Z.F.; Zhang, J.L.; Binosi, D.; de Soto, F.; Mezrag, C.; Papavassiliou, J.; Roberts, C.D.; Rodríguez-Quintero, J.; Segovia, J.; Zafeiropoulos, S. Effective charge from lattice QCD. Chin. Phys. C 2020, 44, 083102. [Google Scholar] [CrossRef]
  80. Zafeiropoulos, S.; Boucaud, P.; De Soto, F.; Rodríguez-Quintero, J.; Segovia, J. Strong Running Coupling from the Gauge Sector of Domain Wall Lattice QCD with Physical Quark Masses. Phys. Rev. Lett. 2019, 122, 162002. [Google Scholar] [CrossRef] [Green Version]
  81. Aguilar, A.C.; De Soto, F.; Ferreira, M.N.; Papavassiliou, J.; Rodríguez-Quintero, J.; Zafeiropoulos, S. Gluon propagator and three-gluon vertex with dynamical quarks. Eur. Phys. J. 2020, C80, 154. [Google Scholar] [CrossRef] [Green Version]
  82. Maas, A.; Vujinović, M. More on the three-gluon vertex in SU(2) Yang-Mills theory in three and four dimensions. SciPost Phys. Core 2022, 5, 019. [Google Scholar] [CrossRef]
  83. Kızılersü, A.; Oliveira, O.; Silva, P.J.; Skullerud, J.I.; Sternbeck, A. Quark-gluon vertex from Nf = 2 lattice QCD. Phys. Rev. D 2021, 103, 114515. [Google Scholar] [CrossRef]
  84. Aguilar, A.C.; De Soto, F.; Ferreira, M.N.; Papavassiliou, J.; Rodríguez-Quintero, J. Infrared facets of the three-gluon vertex. Phys. Lett. B 2021, 818, 136352. [Google Scholar] [CrossRef]
  85. Aguilar, A.C.; Ambrósio, C.O.; De Soto, F.; Ferreira, M.N.; Oliveira, B.M.; Papavassiliou, J.; Rodríguez-Quintero, J. Ghost dynamics in the soft gluon limit. Phys. Rev. D 2021, 104, 054028. [Google Scholar] [CrossRef]
  86. Pinto-Gómez, F.; De Soto, F.; Ferreira, M.N.; Papavassiliou, J.; Rodríguez-Quintero, J. Lattice three-gluon vertex in extended kinematics: Planar degeneracy. Phys. Lett. B 2023, 838, 137737. [Google Scholar] [CrossRef]
  87. Pinto-Gomez, F.; de Soto, F. Three-gluon vertex in Landau-gauge from quenched-lattice QCD in general kinematics. In Proceedings of the 15th Conference on Quark Confinement and the Hadron Spectrum, Stavanger, Norway, 1–6 August 2022. [Google Scholar]
  88. Pinto-Gómez, F.; de Soto, F.; Ferreira, M.N.; Papavassiliou, J.; Rodríguez-Quintero, J. General kinematics of the three-gluon vertex from quenched lattice QCD. arXiv 2022, arXiv:2212.11894. [Google Scholar]
  89. Roberts, C.D.; Schmidt, S.M. Reflections upon the emergence of hadronic mass. Eur. Phys. J. ST 2020, 229, 3319–3340. [Google Scholar] [CrossRef]
  90. Roberts, C.D. On Mass and Matter. AAPPS Bull. 2021, 31, 6. [Google Scholar] [CrossRef]
  91. Roberts, C.D.; Richards, D.G.; Horn, T.; Chang, L. Insights into the emergence of mass from studies of pion and kaon structure. Prog. Part. Nucl. Phys. 2021, 120, 103883. [Google Scholar] [CrossRef]
  92. Binosi, D. Emergent Hadron Mass in Strong Dynamics. Few Body Syst. 2022, 63, 42. [Google Scholar] [CrossRef]
  93. Papavassiliou, J. Emergence of mass in the gauge sector of QCD*. Chin. Phys. C 2022, 46, 112001. [Google Scholar] [CrossRef]
  94. Cornwall, J.M. Quark Confinement and Vortices in Massive Gauge Invariant QCD. Nucl. Phys. 1979, B157, 392. [Google Scholar] [CrossRef]
  95. Parisi, G.; Petronzio, R. On Low-Energy Tests of QCD. Phys. Lett. 1980, B94, 51. [Google Scholar] [CrossRef] [Green Version]
  96. Cornwall, J.M. Dynamical Mass Generation in Continuum QCD. Phys. Rev. D 1982, 26, 1453. [Google Scholar] [CrossRef]
  97. Bernard, C.W. Monte Carlo Evaluation of the Effective Gluon Mass. Phys. Lett. B 1982, 108, 431–434. [Google Scholar] [CrossRef]
  98. Bernard, C.W. Adjoint Wilson Lines and the Effective Gluon Mass. Nucl. Phys. B 1983, 219, 341–357. [Google Scholar] [CrossRef]
  99. Donoghue, J.F. The Gluon ’Mass’ in the Bag Model. Phys. Rev. D 1984, 29, 2559. [Google Scholar] [CrossRef]
  100. Cornwall, J.M.; Papavassiliou, J. Gauge Invariant Three Gluon Vertex in QCD. Phys. Rev. D 1989, 40, 3474. [Google Scholar] [CrossRef] [Green Version]
  101. Lavelle, M. Gauge invariant effective gluon mass from the operator product expansion. Phys. Rev. D 1991, 44, 26–28. [Google Scholar] [CrossRef]
  102. Halzen, F.; Krein, G.I.; Natale, A.A. Relating the QCD pomeron to an effective gluon mass. Phys. Rev. 1993, D47, 295–298. [Google Scholar] [CrossRef] [Green Version]
  103. Wilson, K.G.; Walhout, T.S.; Harindranath, A.; Zhang, W.M.; Perry, R.J.; Glazek, S.D. Nonperturbative QCD: A Weak coupling treatment on the light front. Phys. Rev. 1994, D49, 6720–6766. [Google Scholar] [CrossRef] [Green Version]
  104. Mihara, A.; Natale, A.A. Dynamical gluon mass corrections in heavy quarkonia decays. Phys. Lett. 2000, B482, 378–382. [Google Scholar] [CrossRef] [Green Version]
  105. Philipsen, O. On the nonperturbative gluon mass and heavy quark physics. Nucl. Phys. 2002, B628, 167–192. [Google Scholar] [CrossRef] [Green Version]
  106. Kondo, K.I. Vacuum condensate of mass dimension 2 as the origin of mass gap and quark confinement. Phys. Lett. 2001, B514, 335–345. [Google Scholar] [CrossRef]
  107. Aguilar, A.C.; Natale, A.A.; Rodrigues da Silva, P.S. Relating a gluon mass scale to an infrared fixed point in pure gauge QCD. Phys. Rev. Lett. 2003, 90, 152001. [Google Scholar] [CrossRef] [Green Version]
  108. Aguilar, A.C.; Natale, A.A. A Dynamical gluon mass solution in a coupled system of the Schwinger-Dyson equations. J. High Energy Phys. 2004, 8, 057. [Google Scholar] [CrossRef] [Green Version]
  109. Aguilar, A.C.; Papavassiliou, J. Gluon mass generation in the PT-BFM scheme. J. High Energy Phys. 2006, 12, 012. [Google Scholar] [CrossRef] [Green Version]
  110. Epple, D.; Reinhardt, H.; Schleifenbaum, W.; Szczepaniak, A.P. Subcritical solution of the Yang-Mills Schroedinger equation in the Coulomb gauge. Phys. Rev. 2008, D77, 085007. [Google Scholar] [CrossRef]
  111. Aguilar, A.C.; Papavassiliou, J. On dynamical gluon mass generation. Eur. Phys. J. 2007, A31, 742–745. [Google Scholar] [CrossRef]
  112. Aguilar, A.C.; Binosi, D.; Papavassiliou, J. Gluon and ghost propagators in the Landau gauge: Deriving lattice results from Schwinger-Dyson equations. Phys. Rev. 2008, D78, 025010. [Google Scholar] [CrossRef] [Green Version]
  113. Aguilar, A.C.; Papavassiliou, J. Gluon mass generation without seagull divergences. Phys. Rev. 2010, D81, 034003. [Google Scholar] [CrossRef] [Green Version]
  114. Campagnari, D.R.; Reinhardt, H. Non-Gaussian wave functionals in Coulomb gauge Yang–Mills theory. Phys. Rev. 2010, D82, 105021. [Google Scholar] [CrossRef] [Green Version]
  115. Fagundes, D.A.; Luna, E.G.S.; Menon, M.J.; Natale, A.A. Aspects of a Dynamical Gluon Mass Approach to Elastic Hadron Scattering at LHC. Nucl. Phys. A 2012, 886, 48–70. [Google Scholar] [CrossRef] [Green Version]
  116. Aguilar, A.C.; Binosi, D.; Papavassiliou, J. The dynamical equation of the effective gluon mass. Phys. Rev. 2011, D84, 085026. [Google Scholar] [CrossRef]
  117. Aguilar, A.C.; Ibanez, D.; Mathieu, V.; Papavassiliou, J. Massless bound-state excitations and the Schwinger mechanism in QCD. Phys. Rev. 2012, D85, 014018. [Google Scholar] [CrossRef] [Green Version]
  118. Aguilar, A.C.; Binosi, D.; Papavassiliou, J. Gluon mass through ghost synergy. J. High Energy Phys. 2012, 01, 050. [Google Scholar] [CrossRef] [Green Version]
  119. Aguilar, A.C.; Binosi, D.; Papavassiliou, J. Gluon mass generation in the presence of dynamical quarks. Phys. Rev. 2013, D88, 074010. [Google Scholar] [CrossRef] [Green Version]
  120. Glazek, S.D.; Gómez-Rocha, M.; More, J.; Serafin, K. Renormalized quark–antiquark Hamiltonian induced by a gluon mass ansatz in heavy-flavor QCD. Phys. Lett. 2017, B773, 172–178. [Google Scholar] [CrossRef]
  121. Binosi, D.; Papavassiliou, J. Coupled dynamics in gluon mass generation and the impact of the three-gluon vertex. Phys. Rev. 2018, D97, 054029. [Google Scholar] [CrossRef] [Green Version]
  122. Aguilar, A.C.; Ferreira, M.N.; Figueiredo, C.T.; Papavassiliou, J. Gluon mass scale through nonlinearities and vertex interplay. Phys. Rev. D 2019, 100, 094039. [Google Scholar] [CrossRef] [Green Version]
  123. Eichmann, G.; Pawlowski, J.M.; Silva, J.M. Mass generation in Landau-gauge Yang-Mills theory. Phys. Rev. D 2021, 104, 114016. [Google Scholar] [CrossRef]
  124. Aguilar, A.C.; Ferreira, M.N.; Papavassiliou, J. Exploring smoking-gun signals of the Schwinger mechanism in QCD. Phys. Rev. D 2022, 105, 014030. [Google Scholar] [CrossRef]
  125. Horak, J.; Ihssen, F.; Papavassiliou, J.; Pawlowski, J.M.; Weber, A.; Wetterich, C. Gluon condensates and effective gluon mass. SciPost Phys. 2022, 13, 042. [Google Scholar] [CrossRef]
  126. Aguilar, A.C.; De Soto, F.; Ferreira, M.N.; Papavassiliou, J.; Pinto-Gómez, F.; Roberts, C.D.; Rodríguez-Quintero, J. Schwinger mechanism for gluons from lattice QCD. arXiv 2022, arXiv:2211.12594. [Google Scholar]
  127. Schwinger, J.S. Gauge Invariance and Mass. Phys. Rev. 1962, 125, 397–398. [Google Scholar] [CrossRef]
  128. Schwinger, J.S. Gauge Invariance and Mass. 2. Phys. Rev. 1962, 128, 2425–2429. [Google Scholar] [CrossRef]
  129. Watson, N.J. The gauge-independent QCD effective charge. Nucl. Phys. 1997, B494, 388–432. [Google Scholar] [CrossRef] [Green Version]
  130. Binosi, D.; Papavassiliou, J. The QCD effective charge to all orders. Nucl. Phys. Proc. Suppl. 2003, 121, 281–284. [Google Scholar] [CrossRef] [Green Version]
  131. Aguilar, A.C.; Binosi, D.; Papavassiliou, J.; Rodriguez-Quintero, J. Non-perturbative comparison of QCD effective charges. Phys. Rev. 2009, D80, 085018. [Google Scholar] [CrossRef] [Green Version]
  132. Gell-Mann, M.; Low, F.E. Quantum electrodynamics at small distances. Phys. Rev. 1954, 95, 1300–1312. [Google Scholar] [CrossRef] [Green Version]
  133. Itzykson, C.; Zuber, J.B. Quantum Field Theory; International Series in Pure and Applied Physics; Mcgraw-Hill: New York, NY, USA, 1980; 705p. [Google Scholar]
  134. Nambu, Y.; Jona-Lasinio, G. Dynamical model of elementary particles based on an analogy with superconductivity. I. Phys. Rev. 1961, 122, 345–358. [Google Scholar] [CrossRef] [Green Version]
  135. Lane, K.D. Asymptotic Freedom and Goldstone Realization of Chiral Symmetry. Phys. Rev. 1974, D10, 2605. [Google Scholar] [CrossRef]
  136. Politzer, H.D. Effective Quark Masses in the Chiral Limit. Nucl. Phys. 1976, B117, 397. [Google Scholar] [CrossRef]
  137. Miransky, V.A.; Fomin, P.I. Chiral symmetry breakdown and the spectrum of pseudoscalar mesons in quantum chromodynamics. Phys. Lett. 1981, B105, 387–391. [Google Scholar] [CrossRef]
  138. Atkinson, D.; Johnson, P.W. Chiral Symmetry Breaking in QCD. 2. Running Coupling Constant. Phys. Rev. 1988, D37, 2296–2299. [Google Scholar] [CrossRef]
  139. Brown, N.; Pennington, M.R. Studies of confinement: How quarks and gluons propagate. Phys. Rev. 1988, D38, 2266. [Google Scholar] [CrossRef]
  140. Williams, A.G.; Krein, G.; Roberts, C.D. Quark propagator in an Ansatz approach to QCD. Ann. Phys. 1991, 210, 464–485. [Google Scholar] [CrossRef] [Green Version]
  141. Papavassiliou, J.; Cornwall, J.M. Coupled fermion gap and vertex equations for chiral symmetry breakdown in QCD. Phys. Rev. 1991, D44, 1285–1297. [Google Scholar] [CrossRef]
  142. Hawes, F.T.; Roberts, C.D.; Williams, A.G. Dynamical chiral symmetry breaking and confinement with an infrared vanishing gluon propagator. Phys. Rev. 1994, D49, 4683–4693. [Google Scholar] [CrossRef] [Green Version]
  143. Natale, A.A.; Rodrigues da Silva, P.S. Critical coupling for dynamical chiral-symmetry breaking with an infrared finite gluon propagator. Phys. Lett. 1997, B392, 444–451. [Google Scholar] [CrossRef] [Green Version]
  144. Fischer, C.S.; Alkofer, R. Nonperturbative propagators, running coupling and dynamical quark mass of Landau gauge QCD. Phys. Rev. 2003, D67, 094020. [Google Scholar] [CrossRef] [Green Version]
  145. Maris, P.; Roberts, C.D. Dyson-Schwinger equations: A Tool for hadron physics. Int. J. Mod. Phys. 2003, E12, 297–365. [Google Scholar] [CrossRef] [Green Version]
  146. Aguilar, A.C.; Nesterenko, A.; Papavassiliou, J. Infrared enhanced analytic coupling and chiral symmetry breaking in QCD. J. Phys. 2005, G31, 997. [Google Scholar] [CrossRef]
  147. Bowman, P.O.; Heller, U.M.; Leinweber, D.B.; Parappilly, M.B.; Williams, A.G.; Zhang, J.b. Unquenched quark propagator in Landau gauge. Phys. Rev. 2005, D71, 054507. [Google Scholar] [CrossRef] [Green Version]
  148. Sauli, V.; Adam, J., Jr.; Bicudo, P. Dynamical chiral symmetry breaking with integral Minkowski representations. Phys. Rev. 2007, D75, 087701. [Google Scholar] [CrossRef]
  149. Cornwall, J.M. Center vortices, the functional Schrodinger equation, and CSB. In Proceedings of the 419th WE-Heraeus-Seminar: Approaches to Quantum Chromodynamics, Oberwoelz, Austria, 7–13 September 2022. [Google Scholar]
  150. Alkofer, R.; Fischer, C.S.; Llanes-Estrada, F.J.; Schwenzer, K. The Quark-gluon vertex in Landau gauge QCD: Its role in dynamical chiral symmetry breaking and quark confinement. Ann. Phys. 2009, 324, 106–172. [Google Scholar] [CrossRef] [Green Version]
  151. Aguilar, A.C.; Papavassiliou, J. Chiral symmetry breaking with lattice propagators. Phys. Rev. 2011, D83, 014013. [Google Scholar] [CrossRef] [Green Version]
  152. Rojas, E.; de Melo, J.; El-Bennich, B.; Oliveira, O.; Frederico, T. On the Quark-Gluon Vertex and Quark-Ghost Kernel: Combining Lattice Simulations with Dyson-Schwinger equations. J. High Energy Phys. 2013, 10, 193. [Google Scholar] [CrossRef] [Green Version]
  153. Mitter, M.; Pawlowski, J.M.; Strodthoff, N. Chiral symmetry breaking in continuum QCD. Phys. Rev. 2015, D91, 054035. [Google Scholar] [CrossRef] [Green Version]
  154. Braun, J.; Fister, L.; Pawlowski, J.M.; Rennecke, F. From Quarks and Gluons to Hadrons: Chiral Symmetry Breaking in Dynamical QCD. Phys. Rev. 2016, D94, 034016. [Google Scholar] [CrossRef] [Green Version]
  155. Heupel, W.; Goecke, T.; Fischer, C.S. Beyond Rainbow-Ladder in bound state equations. Eur. Phys. J. 2014, A50, 85. [Google Scholar] [CrossRef] [Green Version]
  156. Binosi, D.; Chang, L.; Papavassiliou, J.; Qin, S.X.; Roberts, C.D. Natural constraints on the gluon-quark vertex. Phys. Rev. 2017, D95, 031501. [Google Scholar] [CrossRef] [Green Version]
  157. Aguilar, A.C.; Cardona, J.C.; Ferreira, M.N.; Papavassiliou, J. Quark gap equation with non-abelian Ball-Chiu vertex. Phys. Rev. 2018, D98, 014002. [Google Scholar] [CrossRef] [Green Version]
  158. Gao, F.; Papavassiliou, J.; Pawlowski, J.M. Fully coupled functional equations for the quark sector of QCD. Phys. Rev. D 2021, 103, 094013. [Google Scholar] [CrossRef]
  159. Eichten, E.; Feinberg, F. Dynamical Symmetry Breaking of Nonabelian Gauge Symmetries. Phys. Rev. D 1974, 10, 3254–3279. [Google Scholar] [CrossRef]
  160. Smit, J. On the Possibility That Massless Yang-Mills Fields Generate Massive Vector Particles. Phys. Rev. D 1974, 10, 2473. [Google Scholar] [CrossRef]
  161. Binosi, D.; Iba nez, D.; Papavassiliou, J. The all-order equation of the effective gluon mass. Phys. Rev. 2012, D86, 085033. [Google Scholar] [CrossRef] [Green Version]
  162. Tissier, M.; Wschebor, N. Infrared propagators of Yang-Mills theory from perturbation theory. Phys. Rev. D 2010, 82, 101701. [Google Scholar] [CrossRef] [Green Version]
  163. Serreau, J.; Tissier, M. Lifting the Gribov ambiguity in Yang-Mills theories. Phys. Lett. 2012, B712, 97–103. [Google Scholar] [CrossRef]
  164. Peláez, M.; Tissier, M.; Wschebor, N. Two-point correlation functions of QCD in the Landau gauge. Phys. Rev. D 2014, 90, 065031. [Google Scholar] [CrossRef] [Green Version]
  165. Siringo, F. Analytical study of Yang–Mills theory in the infrared from first principles. Nucl. Phys. 2016, B907, 572–596. [Google Scholar] [CrossRef] [Green Version]
  166. Aguilar, A.C.; Binosi, D.; Figueiredo, C.T.; Papavassiliou, J. Unified description of seagull cancellations and infrared finiteness of gluon propagators. Phys. Rev. 2016, D94, 045002. [Google Scholar] [CrossRef] [Green Version]
  167. Osterwalder, K.; Schrader, R. Axioms for Euclidean Green’S Functions. Commun. Math. Phys. 1973, 31, 83–112. [Google Scholar] [CrossRef]
  168. Osterwalder, K.; Schrader, R. Axioms for Euclidean Green’s Functions. 2. Commun. Math. Phys. 1975, 42, 281. [Google Scholar] [CrossRef]
  169. Glimm, J.; Jaffe, A.M. Quantum Physics. A Functional Integral Point of View; Springer: New York, NY, USA, 1981. [Google Scholar]
  170. Krein, G.; Roberts, C.D.; Williams, A.G. On the implications of confinement. Int. J. Mod. Phys. 1992, A7, 5607–5624. [Google Scholar] [CrossRef]
  171. Cornwall, J.M. Positivity violations in QCD. Mod. Phys. Lett. 2013, A28, 1330035. [Google Scholar] [CrossRef] [Green Version]
  172. Aguilar, A.C.; Binosi, D.; Iba nez, D.; Papavassiliou, J. Effects of divergent ghost loops on Green’s functions of QCD. Phys. Rev. 2014, D89, 085008. [Google Scholar] [CrossRef] [Green Version]
  173. Pelaez, M.; Tissier, M.; Wschebor, N. Three-point correlation functions in Yang-Mills theory. Phys. Rev. 2013, D88, 125003. [Google Scholar] [CrossRef] [Green Version]
  174. Blum, A.; Huber, M.Q.; Mitter, M.; von Smekal, L. Gluonic three-point correlations in pure Landau gauge QCD. Phys. Rev. 2014, D89, 061703. [Google Scholar] [CrossRef] [Green Version]
  175. Eichmann, G.; Williams, R.; Alkofer, R.; Vujinovic, M. The three-gluon vertex in Landau gauge. Phys. Rev. 2014, D89, 105014. [Google Scholar] [CrossRef] [Green Version]
  176. Williams, R.; Fischer, C.S.; Heupel, W. Light mesons in QCD and unquenching effects from the 3PI effective action. Phys. Rev. 2016, D93, 034026. [Google Scholar] [CrossRef] [Green Version]
  177. Blum, A.L.; Alkofer, R.; Huber, M.Q.; Windisch, A. Unquenching the three-gluon vertex: A status report. Acta Phys. Pol. Suppl. 2015, 8, 321. [Google Scholar] [CrossRef] [Green Version]
  178. Cyrol, A.K.; Fister, L.; Mitter, M.; Pawlowski, J.M.; Strodthoff, N. Landau gauge Yang-Mills correlation functions. Phys. Rev. 2016, D94, 054005. [Google Scholar] [CrossRef] [Green Version]
  179. Aguilar, A.C.; Ferreira, M.N.; Figueiredo, C.T.; Papavassiliou, J. Nonperturbative Ball-Chiu construction of the three-gluon vertex. Phys. Rev. 2019, D99, 094010. [Google Scholar] [CrossRef] [Green Version]
  180. Barrios, N.; Peláez, M.; Reinosa, U. Two-loop three-gluon vertex from the Curci-Ferrari model and its leading infrared behavior to all loop orders. Phys. Rev. D 2022, 106, 114039. [Google Scholar] [CrossRef]
  181. Rivers, R.J. Path Integral Methods in Quantum Field Theory; Cambridge Monographs on Mathematical Physics. Cambridge University Press: Cambridge, UK, 1988. [Google Scholar] [CrossRef]
  182. Fujikawa, K.; Lee, B.W.; Sanda, A.I. Generalized Renormalizable Gauge Formulation of Spontaneously Broken Gauge Theories. Phys. Rev. D 1972, 6, 2923–2943. [Google Scholar] [CrossRef]
  183. Aguilar, A.C.; Papavassiliou, J. Infrared finite ghost propagator in the Feynman gauge. Phys. Rev. 2008, D77, 125022. [Google Scholar] [CrossRef] [Green Version]
  184. Huber, M.Q.; Schwenzer, K.; Alkofer, R. On the infrared scaling solution of SU(N) Yang-Mills theories in the maximally Abelian gauge. Eur. Phys. J. C 2010, 68, 581–600. [Google Scholar] [CrossRef]
  185. Siringo, F. Gluon propagator in Feynman gauge by the method of stationary variance. Phys. Rev. D 2014, 90, 094021. [Google Scholar] [CrossRef] [Green Version]
  186. Aguilar, A.C.; Binosi, D.; Papavassiliou, J. Yang-Mills two-point functions in linear covariant gauges. Phys. Rev. 2015, D91, 085014. [Google Scholar] [CrossRef] [Green Version]
  187. Huber, M.Q. Gluon and ghost propagators in linear covariant gauges. Phys. Rev. 2015, D91, 085018. [Google Scholar] [CrossRef] [Green Version]
  188. Capri, M.A.L.; Fiorentini, D.; Guimaraes, M.S.; Mintz, B.W.; Palhares, L.F.; Sorella, S.P.; Dudal, D.; Justo, I.F.; Pereira, A.D.; Sobreiro, R.F. Exact nilpotent nonperturbative BRST symmetry for the Gribov-Zwanziger action in the linear covariant gauge. Phys. Rev. D 2015, 92, 045039. [Google Scholar] [CrossRef] [Green Version]
  189. Aguilar, A.C.; Binosi, D.; Papavassiliou, J. Schwinger mechanism in linear covariant gauges. Phys. Rev. 2017, D95, 034017. [Google Scholar] [CrossRef] [Green Version]
  190. De Meerleer, T.; Dudal, D.; Sorella, S.P.; Dall’Olio, P.; Bashir, A. Landau-Khalatnikov-Fradkin Transformations, Nielsen Identities, Their Equivalence and Implications for QCD. Phys. Rev. D 2020, 101, 085005. [Google Scholar] [CrossRef]
  191. Napetschnig, M.; Alkofer, R.; Huber, M.Q.; Pawlowski, J.M. Yang-Mills propagators in linear covariant gauges from Nielsen identities. Phys. Rev. D 2021, 104, 054003. [Google Scholar] [CrossRef]
  192. Binosi, D.; Papavassiliou, J. Gauge-invariant truncation scheme for the Schwinger-Dyson equations of QCD. Phys. Rev. 2008, D77, 061702. [Google Scholar] [CrossRef] [Green Version]
  193. Pilaftsis, A. Generalized pinch technique and the background field method in general gauges. Nucl. Phys. B 1997, 487, 467–491. [Google Scholar] [CrossRef] [Green Version]
  194. Binosi, D.; Papavassiliou, J. The Pinch technique to all orders. Phys. Rev. D 2002, 66, 111901. [Google Scholar] [CrossRef] [Green Version]
  195. Binosi, D.; Papavassiliou, J. Pinch technique selfenergies and vertices to all orders in perturbation theory. J. Phys. G 2004, G30, 203. [Google Scholar] [CrossRef] [Green Version]
  196. DeWitt, B.S. Quantum Theory of Gravity. 2. The Manifestly Covariant Theory. Phys. Rev. 1967, 162, 1195–1239. [Google Scholar] [CrossRef]
  197. ’t Hooft, G. Renormalizable Lagrangians for Massive Yang-Mills Fields. Nucl. Phys. B 1971, 35, 167–188. [Google Scholar] [CrossRef] [Green Version]
  198. Honerkamp, J. The Question of invariant renormalizability of the massless Yang-Mills theory in a manifest covariant approach. Nucl. Phys. B 1972, 48, 269–287. [Google Scholar] [CrossRef] [Green Version]
  199. Kallosh, R.E. The Renormalization in Nonabelian Gauge Theories. Nucl. Phys. B 1974, 78, 293–312. [Google Scholar] [CrossRef]
  200. Kluberg-Stern, H.; Zuber, J.B. Renormalization of Nonabelian Gauge Theories in a Background Field Gauge. 1. Green Functions. Phys. Rev. D 1975, 12, 482–488. [Google Scholar] [CrossRef]
  201. Arefeva, I.Y.; Faddeev, L.D.; Slavnov, A.A. Generating Functional for the s Matrix in Gauge Theories. Teor. Mat. Fiz. 1974, 21, 311–321. [Google Scholar] [CrossRef]
  202. Abbott, L. The Background Field Method Beyond One Loop. Nucl. Phys. B 1981, 185, 189–203. [Google Scholar] [CrossRef] [Green Version]
  203. Weinberg, S. Effective Gauge Theories. Phys. Lett. B 1980, 91, 51–55. [Google Scholar] [CrossRef]
  204. Abbott, L.F. Introduction to the Background Field Method. Acta Phys. Polon. 1982, B13, 33. [Google Scholar]
  205. Shore, G.M. Symmetry Restoration and the Background Field Method in Gauge Theories. Ann. Phys. 1981, 137, 262. [Google Scholar] [CrossRef]
  206. Abbott, L.F.; Grisaru, M.T.; Schaefer, R.K. The Background Field Method and the S Matrix. Nucl. Phys. B 1983, 229, 372–380. [Google Scholar] [CrossRef]
  207. Taylor, J. Ward Identities and Charge Renormalization of the Yang-Mills Field. Nucl. Phys. B 1971, 33, 436–444. [Google Scholar] [CrossRef]
  208. Slavnov, A. Ward Identities in Gauge Theories. Theor. Math. Phys. 1972, 10, 99–107. [Google Scholar] [CrossRef]
  209. Grassi, P.A.; Hurth, T.; Steinhauser, M. Practical algebraic renormalization. Ann. Phys. 2001, 288, 197–248. [Google Scholar] [CrossRef] [Green Version]
  210. Grassi, P.A.; Hurth, T.; Steinhauser, M. The Algebraic method. Nucl. Phys. B 2001, 610, 215–250. [Google Scholar] [CrossRef] [Green Version]
  211. Binosi, D.; Papavassiliou, J. Pinch technique and the Batalin-Vilkovisky formalism. Phys. Rev. 2002, D66, 025024. [Google Scholar] [CrossRef]
  212. Aguilar, A.C.; Binosi, D.; Papavassiliou, J. Indirect determination of the Kugo-Ojima function from lattice data. J. High Energy Phys. 2009, 11, 066. [Google Scholar] [CrossRef] [Green Version]
  213. Binosi, D.; Quadri, A. AntiBRST symmetry and Background Field Method. Phys. Rev. 2013, D88, 085036. [Google Scholar] [CrossRef] [Green Version]
  214. Iba nez, D.; Papavassiliou, J. Gluon mass generation in the massless bound-state formalism. Phys. Rev. 2013, D87, 034008. [Google Scholar] [CrossRef] [Green Version]
  215. Aguilar, A.C.; Binosi, D.; Figueiredo, C.T.; Papavassiliou, J. Evidence of ghost suppression in gluon mass scale dynamics. Eur. Phys. J. 2018, C78, 181. [Google Scholar] [CrossRef] [Green Version]
  216. Munczek, H. Dynamical chiral symmetry breaking, Goldstone’s theorem and the consistency of the Schwinger-Dyson and Bethe–Salpeter Equations. Phys. Rev. 1995, D52, 4736–4740. [Google Scholar] [CrossRef] [Green Version]
  217. Bender, A.; Roberts, C.D.; Von Smekal, L. Goldstone theorem and diquark confinement beyond rainbow ladder approximation. Phys. Lett. 1996, B380, 7–12. [Google Scholar] [CrossRef] [Green Version]
  218. Maris, P.; Roberts, C.D.; Tandy, P.C. Pion mass and decay constant. Phys. Lett. 1998, B420, 267–273. [Google Scholar] [CrossRef] [Green Version]
  219. Maris, P.; Roberts, C.D. Pi- and K meson Bethe–Salpeter amplitudes. Phys. Rev. C 1997, 56, 3369–3383. [Google Scholar] [CrossRef] [Green Version]
  220. Chang, L.; Roberts, C.D. Sketching the Bethe–Salpeter kernel. Phys. Rev. Lett. 2009, 103, 081601. [Google Scholar] [CrossRef] [Green Version]
  221. Chang, L.; Roberts, C.D.; Tandy, P.C. Selected highlights from the study of mesons. Chin. J. Phys. 2011, 49, 955–1004. [Google Scholar]
  222. Qin, S.X.; Roberts, C.D. Resolving the Bethe–Salpeter Kernel. Chin. Phys. Lett. 2021, 38, 071201. [Google Scholar] [CrossRef]
  223. Huber, M.Q. Correlation functions of three-dimensional Yang-Mills theory from Dyson-Schwinger equations. Phys. Rev. D 2016, 93, 085033. [Google Scholar] [CrossRef] [Green Version]
  224. Schleifenbaum, W.; Maas, A.; Wambach, J.; Alkofer, R. Infrared behaviour of the ghost–gluon vertex in Landau gauge Yang-Mills theory. Phys. Rev. D 2005, 72, 014017. [Google Scholar] [CrossRef] [Green Version]
  225. Boucaud, P.; Leroy, J.; A., L.Y.; Micheli, J.; Pène, O.; Rodríguez-Quintero, J. On the IR behaviour of the Landau-gauge ghost propagator. J. High Energy Phys. 2008, 06, 099. [Google Scholar] [CrossRef] [Green Version]
  226. Huber, M.Q.; von Smekal, L. On the influence of three-point functions on the propagators of Landau gauge Yang-Mills theory. J. High Energy Phys. 2013, 04, 149. [Google Scholar] [CrossRef] [Green Version]
  227. Aguilar, A.C.; Iba nez, D.; Papavassiliou, J. Ghost propagator and ghost–gluon vertex from Schwinger-Dyson equations. Phys. Rev. 2013, D87, 114020. [Google Scholar] [CrossRef] [Green Version]
  228. Aguilar, A.C.; Ferreira, M.N.; Figueiredo, C.T.; Papavassiliou, J. Nonperturbative structure of the ghost–gluon kernel. Phys. Rev. 2019, D99, 034026. [Google Scholar] [CrossRef] [Green Version]
  229. Dudal, D.; Gracey, J.A.; Sorella, S.P.; Vandersickel, N.; Verschelde, H. A refinement of the Gribov-Zwanziger approach in the Landau gauge: Infrared propagators in harmony with the lattice results. Phys. Rev. 2008, D78, 065047. [Google Scholar] [CrossRef] [Green Version]
  230. Boucaud, P.; Leroy, J.P.; Le Yaouanc, A.; Micheli, J.; Pene, O.; Rodríguez-Quintero, J. IR finiteness of the ghost dressing function from numerical resolution of the ghost SD equation. J. High Energy Phys. 2008, 06, 012. [Google Scholar] [CrossRef]
  231. Kondo, K.I. Infrared behavior of the ghost propagator in the Landau gauge Yang-Mills theory. Prog. Theor. Phys. 2010, 122, 1455–1475. [Google Scholar] [CrossRef] [Green Version]
  232. Pennington, M.R.; Wilson, D.J. Are the Dressed Gluon and Ghost Propagators in the Landau Gauge presently determined in the confinement regime of QCD? Phys. Rev. D 2011, 84, 094028, Erratum in Phys. Rev. D 2011, 84, 119901. [Google Scholar] [CrossRef]
  233. Dudal, D.; Oliveira, O.; Rodriguez-Quintero, J. Nontrivial ghost–gluon vertex and the match of RGZ, DSE and lattice Yang-Mills propagators. Phys. Rev. 2012, D86, 105005, Erratum in Phys. Rev. D 2012, 86, 109902. [Google Scholar] [CrossRef] [Green Version]
  234. Papavassiliou, J.; Aguilar, A.C.; Ferreira, M.N. Theory and phenomenology of the three-gluon vertex. Rev. Mex. Fis. Suppl. 2022, 3, 0308112. [Google Scholar] [CrossRef]
  235. Meyers, J.; Swanson, E.S. Spin Zero Glueballs in the Bethe–Salpeter Formalism. Phys. Rev. 2013, D87, 036009. [Google Scholar] [CrossRef] [Green Version]
  236. Sanchis-Alepuz, H.; Fischer, C.S.; Kellermann, C.; von Smekal, L. Glueballs from the Bethe–Salpeter equation. Phys. Rev. 2015, D92, 034001. [Google Scholar] [CrossRef] [Green Version]
  237. Xu, S.S.; Cui, Z.F.; Chang, L.; Papavassiliou, J.; Roberts, C.D.; Zong, H.S. New perspective on hybrid mesons. Eur. Phys. J. 2019, A55, 113. [Google Scholar] [CrossRef] [Green Version]
  238. Souza, E.V.; Ferreira, M.N.; Aguilar, A.C.; Papavassiliou, J.; Roberts, C.D.; Xu, S.S. Pseudoscalar glueball mass: A window on three-gluon interactions. Eur. Phys. J. A 2020, 56, 25. [Google Scholar] [CrossRef] [Green Version]
  239. Huber, M.Q.; Fischer, C.S.; Sanchis-Alepuz, H. Spectrum of scalar and pseudoscalar glueballs from functional methods. Eur. Phys. J. C 2020, 80, 1077. [Google Scholar] [CrossRef] [PubMed]
  240. Huber, M.Q.; Fischer, C.S.; Sanchis-Alepuz, H. Higher spin glueballs from functional methods. Eur. Phys. J. C 2021, 81, 1083, Erratum in Eur. Phys. J. C 2022, 82, 38. [Google Scholar] [CrossRef]
  241. Aguilar, A.C.; Ferreira, M.N.; Papavassiliou, J. Novel sum rules for the three-point sector of QCD. Eur. Phys. J. C 2020, 80, 887. [Google Scholar] [CrossRef]
  242. Aguilar, A.C.; Ferreira, M.N.; Papavassiliou, J. Gluon dynamics from an ordinary differential equation. Eur. Phys. J. C 2021, 81, 54. [Google Scholar] [CrossRef]
  243. Cornwall, J.M.; Papavassiliou, J.; Binosi, D. The Pinch Technique and its Applications to Non-Abelian Gauge Theories; Cambridge University Press: Cambridge, UK, 2010; Volume 31. [Google Scholar]
  244. Binosi, D.; Papavassiliou, J. New Schwinger-Dyson equations for non-Abelian gauge theories. J. High Energy Phys. 2008, 11, 063. [Google Scholar] [CrossRef] [Green Version]
  245. Grassi, P.A.; Hurth, T.; Quadri, A. On the Landau background gauge fixing and the IR properties of YM Green functions. Phys. Rev. 2004, D70, 105014. [Google Scholar]
  246. Kugo, T. The Universal renormalization factors Z(1) / Z(3) and color confinement condition in nonAbelian gauge theory. arXiv 1995, arXiv:9511033. [Google Scholar]
  247. Kondo, K.I. Kugo-Ojima color confinement criterion and Gribov-Zwanziger horizon condition. Phys. Lett. B 2009, 678, 322–330. [Google Scholar] [CrossRef]
  248. Becchi, C.; Rouet, A.; Stora, R. Renormalization of the Abelian Higgs-Kibble Model. Commun. Math. Phys. 1975, 42, 127–162. [Google Scholar] [CrossRef] [Green Version]
  249. Becchi, C.; Rouet, A.; Stora, R. Renormalization of Gauge Theories. Ann. Phys. 1976, 98, 287–321. [Google Scholar] [CrossRef] [Green Version]
  250. Tyutin, I.V. Gauge Invariance in Field Theory and Statistical Physics in Operator Formalism. arXiv 1975, arXiv:0812.0580. [Google Scholar]
  251. Ball, J.S.; Chiu, T.W. Analytic Properties of the Vertex Function in Gauge Theories. 2. Phys. Rev. D 1980, 22, 2550, Erratum in Phys. Rev. D 1981, 23, 3085. [Google Scholar] [CrossRef]
  252. Davydychev, A.I.; Osland, P.; Tarasov, O. Three gluon vertex in arbitrary gauge and dimension. Phys. Rev. D 1996, 54, 4087–4113, Erratum in Phys. Rev. D 1999, 59, 109901. [Google Scholar] [CrossRef] [Green Version]
  253. von Smekal, L.; Hauck, A.; Alkofer, R. A Solution to Coupled Dyson–Schwinger Equations for Gluons and Ghosts in Landau Gauge. Ann. Phys. 1998, 267, 1–60, Erratum in Ann. Phys. 1998, 269, 182. [Google Scholar] [CrossRef] [Green Version]
  254. Binosi, D.; Papavassiliou, J. Gauge invariant Ansatz for a special three-gluon vertex. J. High Energy Phys. 2011, 03, 121. [Google Scholar] [CrossRef] [Green Version]
  255. Gracey, J.; Kißler, H.; Kreimer, D. Self-consistency of off-shell Slavnov-Taylor identities in QCD. Phys. Rev. D 2019, 100, 085001. [Google Scholar] [CrossRef] [Green Version]
  256. Collins, J.C. Renormalization. An Introduction To Renormalization, The Renormalization Group, And The Operator Product Expansion; Cambridge University Press: Cambridge, UK, 1986. [Google Scholar]
  257. Brodsky, S.J.; Shrock, R. Maximum Wavelength of Confined Quarks and Gluons and Properties of Quantum Chromodynamics. Phys. Lett. 2008, B666, 95–99. [Google Scholar] [CrossRef]
  258. Braun, J.; Gies, H.; Pawlowski, J.M. Quark Confinement from Color Confinement. Phys. Lett. 2010, B684, 262–267. [Google Scholar] [CrossRef] [Green Version]
  259. Gao, F.; Qin, S.X.; Roberts, C.D.; Rodriguez-Quintero, J. Locating the Gribov horizon. Phys. Rev. 2018, D97, 034010. [Google Scholar] [CrossRef] [Green Version]
  260. Jackiw, R.; Johnson, K. Dynamical Model of Spontaneously Broken Gauge Symmetries. Phys. Rev. D 1973, 8, 2386–2398. [Google Scholar] [CrossRef]
  261. Jackiw, R. Dynamical Symmetry Breaking. In Proceedings of the 11th International School of Subnuclear Physics: Laws of Hadronic Matter, Erice, Italy, 8–26 July 1973. [Google Scholar]
  262. Papavassiliou, J. Gauge Invariant Proper Selfenergies and Vertices in Gauge Theories with Broken Symmetry. Phys. Rev. D 1990, 41, 3179. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  263. Nakanishi, N. A General survey of the theory of the Bethe–Salpeter equation. Prog. Theor. Phys. Suppl. 1969, 43, 1–81. [Google Scholar] [CrossRef] [Green Version]
  264. Blank, M.; Krassnigg, A. Matrix algorithms for solving (in)homogeneous bound state equations. Comput. Phys. Commun. 2011, 182, 1391–1401. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  265. Wilson, K.G. Quantum field theory models in less than four-dimensions. Phys. Rev. 1973, D7, 2911–2926. [Google Scholar]
  266. Poggio, E.C.; Tomboulis, E.; Tye, S.H.H. Dynamical Symmetry Breaking in Nonabelian Field Theories. Phys. Rev. 1975, D11, 2839. [Google Scholar] [CrossRef]
  267. Denner, A.; Weiglein, G.; Dittmaier, S. Gauge invariance of green functions: Background field method versus pinch technique. Phys. Lett. 1994, B333, 420–426. [Google Scholar] [CrossRef] [Green Version]
  268. Hashimoto, S.; Kodaira, J.; Yasui, Y.; Sasaki, K. The Background field method: Alternative way of deriving the pinch technique’s results. Phys. Rev. D 1994, 50, 7066–7076. [Google Scholar] [CrossRef] [Green Version]
  269. Papavassiliou, J. On the connection between the pinch technique and the background field method. Phys. Rev. 1995, D51, 856–861. [Google Scholar] [CrossRef]
  270. Aguilar, A.C.; Binosi, D.; Iba nez, D.; Papavassiliou, J. New method for determining the quark-gluon vertex. Phys. Rev. 2014, D90, 065027. [Google Scholar] [CrossRef] [Green Version]
  271. Alkofer, R.; Fischer, C.S.; Llanes-Estrada, F.J. Vertex functions and infrared fixed point in Landau gauge SU(N) Yang-Mills theory. Phys. Lett. B 2005, 611, 279–288, Erratum in Phys. Lett. B 2009, 670, 460–461. [Google Scholar] [CrossRef] [Green Version]
  272. Kugo, T.; Ojima, I. Local Covariant Operator Formalism of Nonabelian Gauge Theories and Quark Confinement Problem. Prog. Theor. Phys. Suppl. 1979, 66, 1–130. [Google Scholar] [CrossRef] [Green Version]
  273. Nakanishi, N.; Ojima, I. Covariant Operator Formalism of Gauge Theories and Quantum Gravity; World Scientific Lectures Notes in Physics: Singapore, 1990; Volume 27. [Google Scholar] [CrossRef]
  274. Alkofer, R.; Huber, M.Q.; Schwenzer, K. Infrared singularities in Landau gauge Yang-Mills theory. Phys. Rev. D 2010, 81, 105010. [Google Scholar] [CrossRef] [Green Version]
  275. Huber, M.Q. On non-primitively divergent vertices of Yang–Mills theory. Eur. Phys. J. 2017, C77, 733. [Google Scholar] [CrossRef] [Green Version]
  276. Aguilar, A.C.; Ferreira, M.N.; Oliveira, B.M.; Papavassiliou, J. Schwinger-Dyson truncations in the all-soft limit: A case study. Eur. Phys. J. C 2022, 82, 1068. [Google Scholar] [CrossRef]
  277. Celmaster, W.; Gonsalves, R.J. The Renormalization Prescription Dependence of the QCD Coupling Constant. Phys. Rev. 1979, D20, 1420. [Google Scholar]
  278. Mintz, B.W.; Palhares, L.F.; Sorella, S.P.; Pereira, A.D. Ghost-gluon vertex in the presence of the Gribov horizon. Phys. Rev. 2018, D97, 034020. [Google Scholar] [CrossRef] [Green Version]
  279. Barrios, N.; Peláez, M.; Reinosa, U.; Wschebor, N. The ghost-antighost–gluon vertex from the Curci-Ferrari model: Two-loop corrections. Phys. Rev. D 2020, 102, 114016. [Google Scholar] [CrossRef]
  280. Cyrol, A.K.; Pawlowski, J.M.; Rothkopf, A.; Wink, N. Reconstructing the gluon. SciPost Phys. 2018, 5, 065. [Google Scholar] [CrossRef]
  281. Binosi, D.; Tripolt, R.A. Spectral functions of confined particles. Phys. Lett. B 2020, 801, 135171. [Google Scholar] [CrossRef]
  282. Kern, W.; Huber, M.Q.; Alkofer, R. The spectral dimension as a tool for analyzing non-perturbative propagators. Phys. Rev. 2019, D100, 094037. [Google Scholar] [CrossRef] [Green Version]
  283. Horak, J.; Papavassiliou, J.; Pawlowski, J.M.; Wink, N. Ghost spectral function from the spectral Dyson-Schwinger equation. Phys. Rev. D 2021, 104, 074017. [Google Scholar] [CrossRef]
  284. Horak, J.; Pawlowski, J.M.; Rodríguez-Quintero, J.; Turnwald, J.; Urban, J.M.; Wink, N.; Zafeiropoulos, S. Reconstructing QCD spectral functions with Gaussian processes. Phys. Rev. D 2022, 105, 036014. [Google Scholar] [CrossRef]
  285. Horak, J.; Pawlowski, J.M.; Wink, N. On the complex structure of Yang-Mills theory. arXiv 2022, arXiv:2202.09333. [Google Scholar]
  286. Kallen, G. On the definition of the Renormalization Constants in Quantum Electrodynamics. Helv. Phys. Acta 1952, 25, 417. [Google Scholar] [CrossRef]
  287. Lehmann, H. On the Properties of propagation functions and renormalization contants of quantized fields. Nuovo Cim. 1954, 11, 342–357. [Google Scholar] [CrossRef]
  288. Gross, D.J.; Pisarski, R.D.; Yaffe, L.G. QCD and Instantons at Finite Temperature. Rev. Mod. Phys. 1981, 53, 43. [Google Scholar] [CrossRef]
  289. Jackiw, R.; Templeton, S. How Superrenormalizable Interactions Cure their Infrared Divergences. Phys. Rev. D 1981, 23, 2291. [Google Scholar] [CrossRef] [Green Version]
  290. Appelquist, T.; Pisarski, R.D. High-Temperature Yang-Mills Theories and Three-Dimensional Quantum Chromodynamics. Phys. Rev. D 1981, 23, 2305. [Google Scholar] [CrossRef]
  291. Deser, S.; Jackiw, R.; Templeton, S. Three-Dimensional Massive Gauge Theories. Phys. Rev. Lett. 1982, 48, 975–978. [Google Scholar] [CrossRef]
  292. Cornwall, J.M. HOW d = 3 QCD RESEMBLES d = 4 QCD. Physica 1989, A158, 97–110. [Google Scholar] [CrossRef]
  293. Cornwall, J.M. Exact zero momentum sum rules in d = 3 gauge theory. Nucl. Phys. B 1994, 416, 335–350. [Google Scholar] [CrossRef]
  294. Alexanian, G.; Nair, V.P. A Selfconsistent inclusion of magnetic screening for the quark - gluon plasma. Phys. Lett. 1995, B352, 435–439. [Google Scholar] [CrossRef] [Green Version]
  295. Cornwall, J.M.; Yan, B. String tension and Chern-Simons fluctuations in the vortex vacuum of d = 3 gauge theory. Phys. Rev. D 1996, 53, 4638–4649. [Google Scholar] [CrossRef] [Green Version]
  296. Cornwall, J.M. On the phase transition in D = 3 Yang-Mills Chern-Simons gauge theory. Phys. Rev. D 1996, 54, 1814–1825. [Google Scholar] [CrossRef] [Green Version]
  297. Buchmuller, W.; Philipsen, O. Magnetic screening in the high temperature phase of the standard model. Phys. Lett. 1997, B397, 112–118. [Google Scholar] [CrossRef] [Green Version]
  298. Jackiw, R.; Pi, S.Y. Seeking an even-parity mass term for 3-D gauge theory. Phys. Lett. 1997, B403, 297–303. [Google Scholar] [CrossRef] [Green Version]
  299. Cornwall, J.M. On one-loop gap equations for the magnetic mass in d = 3 gauge theory. Phys. Rev. 1998, D57, 3694–3700. [Google Scholar]
  300. Karabali, D.; Kim, C.j.; Nair, V. On the vacuum wave function and string tension of Yang-Mills theories in (2+1)-dimensions. Phys. Lett. 1998, B434, 103–109. [Google Scholar] [CrossRef] [Green Version]
  301. Eberlein, F. Two loop gap equations for the magnetic mass. Phys. Lett. 1998, B439, 130–136. [Google Scholar] [CrossRef] [Green Version]
  302. Aguilar, A.C.; Binosi, D.; Papavassiliou, J. Nonperturbative gluon and ghost propagators for d=3 Yang-Mills. Phys. Rev. 2010, D81, 125025. [Google Scholar] [CrossRef]
Figure 1. Upper panel: the diagrammatic representation of the conventional gluon self-energy, Π μ ν ( q ) . Bottom panel: the diagrammatic representation of the Q μ a ( q ) B ν b ( q ) , self-energy δ a b Π ˜ μ ν ( q ) ; the grey circles at the end of the gluon lines indicate a background gluon. The corresponding Feynman rules are given in Appendix B of [14].
Figure 1. Upper panel: the diagrammatic representation of the conventional gluon self-energy, Π μ ν ( q ) . Bottom panel: the diagrammatic representation of the Q μ a ( q ) B ν b ( q ) , self-energy δ a b Π ˜ μ ν ( q ) ; the grey circles at the end of the gluon lines indicate a background gluon. The corresponding Feynman rules are given in Appendix B of [14].
Particles 06 00017 g001
Figure 2. Diagrammatic definition of the ghost–gluon scattering kernel, H ν μ ( r , p , q ) . At the tree level, H ν μ 0 = g ν μ .
Figure 2. Diagrammatic definition of the ghost–gluon scattering kernel, H ν μ ( r , p , q ) . At the tree level, H ν μ 0 = g ν μ .
Particles 06 00017 g002
Figure 3. The diagrammatic representation of the three-gluon and ghost–gluon vertices introduced in Equation (33): I Γ α μ ν ( q , r , p ) (first row) and I Γ α ( r , p , q ) (second row). The first term on the r.h.s. indicates the pole-free part, Γ α μ ν ( q , r , p ) or Γ α ( r , p , q ) , while the second denotes the pole term V α μ ν ( q , r , p ) or V α ( r , p , q ) .
Figure 3. The diagrammatic representation of the three-gluon and ghost–gluon vertices introduced in Equation (33): I Γ α μ ν ( q , r , p ) (first row) and I Γ α ( r , p , q ) (second row). The first term on the r.h.s. indicates the pole-free part, Γ α μ ν ( q , r , p ) or Γ α ( r , p , q ) , while the second denotes the pole term V α μ ν ( q , r , p ) or V α ( r , p , q ) .
Particles 06 00017 g003
Figure 4. The coupled system of Schwinger–Dyson equations (SDEs) for the three-gluon and ghost–gluon vertices, I Γ α μ ν ( q , r , p ) and I Γ α ( r , p , q ) , respectively. The orange ellipses represent four-point scattering kernels, denoted by K i j . We omit diagrams containing five-point scattering kernels.
Figure 4. The coupled system of Schwinger–Dyson equations (SDEs) for the three-gluon and ghost–gluon vertices, I Γ α μ ν ( q , r , p ) and I Γ α ( r , p , q ) , respectively. The orange ellipses represent four-point scattering kernels, denoted by K i j . We omit diagrams containing five-point scattering kernels.
Particles 06 00017 g004
Figure 5. The one-particle exchange approximations, K i j 0 , of the kernels K i j appearing in Figure 4.
Figure 5. The one-particle exchange approximations, K i j 0 , of the kernels K i j appearing in Figure 4.
Particles 06 00017 g005
Figure 6. The diagrammatic representation of the coupled system of Bethe–Salpeter equations (BSEs) that governs the evolution of the functions C ( r 2 ) and C ( r 2 ) .
Figure 6. The diagrammatic representation of the coupled system of Bethe–Salpeter equations (BSEs) that governs the evolution of the functions C ( r 2 ) and C ( r 2 ) .
Particles 06 00017 g006
Figure 7. The solutions for C 🟉 ( r ) (purple dot-dashed) and C 🟉 ( r ) (red dashed) obtained from the coupled BSE system of Equation (46).
Figure 7. The solutions for C 🟉 ( r ) (purple dot-dashed) and C 🟉 ( r ) (red dashed) obtained from the coupled BSE system of Equation (46).
Particles 06 00017 g007
Figure 8. Diagrammatic representation of the basic PT rearrangement in the case of quark–antiquark scattering, corresponding to the S-matrix element T q q ¯ q q ¯ of Equation (70) (left), and gluon–gluon scattering, corresponding to T g g g g of Equation (71) (right).
Figure 8. Diagrammatic representation of the basic PT rearrangement in the case of quark–antiquark scattering, corresponding to the S-matrix element T q q ¯ q q ¯ of Equation (70) (left), and gluon–gluon scattering, corresponding to T g g g g of Equation (71) (right).
Particles 06 00017 g008
Figure 9. Top left: Gluon propagator, Δ ( q ) , from lattice simulations of Reference [85] (points) and a fit given by Equation (C11) of [124] (blue continuous). Top right: The auxiliary function 1 + G ( q ) , defined in Equation (12). Bottom left: The renormalization group invariant (RGI) running interaction strength d ^ ( q ) defined in Equation (69), computed using the Δ ( q ) and 1 + G ( q ) shown in the top panels, with α s = 0.27 [71] and Z 1 = 0.9333 [see Section 8]. Bottom right: The corresponding dimensionless RGI interaction I ( q ) , defined in Equation (72).
Figure 9. Top left: Gluon propagator, Δ ( q ) , from lattice simulations of Reference [85] (points) and a fit given by Equation (C11) of [124] (blue continuous). Top right: The auxiliary function 1 + G ( q ) , defined in Equation (12). Bottom left: The renormalization group invariant (RGI) running interaction strength d ^ ( q ) defined in Equation (69), computed using the Δ ( q ) and 1 + G ( q ) shown in the top panels, with α s = 0.27 [71] and Z 1 = 0.9333 [see Section 8]. Bottom right: The corresponding dimensionless RGI interaction I ( q ) , defined in Equation (72).
Particles 06 00017 g009
Figure 10. Lattice data from Reference [86] for the classical form factor, Γ ˜ 1 ( q 2 , r 2 , p 2 ) , of the transversely projected three-gluon vertex in three different kinematic configurations: the soft-gluon ( q = 0 , p 2 = r 2 , orange stars), the symmetric limit ( q 2 = p 2 = r 2 , green diamonds), and the case p 2 = r 2 with q 2 = 2 r 2 (purple circles). In the left panel Γ ˜ 1 ( q 2 , r 2 , p 2 ) is plotted as a function of r, while in the right it is plotted as a function of the Bose symmetry variable s defined in Equation (75).
Figure 10. Lattice data from Reference [86] for the classical form factor, Γ ˜ 1 ( q 2 , r 2 , p 2 ) , of the transversely projected three-gluon vertex in three different kinematic configurations: the soft-gluon ( q = 0 , p 2 = r 2 , orange stars), the symmetric limit ( q 2 = p 2 = r 2 , green diamonds), and the case p 2 = r 2 with q 2 = 2 r 2 (purple circles). In the left panel Γ ˜ 1 ( q 2 , r 2 , p 2 ) is plotted as a function of r, while in the right it is plotted as a function of the Bose symmetry variable s defined in Equation (75).
Particles 06 00017 g010
Figure 11. Lattice data from Reference [84] for L s g ( q ) , compared to the fit for it given by Equation (C12) of [124] (blue continuous curve).
Figure 11. Lattice data from Reference [84] for L s g ( q ) , compared to the fit for it given by Equation (C12) of [124] (blue continuous curve).
Particles 06 00017 g011
Figure 12. Top: SDE governing the momentum evolution of the ghost propagator. Bottom: SDE for the ghost–gluon vertex, I Γ ν ( r , q , p ) .
Figure 12. Top: SDE governing the momentum evolution of the ghost propagator. Bottom: SDE for the ghost–gluon vertex, I Γ ν ( r , q , p ) .
Particles 06 00017 g012
Figure 13. (Left): ghost dressing function F ( q ) obtained from the coupled system of SDEs of Equations (80) and (83) (blue continuous line) compared to the lattice data of References [73,85]. (Right): The corresponding result for B 1 ( r 2 , p 2 , θ ) for arbitrary magnitudes of the anti-ghost and gluon momenta, r and p, respectively, and a representative value of θ = 2 π / 3 for the angle between them. The red dot-dashed curve highlights the soft gluon limit ( p = 0 ).
Figure 13. (Left): ghost dressing function F ( q ) obtained from the coupled system of SDEs of Equations (80) and (83) (blue continuous line) compared to the lattice data of References [73,85]. (Right): The corresponding result for B 1 ( r 2 , p 2 , θ ) for arbitrary magnitudes of the anti-ghost and gluon momenta, r and p, respectively, and a representative value of θ = 2 π / 3 for the angle between them. The red dot-dashed curve highlights the soft gluon limit ( p = 0 ).
Particles 06 00017 g013
Figure 14. Soft gluon limit, B 1 T ( r 2 , 0 , θ ) , of the classical form factor of the ghost–gluon vertex in the Taylor scheme. The points correspond to the lattice data of References [42,43]. The red dot-dashed line shows the SDE solution with the three-gluon vertex dressed according to Equation (78), while the green dashed represents the SDE solution with tree-level three-gluon vertex.
Figure 14. Soft gluon limit, B 1 T ( r 2 , 0 , θ ) , of the classical form factor of the ghost–gluon vertex in the Taylor scheme. The points correspond to the lattice data of References [42,43]. The red dot-dashed line shows the SDE solution with the three-gluon vertex dressed according to Equation (78), while the green dashed represents the SDE solution with tree-level three-gluon vertex.
Particles 06 00017 g014
Figure 15. Lattice data for the gluon propagator in the deep infrared. The data displayed correspond to the two lattice setups with the largest volumes of [49], namely, V = 72 4 (left) and V = 80 4 (right). The red dashed lines are smooth fits from which the position of the maximum can be estimated.
Figure 15. Lattice data for the gluon propagator in the deep infrared. The data displayed correspond to the two lattice setups with the largest volumes of [49], namely, V = 72 4 (left) and V = 80 4 (right). The red dashed lines are smooth fits from which the position of the maximum can be estimated.
Particles 06 00017 g015
Figure 16. The SDE of the three-gluon vertex at the one-loop dressed level. Diagrams ( e 1 ) and ( e 2 ) are the gluon and the ghost triangle contributions entering the skeleton expansion of the three-gluon vertex.
Figure 16. The SDE of the three-gluon vertex at the one-loop dressed level. Diagrams ( e 1 ) and ( e 2 ) are the gluon and the ghost triangle contributions entering the skeleton expansion of the three-gluon vertex.
Particles 06 00017 g016
Figure 17. The ghost triangle present in the three-gluon vertex SDE (top) and the ghost loop contributing to the gluon propagator in the corresponding equation (middle). The infrared divergences arising from these diagrams are connected through the Slavnov–Taylor identity (STI) of Equation (23), as shown schematically in the bottom panel.
Figure 17. The ghost triangle present in the three-gluon vertex SDE (top) and the ghost loop contributing to the gluon propagator in the corresponding equation (middle). The infrared divergences arising from these diagrams are connected through the Slavnov–Taylor identity (STI) of Equation (23), as shown schematically in the bottom panel.
Particles 06 00017 g017
Figure 18. Comparison of L s g ( r ) (blue continuous) from Figure 11 and d Δ 1 ( r ) / d r 2 (red dashed) resulting from the fit for Δ ( r ) of Figure 9. Note that both display the characteristic features of infrared suppression with respect to their tree-level values (which is 1 for both quantities), zero-crossing, and logarithmic divergence at the origin.
Figure 18. Comparison of L s g ( r ) (blue continuous) from Figure 11 and d Δ 1 ( r ) / d r 2 (red dashed) resulting from the fit for Δ ( r ) of Figure 9. Note that both display the characteristic features of infrared suppression with respect to their tree-level values (which is 1 for both quantities), zero-crossing, and logarithmic divergence at the origin.
Particles 06 00017 g018
Figure 19. SDE for the ghost–gluon scattering kernel, H μ ν ( r , q , p ) . We omit a diagram containing a 1PI four-point function.
Figure 19. SDE for the ghost–gluon scattering kernel, H μ ν ( r , q , p ) . We omit a diagram containing a 1PI four-point function.
Particles 06 00017 g019
Figure 20. Diagrammatic representation of the cancellation of the infrared divergences originating from massless ghost loops in Equation (109) to yield a finite C ( 0 ) . The red cross indicates that the overall ghost momentum is factored out before being set to zero.
Figure 20. Diagrammatic representation of the cancellation of the infrared divergences originating from massless ghost loops in Equation (109) to yield a finite C ( 0 ) . The red cross indicates that the overall ghost momentum is factored out before being set to zero.
Particles 06 00017 g020
Figure 21. W ( r ) obtained using the approximation Equation (117) based on the observed planar degeneracy of the three-gluon vertex in its kernel (blue solid curve) together with uncertainty estimate (blue band).
Figure 21. W ( r ) obtained using the approximation Equation (117) based on the observed planar degeneracy of the three-gluon vertex in its kernel (blue solid curve) together with uncertainty estimate (blue band).
Particles 06 00017 g021
Figure 22. Left: Result for C ( r ) (black continuous line) obtained from Equation (109) using the W ( r ) shown in Figure 21, the fits to lattice data for Δ ( r ) and L s g ( r ) are shown in Figure 9 and Figure 18, respectively, and the SDE solution for F ( r ) shown in Figure 13. The points are obtained using for L s g ( r ) the data in Reference [84], with error bars denoting the error propagated from L s g and W . The green band is obtained by fitting the upper and lower bounds of the points and guiding the eye to the typical error associated with C ( r ) . Right: The C ( r ) of the left panel is compared to the BSE prediction C 🟉 ( r ) (purple dot-dashed and error band) of Figure 7.
Figure 22. Left: Result for C ( r ) (black continuous line) obtained from Equation (109) using the W ( r ) shown in Figure 21, the fits to lattice data for Δ ( r ) and L s g ( r ) are shown in Figure 9 and Figure 18, respectively, and the SDE solution for F ( r ) shown in Figure 13. The points are obtained using for L s g ( r ) the data in Reference [84], with error bars denoting the error propagated from L s g and W . The green band is obtained by fitting the upper and lower bounds of the points and guiding the eye to the typical error associated with C ( r ) . Right: The C ( r ) of the left panel is compared to the BSE prediction C 🟉 ( r ) (purple dot-dashed and error band) of Figure 7.
Particles 06 00017 g022
Table 1. The different types of gluon propagators of the background field method (BFM), together with their diagrammatic representations, symbols, corresponding self-energies, and the background quantum identities (BQIs) that relate them to the conventional propagator.
Table 1. The different types of gluon propagators of the background field method (BFM), together with their diagrammatic representations, symbols, corresponding self-energies, and the background quantum identities (BQIs) that relate them to the conventional propagator.
  External LegsDiagrammatic RepresentationSymbolSelf-EnergyBQI
Q μ a ( q ) Q ν b ( q ) Particles 06 00017 i001
i δ a b Δ μ ν ( q )
Π μ ν ( q )
Q μ a ( q ) B ν b ( q ) Particles 06 00017 i002
i δ a b Δ ˜ μ ν ( q )
Π ˜ μ ν ( q )
Δ ˜ ( q ) = Δ ( q ) 1 + G ( q )
B μ a ( q ) B ν b ( q ) Particles 06 00017 i003
i δ a b Δ ^ μ ν ( q )
Π ^ μ ν ( q )
Δ ^ ( q ) = Δ ( q ) [ 1 + G ( q ) ] 2
Table 2. Statistical results of the χ 2 analysis for the three hypotheses given in Equation (92) for the form factor B 1 . For each case (first column), we give the corresponding χ j 2 computed from Equation (91) (second column), probability P j computed from Equation (93) (third row), and the same P j expressed in terms of confidence levels σ (fourth row).
Table 2. Statistical results of the χ 2 analysis for the three hypotheses given in Equation (92) for the form factor B 1 . For each case (first column), we give the corresponding χ j 2 computed from Equation (91) (second column), probability P j computed from Equation (93) (third row), and the same P j expressed in terms of confidence levels σ (fourth row).
Case (j) χ j 2 P j Confidence Level in σ
171.37 4.0 × 10 7 5.1
23.3991 − 1.8 × 10 6 2.2 × 10 6
350.03 5.8 × 10 4 3.4
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ferreira, M.N.; Papavassiliou, J. Gauge Sector Dynamics in QCD. Particles 2023, 6, 312-363. https://doi.org/10.3390/particles6010017

AMA Style

Ferreira MN, Papavassiliou J. Gauge Sector Dynamics in QCD. Particles. 2023; 6(1):312-363. https://doi.org/10.3390/particles6010017

Chicago/Turabian Style

Ferreira, Mauricio Narciso, and Joannis Papavassiliou. 2023. "Gauge Sector Dynamics in QCD" Particles 6, no. 1: 312-363. https://doi.org/10.3390/particles6010017

APA Style

Ferreira, M. N., & Papavassiliou, J. (2023). Gauge Sector Dynamics in QCD. Particles, 6(1), 312-363. https://doi.org/10.3390/particles6010017

Article Metrics

Back to TopTop