arxiv: v2 [hep-lat] 13 Dec 2007

Transcription

1 Tests of the lattice inde theorem Gerald Jordan, Roman Höllwieser, and Manfried Faber Atomic Institute, Vienna Universit of Technolog, Wiedner Hauptstr. 8-, 4 Vienna, Austria Urs M. Heller American Phsical Societ, One Research Road, Bo 9, Ridge, NY 9-9, USA (Dated: August 9, 8) arxiv:7.5445v [hep-lat] 3 Dec 7 We investigate the lattice inde theorem and the localiation of the ero-modes for thick classical center vortices. For non-orientable spherical vortices, the inde of the overlap Dirac operator differs from the topological charge although the traces of the plaquettes deviate onl b a maimum of.5% from trivial plaquettes. This ma be related to the fact that even in Landau gauge some links of these configuration are close to the non-trivial center elements. PACS numbers:.5.ha,.38.gc I. INTRODUCTION At large separations static color sources of non-vanishing N-alit give rise to gluonic flu tubes and a linearl rising potential. This behavior is caused b the non-perturbative QCD-vacuum which compresses the gluonic flu. Center vortices [,, 3, 4, 5, ], quantied magnetic flu lines, pla an important role for the confinement of color charges, as has been shown b numerical calculations[7, 8]. In addition, simulations have indicated that vortices could also account for phenomena related to chiral smmetr, such as topological charge and spontaneous chiral smmetr breaking (SCSB) [9,,, ]. It is well-known that these non-perturbative features of QCD are intimatel linked to the properties of the low-ling spectrum of the Dirac operator. The topological charge of a gauge field equals the inde of the Dirac operator, while the chiral condensate, the order parameter for SCSB, is proportional to the spectral densit of the near-ero modes. With the advent of the overlap operator [3, 4], the fundamental problems of investigating chiral features on the lattice have been overcome. The overlap operator has an eact chiral smmetr [5], implements a lattice version of the inde theorem [3], and ma even be used for the definition of a local topological charge densit [9]. In this paper, we report on our calculations with the overlap operator applied to thick classical center vortices. We investigate the localiation of ero-modes with respect to the position of the thick vortices and find an interesting discrepanc in the topological charge determined b different methods. II. STRATEGY We workwith SU() gaugefields onlattices with N s sites in the -, -and-direction, andn t sites in the t-direction. In four dimensions vortices form closed two-dimensional surfaces. We investigate thick classical center vortices in the shapes of planes (closed b lattice periodicit) and spheres. The details of the individual vorte tpes will be discussed along with our results. We compare different definitions of the lattice topological charge: () The topological charge determined from the shape of P-vortices, that means from vorte intersections and writhing points according to ref. [9]. () The inde of the overlap Dirac operator [3, 4]. According to the Atiah-Singer inde theorem, which holds for continuous fields, the topological charge is given b the inde indd[a] = n n + = Q () Also at Photonics, Vienna Universit of Technolog, Gusshausstrasse 5-9, 4 Vienna, Austria. Electronic address: faber@kph.tuwien.ac.at

2 where n and n + are the number of left- and right-handed ero modes [, 7, 8]. It has been shown that this relation is also valid on the lattice for the overlap Dirac operator: its inde coincides with the discretied gluonic definition in the continuum limit [9]. The overlap Dirac operator is defined b [] Here, ǫ is the sign function, D = [+γ 5ǫ(H + L )] () H + L = γ 5D w ( m ), (3) and D w is the usual lattice Wilson Dirac operator with mass m (we use m =.). The overlap operator can be used to define the topological charge onl on so-called admissible gauge fields. This restriction assures that H + L has no ero eigenvalues so that the sign-function is well-defined. A sufficient, but not necessar condition for this is [, ] tr(½ U µν ) <.3 (4) The diagonaliation of the overlap operator ields the corresponding Dirac ero modes with well-defined chiralit properties. We further look at the distribution of the scalar fermionic densit ρ() = c,d v() cd, (5) where the summation indices c and d refer to color and Dirac indices of the eigenvectors v, respectivel, and the left- or right-handed chiral densities ρ ± () = v() ( γ 5 ) d,d cd v() cd, () c,d,d where onl ero-modes of the corresponding chiralit contribute. When the ero modes of a specific chiralit are degenerate, it is useful to investigate the square root of the sum of the scalar densities of these modes. This average densit respects eventual smmetries of the field configuration. (3) The integral (sum, on the lattice) of the gluonic charge densit q() = tr(f µν Fµν ) in the plaquette and/or hpercube definitions on the lattice, see ref. [3, 4]. Since Monte-Carlo configurations are in general too coarse, these definitions are usuall onl applied after cooling. In the continuum, since the topological charge densit q() is the total derivative of the topological current k µ [5], q() = µ k µ, k µ = 8 ǫ µνρσtr [A ν ρ A σ +i 3 ] A νa ρ A σ, (7) one can appl Gauss theorem to transform the epression for Q into a surface integral, Q = d 4 q() = dσ µ k µ, (8) S 3 provided the gauge field is smooth. Assuming a sufficientl fast decaing field strength for, the gauge field at infinit can be written as a pure gauge A µ ( ) = i µ ΩΩ. (9) It can be shown that the surface integral Eq. (8), and hence Q for smooth fields, measures the winding number of Ω, which maps the group space SU() onto the space-time surface S 3 at infinit.

3 3 Geometr Topological charge densit Fermionic densit frag replacements FIG. : Plane vortices on a 4 -lattice in - and t-plane intersect in four points giving topological charge and fermionic densit. III. DESCRIPTION OF CONFIGURATIONS AND RESULTS A. Plane vortices The construction of planar vortices is eplained in []. We use periodic boundar conditions where plane center vortices alwas come in pairs. All three definitions of Q ield identical results for all configurations containing onl plane vortices. As an eample, we show in Fig. some diagrams for two orthogonal pairs of plane vortices on a 4 lattice, which intersect in 4 places. The left picture shows the position of the P-plaquettes after going to maimal center gauge and center projection. One vorte pair etends in --planes, the other pair in -t-planes. The short pieces of lines attached to the vertical line smbolie the etent of the -t-vortices in the time-direction. The diagram in the middle of Fig. shows an equi-densit surface of the topological charge densit determined in the plaquette definition and the right diagram shows an equi-densit surface of the scalar fermionic densit. EachofthefourintersectionpointsgivesrisetoalumpoftopologicalchargeQ = ± [9]. Thesignofthecontribution at a given intersectionpoint canbe changedb a flip ofthe orientationofthe vortesurface, this means b a transition from links U µ () creating one of the vortices to U µ() in the region of the intersection. After abelian projection a smooth transition of the vorte orientation in color space produces a monopole line. In this wa one can produce configurationswith two vortepairs oftopological chargeq =,± or ±. In all these cases the lattice inde theorem is fulfilled and our results agree with the analtical solution for the ero modes presented in [7]. B. A non-orientable spherical vorte A non-orientable spherical vorte can be constructed with an appropriate profile -function α( r r ) b time-like links which form a hedgehog in one time-slice of the lattice, where r is the center of a spatial sphere { ep{iα( r r ) n() σ} t =, µ = 4, U µ () = () ½ elsewhere. A characteristic propert of a hedgehog configuration is the agreement between the color direction n and the spatial direction of r r n( r,t) = r r r r, () see left diagram in Fig.. For the choice () of the links the time-like Wilson lines (holonomies) agree with the time-like links in the t = time-slice. At r = r and at the largest possible distances from r on a periodic lattice the direction of the color-vectors () are undefined. This does not lead to a singularit of the gauge field if the profile-function α(r) is appropriatel

4 4 link phase Sfrag replacements link phase / / / link phase 3 4 FIG. : Thick Spherical SU()-vorte (hedgehog, non-orientable) and change of its link phase α(r) = α (r) according to Eq. (). chosen such that the Wilson lines at these special points become center elements of the gauge group. A good choice is α(r) = α (r) with ( ) r < R α (r) = r R R < r < R+ () R+ < r. As shown in the right diagram of Fig. for a 4 3 lattice the phase α of the links in the first time slice varies from at the center of the sphere to outside of the vorte core of thickness, chosen as = 8. The check that this configuration is a vorte is done with maimal center gauge and center projection and results in a P-vorte forming a lattice representation of a sphere of radius in the first time slice. The color structure of the thick vorte which is smbolicall indicated in the left diagram of Fig. 3 leads to a monopole loop on a great circle of the P-vorte after maimal abelian gauge and abelian projection. The direction of the loop depends on the U() subgroup chosen as abelian degrees of freedom. For the subgroup defined b the Pauli matrices σ,σ or σ 3 the monopole loops are in the --, -- and --plane, respectivel. This is indicated schematicall in Fig. 3 with three colors. FIG. 3: The color structure of the vorte surface (l) leads to monopole lines after abelian projection (r). The P-vorte surface which is a closed surface in dual space consists onl of the duals of space-time plaquettes. There are no vorte intersections and no writhing points. Therefore, the topological charge Q determined from the P-vorte surface [9] vanishes, Q =. The traces of all plaquettes are close to unit, tr(½ U µν ) cos 8 =.5,

5 5 Polakov phase Sfrag replacements / 3 4 Polakov phase / FIG. 4: Phases of Wilson lines along a line through the center of a spherical vorte (l), normalied scalar densities for all ero-modes of one chiralit in the planes with = close to the center of the sphere (r). The densit vanishes in the region where the phase of the Polakov loop is. Both chiralities and both time-slices on the 4 3 -lattice give the same picture. compare the link-phases in Fig.. The topological charge in the plaquette definition or in the hpercube definition can be used safel without cooling. According to Eq. () the onl non-trivial links are time-like, therefore all space-space plaquettes are trivial, which confirms the vanishing of the topological charge. This configuration has three ero-modes [9] of positive chiralit and four ero-modes of negative chiralit. Normalied densities of positive and negative chiralit modes have the same spatial distribution. For the spherical vorte with center at (9.5,9.5,9.5), we show the scalar densit in the plane = in Fig. 4. The non-trivial links connect = and =, and the densities in both time-slices are identical. The ero-modes are localied mainl in the region of trivial Wilson lines (for r r ) and seem compressed in the regions where the mirror pictures of the spherical vorte on the periodic lattice approach each other. In the region where the Wilson lines are non-trivial center elements (for r r around ), the densit vanishes. According to the inde theorem the topological charge is given b Q = n n + =. (3) This value is in disagreement with the values determined above (Q = ). But it agrees with the gluonic result which one gets after cooling. The cooling histor is shown in Fig. 5. The left scale shows the action S = ( TrU ) in units of the one-instanton action S inst and the right scale the value of the topological charge in the plaquette definition. The spherical vorte is not a minimum of the action, i.e. not an instanton. During cooling its action is first rapidl decreasing, whereas its topological charge starts with ero and reaches a value close to one after a few cooling steps. On the 4 3 -lattice the action does not reall approach a plateau value. Since a plateau appears for a 4 -lattice this seems to be related to the strong asmmetr of the lattice. Hence one would epect to see a clear plateau on a smmetric, but finer 4 4 lattice, as well. The number of ero-modes is the same for various lattice sies we used: 4, and 3 up to 4 3, increasing the spacial etent L s in steps of 4. We conclude that cooling and overlap fermions give a different topological charge than the discretied integral d 4 tr(f µν Fµν ) even for ver smooth field configurations. The hedgehog configuration is characteried b links in the center of the spatial sphere which are non-trivial center elements and Wilson lines L( ), = (,,), with a topological winding number ν =, ν({l}) = 4 d 3 rǫ ijk Tr(V i V j V k ), V i = i LL. (4) The topological charge that we get in the cooling process and from the inde theorem agrees with this winding number. In the continuum limit, the gauge field at the center of the hedgehog becomes singular. This singularit is due to the topological obstruction and can t be removed b a gauge transformation. With Landau-gauge, which we did for gauge copies on a 4 3 -lattice, it is not possible to get rid of links close to non-trivial center elements. For the lattice spacing a approaching ero, the gauge field for these links, A () σ 3 /(ga), diverges. In all gauge copies we found the same links with negative traces. The number of such negative links and the values of their traces are listed in Table I.

6 action S/Sinst 5 5 Q S.5 topological charge Q number of cooling steps FIG. 5: B cooling the topological charge rises from ero close to one for α = α (right scale) while the action S (in units of the one-instanton action S inst) decreases slowl (left scale). The data is for a 4 3 -lattice. number of links truµ() TABLE I: Number of links with negative traces and trace values after Landau gauge for the non-orientable spherical vorte of Eq. () on a 4 3 -lattice, e.g. there are eight links with truµ() =.93. C. Other configurations with spherical vortices An orientable spherical vorte we construct with the prescription { ep{iα( r r ) n i () σ i } t =,µ = 4, U µ () = (5) ½ elsewhere. Here, n i ( r,t) = ( r r)i r r is the absolute value of the coordinates of the direction vector n. These links cover onl half of the SU() group space and this half the cover twice. The time-like Wilson lines have topological winding number ero. This configuration has no ero-modes and all was to determine the topological charge lead to the same result, Q =. With maimal Abelian gauge no monopole lines appear. With Landau-gauge, which we did for gauge copies on a 4 3 -lattice, one can get rid of links close to non-trivial center elements. The number of such negative links and the values of their traces are listed in Table II. number of links truµ() TABLE II: Number of links with negative traces and trace values after Landau gauge for the orientable spherical vorte of Eq. (5) on a 4 3 -lattice. A second tpe of a non-orientable spherical vorte we construct with the prescription () and the profile function α(r) = α + (r), ( ) r < R α + (r) = + r R R < r < R+ () R+ < r The phases of the links in the time slice agree again with the phases of the Wilson lines in time direction and are shown in Fig.. In the center of the spherical vorte the temporal Wilson lines are trivial and at infinit the are

7 non-trivial center elements. This configuration gives onl one ero-mode of negative chiralit and none of positive chiralit. The ero-mode is localied in the center of the spatial sphere and in both time-slices occupied b the vorte. Again, the scalar densit prefers regions of trivial Wilson lines over those of non-trivial values. According to the inde theorem the topological charge is - and agrees with the winding number (4) of the hedgehog of Wilson lines. The topological charge computed from plaquettes is Q = while the plaquette traces differ b less than.5 % from trivial ones. Since this configuration is related to the configuration in Fig. 4 b a non-periodic gauge transformation, it is not surprising that again the different determinations of the topological charge do not agree. 7 Polakov phase / 3 4 Polakov phase / FIG. : Phase-diagram of the Polakov loop (l) along a line through the center of a spherical vorte with the profile function α = α + (l) and scalar densit of the onl fermionic ero mode (r) which has positive chiralit. For two equal vortices of this tpe, α = α +, one in time slice t = and another at t = of the 4 3 -lattice we get three ero-modes of positive chiralit with a scalar densit like the diagram in Fig. 4 and five of negative chiralit with the densities shown in Fig. 7. Since the modes are degenerate, the can mi and the plots displa the sum of the densities of all ero-modes. The left diagram displas the t = slice, where onl four of these modes contribute significantl, while the fifth mode leads onl to the small peak in the middle. Conversel, the densit in the t = slice (shown in the right picture) is dominated b the fifth mode in the center, the other modes being of a magnitude invisible on this scale. The phase of the Polakov line for this configuration is ero in the center of the spherical vortices and due to the -periodicit in the region between the mirror pictures on the periodic lattice. Again the ero-modes are localied in the region with trivial time-like Wilson lines PSfrag replacements t = t = FIG. 7: Scalar densities of fermionic ero modes of negative chiralit for two equal vortices with the profile function α = α +, in the time slices t = and t = of the 4 3 -lattice.

8 8 A configuration with two vortices of different tpe in the two time-slices, one with α = α at t = and another with α = α + at t = does not show an ero-mode. In this case there is no discrepanc in the determinations of the topological charge. The phase of the Wilson line, which is the sum of the phases from the two vortices, is constant and equal to. A summar of the the number of ero-modes for the configurations described above with spherical vortices is given in Table III. profile α n + n α at t = 3 4 α + at t = α at t = and α + at t = 3 5 α + at t = and α at t = 5 3 α + at t = and α + at t = TABLE III: Number of ero-modes for various smooth field configurations with spherical vortices. The profile functions α ± are defined in Eq. () and (). The above table can be summaried in the following empirical rule: Spherical vortices in slices (3D volumes) of the lattice contribute to the Dirac operator inde with an integer given b the winding number (4) of the corresponding Wilson lines, mapping the 3D volume of the slice to the SU() manifold of the Wilson lines. IV. CONCLUSION We have investigated various configurations with plaquettes close to unit. For configurations with a non-ero winding number of the Wilson lines we get disagreement between the topological charge determined b the plaquette and hpercube definition, and from the vorte intersections and writhing points on one side, and the topological charge from the inde theorem and cooling on the other side. The transition in Eq. (8), from the 4D integral over the topological charge densit, to the 3D integral of the SU()-valued map Ω, is onl possible for smooth gauge fields. Vorte configurations with non-trivial winding numbers have singularities in the gauge fields which can t be removed b gauge transformations or smooth deformations of the field configuration. According to the vorte model of confinement vortices form closed two-dimensional surfaces which in the confined phase penetrate the whole space-time. In ref. [8] it was shown that P-vorte surfaces are mostl non-orientable. It is ver likel that this non-orientabilit is related to a non-ero winding number of the gauge field and the artificial configurations which we discussed in this article reflect properties of those field configuration which saturate the path integral. The discrepancies which we have indicated are not due to the coarseness of the lattice, the traces of all plaquettes are close to unit, tr(½ U µν ) cos 8 =.5 but probabl due to singularities in the gauge field related to non-trivial mappings of Wilson lines. Acknowledgment This stud was partiall supported b the Austria Science Fund (FWF) under grant P-N. [] G. t Hooft, Nucl. Phs. B38, (978). [] P. Vinciarelli, Phs. Lett. B78, 485 (978). [3] T. Yonea, Nucl. Phs. B44, 95 (978). [4] J. M. Cornwall, Nucl. Phs. B57, 39 (979). [5] G. Mack and V. B. Petkova, Ann. Phs. 3, 44 (979). [] H. B. Nielsen and P. Olesen, Nucl. Phs. B, 38 (979). [7] L. Del Debbio, M. Faber, J. Greensite, and S. Olejnik, Phs. Rev. D55, 98 (997), hep-lat/95. [8] T. G. Kovacs and E. T. Tomboulis, Phs. Rev. D57, 454 (998), hep-lat/979. [9] H. Reinhardt and M. Engelhardt, Center vortices in continuum ang-mills theor, in Quark Confinement and the Hadron Spectrum IV, edited b W. Lucha and K. M. Maung, pp. 5, World Scientific,, hep-th/3. [] P. de Forcrand and M. D Elia, Phs. Rev. Lett. 8, 458 (999), hep-lat/99. [] C. Aleandrou, P. de Forcrand, and M. D Elia, Nucl. Phs. A3, 3 (), hep-lat/9995.

9 [] M. Engelhardt, Nucl. Phs. B38, 8 (), hep-lat/4. [3] R. Naraanan and H. Neuberger, Nucl. Phs. B443, 35 (995), hep-th/948. [4] H. Neuberger, Phs. Lett. B47, 4 (998), hep-lat/977. [5] M. Luscher, Phs. Lett. B48, 34 (998), hep-lat/98. [] M. F. Atiah and I. M. Singer, Annals Math. 93, 39 (97). [7] A. S. Schwar, Phs. Lett. B7, 7 (977). [8] L. S. Brown, R. D. Carlit, and C.-k. Lee, Phs. Rev. D, 47 (977). [9] D. H. Adams, J. Math. Phs. 4, 55 (), hep-lat/9. [] R. G. Edwards, U. M. Heller, and R. Naraanan, Nucl. Phs. B54, 457 (999), hep-lat/9877. [] M. Luscher, Nucl. Phs. B549, 95 (999), hep-lat/983. [] H. Neuberger, Phs. Rev. D, 855 (), hep-lat/994. [3] P. Di Vecchia, K. Fabricius, G. C. Rossi, and G. Veneiano, Nucl. Phs. B9, 39 (98). [4] P. Di Vecchia, K. Fabricius, G. C. Rossi, and G. Veneiano, Phs. Lett. B8, 33 (98). [5] P. van Baal, Commun. Math. Phs. 85, 59 (98). [] J. Greensite, Prog. Part. Nucl. Phs. 5, (3), hep-lat/33. [7] H. Reinhardt, O. Schroeder, T. Tok, and V. C. Zhukovsk, Phs. Rev. D, 854 (), hep-th/3. [8] R. Bertle, M. Faber, J. Greensite, and S. Olejnik, JHEP 3, 9 (999), hep-lat/9933. [9] In this work, we denote as ero-mode an eigenmode of D D of the overlap Dirac operator () with eigenvalue smaller than 5. If such modes eist of both positive and negative chiralit the are most likel not eact ero-modes, but lifted slightl above ero pairwise. This numerical inaccurac does not affect the value of the inde (). 9