Fermion determinants in matrix models of QCD at nonzero chemical potential

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 a nonzero chemical potential on the Dirac spectrum by means of the resolvent [Formula presented] 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 [Formula presented] of the resolvent is not equal to the mass [Formula presented] in the fermion determinant, the resolvent diverges in the thermodynamic limit. However, for [Formula presented] the resolvent in both models is well defined. In particular, the nature of the limit [Formula presented] is illuminated in the Gibbs model. The phase structure of the random matrix model in the complex [Formula presented] 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.