Progress on the Microscopic Spectrum of the Dirac Operator for QCD with Wilson Fermions
Starting from the chiral Lagrangian for Wilson fermions at nonzero lattice spacing we have obtained compact expressions for all spectral correlation functions of the Hermitian Wilson Dirac operator in the -domain of QCD with dynamical quarks. We have also obtained the distribution of the chiralities over the real eigenvalues of the Wilson Dirac operator for any number of flavors. All results have been derived for a fixed index of the Dirac operator. An important effect of dynamical quarks is that they completely suppress the inverse square root singularity in the spectral density of the Hermitian Wilson Dirac operator. The analytical results are given in terms of an integral over a diffusion kernel for which the square of the lattice spacing plays the role of time. This approach greatly simplifies the expressions which we here reduce to the evaluation of two-dimensional integrals.
Progress on Wilson Dirac Spectrum
\FullConference The XXIX International Symposium on Lattice Field Theory - Lattice 2011
July 10-16, 2011
Squaw Valley, Lake Tahoe, California
1 Introduction and Motivation
The low-lying spectrum of the continuum Dirac operator has been described successfully in terms of chiral Lagrangians or Random Matrix Theory [1, 2]. The reason for this success is that the smallest eigenvalues are of order , which is deep inside the -domain. In this domain the static part of the chiral Lagrangian, which is equivalent to Random Matrix Theory, factorizes from the partition function, and determines properties of the low-lying Dirac eigenvalues. In this lecture we discuss the extension of these results to nonzero lattice spacing. We can distinguish two different types of discretization effects: those that do not affect the symmetries of the Dirac operator and those that do affect its symmetries. In the first case, discretization effects can be absorbed into a redefinition of the low-energy constants. In the second case, the discretization errors take us outside the universality class which requires additional terms in the chiral Lagrangian. The two main discretization schemes are staggered fermions and Wilson fermions and both break the symmetries of the continuum Dirac operator. In this talk we will discuss recent results [3, 4, 5, 6] for the the Wilson Dirac operator including its dependence on the topological index and the number of flavors. The confirmation of these results by lattice simulations [23, 8] will be discussed elsewhere in these Proceedings. For recent progress on the staggered Dirac operator see .
We will consider discretization effects of the Wilson Dirac operator
Since is nonhermitian, its eigenvalues are complex, but because of the Hermiticity, they occur in complex conjugate pairs or are real. Therefore the complex eigenvalues can collide with the real axis and turn into a pair of real eigenvalues. This is different from QCD at nonzero chemical potential  where the nonzero Dirac eigenvalues are complex and occur in pairs of opposite sign. Then both their real and imaginary parts have to vanish for a pair to collide with the real axis.
In this talk we report on results for discretization effects on spectra of the Wilson Dirac operator of QCD with quarks. A direct evaluation of the generating function is complicated, but using a diffusion method  we obtain compact expressions for correlation functions of Dirac eigenvalues which we here further simplified by exploiting an underlying Pfaffian structure .
2 Spectrum of the Dirac Operator
In order to access the spectrum of , we introduce an axial mass
In the microscopic domain, where the combinations , and are kept fixed in the thermodynamic limit, the , and -dependence of the partition function resides in its zero momentum part that factorizes from the nonzero momentum part 
To this order the chiral Lagrangian also contains but they can be added to the mass term after linearizing the squares at the expense of a Gaussian integral and will not be considered below.
Since small Dirac eigenvalues are very sensitive to the index, , of the Dirac operator we will work at fixed . This index is determined from the spectral flow of the eigenvalues of [22, 23] (see Fig. 1). At the crossing point, , with the real axis we have (for small , the physical part of the lattice Dirac spectrum can be separated from the unphysical part)
Note that multiple crossings may occur. Therefore the number of real eigenvalues of is at least equal to its index. The total number of flow lines with a net flow across the real axis is a topological invariant of the Dirac operator.
It can be written as a sum over eigenvalues of 
The -dependence of the QCD partition function is given by
In the continuum limit, because of the anomaly, the factor can be written as the Jacobian of an axial transformation on the quark fields. At nonzero the natural extension is (See .)
As is the case in the continuum limit, the partition function at fixed is defined by inverting the decomposition (8), which also can be applied to the chiral Lagrangian
where the -dependence of the chiral Lagrangian is induced by the transformation (9).
The generating functions for the Wilson Dirac spectrum are given by 
The resolvent and the spectral density are given by
For the -dependence of the generating function is given by a chiral Lagrangian that is uniquely determined by symmetries that are compatible with the convergence of the bosonic integrals. In the microscopic domain the generating function reduces to the graded integral [3, 2]
where, . The transformation is required to get convergent integrals in the noncompact sector for (for which the Aoki phase  can occur).
A direct evaluation of the integrals over the supergroup becomes very laborious for but they still can be worked out for one flavor [5, 26]. However, for degenerate quark masses the expressions can be simplified by starting from the identity
where is an graded “Hermitian” matrix. In a diagonal representation of the partition function can be expressed in terms of an integral over shifted quark masses. The result for the resolvent is given by 
where is the Jacobian of transforming to the diagonal representation .
2.1 Effect of Light Flavors
For the spectral density of is obtained by a simple variable transformation of the spectral density of : (See Fig. 2 (right)). For the result for the spectral density follows from the resolvent (LABEL:G_N_f+p|p) (see Fig. 2 (left)). Although, for and the integrals in (LABEL:G_N_f+p|p) can still be evaluated directly, it becomes computationally expensive. However, because of an underlying Pfaffian structure , the integrals can be written as the sum of products of two and one dimensional integrals.
The three dimensional integral for one flavor resolvent (LABEL:G_N_f+p|p) can be rewritten as
where the first term is the quenched resolvent with
The various partition functions that enter in the expression are given by
Likewise the four-dimensional integral (LABEL:G_N_f+p|p) for the resolvent be re-expressed as
These simplified expressions for the unquenched resolvents have not been presented previously.
2.2 Density of Real Eigenvalues
The Wilson Dirac operator has as least real eigenvalues in the physical sector of the Dirac spectrum. Because the eigenvalues of the Wilson Dirac operator occur in complex conjugate pairs, additional real eigenvalues may be produced when a pair collides with the real axis. Depending on the expectation value of the density of the real eigenvalues may be decomposed as 
The distribution of the chiralities over real part of the spectrum of
The expression for can again be simplified to at most two-dimensional integrals.
We have obtained the lattice spacing dependence of the microscopic spectral density and all spectral correlators of the Hermitian Wilson Dirac operator with degenerate dynamical quark masses. This makes it possible to extract the low-energy constant from the distribution of the smallest eigenvalues. It should satisfy the constraint due to -Hermiticity that (see also the recent discussion in ). These results can be extended to include the two other low-energy constants to this order. Our results also explain the lattice result that the width of the distribution of the smallest eigenvalue scales as . More results can be obtained using a random matrix formulation of the chiral Lagrangian. Then powerful random matrix methods make it possible to obtain all spectral correlation functions of in the microscopic domain.
Acknowledgments: This work was supported by U.S. DOE Grant No. DE-FG-88ER40388 (JV) and the Sapere Aude program of The Danish Council for Independent Research (KS). Gernot Akemann, Poul Damgaard, Thomas Guhr and Mario Kieburg are thanked for useful discussions.
- E. V. Shuryak and J. J. M. Verbaarschot, Nucl. Phys. A 560, 306 (1993) [hep-th/9212088]; J. J. M. Verbaarschot and I. Zahed, Phys. Rev. Lett. 70, 3852 (1993) [hep-th/9303012]; J. J. M. Verbaarschot, Phys. Rev. Lett. 72, 2531 (1994) [hep-th/9401059].
- P. H. Damgaard, J. C. Osborn, D. Toublan and J. J. M. Verbaarschot, Nucl. Phys. B 547, 305 (1999) [hep-th/9811212].
- P. H. Damgaard, K. Splittorff and J. J. M. Verbaarschot, Phys. Rev. Lett. 105, 162002 (2010). [arXiv:1001.2937 [hep-th]].
- G. Akemann, P. H. Damgaard, K. Splittorff, J. J. M. Verbaarschot, Phys. Rev. D 83, 085014 (2011) [arXiv:1012.0752 [hep-lat]].
- G. Akemann, P. H. Damgaard, K. Splittorff, J. J. M. Verbaarschot, PoS LATTICE2010, 079 (2010). [arXiv:1011.5121 [hep-lat]].
- K. Splittorff, J. J. M. Verbaarschot, Phys. Rev. D84, 065031 (2011). [arXiv:1105.6229 [hep-lat]].
- P. H. Damgaard, U. M. Heller, K. Splittorff, [arXiv:1110.2851 [hep-lat]].
- A. Deuzeman, U. Wenger, J. Wuilloud, [arXiv:1110.4002 [hep-lat]].
- J. C. Osborn, Phys. Rev. D83, 034505 (2011). [arXiv:1012.4837 [hep-lat]].
- L. Del Debbio, L. Giusti, M. Lüscher, R. Petronzio and N. Tantalo, JHEP 0602, 011 (2006) [hep-lat/0512021]; JHEP 0702, 056 (2007) [hep-lat/0610059].
- M. Golterman, S. R. Sharpe, R. L. Singleton, Jr., Phys. Rev. D71, 094503 (2005). [hep-lat/0501015].
- S. R. Sharpe, Phys. Rev. D 74, 014512 (2006) [arXiv:hep-lat/0606002].
- S. Necco, A. Shindler, [arXiv:1101.1778 [hep-lat]].
- G. Akemann, J. C. Osborn, K. Splittorff, J. J. M. Verbaarschot, Nucl. Phys. B712, 287-324 (2005). [hep-th/0411030].
- T. Guhr, Annals Phys. 250, 145-192 (1996).
- M. Kieburg, [arXiv:1109.5109 [math-ph]].
- M. Kieburg, J. J. M. Verbaarschot, S. Zafeiropoulos, [arXiv:1109.0656 [hep-lat]].
- G. Akemann, T. Nagao, JHEP 1110, 060 (2011). [arXiv:1108.3035 [math-ph]].
- S. R. Sharpe and R. L. Singleton, Phys. Rev. D 58, 074501 (1998) [hep-lat/9804028].
- G. Rupak and N. Shoresh, Phys. Rev. 66, 054503 (2002), [arXiv:hep-lat/0201019].
- O. Bär, G. Rupak and N. Shoresh, Phys. Rev. D 70, 034508 (2004), [arXiv:hep-lat/0306021].
- S. Itoh, Y. Iwasaki and T. Yoshie, Phys. Rev. D 36, 527 (1987).
- K. M. Bitar, U. M. Heller and R. Narayanan, Phys. Lett. B 418, 167 (1998). [arXiv:hep-th/9710052].
- S. R. Sharpe, Phys. Rev. D 79, 054503 (2009) [arXiv:0811.0409 [hep-lat]].
- S. Aoki, Phys. Rev. D 30 (1984) 2653.
- R. N. Larsen, arXiv:1110.5744 [hep-th].
- M. Kieburg, J. J. M. Verbaarschot, S. Zafeiropoulos, [arXiv:1110.2690 [hep-lat]].
- M. T. Hansen, S. R. Sharpe, [arXiv:1111.2404 [hep-lat]].