QCD Dirac operator at nonzero chemical potential: lattice data and matrix model

Recently, a non-Hermitian chiral random matrix model was proposed to describe the eigenvalues of the QCD Dirac operator at nonzero chemical potential. This matrix model can be constructed from QCD by mapping it to an equivalent matrix model which has the same symmetries as QCD with chemical potential. Its microscopic spectral correlations are conjectured to be identical to those of the QCD Dirac operator. We investigate this conjecture by comparing large ensembles of Dirac eigenvalues in quenched SU(3) lattice QCD at nonzero chemical potential to the analytical predictions of the matrix model. Excellent agreement is found in the two regimes of weak and strong non-Hermiticity, for several different lattice volumes

Individual complex Dirac eigenvalue distributions from random matrix theory and lattice QCD at nonzero chemical potential

We analyze how individual eigenvalues of the QCD Dirac operator at nonzero chemical potential are distributed in the complex plane. Exact and approximate analytical results for such distributions are derived from non-Hermitian random matrix theory. When comparing these to lattice QCD spectra close to the origin, excellent agreement is found for zero and nonzero topology at several values of the chemical potential. Our analytical results are also applicable to other physical systems in the same symmetry class

Distributions of individual Dirac eigenvalues for QCD at non-zero chemical potential: RMT predictions and lattice results

For QCD at non-zero chemical potential $\mu$, the Dirac eigenvalues are scattered in the complex plane. We define a notion of ordering for individual eigenvalues in this case and derive the distributions of individual eigenvalues from random matrix theory (RMT). We distinguish two cases depending on the parameter $\alpha=\mu^2 F^2 V$, where $V$ is the volume and $F$ is the familiar low-energy constant of chiral perturbation theory. For small $\alpha$, we use a Fredholm determinant expansion and observe that already the first few terms give an excellent approximation. For large $\alpha$, all spectral correlations are rotationally invariant, and exact results can be derived. We compare the RMT predictions to lattice data and in both cases find excellent agreement in the topological sectors $\nu=0,1,2$

Random matrix theory of unquenched two-colour QCD with nonzero chemical potential

We solve a random two-matrix model with two real asymmetric matrices whose primary purpose is to describe certain aspects of quantum chromodynamics with two colours and dynamical fermions at nonzero quark chemical potential mu. In this symmetry class the determinant of the Dirac operator is real but not necessarily positive. Despite this sign problem the unquenched matrix model remains completely solvable and provides detailed predictions for the Dirac operator spectrum in two different physical scenarios/limits: (i) the epsilon-regime of chiral perturbation theory at small mu, where mu^2 multiplied by the volume remains fixed in the infinite-volume limit and (ii) the high-density regime where a BCS gap is formed and mu is unscaled. We give explicit examples for the complex, real, and imaginary eigenvalue densities including Nf=2 non-degenerate flavours. Whilst the limit of two degenerate masses has no sign problem and can be tested with standard lattice techniques, we analyse the severity of the sign problem for non-degenerate masses as a function of the mass split and of mu. On the mathematical side our new results include an analytical formula for the spectral density of real Wishart eigenvalues in the limit (i) of weak non-Hermiticity, thus completing the previous solution of the corresponding quenched model of two real asymmetric Wishart matrices.Comment: 45 pages, 31 figures; references added, as published in JHE

Overlap Dirac operator at nonzero chemical potential and random matrix theory

We show how to introduce a quark chemical potential in the overlap Dirac operator. The resulting operator satisfies a Ginsparg-Wilson relation and has exact zero modes. It is no longer gamma_5-hermitian, but its nonreal eigenvalues still occur in pairs. We compute the spectral density of the operator on the lattice and show that, for small eigenvalues, the data agree with analytical predictions of nonhermitian chiral random matrix theory for both trivial and nontrivial topology.Comment: 4 pages, 2 figure

Dirac spectrum and chiral condensate for QCD at fixed Θ\Theta angle

Kieburg M, Verbaarschot JJM, Wettig T. Dirac spectrum and chiral condensate for QCD at fixed $\Theta$ angle. PHYSICAL REVIEW D. 2019;99(7): 074515.We analyze the mass dependence of the chiral condensate for QCD at nonzero theta angle and find that in general the discontinuity of the chiral condensate is not on the support of the Dirac spectrum. To understand this behavior we decompose the spectral density and the chiral condensate into contributions from the zero modes, the quenched part, and a remainder which is sensitive to the fermion determinant and is referred to as the dynamical part. We obtain general formulas for the contributions of the zero modes. Expressions for the quenched part, valid for an arbitrary number of flavors, and for the dynamical part, valid for one and two flavors, are derived in the microscopic domain of QCD. We find that at nonzero theta angle the quenched and dynamical parts of the Dirac spectral density are strongly oscillating with an amplitude that increases exponentially with the volume V and a period of order of 1/V. The quenched part of the chiral condensate becomes exponentially large at theta not equal 0, but this divergence is canceled by the contribution from the zero modes. The oscillatory behavior of the dynamical part of the density is essential for moving the discontinuity of the chiral condensate away from the support of the Dirac spectrum. As important by-products of this work we obtain analytical expressions for the microscopic spectral density of the Dirac operator at nonzero theta angle for both one- and two-flavor QCD with nonzero quark masses