An accelerated conjugate gradient algorithm to compute low-lying eigenvalues -- a study for the Dirac operator in SU(2) lattice QCD

DOI10.1016/0010-4655(95)00126-3MaRDI QIDQ1276552zbMATH OpenOpenAlexFDO

Authors Thomas Kalkreuter, Hubert Simma

Publication date 26 May 1999

Published in Computer Physics Communications (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/hep-lat/9507023

zbMATH Keywords

eigenvectorsconvergence accelerationconjugate gradient methodDirac operatorgauge fieldsWilson fermionslow-lying eigenvaluessparse Hermitian matrix

Mathematics Subject Classification ID

Numerical computation of eigenvalues and eigenvectors of matrices (65F15) PDEs in connection with quantum mechanics (35Q40) Closed and approximate solutions to the Schrödinger, Dirac, Klein-Gordon and other equations of quantum mechanics (81Q05) Yang-Mills and other gauge theories in quantum field theory (81T13) Strong interaction, including quantum chromodynamics (81V05) Quantum field theory on lattices (81T25)

Abstract: The low-lying eigenvalues of a (sparse) hermitian matrix can be computed with controlled numerical errors by a conjugate gradient (CG) method. This CG algorithm is accelerated by alternating it with exact diagonalisations in the subspace spanned by the numerically computed eigenvectors. We study this combined algorithm in case of the Dirac operator with (dynamical) Wilson fermions in four-dimensional $SUtwo$ gauge fields. The algorithm is numerically very stable and can be parallelized in an efficient way. On lattices of sizes $4^4-16^4$ an acceleration of the pure CG method by a factor of~$4-8$ is found.

Recommendations

• Computing and deflating eigenvalues while solving multiple right-hand side linear systems with an application to quantum chromodynamics
• Acceleration of the Arnoldi method and real eigenvalues of the non-Hermitian Wilson-Dirac operator
• Quark propagator on the Connection Machine
• scientific article; zbMATH DE number 1931072

Cites work

• scientific article; zbMATH DE number 3584841 (Why is no real title available?)
• scientific article; zbMATH DE number 3622441 (Why is no real title available?)
• scientific article; zbMATH DE number 274379 (Why is no real title available?)
• scientific article; zbMATH DE number 3225665 (Why is no real title available?)
• scientific article; zbMATH DE number 3407464 (Why is no real title available?)
• Almost diagonal matrices with multiple or close eigenvalues
• Computing eigenvalues of very large symmetric matrices. An implementation of a Lanczos algorithm with no reorthogonalization
• Optimal gradient minimization scheme for finite element eigenproblems
• The Jacobi method for real symmetric matrices
• The computational efficiency of a new minimization algorithm for eigenvalue analysis

Cited in

(24)

• Quark Spectra, Topology, and Random Matrix Theory
• Ultraviolet filtering of lattice configurations and applications to Monte Carlo dynamics
• SU(2) vortex configuration in Laplacian center gauge
• Locality properties of Neuberger's lattice Dirac operator
• A quadrature-based eigensolver with a Krylov subspace method for shifted linear systems for Hermitian eigenproblems in lattice QCD
• Global anomalies in chiral lattice gauge theory
• Monopole and vortex content of a meron pair
• Numerical techniques for lattice QCD in the \(\epsilon \)-regime
• Low-lying Wilson Dirac operator eigenvector mixing in dynamical overlap hybrid Monte Carlo
• Numerical methods for the QCD overlap operator IV: Hybrid Monte Carlo
• Chiral symmetry breaking at large \(N_c\)
• Analyzing the spectrum of general, non-hermitian Dirac operators
• The low-lying Dirac eigenmodes from domain wall fermions
• Understanding chiral symmetry breaking with the overlap action
• Overlap fermions, improved cooling and the lattice index theorem
• Triangularized orthogonalization-free method for solving extreme eigenvalue problems
• Exact zero-modes of the compact QED Dirac operator
• Low-lying eigenmodes of the Wilson-Dirac operator and correlations with topological objects
• On Witten's global anomaly for higher SU(2) representations
• Effective noise reduction techniques for disconnected loops in lattice QCD
• On Witten's global anomaly for higher SU(2) representations
• Spectrum of the Hermitian Wilson Dirac operator
• Dynamical fermions as a global correction
• Numerical determination of partial spectrum of Hermitian matrices using a Lánczos method with selective reorthogonalization

This page was built for publication: An accelerated conjugate gradient algorithm to compute low-lying eigenvalues -- a study for the Dirac operator in SU(2) lattice QCD

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q1276552)