[Paper Review] Dynamical overlap fermions: techniques and results

Stefan Schaefer, Thomas DeGrand|ArXiv.org|Aug 26, 2005

Quantum Chromodynamics and Particle Interactions18 references3 citations

TL;DR

This paper presents a stable, efficient algorithm for dynamical overlap fermion simulations in lattice QCD using stout-smeared gauge links and a multiple time-scale integrator, enabling accurate computation of the chiral condensate via Random Matrix Theory analysis of Dirac eigenmodes. The method achieves a tenfold speedup and successfully measures Σ ≈ (280 MeV)³ at light quark masses, validating the approach against RMT predictions.

ABSTRACT

We summarize our recent investigations of lattice QCD with dynamical overlap fermions. We sketch algorithmic issues and our approach to solving them. We show our measurement of the topological susceptibility. We describe a computation of the chiral condensate using an analysis of the distribution of eigenmodes of the Dirac operator and Random Matrix Theory.

Motivation & Objective

• To develop a stable and efficient algorithm for dynamical overlap fermion simulations in lattice QCD.
• To address the two major challenges: slow simulation speed and difficulty in changing topological charge.
• To measure the chiral condensate using the distribution of low-lying Dirac eigenmodes and Random Matrix Theory (RMT).
• To validate the method by comparing RMT predictions with simulated eigenmode distributions across different topological sectors.
• To explore the behavior of the system at light quark masses and confirm ergodicity in restricted topological sectors.

Proposed method

• Use of stout-smeared gauge links to improve the condition number of the kernel operator, reducing the number of small eigenmodes and speeding up the overlap fermion action evaluation.
• Adoption of a multiple time-scale molecular dynamics integrator (Sexton-Weingarten) with a 1/12 ratio of gauge to fermion time steps to handle the large difference in force scales.
• Implementation of a refraction/reflection procedure to handle discontinuities in the effective action when eigenvalues of the kernel operator cross zero, preserving detailed balance in the HMC algorithm.
• Application of the Sherman–Morrison formula to compute the step change in the effective action and fermionic determinant at topology-changing surfaces.
• Use of the spectral representation of the sign function in the overlap Dirac operator to define the Hermitian overlap operator and its squared form.
• Fitting the RMT prediction for the distribution of dimensionless eigenmodes (ζ = ρλₖΣV) to measured data to extract the chiral condensate Σ.

Experimental results

Research questions

• RQ1 Can the use of stout-smeared gauge links significantly reduce the computational cost of dynamical overlap fermion simulations?
• RQ2 How accurately can the chiral condensate be extracted from the distribution of low-lying Dirac eigenmodes using Random Matrix Theory?
• RQ3 Does the RMT prediction for the eigenmode distribution hold in finite-volume simulations with dynamical overlap fermions?
• RQ4 What is the behavior of the system in restricted topological sectors, and is the simulation ergodic when topology changes are forbidden?
• RQ5 How do the results for the chiral condensate vary with quark mass, and what is the value in the chiral limit?

Key findings

• The use of stout-smeared gauge links reduced the number of Dirac operator matrix-vector multiplications per trajectory by approximately an order of magnitude.
• The fermion force was reduced by about an order of magnitude, enabling efficient multiple time-scale integration.
• The RMT prediction for the distribution of the two lowest Dirac eigenmodes matched the measured data well, particularly in the ν = 0 and ν = ±1 sectors.
• The chiral condensate was extracted as Σ ≈ (280 MeV)³ using the physical pion decay constant and a finite-volume correction factor ρ ≈ 1.4.
• Simulations at amₚ = 0.01 and amₚ = 0.005 ran stably in the ν = 0 sector, with no signs of non-ergodicity or systematic biases.
• The distribution of the third eigenmode showed good agreement with RMT predictions, though slight deviations were observed in the |ν| = 1 sector, likely due to the Thouless energy scale being approached.