[论文解读] Improvement of algorithms for dynamical overlap fermions

JLQCD Collaboration, Hideo Matsufuru+8 more|ArXiv.org|Oct 3, 2006

Quantum Chromodynamics and Particle Interactions参考文献 7被引用 6

一句话总结

该论文通过结合5D共轭梯度求解器、多时间步积分和噪声Metropolis测试，改进了动态重叠费米子的混合蒙特卡洛（HMC）算法，在Blue Gene系统上实现了相较于标准方法3倍的性能提升。该方法通过引入额外的威尔逊费米子抑制近零模，并优化求解器精度，在保持精度的同时降低了数值成本。

ABSTRACT

We investigate the algorithms for dynamical overlap fermions aiming at improving the performance for large-scale simulations. We look for the best combination of Hybrid Monte Carlo options and iterative quark solvers with respect to the numerical costs. Our main target is a $N_f=2$ simulation with overlap fermion on a $16^3 imes 32$ lattice at lattice spacing around 0.12 fm.

研究动机与目标

• 提升大规模格点QCD中动态重叠费米子模拟的性能。
• 降低HMC模拟中最昂贵的部分——重叠狄拉克算符求逆的计算成本。
• 优化HMC算法参数，包括时间步长和求解器精度，以提高效率。
• 使现代超级计算机上实现包含两个轻夸克和小夸克质量的大体积模拟成为可能。
• 通过引入额外的威尔逊费米子抑制重叠算符的近零本征值，确保算法的稳定性和精度。

提出的方法

• 采用带有部分分式近似的5D共轭梯度求解器，计算重叠狄拉克算符中符号函数的值。
• 实施多时间步积分，根据力的大小层次，为规范项、$S_{PF1}$、$S_{PF2}$和$S_E$项分别使用不同的时间步长。
• 使用更重的夸克质量$m'$进行质量预条件处理，以加速费米子求解器的收敛。
• 引入具有扭曲质量$\mu=0.2$的额外威尔逊费米子，抑制$H_W$的近零模，从而消除对反射/折射检查的需求。
• 引入噪声Metropolis测试以校正5D求解器中的不准确性，以概率$P = \min\{1, e^{-dS}\}$接受新配置。
• 优化求解器参数，包括极点数$N$、$m'$以及步长比值，以在精度与性能之间取得平衡。

实验结果

研究问题

• RQ1 在动态重叠费米子模拟中，如何最优组合HMC参数与迭代求解器，以最小化计算成本？
• RQ2 引入额外威尔逊费米子对HMC算法的稳定性和效率有何影响？
• RQ3 5D共轭梯度求解器在速度和精度方面是否优于标准CG和松弛CG求解器？
• RQ4 多时间步积分对具有层级力项的HMC算法的稳定性和性能有何影响？
• RQ5 噪声Metropolis测试在多大程度上可补偿5D求解器的不准确性，同时保持物理精度？

主要发现

• 使用$N=20$个极点的5D共轭梯度求解器，其速度比相同精度下的松弛CG方法快2至3倍。
• 采用$\Delta\tau_{(PF2)}/\Delta\tau_{(PF1)}=5$和$\Delta\tau_{(PF1)}/\Delta\tau_{(G,E)}=6$的多时间步方法显著提升了效率。
• 噪声Metropolis测试使得可以使用精度较低的5D求解器，同时保持物理精度，从而减少了对高精度求逆的需求。
• 5D求解器、多时间步和噪声Metropolis测试的组合，在Blue Gene（512节点）上相比标准嵌套CG实现了3倍提速。
• 具有$\mu=0.2$的额外威尔逊费米子项能有效抑制$H_W$的近零模，从而消除对反射/折射过程的需求。
• 在Blue Gene（512节点）上的性能表现显示，使用改进算法后，轨迹时间从标准嵌套CG的112分钟降低至22分钟。