Spectrum of the SU(3) Dirac operator on the lattice: Transition from random matrix theory to chiral perturbation theory

We calculate complete spectra of the Kogut-Susskind Dirac operator on the lattice in quenched SU(3) gauge theory for various values of coupling constant and lattice size. From these spectra we compute the connected and disconnected scalar susceptibilities and find agreement with chiral random matrix theory up to a certain energy scale, the Thouless energy. The dependence of this scale on the lattice volume is analyzed. In the case of the connected susceptibility this dependence is anomalous, and we explain the reason for this. We present a model of chiral perturbation theory that is capable of describing the data beyond the Thouless energy and that has a common range of applicability with chiral random matrix theory.Comment: 8 pages, RevTeX, 15 .eps figure

Equivalent of a Thouless energy in lattice QCD Dirac spectra

Random matrix theory (RMT) is a powerful statistical tool to model spectral fluctuations. In addition, RMT provides efficient means to separate different scales in spectra. Recently RMT has found application in quantum chromodynamics (QCD). In mesoscopic physics, the Thouless energy sets the universal scale for which RMT applies. We try to identify the equivalent of a Thouless energy in complete spectra of the QCD Dirac operator with staggered fermions and $SU_c(2)$ lattice gauge fields. Comparing lattice data with RMT predictions we find deviations which allow us to give an estimate for this scale.Comment: LATTICE99 (theor. devel.), 3 pages, 4 figure

Can we do better than Hybrid Monte Carlo in Lattice QCD?

The Hybrid Monte Carlo algorithm for the simulation of QCD with dynamical staggered fermions is compared with Kramers equation algorithm. We find substantially different autocorrelation times for local and nonlocal observables. The calculations have been performed on the parallel computer CRAY T3D.Comment: Talk presented at LATTICE96(algorithms), LaTeX 3 pages, uses espcrc2, epsf, 2 postscript figure

Small eigenvalues of the SU(3) Dirac operator on the lattice and in Random Matrix Theory

We have calculated complete spectra of the staggered Dirac operator on the lattice in quenched SU(3) gauge theory for \beta = 5.4 and various lattice sizes. The microscopic spectral density, the distribution of the smallest eigenvalue, and the two-point spectral correlation function are analyzed. We find the expected agreement of the lattice data with universal predictions of the chiral unitary ensemble of random matrix theory up to a certain energy scale, the Thouless energy. The deviations from the universal predictions are determined using the disconnected scalar susceptibility. We find that the Thouless energy scales with the lattice size as expected from theoretical arguments making use of the Gell-Mann--Oakes--Renner relation.Comment: REVTeX, 5 pages, 4 figure

Universal and non-universal behavior in Dirac spectra

We have computed ensembles of complete spectra of the staggered Dirac operator using four-dimensional SU(2) gauge fields, both in the quenched approximation and with dynamical fermions. To identify universal features in the Dirac spectrum, we compare the lattice data with predictions from chiral random matrix theory for the distribution of the low-lying eigenvalues. Good agreement is found up to some limiting energy, the so-called Thouless energy, above which random matrix theory no longer applies. We determine the dependence of the Thouless energy on the simulation parameters using the scalar susceptibility and the number variance.Comment: LATTICE98(confine), 9 pages, 11 figure

Statistical analysis and the equivalent of a Thouless energy in lattice QCD Dirac spectra

Random Matrix Theory (RMT) is a powerful statistical tool to model spectral fluctuations. This approach has also found fruitful application in Quantum Chromodynamics (QCD). Importantly, RMT provides very efficient means to separate different scales in the spectral fluctuations. We try to identify the equivalent of a Thouless energy in complete spectra of the QCD Dirac operator for staggered fermions from SU(2) lattice gauge theory for different lattice size and gauge couplings. In disordered systems, the Thouless energy sets the universal scale for which RMT applies. This relates to recent theoretical studies which suggest a strong analogy between QCD and disordered systems. The wealth of data allows us to analyze several statistical measures in the bulk of the spectrum with high quality. We find deviations which allows us to give an estimate for this universal scale. Other deviations than these are seen whose possible origin is discussed. Moreover, we work out higher order correlators as well, in particular three--point correlation functions.Comment: 24 pages, 24 figures, all included except one figure, missing eps file available at http://pluto.mpi-hd.mpg.de/~wilke/diff3.eps.gz, revised version, to appear in PRD, minor modifications and corrected typos, Fig.4 revise

Small eigenvalues of the staggered Dirac operator in the adjoint representation and Random Matrix Theory

The low-lying spectrum of the Dirac operator is predicted to be universal, within three classes, depending on symmetry properties specified according to random matrix theory. The three universal classes are the orthogonal, unitary and symplectic ensemble. Lattice gauge theory with staggered fermions has verified two of the cases so far, unitary and symplectic, with staggered fermions in the fundamental representation of SU(3) and SU(2). We verify the missing case here, namely orthogonal, with staggered fermions in the adjoint representation of SU(N_c), N_c=2, 3.Comment: 3 pages, revtex, 2 postscript figure

Fake symmetry transitions in lattice Dirac spectra

In a recent lattice investigation of Ginsparg-Wilson-type Dirac operators in the Schwinger model, it was found that the symmetry class of the random matrix theory describing the small Dirac eigenvalues appeared to change from the unitary to the symplectic case as a function of lattice size and coupling constant. We present a natural explanation for this observation in the framework of a random matrix model, showing that the apparent change is caused by the onset of chiral symmetry restoration in a finite volume. A transition from unitary to symplectic symmetry does not occur.Comment: 6 pages, 3 figures, REVTe

Lectures on Chiral Disorder in QCD

I explain the concept that light quarks diffuse in the QCD vacuum following the spontaneous breakdown of chiral symmetry. I exploit the striking analogy to disordered electrons in metals, identifying, among others, the universal regime described by random matrix theory, diffusive regime described by chiral perturbation theory and the crossover between these two domains.Comment: Lectures given at the Cargese Summer School, August 6-18, 200

Level Spacing Distribution of Critical Random Matrix Ensembles

We consider unitary invariant random matrix ensembles which obey spectral statistics different from the Wigner-Dyson, including unitary ensembles with slowly (~(log x)^2) growing potentials and the finite-temperature fermi gas model. If the deformation parameters in these matrix ensembles are small, the asymptotically translational-invariant region in the spectral bulk is universally governed by a one-parameter generalization of the sine kernel. We provide an analytic expression for the distribution of the eigenvalue spacings of this universal asymptotic kernel, which is a hybrid of the Wigner-Dyson and the Poisson distributions, by determining the Fredholm determinant of the universal kernel in terms of a Painleve VI transcendental function.Comment: 5 pages, 1 figure, REVTeX; restriction on the parameter stressed, figure replaced, refs added (v2); typos (factors of pi) in (35), (36) corrected (v3); minor changes incl. title, version to appear in Phys.Rev.E (v4