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-3zbMATH Open0917.65036arXivhep-lat/9507023OpenAlexW3124130186MaRDI QIDQ1276552FDOQ1276552

Thomas Kalkreuter, Hubert Simma

Publication date: 26 May 1999

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

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.

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

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

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)

Cites Work

• Title not available (Why is that?)
• Title not available (Why is that?)
• Title not available (Why is that?)
• Computing eigenvalues of very large symmetric matrices. An implementation of a Lanczos algorithm with no reorthogonalization
• The computational efficiency of a new minimization algorithm for eigenvalue analysis
• The Jacobi method for real symmetric matrices
• Title not available (Why is that?)
• Almost diagonal matrices with multiple or close eigenvalues
• Title not available (Why is that?)
• Optimal gradient minimization scheme for finite element eigenproblems

Cited In (23)

• Low-lying eigenmodes of the Wilson-Dirac operator and correlations with topological objects
• Triangularized orthogonalization-free method for solving extreme eigenvalue problems
• Effective noise reduction techniques for disconnected loops in lattice QCD
• Locality properties of Neuberger's lattice Dirac operator
• Overlap fermions, improved cooling and the lattice index theorem
• Spectrum of the Hermitian Wilson Dirac operator
• Monopole and vortex content of a meron pair
• Quark Spectra, Topology, and Random Matrix Theory
• Low-lying Wilson Dirac operator eigenvector mixing in dynamical overlap hybrid Monte Carlo
• Numerical methods for the QCD overlap operator IV: Hybrid Monte Carlo
• Exact zero-modes of the compact QED Dirac operator
• SU(2) vortex configuration in Laplacian center gauge
• On Witten's global anomaly for higher SU(2) representations
• Numerical techniques for lattice QCD in the \(\epsilon \)-regime
• Ultraviolet filtering of lattice configurations and applications to Monte Carlo dynamics
• Dynamical fermions as a global correction
• On Witten's global anomaly for higher SU(2) representations
• 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
• Global anomalies in chiral lattice gauge theory
• 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)