Numerical Methods for the QCD Overlap Operator: I. Sign-Function and Error Bounds

Abstract

The numerical and computational aspects of the overlap formalism in lattice quantum chromodynamics are extremely demanding due to a matrix-vector product that involves the sign function of the hermitian Wilson matrix. In this paper we investigate several methods to compute the product of the matrix sign-function with a vector, in particular Lanczos based methods and partial fraction expansion methods. Our goal is two-fold: we give realistic comparisons between known methods together with novel approaches and we present error bounds which allow to guarantee a given accuracy when terminating the Lanczos method and the multishift-CG solver, applied within the partial fraction expansion methods.

30 pages, 2 figures

Cited by in corpus (68)

• Full nonperturbative QCD simulations with 2+1 flavors of improved staggered quarks
• Domain wall QCD with physical quark masses
• Light Hadron Masses from Lattice QCD
• Numerical techniques for lattice QCD in the $\epsilon$--regime
• Dynamical overlap fermions, results with hybrid Monte-Carlo algorithm
• Scale-invariance of parity-invariant three-dimensional QED
• Non-equilibrium study of the Chiral Magnetic Effect from real-time simulations with dynamical fermions
• Overlap Valence on 2+1 Flavor Domain Wall Fermion Configurations with Deflation and Low-mode Substitution
• Two-flavor lattice QCD in the epsilon-regime and chiral Random Matrix Theory
• Deflated GMRES for Systems with Multiple Shifts and Multiple Right-Hand Sides
• A note on Zolotarev optimal rational approximation for the overlap Dirac operator
• Axial Correlation Functions in the epsilon-Regime: a Numerical Study with Overlap Fermions
• Chiral symmetry breaking at large N_c
• QCD thermodynamics with dynamical overlap fermions
• Chiral Logs in Quenched QCD
• Quark propagator in Landau and Laplacian gauges with overlap fermions
• Numerical Methods for the QCD Overlap Operator:III. Nested Iterations
• Algorithms for Lattice QCD with Dynamical Fermions
• Numerical Methods for the QCD Overlap Operator IV: Hybrid Monte Carlo
• Filtered overlap: speedup, locality, kernel non-normality and Z_A~1
• Topological susceptibility from the overlap
• Scaling behavior of the overlap quark propagator in Landau gauge
• An iterative method to compute the overlap Dirac operator at nonzero chemical potential
• Light hadron and diquark spectroscopy in quenched QCD with overlap quarks on a large lattice
• First-order restoration of SU(Nf) x SU(Nf) chiral symmetry with large Nf and Electroweak phase transition
• Lattice QCD at the end of 2003
• A nested Krylov subspace method to compute the sign function of large complex matrices
• Low-lying Wilson Dirac operator eigenvector mixing in dynamical overlap Hybrid Monte-Carlo
• Lattice QCD at finite T and μ
• Topological tunneling with Dynamical overlap fermions
• A Numerical Study of the 2-Flavour Schwinger Model with Dynamical Overlap Hypercube Fermions
• Short-recurrence Krylov subspace methods for the overlap Dirac operator at nonzero chemical potential
• Limited-memory polynomial methods for large-scale matrix functions
• Screening correlators with chiral Fermions
• Nonperturbative renormalization of composite operators with overlap fermions
• The Padé iterations for the matrix sign function and their reciprocals are optimal
• 2-norm error bounds and estimates for Lanczos approximations to linear systems and rational matrix functions
• Magnetic catalysis in the (2+1)-dimensional Gross-Neveu model
• Current status of Dynamical Overlap project
• Nonperturbative renormalisation of composite operators with overlap quarks
• Upper and lower Higgs boson mass bounds from a chirally invariant lattice Higgs-Yukawa model
• Error bounds for Lanczos-based matrix function approximation
• QCD thermodynamics with continuum extrapolated dynamical overlap fermions
• Principal Component Projection Without Principal Component Analysis
• A Posteriori Error Estimates of Krylov Subspace Approximations to Matrix Functions
• Optimised Dirac Operators on the Lattice: Construction, Properties and Applications
• Spectral properties and chiral symmetry violations of (staggered) domain wall fermions in the Schwinger model
• Dynamical overlap fermions, results with HMC algorithm
• An Exact Algorithm for Any-flavor Lattice QCD with Kogut-Susskind Fermion
• A numerical solution to the local cohomology problem in U(1) chiral gauge theories
• Improvement of algorithms for dynamical overlap fermions
• Low-memory Krylov subspace methods for optimal rational matrix function approximation
• Exploring the chiral regime with dynamical overlap fermions
• Faster Principal Component Regression and Stable Matrix Chebyshev Approximation
• Speed and Adaptability of Overlap Fermion Algorithms
• A numerical method to compute derivatives of functions of large complex matrices and its application to the overlap Dirac operator at finite chemical potential
• A lattice Dirac operator for QCD with light dynamical quarks
• Preconditioning the Restarted and Shifted Block FOM Algorithm for Matrix Exponential Computation
• Multigrid Preconditioning for the Overlap Operator in Lattice QCD
• The eigSUMR inverter for overlap fermion
• A curious behavior of three-dimensional lattice Dirac operators coupled to monopole background
• Improving Inversions of the Overlap Operator
• Error Bounds for the Sign Function
• A posteriori error bounds for the block-Lanczos method for matrix function approximation
• A SQMRCGstab Algorithm for Families of Shifted Linear Systems
• Dynamical fermions in lattice quantum chromodynamics
• Continuum reduction in large N gauge theories
• Lattice QCD and the Balkan physicists contribution