BRST invariant d = condensates in Gribov-Zwanziger theory D. Dudal,, C. P. Felix, L. Palhares, F. Rondeau,, D. Vercauteren + KU Leuven Campus Kortrijk Kulak, Department of Physics, Etienne Sabbelaan 53 bus 7657, 8500 Kortrijk, Belgium. Ghent University, Department of Physics and Astronomy, Krijgslaan 8-S9, 9000 Gent, Belgium Instituto de Física Teórica, Rua São Francisco Xavier 54, 0550-03, Maracanã, Rio de Janeiro, Brasil Ecole Normale Supérieure Paris-Saclay, Avenue du Président Wilson 6, 9435 Cachan Cedex, France + Duy Tân University, Institute of Research and Development, P809, 3 Quang Trung, Hải Châu, Đà Nẵng, Vietnam arxiv:8.54v [hep-ph] 9 Nov 08 In this proceeding, SU(N Yang-Mills theory is quantized in the linear covariant gauges, while taking into account the issue of Gribov copies and we construct the one-loop effective potential for a set of mass dimension condensates, including the Gribov parameter, that refines the infrared region of the Gribov-Zwanziger theory, whilst respecting renormalization group invariance and BRST symmetry. PACS numbers: I. INTRODUCTION The analytical explanation of quark and gluon confinement has been a big challenge in recent decades. At low temperatures (or in the infrared regime, IR, where confinement happens, the coupling constant g is large, hence the perturbative formalism can not be used in this regime. This research main focus has been on the gluon propagator and also ghost propagator in the IR. The gluon propagator is suppressed to a nonvanishing value at zero momentum violating reflection positivity and the ghost propagator is not enhanced at large volume according to the lattice data []. A possible analytical explanation for this behavior is obtained through adding dimension condensates to the Gribov-Zwanziger (GZ formalism, yielding to so-called Refined Gribov-Zwanziger (RGZ framework [, 3] that fits the lattice data quite well []. In this present study, the analysis of a non-trivial minimum of the effective action, which leads us to a dynamical transformation of the GZ action into the RGZ action, has been done in the presence of the A h A h and ϕϕ condensates at one-loop following earlier steps of [4], suitably generalized to respect BRST invariance following recent developments by some of us in the field. In this proceeding, we highlight a few steps, a comprehensive paper will be presented elsewhere. II. THE GRIBOV-ZWANZIGER ACTION IN THE LINEAR COVARIANT GAUGE In the IR region, the Gribov copies appear. Since the coupling constant g is large, these copies can not be eliminated [5]. A way to work around this problem is to restrict the functional integral to a specific region Ω in field space, a solution proposed by Gribov using the Landau gauge [5]. Moreover, this solution given by Gribov can be generalized to linear covariant gauge [6]: Ω = {A a µ; µ A a µ = iαb a, M ab (A h = µ D ab µ (A h > 0}. ( whereby the Hermitian Faddeev-Popov-related operator, M ab (A h = δ ab + gf abc (A h c µ µ, is positive. In (, A h µ is a non-local power series in the gauge field, gotten from the minimization of the functional f A [u] along the gauge Work presented at XIII Quark Confinement and the Hadron Spectrum - Confinement08, 3 July - 6 August 08, Maynooth University, Ireland Electronic address: david.dudal@kuleuven.be,caroline.felix@kuleuven.be, leticiapalhares@gmail.com, francois.rondeau@ens-parissaclay.fr, vercauterendavid@dtu.edu.vn
orbit of A µ [7 9], f A [u] min Tr d d x A u µa u µ, {u} A u µ = u A µ u + i g u µ u. ( A local minimum is found and given by ( A h µ = δ µν µ ν φ ν, µ A h µ = 0, [ ] φ ν = A ν ig A, A ν + ig [ ] A, ν A + O(A 3. (3 Here we highlight that A h µ is gauge invariant order by order [6]. The field A h µ can be localized by introducing an auxiliary Stueckelberg field ξ a [6, 0], while A h µ = (A h a µt a = h A a µt a h + i g h µ h, (4 h = e ig ξa T a, (5 Now, the local gauge invariance of A h µ under a gauge transformation u SU(N can be obtained from h u h, h h u, A µ u A µ u + i g u µ u. (6 Now, considering the BRST invariance, the Gribov-Zwanziger action in the linear covariant gauges, the total action is given by S = S YM + S GF + S GZ + S ε, (7 whereby S YM is the Yang-Mills action, S YM = d 4 xf a 4 µνf a µν, (8 S GF is the Faddeev-Popov gauge-fixing in linear covariant gauges, ( α S GF = d 4 x ba b a + ib a µ A a µ + c a µ D ab µ (Ac b, (9 α being the gauge parameter and α = 0 in Landau gauge; S GZ is the Gribov-Zwanziger action in its local form, S GZ = d 4 x [ ϕ ac µ ν D ab ν (A h ϕ bc µ ω ac µ ν (D ab ν (A h ω bc µ ] [ γ g d 4 x f abc (A h a µϕ bc µ + f abc (A h a µ ϕbc µ + d ] g (N c γ, (0 with ( ϕ ac µ, ϕ ac µ a pair of complex-conjugate bosonic fields, ( ω ac µ, ω ac µ a pair of anti-commuting complex-conjugate fields; and γ the Gribov parameter which is dynamically fixed by a gap equation that gives us the horizon function [, ], which can also be rewritten as [4] f abc (A h a µ(ϕ bc µ + ϕ bc µ = d(n γ g, ( Γ = 0, ( γ
3 whereby Γ is the quantum action defined by The last term from (7, e Γ = [dφ]e S. (3 S ε = d 4 x ε a µ (A h a µ (4 ensures, through the Lagrange multiplier ε, the transversality of the composite operator (A h a µ, µ (A h a µ = 0. The action S, (7, enjoys an exact BRST invariance, ss = 0 and s = 0 [6] sa a µ = D ab µ c b, sc a = g fabc c b c c, s c a = ib a, sb a = 0, sϕ ab µ = 0, sω ab µ = 0, s ω ab µ = 0, s ϕ ab µ = 0, sε a = 0, s(a h a µ = 0, sh ij = igc a (T a ik h kj. (5 III. REFINED GRIBOV-ZWANZIGER ACTION The BRST invariant d = condensates, A a h,µ Aa h,µ and ϕab µ, cause non-perturbative dynamical instabilities disturbing the Gribov-Zwanziger formalism [ 4]. ϕ ab µ ϕ ab µ guarantees that the gluon propagator is non-vanishing at zero momentum, and A a µa a µ is crucial to fit the lattice data [4]. The refined Gribov-Zwanziger action (RGZ is obtained adding these condensates to the GZ action via the local composite operator (LCO formalism, see [4]. The operators A h A h and ϕϕ will be added to the action via two BRST invariant bosonic sources τ and Q, µ ϕ ab sτ = 0 and sq = 0. (6 From here, we opted for the Landau gauge A = 0 for convenience, so that we can work with A h = A, as formally proven in [6]. Then, the action with these operators is written as Σ = S + S A + S ϕ ϕ + S vac, (7 whereby S is given by (7 and we also have S A = d d x τ Aa µa a µ, S ϕϕ = d d xq ϕ ac µ ϕ ac µ, ( ζ S vac = d d x τ + αqq + χqτ. (8 The parameters α, χ and ζ are the LCO parameters which guarantee that the divergences of the kind A (xa (y, etc. can be properly dealt with, see [4]. x y IV. THE EFFECTIVE ACTION CALCULUS In order to get the effective action, we have written the energy functional as e W(J = [dφ]e Σ, (9
4 where Σ is given by (7. The action Σ, (7, has three terms quadratic in the sources and they should be removed to facilitate calculations and interpret our results in terms of the vacuum energy. The easiest way to remove these terms is by introducing two auxiliary fields σ and σ via two identities written in terms of renormalized fields = [Dσ ] e d d x(σ + ā A + bq+ cτ, = [Dσ ] e + d d x(σ + dϕϕ+ēq+ f A, (0 multiplying the integral in (9. If we choose in dimensional regularization (d = 4 ɛ/ Z A ā = Zζζ ζ µɛ/ b = Z QQZ χχ χ Zζζ ζ µ ɛ/, c = Z ττ Zζζ ζµ ɛ/, Z ϕ d = µ ɛ/, ( Z χχ χ Z ζζ ζ Z ααα Z χχχ ē = Z QQ Z ζζ ζ Z ααα µ ɛ/, Z A f = Z τqz ζζ ζ Z QQ Z χχ χ µ Zζζ ζ Z QQ Z ɛ/, χχχ Z αα Z ζζ ζα we can remove the quadratic terms in sources. In the MS scheme and at one loop, the Z factors are given by [4] Z A = + 3 Ng 3 6π ɛ, Zζ = Z ζζ Z ττ = 3 Ng 3 6π ɛ, Z ζζ = + Ng 3 6π ɛ, Z g = Ng 3 6π ɛ, Z τ = 35 Ng 6 6π ɛ, Z QQ = Z g Z / A = 3 Ng 6π ɛ, Z χχ =, Z τq = 0, Zα = Z αα Z QQ = + 35 Ng 6 6π ɛ, Z αα = + 53 Ng 6 6π ɛ, Z ϕ = Z ϕ = Z g Z / A = + 3 Ng 6π ɛ. ( Therefore, (9 becomes where e W(Q,τ = [DΦ][Dσ,3 ] exp [ S GZ ( d d x cσ τ + σ 3 Q (( ē (σ 3 bσ σ 3 + ā f b ē ( b ē σ σ + f ē σ 3 d ] ē ( b σ σ 3 ϕϕ, (3 σ 3 = σ b σ ē. (4 So far, all LCO parameters, sources and fields have been renormalized, except the auxiliary fields σ s. Analyzing the term cσ τ = Z τ Zζ ζµ ɛ/ σ τ that appears in Eq.(3, it is easy to see that this field is indeed infinite, thus it also must be renormalized. As the Z-factors are infinite, if σ would be finite, then the quantity multiplying τ would be infinite. This would not make sense as a physical (and thus finite local operator O. The original σ field should then be infinite in order to get a finite quantity multiplying the finite source τ. It is then natural to define a renormalized A
finite field σ by σ Z ττ Zζζ σ Z ζ σ. Regarding σ 3, the term σ 3 Q in Eq.(3 teaches us, knowing that Q is finite, that σ 3 is already finite and thus should not be renormalized. In terms of the finite fields σ and σ 3, the energy functional now reads : e W(Q,τ S = [DΦ][Dσ,3 ] exp S GZ d d x ζσ τ + σ 3 Q ( b σ ē Z ζ ( σ ē 3 b σ σ 3 ā Z f b σ + f ē ē ζ Z σ 3 A ζ d σ b σ 3 ϕϕ. (5 ē Z ζ In this expression, all LCO parameters, sources and fields are now finite, and infinities are only present in the renormalization factors Z s, explicitly written or present in the coefficients ā,..., f. At one loop, χ = 0, Z τq = 0 [4] which implies that b = f = 0, then σ 3 = ēσ. Now, by analysis of the term σ 3 Q = ēσ Q in (5, a similar reasoning as above for σ shows that σ is infinite and should be renormalized defining a new finite field σ through σ Z QQ Zαα σ Z α σ. Hence, the energy functional in terms of the finite fields σ and σ and with one-loop coefficients reads : e W(Q,τ = [DΦ][Dσ, ] exp S GZ d d x σ σ + ā σ A d σ ϕϕ Z ζ Z α Z ζ Z α ] ζσ τ + ασ Q 5 (6 In this expression, infinities are now localized, only in the renormalization factors Z ζ, Z α, and in those hidden in ā and d. In order to have an expression of the form m A M ϕϕ, we defined the effective masses, m and M, respectively linked to AA and ϕϕ by: m Z ā ( σ = + 7 6 ζ Ng 3Ng 6π ɛ 9(N σ + O(g4, (7 M d ( σ = 35 6 Z α Ng 35Ng 6π ɛ 48(N σ + O(g4. (8 where the last equalities follow from considering the first order term of the Z-factors in MS scheme and α = α 0 = g 4(N and ζ = ζ 35Ng 0 = 9(N [4]. g 3Ng The ghost fields c, c, ω, ω give us just an overall factor. Now, to integrate over the ϕ and ϕ fields, we proceed as follows ϕ ab µ = U ab µ + iv ab µ, ϕ ab µ = U ab µ iv ab µ. (9 Then [DU, V] e d d x[vµ ab Pµν ac,bd Vν cd +U ab µ Pµν ac,bd U cd ν gγ f abc A a µu bc µ ] = = det Pµν ac,bd e d d xng γ 4 A a µpµν δ ab A b ν, (30
6 with Therefore, the first contribution Γ a to the effective potential is obtained by : P ac,bd µν ( M δ ac δ bd δ µν. (3 ΩΓ a = ln det P ac,bd µν = Tr ln P ac,bd µν, (3 resulting in [ Γ a = (N M 4 ] ɛ 4π + M4 M ln 8π µ M4 8π ( = 35Ng σ ( ɛ 48 4π + ln 35Ng σ 48(N µ + O(g 4. (33 The second contribution Γ b to the effective potential comes from the gluon field A µ. The quadratic part of the action containing A µ is e d d xa a µr ab µνa b ν (34 where [( R ab µν δ ab + m Nγ4 g ( ] α M δ µν µ ν. (35 Therefore, ΩΓ b = ln det Rac,bd µν = Tr ln Rac,bd µν, (36 resulting in ( Γ b = (N 3 (m (4π ɛ +5 4 4γ 4 g N + 3(N 4 4(4π ( ] M M 4 ln µ, [ x ln ( ( x µ + x x ln µ (37 where x and x are the solutions of the equation x + (M + m x + M m + λ 4 = 0, x = ( m + M + (m M 4λ 4, x = ( m + M (m M + 4λ 4. (38 The third part of effective potential Γ c is from the Gribov-Zwanziger action, Knowing that Z γ = Z / g Z /4 A, we get Γ c = dγ 4 0(N. (39 hence γ 4 0 = Z γ γ4, with Z γ = + 3 g N 6π ɛ, (40 ( Γ c = (N γ 4 4 + 3Ng 3π 3Ng 8π + O(g 4. (4 ɛ
And the last contribution comes from σ and σ : Γ d = ( σ Zζ σ Z α = σ σ + 3 6 Ng 6π ɛ σ + 35 Ng 6π ɛ σ + O(g 4. (4 The full effective potential given by Γ = Γ a + Γ b + Γ c + Γ d is finite when ɛ 0 at first order in g. Therefore, it can be written as Γ(m, M, λ 4 = 9(N m 4 3Ng 48(N M 4 35Ng λ4 (N Ng N 6π { λ 4 + 5 8 m4 + (N M 4 + 3 N 8 6π + 3 N 8 6π ( (N 3 ( M ln 4 µ {m 4 + M 4 λ 4 + ( m + M } (m M 4λ 4 [ ( ] ln µ m + M + (m M 4λ 4 {m 4 + M 4 λ 4 ( m + M } (m M 4λ 4 M 4 } [ ( ] ln µ m + M (m M 4λ 4. (43 7 with m 3Ng = 9(N σ, (44 M 35Ng = 48(N σ. (45 and λ 4 Ng γ 4. The next step is to analyze the gap equations given by Γ M = 0, Γ m = 0, Γ = 0. (46 λ4 Unfortunately, an acceptable result was not obtained in this particular scheme. The resolution was to not fix a scheme, we rather rewrote the effective potential in a general scheme as has been done in e.g. [3]. Details of this procedure will be published elsewhere. The effective potential (43, in general scheme, becomes Γ gen (m, M, λ 4, b 0 = 9(N 6Ng m4 4(N 35Ng M 4 (N M 4 ( ( M 6π ln µ λ 4 N Ng λ 4 N 6π (b 0 + 3 N { 5 4 6π 4 (m4 + M 4 λ 4 m + M λ 4 [ m M + λ 4 ] ln µ 4 [ ] [ ] } +(m + M 4λ4 (m 4λ 4 (m M arctan M M m + M + ln µ M 4, (47 b 0 being a parameter related to the chosen scheme for the coupling. It was fixed, at the end, by matching our values for the complex conjugate poles masses of the transverse gluon propagator to those estimated from lattice data [] when the gap equations are solved for. The effective masses m and M and the Gribov parameter γ were obtained as functions of the parameter b 0 and µ in units Λ =, N = 3 and also for N =. Notice that these poles masses are gauge and scheme independent [6], so we benefitted from this to fix the parameters b 0 and µ by using a minimal
8 external lattice input to determine the optimum scheme. We got b 0 = 3.643 and µ =.49. With this procedure, we obtained a reasonable value for the coupling constant, namely 0.38. Therefore, in this case, the perturbative result is relatively trustworthy. The Gribov parameter γ is 0.637 and the vacuum energy is 6.955. The Hessian determinant is positive and also the second derivatives, Γ gen Γ gen =.668 Γ gen = 0.6 solved solved M m = 0.0. (48 solved M m Then, the solution does correspond to a minimum. The future step will be to extend this research to finite temperatures and to study if the deconfinement transitions reflects itself in a change in the propagator behavior and to check if, with the Polyakov loop added to the game, we observe the transition also in that order parameter. [] D. Dudal, O. Oliveira and P. J. Silva, Annals Phys. 397 (08 35; A. Cucchieri, D. Dudal, T. Mendes, O. Oliveira, M. Roelfs and P. J. Silva, arxiv:809.084 [hep-lat]; A. Cucchieri, T. Mendes and E. M. S. Santos, Phys. Rev. Lett. 03 (009 460; P. Bicudo, D. Binosi, N. Cardoso, O. Oliveira and P. J. Silva, Phys. Rev. D 9 (05 no., 454; A. Cucchieri, D. Dudal, T. Mendes and N. Vandersickel, Phys. Rev. D 85 (0 09453. [] D. Dudal, S. P. Sorella, N. Vandersickel and H. Verschelde, Phys. Rev. D 77 (008 0750. [3] D. Dudal, J. A. Gracey, S. P. Sorella, N. Vandersickel and H. Verschelde, Phys. Rev. D 78 (008 065047. [4] D. Dudal, S. P. Sorella and N. Vandersickel, Phys. Rev. D 84 (0 065039. [5] V. N. Gribov, Nucl. Phys. B 39 (978. [6] M. A. L. Capri et al., Phys. Rev. D 9 (05 045039; M. A. L. Capri et al., Phys. Rev. D 93 (06 06509; M. A. L. Capri et al., Phys. Rev. D 94 (06 05035; M. A. L. Capri, D. Fiorentini, A. D. Pereira, R. F. Sobreiro, S. P. Sorella and R. C. Terin, Annals Phys. 376 (07 40; M. A. L. Capri, D. Dudal, A. D. Pereira, D. Fiorentini, M. S. Guimaraes, B. W. Mintz, L. F. Palhares and S. P. Sorella, Phys. Rev. D 95 (07 0450. [7] G. Dell Antonio and D. Zwanziger, Nucl. Phys. B 36 (989 333. [8] P. van Baal, Nucl. Phys. B 369 (99 59. [9] M. Lavelle and D. McMullan, Phys. Rept. 79 (997. [0] M. A. L. Capri, D. Fiorentini, M. S. Guimaraes, B. W. Mintz, L. F. Palhares and S. P. Sorella, Phys. Rev. D 94 (06 065009. [] D. Zwanziger, Nucl. Phys. B 33 (989 53. [] D. Zwanziger, Nucl. Phys. B 399 (993 477. [3] D. Dudal, R. F. Sobreiro, S. P. Sorella and H. Verschelde, Phys. Rev. D 7 (005 0406.