Bosonic stringlike behavior and the ultraviolet filtering of QCD
Abstract
The gluonic action density is calculated in static mesons at finite temperature just below the deconfinement point. Our focus is to elucidate the role of vacuum ultraviolet fluctuations which are filtered using an improved smearing algorithm. In the intermediate source separation distance, where the free string picture poorly describes the flux tube width profile, we find upon reducing the vacuum action towards the classical instanton vacuum, the characteristics of the flux tube converge and compare favorably with the predictions of the free bosonic string. This result establishes a connection between the free string action and vacuum gauge fields and reveals the important role of ultraviolet physics in understanding the lattice data at this temperature scale. As a byproduct of these calculations, we find the broadening of the QCD flux tube to be independent of the ultraviolet filtering at large distances. Our results exhibit a linearly divergent pattern in agreement with the string picture predictions.
pacs:
12.38.Gc, 12.38.Lg, 12.38.AwI Introduction
In the dual superconductor scenario of quark confinement, the quantum chromodynamical (QCD) vacuum squeezes the color fields into a confining string dual to the Abrikosov line by the dual Meissner effect. The string conjecture Luscher et al. (1980) follows as an intuitive realization of this squeezed color field with the major objective of deriving the leading and subleading properties of the flux tube in the infrared region of confining gauge theories. This effective description is expected to hold on distance scales larger than the intrinsic thickness of the flux tube Caselle et al. (1996) in the rough phase of lattice gauge theories (LGT). The linearly rising potential part arises from the classical configuration of the string, and the quantum fluctuations of the string lead to the presence of a long distance term in the potential well known as the Lüscher term. The existence of the subleading term has been verified in high precision measurements of Polyakov loop correlators in the SU(3) gauge group at zero temperature Luscher and Weisz (2002). The fluctuations of the string render an effective width for the fluxtube which grows logarithmically Luscher et al. (1981) as the color sources are pulled apart. The logarithmic divergence has been verified in many lattice simulations corresponding to a variety of confining gauge models Bali et al. (1995); Pennanen et al. (1997); Caselle et al. (1996); Gliozzi et al. (2010a).
At high temperature, higherorder gluonic modes are present. The corresponding free bosonic string predicts a new set of measurable thermal effects. These include a decrease in the effective string tension Gao (1989); de Forcrand et al. (1985); Pisarski and Alvarez (1982), a change in the pattern of the tube’s growth in width from a logarithmic divergence into a linear divergence Allais and Caselle (2009), and a nonconstant width profile Allais and Caselle (2009); Bakry et al. (2010) along the line.
Unlike the situation at zero temperature, the thermal behavior of the free string manifests only at source separation distance scales larger than what one expects normally in the zero temperature regime Kaczmarek et al. (2000); Allais and Caselle (2009); Bakry et al. (2010); Dass and Majumdar (2008). The fact that the lattice data are poorly described by the free theory in the intermediate distance regime has been a subject of analytic and numerical studies which include higherorder terms of the effective string’s action Aharony and Karzbrun (2009); Luscher and Weisz (2004) into the corresponding partition function. The consequences of such an approach have been studied on the level of the potential Caselle et al. (2005, 2003) and, recently, has been extended to the fluxtube width profile Gliozzi et al. (2010b, a); Caselle (2010). Other studies investigate a possible finite intrinsic thickness of the QCD fluxtube Vyas (2010).
Apart from the linearly rising potential, the interesting physics of the effective confining string is mainly due to its quantum fluctuations. As we will see, remarkable features arise when the ultraviolet (UV) part of the fluctuations of the stringlike flux tube is filtered out for intermediate quark separations at high temperatures. At this distance scale it is not yet clear if the deviations from the string picture are due to a nonNambuGoto action or the fact that a stringlike behavior has not yet set in. Thus, it is interesting to address this problem in a variant context by reporting an observation regarding the role played by the UV fluctuations of the vacuum in these discrepancies. We do this by tracking the response of the QCD vacuum, which is subject to UV filtering, to the presence of external static color sources. This work extends the region for which the free string picture is of utility.
In the following we measure the gluonic actiondensity distribution by correlating an action density operator to Polyakov loop correlators. Measurements are taken on a set of SU(3) pure gauge configurations. The configurations are generated using the standard Wilson gauge action on two lattices of a spatial volume of and temporal extents of and , corresponding to temperatures and respectively. The simulations are performed for coupling value . At this value the lattice spacing is fm to reproduce the standard value of the string tension Mev Bali and Schilling (1993). The Monte Carlo updates are implemented with a pseudoheatbath algorithm Cabibbo and Marinari (1982) using FabriciusHaan and KennedyPendelton (FHKP) Fabricius and Haan (1984); Kennedy and Pendleton (1985) updating. Each update step consists of one heat bath and 5 overrelaxations. The measurements are taken on 500 bins separated with 2000 updating sweeps. Averaging inside each bin is performed by taking 5 measurements separated by 70 updating sweeps. This leads to a hierarchical integration, that is apart from updating the last time slice, similar to implementing a onelevel Lüscher Weisz (LW) algorithm Luscher and Weisz (2002).
The measurements are taken after smoothing the gauge field by an overimproved stoutlink smearing algorithm Moran and Leinweber (2008). The value of the smearing parameters used are and . Smoothing the gauge field reduces the action towards the action minimum or the classical instanton solution de Forcrand (2007). The UV characteristics of the gauge fields can be characterized in terms of Dirac eigenmodes. For example, the number of overimproved stoutlink smearing sweeps used here has been calibrated to a given spectral cutoff in the spectral representation of the Dirac operator Ilgenfritz et al. (2008). The measurements are taken on sets of smeared gauge configurations with increasing levels of smearing. This way we are able to set the limit where the QCD vacuum response to the presence of an external static color sources asymptotically approaches the low energy free effective theory behavior.
Ii Quark–antiquark potential
At fixed temperature , the Monte Carlo evaluation of the quark–antiquark potential at each is calculated through the Polyakov loop correlators
(1) 
with the Polyakov loop given by
(2) 
In the string picture, the Polyakov loop correlator assumes the functional form of the partition function of the twodimensional bosonic string
(3) 
The vector maps the region into , with Dirichlet boundary condition , and periodic boundary condition along the time direction , . is the string action in the physical gauge Luscher and Weisz (2002)
(4) 
The action decomposes into the classical configuration, the fluctuation part, and the string higherorder self interactions. A leading order approximation can be made by neglecting the selfinteraction terms. is a complex parametrization of the world sheet, such that , and is the modular parameter of the cylinder. Solving the path integral of Eq. (4) Dietz and Filk (1983), and using Eq. (II), the quark antiquark potential reads
(5) 
is the Dedekind eta function
(6) 
and is a renormalization parameter.
The numerical evaluation of the quark antiquark potential, Eq. (II), using a four dimensional smearing scheme leads to a systematic ambiguity in regard to the transfer matrix interpretation which allows one to identify the expectation values of the Polyakov loop correlators with . We recourse, instead, to three dimensional smearing keeping the temporallinks unsmeared. The same smearing parameters as above are used.
Our approach is as follows:

We start with 500 configurations. Each configuration is smeared only altering the spatial directions of the links, using spatiallyoriented staples. This is depicted by shading of spatial links in Fig. 1.

For every 3D smeared configuration, an update sweep is applied. Each update sweep consists of a heatbath step and four overrelaxation steps on the lattice. Four or more overrelaxation steps provide the same result within errors. Each step proceeds as an update of every single link based on its neighbors, effectively one at a time as one sweeps across the lattice in all four directions; spatial and temporal as indicated in Fig. 1.

An update sweep and subsequent measurement is repeated three more times on each configuration of the ensemble. This results in four measurements one over each of the four ensembles, where the ensembles differ in how many update sweeps they have had after spatial smearing, i.e, 1 to 4.
The measurements proceed as follows:

The temporal links in each of theses newly created configurations are integrated out using a source sum of staples. The temporal link variables are replaced with the new link variable
(7) using the numerical linkintegration method of Ref. de Forcrand and Roiesnel (1985). It should be noted, that unlike the link integration implementation of Ref. Parisi et al. (1983) where pseudoheat bath hits are performed only on the temporal links, the pseudoheat bath updates we described above are performed on both the smeared space oriented links and the time oriented links. The temporal links, however, are integrated out by the numerical evaluation of the equivalent contour integral of Eq. (7) as detailed in Ref. de Forcrand and Roiesnel (1985).

Finally, the Polyakov loop correlators are calculated on each of the four one update sweep separated configurations. The result is averaged and binned as a single jackknife entry to avoid artificial error reduction. The 500 decorrelated bins are then averaged.
The Monte Carlo update step starts from a low action configuration in the spatial directions due to smearing. The above described update procedure brings in a newly updated timelink such that the effects of local action reduction that was only in the spatial torus takes place in the fourdimensional lattice. In this way, the UV filtering is implemented keeping the integration over the path integral Eq. (II) systematic, thus, preserving the transfer matrix interpretation.
In Fig. 2 the value of the potential measured on various levels of spatially smeared configurations, normalized to its value at fm, are plotted. Fig. 2(a) shows the numerical behavior of the data using the above described 3D smeared heat bath/overrelaxation driven updates. On the other hand, Fig. 2(b) shows the corresponding numerical behavior of the data measured on standard fourdimensional smeared configurations. The data corresponding to the unsmeared lattice and the string model predictions of Eq. (5) at are also included.
The discrepancies between the unsmeared lattice data and the free string model occur in the intermediate distances fm. The numerical results for the quark–antiquark potential evaluated on the 3D smeared updated configurations show an interesting behavior with respect to the number of smearing sweeps (see Fig. 2(a)). The data at large distances show no response to the filtering of the UV fluctuations of the gauge field. The data at intermediate separation distances converge for a large number of smearing sweeps. Moreover, the data approach the free string model predictions.
It is interesting to compare these results to those obtained from the fourdimensional smearing illustrated in Fig. 2(b), where 20 and 100 sweeps of smearing are compared. The results for 100 sweeps of smearing coincide very well with the convergence toward the string model predictions observed in Fig. 2(a). Thus, the smearing approach can be used as an efficient method for exploring the more demanding threepoint functions required to determine the distribution of gluon flux, and this is used in the following. Similar results are observed at .
Iii The Gluonic Profile
The transverse degrees of freedom of the stringlike flux tube render an effective width for the tube. The meansquare width of the free bosonic string is defined as
(8) 
The above integral can be solved analytically Allais and Caselle (2009); Gliozzi et al. (2010b). The width of the tube in dimensions reads
(9) 
where are Jacobi elliptic functions, and is the UV cutoff which has been generalized to be dependent on distances from the sources. This solution gives the mean square width at all the planes transverse to the quarkantiquark line, and hence, is describing the topological shape of the fluctuating flux tube and its dependence on the temperature as well as its evolution with the increase of color source separation. In the above formula, we assume a dependence of the ultraviolet cutoff, , on the position of the transverse planes, since as we will see, this quantity assumes different values near the quark sources as we fit the above formula to lattice data.
Using a modular transform Allais and Caselle (2009); Gliozzi et al. (2010b), Eq.(9), at , in four dimensions Caselle (2010) becomes
(10) 
which indicates linear growth of the tube’s width at large distance.
The width of the action density of the free bosonic string can be compared to the width of the action density of the corresponding flux tube of the lattice gauge theory. After constructing the coloraveraged infinitelyheavy staticmesonic state
subsequent measurement by an action density operator is taken at each point of the threedimensional torus at each corresponding Euclidean time slice for every source configuration. The action density operator is constructed via a highlyimproved threeloop improved lattice fieldstrength tensor BilsonThompson et al. (2003). The measurements taken are averaged over the time slices. A scalar field that characterizes the gluonic actiondensity distribution field can be then measured using the definition Bissey et al. (2007)
(11) 
where denotes averaging over configurations and lattice symmetries, and the vector refers to the spatial position of the flux probe with respect to some origin. To further suppress the statistical fluctuations, the density distributions have been symmetrised around all the symmetry planes of the tube.
Unless otherwise indicated, the measurements presented through out this section are taken on the 4D smeared configurations. The decrease in with the increase of the smearing sweeps is depicted in Fig. 3. The actiondensity asymptotically converges to a minimum for values around to sweeps of smearing. Cluster decomposition of the operators leads to away from the quarks.
A measurement of the width of the fluxtube’s action density may be taken through fitting the density distribution , Eq. (11), to a Gaussian of the form
(12) 
with in each selected transverse plane to quark axis by making use of the cylindrical symmetry of the tube.
The mean square width of the flux tube is defined as the second moment of the flux density with respect to the central line connecting the two quarks, ,
(13) 
with or denoting the position of the quark source. The width is normalized with respect to the ultraviolet cutoff of Eq. (9) according to
(14) 
The measured values of the meansquare width in the middle plane of the tube versus the source separation are plotted in Fig. 4. Similar plots at three consecutive transverse planes , , and to the line joining the two color sources are illustrated in Figs. 5 and 6. has been measured for each smearing level and value by fitting Eq. (9) to data points having fm. Good is obtained.
At large distances, in the middle plane of the flux tube, the tube shows a width broadening pattern for increasing that does not depend strongly on the corresponding smearing level. The data at large distances are increasing linearly in agreement with the string model predictions Gliozzi et al. (2010b, a); Caselle (2010). The UV effects, on the other hand, are manifest in the data points at shorter distances. The width of the fluxtube measured on the lowest smearing level, where the short distance physics is best preserved, is poorly described by the free bosonic string model at short distances. As higher smearing levels are considered, the subsequent removal of the short distance physics from the gauge sector regulates the fast rate of growth of the flux tube width. However, this does not continue uncontrollably. The data ultimately converge near sweeps of smearing, in accord with the saturation in the action density of Fig. 3 and the potential in Fig. 2(a). Moreover, the UVfiltered results converge to the free string predictions.
Table 1 summarizes the measured for fits of Eq. (9) for the fit range fm. With the increase of the number of smearing sweeps the returned values of the improves and becomes stable near the regime of the action saturation of approximately 100 sweeps.
No.Sweeps  

3.2  1.6  1.20  0.98  0.96 
To clarify this point further, we investigate the response of the QCD vacuum to the presence of infinitely heavy sources that are not constructed using smeared temporal links. Instead, the Polyakov lines are evaluated in the calculations of the flux strength, Eq. (11), using the numerical link integration Parisi et al. (1983) procedure of Ref. de Forcrand and Roiesnel (1985) for noise reduction. The temporal links have not been smeared in the evaluation of Polyakov loops, rather the temporal links have been integrated out. This time, the Polyakov loops are taken from the unsmeared configurations and correlated with the smeared action density. For our analysis, performed on 500 configurations, we observe that the data corresponding to the width profile of the flux tube measured on high levels of vacuum UV filtering do display similar behavior to the results in Fig. 5. Fig. 6 presents results for the plane . Again the results systematically converge with a large number of smearing sweeps and approach the string model predictions. Note that the evaluation of the correlation function using this method involves a three point correlation function which becomes noisy at large distances.
The shape aspects of the fluctuating free string are contained mainly in the second term of Eq. (9) and can be isolated by considering the difference in the mean square width at a given plane with respect to the central plane, . The measured value of , however, depends on the corresponding plane at which the lattice data is fit to Eq. (9). The value of the fit parameter is fixed for each plane using lattice data at large separations fm and fm.
The changes in with respect to central plane is plotted separately in Fig. 7. The measurements performed at adjacent planes of the fluxtube reveal that this quantity varies along different planes orthogonal to the string. Indeed, we don’t get a perfect match with the freestring profile unless we take such changes into account. The increase in the value of the UV cutoff is mostly obvious near the quarks and may indicate the importance of string self interactions near the boundary or the string interaction with the quark source itself. These interactions are switchedoff in the NambuGoto (NG) action and would appear when fitting the lattice data to the free string profile as a variation in the value of the UV cutoff along the tube.
This variation in has been included as part of the generalized string model solution of Eq. (9) illustrated by the curves in the following figures. A comparison of generalized string model predictions with the corresponding change in the flux tube’s width on the lattice (obtained from fits of Eq. (11)) is shown in Fig. 8 for a source separation in the intermediate distance fm, fm and at large distance fm.
Lattice data for each gauge smoothing level is also depicted in Fig. 8. The jackknife uncertainties associated with the change in the meansquared width of the tube reveal correlated errors between the adjacent planes. Only subtle changes are observed in the tube’s width along the transverse planes with respect to the central plane for the analysis performed on smeared gauge configurations of the lowest smearing level sweeps. The tube tends to exhibit larger curvatures as higher levels of gauge smoothing are considered. At large values of UV filtering, the tube profile converges and approaches the geometrical shape of the freebosonic string. At large distances, on the other hand, the flux tube displays a curved width profile which compares well with the bosonic string profile and is not affected by smearing, as is evident in Fig. 8 at fm. The string model curve has been calculated based on the changes in the second term of Eq. (9).
Iv Conclusion
The presence of a pair of static external sources in the QCD vacuum induces a response of an effective free bosonic string for source separations in the intermediate separation region, provided the short distance vacuum fluctuations are filtered out. The flux tube, measured as a correlation between the mesonic operator and the vacuum action density, is found to exhibit a broadening pattern and a transverse structure similar to the freebosonic string for measurements taken near the saturation in action density minimization under smearing even at intermediate distances. At large distances, the UV fluctuations do not affect the tube growth, which exhibits a linear divergent pattern consistent with the string model predictions.
In carrying out the potential calculations, we introduced a novel method for studying the effects of UV filtering of the QCD vacuum. This method avoids the ambiguities of performing unsystematic integrations due to smearing the temporal links, thus, preserving the transfer matrix interpretation. Instead threedimensional spatial smearing is combined with single pseudoheatbath driven updates. The numerical data of the potential obtained this way converge towards the string model predictions at large number of smearing sweeps. A comparison with fourdimensional smearing results reveals that any systematic effects associated with smearing the temporal links are subtle.
The analysis performed at short distances provides an extension of the QCD vacua where the free string picture is of utility. The infrared region of the UV filtered vacuum can be described merely on the basis of a free string picture in the intermediate distances as well as large distances. This fact is relevant and complements recent investigations including higherorder self interactions to match lattice results.
Acknowledgments
This research was undertaken on the NCI National Facility in Canberra, Australia, which is supported by the Australian Commonwealth Government. We also thank eResearch SA for generous grants of supercomputing time which have enabled this project. This research is supported by the Australian Research Council.
References
 Luscher et al. (1980) M. Luscher, K. Symanzik, and P. Weisz, Nucl. Phys. B173, 365 (1980).
 Caselle et al. (1996) M. Caselle, F. Gliozzi, U. Magnea, and S. Vinti, Nucl. Phys. B460, 397 (1996), eprint heplat/9510019.
 Luscher and Weisz (2002) M. Luscher and P. Weisz, JHEP0207 049 (2002), URL http://arXiv.org:heplat/0207003.
 Luscher et al. (1981) M. Luscher, G. Munster, and P. Weisz, Nucl. Phys. B180, 1 (1981).
 Bali et al. (1995) G. S. Bali, C. Schlichter, and K. Schilling, Phys. Rev. D 51, 5165 (1995).
 Pennanen et al. (1997) P. Pennanen, A. M. Green, and C. Michael, Phys. Rev. D56, 3903 (1997), eprint heplat/9705033.
 Gliozzi et al. (2010a) F. Gliozzi, M. Pepe, and U. J. Wiese (2010a), eprint 1006.2252.
 Gao (1989) M. Gao, Phys. Rev. D40, 2708 (1989).
 de Forcrand et al. (1985) P. de Forcrand, G. Schierholz, H. Schneider, and M. Teper, Phys. Lett. B160, 137 (1985).
 Pisarski and Alvarez (1982) R. D. Pisarski and O. Alvarez, Phys. Rev. D 26, 3735 (1982).
 Allais and Caselle (2009) A. Allais and M. Caselle, JHEP 01, 073 (2009), eprint 0812.0284.
 Bakry et al. (2010) A. S. Bakry et al., Phys. Rev. D 82, 094503 (2010).
 Dass and Majumdar (2008) N. H. Dass and P. Majumdar, Physics Letters B 658, 273 (2008), ISSN 03702693.
 Kaczmarek et al. (2000) O. Kaczmarek, F. Karsch, E. Laermann, and M. Lutgemeier, Phys. Rev. D 62, 034021 (2000).
 Aharony and Karzbrun (2009) O. Aharony and E. Karzbrun, JHEP 06, 012 (2009), eprint 0903.1927.
 Luscher and Weisz (2004) M. Luscher and P. Weisz, JHEP 07, 014 (2004), eprint hepth/0406205.
 Caselle et al. (2005) M. Caselle, M. Hasenbusch, and M. Panero, J. High Energy Phys. 03, 026. 25 p (2005).
 Caselle et al. (2003) M. Caselle, M. Panero, P. Provero, and M. Hasenbusch, Nucl. Phys. Proc. Suppl. 119, 499 (2003), eprint heplat/0210023.
 Gliozzi et al. (2010b) F. Gliozzi, M. Pepe, and U. J. Wiese (2010b), eprint 1002.4888.
 Caselle (2010) M. Caselle (2010), eprint 1004.3875.
 Vyas (2010) V. Vyas (2010), eprint 1004.2679.
 Bali and Schilling (1993) G. S. Bali and K. Schilling, Phys. Rev. D 47, 661 (1993).
 Cabibbo and Marinari (1982) N. Cabibbo and E. Marinari, Phys. Lett. B119, 387 (1982).
 Fabricius and Haan (1984) K. Fabricius and O. Haan, Phys. Lett. B143, 459 (1984).
 Kennedy and Pendleton (1985) A. D. Kennedy and B. J. Pendleton, Phys. Lett. B156, 393 (1985).
 Moran and Leinweber (2008) P. J. Moran and D. B. Leinweber, Phys. Rev. D77, 094501 (2008), eprint 0801.1165.
 de Forcrand (2007) P. de Forcrand, AIP Conf. Proc. 892, 29 (2007), eprint heplat/0611034.
 Ilgenfritz et al. (2008) E.M. Ilgenfritz, D. Leinweber, P. Moran, K. Koller, G. Schierholz, and V. Weinberg, Phys. Rev. D 77, 074502 (2008).
 Dietz and Filk (1983) K. Dietz and T. Filk, Phys. Rev. D 27, 2944 (1983).
 de Forcrand and Roiesnel (1985) P. de Forcrand and C. Roiesnel, Physics Letters B 151, 77 (1985).
 Parisi et al. (1983) G. Parisi, R. Petronzio, and F. Rapuano, Phys. Lett. B128, 418 (1983).
 BilsonThompson et al. (2003) S. O. BilsonThompson, D. B. Leinweber, and A. G. Williams, Ann. Phys. 304, 1 (2003), eprint heplat/0203008.
 Bissey et al. (2007) F. Bissey et al., Phys. Rev. D76, 114512 (2007), eprint heplat/0606016.