Subleading properties of the QCD fluxtube in 3d lattice gauge theory
Abstract
We study the continuum limit of the stringlike behaviour of flux tubes formed between static quarks and antiquarks in three dimensional lattice gauge theory. We compare our simulation data with the predictions of both effective string models as well as perturbation theory. On the string side we obtain clear evidence for convergence of data to predictions of NambuGoto theory. We comment on the scales at which the static potential starts departing from one loop perturbation theory and then again being well described by effective string theories. We also estimate the leading corrections to the oneloop perturbative potential as well as the NambuGoto effective string. In the intermediate regions we find an empirical formula which gives surprisingly good fits.
Subleading properties of the QCD fluxtube in 3d lattice gauge theory
Pushan Majumdar^{†}^{†}thanks: Speaker.
Institut für Theoretische Physik, WestfälischeWilhelms Universität, Münster
Email: pushan@unimuenster.de
\abstract@cs
1 Introduction
Confinement of quarks, at least on the lattice seems to be due to the formation of a flux tube between a quark and an antiquark in the QCD vacuum [1]. It has been conjectured that the properties of this flux tube can be described by an effective hadronic string [2].
On the lattice, this flux tube can be observed by measuring the potential between a static quark and an antiquark. One of the characteristics of the string like behaviour of the flux tube is the presence of a long distance term in the potential in all dimensions and with a universal coefficient. It is known as the Lüscher term [3, 4].
The Lüscher term has been looked at in lattice simulations since the eighties [5, 6]. However in recent times, increase in computing power and improvement in algorithms have allowed really precise measurements of that term and also the next subleading term. See [7, 8, 9, 10, 11, 12, 13] for example.
Another aspect of the string behaviour is the excitation spectrum of the flux tube. A recent review on this topic can be found in [14].
In this article we present results for the continuum limit of our simulations of the Polyakov loop correlators for YangMills theory and compare the resulting static potential with both perturbation theory and string model predictions. This allows us to narrow down bounds on the distance beyond which we can say the flux tube indeed shows a string like behaviour.
2 Simulation parameters
We carried out simulations of three dimensional lattice gauge theory on lattices at four different lattice spacings. On these lattices the scale was set by the Sommer parameter , implicitly defined by , where is the force between the static quark and the antiquark. Our measured values were at respectively. We also observed that ( string tension) is a constant to a very good approximation as expected in the continuum limit. With this scale our coarsest lattice had a spacing of slightly below 0.13 fm and our finest lattice spacing is about 0.045 fm. We used symmetric cubic lattices and the Wilson gauge action. On all these lattices, we computed Polyakov loop correlators for various spatial separations , with separations being taken along the axes only.
3 Results
From the correlator one can extract the static quarkantiquark potential by We prefer to look at the first and the second derivative of this potential for the force between the quark and the antiquark and information about subleading terms. To facilitate our comparison with string models we actually compute a scaled second derivative which we call . This quantity is expected to become the Lüscher term () asymptotically.
On the lattice these are given by
(3.0) 
where and are defined as in [9] to reduce lattice artifacts.
For large distances, we are going to concentrate on the potential due to the NG string, the so called Arvis potential [16] given by We will define the L.O. and N.L.O. approximations by retaining and terms in the potential respectively. We will compare our lattice data on force and with predictions from leading order, NLO and from the full Arvis potential.
At short distances we compare with the perturbative potential obtained by Schröder [17] as with . For , . The perturbative force and can be computed by and
In Fig. 1 we plot versus . The horizontal line is and defines the Sommer scale . The N.L.O. curve lies in between the L.O. and Arvis curves and has been omitted for the sake of clarity. The data starts departing from the oneloop perturbative curve around 0.22 or 1.8 GeV and joins onto the string curves around 1.5 or 260 MeV. The scaling exhibited by the data is very good with all the four different beta values falling on the same curve. Beyond 1.5 it is virtually impossible to distinguish the different theoretical curves as they are all dominated by the string tension. The force data in fact gives the impression of the string description being good even at distances as small as .
does not contain the string tension, and has a universal value in the L.O. It is therefore more sensitive to the subleading behaviour of the flux tube. In Fig. 2 we plot in units of . We look at a wide range of starting from where the data almost touches the perturbative curves going all the way to the region where the string predictions hold.
3  2.808  0.10274 (3)  0.06172 (1)  0.05095 (1)  

4  3.838  0.11886 (8)  0.07737 (3)  0.06530 (3)  
5  4.876  0.1265 (2)  0.1052 (1)  0.08937 (6)  0.07714 (5) 
6  5.903  0.1299 (4)  0.1128 (2)  0.0985 (1)  0.0868 (1) 
7  6.920  0.1322 (9)  0.1183 (3)  0.1060 (5)  0.0946 (2) 
8  7.932  0.1342 (10)  0.1222 (3)  0.1118 (7)  0.1011 (2) 
9  8.941  0.1333 (6)  0.1251 (6)  0.1156 (2)  0.1064 (2) 
10  9.948  0.1340 (10)  0.1273 (8)  0.1190 (3)  0.1108 (3) 
11  10.95  0.1375 (34)  0.1288 (11)  0.1217 (4)  0.1145 (4) 
12  11.96  0.1296 (7)  0.1241 (5)  0.1174 (5)  
13  12.96  0.1308 (9)  0.1258 (7)  0.1201 (6)  
14  13.96  0.1323 (9)  0.1270 (6)  
15  14.97  0.1332 (13)  0.1282 (9)  
16  15.97  0.1345 (21)  0.1288 (11)  
17  16.97  0.1303 (15)  
18  17.97  0.1308 (23) 
The data almost lie on top of each other exhibiting nice scaling behaviour as one goes to larger values of . The and data come together already in the range 0.5 and 1.25 . The set joins onto this at around 1.5 and even the data joins up at around 2.25. This points to the possibility that the continuum limit of the scale where the flux tube is well described by the Arvis curve can be obtained even on relatively coarse lattices.
We estimate a temporal extent correction factor of about 0.1% for for our largest values. This is about an order of magnitude lower than our statistical errors at such values. Corrections due to finite spatial extents are of similar magnitude.
In , IR divergences prevent computation of the perturbative potential beyond one loop^{1}^{1}1NDH wishes to thank G. ‘t Hooft for an illuminating discussion on this.. We tried to obtain these terms from our data. Assuming the perturbative potential to be of the form
(3.0) 
we obtain and . Also we estimate the range of validity of first order perturbation theory to be about fermi (consistent with our estimate from the force data).
From our data we find that the string tension constitutes 95% of the force at around 1.3, 98% at around 2.1 and 99% at around 2.9. The relative difference between the Arvis and the leading order force, which gives us an idea about the importance of the subleading behaviour, is about 2% at 1.02, 1% at 1.2 and goes down to 0.1% at about 1.9.
The type of the string is determined by the subleading behaviour of the flux tube as at leading order, a variety of theories yield the universal Lüscher term [18, 19]. It is clear from the data that the approach to the Lüscher term is from below, consistent with effective bosonic string model predictions while at short distances the data matches perturbation theory. Beyond a distance of about 2.75, the data seems to be well described by the Arvis curve.
An interesting question is what happens to this distance for as increases. For on coarsest lattices [9] it is about 1.8. Similar indications are also there for [12]. However this scale shifts towards larger as one approaches the continuum limit.
Beyond N.L.O., effective string theories parametrise to be of the type . Our data gives and , which is quite different from the Arvis values. It is therefore interesting to know the effective string predictions for these coefficients.
At intermediate distances over a wide region of varying from to , the data is very well described by a formula of the type
(3.0) 
Existing data in 3d [9] also admit a similar description. The curves are shown in Fig. 3. It is not clear if there is any theoretical basis for such a description, but at least they provide accurate interpolation formulae.
4 Conclusions
In this article we have looked at the continuum limit behaviour of the flux tube at intermediate distances by measuring the static potential. Starting from a distance of about 0.1 fermi where the potential starts breaking away from 1loop perturbation theory, we go to distances of about 1.4 fermi where the data begins to be well described by the Arvis potential.
Our data on seems to approach the Lüscher term from below as expected in bosonic string models. At distances below 0.15 fermi the data joins onto the perturbative values. An empirical formula describes the data well in the intermediate region.
Acknowledgments
One of the authors, PM, gratefully acknowledges the numerous discussions with Peter Weisz. NDH thanks the Hayama Centre for Advanced Studies for hospitality. The simulations were carried out on the teraflop Linux cluster KABRU at IMSc as part of the Xth plan project ILGTI.
References
 [1] G. S. Bali, QCD forces and heavy quark bound states, Phys. Rep. 343 (2001) 1 [hepph/0110254].

[2]
P. Goddard et. al., Quantum dynamics of a massless relativistic string,
Nucl. Phys. B 56 (1973) 109 ;
J. Polchinski, A. Strominger, Effective string theory, Phys. Rev. Lett. 67 (1991) 1681;
J.M. Drummond, Universal subleading spectrum of effective string theory, hepth/0411017;
N.D. Hari Dass, P. Matlock, Universality of correction to Luescher term in PolchinskiStrominger effective string theories, hepth/0606265;
N.D. Hari Dass, P. Matlock, Field Definitions, Spectrum and Universality in Effective String Theories, hepth/0612291.
 [3] M. Lüscher, K. Symanzik, P. Weisz Anomalies of the Free Loop Wave Equation in the WKB Approximation, Nucl. Phys. B 173 (1980) 365 .
 [4] M. Lüscher, Symmetry Breaking Aspects of the Roughening Transition in Gauge Theories, Nucl. Phys. B 180 (1981) 317 .
 [5] J. Ambjorn, P. Olesen, C. Peterseon, Observation of a string in threedimensional SU(2) lattice gauge theory, Phys. Lett. B 142 (1984) 410 .
 [6] P. de Forcrand, G. Schierholz, H. Schneider, M. Teper, The string and its tension in SU(3) lattice gauge theory: Towards definitive results Phys. Lett. B 160 (1985) 137 .

[7]
M. Caselle, R. Fiore, F. Gliozzi, M. Hasenbusch, P. Provero, String effects in the
Wilson loop: A High precision numerical test, Nucl. Phys. B 486 (1997) 245
[heplat/9609041];
M. Caselle, M. Panero, P. Provero, String effects in Polyakov loop correlators, JHEP 0206 (2002) 061 [heplat/0205008];
M. Caselle, M. Hasenbusch, M. Panero, String effects in the 3d gauge Ising model, JHEP 0301 (2003) 057 [heplat/0211012];
M. Caselle, M. Hasenbusch, M. Panero, Short distance behavior of the effective string, JHEP 0405 (2004) 032 [heplat/0403004];
M. Caselle, M. Pepe, A. Rago, String effects in SU(2) lattice gauge theory, Nucl. Phys. Proc. Suppl. 129 (2004) 721 [heplat/0310005];
M. Caselle, M. Pepe, A. Rago, Static quark potential and effective string corrections in the (2+1)d SU(2) YangMills theory., JHEP 0410 (2004) 005 [heplat/0406008].  [8] M.Lüscher, P. Weisz, Locality and exponential error reduction in numerical lattice gauge theory, JHEP 0109 (2001) 010 [heplat/0108014].
 [9] M. Lüscher, P. Weisz, Quark confinement and the bosonic string, JHEP 0207 (2002) 049 [heplat/0207003].
 [10] N.D. Hari Dass, P. Majumdar, Stringlike behaviour of 4D SU(3) YangMills flux tubes, JHEP 0610 (2006) 020 [heplat/0608024].

[11]
P. Majumdar, The String spectrum from large Wilson loops,
Nucl. Phys. B 664 (2003) 213 [heplat/0211038];
P. Majumdar, Continuum limit of the spectrum of the hadronic string, heplat/0406037.  [12] H. B. Meyer, Static forces in d=2+1 SU(N) gauge theories, Nucl. Phys. B 758 (2006) 204 [heplat/0607015].

[13]
J. Juge, J. Kuti, C. Morningstar,
Excitations of the static quark antiquark system in several gauge theories,
heplat/0312019;
J. Juge, J. Kuti, C. Morningstar, QCD string formation and the Casimir energy, heplat/0401032 .  [14] J. Kuti, Lattice QCD and String Theory, PoS(LAT2005)001 [heplat/0511023].
 [15] N.D. Hari Dass, P. Majumdar, Continuum limit of string formation in 3d SU(2) LGT, heplat/0702019.
 [16] J. F. Arvis, The exact q antiq potential in Nambu string theory, Phys. Lett. B 127 (1983) 106 .
 [17] Y. Schröder, The Static potential in QCD., Ph.D. Thesis, DESYTHESIS1999021, Jun 1999.
 [18] K. Dietz, T. Filk, On The Renormalization Of String Functionals, Phys. Rev. D 27 (1983) 2944 .
 [19] S. Naik, Improved heavy quark potential at finite temperature from Antide Sitter supergravity , Phys. Lett. B 464 (1999) 73 [hepth/9904147].