Effect of the Gribov horizon on the Polyakov loop and vice versa

We consider finite temperature SU(2) gauge theory in the continuum formulation, which necessitates the choice of a gauge fixing. Choosing the Landau gauge, the existing gauge copies are taken into account by means of the Gribov-Zwanziger (GZ) quantization scheme, which entails the introduction of a dynamical mass scale (Gribov mass) directly influencing the Green functions of the theory. Here, we determine simultaneously the Polyakov loop (vacuum expectation value) and Gribov mass in terms of temperature, by minimizing the vacuum energy w.r.t. the Polyakov loop parameter and solving the Gribov gap equation. Inspired by the Casimir energy-style of computation, we illustrate the usage of Zeta function regularization in finite temperature calculations. Our main result is that the Gribov mass directly feels the deconfinement transition, visible from a cusp occurring at the same temperature where the Polyakov loop becomes nonzero. In this exploratory work we mainly restrict ourselves to the original Gribov-Zwanziger quantization procedure in order to illustrate the approach and the potential direct link between the vacuum structure of the theory (dynamical mass scales) and (de)confinement. We also present a first look at the critical temperature obtained from the Refined Gribov-Zwanziger approach. Finally, a particular problem for the pressure at low temperatures is reported.


Introduction
Within SU (N ) Yang-Mills gauge theories, it is well accepted that the asymptotic particle spectrum does not contain the elementary excitations of quarks and gluons.These color charged objects are confined into color neutral bound states: this is the so-called color confinement phenomenon.It is widely believed that confinement arises due to non-perturbative infrared effects.Many criteria for confinement have been proposed (see the nice pedagogical introduction [1]).A very natural observation is that gluons (due to the fact that they are not observed experimentally) should not belong to the physical spectrum in a confining theory.On the other hand, the perturbative gluon propagator satisfies the criterion to belong to the physical spectrum (namely, it has a Källén-Lehmann spectral representation with positive spectral density).Hence, non-perturbative effects must dress the perturbative propagator in such a way that the positivity conditions are violated, such that it does not belong to the physical spectrum anymore.A well-known criterion is related to the fact that the Polyakov loop [2] is an order parameter for the confinement/deconfinement phase transition via its connection to the free energy of a (very heavy) quark.The importance to clarify the interplay between these two different points of view (nonperturbative Green function's behaviour vs. Polyakov loop) can be understood by observing that while there are, in principle, infinitely many different ways to write down a gluon propagator which violates the positivity conditions, it is very likely that only few of these ways turns out to be compatible with the Polyakov criterion.
One of the most fascinating non-perturbative infrared effects is related to the appearance of Gribov copies [3] which represent an intrinsic overcounting of the gauge-field configurations which the perturbative gauge-fixing procedure is unable to take care of.Soon after Gribov's seminal paper, Singer showed that any true gauge condition, as the Landau gauge1 , presents this obstruction [4] (see also [5]).The presence of Gribov copies close to the identity induces the existence of non-trivial zero modes of the Faddeev-Popov operator, which make the path integral ill defined.Even when perturbation theory around the vacuum is not affected by Gribov copies close to the identity (in particular, when YM-theory is defined over a flat space-time2 with trivial topology [10]), Gribov copies have to be taken into account when considering more general cases (such as with toroidal boundary conditions on flat space-time [11,12]).Thus, in the following only the standard boundary conditions will be considered.
The most effective method to eliminate Gribov copies, at leading order proposed by Gribov himself, and refined later on by Zwanziger [3,13,14,15] corresponds to restricting the path integral to the so-called Gribov region, which is the region in the functional space of gauge potentials over which the Faddeev-Popov operator is positive definite.The Faddeev-Popov operator is Hermitian in the Landau gauge, so it makes sense to discuss its sign.In [16,17] Dell'Antonio and Zwanziger showed that all the orbits of the theory intersect the Gribov region, indicating that no physical information is lost when implementing this restriction.Even though this region still contains copies which are not close to the identity [18], this restriction has remarkable effects.In fact, due to the presence of a dynamical (Gribov) mass scale, the gluon propagator is suppressed while the ghost propagator is enhanced in the infrared.More general, an approach in which the gluon propagator is "dressed" by non-perturbative corrections which push the gluon out of the physical spectrum leads to propagators and glueball masses in agreement with the lattice data [19,20].With the same approach, one can also solve the old problem of the Casimir energy in the MIT-bag model [21].Moreover, the extension of the Gribov gap equation at finite temperature provides one with a good qualitative understanding, already within the semiclassical approximation, of the deconfinement temperature as well as of a possible intermediate phase in which features of the confining phase coexist with features of the fully deconfined phase in agreement with different approaches (see [22] and references therein).Furthermore, within this framework the presence of the Higgs field [23,24] as well of a Chern-Simons term in 2+1 dimensions [25] can be accounted for as well.
For all these reasons, it makes sense to compute the vacuum expectation value of the Polyakov loop when we eliminate the Gribov copies using the Gribov-Zwanziger (GZ) approach.Related computations are available using different techniques to cope with nonperturbative propagators at finite temperature, see e.g.[26,27,28,29,30,31,32,33,34,35,36].In the present paper, we will perform for the first time (to the best of the authors' knowledge) this computation, using two different techniques, to the leading oneloop approximation.In [37,38,39], it was already pointed out that the Gribov-Zwanziger quantization offers an interesting way to illuminate some of the typical infrared problems for finite temperature gauge theories.
In Section 2, we provide a brief technical overview of the Gribov-Zwanziger quantization process and eventual effective action.In the following Section 3, the Polyakov loop is introduced into the GZ theory via the background field method, building on work of other people [27,28,32].Next, Section 4 handles the technical computation of the leading order finite temperature effective action, while in Section 5 we discuss the gap equations, leading to our estimates for both Polyakov loop and Gribov mass.The key finding is a deconfinement phase transition at the same temperature at which the Gribov mass develops a cuspy behaviour.We subsequently also discuss the pressure and energy anomaly.Due to a problem with the pressure in the GZ formalism (regions of negativity), we take a preliminary look at the situation upon invoking the more recently developed Refined Gribov-Zwanziger approach.We summarize in Section 7.

Mills theories
Let us start by giving a short overview of the Gribov-Zwanziger framework [3,13,14,15].As already mentioned in the Introduction, the Gribov-Zwanziger action arises from the restriction of the domain of integration in the Euclidean functional integral to the Gribov region Ω, which is defined as the set of all gauge field configurations fulfilling the Landau gauge, ∂ µ A a µ = 0, and for which the Faddeev-Popov operator The boundary ∂Ω of the region Ω is the (first) Gribov horizon.
One starts with the Faddeev-Popov action in the Landau gauge where S YM and S gf denote, respectively, the Yang-Mills and the gauge-fixing terms, namely and where (c a , c a ) stand for the Faddeev-Popov ghosts, b a is the Lagrange multiplier implementing the Landau gauge, D ab µ = (δ ab ∂ µ − gf abc A c µ ) is the covariant derivative in the adjoint representation of SU (N ), and F a µν denotes the field strength: Following [3,13,14,15], the restriction of the domain of integration in the path integral is achieved by adding to the Faddeev-Popov action S FP an additional term H(A), called the horizon term, given by the following non-local expression where M −1 stands for the inverse of the Faddeev-Popov operator.The partition function can then be written as [3,13,14,15]: where V is the Euclidean space-time volume.The parameter γ has the dimension of a mass and is known as the Gribov parameter.It is not a free parameter of the theory.It is a dynamical quantity, being determined in a self-consistent way through a gap equation called the horizon condition [3,13,14,15], given by where the notation H(A, γ) GZ means that the vacuum expectation value of the horizon function H(A, γ) has to be evaluated with the measure defined in Eq.( 6).An equivalent all-order proof of eq.( 7) can be given within the original Gribov no-pole condition framework [3], by looking at the exact ghost propagator in an external gauge field [40].
Although the horizon term H(A, γ), eq.( 5), is non-local, it can be cast in local form by means of the introduction of a set of auxiliary fields (ω ab µ , ω ab µ , φab µ , ϕ ab µ ), where ( φab µ , ϕ ab µ ) are a pair of Bosonic fields, while (ω ab µ , ω ab µ ) are anti-commuting.It is not difficult to show that the partition function Z GZ in eq.( 6) can be rewritten as [13,14,15] where Φ accounts for the quantizing fields, A, c, c, b, ω, ω, φ, and ϕ, while S GZ [Φ] is the Yang-Mills action plus gauge fixing and Gribov-Zwanziger terms, in its localized version, with and It can be seen from ( 6) that the horizon condition (7) takes the simpler form which is called the gap equation.The quantity E v (γ) is the vacuum energy defined by The local action S GZ in eq.( 9) is known as the Gribov-Zwanziger action.Remarkably, it has been shown to be renormalizable to all orders [13,14,15,41,42,43,44,45].This important property of the Gribov-Zwanziger action is a consequence of an extenstive set of Ward identities constraining the quantum corrections in general and possible divergences in particular.In fact, introducing the nilpotent BRST transformations it can immediately be checked that the Gribov-Zwanziger action exhibits a soft breaking of the BRST symmetry, as summarized by the equation where Notice that the breaking term ∆ is of dimension two in the fields.As such, it is a soft breaking and the ultraviolet divergences can be controlled at the quantum level.The properties of the soft breaking of the BRST symmetry of the Gribov-Zwanziger theory and its relation with confinement have been object of intensive investigation in recent years, see [46,47,48,49,50,51,52,53,54,55].Here, it suffices to mention that the broken identity ( 15) is connected with the restriction to the Gribov region Ω.However, a set of BRST invariant composite operators whose correlation functions exhibit the Källén-Lehmann spectral representation with positive spectral densities can be consistently introduced [56].These correlation functions can be employed to obtain mass estimates on the spectrum of the glueballs [19,20].
Let us conclude this brief review of the Gribov-Zwanziger action by noticing that the terms S gf and S 0 in expression ( 9) can be rewritten in the form of a pure BRST variation, i.e.
so that from which eq.( 15) becomes apparent.

The Polyakov loop and the background field formalism
In this section we shall investigate the confinement/deconfinement phase transition of the SU (2) gauge field theory in the presence of two static sources of (heavy) quarks.The standard way to achieve this goal is by probing the Polyakov loop order parameter, with P denoting path ordering, needed in the non-Abelian case to ensure the gauge invariance of P.This path ordering is not relevant at one-loop order, which will considerably simplify the computations of the current work.In analytical studies of the phase transition involving the Polyakov loop, one usually imposes the so-called "Polyakov gauge" on the gauge field, in which case the time-component A 0 becomes diagonal and independent of (imaginary) time.This means that the gauge field belongs to the Cartan subalgebra.More details on Polyakov gauge can be found in [28,57,58].Besides the trivial simplification of the Polyakov loop, when imposing the Polyakov gauge it turns out that the quantity A 0 becomes a good alternative choice for the order parameter instead of P.This extra benefit can be proven by means of Jensen's inequality for convex functions and is carefully explained in [28], see also [27,29,30,31,32].For example, for the SU (2) case we have the following: if 1 2 gβ A 0 = π 2 then we are in the "unbroken symmetry phase" (confined or disordered phase), equivalent to P = 0; otherwise, if 1  2 gβ A 0 < π 2 , we are in the "broken symmetry phase" (deconfined or ordered phase), equivalent to P = 0. Since P ∝ e −F T with T the temperature and F the free energy of a heavy quark, it is clear that in the confinement phase, an infinite amount of energy would be required to actually get a free quark.The broken/restored symmetry referred to is the Z N center symmetry of a pure gauge theory (no dynamical matter in the fundamental representation).
A slightly alternative approach to access the Polyakov loop was worked out in [32].In order to probe the phase transition in a quantized non-Abelian gauge field theory, we use, following [32], the Background Field Gauge (BFG) formalism, detailed in general in e.g.[65].Within this framework, the effective gauge field will be defined as the sum of a classical field Āµ and a quantum field A µ : a µ (x) = a a µ (x)t a = Āµ +A µ , with t a the infinitesimal generators of the SU (N ) symmetry group.The BFG method is a convenient approach, since the tracking of breaking/restoration of the Z N symmetry becomes easier by choosing the Polyakov gauge for the background field.
Within this framework, it is convenient to define the gauge condition for the quantum field, known as the Landau-DeWitt (LDW) gauge fixing condition, where Dab µ = δ ab ∂ µ − gf abc Āc µ is the background covariant derivative.After integrating out the (gauge fixing) auxiliary field b a , we end up with the following Yang-Mills action, Notice that, concerning the quantum field A µ , the condition ( 20) is equivalent to the Landau gauge, yet the action still has background center symmetry.The LDW gauge is actually recovered in the limit ξ → 0, taken at the very end of each computation.
It is perhaps important here to stress that we are restricting our analysis to the (background) Landau gauge, for which a derivation argument in favor of the action ( 21) can be provided.For a vanishing background, this is precisely the original Gribov-Zwanziger construction [3,14,15], also applicable to the Coulomb gauge.More recently, it was also generalized to the SU (2) maximal Abelian gauge in [59].Intuitively, it might be clear that the precise influence on the quantum dynamics by Gribov copies can strongly depend on the chosen background, given that Gribov copies are defined via the zero modes of the Faddeev-Popov operator of the chosen gauge condition, which itself explicitly depends on the chosen background.This is open to further research, as it has not been pursued in the literature yet.Though, for a constant background as relevant for the current purposes, it will be discussed elsewhere that the action is indeed obtainable via a suitable extension of the arguments of [3,14,15].
In the absence of a background, a proposal for a generalization to the linear covariant gauges was put forward in [60,61], albeit leading to a very complicated nonlocal Lagrangian structure, containing e.g.reciprocals and exponentials of fields.To our knowledge, no practical computations were done so far with this formalism.Nonetheless, potential problems with gauge parameter dependence of physical quantities were discussed in [60,61], not surprisingly linked to the softly broken BRST symmetry, see also our Section 2 for more on this and relevant references.
A very recent alternative for the linear covariant gauges was worked out in [62], partially building on earlier work of [63].With this proposal, it was explicitly checked at one loop that the Gribov parameter γ 2 and vacuum energy are gauge parameter independent.This at least suggests that in this class of covariant gauges, an approach to Gribov copies can be worked out that is compatible with gauge parameter independence [64].
As explained for the simple Landau gauge in the previous section, the Landau background gauge condition is also plagued by Gribov ambiguities, and the Gribov-Zwanziger procedure is applicable also in this instance.The starting point of our analysis is, therefore, the GZ action modified for the BFG framework (see [66]): As mentioned before, with the Polyakov gauge imposed to the background field Āµ , the time-component becomes diagonal and time-independent.In other words, we have Āµ (x) = Ā0 δ µ0 , with Ā0 belonging to the Cartan subalgebra of the gauge group.For instance, in the Cartan subalgebra of SU (2) only the t 3 generator is present, so that Āa 0 = δ a3 Ā3 0 ≡ δ a3 Ā0 .As explained in [32], at leading order we then simply find, using the properties of the Pauli matrices, where we defined r = gβ Ā0 , with β the inverse temperature.Just like before, r = π corresponds to the confinement phase, while 0 ≤ r < π corresponds to deconfinement.With a slight abuse of language, we will refer to the quantity r as the Polyakov loop hereafter.
Since the scope of this work is limited to one-loop order, only terms quadratic in the quantum fields in the action ( 22) shall be considered.One then immediately gets an action that can be split in term coming from the two color sectors: the 3rd color direction, called Cartan direction, which does not depend on the parameter r; and one coming from the 2 × 2 block given by the 1st and 2nd color directions.This second 2 × 2 color sector is orthogonal to the Cartan direction and does depend on r.The scenario can then be seen as a system where the vector field has an imaginary chemical potential irT and has isospins +1 and −1 related to the 2 × 2 color sector and one isospin 0 related to the 1 × 1 color sector. 4 The finite temperature effective action at leading order Considering only the quadratic terms of ( 22), the integration of the partition function gives us the following vacuum energy at one-loop order, defined according to (13), where V is now just the spacial volume.Here, D is the covariant derivative in the adjoint representation in the presence of the background A3 0 field and λ 4 = 2N g 2 γ 4 .Throughout this work, it is always tacitly assumed we are working with N = 2 colors, although we will frequently continue to explicitly write N dependence for generality.Using the usual Matsubara formalism, we have that D 2 = (2πnT + rsT ) 2 + q 2 , where n is the Matsubara mode, q is the spacelike momentum component, and s is the isospin, given by −1, 0, or +1 for the SU (2) case 3 .
The general trace is of the form which will be computed immediately below.

The sum-integral: 2 different computations
We want to compute the following expression: One way to proceed is to start by deriving the previous expression with respect to m 2 .Then, one can use the well-known formula from complex analysis where the sum is over the poles z 0 of the function f (z).Subsequently we integrate with respect to m 2 (and determine the integration constant by matching the result with the known T = 0 case).Finally one can split off the analogous T = 0 trace (which does not depend on the background field) to find where the limit → 0 was taken in the (convergent) second integral.The first term in the r.h.s. is the (divergent) zero temperature contribution.
Another way to compute the above integral is by making use of Zeta function regularization techniques, which are particularly useful in the computation of the Casimir energy in various configurations see [68,69].The advantage of this second technique is that, although it is less direct, it provides one with an easy way to analyze the high and low temperature limits as well as the small mass limit, as we will now show.Moreover, within this framework, the regularization procedures are often quite transparent.One starts by writing the logarithm as ln x = − lim s→0 ∂ s x −s , after which the integral over the momenta can be performed: where the renormalization scale µ has been introduced to get dimensional agreement for s = 0, and where we already put = 0, as s will function as a regulator -i.e.we assume s > 3/2 and analytically continuate to bring s → 0. Using the integral representation of the Gamma function, the previous expression can be recast to where the variable of integration was transformed as y = 4π 2 T 2 t ≥ 0 in the second line.Using the Poisson rule (valid for positive ω): we obtain that where K ν (z) is the modified Bessel function of the second kind.Simplifying this, we find where the first term is the T = 0 contribution, and the sum is the finite-temperature correction.Using numerical integration and series summation, it can be checked that both results ( 29) and ( 34) are indeed identical.Throughout this paper, we will mostly base ourselves on the expression (29).Nonetheless the Bessel series is quite useful in obtaining the limit cases m = 0, T → ∞, and T → 0 by means of the corresponding behaviour of K 2 (z).Observing that we obtain where Li n s is the polylogarithm or Jonquière's function. Analogously, so that Li 2 e −irs + Li 2 e irs − 4T 2 Li 4 e −irs + Li 4 e irs .
(37) Finally for T → 0 we can use the asymptotic expansion of the Bessel function [70]: where a k (ν) are finite factors.So, at first order (k = 0),

The result for further usage
Making use of the result (29) we may define so that the vacuum energy ( 25) can be rewritten as where tr T =0 denotes the trace taken at zero temperature.
5 Minimization of the effective action, the Polyakov loop and the Gribov mass

Warming-up exercise: assuming a T -independent Gribov mass λ
As a first simpler case, let us simplify matters slightly by assuming that the temperature does not influence the Gribov parameter λ.This means that λ will be supposed to assume its zero-temperature value, which we will call λ 0 , given by the solution of the gap equation ( 7) at zero temperature.In this case, only the terms with the function I matter in (41), since the other terms do not explicitely depend on the Polyakov line r.Plotting this part of the potential (see Figure 1), one finds by visual inspection that a second-order phase transition occurs from the minimum with r = π to a minimum with r = π.
The transition can be identified by the condition Using the fact that when s = ±1 and zero when s = 0, the equation ( 42) can be straightforwardly solved numerically for the critical temperature.We find

The T -dependence of the Gribov mass λ
Let us now investigate what happens to the Gribov parameter λ when the temperature is nonzero.Taking the derivative of the effective potential (41) with respect to λ 2 and dividing by d(N 2 − 1)λ 2 /N g 2 (as we are not interested in the solution λ 2 = 0) yields the gap equation for general number of colors N : where the notation ∂I/∂m 2 denotes the derivative of I with respect to its first argument (written m 2 in ( 40)).If we now define λ 0 to be the solution to the gap equation at T = 0: then we can subtract this equation from the general gap equation (45).After dividing through (d − 1)N g 2 /2d and setting d = 4 and N = 2, the result is where now all integrations are convergent.This equation can be easily solved numerically to yield λ as a function of temperature T and background r, in units λ 0 .This is shown in Figure 2.

Absolute minimum of the effective action
As λ does not change much when including its dependence on temperature and background, the transition is still second order and its temperature is, therefore, still given by the condition (42).Now, however, the potential depends explicitely on r, but also implicitely due to the presence of the r-dependent λ.We therefore have Now, dλ/dr| r=π = 0 due to the symmetry at that point.Furthermore, as we are considering λ = 0, ∂E v /∂λ = 0 is the gap equation and is solved by λ(r).Therefore, we find for the condition of the transition: where the derivative is taken with respect to the explicit r only.
We already solved equation ( 49) in section 5.1, giving (44): As we computed λ as a function of r and T in section 5.2 already, it is again straightforward to solve this equation to give the eventual critical temperature: as expected only slightly different from the simplified estimate (44) found before.

5.4
The T -dependence of the Polyakov loop r and the equation of state

Deconfinement transition and its imprint on the Gribov mass
Let us now investigate the temperature dependence of r.The physical value of the background field r is found by minimizing the vacuum energy: From the vacuum energy (41) we have The expression (53) was obtained after summation over the possible values of s.Furthermore, we used the fact that I(m 2 , r, +1, T ) = I(m 2 , r, −1, T ) and that s = 0 accounts for terms independent of r, which are cancelled by the derivation w.r.t.r.From (40) one can get, whenever s = ±1: Since ( 53) is finite, we can numerically obtain r as a function of temperature.From the dotted curve in Figure 3 one can easily see that, for T > T crit ≈ 0.40λ 0 , we have r = π, pointing to a deconfined phase, confirming the computations of the previous section.In the same figure, λ(T ) is plotted in a continuous line.We observe very clearly that the Gribov mass λ(T ) develops a cusp-like behaviour exactly at the critical temperature T = T crit .

Equation of state
Following [71], we can also extract an estimate for the (density) pressure p and the interaction measure I/T 4 , shown in Figure 4 (left and right respectively).As usual the (density) pressure is defined as which is related to the free energy by p = −E v .Here the plot of the pressure is given relative to the Stefan-Boltzmann limit pressure: p SB = κT 4 , where κ = (N 2 − 1)πT 4 /45 is the Stefan-Boltzmann constant accounting for all degrees of freedom of the system at high temperature.We subtract the zero-temperature value, such that the pressure becomes zero at zero temperature: Namely, after using the MS renormalization prescription and choosing the renormalization parameter μ so that the zero temperature gap equation is satisfied, we have the following expression for the pressure (in units of λ 4 0 ), In (57) prime quantities stand for quantities in units of λ 0 , while λ and λ 0 satisfy their gap equation.The last term of (57) accounts for the zero temperature subtraction, so that p(0) = 0, according to the definition of I(m 2 , r, T ) in (40).Note that the coupling constant does not explicitly appear in (57) and that λ 0 stands for the Gribov parameter at T = 0.
The interaction measure I is defined as the trace anomaly in units of T 4 , and I is nothing less than the trace of the of the stress-energy tensor, given by with being the internal energy density, which is defined as = E v + T s (with s the entropy density), u = (1, 0, 0, 0) and η µν the (Euclidean) metric of the space-time.Given the thermodynamic definitions of each quantity (energy, pressure and entropy), we obtain Both quantities display a behavior similar to that presented in [39] (but note that in they plot the temperature in units of the critical temperature (T c in their notation), while we use units λ 0 ).Besides this, and the fact that we included the effect of Polyakov loop on the Gribov parameter, in [39] a lattice-inspired effective coupling was introduced at finite temperature while we used the exact one-loop perturbative expression, which is consistent with the order of all the computations made here.
However, we notice that at temperatures relatively close to our T c , the pressure becomes negative.This is clearly an unphysical feature, possibly related to some missing essential physics.For higher temperatures, the situation is fine and the pressure moreover displays a similar to what is seen in lattice simulations for the nonperturbative pressure (see [72] for the SU (3) case).A similar problem is present in one of the plots presented in [39,Fig. 4], although no comment is made about it.Another strange feature is the oscillating behaviour of both pressure and interaction measure at low temperatures.Something similar was already observed in [73] where a quark model was employed with complex conjugate quark mass.It is well-known that the gluon propagator develops two complex conjugate masses in Gribov-Zwanziger quantization, see e.g.[19,20,56,77] for some more details, so we confirm the findings of [73] that, at least at leading order, the thermodynamic quantities develop an oscillatory behaviour.We expect this oscillatory behaviour would in principle also be present in [39] if the pressure and interaction energy were to be computed at lower temperatures than shown there.In any case, the presence of complex masses and their consequences gives us a warning that a certain care is needed when using GZ dynamics, also at the level of spectral properties as done in [74,75], see also [46,76].
These peculiarities justify giving an outline in the next Section of the behaviours of the pressure and interaction measure in an improved formalism, such as in the Refined Gribov-Zwanziger one.6 Outlook to the Refined Gribov-Zwanziger formalism The previous results can be slightly generalized to the case of the Refined Gribov-Zwanziger (RGZ) formalism studied in [42,43,45,79,80].In this refined case, additional nonperturbative vacuum condensates such as A 2 µ and φab µ ϕ ab µ are to be introduced.The corresponding mass dimension two operators get a nonzero vacuum expectation value (thereby further lowering the vacuum energy) and thus influence the form of the propagator and effective action computation.The predictions for the RGZ propagators, see also [81,82,83], are in fine agreement with ruling T = 0 lattice data, see e.g. also [84,85,86,87,88,89,90,91,92].This is in contrast with the original GZ predictions, such that it could happen that the finite temperature version of RGZ is also better suited to describe the phase transition and/or thermodynamical properties of the pure gauge theory.
Due to the more complex nature of the RGZ effective action (more vacuum condensates), we will relegate a detailed (variational) analysis of their finite temperature counterparts4 to future work, as this will require new tools.Here, we only wish to present a first estimate of the deconfinement critical temperature T c using as input the T = 0 RGZ gluon propagator where the nonperturbative mass parameters are fitted to lattice data for the same propagator.More precisely, we use [77] ∆ ab µν (p) = δ ab where we omitted the global normalization factor Z which drops out from our leading order computation 5 .In this expression, we have that The free energy associated to the RGZ framework can be obtained by following the same steps as in section 4, leading to with r 2 ± standing for minus the roots of the denominator of the gluon propagator (60), N 2 = M 2 + ρ 1 , and I(m 2 , r, T ) given by (40).Explicitly, the roots are The (central) condensate values were extracted from [77]: Once again the vacuum energy will be minimized with respect to the Polyakov loop expectation value r.For the analysis of thermodynamic quantities, only contributions coming from terms proportional to I(m 2 , r, T ) will be needed.Therefore, we will always consider the difference E v (T ) − E v (T = 0).Since in the present (RGZ) prescription the condensates are given by the zero temperature lattice results (64) instead of satisfying gap equations, the divergent contributions to the free energy are subtracted, and no specific choice of renormalization scheme is needed.Furthermore, explicit dependence on the coupling constant seems to drop out of the one-loop expression, such that no renormalization scale has to be chosen.Following the steps taken in Section 5.1, we find a second order phase transition at the temperature: which is not that far from the value of the SU (2) deconfinement temperature found on the lattice: T c ≈ 0.295 GeV, as quoted in [96,97].
In future work, it would in particular be interesting to find out whether -upon using the RGZ formalism-the Gribov mass and/or RGZ condensates directly feel the deconfinement transition, similar to the cusp we discovered in the Gribov parameter following the exploratory restricted analysis of this paper.This might also allow to shed further light on the ongoing discussion of whether the deconfinement transition should be felt at the level of the correlation functions, in particular the electric screening mass associated with the longitudinal gluon propagator [98,99,100,101].
Let us also consider the pressure and interaction measure once more.The results are shown in Figure 5 and Figure 6, respectively.The oscillating behaviour at low temperature persists at leading order, while a small region of negative pressure is still present -see the right plot of Figure 5.These findings are similar to [102] (low temperature results are not shown there), where two sets of finite temperature RGZ fits to the SU (3) lattice data were used [103,104], in contrast with our usage of zero temperature SU (2) data.In any case, more involved analysis of the RGZ finite temperature dynamics is needed to make firmer statements.As already mentioned before, there is also the possibility that important low temperature physics is missing, as for instance the proposal of [102] related to the possible effect of light electric glueballs near the deconfinement phase transition [105,106].Obviously, these effects are absent in the current treatment (or in most other treatments in fact).

Summary
In this paper we studied the Gribov-Zwanziger (GZ) action for SU (2) gauge theories with the Polyakov loop coupled to it via the background field formalism.Doing so, we were able to compute in a simultaneous fashion the finite temperature value of the Polyakov loop and Gribov mass to the leading one-loop approximation.The latter dynamical scale enters the theory as a result of the restriction of the domain of gauge field integration to avoid (infinitesimally connected) Gribov copies.Our main result is that we found clear evidence of a second order deconfinement phase transition at finite temperature, an occurrence accompanied by a cusp in the Gribov mass, which thus directly feels the transition.It is perhaps worthwhile to stress here that at temperatures above T c , the Gribov mass is nonzero, indicating that the gluon propagator still violates positivity and as such it rather describes a quasi-than a "free" observable particle, see also [26,107] for more on this.
We also presented the pressure and trace anomaly, indicating there is a problem at temperatures around the critical value when using the original GZ formulation.We ended with a first look at the changes a full-fledged analysis with the Refined Gribov-Zwanziger (RGZ) formalism might afflict, given that the latter provides an adequate description of zero temperature gauge dynamics, in contrast to the GZ predictions.This will be studied further in upcoming work.Note that, even not considering finite temperature corrections to the condensates in the RGZ formalism, the region of negative pressure is considerably smaller than the region found with the GZ formalism.
A further result of the present paper, which is interesting from the methodological point of view, is that it shows explicitly that finite-temperature computations (such as the computation of the vacuum expectation value of the Polyakov loop) are very suitable to be analyzed using analytical Casimir-like techniques.The interesting issue of Casimir-style computations at finite temperatures is that, although they can be more involved, they provide one with easy tools to analyze the high and low temperature limits as well as the small mass limit.Moreover, within the Casimir framework, the regularization procedures are often quite transparent.Indeed, in the present paper, we have shown that the computation of the vacuum expectation value of the Polyakov loop is very similar to the computation of the Casimir energy between two plates.We believe that this point of view can be useful in different contexts as well.

)Ev T Ev T 0 Λ 0 4 Figure 1 :
Figure1: The effective potential(41) at the temperatures (from below upwards at r = π) 0.42, 0.44, 0.46, and 0.48 times λ as a function of r, with the simplifying assumption that λ maintains its zero-temperature value λ 0 throughout.It can be seen that the minimum of the potential moves away from r = π in between T = 0.44 λ and 0.46 λ.

Figure 2 :
Figure 2: The Gribov parameter λ as a function of the temperature T at r equals to zero (upper line) and π (lower line), in units of the zero-temperature Gribov parameter λ 0 .

Figure 3 :
Figure 3: The dotted curve represents r(T ), while the continuous line is λ(T ).At T ≈ 0.40λ 0 , both curves clearly have a discontinuous derivative.

Figure 5 :
Figure 5: Right and left plots refer to the RGZ pressure in terms of T /T c and in units of T 4 .In the left plot, a wide temperature range of is shown.In the right plot, a zoom is made for temperatures around 1.10 T c to show the existence of negative pressure.

4 Figure 6 :
Figure 6: The RGZ interaction measure I/T 4 in units T /T c .