A0--condensation in quark-gluon plasma with finite baryon density

In the present paper, we return to the problem on a spontaneous generation of the $A_0$-background field in QCD at finite temperature and include in addition a quark chemical potential, $\mu$. We reproduce the known expressions in terms of Bernoulli's polynomials for the gluons and quarks using the same formalism. Then we calculate the $\mu$-dependence, both for small $\mu$ as expansion, and numerically for finite $\mu$. One result is that the chemical potential only weakly changes the values of the condensed field, but quite strongly deepens the minima of the effective potential. The gauge dependence of the condensation is discussed in terms of Nielsen's identity. It also is expressed in terms of Polyakov's loop. Both these representations guarantee gauge invariance of the condensation phenomenon. Also, we discuss the dependence of the Polyakov loop in the minimum of the effective potential on the chemical potential.


Introduction
Investigations of the deconfinement phase transition (DPT) and quark-gluon plasma (QGP) are in the center of modern high energy physics. One way to access these problems is the consideration of vacuum states in the presence of constant background fields. For instance, a constant A 0 -background allows for an investigation of the deconfinement phase transition in terms of the average of the Polyakov loop, (in SU (3)) as an order parameter. At finite temperature, A 0 cannot be zero and for reasons of gauge invariance there is a number of equivalent choices which provide equivalent vacuua. These states belong to the center Z 3 of the gauge group and the properties of, say the effective potential, or the Polyakov loop, (1), under the corresponding transformations, allow to discriminate confined and deconfined states. These ideas were discussed in literature many times beginning with [1] and evolved over the decades. On the one loop level, the free energy F = −β ln Z, belonging to (2), will be the blackbody radiation of the free gluon and quark fields. It has the known Z 3 -structure for the gluon sector, which is broken in the quark sector. If including the next loop, it was found in the early 90-ies that there are non-trivial, lower lying minima. These imply a spontaneous generation of the A 0 -background, or in other words, a condensate. This is similar to the Savvidy vacuum (for a recent review see [2]). While the magnetic case is plagued by instabilities, the A 0 -background is stable.
It is common knowledge that the A 0 -background appears at finite temperature. Due to the periodicity of the fields, it cannot be gauged away and the matrix A 0 = N 2 −1 a=1 A a 0 t a (for SU (N )) can be chosen to be diagonal (and traceless). The chemical potential µ appears as a shift, A 0 → A 0 + i µ g 1 proportional to the unit matrix, with imaginary chemical potential θ, iµ ≡ θ. In the perturbative approach, on the one loop level this picture was understood in the early 80-ies, for instance in the papers [3], [4], [1] (see also the review [5]). The two loop contribution was calculated in the early 90-ies, for gluons in [6], [7], using the so-called 'mesonic (charged) basis'. In [8], overcoming some shortcomings in the literature, it was shown that in two loops the effective potential has non-trivial minima A 3 0 = g(3−ξ) 4πβ , A 8 0 = 0 in the basic sector. The contribution of quarks in two loops was calculated therein. It was found that the minima known from the gluonic sector stay in place and that the effective potential gets lowered.
These results are gauge dependent which, naturally, caused questions. In [9], the gauge dependence of the condensation was shown to satisfy the Nielsen identities. Recently, this topic was reconsidered in [10], showing that the result, expressed in terms of the eigenvalues of the Polyakov loop, is gauge invariant. Also, it was found that this is equivalent to take the gauge ξ = −1. Another, independent, approach is the use of constraint potentials where the condensation shows up in three loop approximation. For a recent account see [11] and literature cited therein.
In the case of quark chemical potential µ the situation is not so well studied. Mainly, because the partition function becomes complex and special considerations are needed. First of all this concerns the lattice calculations which must be real basically and therefore an analytic continuation from imaginary values of µ to real values is applied. These peculiarities are widely discussed in the literature for many years. In fact, they are both theoretical and computational ones. Some new results and current references can be found, for instance, in [12], [13], where the results on condensation of A 3 0 and A 8 0 in the presence of µ are obtained in analytic and numeric lattice calculations based on a model of interacting Polyakov's loops. These are in agreement and correspondence with the ones obtained below.
In the present paper we investigate in detail how a chemical potential, real and imaginary, influences these features. We repeat some known results and present some new ones. For in-stance, using explicit analytic expressions we discuss for the condensates the relation between imaginary and real chemical potential which might be of relevance for the extrapolation of lattice results.
We also derive an expansion parameter which regulates an applicability of the loop expansion theory for the effective potential W ef f (A 0 ) at high temperature, r = g 2 8π 2 , which is sufficiently small for values of coupling close to the confinement temperature. Using the constructed gauge-fixing independent effective potential W L (A cl 0 ), we calculate a thermodynamic pressure of plasma and the screening Debye's mass of gluons at the considered environment.
The paper is organized as follows. In the next section we provide the necessary information about the Lagrangian, the notations and the results of the two-loop free energy in the background potentials A 3 0 and A 8 0 . We also include the quark chemical potential µ. In sect. 3, for µ = 0, we investigate gauge-fixing dependence of the effective potential within the Nielsen identity approach and the effective potential of order parameter and derive the gauge invariant effective potential. In sects. 4 and 5 we consider in detail the dependence of the condensation on a chemical potential, both, real and imaginary. In the last section we discuss the results and give some outlook. Appendix A is devoted to the derivation of µ dependent effective potential in terms of Bernoulli's polynomials. Appendix B contains Feynman's rules and the definition of the periodic Bernoulli polynomials.

Two-loop free energy
In this section we remind the expressions for the effective action in one and two loop order in a constant A 0 -background, which are known so far, and discuss the inclusion of a finite chemical potential µ in the quark sector. The starting point is the QCD Lagrangian with its gluon fields, G a µν , ghost fields, χ and quark fields, Ψ a , in the background R ξ gauge, Euclidean space-time and the presence of a chemical potential µ, where (t c ) a b are generators of SU (3) group, are the components of the background field and holds. It is useful to introduce a "charged (mesonic) basis" for gluon fields [6], [9], These sets of fields correspond to the SU I (2), SU U (2) and SU V (2) subgroups of the SU (3) group. In terms of them the calculations and the representation of results are considerably simplified.
In particular, in this basis the constant background potentials A 3 0 and A 8 0 enter the Lagrangian, and the Feynman rules, in the form of a shift of the fourth momentum component of the color charged fields in (6), with The explicit expressions read Here a 1 is related to the SU I (2) color isotopic spin subgroup, and a 2 and a 3 correspond to the SU U (2), SU V (2) ones. The remaining three charged gluon lines follow with the property The three quark lines undergo shifts with with the notationμ = µ/(2πT ). The explicit expressions are These parameters follow from the covariant derivative in the fundamental representation and are connected by q i = 2c i + 1 with those used in [5]. The c i include the fermionic '+ 1 2 '. The relation to the parameters a i , which follow from the adjoint representation, reads The parameters a i are physically more instructive than the differences q i − q j used frequently in literature. We mention the symmetry under rotation by π/3 in the (A 3 0 , A 8 0 ) plane, or equivalently, in the x, y-plane, of the six frequency shifts of the charged gluons. There is a corresponding symmetry for the quark lines, however with the double angle, 2π/3. Also we mention that the sums of the parameters a i (accounting for (10)) and of c i (except for the addendum + 1 2 + iμ) are zero.
The corresponding Feynman's rules are the standard ones (see also, for example, equations on p. 356, right column, in [6], for their relations to the fields), with the only modification by iµ. We remind them in App. B.
The first calculations of the effective action with a constant addendum to p 4 were undertaken long ago in [14] and [4]. The results for one and two loop levels, formulated in terms of the Bernoulli polynomials, were obtained in [7] and [6], both for gluons. In [8], the two loop contribution from quarks was calculated. In [9] and [15] it was shown that the results are in agreement with Nielsen's identities. More recently, these results were reviewed and reformulated in terms of the Polyakov loop in [10].
The effective potential, including the two loops contributions of gluons and quarks, can be written in the form as sum of the gluonic part (we follow [15], where some misprints in earlier papers were corrected), and the quark part (eq. (2) in [8]) Here, N f is the number of quark flavors. We also take for simplicity µ to be the same for all flavors. It is seen that the one-loop contribution is additive and the two-loop part (∼ g 2 ) includes interference terms of the corresponding color subgroups. Starting from here, the discussion in this section is for µ = 0. In the representations used in literature so far, the polynomials, initially defined in 0 ≤ x ≤ 1, were continued to −1 ≤ x ≤ 0 to obey B n (−x) = (−1) n B n (x) so that (10) can be exploited. This way it was possible to account by three variables, a i (i = 1, 2, 3), for all six charged gluon fields. In the present paper we use a modified scheme. We define, starting from the initial interval 0 ≤ x ≤ 1, the continuation of the Bernoulli polynomials by periodicity in the real part of the argument. The above formulas (15) stay in place which follows from the property B n (x) = (−1) n B n (1 − x) of the Bernoulli polynomials together with the periodicity. In App. A we re-derive the Bernoulli polynomials, using a polylogarithm, where this property becomes obvious.
The effective potential (14) depends explicitly on the gauge fixing parameter ξ. This circumstance caused a large number of discussions about the physical meaning of W ef f and the whole approach with the A 0 -condensation. However, as it was shown in [8] and [15], this effective potential satisfies Nielsen's identity justifying its ξ-independence. This property is established along special characteristic orbits in the (A 0 , ξ) plane, where the variation in gauge-fixing parameter is compensated by corresponding variation in fields. At the same time, it is desirable to have a function which would be independent of ξ and explicitly written in terms of gauge invariant observable (Polyakov's loop, usually). This idea was expressed by Belyaev [16] and we call this object the "effective potential of order parameter", W L (A cl 0 ). However, Belyaev obtained the value A cl 0 = 0 in the two-loop order calculations. Below, in sect. 3 we show that the choice ξ = −1 delivers effectively the gauge invariant results.
Consider the effective potential (14) as function of x and y. It has the known minima. In the main topological sector, there is one minimum, located at and the effective potential in this minimum reads where the first two terms correspond to the zero field case (these are present also without condensation and represent the blackbody radiation) and the last term is due to the condensation. It is negative and lowers the free energy. The other minima are located on a circle with radius x and form a hexagon (in [4] this was mentioned as permutation symmetry). This is the way how, for the condensates, the above mentioned rotation symmetry is realized, whenever this is not obvious for the quarks since the c i rotate by 2π/3. In Fig. 1, in the right panel, we show these 6 minima, which appear in the two loop approximation. For comparison, in the left panel we show the effective potential for g = 0, where we have only the 'trivial' minimum at x = 0. By periodicity resulting from the Matsubara frequencies, all these minima are replicated on a super lattice given by the shifts (n, m are integers) The replications of the minima in the pure gluonic case realize the Z 3 -symmetry. All these minima have the same depth. As we see, the depth of a minimum increases with increasing number of flavors. As first mentioned in [8], the minima of the pure gluonic effective potential are not moved by the inclusion of the quarks. However, as can be seen from (21), the replication happens in the y-direction with a three times larger period ( √ 3 in place of 1 √ 3 ) such that only a third of them remains. This way, the Z 3 -symmetry is broken by the quarks. This difference is, of course, due to the different group representations; associated for the gluons and fundamental for the quarks. The above discussion illustrates the commonly expected (e.g., [1]) breakdown of the Z 3 -symmetry by the quarks in loop expansion. We demonstrate this feature in Fig. 2.
Finally in this section, we express the Polyakov loop (1) in terms of the condensates. We have to insert, using the known SU(3) generators t c , (4) and (8) into (1) and to get L in terms of x and y, For x and y we insert the condensates (20) and (21) and arrive at As will be shown in sect. 3, one has to use (18) with ξ = −1. It is seen that on the one-loop level, i.e., for x (0) min = 0, we have L = 1 and on the two loop level with x (0) min given by (18), a smaller value, L 1. It is to be mentioned that expression (23) remains unchanged under (20) and (21), i.e., it takes the same value in all minima related by these formulas.

On the gauge independence of the A 0 -condensation
When the minima of the effective potential in the A 0 -background were found, their dependence on the gauge fixing parameter ξ, see (18) and (19), raised questions whether these can be physical. There were two ways found to get an affirmative answer. One way involves the idea to restrict the integration space in the functional integral representing the effective action to such fields, which have a given value of the Polyakov loop (1). This approach rests on the ideas expressed in [17] and is actively developed. We mention the very recent paper [11]. The other, much more direct way, was suggested in [16]. In [9] it was related to the Nielsen identities, which are a special case of the Ward-Takahashi (Slavnov-Taylor) identities. The basic statement is that a change in the gauge fixing parameter can be compensated by a change in the background field, see, for instance, eq. (1.2) in [18] or eq. (1.1) in [19]. The Nielsen identities cast this statement into formulas. We adopt this approach here, in a much simplified form, to get a self contained representation and clarify some deficiencies in the literature.
We start from the equivalent statement in [16] that the effective potential in the minimum (at zero sources) is a gauge invariant object. Now, in our case it depends on the background field A 0 and on the gauge fixing parameter ξ. Since A 0 is a gauge potential, it is not physical (measurable) and may itself depend on ξ, A 0 = A 0 (ξ). But, its dependence must be such that the effective potential is independent on ξ, This way, we expect a set of curves in the (A 0 , ξ)-plane. We mention that this statement was derived in [16] using the Nielsen identities. Now, in order to get the mentioned curves, we apply (24) to the effective potential (14). We restrict ourselves to the order g 2 in the perturbation expansion. Further, for simplicity, we restrict ourselves (for the rest of this section) to the SU I (2)-case, where the effective potential reads  (21) are seen, especially the lifting of minima when including the quarks.
For notation simplicity we consider x = x(ξ) (in place of A 0 ), assume ∂x ∂ξ ∼ g 2 and remind, after (9) and (13), the relations a = x 2 = gA 3 0 2πT and c = x 4 + 1 2 . Then we have, to order g 2 , from (24), Carrying out the derivatives, using (25), we note and Using the recursion relations, B n (a) = nB n−1 (a), for the Bernoulli polynomials, and inserting into (26), we arrive at The expression factorized and the equation for the curve in the (A 0 , ξ)-plane is simply To get the solution of this equation to order g 2 , we may substitute in the right side for a = x 2 some tree value, x cl , and integrate. As a result we get where ζ is an integration constant, characterizing the gauge orbits. The solutions are parallel straight lines (at the given approximation). On the tree level we have x = x cl and, of course, no ξ-dependence. More precisely, inserting this relation into (25) and expanding in a series of g 2 , we could obtain a ξ-independent expression for the effective potential as a function of the characterization ζ of the gauge orbits. A second line of reasoning is based on the circumstance that the average < L > of the Polyakov loop (1) is gauge invariant. Following [16], we calculate this average perturbatively in order g 2 (one loop). On the tree level we have x cl = gA cl 0 πT and holds. Including loop corrections, this relation turns into < L >= cos π 2 x =< L > tree +∆L.
For the correction ∆L, from the known Feynman rules, the expression and k = 2πT ( + a) follows. For details we refer to [16], eq. (9). Such kinds of expression were calculated in the literature repeatedly. However, we found it instructing to give a version in terms of zeta functional regularization. Thus we consider the expressions with with s = 0 at the end. We integrate over the spatial momenta and, in the first integral, we integrate in t by parts. For the factor k in front of the exponential we use a derivative with respect to a. We arrive at and I 2 = Γ(s + 1)(s + 1 2 ) Γ(s + 2) This way, from (36) we get Further we apply Poisson re-summation to I 1 , Further we drop the T = 0 contribution (this is N = 0) and can now put s = 0. We arrive at Carrying out the t-integration, we may use the known formula where B 2 (a) is the periodic continuation of a Bernoulli polynomial. Inserting into (42) gives Finally, inserting into (34) and using (41) and (40) we come to and with (33) to Comparing with (31) we come to the conclusion This formula can be compared with (31) and we find for the integration constant ζ = −1. Just for this value an arbitrary characteristic coincides with the curve corresponding physical "observable" x cl . As a result, to express (25) in terms of Polyakov's loop we have to set ξ = −1 and x = x cl in the initial effective potential.
Another important consequence of (48) is the ratio of two-and one-loop contributions r = W (2) W (1) determining the expansion parameter of the theory, r = g 2 8π 2 . It is small for sufficiently large coupling values g ≤ 4 − 5, and the two-loop effective potential (25) is applicable till the confinement temperature where Polyakov's loop turns to zero.
The relation between the x and the x cl (48) can be immediately extended to U − and V − subgroups (these are the fields K ± ν ,K ± ν in (6)). We have just to replace a 1 = x 2 by a 2 and a 3 (9), correspondingly. This possibility is due to the obvious fact that the calculation procedure resulting in (48) is independent of the actual value of a 1 = x 2 . More detailed derivation of (31), (and hence (48)) can be found in [9] for SU (3) case and in [15] for full QCD.
Thus, to derive the gauge-fixing independent effective potential expressed in terms of the Polyakov loop we have to set in (15) a 1 = a cl 1 , a 2 = a cl 2 , a 3 = a cl 3 and ξ = −1. From the obtained effective potential, the observable values of the condensed fields are determined from their minimum positions.

A 0 condensation in the presence of µ
In this section we investigate how the picture described in the preceding sections changes under inclusion of a chemical potential µ, i.e., of a finite density of quarks. We present the results obtained in two forms, one for the gauge-fixing dependent effective potential, in order to verify Nielsen's identity. The other form is the the gauge-fixing independent effective potential W L (A cl 0 , µ). On the level of the Feynman rules, the inclusion of µ is shown in (12). The structure of the Feynman diagrams remains unchanged and formally we arrive at the expressions (16), where µ enters through the variables c i , (12). These formulas can be viewed as a kind of analytic continuation in the variables A 3,8 0 , or (x, y), in (12). The question arises whether this is justified, especially having in mind that the B n are not really polynomials. The answer is yes, as we demonstrate in Appendix A. For instance, it can be seen that the periodic continuation is to be done in the real part of the argument. In doing so, obviously the symmetries (20) and (21) stay in place.
Including the chemical potential this way, we first consider the µ dependence analytically. As in sect. 2, we find the minimum of the effective potential in the fundamental sector in the (x, y)-plane, namely its position and depth, and expand the resulting expressions for small g. It must be mentioned that the effective potential (14) is a fourth order polynomial in x and y. Its roots are solutions of third order equations and contain third order square roots. First we do an expansion of these expressions in powers of g.
For the position of the minimum, as it turns out, we have now that (18) is kept for N f < 9, i.e., there is no dependence on µ for known values of flavors. However, the effective potential (19) depends on µ. It is a fourth order polynomial in µ, where W 0 is the same as (18) (recallμ = µ/(2πT )). We see that the chemical potential, for small µ, does not change the condensate, but deepens the minimum. There are odd powers of µ (and with them an imaginary part for real chemical potential). A plot of the real part of this potential is shown in Fig. 3. Again we note that the chemical potential enters (14) through the arguments of Bernoulli's polynomials in (12). This does not change the structure of the effective potential that satisfies Nielsen's identity derived in [15], [8], and gauge invariance is preserved. In order to obtain the gauge-fixing independent effective potential corresponding to Polyakov's loop, we have to substitute A 0 → A cl 0 , ξ → ζ = −1, as it is described above. These substitutions are equivalent to put ξ = −1 in (49). The emerging expression can be used for calculating thermodynamic pressure in the plasma which is determined as: p = −Re(β 4 (W q + W gl ) |min ). Hence it follows that chemical potential increases the pressure.
In (23) we have seen that the Polyakov loop (1) depends of the condensate, see (22). Following the µ independence of the condensates, L is also independent on µ.
We complete this section with the calculation of the Debye screening mass for neutral gluons. This mass is defined as It is a function of the background, x (we put y = 0 for simplicity), N f and µ. This function is a polynomial in µ. For x = 0 it reads m 2 D,0 = 1 + In the minimum (18) of the effective potential the Debye mass is The difference between them is m 2 As already, to obtain the Debye mass expressed in the terms of the Polyakov loop we have to set ξ = −1 in the above expressions. Finally, we mention the special case when µ = 0, N f = 0, which is the pure gluonic case.

Condensation with imaginary chemical potential
In this section we consider imaginary chemical potential. In the above formulas this is reached by the substitution µ → iµ (andμ → iμ). In that case the addendum to the frequency (11) is real. As such, it fits to the periodic continuation of the Bernoulli polynomials and the effective potential will be periodic in this µ (with period 1 forμ). In the perturbative results in Section 3 one has to do simply the mentioned substitution. The first order (in µ) correction will dominate and its sign depends on the sign of µ. The position of the minima shows a weak dependence on µ. There is also a small contribution to y. It is shown in Fig. 5. It should be mentioned that this dependence is for finite g, in distinction from the preceding section where we considered the perturbative expansion in g.
We provide also some pictures of the effective potential for imaginary µ, see Fig. 6. As can be seen, the role of the minima changes as function of µ, cycling with period of approximately 1/3 (forμ). We demonstrate this by plotting the depths of three consecutive minima in Fig.  7 (left panel).
The relation between real and imaginary potential is in the above formulas in complete analogy with the relation between trigonometric and hyperbolic functions. This can be most clearly seen by looking on the left side of (61), i.e., before the application of Jonquière inversion formula. Although this is on some stage nothing more than some analytic continuation, the results are quite different and hard to connect post factum.
We turn to the Polyakov loop (1). In the presence of imaginary chemical potential it depends, of course, also on µ. Inserting the numerical µ-dependence, which is shown in Fig.  5, into L, we can represent this dependence graphically. This is done in Fig. 8. We see that the change with respect to eq. (23) is weak. An imaginary part shows up, which is small since it results from the y-dependence shown in Fig. 5. The Fig. 8 is done for the minimum in    (20), for odd k the imaginary part changes sign, the real part stays in place.
We complete this section with the following speculation. Assume, the effective potential with imaginary chemical potential had some physical relevance. Then we had, for different µ, equivalent minima. Hence, there should be some tunneling between them. If we take a simple average of the effective potential over these µ, we get the picture shown in Fig. 7 (right panel). This would be a way to restore the Z 3 symmetry including quarks.

Conclusions and discussion
In QCD, on the base of Nielsen's identity and the effective potential of order parameter W L (A cl 0 ) calculated in sects. 2, 3, we investigated the influence of chemical potential µ on the condensation of color constant potentials A 3 0 and A 8 0 at high temperature after the chiral phase transition when quark masses are m f = 0.
For µ = 0, in two loop order we confirm the known results for condensation, A 3 0 = 0 and A 8 0 = 0 to order g 2 , in the main sector, eqs. (18), and for the minimum of the effective potential (19), to order g 4 . We give explicit formulas for all other sectors, (20) and (21).
For non-zero µ, within the perturbative expansion, we observe no dependence of the condensate on µ. Not expanding in g, we observe a weak dependence of the condensates on µ (see Fig. 5). We see a deepening of the minima on µ (which is not small); see Figs. 3 and 4..
Within the perturbative expansion, we calculated the Debye mass, (50), in the minimum of the effective potential, (52).
Also we investigated imaginary chemical potential and found that, including the fermions, there is no Z 3 -symmetry in our approach, in accordance with general expectations.
We obtain that a real chemical potential decreases the effective action and, to some extend, amplifies its dependence of the condensates. The dependence on imaginary chemical potential is, of course, periodic. Examples are shown in Fig. 6 and Fig. 7. We mention that the expansion of the effective potential in powers of µ has, again in the given approach, odd contributions.
We mention that the effective potential in a minimum is not even in the chemical potential, see for instance the linear term in (49). This could pose an additional problem in attempts of analytic continuation from imaginary to real chemical potential as it was done, for example in [20]. We mention that an A 0 -condensate will be present at least in the perturbative approach, for instance also in a quark-gluon plasma in the parameter regions where this approach is reliable. Since this background violates Z 3 and color symmetries, new effective vertexes and related with them phenomena have to exist. Possibly, these open the possibility for new signals of the deconfinement phase transition. It is interesting and important that an increase of µ increases the corresponding effects. Some of these phenomena (such as decay of gluons into photons, etc.) were discussed in [21], [22].
In [1], the Z 3 -symmetry of the effective potential was discussed. It should have the property Z(µ + 2πT 3 ) = Z(µ), where µ is the imaginary chemical potential (in [1] the notation Θ was used), if only color singlet states contribute to the functional integral (2). In a loop expansion, as we do, this symmetry should be absent. This feature is clearly confirmed by our results, see especially Fig. 6. The mechanism is that the presence of the quarks makes two third of the gluonic minima (Fig. 2, right panel) less deep. Now, increasing µ by 2πT 3 , one cycles through the panels of Fig. 6. Accordingly, the system moves from one minimum to the next. In Fig. 7, left panel, the depths of 3 neighbored minima are shown as function of the imaginary chemical potentialμ. It is seen that at most two of them may be of equal depth. We mention a further symmetry which, possibly, was not observed so far. As function of x, for y = 0, for the effective potential (2) the relation Z(n + 1 3 2πT ) = Z(n + 2 3 2πT ) (n-integer) holds. This property can also be spotted on Fig. 6.
It is only a speculative tunneling between these minima that could restore the Z 3 symmetry as shown in Fig. 7, right panel. We mention that one observes a different picture at lower temperature, as is discussed in many lattice investigations, see [20] for example.
On the base of the obtained gauge-independent effective potential W L (A cl 0 ), we have calculated thermodynamic pressure p and screening Debye's masses of gluons m 3 D and m 8 D at the considered environment. Now, we compare the results obtained with the calculations existing in the current literature. In [23], the two-loop constrained effective potential for QCD at finite temperature was calculated for either massless or massive quarks and chemical potential. Therein the goal to investigate a spontaneous generation of gauge fields was not intended. Nevertheless some comparisons are of interest and possible. First, as concerns µ dependence of the twoloop effective potential calculated in [23], it is similar to that derived in (16). But no gauge field condensation was detected. Second, in [23], [11] it was noted that the methods of Nielsen's identity and effective potential of order parameter W L (A cl 0 ) are an alternative to the constrained effective potential. But in fact they are mutually related. This is because in all of them the key idea is realized -to calculate an addition part for the two-loop ξ-dependent effective potential, which compensates the change of the tree-level Lagrangian followed when the gauge-fixing term is introduced. In Nielsen's identity approach, this is realized in the form of the special variation of gluon fields removing the variation ξ → ξ + δξ of the gauge-fixing term. There is the set of such orbits (parameterized by an integration constant ζ) in the (A 0 , ξ)-plane along which the effective potential is independent of ξ. In case of W L (A cl 0 ), the value of ζ is chosen in such a way that the A cl 0 and A 0 belong to Polyakov's loop. This idea was proposed by Belyaev [16]. It is very important that the characteristic and the one-loop variation of the PL ∆ L are described by the same equation [9]. We have investigated this point in detail in Sect. 3 and explicitly shown that the factor regulating this contribution should be ∼ g 2 (ξ + 1)B 1 (a i ) instead of ∼ g 2 (ξ − 3)B 1 (a i ) as obtained by Belyaev (and also in [23], [11] for the constrained effective potential). Actually, just (ξ − 3) factor results, in particular, in the A cl 0 = 0 in two-loop order [16] (same also in [11]). This is a special gauge where there is no renormalization of fields in order g 2 . Detailed discussion of these (as well as other related with them) peculiarities were given in our papers [15], [8], [9] devoted to the Nielsen identity method. In [11], the gauge field condensation has been detected in order g 3 , as it was assumed in [9], [24].
In [12] and [13], the results on condensation of A 3 0 and A 8 0 in plasma with chemical potential have been obtained in the lattice calculations based on a Polyakov loop model. Thus, physical result is the same in all the calculations discussed. It confirms that the A 0 condensation takes place and that it regulates the infrared dynamics of the fields at high temperature. It worth also to mention that the latter results have been obtained as expansion over g −1 in a large coupling constant approximation which is reliable at temperatures not much above the deconfinement temperature T d . Therefore these calculations are complementary to ours. They support our results when one moves from low to high temperatures.
We would like to complete with the conclusion that the A 0 background has to result in new interesting phenomena proper to the plasma. The presence of chemical potential should amplify these. It is to be mentioned that the physical consequences of the A 0 -condensation, without and with chemical potential, are underrepresented in the modern discussion of quarkgluon plasma. We hope, that the present, detailed investigation will promote such discussions.
As an outlook we mention that it would be interesting to see whether the A 0 -condensation is stable against the inclusion of magnetic background fields. Possibly, it could even stabilize the magnetic background as discussed recently in [10].
In the case µ = 0, i.e., without chemical potential, we have the well known expressions