Short-recurrence Krylov subspace methods for the overlap Dirac operator at nonzero chemical potential

详细信息    查看全文

• 作者：Jacques C.R. Bloch ; Tobias Breu ; Andreas Frommer ; Simon Heybrock ; Katrin Sch&#xe4 ; fer ; Tilo Wettig
• 关键词：Overlap Dirac operator ; Quark chemical potential ; Sign function ; Non-Hermitian matrix ; Iterative methods
• 刊名：Computer Physics Communications
• 出版年：2010
• 出版时间：August 2010
• 年：2010
• 卷：181
• 期：8
• 页码：1378-1387
• 全文大小：1255 K

文摘

The overlap operator in lattice QCD requires the computation of the sign function of a matrix, which is non-Hermitian in the presence of a quark chemical potential. In previous work we introduced an Arnoldi-based Krylov subspace approximation, which uses long recurrences. Even after the deflation of critical eigenvalues, the low efficiency of the method restricts its application to small lattices. Here we propose new short-recurrence methods which strongly enhance the efficiency of the computational method. Using rational approximations to the sign function we introduce two variants, based on the restarted Arnoldi process and on the two-sided Lanczos method, respectively, which become very efficient when combined with multishift solvers. Alternatively, in the variant based on the two-sided Lanczos method the sign function can be evaluated directly. We present numerical results which compare the efficiencies of a restarted Arnoldi-based method and the direct two-sided Lanczos approximation for various lattice sizes. We also show that our new methods gain substantially when combined with deflation.