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, https://bib-pubdb1.desy.de/search?p=id:%22PUBDB-2023-05897%22

zbMATH Keywords

conjugate gradient methodconvergence accelerationDirac operatoreigenvectorsgauge fieldslow-lying eigenvaluessparse Hermitian matrixWilson fermions

Mathematics Subject Classification ID

PDEs in connection with quantum mechanics (35Q40) Numerical computation of eigenvalues and eigenvectors of matrices (65F15) 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) Quantum field theory on lattices (81T25) Strong interaction, including quantum chromodynamics (81V05)

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

• Almost diagonal matrices with multiple or close eigenvalues
• Computing eigenvalues of very large symmetric matrices. An implementation of a Lanczos algorithm with no reorthogonalization
• 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?)
• Optimal gradient minimization scheme for finite element eigenproblems
• The computational efficiency of a new minimization algorithm for eigenvalue analysis
• The Jacobi method for real symmetric matrices

Cited in

(25)

• 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
• Dynamical fermions as a global correction
• Locality properties of Neuberger's lattice Dirac operator
• Monopole and vortex content of a meron pair
• Overlap fermions, improved cooling and the lattice index theorem
• Triangularized orthogonalization-free method for solving extreme eigenvalue problems
• 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
• Quark Spectra, Topology, and Random Matrix Theory
• Ultraviolet filtering of lattice configurations and applications to Monte Carlo dynamics
• Numerical determination of partial spectrum of Hermitian matrices using a Lánczos method with selective reorthogonalization
• Effective noise reduction techniques for disconnected loops in lattice QCD
• On Witten's global anomaly for higher SU(2) representations
• On Witten's global anomaly for higher SU(2) representations
• Exact zero-modes of the compact QED Dirac operator
• SU(2) vortex configuration in Laplacian center gauge
• Low-lying eigenmodes of the Wilson-Dirac operator and correlations with topological objects
• Spectrum of the Hermitian Wilson Dirac operator
• Low-lying Wilson Dirac operator eigenvector mixing in dynamical overlap hybrid Monte Carlo
• Numerical methods for the QCD overlap operator IV: Hybrid Monte Carlo
• Numerical techniques for lattice QCD in the \(\epsilon \)-regime
• Understanding the approach to thermalization from the eigenspectrum of non-abelian gauge theories
• Chiral symmetry breaking at large N_c

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)