  Random matrix theory is a powerful way to describe universal correlations of
eigenvalues of complex systems. It also may serve as a schematic model for
disorder in quantum systems. In this review, we discuss both types of
applications of chiral random matrix theory to the QCD partition function. We
show that constraints imposed by chiral symmetry and its spontaneous breaking
determine the structure of low-energy effective partition functions for the
Dirac spectrum. We thus derive exact results for the low-lying eigenvalues of
the QCD Dirac operator. We argue that the statistical properties of these
eigenvalues are universal and can be described by a random matrix theory with
the global symmetries of the QCD partition function. The total number of such
eigenvalues increases with the square root of the Euclidean four-volume. The
spectral density for larger eigenvalues (but still well below a typical
hadronic mass scale) also follows from the same low-energy effective partition
function. The validity of the random matrix approach has been confirmed by many
lattice QCD simulations in a wide parameter range. Stimulated by the success of
the chiral random matrix theory in the description of universal properties of
the Dirac eigenvalues, the random matrix model is extended to nonzero
temperature and chemical potential. In this way we obtain qualitative results
for the QCD phase diagram and the spectrum of the QCD Dirac operator. We
discuss the nature of the quenched approximation and analyze quenched Dirac
spectra at nonzero baryon density in terms of an effective partition function.
Relations with other fields are also discussed.
