Fermion determinants in matrix models of QCD at nonzero chemical potential

Abstract

The presence of a chemical potential completely changes the analytical structure of the QCD partition function. In particular, the eigenvalues of the Dirac operator are distributed over a finite area in the complex plane, whereas the zeros of the partition function in the complex mass plane remain on a curve. In this paper we study the effects of the fermion determinant at nonzero chemical potential on the Dirac spectrum by means of the resolvent, G(z), of the QCD Dirac operator. The resolvent is studied both in a one-dimensional U(1) model (Gibbs model) and in a random matrix model with the global symmetries of the QCD partition function. In both cases we find that, if the argument z of the resolvent is not equal to the mass m in the fermion determinant, the resolvent diverges in the thermodynamic limit. However, for z =m the resolvent in both models is well defined. In particular, the nature of the limit $z \rightarrow m$ is illuminated in the Gibbs model. The phase structure of the random matrix model in the complex m and μ-planes is investigated both by a saddle point approximation and via the distribution of Yang-Lee zeros. Both methods are in complete agreement and lead to a well-defined chiral condensate and quark number density.