The spatial string tension from the field correlator method

The behaviour of spatial string tension σs(T)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma _s(T)$$\end{document} as a function of the temperature T is found in the framework of the Field Correlator Method (FCM). Here the string tension is calculated using the gluelump Green’s function, where gluons in the gluelump are interacting via the same spatial string tension. The resulting T-dependence was obtained without extra parameters in the region Tc<T<5Tc\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$T_c< T < 5 T_c$$\end{document} using the formalism of elliptic functions θ3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\theta _3$$\end{document}, demonstrating good agreement with available lattice data.


Introduction
The confinement in QCD is a basic phenomenon which ensures more than 90 percent of the visible mass in the Universe and makes the world such as we see it.At zero temperature the theory of confinement in QCD was formulated in the framework of the Field Correlator Method (FCM) [1][2][3][4][5][6], via the vacuum field correlators of the colorelectric (CE) and the colormagnetic (CM) fields E a i , H a i , and at the temperature T = 0 the behaviour of all physical quantities is expressed via the basic nonperturbative parameter -the string tension, which can have different values in the light-like σ E and space-like σ H areas, but at zero temperature T, σ E (T = 0) = σ H (T = 0) = σ.Very important role in FCM plays bilocal correlator (BC) of gluonic fields strength: < tr f Φ(y, x)F µν (x)Φ(x, y)F λρ (y) >≡ D µν,λρ (x, y).
From this moment we use F µν ≡ F a µν T a , a = 1..N f , T a -are generators of fundamental representation of SU(N c ) .In Eq. ( 1) symbol <> means averaging over Yang-Mills .N2 c − 1, Φ(x, y) = P exp(i x y A µ dz µ ), µ = 1..41 is the Wilson line in fundamental representation.We can write BC as follows: D(x−y), D 1 (x−y) -are scalar functions.We also can add index E or H 2 to D, D 1 , because E i = F 0i , H i = ǫ ijk F jk /2, i = 1, 2, 3 and < EH >= 0. Functions D E,H (x), D E,H 1 (x) define all confining QCD dynamics and in particular the string tensions: These functions were calculated in good agreement between the FCM [7] and the lattice data [3,5], while D E (x), D H (x) were also studied in details at T > 0 on the lattice [5].
The most interesting fact is that at T > 0 σ E (T ) and σ H (T ) = σ s (T ) behave differently.Namely: σ E (T ) displays a spectacular drop before T = T c and disappears above T = T c , while in contrast to that σ s (T ) grows almost quadratically at large T , as was found on the lattice [4,[8][9][10] and supported by the studies in the framework of the FCM [11][12][13][14][15][16][17][18][19].Indeed, as it was found in Ref. [9,10] that the dominant part of the spatial string tension σ s (T ) grows quadratically at large where c σ was defined numerically in the lattice calculations [9,10] in case On the theoretical side the quadratic growth of the σ s (T ) was derived in the framework of FCM [7,11,13,14], and the value of c σ was found in Ref. [15] in good agreement with the lattice data of [9,10].
However, in the full FCM expression for the spatial string tension the term in the Eq. ( 4) is only a fast-growing part of the whole expression, which was hitherto not known.
The purpose of this paper is to derive the total expression of the spatial string tension including the linear in T part, to calculate the numerical value of σ s (T ) in the region of [T c ..5T c ]3 and compare it with the lattice data.As will be seen, the results, obtained within the FCM, will provide the values of σ s (T ) in good agreement with the lattice data.In the next section we discuss connection between BC and gluelump Green's function.In section 3 we discuss the general expression for σ s (T ) in terms of the field correlators (and gluelump Green's functions respectively), we formulate its final form and obtain the expression of the string tension, which is discussed and compared with the lattice data in section 4. The renormalized coupling constant g 2 (T ) is given in Appendix A1, the detailed discussion of two-gluelump Green's function is given in Appendix A2.We discuss diagonalization of the two-glulump Hamiltonian in Appendix A3.Calculations of the two-gluelump Hamiltonian eigenvalues and eigenfunctions is given in Appendix A4.

Gluelump and bilocal correlator of gluonic field strength
We need to find a connection between D(x) and so called gluelump Green's function [20].For this purpose we rewrite the expression in Eq. (2) in the form: The integration in the last expression is performed along the straight line connecting the points x and y, so we can rewrite the expression in square brackets in the form [12]: and finally we have: Last equation coincides with the expression from the paper [21].For our purposes we rewrite it as: here all dashed letters stand for operators in the adjoint representation.At the next step we need to find connection between BC and so called gluelump Green's functions .Expanding F µν into abelian (parentheses ) and nonabelian parts: we can write BC as: where the number at the top of the letter D means power minus two of coupling constant g.For D 0 (x, y) we obtain: where ∆ 0 µν,λρ (x, y) contains contribution of higher field cumulants, which we systematically discard.Here G 1g is one-gluelump Green's function: tr adj is a trace over adjoint indices.As shown in Ref. [22], this term is connected with the functions D E,H

1
.From the physical point of view the Eq. ( 13) describes the gluon that is moving in the field of adjoint source (see Fig. 4 in Appendix A.5). Interaction between two objects in the adjoint representation is leading to formation of the string that according to Casimir scaling law found in the framework of FCM in [23] and supported by lattice data in [24,25] have a tension are Casimir operators for adjoint and fundamental representations.This hypothesis give us a chance to calculate one-gluelump mass M 0 .This mass governs the nonconfinig part of the colour Coulomb's potential (we can calculate it from the correlator of Polyakov lines), and from our reasoning it is obvious that we can make the assumption that M 0 ∼ √ σ adj .From direct calculations we have [18]: r-is a distance, T is a temperature, α s is a strong coupling constant.From comparison of the last equation with a simple Debye potential [26] we can say that M 0 is playing the role of Debye mass.
As for D 2 µν,λρ (x, y), it is of basic importance, since ensures confinement via D(x-y) and is expressed via two-gluon gluelump Green's function G 2g (x, y).The expression for D 2 µν,λρ (x, y) reads as: We remind that: Let's consider: We can fix in Eq. ( 17) color indices a, b; d, e and average Green's function over all fields A h µ with h = a, b; d, e.This averaging will produce the white string (of triangle shape at any given moment), and hence it will ensure terms (δ ad δ be + permutations).As a result we can represent G µν,λρ (x, y) in the form: where G 2gl (x, y) is the Green's function of the two-gluon gluelump.Comparison of Eqs. ( 2), ( 15) and ( 18) immediately yields the following expression for D(x-y): Both one-and two-gluon gluelump functions can be written in terms of path integrals [27] and finally expressed via eigenvalues and eigenfunctions of relativistic string Hamiltonian [28].
3 The spatial string tension in the FCM We need to clarify some important moments: at temperatures above deconfinement, T > T c , large spatial Wilson loops still comply with the area law.For pure gauge SU(3) Yang-Mills theory T c = 270MeV .This behaviour of spatial Wilson loops [29] is the main well established nonperturbative phenomenon at T > T c , which is usually called "magnetic " or "spatial" confinement.Of course, it does not contradict true deconfinement of a static quark-antiquark pair [30][31][32][33], because spatial-time Wilson loops indeed lose the exponential damping with the area for T > T c .To illustrate this phenomenon one can calculate polarization operator Π(x, y) in Yang-Mills theory [40].As a standard step we need to decompose gluon field in non-perturbative B a µ and perturbative a b µ parts: usually field B is treated as external.The expression for scalar part of polarization operator reads as: where In Eq. ( 21) averaging is over B a µ .At T > T c , T << σ s R, (σ s = σ H ) one observes dramatic difference for spatial and time-like domains.For the first one (x 0 −y 0 = 0, x−y = R) one obtains at large distance: and at small distance polarization operator tends to non-interacting case: For a time-like domain at zero temperature we have σ E = σ H , but at T > T c σ E = 0 as a manifestation of the confinement-deconfinement phase transition.
With these examples in mind we can focus on calculations of two-gluelump Green's function at non-zero temperature.The Eq. ( 17) describes two gluons moving in the field of adjoint source and interacting nonperturbatively (via σ E,s ) with it and between themselves.
The spatial string tension is proportional to the integral of two-gluon gluelump Green's function in the 3d space, where one of three space coordinates can be treated as an evolution parameter ("'the Euclidean time").Using the technic, developed in Ref. [11,13,14] we can write We neglect the spin interactions in the first approximation.Concerning G (g) 4d we have [14,15]: where and a winding path measure is The important point for the resulting T dependence of the string tension is the integration in the gluon propagator G One can notice that the sum in the Eq. ( 28) is a known function: where the function ϑ 3 (q) is defined as: and thus, the first term in this expression is connected with the vacuum contribution.
Then starting from low temperature there is an expansion: To find the asymptotics at high T one can use the relation: As a result at large T one obtains an equality: Here we use the elliptic functions ϑ 3 (q) defined in Eq. ( 30).Their behaviour as function of q is given in Fig. 1.As a result we can use J 4 at an arbitrary T in the form: In this way starting from low T one obtains an exact expression for J 4 (T ) valid in the whole range of T .One could approximate this behaviour as a sum of linear and constant terms implying a soft transition from T = 0 case to the linear in T behaviour however this approximation fails numerically and actually one observes a sharp transition at some intermediate point T * from the regime T = 0 to the large T behaviour as given by the Eq.(34).For simplification we define: At this point we turn to the general form of the field correlator D H (z) with the aim to express the string tension via the factors f (x).Rewriting the Eq. ( 19) with index "H" we have: where G (2g) (z, T ) is the two-gluelump Green's function.In the path integral representation we can write it (see Appendix A2 for details) as: As a result we obtain σ s (T ) in the following form: One can see in Eq. ( 39) the only T-dependent factors g 4 (T ) and f ( z/2ω 1 T ) which define the dependence of σ s (T ).Therefore one can write σ s (T ) (denoting the z-and ω-integration in Eq. ( 39) with the average sign < ... >) in the following form: where we have denoted the average values of z/2ω (obtained as a result of integration over the T-independent region of parameters with the T-independent kernel.We have also taken into account the symmetries of the Hamiltonian H(ω 1 , ω 2 ) with respect to permutation of ω 1 and ω 2 ) as w = ρ/T c and both ρ, T c are fixed parameters.The appearance of g 4 (T ) which is decreasing with T as (ln T ) −2 defines the T dependence of σ s (T ) to be lower than T 2 , thus confirming the behaviour of σ s (T ) in the lattice data of [9], where the data were fitted as σ s (T ) = constg 4 (T )T 2 .However this fit fails for T < 2T c claiming the necessity of another factor in Eq. ( 41).Correspondingly we are writing the resulting equation for the σ s (T ) denoting the average value of z/(2ω)T as ρT /T c .
In the next sections we try to test our arguments and to demonstrate that this new form with the well-defined factor f (wT ) describes the whole region of T > T c with good accuracy.

General expression for the spatial string tension vs lattice data
For approve our predictions we need to find the parameter ρ that describes the all data from the lattice simulations.For f (ρT /T c) we have: The numerical analysis of the data [9] allows to reproduce well the data with the Eq. ( 41), derived in the previous section Figure 2: Spatial string tension σ s (T )/σ for SU(3) gauge theory as function of T /T c .The lattice data with errors are from Ref. [9].T c =270 MeV The comparison with the lattice data of [9] for Eq. ( 43) is shown in Fig. 2 , and the expression for g 4 (T ) is given in the Appendix 1 and the value of the ρ-parameter ρ = 1/ √ 3.2.The Fig. 2 demonstrates a good agreement between the lattice data and Eq.(43), including the region T < 2.5T c where the lattice fit T 2 g 4 (T ) in [9] starts to disagree with numerical data.
The points on Fig. 2 as fuctions of T /T c completely coincide with the points from Fig. 3 in Ref. [12].This fact means that at sufficiently high temperatures we numerically reproduce the Eq.( 4).We also show the dependence of running coupling as function of temperature (see Fig. 3 in the Appendix A1).As for the dependence of ϑ 3 (e 2 it is equal to one with accuracy of two percent for temperatures higher than 1.3T c .From these data we also can understand that in the absence of running coupling the string tension at high T is proportional to T 2 .Thus one can conclude that running coupling has a small effect in comparison with ϑ 2  3 .From all this facts we can say that temperature of dimensional reduction should be T d ≥ 1.3T c.That doesn't contradict neither the lattice data no earlier FCM predictions [12].

Discussion of results and Conclusions
The main purpose of our work is the construction of the detailed mechanism of the spatial confinement in the whole region of the temperature from T = T c to asymptotically large temperatures.The previous analysis in Ref. [19] has shown that the qualitative behaviour near T = T c can be continuosly connected with the asymptotic behaviour of the σ s (T ).Due to complications of analytic calculations we tried to find the form of dependence of spatial string tension from the main QCD parameters.And for this dependence we found the behaviour that well describes the lattice data.As can be seen in Fig. 2 our resulting curve for the spatial string tension is in a good agreement with the accurate lattice data in the whole measured region T c < T < 5T c .We have exploited the coupling constant depending on the temperature T given in the Appendix A1, which has also allowed the authors of [9] to get agreement with their data in the asymptotic region.Also to describe the region of smaller T we have used the formalism of the elliptic functions ϑ 3 (z) which describe well the sharp transition of the gluon propagator in the two-gluelump Green's function from the constant to the linear behaviour.It should be emphasized that the formalism presented in the paper is standard for the temperature dependence of any Green's functions developing in the spacial or time-like continuum with inclusion of interaction via the field correlators [11,13].In particular the same spatial string tension appears in the expression for the screening mass (we have called it "Debye mass" m D (T ) = 2.06 σ s (T ) in [18,19]) characterizing the spatially oriented parts of the area of the Wilson loops.The inclusion of the temperature via the Matsubara-type formalism with the T-dependent factors I(x 4 − y 4 , T ) = exp(−ip 4 (x 4 − y 4 − n/T )) in the gluon Green's functions is shown in the Eq.(27).For the space-like correlators with x 4 = y 4 the use of the Poisson summation formula 1/(2π) n exp(ip 4 nβ) = k δ(p − 4β − 2πk), β = 1/T , brings about an additional factor of T .This finally leads to the T 2 dependence of the leading term in the σ s (T ) as in the Eq. ( 4).On the contrary for the time-like Green's functions with the nonzero x 4 − y 4 the T dependence is dictated by the corresponding mass parameters and for σ E (T ) the situation is even more dramatic since it drops to zero (deconfinement) at T = T c approximately as (1 − (T /T c ) 4 ) 1/2 [34].One can wonder why these two phenomena -spatial (colormagnetic) confinement and colorelectric confinement are so different and hence disconnected and as follows from the lattice data (see Fig 9,10 in Ref. [5]) the CE gluon condensate < G E 2 (T ) > and the CM condensate < G M 2 (T ) > being equal at T = 0 behave also in a similarly different manner with growing T ?The answer lies in the different active regions of these phenomenathe space-like continuum for CM and the time-like continuum for CE confinement which have a little dynamical intersection as space-like and time-like surfaces, which is evident in the FCM and is an additional argument in favor of its selfconsistency.As it is, we have found good agreement of our FCM approach for CM string tension with lattice data [9] in this paper as well as good agreement of all our CE calculations with the corresponding lattice and experimental data [2-5, 7, 11, 13-15] including the latest CE calculations of the deconfining process [34].Turning back to the CM physics it was found within our approach that an even more important role of the spatial string tension may be in the high T thermodynamics where in the framework of FCM it provides the basic nonperturbative contribution to the pressure and other observables, see e.g.Ref. [17], in good agreement with the lattice data and solving as in Ref. [14] the old "Linde problem" which precludes pure perturbative thermodynamic calculations of interacting systems at large temperature.Another interesting development of this method is the dynamical theory of QCD systems in the external magnetic field where the FCM yields all results in good agreement with lattice data without any parameters, see e.g.[35,36].In this way source with spatial coordinate r = (0, 0, 0) that moves along z-axis entirely.This procedure leads to formation of strings between gluons themselves and glouns and the source.The distance of movement along z-axes is t.This interval must be large enough to avoid contribution of gluon fields fluctuations.With the help of Eq. ( 28) one can integrate out fourth component of Eq. (A2.1) and obtain: S spatial -is a rest part of the action without fourth component containing "spatial" Wilson factor.Changing s = t 2ω 1 , s = t 2ω 2 and Eq.(34) with third coordinate as "Euclidean time" 4 we can rewrite the last equation in the following form: D 3 zD 3 zexp(−S spatial )J(t, ω 1 , ω 2 ), (A2.5) 2 )ϑ 3 (e − ω 2 2tT 2 ), (A2.6) and according to [39] we can average Wilson line over the " euclidean time direction"'.
Thus we obtain: (D 2 z 1 ) xy (D 2 z2 ) xy e − i=1,2 K i (ω i )−V t ϑ 3 (e − ω 1 2tT 2 )ϑ 3 (e − ω 2 2tT 2 ), (A2.7) 4 Such a choice seems a bit surprising but it is a price for possibility to include finite temperatures in the formalism [19].That also means that gluons are moving in the field of the static adjoint source, that is evolving entirely along the z-axis.
For simplicity we assume that ω 1 = ω 2 = ω, ν 1 = ν 2 = ν 3 = ν. 5We calculate all quantities in Eqs.(A3.2), (A3.3) and obtain mass of the lowest two-gluelump: The conditions of minima yields the final result with notation ω, ν for the extremal values.At final step we obtain6 : As for the eigenfunctions, they are given by the product of the eigenfunctions of two noninteracting two-dimensional oscillators.That means that we can obtain all eigenvalues and eigenfunctions that we need in Eq. ( 39).

4 s0 dτ dz 4 dτ 2 ,
over the 4-th direction in Eq. (25) with the exponent K 4 = 1 which gives for the spatial string tension with x 4 = y 4 , and for the temporal string tension with the nonzero x 4 − y 4 completely different behaviour, namely for the σ s case:

Figure 3 :
Figure 3: Behaviour of g(T ) 4 g(Tc) 4 as a function of T Tc .

Figure 4 :
Figure 4: Two-gluelump Green function.Continuous and dashed lines are gluons trajectories.Bold straight line is trajectory of adjoint source.Shaded domain is x-y plane that is perpendicular to z axes.Sides of shaded triangle are fundamental strings.

Figure 5 :
Figure 5: One-gluelump Green function.Bold straight line is trajectory of adjoint source that interacts with gluon (dashed line) through adjoint string).