The Landau gauge gluon and ghost propagators in 4D SU(3) gluodynamics in large lattice volumes
Abstract
We present recent results on the Landau gauge gluon and ghost propagators in pure gauge theory at Wilson for lattice sizes up to corresponding to physical volumes up to . In particular, we focus on finitevolume and Gribovcopy effects. We employ a gaugefixing method that combines a simulated annealing algorithm with finalizing overrelaxation. We find the gluon propagator for the largest volumes to become flat at . Although not excluded by our data, there is still no clear indication of a gluon propagator tending towards zero in the zeromomentum limit. New data for the ghost propagator are reported, too.
HUEP07/47 ADP0711/T651
\ShortTitleGluon and ghost propagators in large volumes
\FullConferenceThe XXV International Symposium on Lattice Field Theory
July 30 – August 4, 2007
Regensburg, Germany
1 Introduction
Presently, there is an intensive exchange of results and opinions between groups carrying out analytical and numerical studies of infrared QCD. This ongoing research focuses in particular on the infrared behavior of the Landaugauge gluon and ghost propagators. The latter is intimately related to the confinement [1, 2, 3]. What makes analytical predictions possible also in the nonperturbative sector of the theory is the possibility to write down a hierarchy of DysonSchwinger equations (DSE) connecting propagators and vertices. Under rather mild assumptions the hierarchy can be truncated. However, not all assumptions have been thoroughly checked. Note that recently the full system of Landau gauge DSE has been solved without any truncations within the asymptotic infrared region with a power ansatz for all Green functions involved [4]. Numerically, the propagators can be studied from first principles in terms of Monte Carlo (MC) simulations of lattice QCD. It is worth to compare the lattice results with the asymptotic powerlike behavior, and with numerical DSE solutions found in finite volumes [5]. It is interesting to see whether there remain differences as the infinitevolume limit is approached.
On the lattice we approximate the gluon propagator as the MC average
(1) 
with the gluon field transformed into Fourier space. The lattice momenta with integer are related to their physical values by .
The ghost propagator in momentum space at nonzero is defined by double Fourier transformation
(2) 
Practically, this is done by an inversion of the FaddeevPopov (FP) matrix using a conjugategradient algorithm with plane waves as sources. The FaddeevPopov operator in terms of the Landau gaugefixed links is
(3) 
with ( are the GellMann matrices). The functions and are called dressing functions of the respective propagator.
In Landau gauge the gluon and ghost dressing functions are predicted to follow the simple power laws [6]
(4) 
in the asymptotic regime . Thereby, both exponents are related to some which, under the assumption that the ghostgluon vertex is infraredregular, takes a value of about [7]. That is, the gluon propagator is predicted to decrease towards lower momenta and to vanish at . At which scale this asymptotic behavior sets in cannot be concluded from those studies, however.
A truncated system of DSE formulated on a 4D torus and numerically solved [5] predicts a specific finitevolume behavior which, at the first glance, looks quite similar to earlier lattice results obtained in particular by some of us [8]. Characteristic deviations for momenta of the gluon and the ghost propagator from the momentum dependence of the respective propagators at infinite volume should be expected. In order to check the DSE predictions we decided to evaluate the gluon and the ghost propagator for increasingly large symmetric lattices. We have measured the gluon and ghost propagators for configurations generated with the Wilson gauge action at fixed on , , and lattices. Note that the latter corresponds to a volume of about . Comparing results from either analytic or numerical approaches for varying 4volume will hopefully allow (i) to conclude for which momenta data on both propagators are reliable, and (ii) to estimate the order of magnitude of distortion of the momentum dependence due to finitesize and Gribovcopy effects. This paper presents first results of this study.
2 Gauge fixing
To fix the Landau gauge, we apply to all links a gauge transformation () mapping , with the aim to maximize a gauge functional
(5) 
or, more exactly, to find the global maximum of [1]. In this work is considered as a periodic field on all lattice sites. To find the global maximum in practice is a complicated problem which becomes exceedingly time consuming with increasing lattice volume. Starting from an initial random gauge transformation one generally arrives at one of many local maxima of . The corresponding gaugefixed configurations are called Gribov copies. They all satisfy the differential gauge condition together with the additional necessary condition that the FaddeevPopov operator has a positive spectrum (apart from its 8 trivial zero modes). With increasing volume, the copies become dense with respect to the value of the functional (5) and the spectral density of the FP operator near zero grows [9].
One way to suppress the effect of the Gribov ambiguity is to find
local maxima of and to choose among them
the “best” one (“bc”), which possesses the largest value of
. The underlying idea is that the maximal value of the
local maxima approaches the global maximum of , and the
distortion of gaugedependent observables, computed on such
copies, vanishes in the limit . Such
studies normally use the overrelaxation (OR) technique to
search for the maximum of . They have been carried out in
[8, 10, 11, 12] and have shown the Gribovcopy
effect to become weaker with growing lattice extension – in
accordance with Zwanziger’s conjecture [1]. This
suggests that it is tolerable to restrict gaugefixing
computations on large lattices to one gauge copy
only
For the case we chose such that it leads to a sufficiently large mobility in the functional space. The final temperature was taken low enough that the subsequent OR was not slowed down while penetrating further functional barriers. This is witnessed by the check that the violation of the differential gauge condition, , monotonously decreases until the machine precision is reached (stopping criterion) in almost all cases. In practice, for and we restricted ourselves to one copy, and carried out from to heatbath (HB) sweeps of SA with 4 microcanonical sweeps after each HB one. We checked that the smaller steps are done in between, the higher the local maxima being reached. Finally, we note that a linear decrease in seems not to be the optimal choice. schedules with smaller steps close to and larger steps at the end (with fixed) lead to higher values (after completing the full SA procedure).
3 Ghost propagator results
The ghost propagator at is shown as a function of in Fig. 1; on the left hand side for a single configuration, simply comparing results after either OR or SA gauge fixing, and on the right hand side as an average over 14 configurations in the case of OR, and over 7 configurations using SA gauge fixing. The influence of the gaugefixing method, here through the emerging copy, can be seen only for the three lowest momenta at this lattice size. In general, one notices that the higher the gauge functional, the lower the estimates of the ghost propagator at the smallest momenta. This comparison (Fig. 1) demonstrates that the problem of Gribov copies does still exist for , resulting in maximally difference of ghost propagators at lowest momenta. Note that this result cannot be directly compared to our previous studies [10, 8, 11, 12] of Gribovcopies effect, in which the “fcbc” comparison was used to assess the Gribov ambiguity. A detailed check of Zwanziger’s conjecture on the weakening of the Gribov problem with an increase of the lattice volume (using the SA vs. OR “onecopy” comparison) requires further studies both for smaller and larger lattices.
In Fig. 2 a scatter plot is shown of the ghost dressing function for a broad range of momentum, combining data obtained for 7 configurations on a lattice, 14 configurations on a lattice, 3 configurations on a and 3 configurations on a lattice, all thermalized at . The gauge field configurations were produced with a heatbath algorithm applying O(1000) thermalization sweeps in between. We consider only momenta surviving a cylinder cut with [16]. For all lattice sizes GF has been carried out with the SA algorithm. Surprisingly, all the values for the ghost propagator fall perfectly on one universal curve (within accuracy), besides those for the 2 smallest momenta. The results, especially those found at , show that a true IR exponent cannot yet be defined or does not exist at all. This is at variance with the asymptotic DSE prediction [5] and also with motivated by thermodynamic considerations in [18]. The latter estimate of is based on the required cancellation of gluon and ghost contributions to the pressure, that otherwise were building up a StefanBoltzmann law, in the confinement phase. Note that downward deviations of the data at lowest momentum for each physical box size from the infinitevolume curve of are predicted by the DSE approach on a finite torus [5]. However, we do not find such deviations.
4 Gluon propagator results
Fig. 3 shows data for the gluon propagator computed for three different lattice sizes at . There, the gauge was fixed with the SA algorithm. At the present stage, the data favor a nonvanishing gluon propagator at zero momentum, as there is no sight of a different behavior even at the largest lattice volume available to us. Also, the decrease of the zeromomentum propagator upon increasing the volume seems to become less with bigger . The data for resemble the pattern of overshooting deviations from an universal function of momentum known from the DSE solutions on a finite torus [5], though. If the lowest two or three momenta were removed from the plot, the picture would be less convincing in favor of an emerging plateau. Better statistics and data on even larger symmetric lattices (with ) will make us more confident in this.
5 Discussion
We conclude that using the SA technique for the purpose of gauge fixing considerably facilitates simulations of ghost and gluon propagators in Landau gauge on large lattices. We find that SAbased computations of the ghost propagator seem to be less affected by statistical fluctuations compared to other calculations where OR is used. Additionally, estimates of the ghost propagator at low momenta are systematically lower than those obtained after simple OR.
A continuous decrease in slope in the ghost dressing function below 0.4 GeV does not conform to a simple powerlaw ansatz. Therefore, any attempts to extract infrared exponents from lattice data seem to be premature at the present stage. Qualitatively, the same behavior is seen for the case of (see [17]) and also in the DSE solutions on a torus [5]. However, the effects of finite volumes are much less than expected from there though. The same we find for the gluon propagator which we cannot confirm to be infrareddecreasing even at volumes larger than .
Future analytical (DSE and renormalization group) studies and lattice simulations, optionally including also flip operations into the GF procedure [11], at even larger volumes will hopefully help to resolve or explain the existing discrepancies between the lattice and analytical findings.
Acknowledgements
Simulations were done on the MVS15000BM and MVS50000BM at the Joint Supercomputer Centre (JSCC) in Moscow and on the IBM pSeries 690 at HLRN. This work was supported by joint grants DFG 436 RUS 113/866/0 and RFBR 060204014. Part of this work is supported by DFG under contract FOR 465 Mu 932/24, and by the Australian Research Council.
Footnotes
 This is called “first copy” (“fc”) in a multicopy approach.
References
 D. Zwanziger, Nucl. Phys. B412, 657 (1994); Phys. Rev. D69, 016002 (2004), [hepph/0303028].
 T. Kugo and I. Ojima, Prog. Theor. Phys. Suppl. 66 (1979) 1; T. Kugo, hepth/9511033.
 R. Alkofer and L. von Smekal, Phys. Rept. 353, 281 (2001), [hepph/0007355].

R. Alkofer, C.S. Fischer, and F.J. LlanesEstrada,
Phys. Lett. B611, 279 (2005), [hepth/0412330];
C.S. Fischer and J.M. Pawlowski, Phys. Rev. D75, 025012 (2007), [hepth/0609009].  C.S. Fischer, A. Maas, J.M. Pawlowski, and L. von Smekal, Ann. Phys. 322, 2916 (2007), [hepph/0701050].
 L. von Smekal, R. Alkofer, and A. Hauck, Phys. Rev. Lett. 79, 3591 (1997), [hepph/9705242].

Ch. Lerche and L. von Smekal,
Phys. Rev. D65, 125006 (2002),
[hepph/0202194];
D. Zwanziger, Phys. Rev. D65, 094039 (2002), [hepph/0109224].  A. Sternbeck, E.M. Ilgenfritz, M. MüllerPreussker, and A. Schiller, Phys. Rev. D72, 014507 (2005), [heplat/0506007].
 A. Sternbeck, E. M. Ilgenfritz, and M. MüllerPreussker, Phys. Rev. D73, 014502 (2006), [heplat/0510109].
 T.D. Bakeev, E.M. Ilgenfritz, V.K. Mitrjushkin, and M. MüllerPreussker, Phys. Rev. D69, 074507 (2004), [heplat/0311041].
 I.L. Bogolubsky, G. Burgio, V.K. Mitrjushkin, and M. MüllerPreussker, Phys. Rev. D74, 034503 (2006), [heplat/0511056]; I.L. Bogolubsky, et al., arXiv:0707.3611 [heplat].
 E.M. Ilgenfritz, M. MüllerPreussker, A. Sternbeck, A. Schiller, and I.L. Bogolubsky, Braz. J. Phys. 37, 193 (2007), [heplat/0609043].
 S. Kirkpatrick, C.D. Gelatt Jr., and M.P. Vecchi, Science 220, 671 (1983).
 V. Cerny, J. Opt. The. Appl. 45, 41 (1985).
 G.S. Bali, V. Bornyakov, M. MüllerPreussker, and F. Pahl, Nucl. Phys. Proc. Suppl. 42, 852 (1995), [heplat/0303028].
 D.B. Leinweber, J.I. Skullerud, A.G. Williams, and C.Parinello, Phys. Rev. D60, 094507 (1999), [heplat/9811027].
 A. Sternbeck, L. von Smekal, D.B. Leinweber, and A.G. Williams, PoS(LATTICE 2007) 340 (2007).
 M.N. Chernodub, V.I. Zakharov, [hepph/0703167].