Landau gauge Yang-Mills correlation functions
Anton K. Cyrol,1 Leonard Fister,2 Mario Mitter,1 Jan M. Pawlowski,1, 3 and Nils Strodthoff1, 4 1Institut f¨urTheoretische Physik, Universit¨atHeidelberg, Philosophenweg 16, 69120 Heidelberg, Germany 2Institut de Physique Th´eorique,CEA Saclay, F-91191 Gif-sur-Yvette, France 3ExtreMe Matter Institute EMMI, GSI, Planckstr. 1, D-64291 Darmstadt, Germany 4Nuclear Science Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA We investigate Landau gauge SU(3) Yang-Mills theory in a systematic vertex expansion scheme for the effective action with the functional renormalisation group. Particular focus is put on the dynamical creation of the gluon mass gap at non-perturbative momenta and the consistent treatment of quadratic divergences. The non-perturbative ghost and transverse gluon propagators as well as the momentum-dependent ghost-gluon, three-gluon and four-gluon vertices are calculated self- consistently with the classical action as only input. The apparent convergence of the expansion scheme is discussed and within the errors, our numerical results are in quantitative agreement with available lattice results.
PACS numbers: 12.38.Aw, 12.38.Gc
I. INTRODUCTION tex expansion scheme. Therefore, the results [9, 33] for the YM correlations functions give no access to system- The past two decades have seen tremendous progress in atic error estimates. In general, many applications of the description of QCD with functional approaches such functional methods to bound states and the QCD phase as the functional renormalisation group (FRG), Dyson- diagram use such mixed approaches, where part of the Schwinger equations (DSE), and n-particle irreducible correlation functions is deduced from phenomenological methods (nPI). These approaches constitute ab initio constraints or other external input. Despite the huge descriptions of QCD in terms of quark and gluon cor- success of mixed approaches, a full ab initio method is relation functions. The full correlation functions satisfy wanted for some of the most pressing open questions of a hierarchy of loop equations that are derived from the strongly-interacting matter. The phase structure of QCD functional FRG, DSE and nPI relations for the respec- at large density is dominated by fluctuations and even a tive generating functionals. By now, systematic compu- partial phenomenological parameter fixing at vanishing tational schemes are available, which can be controlled by density is bound to lead to large systematic errors [34]. apparent convergence. In the present work on pure Yang- The same applies to the details of the hadron spectrum, Mills (YM) theory we complement the work in quenched in particular with regard to the physics of the higher QCD [1], where such a systematic expansion scheme has resonances, which requires knowledge about correlation been put forward within the FRG. Equipped with such a functions deep in the complex plane. controlled expansion, functional approaches to QCD are In the present work we perform a systematic vertex specifically interesting at finite temperatures and large expansion of the effective action of Landau gauge YM density, where reliable ab initio theoretical predictions theory within the functional renormalisation group ap- and experimental results are missing at present. proach, discussed in Sec. II. The current approximation Most progress with functional approaches has been is summarised in Sec. II B, which also includes a compari- made in Landau gauge QCD, which has many convenient son to approximations used in other works. This ab initio properties for non-perturbative numerical computations. approach starts from the classical action. Therefore the Applications of functional methods include the first-ever only parameter is the strong coupling constant αs at a arXiv:1605.01856v2 [hep-ph] 9 Sep 2016 calculation of qualitative non-perturbative Landau gauge large, perturbative momentum scale. The most distinct propagators as well as investigations of the phase struc- feature of YM theory is confinement, which is reflected ture of QCD. For reviews see [2–13], for applications to by the creation of a gluon mass gap in Landau gauge. Yang-Mills theory see e.g. [14–27], and e.g. [28] for re- We discuss the necessity of consistent infrared irregular- lated studies in the Coulomb gauge. The formal, alge- ities as well as mechanisms for the generation of a mass braic, and numerical progress of the past decades sets gap in Sec. III. Numerical results from a parameter-free, the stage for a systematic vertex expansion scheme of self-consistent calculation of propagators and vertices are Landau gauge QCD. Quantitative reliability is then ob- presented in Sec. IV. Particular focus is put on the im- tained with apparent convergence [1] as well as further portance of an accurate renormalisation of the relevant systematic error controls inherent to the method, see e.g. vertices. We compare with corresponding DSE and lat- [5, 29–32]. In the aforementioned quenched QCD in- tice results and discuss the apparent convergence of the vestigation [1], the gluon propagator was taken from a vertex expansion. Finally, we present numerical evidence separate FRG calculation in [9, 33]. This gluon propa- for the dynamic mass gap generation in our calculation. gator shows quantitative agreement with the lattice re- Further details, including a thorough discussion of the sults, but has been obtained within an incomplete ver- necessary irregularities, can be found in the appendices. 2
II. FRG FLOWS FOR YANG-MILLS THEORY IN A VERTEX EXPANSION
Functional approaches to QCD and Yang-Mills theory are based on the classical gauge fixed action of SU(3) Yang-Mills theory. In general covariant gauges in four full. mom. dep. full. mom. dep. sym. point and dimensions it is given by tadpole config. Z Z Z 1 a a 1 a 2 a ab b FIG. 1. Approximation for the effective action. Only clas- Scl = Fµν Fµν + (∂µAµ) c¯ ∂µDµ c . x 4 2ξ x − x sical tensor structures are included. See Fig. 2 for diagrams (1) that contribute to the individual propagators and vertices.
Here, ξ denotes the gauge fixing parameter, which is R R 4 where R = R d4p/(2π)4 . Here we have introduced the taken to zero in the Landau gauge and x = d x. The p field strength tensor and covariant derivative are given RG-time t = ln(k/Λ) with a reference scale Λ , which is by typically chosen as the initial UV cutoff Λ . Although a a a abc b c this flow equation comes in a simple one-loop form, it F = ∂µA ∂ν A + gf A A , µν ν − µ µ ν provides an exact relation due to the presence of the full field-dependent propagator, ab ab abc c D = δ ∂µ gf A , (2) µ − µ 1 G [φ](p, q) = (p, q) , (7) using the fundamental generators T a, defined by k (2) Γk [φ] + Rk 1 T a,T b = if abcT c , tr T aT b = δab . (3) on its right-hand side. Furthermore, the flow is infrared 2 and ultraviolet finite by construction. Via the integration In general, our notation follows the one used in the works of momentum shells in the Wilsonian sense, it connects [1, 35, 36] of the fQCD collaboration [37]. the ultraviolet, bare action Scl = Γk→Λ→∞ with the full quantum effective action Γ = Γk→0 . The flow equations for propagators and vertices are A. Functional Renormalisation Group obtained by taking functional derivatives of (6). At the vacuum expectation values, these derivatives give equa- (n) n n We use the functional renormalisation group approach tions for the 1PI correlation functions Γk = δ Γk/δφ , as a non-perturbative tool to investigate Yang-Mills the- which inherit the one-loop structure of (6). As the cutoff- ory. The FRG is built on a flow equation for the one- derivative of the regulator functions, ∂tRk, decays suffi- particle irreducible (1PI) effective action or free energy ciently fast for large momenta, the momentum integra- of the theory, the Wetterich equation [38]. It is based tion in (6) effectively receives only contributions for mo- 2 2 on Wilson’s idea of introducing an infrared momentum menta p . k . Furthermore, the flow depends solely cutoff scale k. Here, this infrared regularisation of the on dressed vertices and propagators, leading to a consis- gluon and ghost fluctuations is achieved by modifying tent RG and momentum scaling for each diagram result- the action Scl Scl + ∆Sk with ing from derivatives of (6). Despite its simple structure, → the resulting system of equations does not close at a fi- Z Z 1 a ab b a ab b nite number of correlation functions. In general, higher ∆Sk = A R A + c¯ R c . (4) 2 µ k,µν ν k (n+2) x x derivatives up to the order Γk of the effective action appear on the right hand side of the functional relations The regulator functions Rk are momentum-dependent (n) masses that suppress the corresponding fluctuations be- for the correlation functions Γk . low momentum scales p2 k2 and vanish in the ultra- violet for momenta p2 ≈k2. See App. E for details on B. Vertex expansion of the effective action the regulators used in the present work. Consequently, the effective action, Γk[φ] , is infrared regularised, where φ denotes the superfield The structural form of the functional equations dis- cussed in the previous section necessitates the use of ap- φ = (A, c, c¯) . (5) proximations in most practical application. One system- atic expansion scheme is the vertex expansion, i.e. an The fluctuations of the theory are then successively taken expansion of the effective action in terms of 1PI Green’s into account by integrating the flow equation for the ef- functions. This yields an infinite tower of coupled equa- fective action, see e.g. [39, 40], tions for the correlation functions that has to be trun- Z Z cated at a finite order. This expansion scheme allows a 1 µν ba ab ba ∂tΓk[φ] = Gk,ab[φ] ∂tRk,µν Gk [φ] ∂tRk . (6) systematic error estimate in terms of apparent conver- p 2 − p gence upon increasing the expansion order or improving 3
vertices has been summarised as 3 and , which 1 A ,cl Acc,¯ cl ∂t − = + are listed explicitly in App. D. InT the case ofT the transver- sally projected ghost-gluon vertex, Acc,¯ cl represents al- ready a full basis whereas a full basisT for the transversally 1 1 ∂t − = 2 − − 2 projected three-gluon vertex consists of four elements. However, the effect of non-classical tensor structures has been found to be subleading in this case [41]. ∂t = + perm. − − The most important four-point function is given by the four-gluon vertex, which appears already on the classical level. Similarly to the three-gluon vertex, we approxi- ∂t = +2 + + perm. − mate it with its classical tensor structure
(4) abcd abcd [ΓA4 ]µνρσ(p, q, r) = ZA4,⊥(¯p)[ A4,cl]µνρσ(p, q, r) , (10) ∂t = + + 2 + perm. − − T see App. D for details. The dressing function of the four- gluon vertex is approximated from its momentum depen- FIG. 2. Diagrams that contribute to the truncated flow of dence at the symmetric point via the average momen- p propagators and vertices. Wiggly (dotted) lines correspond tump ¯ p2 + q2 + r2 + (p + q + r)2/2 , which has been to dressed gluon (ghost) propagators, filled circles denote shown≡ to be a good approximation of the full momen- dressed (1PI) vertices and regulator insertions are represented by crossed circles. Distinct permutations include not only tum dependence [42, 43]. To improve this approximation (anti-)symmetric permutations of external legs but also per- further, we additionally calculate the momentum depen- mutations of the regulator insertions. dence of the four-gluon dressing function ZA4,⊥( p , q , t) on the special configuration (p, q, r) = (p, q, p) .| We| | | use this special configuration exclusively in the− tadpole di- further approximations for example in the momentum agram of the gluon propagator equation, cf. Sec. IV B. resolution or tensor structures of the included correla- We show the difference between the special configura- tion functions. We discuss the convergence of the vertex tion and the symmetric point approximation in the ap- expansion in Sec. IV B. pendix in Fig. 11. Although the four-gluon vertex has Here we calculate the effective action of SU(3) Yang- been the subject of several studies [42–46], no fully con- Mills theory in Landau gauge within a vertex expansion, clusive statements about the importance of additional see Fig. 1 for a pictorial representation. The diagrams non-classical tensors structures are available. contributing to the resulting equations of the constituents In summary we have taken into account the propa- of our vertex expansion are summarised graphically in gators and the fully momentum-dependent classical ten- Fig. 2. The lowest order contributions in this expan- sor structures of the three-point functions, as well as se- sion are the inverse gluon and ghost propagators param- lected momentum-configurations of the gluon four-point eterised via function, see the paragraph above, and App. D. For a comparison of the current approximation with that used (2) ab 2 ab ⊥ 1 ab [ΓAA]µν (p) = ZA(p) p δ Πµν (p) + δ pµpν , in other functional works one has to keep in mind that ξ FRG, Dyson-Schwinger or nPI equations implement dif-
(2) ab 2 ab ferent resummation schemes. Thus, even on an identical [Γcc¯ ] (p) = Zc(p) p δ , (8) approximation level of a systematic vertex expansion, the included resummations differ. with dimensionless scalar dressing functions 1/ZA and ⊥ 2 In the present work we solve the coupled system of 1/Zc. Here, Πµν (p) = δµν pµpν /p denotes the corre- sponding projection operator.− We use this splitting in all momentum-dependent classical vertex structures and tensor structures with canonical momentum dimension propagators. In former works with functional methods, and dimensionless dressings also for the higher order ver- see e.g. [14–23, 41–53], only subsets of these correlation tices. functions have been coupled back. A notable exception On the three-point level we include the full transverse is [26], where a similar self-consistent approximation has ghost-gluon vertex and the classical tensor structure of been used for three-dimensional Yang-Mills theory. the three-gluon vertex
(3) abc abc C. Modified Slavnov-Taylor identities and [ΓAcc¯ ]µ (p, q) = ZAcc,¯ ⊥( p , q , t)[ Acc,¯ cl]µ (p, q) , | | | | T transversality in Landau gauge (3) abc abc [Γ 3 ]µνρ(p, q) = ZA3,⊥( p , q , t)[ A3,cl]µνρ(p, q) . (9) A | | | | T In Landau gauge, the dynamical system of correlation Here, the momentum p (q) corresponds to the indices a functions consists only of the transversally projected cor- (b) and t denotes the cosine of the angle between the relators [9]. Those with at least one longitudinal gluon momenta p and q . The classical tensor structure of the leg do not feed back into the dynamics. To make these 4 statements precise, it is useful to split correlation func- Consequences of the STIs & mSTIs tions into purely transverse components and their com- plement with at least one longitudinal gluon leg. The For the purposes of this work, the most important ef- (n) purely transverse vertices Γ⊥ , are defined by attach- fect of the modification of the STIs due to the regulator ing transverse projection operators to the corresponding term is that it leads to a non-vanishing gluon mass pa- gluon legs, rameter [54],
h (n)i h i h iab Γ Π⊥ Π⊥ Γ(n) , (11) (2) ab 2 ⊥ µ1ν1 µnA νnA ∆mSTI ΓAA ∝ δ δµν α(k) k . (16) µ1···µnA ≡ ··· ν1···νnA µν where nA is the number of gluon legs and group indices At k = 0, where the regulators vanish, this modification and momentum arguments have been suppressed for the disappears, as the mSTIs reduce to the standard STIs. In sake of brevity. This defines a unique decomposition of particular, this entails that, at k = 0, the inverse longi- (2) n-point functions into tudinal gluon propagator, ΓAA,L, reduces to the classical one, solely determined by the gauge fixing term (n) (n) (n) Γ = Γ + ΓL , (12) ⊥ (2) ab (2) ab pµ [ΓAA,L]µν (p) [SAA,L]µν (p) = 0 . (17) (n) − where the longitudinal vertices ΓL , have at least one longitudinal gluon leg. Consequently, they are always This provides a unique condition for determining the projected to zero by the purely transverse projection op- value of the gluon mass parameter (16) at the ultraviolet erators of (11). initial scale Λ. However, it can only serve its purpose, if Functional equations for the transverse correlation the longitudinal system is additionally solved. functions close in the Landau gauge, leading to the struc- One further conclusion from (15) is that the mSTIs do ture [9], not constrain the transverse correlation functions with- out further input. This fact is not in tension with one of the main applications of STIs in perturbation theory, i.e. Γ(n) = Diag[ Γ(n) ] . (13) ⊥ { ⊥ } relating the running of the relevant vertices of Yang-Mills theory that require renormalisation. As Yang-Mills the- In (13) Diag stands for diagrammatic expressions of ei- ory is renormalisable, only the classical vertex structures ther integrated FRG, Dyson-Schwinger or nPI equations. are renormalised and hence the renormalisation functions Equation (13) follows from the fact that all internal legs of their transverse and longitudinal parts have to be iden- are transversally projected by the Landau gauge gluon tical. propagator. Hence, by using transverse projections for As an instructive example we consider the ghost-gluon the external legs one obtains (13). In contradistinction vertex. For this example and the following discussions to this, the functional equations for the vertices with at (n) we evaluate the STIs within the approximation used in least one longitudinal gluon leg, ΓL , are of the form the present work: on the right hand side of the STIs we only consider contributions from the primitively diver- (n) (n) (n) ΓL = Diag[ ΓL , Γ⊥ ] . (14) gent vertices. In particular, this excludes contributions { } { } from the two-ghost–two-gluon vertex. The ghost-gluon In other words, the solution of the functional equations vertex can be parameterised with two tensor structures, (14) for Γ(n) requires also the solution of the transverse L (3) abc abch i set of equations (13). [ΓAcc¯ ]µ (p, q) = if qµZAcc,¯ cl(p, q) + pµZAcc,¯ ncl(p, q) . In the present setting, gauge invariance is encoded in (18) modified Slavnov-Taylor identities (mSTIs) and Ward- Takahashi identities (mWTIs). They are derived from In (18) we have introduced two dressing functions ZAcc,¯ cl the standard Slavnov-Taylor identities (STIs) by includ- and ZAcc,¯ ncl as functions of the gluon momentum p and ing the gauge or BRST variations of the regulator terms, anti-ghost momentum q . In a general covariant gauge see [5, 14, 54–57] for details. The mSTIs are of the only ZAcc,¯ cl requires renormalisation. Similar splittings schematic form into a classical tensor structure and the rest can be used in other vertices. Trivially, this property relates the per- (n) (n) (n) ΓL = mSTI[ ΓL , Γ⊥ ,Rk] , (15) turbative RG-running of the transverse and longitudi- { } { } nal projections of the classical tensor structures. Then, which reduce to the standard STIs in the limit of van- the STIs can be used to determine the perturbative RG- ishing regulator, Rk 0. The STIs and mSTIs have a running of the classical tensor structures, leading to the similar structure as (14)≡ and can be used to obtain infor- well-known perturbative relations mation about the longitudinal part of the correlators. Al- 2 2 ternatively, they provide a non-trivial consistency check ZAcc,¯ cl ZA3,cl ZA4,cl 2 = 3 = 2 , (19) for approximate solutions of (14). Zc ZA ZA ZA 5 at the renormalisation scale µ. Consequently, (19) al- variant gauges. Furthermore, in Yang-Mills theory for- lows for the definition of a unique renormalised two-loop mulated in covariant gauges, the gapping of the gluon coupling αs(µ) from the vertices. relative to the ghost is necessary and sufficient for pro- The momentum dependent STIs can also be used to ducing a confining potential for the corresponding order make the relation (19) momentum-dependent. Keeping parameter, the Polyakov loop. Hence, understanding the only the classical tensor structures, we are led to the details of the dynamical generation of a gluon mass gap momentum dependent running couplings gives insight into the confinement mechanism. This relation holds for all potential infrared closures of Z2 (p) 1 Acc,¯ ⊥ the perturbative Landau gauge. The standard infrared αAcc¯ (p) = 2 , 4π ZA(p) Zc (p) closure corresponds to a full average over all Gribov re- gions. This leads to the standard Zinn-Justin equation 2 1 ZA3,⊥(p) as used in the literature, e.g. [4]. In turn, the restriction αA3 (p) = 3 , to the first Gribov regime can be implemented within the 4π ZA(p) refined Gribov-Zwanziger formalism, e.g. [63–67], that
1 ZA4,⊥(p) leads to infrared modifications of the STIs. In the follow- α 4 (p) = , (20) A 4π Z2 (p) ing we discuss the consequences of the standard STIs, a A discussion of the refined Gribov-Zwanziger formalism is where the used transverse projection is indicated by the deferred to future work. subscript , for details see App. D. Additionally, the vertices appearing⊥ in (20) are evaluated at the symmet- ric point, see Sec. IV B for the precise definition. The A. Gluon mass gap and irregularities STIs and two-loop universality demand that these run- ning couplings become degenerate at large perturbative In order to study the dynamical generation of the mass momentum scales, where the longitudinal and transverse gap, we first discuss the consequences of the STI for the parts of the vertices agree. longitudinal gluon two point function (17). It states that In Landau gauge, the ghost-gluon vertex is not renor- no quantum fluctuations contribute to the inverse lon- malised on specific momentum configurations, and we gitudinal gluon propagator, i.e. the longitudinal gluon can alternatively define a running coupling from the wave propagator is defined by the gauge fixing term. There- function renormalisation of ghost and gluon [15, 58], fore, the dynamical creation of a gluon mass gap requires different diagrammatic contributions to the longitudinal 1 g2 αs(p) = 2 . (21) and transverse gluon mass parameter. The discussion 4π ZA(p)Zc (p) of the prerequisites for meeting this condition is qualita- Note that the momentum-dependence of the running cou- tively different for the scaling and the decoupling solu- pling (21) does not coincide with that of the correspond- tions. Hence, these two cases are discussed separately. ing running couplings obtained from other vertices, i.e. The scaling solution is characterised by the infrared (20). This is best seen in the ratio αAcc¯ (p)/αs(p) = behaviour [15, 18, 69–75] 2 2 ZAcc,¯ ⊥(p)/g . In this context we also report on an im- 2 2 κ lim Zc(p ) (p ) , portant result for the quark-gluon vertex coupling, p→0 ∝ 2 1 ZAqq,¯ ⊥(p) 2 2 −2 κ αAqq¯ (p) = 2 , (22) lim ZA(p ) ∝ (p ) , (23) 4π ZA(p)Zq(p) p→0 with the dressing function of the classical tensor struc- 2 with the scaling coefficient 1/2 < κ < 1. A simple calcu- ture of the quark-gluon vertex ZAqq,¯ ⊥(p) and the quark lation presented in App. A shows that the ghost loop with dressing function 1/Zq(p) [1]. The solution of the cor- an infrared constant ghost-gluon vertex and scaling ghost responding STI reveals that the quark-gluon vertex cou- propagator is already capable of inducing a splitting in pling αAqq¯ agrees perturbatively with αs(p) in (21), and the longitudinal and transverse gluon mass parameter. hence it differs from the other vertex couplings in (20). Next we investigate the decoupling solution, e.g. [21, Note that the present truncation only considers contri- 22], which scales with butions from primitively divergent vertices. Accordingly, 2 the two-quark–two-ghost vertex contribution in the STI lim Zc(p ) ∝ 1 , for the quark-gluon vertex, see e.g. [4], has been dropped. p→0
2 2 −1 lim ZA(p ) (p ) , (24) p→0 ∝ III. CONFINEMENT, GLUON MASS GAP, AND IRREGULARITIES at small momenta. Assuming vertices that are regular in the limit of one vanishing gluon momentum, one finds It has been shown in [40, 59–62] that a mass gap in that all diagrammatic contributions to the longitudinal the gluon propagator signals confinement in QCD in co- and transverse gluon mass parameter are identical. For 6
3 FRG, scaling FRG, scaling FRG, decoupling FRG, decoupling 10 Sternbeck et al. Sternbeck et al. 2
1
1 ghost propagator dressing gluon propagator dressing
0 0.1 1 10 0.1 1 10 p [GeV] p [GeV]
FIG. 3. Gluon dressing 1/ZA (left) and ghost dressing 1/Zc (right) in comparison to the lattice results from [68]. The scale setting and normalisation procedures are described in App. F. example, if the ghost loop were to yield a non-vanishing pearing in the gluon propagator equation. Consequently, contribution to the gluon mass gap, the ghost-gluon ver- if all vertices were regular, no gluon mass gap would be tex would have to be a function of the angle θ = arccos(t) created. In particular, regular vertices would entail the between the gluon and anti-ghost momenta p and q, absence of confinement. The necessity of irregularities for the creation of a gluon mass gap was already realised (3) abc (3) abc lim [ΓAcc¯ ]µ ( p , q , t) = [ΓAcc¯ ]µ (0, q , t) , (25) by Cornwall [76]. |p|→0 | | | | | | even in the limit of vanishing gluon momentum p 0 . In the light of these findings it is interesting to in- Since the above limit depends on the angle, the vertex| | → is vestigate the consistency of irregularities with further irregular. See App. A for more details on this particular Slavnov-Taylor identities. Therefore, we consider the case. Similar conclusions can be drawn for all vertices ap- Slavnov-Taylor identity of the three-gluon vertex, e.g. [4],
(3) abc abc 1 ˜ 2 2 ⊥ ˜ 2 2 ⊥ irρ[ΓA3 ]µνρ(p, q, r) f 2 Gµσ(p, q)q ZA(q )Πσν (q) Gνσ(q, p)p ZA(p )Πσµ(p) , (26) ∝ Zc(r ) −
where G˜µν relates to the ghost-gluon vertex via three-gluon vertex projected with one non-zero longitudi- nal leg rρ . Although this momentum configuration does (3) abc abc ˜ [ΓAcc¯ ]µ (p, q) = igf qν Gµν (p, q) . (27) not enter the gluon mass gap directly, crossing symme- try implies the necessary irregularity. In summary, these ˜ For a regular Gµν in the limit p 0 in (26), the scaling arguments illustrate that also the three-gluon vertex STI → solution leads to a singular contribution of the type is consistent with the necessity of irregularities for both 2 1−2κ ˜ ⊥ types of solutions. lim(p ) Gνσ(q, 0) Πσµ(p) + regular , (28) p→0 We close the discussion of vertex irregularities with the remark that the infrared modification of the propagator- where κ is the scaling coefficient from (23). This is con- STI in the refined Gribov-Zwanziger formalism may re- sistent with the expected scaling exponent of the three- move the necessity for irregularities in the vertices. gluon vertex in this limit [74]. In the same limit, the decoupling solution leads to a singular contribution of the form B. Origin of irregularities ˜ ⊥ lim Gνσ(q, 0) Πσµ(p) + regular . (29) p→0 As discussed in the previous section, self-consistency ⊥ Since the transverse projector Πσµ(p) introduces an an- in terms of the Slavnov-Taylor identities entails a corre- gular dependence in the limit p 0 , the STI again spondence between the dynamical generation of a gluon demands an irregularity in limit of→ one vanishing mo- mass gap and the presence of irregularities. But the STIs mentum. Note that this is just a statement about the do not provide a mechanism for the creation of irregular- 7
8 ] FRG, scaling 0 -2 10 FRG, decoupling 6 Sternbeck et al. -2 2 10 1.5 α 4 Acc -4 α 3 10 1 A α 4 running couplings A 2
gluon propagator [GeV -6 10 0.5 0 0.1 1 10 0.1 1 10 100 p [GeV] p [GeV]
FIG. 4. Left: Gluon propagator in comparison to the lattice results from [68]. Right: Effective running couplings defined in (20) as obtained from different Yang-Mills vertices as function of the momentum. ities, the gluon mass gap, and in turn confinement. sufficient for creating a physical mass gap in the gluon. In the scaling solution, (23), the irregularities arise nat- This scenario provides a direct relation of confinement urally from the non-trivial scaling. Hence they are tightly and spontaneous symmetry breaking. Therefore it is pos- linked to the original Kugo-Ojima confinement scenario sibly connected to the presence of resonances that are [77], that requires the non-trivial scaling. Note, however, triggered in the longitudinal sector of the theory, where that this simply links different signatures of confinement they do not spoil the gapping of the completely trans- but does not reveal the mechanism at work. verse sector. A purely longitudinal massless mode, as a For the decoupling solution (24), we want to discuss source for irregularities in the gluonic vertices, has been two possible scenarios. In the first scenario, the irregu- worked out in [80, 81], for a short summary see [82]. As larities are generated in the far infrared. A second possi- a consequence, an irregularity appears in the purely lon- bility is that they are triggered via a condensate and/or a gitudinal three-gluon vertex in a way that preserves the resonance, providing a direct connection of confinement corresponding Slavnov-Taylor identity. The creation of and spontaneous symmetry breaking. a purely transverse background and the presence of lon- gitudinal massless mode would then be two sides of the In the first scenario it is sufficient to focus on ghost same coin. Furthermore, the longitudinal resonance has loops as possible sources of such irregularities, since the to occur at about the same scale as the gluon conden- gluonic diagrams decouple from the infrared dynamics sate, in order to trigger the correct gluon mass gap. A due to the gluon mass gap. This is a seemingly appealing more detailed discussion and computation of this scenario scenario as it is the dynamical ghost that distinguishes cannot be assessed in the purely transverse system and confining Yang-Mills theory from e.g. QED. However, in is therefore deferred to future work. the decoupling solution (24) both, the ghost-gluon ver- tex as well as the ghost propagator, have infrared finite quantum corrections: no ghost-loops contribute to their equation and (infrared) constant dressing functions can C. The purely transverse system be assumed for both. As a consequence the ghost loop contributions to correlation functions have the same in- In this work we restrict ourselves to a solution of the frared structure as perturbative ghost-loop contributions. purely transverse system (13), which is closed. The only However, none of these perturbative ghost loops yields relevant UV parameters in this system are the strong the necessary irregularities, see App. B for an explicit coupling and the transverse gluon mass parameter. In calculation. the UV the transverse mass parameter agrees with the In the second scenario, the generation of irregulari- longitudinal one. The latter is fixed by the mSTI for ties can be based on the dynamical generation of a non- the longitudinal gluon propagator. Hence, the only in- a a vanishing transverse background, Fµν Fµν = 0 , in the formation needed from the longitudinal system is the ini- infrared. This gluon condensate is the Savvidi6 vacuum tial value for the transverse gluon mass parameter (16). [78], and its generation in the present approach has been Note also that there is at least one value for the initial a a 4 discussed in [79] with Fµν Fµν 1 GeV . Then, a ver- gluon mass parameter that yields a valid confining solu- tex expansion about this non-trivial≈ IR-solution of the tion. In the following we vary the gluon mass parameter equation of motion introduces an IR-splitting of trans- and discuss the properties of the ensuing solutions. We verse and longitudinal vertices due to the transversality find a confining branch with both scaling and decoupling of the background field. This IR-splitting automatically solutions. In addition, we observe a transition to the implies irregularities as discussed in Sec. III A, and is deconfined Higgs-type branch. No Coloumb branch is 8
3.5 4 FRG, scaling FRG, scaling FRG, decoupling FRG, decoupling 3
2 3 1
0 three-gluon vertex dressing ghost-gluon vertex dressing
2.5 -1 0.1 1 10 0.1 1 10 p [GeV] p [GeV]