Eigenvalue distributions of the QCD Dirac operator

We compute by Monte Carlo methods the individual distributions of the kth smallest Dirac operator eigenvalues in QCD, and compare them with recent analytical predictions. We do this for both massless and massive quarks in an SU(3) gauge theory with staggered fermions. Very precise agreement is found in all cases. As a simple by-product we also extract the microscopic spectral density of the Dirac operator in SU(3) gauge theory with dynamical massive fermions for Nf=1 and 2, and obtain high-accuracy agreement with analytical expressions.