An iterative method to compute the overlap Dirac operator at nonzero chemical potential
نویسندگان
چکیده
The overlap Dirac operator at nonzero quark chemical potential involves the computation of the sign function of a non-Hermitian matrix. In this talk we present an iterative method, first proposed by us in Ref. [1], which allows for an efficient computation of the operator, even on large lattices. The starting point is a Krylov subspace approximation, based on the Arnoldi algorithm, for the evaluation of a generic matrix function. The efficiency of this method is spoiled when the matrix has eigenvalues close to a function discontinuity. To cure this, a small number of critical eigenvectors are added to the Krylov subspace, and two different deflation schemes are proposed in this augmented subspace. The ensuing method is then applied to the sign function of the overlap Dirac operator, for two different lattice sizes. The sign function has a discontinuity along the imaginary axis, and the numerical results show how deflation dramatically improves the efficiency of the method.
منابع مشابه
Comparing iterative methods to compute the overlap Dirac operator at nonzero chemical potential
The overlap Dirac operator at nonzero quark chemical potential involves the computation of the sign function of a non-Hermitian matrix. In this talk we present iterative Krylov subspace approximations, with deflation of critical eigenvalues, which we developed to compute the operator on large lattices. We compare the accuracy and efficiency of two alternative approximations based on the Arnoldi...
متن کاملA nested Krylov subspace method for the overlap operator
We present a novel method to compute the overlap Dirac operator at zero and nonzero quark chemical potential. To approximate the sign function of large, sparse matrices, standard methods project the operator on a much smaller Krylov subspace, on which the matrix function is computed exactly. However, for large lattices this subspace can still be too large for an efficient calculation of the sig...
متن کامل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
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 app...
متن کاملOverlap Dirac operator at nonzero chemical potential and random matrix theory.
We show how to introduce a quark chemical potential in the overlap Dirac operator. The resulting operator satisfies a Ginsparg-Wilson relation and has exact zero modes. It is no longer gamma5 Hermitian, but its nonreal eigenvalues still occur in pairs. We compute the spectral density of the operator on the lattice and show that, for small eigenvalues, the data agree with analytical predictions ...
متن کاملar X iv : h ep - l at / 0 60 90 20 v 1 1 3 Se p 20 06 Nonzero chemical potential in the overlap Dirac operator and comparison to random matrix theory
In this talk we present the results published recently in Ref. [1], where we showed how to introduce a quark chemical potential in the overlap Dirac operator. The resulting operator satisfies a Ginsparg-Wilson relation and has exact zero modes. It is no longer γ5-Hermitian, but its nonreal eigenvalues still occur in pairs. We compute the spectral density of the operator on the lattice and show ...
متن کامل