Sign of Wilson-Fermion Determinant using Contour-Integral Spectral Projection

技术摘要：利用轮廓积分谱投影确定 Wilson-Fermion 行列式符号

问题陈述
在利用 Wilson 费米子的格点量子色动力学（Lattice QCD）中，费米子行列式的符号 $\text{sign}(\det D)$ 是一个关键量。虽然 $\gamma_5$-Hermiticity 确保了行列式是实数，但对于单味（single-flavor）或质量非简并的 Wilson 费米子，其正定性并不能得到保证。该符号完全由非厄米 Wilson-Dirac 算符（$D$）位于负实轴上的特征值数量 $N_-$ 的奇偶性决定。确定这一符号具有计算挑战性，特别是在存在近零模或“异常构型”（exceptional configurations）的参数区域，即算符处于病态条件下的情况。现有方法，如谱流法（通过追踪厄米算符 $Q=\gamma_5 D$ 随质量变化的特征值）或对非厄米矩阵进行直接迭代特征值求解，通常需要解析宽广的低能谱区域，或者依赖于需要重加权因子（$W_{\text{sign}}$）来修正累积误差的多项式/有理近似。

方法论
作者提出了一种基于轮廓积分的有针对性的谱投影方法，旨在直接统计负实特征值，从而绕过对整个低能谱的解析需求。该方法的核心组成部分包括：

1. 轮廓积分谱投影：
   该方法利用在复 Argand 平面中关于实轴对称的轮ло（ellipse）轮廓 $\Gamma$ 来隔离感兴趣的谱区域。谱投影器定义为：
   $$P_\Gamma = \frac{1}{2\pi i} \oint_\Gamma (zI - A)^{-1} dz$$
   其中 $A$ 是非厄米 Wilson-Dirac 算符。该积分通过使用有限的正交规则（包含 $N_q$ 个节点 $z_j$ 和权重 $\omega_j$）进行数值近似。生成的近似投影器 $\hat{P}_\Gamma$ 将算符映射到一个缩减后的子空间。

2. 利用 $\gamma_5$-Hermiticity：
   由于 Wilson-Dirac 算符具有 $\gamma_5$-Hermiticity，左、右特征系是相关的（即 $r_k^\dagger \gamma_5$ 作为 $\lambda_k^*$ 的左特征向量）。因此，该方法仅需求解右投影子空间，避免了独立确定左特征系的需要。

3. 移位线性系统与求解器适配：
   投影过程涉及在每个正交节点处求解移位线性系统 $(z_j I - \hat{D})q_j = Y$，其中 $\hat{D}$ 是 $D$ 的偶 Schur 补。作者对原本为全算符 $D$ 设计的 DFL-SAP-GCR 求解器进行了适配，使其专门作用于移位的偶 Schur 补 $\hat{D} - zI$。这种适配至关重要，因为移位后的全算符的 Schur 补并不等于原算符的移位 Schur 补。

4. 用于谱边界限制的奇异值排除法：
   一个关键挑战是确保轮廓的左边界（$X_L$）既能包含所有负实特征值，又不会因过大而导致计算成本激增。作者引入了一种基于奇异值估计的务实策略：通过评估一系列点 $s_i$ 处移位系统 $\hat{D} - s_i I$ 的最小奇异值 $\sigma_{\min}$，识别出不含特征值的排除区间。这使得能够确定一个保守且高效的 $X_L$，避免了仅基于厄米部分 Rayleigh 商所得到的过于保守的界限。

5. 子空间迭代：
   采用 FEAST 软件包来实现试探子空间向量。试探子空间的维度（$m$）被设定为超过预期包含的特征值数量。该方法表现出快速收敛性，通常在两次 FEAST 迭代内即可解析所有负实模。

关键结果
该方法在来自 CLS 财团生成的各类格点 QCD 集合上进行了测试。主要发现包括：

• 鲁棒性与准确性： 该程序成功地在各种构型下重现了 $\hat{D}$ 的所有负实解，精度达到机器精度。这些构型涵盖了无负解、浅负解、深负解以及多个解之间相差数个数量级的案例。
• 参数效率： 对于大多数构型，仅需 $N_q = 8$ 个正交节点和 $10^{-3}$ 的移位求解器容差即可满足要求。所有负实特征值均在第二次 FEAST 迭代后得到解析。
• 处理极端情况： 在特征值跨越两个数量级的构型中（例如配置 C102r005n56），通过将 $N_q$ 增加到 12 并将求解器容差收紧至 $10^{-5}$，可以在不产生显著计算开销的情况下，忠实地捕获所有模。
• 计算成本： 主要成本源于求解移位线性系统。选择轮廓和正交节点的额外开销不到总评估成本的百分之几。奇异值排除策略有效地最小化了积分区域，减少了所需的正交节点数，相比于使用仅基于厄米部分导出的界限。

意义与主张
论文声称，这种基于轮廓积分的谱投影提供了一种系统且高效的方法，用于确定费米子行列式符号，而无需承担解析无关物理模的计算负担，也不依赖于宽广的谱流分析。通过将特征值求解器限制在控制行列式符号的确切定义域内，该方法为传统的补充了多项式变换的迭代特征值求解器提供了一个鲁棒的替代方案。

除了在评估格点 QCD 中 $W_{\text{sign}}$ 的特定应用外，作者还指出该方法提供了一种通用的策略，用于从稠密的复谱中隔离具有物理意义的特征值。这对于非厄米量子系统、稳定性分析以及大规模特征值问题具有潜在的应用价值，能够作为针对特定谱子集的定向过滤器，而非进行宽泛的搜索。这项工作强调，可靠地确定行列式符号对于通过重加权纠正有理近似误差（用于模拟）或验证单味情况下测度的正定性至关重要。