[论文解读] Disconnected contributions to hadronic structure: a new method for stochastic noise reduction

Sara Collins, Gunnar Bali+1 more|arXiv (Cornell University)|Sep 20, 2007

Particle physics theoretical and experimental studies参考文献 9被引用 18

一句话总结

该论文提出了一种新颖且无偏的方法，通过在共轭梯度求逆未完全收敛时截断，以减少格点QCD计算中非连通夸克圈的随机噪声。该方法通过使用廉价的局部求逆来实现大部分噪声抑制，仅用少数几次完整求逆来校正误差，从而在计算机时间上实现了4–12倍的加速，若结合跳跃参数展开法，甚至可实现高达30倍的加速，显著提升了计算强子结构可观测量的可行性。

ABSTRACT

We present a new method for reducing the stochastic noise of all-to-all propagators based on stopping the inversion of the propagator before convergence. The method is easy to implement, unbiased and independent of the quark action. Applying this method to the calculation of disconnected loops needed for hadronic structure observables we find savings in computer time of factors of 4-12 depending on the operator inserted in the loop. When combined with a hopping parameter expansion technique we obtain combined gains of up to factors of 30 for some operators.

研究动机与目标

• 为解决强子结构可观测量的非连通夸克圈计算中高计算成本和随机噪声的问题。
• 开发一种无偏、适用于所有夸克作用形式且易于实现的噪声抑制技术。
• 提升计算核子形式因子和部分子分布矩的非连通贡献的可行性。
• 将新方法与现有技术（如跳跃参数展开）结合，以进一步提升性能。

提出的方法

• 该方法在共轭梯度求解器迭代次数较少（$ n_t \ll n_c $）时截断，生成近似逆矩阵 $ \mathrm{M}^{-1}_{n_t} $，显著降低计算成本。
• 通过将 $ L_1 $ 次局部求逆与 $ L_2 $ 次校正项 $ \mathrm{M}^{-1}_{n_c} - \mathrm{M}^{-1}_{n_t} $ 的完整求逆相结合，构建全逆的无偏估计。
• 两组噪声向量相互独立，确保整体估计保持无偏，且随机误差按 $ \mathcal{O}(1/\sqrt{L}) $ 缩放。
• 在固定计算成本下，通过最小化总方差来确定最优的 $ n_t $ 和 $ L_1/L_2 $ 比例。
• 将该方法与跳跃参数展开（HPE）结合，进一步降低方差，尤其在重夸克区域能效显著。
• 该方法在共轭梯度（CG）和BiCGStab求解器上均进行了测试，结果表现出一致的性能提升。

实验结果

研究问题

• RQ1 截断求解器方法是否能在不引入偏差的前提下，减少全通量夸克传播子中的随机噪声？
• RQ2 在给定方差下，如何在早期截断与校正求逆之间实现计算成本的最优权衡？
• RQ3 该方法在不同狄拉克结构 $ \Gamma $（如 $ \mathbb{1} $、$ \gamma_5 $ 和 $ \gamma_3\gamma_5 $）下的表现如何？
• RQ4 该方法在多大程度上可与跳跃参数展开结合以进一步降低方差？
• RQ5 随着夸克质量减小，特别是接近奇异夸克和轻夸克区域时，该方法的性能如何变化？

主要发现

• 对于接近奇异夸克质量的夸克质量，该截断求解器方法在非连通圈计算中实现了4–12倍的计算机时间减少，具体取决于算符 $ \Gamma $。
• 当与跳跃参数展开结合使用时，对于 $ \Gamma = \gamma_3\gamma_5 $，该方法可实现高达30倍的综合加速。
• 最优截断点 $ n_t^{\text{opt}} $ 在引入HPE后有所增加，但仍远小于完全收敛所需的 $ n_c $，表明噪声抑制效果显著。
• 该方法在 $ m_{PS} \approx 300 \, \text{MeV} $ 时仍保持有效性，对多个 $ \Gamma $ 算符，其性能较标准随机方法至少提升2倍。
• 使用BiCGStab替代CG求解器可获得相当的性能提升（例如，对于 $ \Gamma = \gamma_3\gamma_5 $，分别实现9倍和24倍加速），表明该方法在不同求解器间具有鲁棒性。
• 该方法计算效率高，应用 $ \not\!D^4 $ 的成本仅占完整求逆的5%，使得基于HPE的校正成本可控。