Dirac spectrum in QCD and quark masses

We use a chiral random matrix model to investigate the effects of massive quarks on the distribution of eigenvalues of QCD inspired Dirac operators. Kalkreuter's lattice analysis of the spectrum of the massive (hermitian) Dirac operator for two colors and Wilson fermions is shown to follow from a cubic equation in the quenched approximation. The quenched spectrum shows a Mott transition from a (delocalized) Goldstone phase softly broken by the current mass, to a (localized) heavy quark phase, with quarks localized over their Compton wavelength. Both phases are distinguishable by the quark density of states at zero virtuality, with a critical quark mass of the order of 100-200 MeV.At the critical point, the quark density of states is given by ν _Q (λ) ∼ |λ| ^{1}/{3}. Using Grassmaniian techniques, we derive an integral representation for the resolvent of the massive Dirac operator with one flavor in the unquenched approximation, and show that near zero virtuality the distribution of eigenvalues is quantitatively changed by a non-zero quark mass. The generalization of our construction to arbitrary flavors is also discussed. Some recommendations for lattice simulations are suggested.