Confinement Mechanism in the Field Correlator Method

Abstract

Confinement in QCD results from special properties of vacuum fluctuations of gluon fields. There are two numerically different scales, characterizing nonperturbative QCD vacuum dynamics: "small" one, corresponding to gluon condensate, critical temperature etc, which is about 0.1-0.3 GeV, and a "large" one, given by inverse confining string width, glueball and gluelump masses, and so forth, which is about 1.5-2.5 GeV. We discuss the origin of this hierarchy in a picture where confinement is ensured by quadratic colorelectric field correlators of the special type. These correlators, on the other hand, can be calculated via gluelump Green's function, whose dynamics is defined by the correlators themselves. In this way one obtains a self-consistent scheme, where string tension can be expressed in terms of _ΛQCD .

(ProQuest: ... denotes non-US-ASCII text omitted.)

Recommended by Frederik G. Scholtz

State Research Center, Institute of Theoretical and Experimental Physics, B. Cheremushkinskaya 25, Moscow 117218, Russia

Received 2 July 2009; Accepted 3 November 2009

1. Introduction

Confinement of color is the most important property of quantum chromodynamics (QCD), ensuring stability of matter in the Universe. Attempts to understand physical mechanism of confinement are incessant since advent of constituent quark model and QCD; for reviews, see, for example, [1-4]. Among many suggestions, one can distinguish three major approaches as follows.

(i) Confinement is due to classical field lumps like instantons or dyons.

(ii) Confinement can be understood as a kind of Abelian-like phenomenon, for example, according to seminal suggestion by Hooft [5] and Mandelstam [6] of dual Meissner scenario.

(iii): Confinement results from properties of quantum stochastic ensemble of nonperturbative (np ) fields filling QCD vacuum.

The first approach has a body of successful phenomenology; see review in [7]. On the other hand, as is well known, instanton gas model lacks confinement. As for the second approach, it has gained much popularity in the lattice community, and studies of various projected objects such as Abelian monopoles and center vortices are going on (see, e.g., [8, 9]). It is worth mentioning that despite tremendous recent progress in lattice calculations, the problem of identifying correct gauge-invariant picture of the QCD vacuum is extremely difficult on the lattice as well as in the continuum. For example, the results of [10, 11]--lowering deconfinement critical temperature with rising magnitude of external Abelian chromomagnetic (not chromoelectric!) field--are not easily incorporated to the conventional dual superconductor picture with space-time independent monopole condensate and presumably suggest the formation of nontrivial inhomogeneous condensates. It seems to be correct to conclude that the exact pattern of the Abelian confinement scenario is not yet fixed.

The development of the third approach in a systematic way started in 1987 [12-15] (see, e.g., [16] as example of earlier investigations concerning stochasticity of QCD vacuum). The nontrivial structure of np vacuum can be described by a set of nonlocal gauge-invariant field strength correlators (FCs). The discussion presented below is in the framework of this scenario.

QCD sum rules [17-19] were suggested as an independent approach to np QCD dynamics, not addressing directly confinement mechanism. The key role is played by gluon condensate _G2 , which is defined as np average of the following type: [figure omitted; refer to PDF] In the sum rule framework, it is a universal quantity characterizing QCD vacuum as it is, while it enters power expansion of current-current correlators with channel-dependent coefficients. This is the cornerstone of QCD sum rules ideology.

It is worth stressing from the very beginning that _G2 and other condensates are usually considered as finite physical quantities. Naively in perturbation theory, one would get _G2 ~^a-4 , where a is a space-time ultraviolet cutoff (e.g., lattice spacing). It is always assumed that this "hard'' contribution is somehow subtracted from (1.1) and the remaining finite quantity results from "soft'' np fields, in some analogy with Casimir effect where modification of the vacuum by boundaries of typical size L yields nonzero shift of energy-momentum tensor by the amount of order ^L-4 . In QCD, the role of L is played by dynamical scale ^ΛQCD-1 .

The idea of np gluon condensate has proved to be very fruitful. However the relation of _G2 to confinement is rather tricky. Let us stress that since there can be no local gauge-invariant order parameter for confinement-deconfinement transition, _G2 is not an order parameter, and, in particular, neither perturbative nor np contributions to this quantity vanish in either phase. (Let us repeat that following [17, 18], we define _G2 in (1.1) as purely np object.) On the other hand, one can show (see details in [3]) that the scale of deconfinement temperature is set by the condensate, or, more precisely, by its electric part: _Tc ~^G21/4 . Let us remind the original estimate for the condensate [17, 18] given by _G2 =0.012 Ge^V4 , with large uncertainties. One clearly sees some tension between this "small'' scale and a "large'' mass scale in QCD given by, for example, the mass of the lightest ⁰⁺⁺ glueball, which is about 1.5 GeV.

The existence of two numerically different np scales in strong interaction physics is a common feature one encounters in many models. As an early example, let us mention seminal discussion [20] concerning the difference between confinement scale and chiral symmetry breaking scale. Another example is given by instanton model [21] where the mean instanton size (taken to be about 0.3 Fm) is assumed to be smaller than mean instanton radius (about 1 Fm) and this fact is nothing but the "diluteness'' of the instanton gas in this model. Dual superconductor picture [8, 9] provides more recent example of this sort: for the QCD vacuum understood as the so-called type-II superconductor, the mass of diagonal field--the "photon'' (inverse London penetration depth)--is parametrically larger than the effective Higgs mass. Correspondingly the nonabelian electric field penetrates the vacuum forming Abrikosov fluxes (in dual analogy with the magnetic fields in conventional metal superconductors).

Of immediate interest are recent lattice studies of relevant nonperturbative scales [22]. The authors project each lattice link variable to specific zone in the momentum space bound by some infrared and ultraviolet cuts--the momentum ring--and study dependence of various nonperturbative quantities of physical interest on the location and width of this ring. This analysis demonstrates numerically that the relevant scales are indeed different for different quantities. It is far from obvious, however, what gauge-invariant statements based on these results can be made (since the approach is gauge dependent).

The connections between these models and the approach discussed here were studied in many papers (see reviews in [12-15] and references therein). For example, direct analogy between λ and effective instanton size was analyzed in [23]. The relation between dual superconductor picture and field correlator method was studied in details in [24, 25]. Notice that in both cases the main advantage of field correlator method is its manifest gauge invariance.

The phenomenon of this hierarchy of scales is of prime importance for the stochastic scenario. It was suggested in [12-15]; see also, [26-28] that qualitative physical explanation is the following: QCD vacuum is filled by short-distance correlated nonperturbative gluon fields . In other words, besides gluon condensate _G2 , there is another important dimensionfull parameter characterizing np dynamics of vacuum fields: correlation length λ (also denoted as _Tg in some papers). The gauge-invariant dimensionless product (_G2^λ4)1/2 is about a few percents (see below). This leads to numerous consequences, some of which are discussed in the rest of this paper.

An immediate question to ask is about the origin of this smallness. To answer it is one of the goals of the present paper. The key observation is the following: correlation length λ is nothing but inverse mass of certain gluelump. Gluelumps (for details, we refer the reader to [29, 30]) are rather unusual objects from perturbative field theory point of view. In some sense they are analogous to heavy-light mesons. The latter can be approximated as a bound state of fundamental color charge (light quark) in the static field of fundamental color anticharge (heavy antiquark). In the limit of heavy mass M going to infinity, the difference between the meson mass and M stays constant and the object as a whole is of course color singlet. Analogously, gluelump is a singlet bound state of adjoint color charge "valent gluons'' in the static field of "infinitely heavy'' adjoint color source. Roughly speaking, gluelump is what replaces gluon propagator in gauge-invariant formulation of the theory. Of course there are no elementary heavy particles with adjoint charge in QCD, and in this sense, gluelump cannot be "detected'' as a real particle. Nevertheless, gluelumps happen to be very useful objects in discussions of the most fundamental properties of confining QCD vacuum (see more on that below).

Let us illustrate the appearance of numerical hierarchy of the kind we are discussing on a simple phenomenological example. Consider leading the following ρ -meson Regge trajectory: [figure omitted; refer to PDF] This corresponds to the product [figure omitted; refer to PDF] where σ=0.18 Ge^V2 is the confining string tension. One can say that the mass squared is a factor "2π '' larger than string tension. On the other hand, in stochastic picture one has, parametrically, σ[proportional, variant]^λ2_G2 . Here the smallness of string tension is controlled by two powers of λ . Thus to get self-consistent picture with the desired hierarchy, roughly speaking, one is to make manifest this "1/2π '' factor in ^λ2 . This is the essence of the phenomenon discussed in the present paper.

Let us come back to the stochastic picture. Its crucial feature found in [31-33] is that the lowest, quadratic, nonlocal FC, ..._Fμν (x)_Fλσ (y)..., describes all np dynamics with very good accuracy. It was also shown that a simple exponential form of quadratic correlators found on the lattice [34-38] allows to calculate all properties of lowest mesons, glueballs, hybrids, and baryons, including Regge trajectories, lepton widths, etc.; see reviews in [2, 27, 28] and references therein. Since the correlation length λ entering these exponents is small, potential relativistic quantum-mechanical picture is applicable and all QCD spectrum is defined mostly by string tension σ (which is an integral characteristic of the nonlocal correlator; see below) and not by its exact profile. This is discussed in details in the next section.

However to establish the confinement mechanism unambiguously, one should be able to calculate vacuum field distributions, that is, field correlators, self-consistently. In the long run, it means that one is to demonstrate that it is essential property of QCD vacuum fields ensemble to be characterized by correlators, which support confinement for temperatures T below some critical value _Tc and deconfinement at T>_Tc .

Attempts to achieve this goal have been undertaken in [39, 40]; however, the resulting chain of equations is too complicated to use in practice.

Another step in this direction is done in [41], where FCs are calculated via gluelump Green's functions and self-consistency of this procedure was demonstrated for the first time. These results were further studied and confirmed in [42].

The main aim of this paper is to present this set of equations as a self-consistent mechanism of confinement and to clarify qualitative details of the FC-gluelump connection. In particular we demonstrate how the equivalence of color-magnetic and colorelectric FCs for T=0 ensures the Gromes relation [43, 44]. We also show, that self-consistency condition for FC as gluelump Green's function allows to connect the mass scale to _αs , and in this way, to express _ΛQCD via string tension σ .

Lattice computations play important role as a source of independent knowledge about vacuum field distributions. Recently a consistency check of this picture was done on the lattice [45] by measurement of spin-dependent potentials. The resulting FCs are calculated and compared with gluelump predictions in [46] demonstrating good agreement with analytic results; in particular small vacuum correlation length λ~0.1 Fm is shown to correspond to large gluelump mass _M0 [approximate]2 GeV.

The paper is organized as follows. In the next section, we remind the basic expressions for Green's functions of qq... and gg systems in terms of Wegner-Wilson loops, and qualitative picture of FCs dynamics is discussed. Section 3 is devoted to the expressions of spin-dependent potentials in terms of FC. We also shortly discuss the lattice computations of FC. In Section 4 FC are expressed in terms of gluelump Green's functions, while in Section 5 the self-consistency of resulting relations is studied. Our conclusions are presented in Section 6.

2. Wegner-Wilson Loop, Field Correlators, and Green's Functions

Since our aim is to study confinement for quarks (both light and heavy) and also for massless gluons, we start with the most general Green's functions for these objects in a proper physical background using Fock-Feynman-Schwinger representation. The formalism is built in such a way that gauge invariance is manifest at all steps. A reader familiar with this technique can skip this section and go directly to Section 3.

First, we consider the case of mesons made of quarks while gluon bound states are discussed below. For quarks, one forms the initial and final state operators [figure omitted; refer to PDF] where ^ψ[dagger] , ψ are quark operators, Γ is a product of Dirac matrices, that is, 1, _γμ ,_γ5 ,(_γμγ5 ),..., and Φ(x,y)=Pexp (ig^∫yx_Aμ d_zμ ) is parallel transporter (also known as Schwinger line or phase factor). The meson Green's function can be written (in quenched case) as [figure omitted; refer to PDF] and the quark Green's function _Sq is given in Euclidean space time (see [47, 48] and references therein) as [figure omitted; refer to PDF] where K is kinetic energy term as [figure omitted; refer to PDF] where m is the pole mass of quark, and quark trajectories _zμ (τ) with the end points x and y are being integrated over in (Dz_)xy .

The factor _Pσ (x,y,s) in (2.3) is generated by the quark spin (color-magnetic moment) and is equal to [figure omitted; refer to PDF] where _σμν =(1/4i)(_γμγν -_γνγμ ), and _PF (_PA ) in (2.5) and (2.3) are, respectively, ordering operators of matrices _Fμν (_Aμ ) along the path _zμ (τ) .

With (2.2) and (2.3), one easily gets [figure omitted; refer to PDF] where Tr means trace operation both in Dirac and color indices, while [figure omitted; refer to PDF] In (2.7) the closed contour C(x,y) is formed by the trajectories of quark _zμ (τ) and antiquark ^{zν[variant prime]} (^{τ[variant prime]} ) , and the ordering _PA and _PF in _Pσ , ^{Pσ[variant prime]} is universal, that is, _Wσ (x,y) is the Wegner-Wilson loop (W-loop) for spinor particle.

The factors (m-D...) and (^{m[variant prime]} -^{D...[variant prime]} ) in (2.6) need a special treatment when being averaged over fields. As shown in Appendix 1 of [49], one can use the simple replacement [figure omitted; refer to PDF]

The representation (2.6) is exact in the limit _Nc [arrow right]∞ , when internal quark loops can be neglected, and is a functional of gluonic fields _Aμ , _Fμν , which contains both perturbative and np contributions, not specified at this level.

The next step is averaging over gluon fields, which yields the physical Green's function _Gqq... : [figure omitted; refer to PDF] The averaging is done with the usual Euclidean weight exp (-Action/...) , containing all necessary gauge-fixing and ghost terms.

To proceed, it is convenient to use nonabelian Stokes theorem (see [2] for references and discussion) for the first factor on the r.h.s. in (2.7) (corresponding to the W-loop for scalar particle) and to rewrite it as an integral over the surface S spanned by the contour C=C(x,y) : [figure omitted; refer to PDF] Here F(_ui )d_si =Φ(_x0 ,_ui )^Fμνa (_ui )^ta Φ(_ui ,_x0 )d_sμν (_ui ) , and _ui , _x0 are the points on the surface S bound by the contour C=∂S . The double brackets ............... stay for irreducible correlators proportional to the unit matrix in the color space (and therefore, only spacial ordering _Px enters (2.10)). Since (2.10) is gauge-invariant, one can make use of generalized contour gauge [50-52], which is defined by the condition Φ(_x0 ,_ui )≡1 . Notice that throughout the paper we normalize trace over color indices as Tr 1=1 in any given representation.

Since (2.10) is an identity, the r.h.s. does not depend on the choice of the surface, which is integrated over in d_sμν (u) . On the other hand, it is clear that each irreducible n -point correlator ......F...F...... integrated over S (these are the functions ^Δ(n) [S] in (2.10)) depends, in general case, on the choice one has made for S . Therefore it is natural to ask the following question: is there any hierarchy of ^Δ(n) [S] as functions of n for a given surface S ?

To get the physical idea behind this question, let us take the limit of small contour C . In this case one has for the np part of the W-loop in fundamental representation the following: [figure omitted; refer to PDF] where S is the minimal area inside the loop C . The np short-distance dynamics is governed by the vacuum gluon condensate _G2 ; see (1.1). Higher condensates and other possible vacuum averages of local operators, in line with Wilson expansion, have been introduced and phenomenologically estimated in [19]; for a review see [53].

Suppose now that the size of W-loop is not small, that is, it is larger than some typical dynamical scale λ=...AA;(1 GeV) to be specified below. Let us now ask the following question: if one still wants to expand the loop formally in terms of local condensates, what would the structure of such expansion be? It is easy to see from (2.10) that there are two subseries in this expansion. The first one is just an expansion in n , that is, in the number of field strength operators. It is just [figure omitted; refer to PDF] The second series is an expansion of each ^Δ(n) [S] in terms of local condensates, that is, expansion in powers of derivatives. Taking for simplicity the lowest n=2 term, this expansion reads [figure omitted; refer to PDF]

It is worth mentioning (but usually skipped in QCD sum rules analysis) that the latter series contain ambiguity related to the choice of contours in parallel transporters. The expression (2.13) is written for the simplest choice of straight contour connecting the points x and y (and not x and _x0 or y and _x0 ). This is the choice we adopt in what follows. In general case each condensate in (2.13) is multiplied by some contour-dependent function of x , y, and _x0 . It is legitimate to ask whether this approximation, that is, replacement of [figure omitted; refer to PDF] by [figure omitted; refer to PDF] affects our final results. Since the only difference between the above two expressions is profiles of the transporters, this is the question about contour dependence of FCs. We will not discuss this question in details in the present paper and refer an interested reader to [12-15]. Here we only mention that the result (weak dependence on transporters' choice) is very natural in Abelian dominance picture since for Abelian fields the transporters exactly cancel.

Let us stress that both expansions (2.12) and (2.13) are formal since the loop is assumed to be large. Moreover, from dimensional point of view, the terms of these two expansions mix with each other. For example, one has nonzero v.e.v. of two operators of dimension eight: [figure omitted; refer to PDF] with no a priori reason (one could have such reasons if there is some special kinematics in the problem, but this is not the case.) to drop any of them. Moreover, even if one assumes that all these condensates of high orders can be self-consistently defined analogously to _G2 , it remains to be seen whether the whole series converges or not.

Here the physical picture behind FC method comes into play. It is assumed (and confirmed a posteriori in several independent ways) that if the surface S is chosen as the minimal one, the dominant contribution to (2.10) results from v.e.v.s of the operators F^Dk F, and this subseries can be summed up as (2.13). In this language, nonlocal two-point FC can be understood as a representation for the sum of infinite sequence of local terms (2.13).

Physically the dominance of ...F^Dn F... terms (and hence, of two-point FC) corresponds to the fact that the correlation length is small yielding a small expansion parameter ξ=F...^λ2 <<1 (F... is average modulus of np vacuum fields), or, in other words, that typical inverse correlation length ^λ-1 characterizing ensemble of QCD vacuum fields is parametrically larger than fields themselves. The dimensionless parameter ξ in terms of condensate is given by (_G2^λ4)1/2 , and it is indeed small: of the order of 0.05 according to lattice estimates. One can say that often repeated statement "QCD vacuum is filled by strong and strongly fluctuating chromoelectric and chromomagnetic fields'' is only partly correct, namely, in reality the fields are not so strong as those of ^λ-1 . Technically one can say that we sum up the leading subseries in (2.10) in the same spirit as one does in covariant perturbation theory [54] for weak but strongly varying potentials.

Thus, for understanding of confining properties of QCD vacuum dynamics, not only scale of np fields given by _G2 but also another quantity--the vacuum correlation length λ -- which defines the nonlocality of gluonic excitations is important. It is worth repeating here that physically λ encodes information about properties of the series (2.13) and has the same theoretical status as that of the corresponding np condensates in the following sense: knowledge of all terms in (2.13) would allow to reconstruct λ . Of course, it is impossible in practice and one has to use other methods to study large distance asymptotics of FC. On more phenomenologic level, it was discussed in [55-59], while rigorous definition was given in the framework of the FC method (often referred to as Stochastic Vacuum Model.) (see [12-15, 26, 60] for reviews).

It is important to find also the world-line representation for gluon Green's function, which is done in the framework of background perturbation theory (see [61-65] for details). In this way one finds (where a , b are adjoint color indices and dash sign denotes adjoint representation) that [figure omitted; refer to PDF] where [figure omitted; refer to PDF] All the above reasoning about asymptotic expansions of the corresponding W-loops is valid here as well.

Thus the central role in the discussed method is played by quadratic (Gaussian) FC of the form [figure omitted; refer to PDF] where _Fμν is the field strength and Φ(x,y) is the parallel transporter. Correlation lengths _λi for different channels are defined in terms of asymptotics of (2.19) at large distances: exp (-|x-y|/_λi ) . The physical role of _λi is very important since it distinguishes two regimes: one expects validity of potential-type approach describing the structure of hadrons of spatial size R and at temporal scale _Tq for _λi <<R,_Tq , while in the opposite case, when _λi >>R,_Tq the description in terms of spatially homogeneous condensates can be applied.

At zero temperature, the O(4) invariance of Euclidean space-time holds and FC (2.19) is represented through two scalar functions D(z),_D1 (z) (where z≡x-y ) as follows(all treatment in this paper as well as averaging over vacuum fields is done in the Euclidean space. Notice that only after all Green's functions are computed, analytic continuation to Minkowskii space-time can be accomplished): [figure omitted; refer to PDF]

One has to distinguish from the very beginning perturbative and np parts of the correlators D(z) , _D1 (z) . Beginning with the former, one easily finds at tree level that [figure omitted; refer to PDF] where _C2 (f) is fundamental Casimir _C2 (f)=(^Nc2 -1)/2_Nc . At higher orders, situation becomes more complicated. Namely, one has at n -loop order [figure omitted; refer to PDF] where the gauge-invariant functions ^G(n) (z) , ^G1(n) (z) have the following general structure: [figure omitted; refer to PDF] where _cn is numerical coefficient and μ is renormalization scale, and we have omitted subleading logarithms and constant terms in the r.h.s.; explicit expressions for the case n=1 can be found in [66, 67].

Naively one could take perturbative functions ^Dp,n (z) , ^D1p,n (z) written above and use them for computations of W-loop or static potentials in Gaussian approximation. This, however, would be incorrect. The reason is that at any given order in perturbation expansion over _αs one should take into account perturbative terms of the given orders coming from all FCs, and not only from the Gaussian one. For example, at one loop level, the function ^Dp,1 (z) is nonzero which naively would correspond to area law at perturbative level. Certainly this cannot be the case since self-consistent renormalization program for W-loops is known not to admit any terms of this sort [68, 69]. Technically in FC language, the correct result is restored by cancelation of contributions proportional to ^Dp,1 (z) to all observables by the terms coming from triple FC (see details in [70]).

In view of this general property of cluster expansion, it is more natural just to take relation ^Dp,n (z)≡0 as valid at arbitrary n , having in mind that proper number of terms from (2.12) has to be included to get this result. Thus we assume that the following decomposition takes place: [figure omitted; refer to PDF] and D(z) has a smooth limit when z[arrow right]0 . As for _D1 (z), its general asymptotic at z[arrow right]0 reads as [figure omitted; refer to PDF] where c and _a2 weakly (logarithmically) depend on z . Notice that by _D1 (0), we always understand np "condensate'' part in what follows: [figure omitted; refer to PDF]

The term 1/^z2 deserves special consideration. One may argue that at large distances the simple additivity of perturbative and np contributions to _VQQ... potential is violated [71]. It is interesting whether this additivity holds at small, distances and, in particular, mixed condensate _a2 was suggested in [72] and argued to be phenomenologically desirable.

As we demonstrate below, our approach is self-consistent with _a2 =0 , that is, when perturbative-np additivity holds at small distances and there is no dimension-two condensate. However, at intermediate distances, a complicated interrelation occurs, which does not exclude (2.25) being relevant in this regime.

We proceed with general analysis of (2.19). As proved in [12-15], D(z) , _D1 (z) do not depend on relative orientation of plane (μν) , plane (λσ) , and vector _zα . This orientation can be of the following 3 types: (a) planes (μν) and (λσ) are perpendicular to the vector _zα ; (b) planes (μν) and (λσ) are parallel or intersecting along one direction and _zα lies in one of planes; (c) planes (μν) and (λσ) are perpendicular and _zα lies necessarily in one of them.

It is easy to understand that this classification is O(4) invariant, and therefore, one can assign indices a,b,c to _Dμν,λσ , resulting in general case in three different functions. Physically speaking, case (a) refers to, for example, color magnetic fields _Fik , _Flm with _zα in the 4th direction, and the corresponding FCs will be denoted as _{D[perpendicular]} (z) . Case (b) refers to, for example, the color electric fields _Ei (x) , _Ek (y) , connected by the same temporal links along 4th axis, and the corresponding FCs are denoted as _D|| (z) . Finally, in case (c) only _D1 part survives in (2.20) and it will be given by the subscript EH , ^D1EH .

In case of zero temperature, when O(4) invariance holds, these three functions _{D[perpendicular]} (z) , _D|| (z) , and ^D1EH (z) are easily expressed via D(z) , _D1 (z):

[figure omitted; refer to PDF] (note that just _{D[perpendicular]} , _D|| were measured on the lattice in [34-38]). For nonzero temperature, the correlators of color-electric and color-magnetic fields can be different and in general there are five FCs: ^DE (z) , ^D1E (z) , ^DH (z) , ^D1H (z) , and ^D1EH (z) . As a result, one obtains the full set of five independent quadratic FCs which can be defined as follows: [figure omitted; refer to PDF] It is interesting that the structure (2.28) survives for nonzero temperature, where O(4) invariance is violated: it tells that (2.28) give the most general form of the FCs. Note that at zero temperature, ^DH =^DE , ^D1H =^D1E at coinciding point; the O(4) invariance requires that colorelectric and colormagnetic condensates coincide: ^g2 ...Tr _Hi (x)_Hi (x)...=^g2 ...Tr _Ej (x)_Ej (x)...; hence, [figure omitted; refer to PDF] Here we assume that the np part of all FCs is finite.

At z≠0 for zero temperature, FCs can be expressed through _{D[perpendicular]} (z) , _D|| (z) which do not depend on whether _Ei or _Hk enters in them, since, for example, ^D||E can be transformed into ^D||H by an action of O(4) group elements (the same for ^{D[perpendicular]E} , ^{D[perpendicular]H} ). From this, one can deduce that both _D1 and D do not depend on subscripts E , H for zero temperature and ^D1EH coincides with _D1 .

Below we will illustrate this coincidence by concrete calculations of ^DE , ^DH , and so forth, through the gluelump Green's functions and show that the corresponding correlation lengths satisfy the relations (with obvious notations) [figure omitted; refer to PDF] As shown below, all correlation lengths _λj appear to be just inverse masses of the corresponding gluelumps _λj =1/_Mj .

Having analytic expressions for FCs, one might ask that how to check them versus experimental and lattice data. On experimental side in hadron spectroscopy, one measures masses and transition matrix elements, which are defined by dynamical equation, and the latter can be used of potential type due to smallness of _λj . In static potential, only integrals over distance enter and the spin-independent static potential can be written as [12-15, 73, 74] [figure omitted; refer to PDF]

At this point one should define how perturbative and np contributions combine in ^DE , ^D1E , and this analysis was done in [70]. Making use of (2.21) together with definition of the string tension [12-15], one yeilds [figure omitted; refer to PDF] One obtains from (2.31) the standard form of the static potential for N=3 at distances r>>^λE as [figure omitted; refer to PDF] As argued below, ^λE [approximate]0.1 Fm, so that the form (2.33) is applicable in the whole range of distances r>0.1 Fm provided that asymptotic freedom is taken into account in _αs (r) in (2.33).

At the smallest values of r , r[<, ~]^λE , one would obtain a softening of confining term, σr[arrow right]c^r2 [12-15], which is not seen in accurate lattice data at r[> ~]0.2 Fm [75, 76], imposing a stringent limit on the value of ^λE ; ^λE [<, ~]0.1 Fm; see discussion in [76].

3. Spin-Dependent Potentials and FC

We now turn to the spin-dependent interactions to demonstrate that FC can be extracted from them [73, 74, 77-79]. ^{V1[variant prime]} (r),^{V2[variant prime]} (r),_V3 (r),_V4 (r) , plus static term ^{VQQ...[variant prime]} (r) contain in integral form all five FCs: ^DE ,^D1E ,^DH ,^D1H ,^D1HE , and one can extract the properties of the latter from the spin-dependent potentials. This procedure is used recently for comparison of analytic predictions [46] with the lattice data [45].

To express spin-dependent potentials in terms of FC, one can start with the W-loop (2.7), entering as a kernel in the qq... Green's function, where quarks q , q... can be light or heavy. Since both kernels _Pσ , ^{Pσ[variant prime]} contain the matrix _σμνFμν , the terms of spin-spin and spin-orbit types appear. As before, we keep Gaussian approximation, and in this case spin-dependent potentials can be computed not only for heavy but also for light quarks. They have the following Eichten-Feinberg form [80]: [figure omitted; refer to PDF]

Spin-spin interaction appears in _Wσ in (2.7) which can be written as [figure omitted; refer to PDF] where _z1,2 =z(_τ1,2 ), and spin-orbit interaction arises in (2.19) from the products ..._σμνFμν d_sρλFρλ ... .

It is clear that the resulting interaction will be of matrix form (2×2)×(2×2) (not accounting for Pauli matrices). If one keeps only diagonal terms in _σμνFμν (as the leading terms for large _μi [approximate]M ) then one can write for the spin-dependent potentials the representation of the Eichten-Feinberg form (3.1).

At this point one should note that the term with d[straight epsilon]/dR in (3.1) was obtained from the diagonal part of the matrix (m-D...)_σμνFμν , namely, as product i_σkDk ·_σiEi ; see [73, 74, 77, 78] for details of derivation, while all other potentials _Vi , i=1,2,3,4 are proportional to FCs ..._Hi Φ_Hk Φ... . One can relate FCs of colorelectric and colormagnetic fields to ^DE , ^D1E , ^DH , ^D1H defined by (2.28) with the following result: [figure omitted; refer to PDF]

One can check that Gromes relation [43, 44] acquires the form [79]

[figure omitted; refer to PDF]

For T=0 , when ^DE =^DH , ^D1E =^D1H , the Gromes relations are satisfied identically; however, for T>0 electric and magnetic correlators are certainly different and Gromes relation is violated, as one could tell beforehand; since for T>0, the Euclidean O(4) invariance is violated.

4. Gluelumps and FC

In this section we establish a connection of FC D(z) and _D1 (z) with the gluelump Green's functions.

To proceed, one can use the background field formalism [61-65, 81-83], where the notions of valence gluon field _aμ and background field _Bμ are introduced, so that total gluonic field _Aμ is written as [figure omitted; refer to PDF] The main idea we are going to adopt here is suggested in the second paper of [39, 40]. To be self-contained, let us explain it here in simple form. First, we single out some color index a and fix it at a given number. Then in color components we have, by definition, that [figure omitted; refer to PDF] and there is summation neither over a in the first term nor over c in the second. As a result, the integration measure factorizes [figure omitted; refer to PDF] so one first averages over the fields ^Bμc , and after that one integrates over the fields ^aμa . The essential point is that the integration over ...9F;^Bμc provides colorless adjoint string attached to the gluon ^aμa , which keeps the color index "a '' unchanged in the course of propagation of gluon a in background made of B . This is another way of saying that gluon field of each particular color enters QCD Lagrangian quadratically. Despite the fact that such separation of the measure is a kind of trick, the formation of adjoint string is the basic physical mechanism behind this background technic, and it is related to the properties of gluon ensemble: even for _Nc =3, one has one color degree of freedom ^aμa and 7 fields ^Bμc ; also confining string is a colorless object, and therefore, the singled-out adjoint index "a '' can be preserved during interaction process of the valence gluon ^aμa with the background (^Bμc ) . These remarks to be used in what follows make explicit the notions of the valence gluon and background field.

We assume the background Feynman gauge _Dμaμ =0 and assign field transformations as follows: [figure omitted; refer to PDF] As a result, the parallel transporter Φ(x,y)=Pexp ig^∫yx_Bμ (z)dz keeps its transformation property, and every insertion of _aμ between Φ transforms gauge covariantly. In what follows we will assume that Φ is made entirely of the field _Bμ . Note that in the original cluster expansion of the W-loop W(B+a) in FC [81-83] [figure omitted; refer to PDF] it does not matter whether (_Bμ +_aμ ) or _Bμ enters Φ' s since those factors cancel in the sum. Thus we assume from the beginning that Φ' s are not renormalized, since such renormalization is unphysical anyway according to what have been said above.

As for renormalization having physical meaning, it is known [68, 69] to reduce to the charge renormalization and renormalization of perimeter divergences on the contour C . For static quarks, the latter reduces to the mass renormalization.

In background field formalism, one has an important simplification: the combination g_Bμ is renorminvariant [61-65, 68, 69], and so is any expression made of _Bμ only. This point is important for comparison of FC with those of the lattice correlators in [45, 46]. Renormalization constants for field strengths _Zb used there account for finite size of plaquettes and they are similar to the lattice tadpole terms.

Now we turn to the analytic calculation of FC in terms of gluelump Green's function. To this end, we insert (4.1) into (2.19) and have for _Fμν (x) the following: [figure omitted; refer to PDF] Here, the term ^Fμν(B) contains only the field ^Bμb . It is clear that when one averages over field ^aμa and sums finally over all color indices a , one actually exploits all the fields with color indices from ^Fμν(B) , so that the term ^Fμν(B) can be omitted, if summing over all indices a is presumed to be done at the end of calculation.

As a result, _Dμν,λσ can be written as [figure omitted; refer to PDF] where the superscripts 0, 1, 2 refer to the power of g coming from the term ig[_aμ ,_aν ] .

We can address now an important question about the relation between FCs and gluelump Green's functions. We begin with 1-gluon gluelump, whose Green's function reads [figure omitted; refer to PDF] According to (4.4), ^Gμν(1g) (x,y) is a gauge-invariant function.

As shown in [41], the first term in (4.7) is connected to the functions ^D1E , ^D1H and it can be written as follows: [figure omitted; refer to PDF] where ^{Δμν,λσ(0)} contains contribution of higher FCs, which we systematically discard.

On the other hand, one can find ^Gμν(1g) (x,y) from the expression [47, 48] written as [figure omitted; refer to PDF] where the spin-dependent W-loop is [figure omitted; refer to PDF] and the gluon spin factor is exp F≡exp (2ig^∫0s dτ_F...B (z(τ))) with _F...B made of the background field _Bμ only.

Analogous expression can be constructed for Green's function of 2-gluon gluelump. It is given by the following expression : [figure omitted; refer to PDF] At small, distances ^G(2gl) (x,y) is dominated by the perturbative expansion terms [figure omitted; refer to PDF] however, as we already discussed in details, all these perturbative terms are canceled by those from higher FCs (triple, quartic, etc.); therefore, expansion in fact starts with np terms of dimension four.

To identify the np contribution to ^G(2gl) (z), we rewrite it as follows: [figure omitted; refer to PDF] where the two-gluon gluelump W-loop _WΣ (_C1 ,_C2 ) is depicted in Figure 6 and can be written as (in the Gaussian approximation) [figure omitted; refer to PDF] and the total surface S consists of 3 pieces, as shown in Figure 4, as follows: [figure omitted; refer to PDF] where (F...) has Lorentz tensor and adjoint color indices ^Fijab and lives on gluon trajectories, that is, on the boundaries of _Si . In full analogy with (4.9), we have that [figure omitted; refer to PDF]

The crucial point is that Green's functions (4.9), (4.14) can be calculated in terms of the same FCs D(z) , _D1 (z) . Indeed one has that [figure omitted; refer to PDF] where the index w stays for 1-gluon or 2-gluon gluelump Hamiltonians. The latter are expressed via the same FCs D(z) , _D1 (z) (see [84-86] and references therein): [figure omitted; refer to PDF] where the last term _Hstring [D,ν] depends on D(z) via the adjoint string tension [figure omitted; refer to PDF] and Hamiltonian depends on einbein fields μ and ν . The term corresponding to perimeter Coulomb-like interaction Δ_HCoul [_D1 ] depends on _D1 (z) (compare with (2.31), (3.1)).

The self-consistent regimes correspond to different asymptotics of the solutions to these equations. In Coulomb phase of a gauge theory, both _H1g , _H2g exhibit no mass gap, that is, large-z asymptotic of (4.18) is power like. The function D(z) vanishes in this phase and W-loop obeys perimeter law. In the confinement phase realized in Yang-Mills theory at low temperatures, a typical large-z pattern is given by [figure omitted; refer to PDF] that is, there is a mass gap for both _H1g , _H2g .

Confining solutions are characterized by W-loops obeying area law. In other words, Hamiltonians expressed in terms of interaction kernels depending on D(z) , _D1 (z) exhibit mass gap if these kernels are confining. On the other hand, the same mass gap plays a role of inverse correlation length of the vacuum. This should be compared with well-known mean-field technique. In our case the role of mean-field is played by quadratic FC, which develops nontrivial Gaussian term D(z) . Notice that the exponential form of its large distance asymptotics exp (-|z|/λ) (and not, let's say, exp (-^z2 /^λ2 ) ) is dictated by spectral expansion of the corresponding Green's function at large distances.

Full solution of the above equations is a formidable task not addressed by us here. Instead, as a necessary prerequisite, we check below different asymptotic regimes and demonstrate self-consistency of the whole picture.

We begin with small distance region. The np part of contribution to _D1 (z) at small distances comes from two possible sources: the area law term (first exponent and the exp F term in (4.11)). Both contributions are depicted in Figures 1 and 2, respectively. We will disregard the term exp F in this case, since for the one-gluon gluelump it does not produce hyperfine interaction and only gives rise to the np shift of the gluon mass, which anyhow is eliminated by the renormalization (see the appendix for more details). This is in contrast to the two-gluon gluelump Green's function generating ^DE , ^DH , where the hyperfine interaction between two gluons is dominating at small distances.

Figure 1: [figure omitted; refer to PDF]

Figure 2: [figure omitted; refer to PDF]

As a result, (4.8) without the exp F term yields [figure omitted; refer to PDF] It is remarkable that the sign of the np correction is positive.

At large distances, one can use the gluelump Hamiltonian for one-gluon gluelump from [30] to derive the asymptotics [41] [figure omitted; refer to PDF] where ^M0(1) =(1.2÷1.4) GeV for _σf =0.18 GeV² [29, 30].

We now turn to the FC D(z) as was studied in [41]. The relation (4.17) connects D(z) to the two-gluon gluelump Green's function, studied in [30] at large distances. Here we need its small-z behavior and we will write it in the form [figure omitted; refer to PDF] where ^Gp(2gl) (z) contains purely perturbative contributions which are subtracted by higher-order FCs, while ^Gnp(2gl) (z) contains np and possible perturbative-np interference terms. We are interested in the contribution of the FC ...FF... to ^G(2gl) (z) , when z tends to zero.

One can envisage three types of contributions as follows.

(a) The first one is due to product of surface elements d_sμν d_sλσ , which gives for small surface the factor exp (-^g2 ...FF...^Si2 /24_Nc ) , similar to the situation discussed for _D1 . This is depicted in Figure 4.

(b) The second one is contribution of the type d_τ1 d_τ2 ...FF... , where dτ and d_τ2 belong to different gluon trajectories, 1 and 2, as depicted in in Figure 5. This is the hyperfine gluon-gluon interaction, taken into account in [30] in the course of the gluelump mass calculations (regime of large z ); however, in that case mostly the perturbative part of ...FF... contributes (due to _D1 ). Here we keep in ...FF... only the NP part and consider the case of small z .

(c) Again the third is the term d^τi d_τj , but now i=j . This is actually a part of the gluon selfenergy correction, which should be renormalized to zero, when all (divergent) perturbative contributions are added. As we argued in Appendix of [41], we disregard this contribution here, as well as in the case of _D1 (z) (one-gluon gluelump case). Now we treat both contributions (a) and (b).

(a) In line with the treatment of _D1 , (4.22), one can write the term representing two-gluon gluelump as two nearby one-gluon gluelumps which yield for D(z) the following: [figure omitted; refer to PDF]

(b) In this case one should consider the diagram given in Figure 5. which yields the answer (for details see the appendix) [figure omitted; refer to PDF] which contributes to D(z) as [figure omitted; refer to PDF] where at small z , h(z)[approximate](D(_λ0 )/64^π4 )^{log 2} (_λ0 e/z) and _λ0 is of the order of correlation length λ .

At large distances, one uses the two-gluon gluelump Hamiltonian as in [29] and finds the corresponding spectrum and wave functions; see [29] and Appendix 5 of [41] for details. As a result, one obtains in this approximation [figure omitted; refer to PDF] where ^M0(2) is the lowest two-gluon gluelump mass found in [30] to be about ^M0(2) =(2.5÷2.6) GeV.

We will discuss the resulting properties of D(x) and _D1 (x) in the next section.

5. Discussion of Consistency

We start with the check consistency for D(z) . As is shown above, D(z) has the following behavior at small z: [figure omitted; refer to PDF] Since _αs (μ(z))~2π/_β0 log (Λz^)-1 , the first term is subleading at z[arrow right]0 and the last term on the r.h.s. tends to a constant [figure omitted; refer to PDF] From (5.2) one can infer that D(0)[approximate]0.15D(_λ0 ) for _Nc =3, where _λ0 [> ~]^λE . So D(z) is an increasing function of z at small z , z[<, ~]_λ0 , and for z>>λ, one observes exponential falloff. The qualitative picture illustrating this solution for D(z) is shown in Figure 5.

This pattern may solve qualitatively the contradiction between the values of D(0) estimated from the string tension _Dσ (0)≈ σ/π^λ2 [approximate]0.35 GeV⁴ and the value obtained in naive way from the gluon condensate _DG2 (0)=(^π2 /18)_G2 [approximate](0.007÷0.012) GeV⁴ . One can see that _Dσ (0)[approximate](30÷54)_DG2 (0) . This seems to be a reasonable explanation of the mismatch discussed in Section 1.

As shown in [41], the large distance exponential behavior is self-consistent, since (assuming that it persists for all z , while small-z region contributes very little) from the equality σ=π^λ2_Dσ (0) , comparing with (4.28), one obtains [figure omitted; refer to PDF] where in _αs (μ) the scale μ corresponds roughly to the gluelump average momentum (inverse radius) _μ0 [approximate]1 GeV. Thus (5.3) yields _αs (_μ0 )[approximate]0.4 which is in reasonable agreement with _αs from other systems [87].

We end this section with discussion of three points as follows.

(1 ) D(z) and _D1 (z) have been obtained here in the leading approximation, when gluelumps of minimal number of gluons contribute 2 for D(z) and 1 for _D1 (z) . In the higher orders of O(_αs ), one has an expansion of the type [figure omitted; refer to PDF] Hence the asymptotic behavior for D(z) will contain exponent of ^M0(1) |z| too, but with a small preexponent coefficient.

(2 ) The behavior of D(z) , ^D1(np) (z) at small z is defined by NP terms of dimension four, which are condensate _G2 as in (2.11) and the similar term from the expansion of exp F , namely, [figure omitted; refer to PDF] therefore, one does not encounter mixed terms like O(^m2 /^z2 ); however, if one assumes that expansion of ...F(0)F(z)... starts with terms of this sort as it was suggested in [72], then one will have a self-consistent condition for the coefficient in front of this term.

(3 ) To study the difference between ^D1E , ^D1H at T≤_Tc , one should look at (4.9) and compare the situation, when μ=λ=4 , ν=σ=i and take z=x-y along the 4th axis (for ^D1E ), while for ^D1H one takes μ=λ=i , ν=σ=k, and the same z . One can see that in both cases one ends up with the one-gluon gluelump Green's function _Gμν (4.10) which in the lowest (Gaussian) approximation is the same _Gμν =_δμν f(z) , and f(z) is the standard lowest mass (and lowest angular momentum) Green's function.

Hence ^D1E =^D1H in this approximation, and for T=0, this is an exact relation, as discussed above. Therefore,

[figure omitted; refer to PDF] The value ^M0(1) in (5.6) is taken from the calculations in [30]. The same is true for ^DE , ^DH , as it is seen from (4.12), where ^G(2gl) corresponds to two-gluon subsystem angular momentum L=0 independently of μν , λσ .

Hence one obtains [figure omitted; refer to PDF] where the value ^M0(2) [approximate]2.5 GeV is taken from [30].

6. Conclusions

We have derived, following the method of [61-65], the expressions for FC D(z) , _D1 (z) in terms of gluelump Green's functions. This is done in Gaussian approximation. The latter are calculated using Hamiltonian where np dynamics is given by ^Dnp (z) , ^D1np (z) . In this way one obtains self-coupled equations for these functions, which allow two types of solutions: (1 ) ^Dnp (z)=0 , ^D1np (z)=0 , that is, no np effects at all; (2 ) ^Dnp (z) , ^D1np (z) are nonzero and defined by the only scale, which should be given in QCD, for example, string tension σ or _ΛQCD . All other quantities are defined in terms of these basic ones. We have checked consistency of self-coupled equations at large and small distances and found that to the order ...AA;(_αs ) no mixed perturbative-np terms appear. The function _D1 (z) can be represented as a sum of perturbative and np terms, while D(z) contains only np contributions.

We have found a possible way to explain the discrepancy between the average values of field strength taken from _G2 and from σ by showing that D(z) has a local minimum at z=0 and grows at z~λ . Small value of λ and large value of the gluelump mass _Mgl =1/λ[approximate]2.5 GeV explain the lattice data for λ[approximate]0.1 Fm. Thus the present paper argues that relevant degrees of freedom ensuring confinement are gluelumps, described self-consistently in the language of FCs.

Acknowledgments

This work is supported by the grant for support of scientific schools NS-4961.2008.2. The work of the second author is also supported by INTAS-CERN fellowship 06-1000014-6576.

References

[1] Yu. A. Simonov, "The confinement," Physics-Uspekhi , vol. 39, no. 4, pp. 313-336, 1996.

[2] D. S. Kuz'menko, V. I. Shevchenko, Yu. A. Simonov, "Vacuum, confinement, and QCD strings in the vacuum correlator method," Physics-Uspekhi , vol. 47, no. 1, pp. 1-15, 2004.

[3] Yu. A. Simonov, "Dyons in QCD: confinement and chiral symmetry breaking," http://arxiv.org/abs/hep-ph/9509403v2 in Lecture at the International School of Physics "Enrico Fermi"; Course 80: Selected Topics in Nonperturbative QCD, pp. 339-364, Varenna, Italy, Jun-Jul 1995.

[4] A. Di Giacomo, "Confinement of color: open problems and perspectives," in Proceedings of the 17th Internatinal Conference on Particles and Nuclei, vol. 842, of AIP Conference Proceedings, pp. 216-221, Santa Fe, NM, USA, October 2006.

[5] G. 't Hooft, A. Zichichi, "On the phase transition towards permanent quark confinement," Nuclear Physics B. , vol. 138, pp. 1-43, 1978.

[6] S. Mandelstam, "Vortices and quark confinement in non-abelian gauge theories," Physics Letters B , vol. 53, no. 5, pp. 476-478, 1975.

[7] T. Schäfer, E. V. Shuryak, "Instantons in QCD," Reviews of Modern Physics , vol. 70, no. 2, pp. 323-425, 1998.

[8] K. Amemiya, H. Suganuma, "Off-diagonal gluon mass generation and infrared Abelian dominance in the maximally Abelian gauge in lattice QCD," Physical Review D , vol. 60, no. 11, 1999., [email protected]

[9] V. G. Bornyakov, M. N. Chernodub, F. V. Gubarev, S. M. Morozov, M. I. Polikarpov, "Abelian dominance and gluon propagators in the maximally Abelian gauge of SU(2) lattice gauge theory," Physics Letters B , vol. 559, no. 3-4, pp. 214-222, 2003.

[10] P. Cea, L. Cosmai, "Abelian chromomagnetic fields and confinement," Journal of High Energy Physics , vol. 7, no. 2, pp. 655-665, 2003., [email protected]; [email protected]

[11] L. Campanelli, M. Ruggieri, "Probing the QCD vacuum with an Abelian chromomagnetic field: a study within an effective model," Physical Review D , vol. 80, no. 3, 2009., [email protected]; [email protected]

[12] H. G. Dosch, "Gluon condensate and effective linear potential," Physics Letters B , vol. 190, no. 1-2, pp. 177-181, 1987.

[13] H. G. Dosch, Yu. A. Simonov, "The area law of the Wilson loop and vacuum field correlators," Physics Letters B , vol. 205, no. 2-3, pp. 339-344, 1988.

[14] Yu. A. Simonov, "Vacuum background fields in QCD as a source of confinement," Nuclear Physics B , vol. 307, no. 3, pp. 512-530, 1988.

[15] A. Di Giacomo, H. G. Dosch, V. I. Shevchenko, Yu. A. Simonov, "Field correlators in QCD. Theory and applications," Physics Reports A , vol. 372, no. 4, pp. 319-368, 2002.

[16] J. Ambjorn, P. Olesen, "On the formation of a random color magnetic quantum liquid in QCD," Nuclear Physics B , vol. 170, no. 1, pp. 60-78, 1980.

[17] M. A. Shifman, A. I. Vainshtein, V. I. Zakharov, "QCD and resonance physics. theoretical foundations," Nuclear Physics B , vol. 147, no. 5, pp. 385-447, 1979.

[18] M. A. Shifman, A. I. Vainshtein, V. I. Zakharov, "QCD and resonance physics. applications," Nuclear Physics B , vol. 147, no. 5, pp. 448-518, 1979.

[19] V. A. Novikov, M. A. Shifman, A. I. Vainshtein, V. I. Zakharov, " ^{η[variant prime]} Meson as a pseudoscalar gluonium," Physics Letters B , vol. 86, no. 3-4, pp. 347-350, 1979.

[20] A. Manohar, H. Georgi, "Chiral quarks and the non-relativistic quark model," Nuclear Physics B , vol. 234, no. 1, pp. 189-212, 1984.

[21] D. I. Dyakonov, V. Yu. Petrov, "A theory of light quarks in the instanton vacuum," Nuclear Physics B , vol. 272, no. 2, pp. 457-489, 1986.

[22] A. Yamamoto, H. Suganuma, "Lattice analysis for the energy scale of QCD phenomena," Physical Review Letters , vol. 101, no. 24, 2008.

[23] N. O. Agasian, Yu. A. Simonov, "Instantons in the stochastic QCD vacuum," Modern Physics Letters A , vol. 10, no. 24, pp. 1755-1760, 1995.

[24] D. Antonov, "A possible universal treatment of the field strength correlator in the Abelian-projected SU (2)-theory," Physics Letters B , vol. 506, no. 3-4, pp. 303-310, 2001., [email protected]

[25] D. Antonov, "Field correlators in abelian-projected theories and stochastic vacuum model," Journal of High Energy Physics , vol. 4, no. 7, pp. 18-32, 2000., [email protected]

[26] A. M. Badalian, Yu. A. Simonov, V. I. Shevchenko, "Method of vacuum correlation functions: results and prospects," Physics of Atomic Nuclei , vol. 69, no. 10, pp. 1781-1802, 2006., [email protected]

[27] Yu. A. Simonov, "A new nonperturbative approach to QCD and its application to hadronic physics," Yadernaya Fizika , vol. 54, no. 1, pp. 192-223, 1991.

[28] Yu. A. Simonov, L. Ferreira, P. Nogueira, J. I. Silva-Marcos, "QCD and theory of hadrons," QCD: Perturbative or Nonperturbative? , World Scientific, Singapore, 2001.

[29] M. Foster, C. Michael, "Hadrons with a heavy color-adjoint particle," Physical Review D , vol. 59, no. 9, 1999.

[30] Yu. A. Simonov, "Gluelump spectrum in the QCD string model," Nuclear Physics B , vol. 592, no. 1-2, pp. 350-368, 2000., [email protected]

[31] V. I. Shevchenko, Yu. A. Simonov, "Casimir scaling as a test of QCD vacuum models," Physical Review Letters , vol. 85, no. 9, pp. 1811-1814, 2000.

[32] V. I. Shevchenko, Yu. A. Simonov, "Casimir scaling and models of confinement in QCD," International Journal of Modern Physics A , vol. 18, no. 1, pp. 127-157, 2003., [email protected]

[33] G. S. Bali, "Casimir scaling of SU (3) static potentials," Physical Review D , vol. 62, no. 11, 2000., [email protected]

[34] A. Di Giacomo, H. Panagopoulos, "Field strength correlations in the QCD vacuum," Physics Letters B , vol. 285, no. 1-2, pp. 133-136, 1992.

[35] M. D'Elia, A. Di Giacomo, E. Meggiolaro, "Field strength correlators in full QCD," Physics Letters B , vol. 408, no. 1-4, pp. 315-319, 1997., [email protected]

[36] A. Di Giacomo, E. Meggiolaro, H. Panagopoulos, "Gauge-invariant field strength correlations in QCD at zero and non-zero temperature," Nuclear Physics B , vol. 483, no. 1-2, pp. 371-382, 1997., [email protected]

[37] M. D'Elia, A. Di Giacomo, E. Meggiolaro, "Gauge-invariant field-strength correlators in pure Yang-Mills theory and full QCD at finite temperature," Physical Review D , vol. 67, no. 11, 2003.

[38] G. S. Bali, N. Brambilla, A. Vairo, "A lattice determination of QCD field strength correlators," Physics Letters B , vol. 421, no. 1-4, pp. 265-272, 1998.

[39] D. V. Antonov, Yu. A. Simonov, "Exact equations for vacuum correlators in field theory," Physics of Atomic Nuclei , vol. 60, no. 3, pp. 478-489, 1997.

[40] Yu. A. Simonov, "Equations for field correlation functions," Physics of Atomic Nuclei , vol. 61, no. 5, pp. 855-864, 1998.

[41] Yu. A. Simonov, "Analytic calculation of field-strength correlators," Physics of Atomic Nuclei , vol. 69, no. 3, pp. 528-540, 2006.

[42] D. Antonov, H.-J. Pirner, "Following gluonic world lines to find the QCD coupling in the infrared," European Physical Journal C , vol. 51, no. 3, pp. 633-644, 2007., [email protected]

[43] D. Gromes, "Spin-dependent potentials in QCD and the correct long-range spin orbit term," Zeitschrift für Physik C , vol. 26, no. 3, pp. 401-406, 1984.

[44] N. Brambilla, D. Gromes, A. Vairo, "Poincaré invariance and the heavy-quark potential," Physical Review D , vol. 64, no. 7, 2001., [email protected]; [email protected]; [email protected]

[45] Y. Koma, M. Koma, "Spin-dependent potentials from lattice QCD," Nuclear Physics B , vol. 769, no. 1-2, pp. 79-107, 2007.

[46] A. M. Badalian, A. V. Nefediev, Yu. A. Simonov, "Dynamical suppression of the spin-orbit interaction in hadrons," JETP Letters , vol. 88, no. 3, pp. 151-156, 2008.

[47] Yu. A. Simonov, J. A. Tjon, "The Feynman-Schwinger representation in QCD," Annals of Physics , vol. 300, no. 1, pp. 54-87, 2002.

[48] Yu. A. Simonov, J. A. Tjon, "The Feynman-Schwinger representation for the relativistic two-particle amplitude in field theory," Annals of Physics , vol. 228, no. 1, pp. 1-18, 1993.

[49] Yu. A. Simonov, "Heavy-light mesons and quark constants _fB , _fD in the vacuum correlator method," Zeitschrift für Physik C , vol. 53, no. 3, pp. 419-431, 1992.

[50] S. V. Ivanov, G. P. Korchemskij, "Some supplements of the nonperturbative gauges," Physics Letters B , vol. 154, no. 2-3, pp. 197-201, 1985.

[51] V. I. Shevchenko, Yu. A. Simonov, "Generalized coordinate gauge, nonabelian Stokes theorem and the dual QCD Lagrangian," Physics Letters B , vol. 437, no. 1-2, pp. 146-152, 1998., [email protected]; [email protected]

[52] L. Lukaszuk, E. Leader, A. Johansen, "Contour gauges, canonical formalism and flux algebras," Nuclear Physics B , vol. 562, no. 1-2, pp. 291-314, 1999.

[53] B. L. Ioffe, "QCD (Quantum chromodynamics) at low energies," Progress in Particle and Nuclear Physics , vol. 56, no. 1, pp. 232-277, 2006.

[54] A. O. Barvinsky, G. A. Vilkovisky, "The generalized Schwinger-DeWitt technique in gauge theories and quantum gravity," Physics Reports A , vol. 119, no. 1, pp. 1-74, 1985.

[55] M. B. Voloshin, "On dynamics of heavy quarks in a non-perturbative QCD vacuum," Nuclear Physics B , vol. 154, no. 3, pp. 365-380, 1979.

[56] M. B. Voloshin, "Precoulombic asymptotics for energy levels of heavy quarkonium," Sov. J. Nucl. Phys. , vol. 36, pp. 162-143, 1982.

[57] D. Gromes, "Space-time dependence of the gluon condensate correlation function and quarkonium spectra," Physics Letters B , vol. 115, no. 6, pp. 482-486, 1982.

[58] U. Marquard, H. G. Dosch, "Potential and sum-rule approach in QCD," Physical Review D , vol. 35, no. 7, pp. 2238-2243, 1987.

[59] M. Campostrini, A. Di Giacomo, S. Olejník, "On the possibility of detecting gluon condensation from the spectra of heavy quarkonia," Zeitschrift für Physik C , vol. 31, no. 4, pp. 577-582, 1986.

[60] Yu. A. Simonov, S. Titard, F. J. Ynduráin, "Heavy quarkonium systems and nonperturbative field correlators," Physics Letters B , vol. 354, no. 3-4, pp. 435-441, 1995., [email protected]; [email protected]; [email protected]

[61] B. S. DeWitt, "Quantum theory of gravity. II. The manifestly covariant theory," Physical Review , vol. 162, no. 5, pp. 1195-1239, 1967.

[62] J. Honerkamp, "The question of invariant renormalizability of the massless Yang-Mills theory in a manifest covariant approach," Nuclear Physics B , vol. 48, no. 1, pp. 269-287, 1972.

[63] G. 't Hooft, "An algorithm for the poles at dimension four in the dimensional regularization procedure," Nuclear Physics B , vol. 62, pp. 444-460, 1973.

[64] G. 't Hooft, "Lectures at Karpacz 1975," in Proceedings, Acta Universitatis Wratislaviensis No. 368, Vol. 1, pp. 345-369, Wroclaw, Poland, 1976.

[65] L. F. Abbott, "The background field method beyond one loop," Nuclear Physics B , vol. 185, no. 1, pp. 189-203, 1981.

[66] M. Eidemüller, M. Jamin, "QCD field strength correlator at the next-to-leading order," Physics Letters B , vol. 416, no. 3-4, pp. 415-420, 1998.

[67] V. I. Shevchenko, "QCD field strength correlator at the one-loop order," http://arxiv.org/abs/hep-ph/9802274

[68] A. M. Polyakov, "String representations and hidden symmetries for gauge fields," Physics Letters B , vol. 82, no. 2, pp. 247-250, 1979.

[69] V. S. Dotsenko, S. N. Vergeles, "Renormalizability of phase factors in nonabelian gauge theory," Nuclear Physics B , vol. 169, no. 5-6, pp. 527-546, 1980.

[70] V. I. Shevchenko, Yu. A. Simonov, "Generalized coordinate gauge, nonabelian Stokes theorem and the dual QCD Lagrangian," Physics Letters B , vol. 437, no. 1-2, pp. 146-152, 1998.

[71] N. O. Agasian, Yu. A. Simonov, "New nonperturbative approach to the Debye mass in hot QCD," Physics Letters B , vol. 639, no. 2, pp. 82-87, 2006.

[72] V. I. Zakharov, "From confining fields back to power corrections," Nuclear Physics B--Proceedings Supplements , vol. 164, pp. 240-247, 2007.

[73] Yu. A. Simonov, "Nonperturbative dynamics of heavy quarkonia," Nuclear Physics B , vol. 324, no. 1, pp. 67-90, 1989.

[74] M. Schiestl, H. G. Dosch, "Gluon condensate and spin dependent potentials," Physics Letters B , vol. 209, no. 1, pp. 85-89, 1988.

[75] G. S. Bali, "Are there short distance non-perturbative contributions to the QCD static potential?," Physics Letters B , vol. 460, no. 1-2, pp. 170-175, 1999.

[76] A. M. Badalian, D. S. Kuzmenko, "Freezing of QCD coupling affects the short distance static potential," Physical Review D , vol. 65, no. 1, 2002., [email protected]; [email protected]

[77] A. M. Badalian, Yu. A. Simonov, "Spin-dependent potentials and field correlation functions in QCD," Physics of Atomic Nuclei , vol. 59, no. 12, pp. 2164-2180, 1996.

[78] A. M. Badalian, Yu. A. Simonov, "Spin dependent potentials and field correlation functions in QCD," Yadernaya Fizika , vol. 59, pp. 2247, 1996.

[79] Yu. A. Simonov, "Spin-dependent interaction in the deconfined phase of QCD," Sense of Beauty in Physics, A. Di Giacomo's Festschrift , pp. 29, Pisa University Press, 2006.

[80] E. Eichten, F. Feinberg, "Spin-dependent forces in quantum chromodynamics," Physical Review D , vol. 23, no. 11, pp. 2724-2744, 1981.

[81] Yu. A. Simonov, "Perturbation theory in the nonperturbative QCD vacuum," Physics of Atomic Nuclei , vol. 58, pp. 107-123, 1995.

[82] Yu. A. Simonov JETP Letters , vol. 75, pp. 525, 1993.

[83] Yu. A. Simonov, H. Latal, W. Schweiger, "Perturbation theory in a nonperturbative QCD background," Perturbative and Nonperturbative Aspects of Quantum Field Theory (Schladming, 1996) , vol. 479, of Lecture Notes in Physics, pp. 139-207, Springer, Berlin, Germany, 1997.

[84] A. M. Badalian, "Is the _a0 (1450) a candidate for the lowest qq¯ ³_P0 state?," Physics of Atomic Nuclei , vol. 66, no. 7, pp. 1342-1348, 2003.

[85] A. M. Badalian, B. L. G. Bakker, " _ηc (3654) and hyperfine splitting in charmonium," Physical Review D , vol. 67, no. 7, 2003.

[86] A. M. Badalian, V. L. Morgunov, B. L. G. Bakker, "Fine-structure splittings of excited P- and D-wave states in charmonium," Physics of Atomic Nuclei , vol. 63, no. 9, pp. 1635-1639, 2000.

[87] A. M. Badalian, V. L. Morgunov, "Determination of _αs (1 GeV) from the charmonium fine structure," Physical Review D , vol. 60, no. 11, 1999., [email protected]; [email protected]

Appendix

Appendix

Nonperturbative Contributions to _D1 (z) and D(z) at Small z

We start with _D1 (z) expressed via gluelump Green's function in (4.8), (4.10) _Gμν (z)=_δμν f(z) .

The leading NP behavior at small z , proportional to ...F(x)F(y)... , is obtained from the amplitudes, shown in Figures 1 and 2. For the amplitude of Figure 1, one can use the vacuum average of the W-loop (2.11) in the limit of small contours and neglect the factor exp F in (4.11). One obtains the adjoint loop [figure omitted; refer to PDF] where γ=_C2 (adj)/_C2 (f)=2^Nc2 /(^Nc2 -1) (original derivation in [41] referred to the fundamental loop). Proceeding as in Appendix 2 of [41], one obtains the NP correction Δ_D1 as [figure omitted; refer to PDF] where we have introduced the standard definition [17, 18] _G2 ≡(_αs /π)...^Fμνa (0)^Fμνa (0)... . One can check that Δ_D1 =...AA;(^Nc0 ) .

We turn now to the amplitude of Figure 2. The corresponding f(x) for x[arrow right]0 can be written as follows: [figure omitted; refer to PDF] where the gluon Green's function G(x,y) for small |x-y| can be replaced by the perturbative part: G(x,y)[arrow right]_G0 (x,y)=1/4^π2 (x-y⁾² .

Taking into account that Δ_D1 (x)=-(2^g2 /^Nc2 )(d/d^x2 )Δf(x) , one obtains for x[arrow right]0 the following: [figure omitted; refer to PDF] where we have defined [figure omitted; refer to PDF]

The integral I(0) diverges at small v and large u,v . The latter divergence is removed since at large arguments G(x,y) is damped by confining force, keeping the propagating gluon nearby the 4th axis _x4 =_y4 =0 where the static gluonic source resides. One can estimate I(x[arrow right]0)~...AA;((ln (^λ2 /^x2 )^)α ) . Since _G2 is connected to D , _D1 , then [figure omitted; refer to PDF] and one can see in (A.4) that coefficients of ^D1E on both sides of (A.4) have the same order of magnitude and the same sign, suggesting a self-consistency on this preliminary level. However, as we argued before and in [41], this contribution is actually gluon mass renormalization, which is zero. This is especially clear when ...FF... in Figure 2 is replaced by its perturbative part.

We turn now to the most important case of the confining FC ^DE (x) . The corresponding two-gluon gluelump amplitude is depicted in Figure 3, and the np contributions to ^DE (x) are given in Figures 4 and 5.

Figure 3: [figure omitted; refer to PDF]

Figure 4: [figure omitted; refer to PDF]

Figure 5: [figure omitted; refer to PDF]

Figure 6: [figure omitted; refer to PDF]

We start with the amplitude of Figure 4, which can be represented as a doubled diagram of Figure 1, and therefore, the contribution to the two-gluon gluelump Green's function ^G(2gl) (x) will be [figure omitted; refer to PDF] Now using relation (2.11), one obtains for the np contribution of Figure 4 [figure omitted; refer to PDF] Note that this contribution is finite at x[arrow right]0 , however, with the negative sign. We now turn to the np contribution of Figure 5, which can be obtained from (4.12) inserting there the product of operators _F...fg (w)^d4 w_{F...^{f[variant prime]g[variant prime]}} (^{w[variant prime]} )^{d4w[variant prime]} , where fg(^{f[variant prime]g[variant prime]} ) are adjoint color indices, _F...fg =^FaTfga =^Fa (-i)_fafg . Using the formulas [figure omitted; refer to PDF] one arrives at the expression [figure omitted; refer to PDF] where G(y) is the gluon Green's function in the two-gluon gluelump; the gluon is connected at large distances by two strings to another gluon and to the static gluon source; see Figure 5. At small x (we take it along axis 4 for convenience), the integral in h(x) grows when G(y) is close to the 4th axis and becomes the free gluon G(y)[arrow right]_G0 (y)=1/4^π2y2 . As a result, one obtains [figure omitted; refer to PDF] and one should have in mind that the integral for x[arrow right]0 diverges both at small and at large w , ^{w[variant prime]} . The small w , ^{w[variant prime]} region will give the terms O(ln (^λ02 /^x2 )) , and at large w , ^{w[variant prime]} , the integral is protected by the falloff of G(y) at large y due to confinement; therefore, we must imply in the integrals in (A.11) the upper limits for w , ^{w[variant prime]} at some _λ0 [> ~]1/_σadj . Then introducing for w(w,_w4 ) polar coordinates |w|=ρsinθ , _w4 =ρcos θ (and the same for ^{w[variant prime]} ), one arrives at the integral [figure omitted; refer to PDF] or (π/2^x2 ) (x≥ρ) , and finally one has the estimate of contribution to (A.10) from the region of small w,^{w[variant prime]} (w|w|,|^{w[variant prime]} |[<, ~]_λ0 ) as [figure omitted; refer to PDF] And finally for D(x), one obtains (omitting the perturbative term in (A.7)) [figure omitted; refer to PDF]

Note that the last term on the r.h.s. dominates at small x and compensates the decrease of the ^αs2 term since [figure omitted; refer to PDF] Hence,

[figure omitted; refer to PDF]