Confinement in a three-dimensional Yang-Mills theory

We show that, starting from known exact classical solutions of the Yang-Mills theory in three dimensions, the string tension is obtained and the potential is consistent with a marginally confining theory. The potential we obtain agrees fairly well with preceding findings in literature but here we derive it analytically from the theory without further assumptions. The string tension is in strict agreement with lattice results and the well-known theoretical result by Karabali-Kim-Nair analysis. Classical solutions depend on a dimensionless numerical factor arising from integration. This factor enters into the determination of the spectrum and has been arbitrarily introduced in some theoretical models. We derive it directly from the solutions of the theory and is now fully justified. The agreement obtained with the lattice results for the ground state of the theory is well below 1% at any value of the degree of the group.


Introduction
A deep understanding of Yang-Mills theory in all the range of the coupling represents a fundamental aspect of our comprehension of strong interactions. The reason is that this would open the possibility to accomplish computations of the behavior of the theory in the low-energy limit where the theory displays bound states. Currently, the only way to obtain results that are derived directly from the theory is through extensive use of lattice computations on large computer facilities. This has permitted for the Yang-Mills theory to obtain both the spectrum and the behavior of propagators in several gauges also at finite temperature [1,2,3,4] in four dimensions and similarly for the case Marco Frasca Via Erasmo Gattamelata, 3 00176 Rome (Italy) E-mail: marcofrasca@mclink.it d = 2 + 1 [5,6,7,8,9,10]. Specially in this latter case, very precise results exist for the string tension and the spectrum.
From a theoretical standpoint, the situation appears decisively better for the three-dimensional case where some analysis have been performed producing excellent agreement with lattice computations for the string tension [11,12,13,14] and the spectrum [15,16]. The starting point was a work by Karabali, Kim and Nair that proposed a proper set of matrix variables to work with in this case to put forward a wavefunction and derive fundamental results of the theory [17,11,12]. Karabali, Kim and Nair approach appears greatly successful in the derivation of the string tension and higher order corrections [14]. For the spectrum, a different wavefunction was postulated [15,16] always in the framework of Karabali, Kim and Nair formalism. Again, the agreement with lattice data was impressive.
In a recent paper of ours we were able to build a quantum field theory for the self-interacting scalar field in the limit of the coupling running to infinity [18]. We displayed a set of classical solutions that, notwithstanding we started from a massless equation, showed a massive dispersion relation. These solutions were already proposed in [19] but the idea in [18] was to consider them as the vacuum expectation value of the field and build the quantum theory from them. In this way one has that conformal invariance is broken and a zero mode appears. The particles get a mass and a tower of excited states described by the spectrum of a harmonic oscillator. This theory shares a trivial infrared fixed point and an ultraviolet trivial fixed point making the theory overall trivial but with a mass gap. On this ground it is a natural question to ask if also a Yang-Mills theory can share such classical solutions and a corresponding quantum field theory built upon them. The answer was affirmative as we showed in [20] but this is true asymptotically in the general case while the result holds exactly just in the Lorenz (Landau) gauge. The corresponding quantum field theory develops a mass gap but is trivial at both sides of the range due to the trivial infrared fixed point and asymptotic freedom on the other side. This scenario has received some confirmations recently in a work by Deur [21]. The idea is to use the propagator of the scalar theory and compute the potential with a running coupling as expected in a Yang-Mills theory. The agreement with a confining potential obtained from lattice computations is striking. The running coupling in Yang-Mills theories has been widely discussed in [22,23,24,25,26,27,28,29].
In this paper we develop this approach, alternative to the Nair and Karabali formalism, deriving all the properties of the theory. That is, we solve the Yang-Mills theory in 2+1 dimensions in another way and we will get results in strikingly good agreement with lattice data, validating this approach. As a byproduct we will get an important hint that the flux tube description of hadron emerging in AdS/CFT approach [30] is a successful one as we are able to get the right ground state of the theory by introducing the same factor as in the Isgur-Paton theory [31] as demonstrated by Teper and Johnson [32]. Anyhow, it is important to point out that we show how the correction factor, arbitrarily introduced in the aforementioned works, is here properly obtained by solving the classical equations of motions, being this just an integration constant. This kind of arbitrariness enters into Yang-Mills theory and the scalar field theory, to which it maps, due to the properties of the differential equations of the theories. Another result that is really interesting with our approach is that our confining potential is almost identical to the one obtained in [15,16]. These authors obtained it after some hypotheses to be verified and we show here that they were correct. The three-dimensional theory is so proved to be marginally confining.
The value of this analysis can be just of mathematical interest to study the structure of a quantum field theory in lower dimensions and to get some understanding of the behavior of the four dimensional theory starting from a simpler case. Indeed, as stated above, a lot of effort has been spent, both through lattice and theoretical study, to analyze the three dimensional case. Indeed, some of the mathematical techniques devised for this case could be successfully applied to the more realistic four dimensional case.
The paper is so structured. In Sec. 2 we introduce a set of exact classical solutions to work with for the quantum field theory. In Sec. 3 we derive the gluon propagator for the classical theory. In Sec. 4 we discuss the ghost sector of the theory. In Sec. 5 we evaluate the quantum corrections. In Sec. 6 we compute the Wilson loop obtaining a confining potential in agreement with literature and the string tension in agreement with lattice data. In Sec. 7 we comnpute the ground state of the theory giving the lowest glueball state. Finally, in Sec. 8 the conclusions are presented.

Classical solutions
Motion equation for Yang-Mills theory can be straightforwardly written down for any number of dimensions and SU(N) group in the form [33] the covariant derivative, T a the generators of the group and A a µ the potentials (a, b, c, . . . are color indexes running from 1 to N 2 − 1 ), and the field components with F µν = T a F a µν and f abc the structure constants of the group. As our aim is to work out a result in quantum field theory, we add a term into eq.(1) to fix the gauge in the form with ξ a free parameter determining the gauge choice.
Using perturbation theory, one can show that there exists a set of solutions of Yang-Mills equations of motion that can be cast in the form [20] being η a µ a set of constants to be determined depending on the problem at hand (e.g., for SU (2) in the Landau gauge, one can take η 1 1 = η 2 2 = η 3 3 = 1, all other components being zero). Putting these potentials into the equations of motion yields [20] ∂ These solutions become exact and not just perturbative for the Lorenz (Landau) gauge. An interesting aspect of these solutions is that hold in any dimensions d > 2. For d = 2 Yang-Mills equations of motion are trivial and no such solutions can be found. Without exploiting all the possible solutions of eq.(6) we limit our interest to a subclass of solutions that have the property to be massive even if we started from massless equations of motion. We have fully exploited this case in Ref. [18]. In this paper we consider such exact solutions as a ground state of the quantum field theory of a scalar field. In 3+1 dimensions this can be written down as [18] being sn a Jacobi elliptic function, φ an arbitrary phase, µ an arbitrary constant having the dimension of a mass and provided that So, if we interpret k as a four-vector of momenta, this can be seen as the dispersion relation of a massive wave. These solutions are rather counterintuitive as we started from a pure massless theory. A mass term can be seen to arise from the nonlinearities of the equations we started from. In the following we will assume that such solutions are just the ground state for the quantum field theory we aim to study. In 2+1 dimensions Yang-Mills equations have a coupling g 2 having the dimension of a mass or inverse of a length. This means that our solution takes the form being a an arbitrary dimensionless constant to be fixed in the quantum theory and φ an arbitrary phase. This holds provided the following dispersion relation holds Confinement in a three-dimensional Yang-Mills theory 5 3 Gluon propagator We need to introduce the propagator of Yang-Mills theory in the infrared limit. This is generally accomplished by a current expansion [18,34]. Instead to start from the action, we prefer the equations of motion [35] Then, we assume a functional form A a ν = A a ν [j] and perform a Taylor expansion around the asymptotic solution (5). We have to take in mind that, for the Landau gauge, these solutions are exact but just asymptotic for whatever other gauge choice. So, we take in general We are assuming here that eq.(5) represents the ground state of the theory i.e. A a ν [0] = η a ν χ(x). These describe oscillations around a vacuum expectation value of the fields as seen from our solutions. Then, the propagator of the theory will be We can obtain the corresponding equation by doing the functional derivative on the equation of motion. We get Imposing j = 0 one obtains the following equation for the Green function of Yang-Mills theory In order to compute the propagator, we perform a gauge's choice. The most common is the Landau gauge (α = 1) that also grants that we are using exact formulas rather than asymptotic ones. So, we write as usual for this gauge being p µ the momentum vector. This yields for the above equation that is the equation we were looking for. This equation coincides with that of the Green function of the scalar field obtained in [18] in agreement with the mapping we derived in [20] provided λ ↔ N g 2 , being λ the corresponding coupling for the scalar field theory. We now limit our analysis to the case d = 2+1 and compute the exact Green function for this problem. The technique we follow is that outlined in Ref. [18]. We just note that we have two independent solutions of the homogeneous equation or One is with cn and dn elliptic Jacobi functions, that holds provided The other one can be obtained by writing it as with Now, we introduce a new variablex = p · x + φ and use the dispersion relation (22) to obtain where the primes mean derivative with respect tox. From eq. (21) we can obtain the solution in the rest reference frame p 1 = p 2 = 0 and p 0 = aN g 2 /2 1 4 . The corresponding Green function is where we have set µ 0 = aN g 2 /2 that fixes the mass scale, and θ(t) is the Heaviside function granting that the solution is different from 0 at t > 0 and 0 for t < 0 and provided that cn(φ, −1) = 0. Similarly, one can define a backward propagating Green function as So, the propagator is When we turn to a Fourier transform, Fourier series of Jacobi functions are well-known [36] giving and so one arrives, back to the moving reference frame, at the result [18] being K(−1) the complete elliptic integral of the first kind and we get the "mass spectrum" At this stage this has just a formal meaning. Moving to quantum field theory, we will prove that this is indeed the spectrum of the theory. So, our final result for the Green function in d = 2 + 1 is This result implies that the Yang-Mills theory shows up a mass gap also in this case. The corresponding spectrum can be used to fit with lattice data.

Ghost sector
As shown in [37] in four dimensions, starting from the exact solutions given in Sec. 2 for the 1-point function in the Dyson-Schwinger set of equations, the ghost propagator reduces just to the one of a free massless theory. This signals that the ghost sector decouples from the physical degrees of freedom. For the sake of completeness, we give here the corresponding Dyson-Schwinger equation that is where we identify the 3-point and 2-point function K bcm 3µ (0, x− y), K cm 2µ (x− y) for the ghost-gluon field propagation and the 1-point function P b 1 (x) for the ghost field. With the given solutions for the 1-point function, this just boils down to the propagator for a free field.

Quantum corrections
We want to see how quantum theory modifies the one-and two-point functions we obtained in the classical theory. This can be accomplished using the Dyson-Schwinger equations. We will stop the analysis to the two-point function as already discussed in our recent work [37]. In that paper it is shown that the mass should be renormalized by adding the term (given in d dimensions) where the dimensions of the coupling g grant that of the squared mass. This integral can be exactly evaluated to give This correction diverges for d = 4 while is finite for d = 3. This should be expected for Yang-Mills theory in three dimensions [38,39]. It evaluates to with S 0 = 2.046970223 . . . the result of the sum. This boils down to add a numerical constant to the arbitrary parameter a. So, one can always redefine the factor a in such a way to compensate the numerical factor obtained in this way and we will have for the spectrum

Wilson loop and potential
In order to compute the potential in a pure Yang-Mills theory at the infrared fixed point, we have to evaluate (40) being S g [c, c] the contribution of the ghost field, T a the anti-hermitian generators of the gauge group and P the path ordering operator. In our case, in the infrared limit, we have a trivial fixed point and the contribution coming from the Yang-Mills field is just a Gaussian one. This implies that our generating functional takes also a Gaussian form and the Wilson loop has the simple form being C 2 (R) the quadratic Casimir operator that for SU(N) in the fundamental representation, R = F , is C 2 (F ) = (N 2 − 1)/2N . This yields being We aim to recover the result given in [15,16]. So, in our case we have to evaluate the integral [40] being J 0 a Bessel function and with the propagator given in eq. (31). The integral can be computed exactly giving being K 0 a Bessel function. We recognize here the potential obtained, after some hypotheses, in [15,16]. What is changing is the mass scale but this should be expected due to our approach that involves exact solutions to the classical equations of motion. This potential grants that the three-dimensional theory is marginally confining [15,16]. The reason relies on the fact that K 0 (m n r) ≈ − ln(m n r/2) at smaller distances and decreases really slow with the distance making the Wilson loop not strictly proportional to the area. Anyhow, a logarithmic potential grants that the Gauss law is satisfied in 2+1 dimensions yielding colorless states. We notice that with σ KKN = N 2 g 4 8π (1 − 1/N 2 ) the Karabali-Kim-Nair string tension. Written in this way, we can compare it with the same result given in [15,16]. We get a renormalized Karabali-Kim-Nair string tension given by The value of the arbitrary factor a is irrelevant here as the square root of the string tension determines the spectrum through the ratio m n / σ R KKN and this factor enters also into m n absorbing it. Then, in the fundamental representation, we can take a = 1 to compare with lattice data and one has √ σ a. This result agrees within an error of about 2% with respect to lattice computations for the leading order in 1/N [10].

Glueball spectrum
The glueball spectrum for 0 ++ is easily obtained through the equation Now, the factor a entering into the string tension simplify with the one in the spectrum giving an overall value that, in some models in literature, was identified as a fudge factor but that here is fully justified by the exact solutions of the theory. Indeed, the spectrum yields m n √ σ = (2n + 1) · 5.032050686 . . . · √ a 1 The agreement with the ground state of the theory is reached for a = 2/3, assuming that the lattice data are affected by errors. In this case, the following The agreement is strikingly good being well below 1% error for the ground state at any N . The factor emerging from the analysis of the ground state of the theory is in agreement with similar factors introduced in literature [31,32] but now theoretically well founded.

Conclusions
We have shown that, in the framework of our formalism, marginal confinement is achieved for QCD in three dimensions. We have found extensive agreement with lattice data and preceding theoretical works. We have also shown that numerical factors arbitrarily introduced in some models are completely justified by the set of classical solutions we have chosen to start with. The exceptionally good agreement between lattice data and theoretical predictions we achieved in the present case can serve as a justification a posteriori for the choice of the solutions to start quantum field theory.