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

Bloch, Jacques and Frommer, A. and Lang, B. and Wettig, T. (2007) 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. COMPUTER PHYSICS COMMUNICATIONS, 177 (12). pp. 933-943. ISSN 0010-4655,

Full text not available from this repository. (Request a copy)

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. (c) 2007 Elsevier B.V. All rights reserved.

Item Type: Article
Uncontrolled Keywords: CHIRAL FERMIONS; SQUARE-ROOT; LATTICE; EQUATION; SYMMETRY; overlap Dirac operator; quark chemical potential; sign function; non-Hermitian matrix; iterative methods
Subjects: 500 Science > 530 Physics
Divisions: Physics > Institute of Theroretical Physics
Depositing User: Dr. Gernot Deinzer
Date Deposited: 25 Nov 2020 06:46
Last Modified: 25 Nov 2020 06:46
URI: https://pred.uni-regensburg.de/id/eprint/31804

Actions (login required)

View Item View Item