An iterative method to compute the sign function of a non-Hermitian matrix and its application to the overlap Dirac operator at nonzero chemical potential

DOI10.1016/J.CPC.2007.07.012zbMATH Open1196.65085arXiv0704.3486OpenAlexW2012926860MaRDI QIDQ710108FDOQ710108

J. C. R. Bloch, Bruno Lang, T. Wettig, Andreas Frommer

Publication date: 18 October 2010

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

Abstract: The overlap Dirac operator in lattice QCD requires the computation of the sign function of a matrix. While this matrix is usually Hermitian, it becomes non-Hermitian in the presence of a quark chemical potential. We show how the action of the sign function of a non-Hermitian matrix on an arbitrary vector can be computed efficiently on large lattices by an iterative method. A Krylov subspace approximation based on the Arnoldi algorithm is described for the evaluation of a generic matrix function. The efficiency of the method is spoiled when the matrix has eigenvalues close to a function discontinuity. This is cured by adding a small number of critical eigenvectors to the Krylov subspace, for which we propose two different deflation schemes. The ensuing modified Arnoldi method is then applied to the sign function, which has a discontinuity along the imaginary axis. The numerical results clearly show the improved efficiency of the method. Our modification is particularly effective when the action of the sign function of the same matrix has to be computed many times on different vectors, e.g., if the overlap Dirac operator is inverted using an iterative method.

Full work available at URL: https://arxiv.org/abs/0704.3486

Recommendations

• Short-recurrence Krylov subspace methods for the overlap Dirac operator at nonzero chemical potential
• 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 nested Krylov subspace method to compute the sign function of large complex matrices
• Numerical methods for the QCDd overlap operator. I: Sign-function and error bounds
• Numerical methods for the QCD overlap operator IV: Hybrid Monte Carlo

zbMATH Keywords

iterative methodsnon-Hermitian matrixsign functionoverlap Dirac operatorquark chemical potential

Mathematics Subject Classification ID

Strong interaction, including quantum chromodynamics (81V05) Quantum field theory on lattices (81T25) Numerical computation of matrix exponential and similar matrix functions (65F60)

Cites Work

• Algorithm 432 [C2]: Solution of the matrix equation AX + XB = C [F4]
• Title not available (Why is that?)
• Title not available (Why is that?)
• Title not available (Why is that?)
• A Hessenberg-Schur method for the problem AX + XB= C
• Linear model reduction and solution of the algebraic Riccati equation by use of the sign function†
• Analysis of Some Krylov Subspace Approximations to the Matrix Exponential Operator
• Using Nonorthogonal Lanczos Vectors in the Computation of Matrix Functions
• Numerical methods for the QCDd overlap operator. I: Sign-function and error bounds
• An iterative solution method for solving \(f(A)x=b\), using Krylov subspace information obtained for the symmetric positive definite matrix A
• Stable iterations for the matrix square root
• Rational Iterative Methods for the Matrix Sign Function
• Overlap Dirac Operator at Nonzero Chemical Potential and Random Matrix Theory
• The matrix sign function and computations in systems
• The matrix sign function
• Formal derivation of algorithms
• A Lanczos approach to the inverse square root of a large and sparse matrix
• Analysis of Augmented Krylov Subspace Methods
• A parallel block cyclic reduction algorithm for the fast solution of elliptic equations

Cited In (14)

• A Comparison of Limited-memory Krylov Methods for Stieltjes Functions of Hermitian Matrices
• Overlap Dirac Operator at Nonzero Chemical Potential and Random Matrix Theory
• Error bounds and estimates for Krylov subspace approximations of Stieltjes matrix functions
• Short-recurrence Krylov subspace methods for the overlap Dirac operator at nonzero chemical potential
• The Radau--Lanczos Method for Matrix Functions
• Low-Rank Updates of Matrix Functions II: Rational Krylov Methods
• Title not available (Why is that?)
• On restarted and deflated block FOM and GMRES methods for sequences of shifted linear systems
• Analysis of Probing Techniques for Sparse Approximation and Trace Estimation of Decaying Matrix Functions
• A nested Krylov subspace method to compute the sign function of large complex matrices
• BiCR-type methods for families of shifted linear systems
• Analyzing the spectrum of general, non-hermitian Dirac operators
• Numerical methods for the QCDd overlap operator. I: Sign-function and error bounds
• Restarted Hessenberg method for solving shifted nonsymmetric linear systems

This page was built for publication: An iterative method to compute the sign function of a non-Hermitian matrix and its application to the overlap Dirac operator at nonzero chemical potential

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