A new Chiral Two-Matrix Theory for Dirac Spectra with Imaginary Chemical Potential

We solve a new chiral Random Two-Matrix Theory by means of biorthogonal polynomials for any matrix size $N$. By deriving the relevant kernels we find explicit formulas for all $(n,k)$-point spectral (mixed or unmixed) correlation functions. In the microscopic limit we find the corresponding scaling functions, and thus derive all spectral correlators in this limit as well. We extend these results to the ordinary (non-chiral) ensembles, and also there provide explicit solutions for any finite size $N$, and in the microscopic scaling limit. Our results give the general analytical expressions for the microscopic correlation functions of the Dirac operator eigenvalues in theories with imaginary baryon and isospin chemical potential, and can be used to extract the tree-level pion decay constant from lattice gauge theory configurations. We find exact agreement with previous computations based on the low-energy effective field theory in the two special cases where comparisons are possible.Comment: 31 pages 2 figures, v2 missing term in partially quenched results inserted, fig 2 update

A new chiral two-matrix theory for dirac spectra with imaginary chemical potential

We solve a new chiral Random Two-Matrix Theory by means of biorthogonal polynomials for any matrix size $N$. By deriving the relevant kernels we find explicit formulas for all $(n,k)$-point spectral (mixed or unmixed) correlation functions. In the microscopic limit we find the corresponding scaling functions, and thus derive all spectral correlators in this limit as well. We extend these results to the ordinary (non-chiral) ensembles, and also there provide explicit solutions for any finite size $N$, and in the microscopic scaling limit. Our results give the general analytical expressions for the microscopic correlation functions of the Dirac operator eigenvalues in theories with imaginary baryon and isospin chemical potential, and can be used to extract the tree-level pion decay constant from lattice gauge theory configurations. We find exact agreement with previous computations based on the low-energy effective field theory in the two special cases where comparisons are possible

Comparing matrix models and QCD lattice data with chemical potential

We present a quantitative analysis of the microscopic Dirac spectrum which is complex in the presence of a non-vanishing quark chemical potential. Data from quenched SU(3) lattice simulations for different volumes V and small values of the chemical potential are compared to analytical predictions from matrix models. We confirm the existence of two distinct limits for weakly and strongly nonhermitian Dirac operators. Good agreement is found in both limits, confirming the different scaling of chemical potential and eigenvalues with the volume

Microscopic universality of complex matrix model correlation functions at weak non-Hermiticity

The microscopic correlation functions of non-chiral random matrix models with complex eigenvalues are analyzed for a wide class of non-Gaussian measures. In the large-N limit of weak non-Hermiticity, where N is the size of the complex matrices, we can prove that all k-point correlation functions including an arbitrary number of Dirac mass terms are universal close to the origin. To this aim we establish the universality of the asymptotics of orthogonal polynomials in the complex plane. The universality of the correlation functions then follows from that of the kernel of orthogonal polynomials and a mapping of massive to massless correlators

On matrix model partition functions for QCD with chemical potential

Partition functions of two different matrix models for QCD with chemical potential are computed for an arbitrary number of quark and complex conjugate anti-quark flavors. In the large-N limit of weak nonhermiticity complete agreement is found between the two models. This supports the universality of such fermionic partition functions, that is of products of characteristic polynomials in the complex plane. In the strong nonhermiticity limit agreement is found for an equal number of quark and conjugate flavours. For a general flavor content the equality of partition functions holds only for small chemical potential. The chiral phase transition is analyzed for an arbitrary number of quarks, where the free energy presents a discontinuity of first order at a critical chemical potential. In the case of nondegenerate flavors there is first order phase transition for each separate mass scale

Macroscopic and microscopic (non-)universality of compact support random matrix theory

A random matrix model with a σ-model like constraint, the restricted trace ensemble (RTE), is solved in the large-n limit. In the macroscopic limit the smooth connected two-point resolvent G(z,w) is found to be non-universal, extending previous results from monomial to arbitrary polynomial potentials. Using loop equation techniques we give a closed though non-universal expression for G(z,w), which extends recursively to all higher k-point resolvents. These findings are in contrast to the usual unconstrained one-matrix model. However, in the microscopic large-n limit, which probes only correlations at distance of the mean level spacing, we are able to show that the constraint does not modify the universal sine-law. In the case of monomial potentials V(M)=M2p, we provide a relation valid for finite-n between the k-point correlation function of the RTE and the unconstrained model. In the microscopic large-n limit they coincide which proves the microscopic universality of RTEs

The complex laguerre symplectic ensemble of Non-Hermitian matrices

We solve the complex extension of the chiral Gaussian Symplectic Ensemble, defined as a Gaussian two-matrix model of chiral non-Hermitian quaternion real matrices. This leads to the appearance of Laguerre polynomials in the complex plane and we prove their orthogonality. Alternatively, a complex eigenvalue representation of this ensemble is given for general weight functions. All k-point correlation functions of complex eigenvalues are given in terms of the corresponding skew orthogonal polynomials in the complex plane for finite-N, where N is the matrix size or number of eigenvalues, respectively. We also allow for an arbitrary number of complex conjugate pairs of characteristic polynomials in the weight function, corresponding to massive quark flavours in applications to field theory. Explicit expressions are given in the large-N limit at both weak and strong non-Hermiticity for the weight of the Gaussian two-matrix model. This model can be mapped to the complex Dirac operator spectrum with non-vanishing chemical potential. It belongs to the symmetry class of either the adjoint representation or two colours in the fundamental representation using staggered lattice fermions

Microscopic and bulk spectra of Dirac operators from finite-volume partition functions

The microscopic spectrum of the QCD Dirac operator is shown to obey random matrix model statistics in the bulk region of the spectrum close to the origin using finite-volume partition functions

New critical matrix models and generalized universality

We study a class of one-matrix models with an action containing nonpolynomial terms. By tuning the coupling constants in the action to criticality we obtain that the eigenvalue density vanishes as an arbitrary real power at the origin, thus defining a new class of multicritical matrix models. The corresponding microscopic scaling law is given and possible applications to the chiral phase transition in QCD are discussed. For generic coupling constants off-criticality we prove that all microscopic correlation functions at the origin of the spectrum remain in the known Bessel universality class. An arbitrary number of Dirac mass terms can be included and the corresponding massive universality is maintained as well. We also investigate the critical behavior at the edge of the spectrum: there, in contrast to the behavior at the origin, we find the same critical exponents as derived from matrix models with a polynomial action

Correlations for non-Hermitian Dirac operators: chemical potential in three-dimensional QCD

In the presence of a non-vanishing chemical potential the eigenvalues of the Dirac opera- tor become complex. We use a Random Matrix Model (RMM) approach to calculate ana- lytically all correlation functions at weak and strong non-Hermiticity for three-dimensional QCD with broken flavor symmetry and four-dimensional QCD in the bulk