---
title: Wilson–Dirac Spectrum in Lattice QCD
url: https://www.emergentmind.com/topics/wilson-dirac-spectrum
type: topic
---

# Wilson–Dirac Spectrum in Lattice QCD

The Wilson–Dirac spectrum refers to the statistical properties of the eigenvalues of the Wilson–Dirac operator and its hermitized forms in lattice gauge theory. Detailed analysis of this spectrum provides direct access to nonperturbative aspects of chiral symmetry breaking, the impact of lattice discretization, and critical low-energy constants (LECs) within the Wilson chiral perturbation theory (WχPT) framework. It is foundational to precision extrapolation of lattice QCD results to the continuum and understanding lattice-induced chiral symmetry breaking.

## 1. The Wilson–Dirac Operator and Hermitian Forms

The Wilson–Dirac operator $D_W$ is the standard lattice discretization of the Dirac operator, augmented by a Wilson term to remove doubler modes at finite lattice spacing $a$. In $d=4$, and after inclusion of a bare mass $m_0$:
\[
D_W = m_0 + \frac12 \sum_\mu \left[ \gamma_\mu (\nabla_\mu + \nabla^*_\mu) - a\,\nabla^*_\mu\nabla_\mu \right]
\]
where $\gamma_\mu$ are Euclidean Dirac matrices, and $\nabla_\mu$, $\nabla^*_\mu$ are forward/backward covariant derivatives. 

$D_W$ obeys $\gamma_5$-Hermiticity, $D_W^\dagger = \gamma_5 D_W \gamma_5$, but is nonnormal for $a>0$—hence not all eigenvalues are purely imaginary. For spectral analysis, the Hermitian Wilson–Dirac operator is introduced:
\[
D_5 \equiv \gamma_5 (D_W + m)
\]
which has a real spectrum, making microscopic spectral analysis and comparison with field theory tractable [1012.0752].

In twisted-mass lattice QCD, the Hermitian operator is generalized to
\[
D_5 \equiv \gamma_5 (D_W + m + i\mu \gamma_5 \tau_3)
\]
where $m$ is the untwisted mass and $\mu$ the twisted mass [1510.09169].

## 2. Microscopic Spectral Density and Wilson χPT

The microscopic spectral density refers to the rescaled spectral density in the so-called $\epsilon$-regime ($V\to\infty$ with $m\Sigma V$, $a^2W_iV$, $\lambda^5\Sigma V$ fixed). For index $\nu$:
\[
\rho_s^\nu(\hat\lambda; \hat{m}, \hat{a}_i) = \lim_{V \to \infty} \frac{1}{\Sigma V} \left\langle \sum_k \delta(\hat\lambda - \lambda_k^5 \Sigma V) \right\rangle_\nu
\]
where $\Sigma$ is the infinite-volume chiral condensate, and the $\hat{a}_i^2 = a^2 W_i V$ describe cutoff effects with the LECs $W_6$, $W_7$, $W_8$ entering at $O(a^2)$ in the chiral Lagrangian [1012.0752, 1001.2937, 1112.5160].

At leading order, the zero-momentum sector of WχPT yields a group integral for the partition function at fixed $\nu$:
\[
Z_{N_f}^\nu(m,z;a) = \int_{U(N_f)} dU \det(U)^\nu \exp \left[ S[U] \right]
\]
with
\[
S[U] = \frac12 (m+z)\Sigma \operatorname{Tr} U + \frac12 (m-z)\Sigma \operatorname{Tr} U^\dagger 
- a^2 W_6 [\operatorname{Tr}(U+U^\dagger)]^2
- a^2 W_7 [\operatorname{Tr}(U-U^\dagger)]^2
- a^2 W_8 \operatorname{Tr}(U^2+U^{\dagger2})
\]
Extracting the microscopic spectral density involves constructing a (graded) generating functional and taking imaginary parts of the partially quenched resolvent or its supersymmetric extension [1012.0752, 1001.2937].

## 3. Analytical Results: Diffusion Representation and Random Matrix Theory Matching

For $N_f=2$ (twisted-mass), the microscopic spectral density in sector $\nu$ is given by:
\[
\rho_S^\nu(\zeta; m̂, â) = \frac{1}{\pi} \Im G_{3|1}^\nu(-\zeta, m̂, â)
\]
where $G_{3|1}^\nu$ is derived from the graded chiral partition function [1510.09169, 1201.1361].

A key result is the diffusion representation [1112.0377]:
\[
\rho_s^\nu(\xi; \tau) = \int_0^\infty d\eta\, K_\nu(\xi, \eta; \tau)\,\rho_{s,0}^\nu(\eta)
\]
with $K_\nu$ a Bessel–diffusion kernel reflecting the smearing of zero modes due to $a^2 W_8$ effects.

Wilson chiral Random Matrix Theory (WRMT) provides an alternative, mathematically equivalent formulation for the microscopic spectrum via ensembles of random matrices matching the symmetry class and index structure of the Wilson–Dirac operator [1012.0752, 1307.7251, 1001.2937, 1112.5160].

The mapping between chiral Lagrangian and WRMT parameters is:
\[
N m = \Sigma V m,\quad N z = \Sigma V z,\quad \frac{N \tilde{a}^2}{4} = a^2 W_8 V
\]
This mapping ensures that microscopic observables—eigenvalue densities, spacing distributions, and individual eigenvalue distributions—agree precisely between WχPT and WRMT in the microscopic scaling limit [1012.0752, 1001.2937].

## 4. Extraction of Low-Energy Constants and Spectral Observables

Microscopic spectral properties, particularly near the origin, are sensitive to the chiral condensate $\Sigma$ and the Wilson LECs ($W_6$, $W_7$, $W_8$). Several practical methods have been established:
- The width of the near-zero eigenvalue peak (smeared former zero modes) is $\sim a\sqrt{W_8 V}$, directly fixing $W_8$ [1112.0377, 1510.09169].
- The spectral density of the Hermitian operator $D_5$ and distribution of chirality over real $D_W$ eigenvalues allow for fits of $\Sigma$, $W_8$, and (where resolved) $W_6$, $W_7$ [1112.5160, 1307.7251].
- Volume- and mass-scaling tests of spectral data confirm that fitting the theoretical microscopic density to lattice histograms (usually in fixed-$\nu$ sectors) yields robust extractions of LECs. The mapping of sector scaling for parameters $$\hat{m} = m \Sigma V,~\hat{a}_i = a\sqrt{W_i V}$$ is confirmed by data [1110.2851, 1110.4002, 1510.09169].

Typical values extracted for $\Sigma$ and $W_8$ in $N_f=2$ twisted-mass QCD, using the ETM lattice data, are $\Sigma^{1/3} = 271.1(7.3)_{\rm stat}~{\rm MeV}$ and $W_8 = 0.0064(12)_{\rm stat}$ (in suitably normalized units) [1510.09169]. The analytic determination of all $O(a^2)$ LECs is possible in principle by fitting to sufficiently rich, unquenched spectra.

## 5. Topology, Twisted Mass, and Physical Regime Dependence

The Wilson–Dirac spectrum exhibits strong dependence on the topological index $\nu$. For $|\nu|>0$, the continuum limit shows exact zero modes—at finite $a$, these are broadened into Gaussian peaks whose properties are governed by $W_8$. The peak structure, repulsion from the origin, and overall normalization of the microscopic density encode information about topology and discretization effects [1510.09169, 1112.0377].

Large twisted mass $\mû$ suppresses dynamical fermion contributions, tending the spectrum towards the quenched (determinantless) form. For smaller $\mû$ the dynamical effects are pronounced, influencing the density near the origin [1510.09169, 1201.1361].

The combined scaling regime ($a^2 V \sim O(1)$, $m\Sigma V \sim O(1)$) ensures that the lowest eigenvalues probed are below the Thouless energy and thus dominated by chiral zero-mode dynamics, with finite-volume corrections suppressed [1510.09169, 1110.4002].

## 6. Spectral Features, Continuum Limit, and Lattice Artifacts

Nonzero lattice spacing causes the “hard-edge” structure of the continuum spectrum to be smoothed. For small $a$, $O(a^2)$ corrections parametrized by $W_8$ dominate, while $W_6$, $W_7$ encode double-trace lattice artifacts and are suppressed for $N_c\to\infty$.

Key features:
- The zero-mode peak for $|\nu| > 0$ becomes a finite-width Gaussian for nonzero $a$.
- Tail states appear within the continuum gap as Lifshitz tails, scaling as $\exp[-(x-m)^2/\mathcal{O}(a^2)]$.
- At larger $a^2 V$, the spectrum exhibits square-root edges analogous to Tracy–Widom distributions.
- Dynamical quarks further suppress the singularity at the origin, removing the $1/\sqrt{\xi}$ divergence present in the quenched density [1112.0377].

Fits to the spectrum of $D_5$ (including the lowest several eigenvalues across topological sectors) can thus separate continuum physics from lattice artifacts and facilitate continuum extrapolation. The methodology—analytical prediction from WχPT/WRMT, followed by sector-resolved histogram fitting—is standard [1510.09169, 1110.2851, 1112.5160].

## 7. Impact, Applications, and Extensions

The Wilson–Dirac spectrum is central to several directions in lattice gauge theory:
- It enables precision extraction of the chiral condensate and $O(a^2)$ LECs directly from spectral data.
- Control of lattice artifacts through analysis of the microscopic spectrum informs improved action development, tuning of simulation parameters (e.g., fixing $m$, $a$ to avoid the Aoki phase), and direct assessment of universality [1301.3099].
- The approach extends to enriched lattice actions (e.g., twisted-mass, clover improvement), different gauge groups, and probes of topological charge distributions.
- The analytic realization and matching with WRMT provide efficient computational tools for analysis and interpolation, now standard in extracting nonperturbative low-energy constants [1510.09169, 1001.2937, 1112.0377, 1307.7251].

In summary, the detailed understanding of the Wilson–Dirac spectrum at the microscopic level—grounded in WχPT and confirmed via lattice data—forms the backbone of modern spectral analyses in lattice QCD, chiral dynamics, and the study of discretization effects in strongly coupled gauge theories.

Source: https://www.emergentmind.com/topics/wilson-dirac-spectrum