Chiba Univ./KEK Preprint CHIBAEP165
KEK Preprint 200719
July 2007
Compact lattice formulation of
ChoFaddeevNiemi decomposition:
gluon mass generation and infrared Abelian dominance
A. Shibata, S. Kato, K.I. Kondo, T. Murakami, T. Shinohara S. Ito
Computing Research Center, High Energy Accelerator Research Organization (KEK),
& Graduate Univ. for Advanced Studies (Sokendai), Tsukuba 3050801, Japan
Takamatsu National College of Technology, Takamatsu 7618058, Japan
Department of Physics, Faculty of Science, Chiba University, Chiba 2638522, Japan
Graduate School of Science and Technology, Chiba University, Chiba 2638522, Japan
Nagano National College of Technology, 716 Tokuma, Nagano 3818550, Japan
This paper complements a new lattice formulation of SU(2) YangMills theory written in terms of new variables in a compact form proposed in the previous paper. The new variables used in the formulation were once called the Cho–Faddeev–Niemi or Cho–Faddeev–Niemi–Shabanov decomposition. Our formulation enables us to explain the infrared “Abelian” dominance, in addition to magnetic monopole dominance shown in the previous paper, in the gauge invariant way without relying on the specific gauge fixing called the maximal Abelian gauge used in the conventional investigations. In this paper, especially, we demonstrate by numerical simulations that gluon degrees of freedom other than the “Abelian” part acquire the mass to be decoupled in the lowenergy region leading to the infrared Abelian dominance.
Key words: lattice gauge theory, Abelian dominance, magnetic monopole, quark confinement
PACS: 12.38.Aw, 12.38.Lg Email: akihiro.shibata@kek.jp
Email: kato@takamatsunct.ac.jp
Email: kondok@faculty.chibau.jp
Email: tom@cuphd.nd.chibau.ac.jp
Email: sinohara@graduate.chibau.jp
Email: shoichi@ei.naganonct.ac.jp
1 Introduction
The dual superconductivity in YangMills vacuum [1] is believed to be a promising mechanism [2] for quark confinement. For this picture to be valid, it is supposed that one can extract the relevant Abelian component responsible for duality from nonAbelian gauge theory, i.e., YangMills theory, since the duality is understood as the electric–magnetic dual in the Abelian gauge theory represented by the Maxwell theory. Therefore, it is important to give a prescription to extract a variable playing the role of such an Abelian part from the original nonAbelian gauge theory in some sense. However, such a variable is not necessarily the Abelian gauge field in its appearance. To emphasize this situation, we use double quotation marks as “Abelian”.
In the previous paper [3], we have proposed a formulation of YangMills theory on a lattice in a compact form (say, compact lattice formulation). The compact lattice formulation was constructed so that it reduces in the (naive) continuum limit to the formulation of YangMills theory written in terms of new variables obtained through nonlinear change of variables (NLCV) from the original gauge field [4]. The NLCV generates the new variables which have the same form as the Cho–Faddeev–Niemi (CFN) or Cho–Faddeev–Niemi–Shabanov (CFNS) decomposition [5, 6, 7].
Prior to the compact lattice formulation, we have already given another lattice formulation in a noncompact form (referred to as the noncompact lattice formulation hereafter) [8]. The noncompact lattice formulation enabled one to define the magnetic monopole in the gauge invariant way in YangMills theory on a lattice without introducing fundamental scalar fields. This is a remarkable result, since the conventional approach of defining the magnetic monopole in YangMills theory without fundamental scalar fields heavily relies on a specific choice of gauge fixing, the socalled the maximal Abelian gauge (MAG) [9], which breaks the color symmetry explicitly in addition to the local gauge symmetry. The gaugeinvariant magnetic monopole in pure YangMills theory has been constructed by introducing the unit vector field based on NLCV which plays the role of recovering color symmetry broken by a specific choice of color direction according to a Cartan decomposition in MAG. However, the magnetic charge resulting from the magnetic monopole defined in this way is not integervalued in the noncompact formulation.
This drawback was remedied by the subsequent compact formulation [3] which guarantees that the magnetic charge is integervalued and obeys the Dirac quantisation condition. Moreover, the infrared “Abelian” dominance and magnetic monopole dominance in the string tension were demonstrated by numerical simulations in the compact formulation, although such phenomena were found for the first time in the MAG [10, 11, 12, 13]. These results strongly support the dual superconductor picture of QCD vacuum as a promising mechanism of quark confinement. However, what is the mechanism for the infrared “Abelian” dominance or magnetic monopole dominance is an unanswered question in this investigation. From a theoretical point of view [14], on the other hand, it was clarified which variables should be identified with the “Abelian”part which is responsible for quark confinement within the continuum formulation [4] so that the variable gives the dominant contribution to the string tension to be calculated from the Wilson loop average (infrared “Abelian” dominance). This means that the remaining variable decouples in the lowenergy or longdistance region to become irrelevant for the string tension, once such an identification of the “Abelian” part is achieved.
Our continuum formulation of YangMills theory allows one to introduce the mass term for the remaining field without breaking the local gauge invariance. Therefore, the dynamical mass generation for is not prohibited in this formulation. If such gluon mass for is generated, it could be a mechanism for infrared “Abelian” dominance as pointed out in [14]. The dynamical generation of the gluon mass for yields the decoupling of these degrees of freedom in the lowenergy region leaving the “Abelian” part as the lowenergy modes relevant to quark confinement. To confirm this scenario for dynamical Abelian projection is a main motivation of this paper. For MAG, mass generation was so far reported for the offdiagonal gluon component [15]. In this paper, we have measured the mass for the gluon directly by numerical simulations based on our lattice formulation [3]. In this paper, moreover, we discuss in detail how to define the lattice variable corresponding especially to the remaining part in the compact lattice formulation to obtain the continuum counterparts to the errors of lattice spacing . Preliminary results have already been reported in [16]. Finally, it is worth mentioning that the gauge invariance of the Abelian confinement mechanism was also discussed recently in many publications by other group, e.g., [17].
2 Compact formulation
In order to consider a compact lattice version for the new formulation of YangMills theory, we recall the continuum formulation presented in [4]. In the continuum formulation [5, 4], we have introduced a color vector field of a unit length, i.e., . In what follows, we use the arrow to denote the vector and use the boldface letter to express the Liealgebra valued field, e.g., with generators of , where () are Pauli matrices. The valued gluon field (gauge potential) is decomposed into two parts:
(1) 
in such a way that the color vector field is covariantly constant in the background field :
(2) 
and that the remaining field is perpendicular to :
(3) 
Here we have introduced the gauge coupling and we have adopted the normalization for generators: . Note that , , and are realvalued fields and their Liealgebra forms are Hermitian due to Hermiticity of the generators .
By solving the defining equation (2), the field is obtained in the form:
(4) 
where the second term is perpendicular to , i.e., . Here it should be remarked that the parallel part , proportional to can not be determined uniquely from the defining equation (2). Imposing the perpendicular condition (3) determines and the remaining part as
(5) 
It is easy to check that the sum of and specified respectively by (4) and (5) agrees with the original field according to (1).
On a lattice, on the other hand, we introduce the site variable constructed according to [8], in addition to the original link variable . Note that we define a color vector field on the lattice corresponding to the continuum notation . In this paper, we define the link variable as the exponential of the line integral of a gauge potential along a link from to :
(6) 
where denotes the lattice spacing and denotes the path ordering. In the explicit estimation of the naive continuum limit, we adopt in this paper the midpoint definition for the link variable:
(7) 
using the midpoint of the link running from to . This prescription is adopted to suppress as much as possible lattice artifacts coming from a finite (nonzero) lattice spacing, in contrast to our previous paper [3] where we have adopted the very naive definition:
The link variable and the site variable transform under the gauge transformation II [4] as
(8) 
Note that is Hermitian, , and is unitary, . It should be remarked that this transformation property follows from the most general form (6) for the link variable , irrespective of the prescription for the discrete lattice approximation.
The lattice variables and corresponding to and should be expressed in terms of the site variable and the original link variable , just as the continuum variables and are expressed in terms of and , However, the definition of lattice variables and is not unique. They must be defined in a consistent way with the defining equation on a lattice respecting the transformation property. We achieve this by solving a lattice version [3] of (2) and (3):
(9) 
(10) 
The defining equation must be invariant under the gauge transformation, namely, they are forminvariant: and . We identify the lattice variable with a link variable which transforms in the same way as the original link variable :
(11) 
This requirement guarantees that the defining equation (9) is gauge invariant. On the other hand, we define the lattice variable so that it transforms in just the same way as the site variable :
(12) 
to realize the adjoint color rotation at the site suggested from the transformation property of the continuum variable. By this choice, indeed, the orthogonality condition (10) is kept gauge invariant.
Explicit construction of the new lattice variables are as follows. We define as a link variable which is a group element of related to the valued background field through
(13) 
In the midpoint definition the link variable reads
(14) 
where is to be identified with the continuum variable defined by (4) in the continuum limit. Hence must be unitary . The same remark as the the link variable for the naive continuum limit holds also for the link variable . In the previous paper [3], the lattice version (9) of the defining equation (2) has been solved and the resulting link variable is of the form (up to the normalization) [18]:
(15) 
and the unitary link variable has been obtained after the normalization:
(16) 
Indeed, the naive continuum limit of the link variable (16) reduces to the continuum expression (4).
A naive choice for the lattice variable is given by or These are suggested from the relation . In fact, they satisfy the desired transformation property (12). Note that or is excluded, since it obeys the adjoint rotation at or , not at . Then we can construct a lattice variable as the linear combination:
(17) 
to satisfy the desired transformation property (12). Now we can see that it is reasonable to adopt (10) as a lattice version of the orthogonality equation (3). In fact, implies (See Appendix for the derivation of (18))
(18) 
In this way, the lattice variables is expressed in terms of the site variable and the original link variable as
(19) 
In particular, a good choice is obtained for the symmetric case, i.e., , since this choice enables us to define the lattice variable so as to reproduce the naive continuum limit of the orthogonality equation up to .
Finally, we obtain the unitary lattice variable after the normalization:
(20) 
In the numerical simulations, we have adopted the choice by the reason mentioned above. Moreover, there are some arbitrariness for extracting the Liealgebra valued variable from the compact lattice variable . This issue will be examined by comparing the results of numerical simulations.
3 Abelian dominance and gluon mass generation
3.1 Identifying the Abelian part and mass term for the remaining part
The field can be regarded as the “Abelian” part in the reformulated YangMills theory by the following reasons.

The “Abelian” part corresponds to the diagonal part of the gauge potential in the context of the conventional MAG which is reproduced when the color vector is aligned in the same direction over the whole spacetime, for example,
(21) 
The Wilson loop average in YangMills theory written in terms of is rewritten into the reduced Wilson loop average which is entirely rewritten in terms of in the reformulated YangMills theory, as demonstrated in [14].

The mass term for can be introduced without breaking gauge invariance in this reformulation [4]. In fact, it has been shown to oneloop order [14] that such an effective mass term is generated due to the gaugeinvariant dimension two condensate thanks to the gauge invariant selfinteraction term among gluons, in sharp contrast to the ordinary selfinteraction term which is not gaugeinvariant.
Therefore, in the energy region lower than the mass of the field , the remaining components should decouple or negligible and the field could be dominant. This leads to the infrared Abelian dominance (in the string tension) in our reformulation.
Keeping these facts in mind, we proceed to obtain a fitting function of the two–point correlation function. Suppose that the YangMills theory has the effective mass term:
(22) 
An additional quadratic term in of the following type could be generated from gauge fixing conditions in the differential form [4].
(23) 
This can be understood as follows. Recall that we impose an constraint called the new Maximal Abelian gauge (nMAG) to obtain the reformulated YangMills theory with the original gauge symmetry SU(2) even after introducing the color field which apparently increases gauge degrees of freedom [4]. Then we introduce a gaugefixing parameter for nMAG of the form: . This term does not fix the SU(2) gauge invariance. Therefore, we adopt the Landau gauge for the overall gauge fixing of whose differential form is . This gives an additional quadratic term: coming from the GF term: . Therefore, combining two terms yields an additional term quadratic in : with . Thus we assume the effective propagator for gluon of the form:
(24) 
In particular, the limit reproduces the Proca case:
(25) 
This form was adopted in the study of offdiagonal gluon mass generation in MAG [15] where the mass term was introduced by hand without preserving the gauge invariance. Note that both nMAG and Landau gauge conditions are exactly satisfied only at and . This is realized at limit:
(26) 
Therefore, the limit differs from the previous Proca case used in MAG. However, it will turn out below that the constant shift of the propagator gives the same decay rate and hence the same mass of gluon.
3.2 Numerical simulations
We have generated configurations of link variables based on the standard heat bath method for the standard Wilson action. The numerical simulation are performed at , on lattice, at , , on lattice, at , , on lattice, and at on lattice by thermalizing 15000 sweeps. Here 200 configurations are stored every 300 sweeps. Other settings of numerical simulations are the same as those in the previous paper [3].
We are now ready to study characteristic features of the reformulated YangMills theory written in terms of new variables defined through NLCV of the original field variable : infrared Abelian dominance, magnetic monopole dominance and nonvanishing gluon mass. Among them, the magnetic monopole dominance in the string tension has already been confirmed in the previous paper [3] using the gaugeinvariant magnetic monopole which is guaranteed to have integervalued magnetic charge subject to the Dirac quantization condition according to our construction of magnetic current based on NLCV. ^{1}^{1}1 The proposed NLCV enables one to extract the “Abelian part” irrespective of the choice of the gauge fixing preserving the color symmetry. The YangMills theory in the conventional MAG is reproduced as a very special limit (21) of our reformulated YangMills theory based on NLCV. An advantage of our formulation is that we can confirm such characteristic features for any choice of gauge fixing, not restricted to MAG, since our formulation allows us to take arbitrary type of gauge fixing for the original variable .
To study the infrared Abelian dominance and the gluon mass generation in the reformulated YangMills theory, we first define the twopoint correlation functions (full propagators) for the independent variables in the new formulation on a lattice, i.e., , and , in addition to the original variable . For simplicity, we examine just the contracted scalartype propagator simplified by avoiding the complicated tensor structure:
(27) 
and
(28) 
Here the Liealgebra valued gauge potential or is defined from the respective link variable by
(29) 
For the variable , on the other hand, we examined two options: one is extracted from decomposing the gauge potential (groupvalued):
(30) 
and the other is from the definition of the decomposition (Liealgebravalued):
(31) 
The field is defined by
(32) 
The numerical results are presented in Fig. 1. As is quickly observed from the left panel of Fig. 1, and exhibit quite similar behaviors in the measured range of the Euclidean distance . In order to determine the physical scale, we have used the relationship between the (inverse) gauge coupling and lattice spacing given in [19] ^{2}^{2}2 We use the relationship between the physical units, or . This comes from . which is summarized in Table 1.
lattice spacing  lattice size  

From the right panel of Fig. 1, ( or ) is dominant compared to which decreases more rapidly than other correlation functions in . This implies the infrared “Abelian” dominance, provided that the components composed of and are identified with the “Abelian” part of . As is seen from the left panel of Fig. 1, a nontrivial mixed correlation function exists, since includes a perpendicular component to .
Fig. 1 demonstrates nice independence of our results against variations of the ultraviolet cutoff (the lattice spacing ). The propagators calculated at the lattices with different follow the same curve if plotted in the physical units. These accurate plots provide an additional support that the results presented in this paper are definitely not lattice artifacts.
Note that we must impose the gauge fixing condition for the original variable to obtain the correlation function. In our simulations, we have chosen the lattice Landau gauge (LLG) for the original field for this purpose. Thus we have confirmed the infrared “Abelian” dominance with color symmetry being kept, since the Landau gauge keeps the color symmetry. This is one of our main results. The infrared Abelian dominance was so far obtained only for the MAG which breaks the color symmetry explicitly. As already mentioned, moreover, we can choose any other gauge and we can study using this formulation if the infrared “Abelian” dominance can be observed in any other gauge. We hope we can report the results in the other gauge in future investigations.
Next, we determine the gluon mass generated in the nonperturbative way by examining the correlation functions in more detail. The gauge boson propagator is related to the Fourier transform of the massive propagator :
(33) 
Then the scalartype propagator as a function of should behave for large as (See [15] for details of the integral calculation.)
(34) 
Therefore, the scaled propagator should be proportional to for with being the damping rate of . In other words, the mass of the gauge field can be estimated from the slope in the logarithmic plot of the scaled propagator as a function of . ^{3}^{3}3 Here we have assumed that the anomalous dimension is sufficiently small so that the exponent of the power of is the same as the tree value.
Fig. 2 shows the logarithmic plots of the scaled scalartype propagator for , and as a function of the distance measured in the physical unit [fm] and in unit of square root of the string tension MeV. According to Fig. 2, we find just small difference between two types of defined by (30) or (31) over several choices of lattice spacing (i.e., several values of , ). Therefore, we can use either definition of the lattice variable to obtain in the consistent manner.
In Fig. 3, the measured values for the gluon mass are plotted as the function of the inverse lattice volume in the physical unit, to study the finitesize effect on the mass. The finite latticesize effect seems to be small for the gluon mass . Here the error bars originate from the fitting procedure for obtaining the slope, but no systematic errors such as finitevolume are included. In this way, we have estimated the mass for the gluon:
(35) 
Even after the whole gauge fixing, our formulation preserves color symmetry in sharp contrast to the conventional MA gauge. In view of the fact that our reformulated YangMills theory reproduces the YangMills theory in MA gauge as a special limit, the remaining part could correspond to the offdiagonal part in this limit. From this point of view, our result is consistent with the result obtained for the offdiagonal gluon mass in MAG [15].
Moreover, we have simultaneously estimated the decay rate for the new fields , , and the original gauge field by imposing the LLG as the overall gauge fixing. For , the decay rate is extracted according to the fitting: . Fig. 3 indicates not so small finite volume effect for data of . Using the data of , therefore, we have estimated the decay rate (or “mass”) as
(36) 
The decay rate obtained from the correlation function of field is slightly larger than that expected from the result in MAG. It should be remarked that the decay rate for the correlation function of field is extracted according to the fitting function which is not yet justified from the theoretical consideration. This might be an origin of the large value of . More simulations on the larger lattice are expected to eliminate finite volume effect for these values. However, we have no argument for guaranteeing the gauge invariance of these values or for identifying these values with their “masses”. In fact, the field is not gauge invariant. These issues will be checked in further investigations based on our reformulation.
Finally, we comment on the “Abelian”part , since our treatment of the “Abelian”part is different from the conventional approach based on MAG. The above result yields the “mass” of the “Abelian”part : . This value is nearly equal to that of the diagonal gluon mass obtained by imposing the Landau gauge in the conventional approach as reported in the second paper of [15] where the Landau gauge was imposed on the Abelian diagonal part in addition to the MAG for offdiagonal gluon field defined by the Cartan decomposition (). Therefore, the prescription of gauge fixing in [15] is different from ours.
4 Conclusion and discussion
In this paper we have developed a compact lattice formulation of SU(2) YangMills theory proposed in the previous paper [3] as the lattice version of the NLCV which was once called the CFN or CFNS decomposition. This resolves all drawbacks of the previous noncompact lattice formulation of our own [8].^{4}^{4}4 This is done up to specifying the integration measure for the respective new variable. It was not necessary to resolve this issue for obtaining the results reported in this paper. This compact formulation has enabled one to define the gaugeinvariant magnetic monopole with the magnetic charge subject to Dirac quantisation condition and to extract the “Abelian” part yielding the infrared “Abelian” dominance in the string tension for any choice of the gauge fixing for the original gauge field in the original YM theory.
In order to confirm the dynamical mass generation for the remaining part as a mechanism for the infrared “Abelian” dominance, we have measured the twopoint correlation function (the full propagator in real space) in our lattice formulation by imposing LLG for the original gauge field as the whole gauge fixing. We have found the infrared “Abelian” dominance in the sense that the propagator is suppressed in the long distance compared to and (and propagators as an immediate consequence of dynamically generated mass for (which is larger than the decay rate of other gluon field propagators).
Even after the whole gauge fixing, our formulation can preserve color symmetry by choosing the gaugefixing condition which does not break color symmetry, e.g., Landau gauge. This opens a path to examine color confinement in the same framework as quark confinement in the dual superconductivity picture. This feature is in sharp contrast to the conventional MA gauge breaking color symmetry, although our formulation reproduces the MA gauge as a special limit (21). It is important to demonstrate explicitly the gaugefixing independence of our results obtained in this paper for establishing the gaugeinvariant mechanism for quark confinement.
Acknowledgments
The numerical simulations have been done on a supercomputer (NEC SX5 and NEC SX8R) at Research Center for Nuclear Physics (RCNP), Osaka University. This project is also supported by the Large Scale Simulation Program No.0617 (FY2006) and No.0715 (FY2007) of High Energy Accelerator Research Organization (KEK). This work is financially supported by GrantinAid for Scientific Research (C)18540251 from JSPS and in part by GrantinAid for Scientific Research on Priority Areas (B)13135203 from MEXT.
Appendix A The accuracy of the naive continuum limit
For the naive continuum limit , we show that the lattice variable defined by yields for , while for . The repeated use of the BakerCampbellHausdorff formula yields
(37) 
where we have used , , , and . Thus we have
(38) 
The statement follows from the fact that the choice eliminates order terms.
References
 [1] C.N. Yang and R.L. Mills, Phys. Rev. 96, 191195 (1954).

[2]
Y. Nambu,
Phys. Rev. D 10, 4262(1974).
G. ’t Hooft, in: High Energy Physics, edited by A. Zichichi (Editorice Compositori, Bologna, 1975).
S. Mandelstam, Phys. Report 23, 245(1976).
A.M. Polyakov, (1975). Nucl. Phys. B 120, 429(1977).  [3] S. Ito, S. Kato, K.I. Kondo, T. Murakami, A. Shibata and T. Shinohara, [heplat/0604016], Phys. Lett. B 645, 67 (2007).

[4]
K.I. Kondo, T. Murakami and T. Shinohara,
[hepth/0504107],
Prog. Theor. Phys. 115, 201 (2006).
K.I. Kondo, T. Murakami and T. Shinohara, [hepth/0504198], Eur. Phys. J. C 42, 475 (2005).  [5] Y.M. Cho, Phys. Rev. D 21, 1080 (1980). Phys. Rev. D 23, 2415 (1981).
 [6] L. Faddeev and A.J. Niemi, [hepth/9807069], Phys. Rev. Lett. 82, 1624 (1999).

[7]
S.V. Shabanov,
[hepth/9903223],
Phys. Lett. B 458, 322 (1999).
S.V. Shabanov, [hepth/9907182], Phys. Lett. B 463, 263 (1999).  [8] S. Kato, K.I. Kondo, T. Murakami, A. Shibata, T. Shinohara and S. Ito, [heplat/0509069], Phys. Lett. B 632, 326 (2006).
 [9] A. Kronfeld, M. Laursen, G. Schierholz and U.J. Wiese, Phys.Lett. B 198, 516(1987).
 [10] G. ’t Hooft, Nucl.Phys. B 190 [FS3], 455(1981).
 [11] Z.F. Ezawa and A. Iwazaki, Phys. Rev. D 25, 2681(1982).
 [12] T. Suzuki and I. Yotsuyanagi, Phys. Rev. D 42, 4257 (1990).
 [13] J.D. Stack, S.D. Neiman and R. Wensley, [heplat/9404014], Phys. Rev. D50, 3399 (1994). H. Shiba and T. Suzuki, Phys.Lett.B333, 461 (1994).
 [14] K.I. Kondo, [hepth/0609166], Phys. Rev. D 74, 125003 (2006).

[15]
K. Amemiya and H. Suganuma,
[heplat/9811035],
Phys. Rev. D60, 114509 (1999).
V.G. Bornyakov, M.N. Chernodub, F.V. Gubarev, S.M. Morozov and M.I. Polikarpov, [heplat/0302002], Phys. Lett. B559, 214222 (2003).  [16] A. Shibata, S. Ito, S. Kato, K.I. Kondo, T. Murakami and T. Shinohara, Talk given at 24th International Symposium on Lattice Field Theory (Lattice 2006), heplat/0610023, PoS(LAT2006)074.
 [17] T.Suzuki et al, [heplat/0410001], Phys. Rev. Lett. 94, 132001 (2005). heplat/0511025, heplat/0511026, heplat/0511027.
 [18] M.N. Chernodub, F.V. Gubarev and E.M. Ilgenfritz, [heplat/9710011], Phys. Lett. B 424, 106—114 (1998).
 [19] S. Kato, S. Kitahara, N. Nakamura and T. Suzuki, Nucl. Phys. B520, 323344 (1998).