Distributions of individual Dirac eigenvalues
for QCD at
nonzero chemical potential:
RMT predictions and lattice results
Abstract
For QCD at nonzero chemical potential , the Dirac eigenvalues are scattered in the complex plane. We define a notion of ordering for individual eigenvalues in this case and derive the distributions of individual eigenvalues from random matrix theory (RMT). We distinguish two cases depending on the parameter , where is the volume and is the familiar lowenergy constant of chiral perturbation theory. For small , we use a Fredholm determinant expansion and observe that already the first few terms give an excellent approximation. For large , all spectral correlations are rotationally invariant, and exact results can be derived. We compare the RMT predictions to lattice data and in both cases find excellent agreement in the topological sectors .
Distributions of individual Dirac eigenvalues
for QCD at nonzero chemical potential:
RMT predictions and lattice results
Gernot Akemann, Jacques Bloch, Leonid Shifrin^{†}^{†}thanks: Speaker. , and Tilo Wettig
Department of Mathematical Sciences & BURSt Research Centre,
Brunel University West London, Uxbridge UB8 3PH, United Kingdom
Institute for Theoretical Physics, University of Regensburg, 93040 Regensburg, Germany
Email: Leonid.Shifrin@brunel.ac.uk
\abstract@cs
1 Introduction
Studies of the properties of the Dirac operator spectrum in gauge theories, including QCD, have a long history. For example, the lowlying Dirac modes provide information about spontaneous chiral symmetry breaking through the BanksCasher relation. The Dirac operator spectrum is also a natural object to study in lattice QCD. In the deep infrared, QCD in the regime can be described by the chiral random matrix theory (RMT) introduced in Ref. [1]. One of the advantages of RMT is that many exact analytical results can be derived. These results contain the lowenergy constant (LEC) of chiral perturbation theory (chPT), and in some cases also the LEC . These LECs can then be determined by fitting lattice data to the RMT curves.
The observables that are most natural to compute in RMT are the spectral density correlation functions. At zero chemical potential , all of them are known in RMT [1, 2]. On the other hand, one can consider individual Dirac eigenvalue distributions (IED). From the lattice QCD point of view, these are the most natural observables to measure directly. Certain quantities such as the average positions of the eigenvalues are more pronounced in IEDs and therefore require less statistics to be measured reliably.
At , the Dirac operator is antiHermitian and has a purely imaginary spectrum. In this case, all IEDs are known analytically in RMT [3]. They have become a standard tool in lattice QCD to extract in sectors of fixed topology. At , the Dirac operator is no longer antiHermitian, and its eigenvalues are scattered in the complex plane. Our work is based on an RMT for [4] which has an eigenvalue representation and for which all complex density correlations (both quenched and unquenched) have been calculated. The same results can be obtained from the RMT for introduced earlier by Stephanov [5] or from chPT in the regime [6, 7] and are universal in that sense. These results have been compared to data from quenched lattice simulations with staggered [8] and overlap [9] fermions at . A virtue of is that it couples to in leading order of chPT [10] so that a comparison with lattice data allows us to extract [8].
For the IEDs at much less is known. One of the problems here is to define an ordering of complex eigenvalues. Previous work on the repulsion between complex eigenvalues in RMT [11] and on the lattice [12] was done in the bulk of the spectrum, where no link to chPT is apparent. In the present work, we are interested in IEDs for eigenvalues close to the origin since they provide information on topological properties and LECs.
We first define the general notion of IEDs for complex eigenvalues, and then compute the first few IEDs approximately by truncating a socalled Fredholm determinant expansion to the first few terms. It was already observed in Ref. [13] for that this is a very good approximation. For large values of the parameter , where is the volume, we are able to derive all IEDs in closed form [14]. We use these exact results as a consistency check of the Fredholm determinant expansion. Our results are then compared to the lattice data of Ref. [9], in which the overlap Dirac operator for was constructed. This operator has good chiral properties (it satisfies a GinspargWilson relation, has an exact lattice chiral symmetry and exact zero modes, and satisfies the index theorem), which is essential for the present work. The same operator is obtained if a chemical potential is introduced in the domainwall fermion formalism in the limit of infinite extent of the fifth dimension [15]. Due to the sign problem at , the lattice analysis in Ref. [9] was restricted to the quenched case, and this restriction on the lattice data applies to this work as well.
2 Individual eigenvalue distributions for complex eigenvalues
Consider an operator with a finite number of complex eigenvalues, distributed according to a joint probability distribution which is symmetric in all its arguments. We also assume a symmetry and only consider the upper halfplane . The partition function is then given by , and the spectral density correlation functions are defined as
(2.0) 
Now consider any oneparameter family of mutually nonintersecting closed contours which cover . For fixed , is the boundary of a set . Let us parametrize the contour as with . We then define the th eigenvalue distribution as the probability that eigenvalues are inside , one is at the point on the boundary , and the remaining are in the complement ,
(2.0) 
Note that the eigenvalue ordering is induced by the entire contour family. The play the role of generalized polar coordinates. In the following, the argument of will be suppressed. It is possible [14] to express all through the densities Eq. (2). In particular, for the distribution of the first eigenvalue one obtains
(2.0) 
One can show that the integrated distributions are normalized as , where is the range of , and is the Jacobian of the transformation from to . The proof is similar to that for real eigenvalues [13]. We emphasize that in this framework the choice of the contour family becomes part of the definition of the quantities we measure (i.e., the ’s). Different contour families in general lead to different ’s. However, one relation always holds trivially, namely .
3 RMT predictions
The partition function of the matrix model we use [4] reads
(3.0) 
Here, and are complex matrices with no further symmetries, is the topological charge, the are the masses of flavors of dynamical quarks, and is the chemical potential in the matrix model. In the large limit this model describes QCD in the regime. All density correlation functions of this model follow from the kernel of biorthogonal polynomials with respect to the weight
(3.0) 
where (and below) are modified Bessel functions, according to
(3.0) 
We rescale the parameters of the model such that the parameters , , and stay finite in the large limit. The scaling of these parameters in terms of the LECs of chPT is given in parentheses. In the quenched case, the microscopic kernel is given by [4]
(3.0) 
and the microscopic spectral density follows as .
Approximate computations for arbitrary
In order to obtain the distribution , we use Eq. (2) and substitute the densities from Eqs. (3) and (3). Similar formulas can be written down for other ’s. For practical purposes, we truncate the (socalled Fredholm determinant) expansion in Eq. (2) to the first three terms and the corresponding expansion for to the first two terms. The density and the distributions and of the first and second eigenvalue are shown in Fig. 1, together with their counterparts in the case. For exact results are available [3] which facilitate a detailed comparison. As for [13], we see that the expansion converges rapidly. Higherorder terms merely assure that remains zero for large .
Exact results in the large limit
In the limit, the problem becomes radially symmetric. For finite but large , the symmetry is apparent close to the origin. The limiting microscopic spectral density expressed in the new variable reads
(3.0) 
In this limit, we can derive a closed expression for all eigenvalue distributions [14]. Because of the rotational symmetry we choose to be a semicircle in of radius and obtain for
(3.0) 
Here, we have introduced the incomplete Bessel function for , and zero otherwise. Our result is analogous to the result for the nearestneighbor spacing distribution of the Ginibre ensemble [11, 12], which can be interpreted as the distribution of the smallest nonzero eigenvalue if one eigenvalue is fixed at zero. Expressions for with are also available [14, 16]. Fig. 2 shows several , summing up nicely to the spectral density Eq. (3).
4 Lattice calculations
The lattice part of our work is based on the data obtained in Ref. [9]. The overlap Dirac operator introduced there is
(4.0) 
where is the sign function of a nonHermitian matrix and is the Wilson Dirac operator at . This overlap operator was shown to satisfy a GinspargWilson relation and to have good chiral properties [9, 15]. Equation (4) reduces to the standard overlap operator [17] at .
From the computational standpoint, the most demanding part is the computation of the matrix sign function. For the present set of data, this was done exactly using the spectral definition of the sign function. The lattice size is only , since high statistics are needed for a comparison with RMT. The coupling in the Wilson action is in order to stay in the regime (where RMT applies) for the first eigenvalue(s) [9]. The Wilson mass is ( is the lattice spacing), and the quark mass is zero. Data were sampled in the topological sectors , 1, 2 for the values of , 0.2, 0.3, and 1.0, corresponding to , 0.615, 1.42, and 4.51, respectively. The number of configurations varied from about 9000 for to about for . The parameters and were determined by a fit to the spectral density from Eq. (3). For , the data showed rotational invariance up to , and hence only the combination could be determined by a fit to Eq. (3). Our comparisons use these values and are thus parameterfree.
To compute from the lattice, we choose for our contours concentric semicircles with radius for all values of . (Other choices are also possible.) The localized nature of the IEDs allows us to integrate over the phase (i.e., we compute ) rather than to consider cuts as in Ref. [9]. This procedure results in a much better signal. As a consequence, we are able to obtain a better comparison in topological sectors , 1 and, for the first time, to successfully test the RMT predictions in the sector. In Fig. 3 we show the comparison of RMT predictions and lattice results for , 0.3, 1.0 and , 1, 2 (the case is similar [16]). We see that for smaller , the agreement is excellent, whereas there are deviations for and , 2. In these cases we are outside the regime of QCD so that RMT no longer applies. We emphasize that while the ascent of the distributions from zero was in principle already tested in Ref. [9] through the density, their descent represents a new, parameterfree test.
The deviations of the theoretical curves from zero for large are an artifact of the truncation of the Fredholm determinant expansion. As (or ) is increased, the convergence of the approximation becomes slower, i.e., more terms are needed.
5 Conclusions
In this work we have studied individual Dirac eigenvalue distributions in the regime of QCD at nonzero chemical potential. We provided a general framework for ordering complex eigenvalues. Our RMT computation for arbitrary was based on the truncation of a Fredholm determinant expansion. In the limit, we were able to derive all IEDs analytically in closed form. These predictions were then tested against lattice data based on the generalization of the overlap Dirac operator to nonzero chemical potential, in the topological sectors , 1, 2. We found excellent agreement between RMT and lattice results for several values of in the domain of the applicability of RMT. The descent of the IEDs represents a new, parameterfree test of RMT predictions. The much improved signal (resulting from the integration over the phase) allowed us, for the first time, to successfully test the RMT predictions in the topological sector .
Acknowledgments
This work was supported by EPSRC grant EP/D031613/1 (GA & LS), by EU network ENRAGE MRTNCT2004005616 (GA) and by DFG grant FOR 465 (JB & TW).
References

[1]
E. V. Shuryak and J. J. M. Verbaarschot,
Nucl. Phys. A 560 (1993) 306
[hepth/9212088];
J. J. M. Verbaarschot, Phys. Rev. Lett. 72 (1994) 2531 [hepth/9401059]. 
[2]
G. Akemann, P. H. Damgaard, U. Magnea and S. Nishigaki,
Nucl. Phys. B 487 (1997) 721
[hepth/9609174];
P. H. Damgaard and S. M. Nishigaki, Nucl. Phys. B 518 (1998) 495 [hepth/9711023]. 
[3]
T. Wilke, T. Guhr and T. Wettig,
Phys. Rev. D 57 (1998) 6486
[hepth/9711057];
S. M. Nishigaki, P. H. Damgaard and T. Wettig, Phys. Rev. D 58 (1998) 087704 [hepth/9803007];
P. H. Damgaard and S. M. Nishigaki, Phys. Rev. D 63 (2001) 045012 [hepth/0006111]. 
[4]
J. C. Osborn,
Phys. Rev. Lett. 93 (2004) 222001
[hepth/0403131];
G. Akemann, J. C. Osborn, K. Splittorff and J. J. M. Verbaarschot, Nucl. Phys. B 712 (2005) 287 [hepth/0411030].  [5] M. A. Stephanov, Phys. Rev. Lett. 76 (1996) 4472 [heplat/9604003].
 [6] K. Splittorff and J. J. M. Verbaarschot, Nucl. Phys. B 683 (2004) 467 [hepth/0310271].
 [7] F. Basile and G. Akemann, arXiv:0710.0376 [hepth].

[8]
G. Akemann and T. Wettig,
Phys. Rev. Lett. 92 (2004) 102002
[Erratumibid. 96 (2006) 029902]
[heplat/0308003];
J. C. Osborn and T. Wettig, PoS LAT2005 (2006) 200 [heplat/0510115].  [9] J. Bloch and T. Wettig, Phys. Rev. Lett. 97 (2006) 012003 [heplat/0604020].
 [10] D. Toublan and J. J. M. Verbaarschot, Int. J. Mod. Phys. B 15 (2001) 1404 [hepth/0001110].
 [11] R. Grobe, F. Haake and H.J. Sommers, Phys. Rev. Lett. 61 (1988) 1899.
 [12] H. Markum, R. Pullirsch and T. Wettig, Phys. Rev. Lett. 83 (1999) 484 [heplat/9906020].
 [13] G. Akemann and P. H. Damgaard, Phys. Lett. B 583 (2004) 199 [hepth/0311171].
 [14] G. Akemann and L. Shifrin (2007), to be published.
 [15] J. Bloch and T. Wettig, arXiv:0709.4630 [heplat], to appear in Phys. Rev. D.
 [16] G. Akemann, J. C. R. Bloch, L. Shifrin and T. Wettig, arXiv:0710.2865 [heplat].

[17]
R. Narayanan and H. Neuberger,
Nucl. Phys. B 443 (1995) 305
[hepth/9411108];
H. Neuberger, Phys. Lett. B 417 (1998) 141 [heplat/9707022].