一种计算非零化学势下重叠 Dirac 算符的迭代方法

非零夸克化学势下的重叠 Dirac 算符涉及非厄米矩阵符号函数的计算。本文介绍一种迭代方法（最初由我们在文献 [1] 中提出），该方法即使在大格点上也能高效计算该算符。出发点是基于 Arnoldi 算法的 Krylov 子空间近似，用于评估一般矩阵函数。当矩阵具有接近函数不连续点的本征值时，该方法的效率会降低。为解决此问题，我们将少量临界本征向量加入 Krylov 子空间，并在此增广子空间中提出了两种不同的降阶方案。随后将该方法应用于两种不同格点尺寸下的重叠 Dirac 算符符号函数计算。符号函数沿虚轴存在不连续点，数值结果表明降阶技术显著提高了该方法的效率。