[论文解读] Improving DWF Simulations: the Force Gradient Integrator and the Möbius Accelerated DWF Solver

Hantao Yin, Robert D. Mawhinney|arXiv (Cornell University)|Nov 21, 2011

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

一句话总结

本文提出了一种力梯度积分器（FGI）和一种Möbius加速的域墙费米子（DWF）求解器，以提升格点QCD模拟的效率。FGI降低了积分误差并允许使用更大的步长，而Möbius求解器通过减小第五维尺寸，将狄拉克算符应用次数减少约2倍，在32³×64体积、π介子质量约为180 MeV的大规模格点上实现了2倍的加速。

ABSTRACT

We have implemented a variant of the force gradient integrator proposed by Kennedy et.al. and are using it in our production 2+1 flavor DWF simulations with pion masses of 180 MeV in (4.5fm)3 volumes. We find modest speed-ups (\sim 20%) from using the force gradient integrator, compared to our previously used Omelyan integrator. On other ensembles, primarily finite temperature 2+1 flavor DWF QCD, we have extensively tuned the Hasenbusch preconditioning masses and achieved speed-ups of 2-3x. Here we have also switched to the force gradient integrator, but this change has not had any impact on the speed. We also report on an improved solver for DWF, which uses Möbius fermions, with a smaller fifth dimension than the original DWF fermions, as an intermediate step in the generation of solutions of the Dirac equation. This approach cuts the number of effective Dirac applications by approximately a factor of 2 when the conjugate gradient iteration count is large.

研究动机与目标

• 通过优化分子动力学积分方案，提升2+1味域墙费米子（DWF）格点QCD模拟的效率。
• 降低HMC算法中占主导时间的DWF狄拉克方程求解的计算成本。
• 利用Möbius费米子与标准DWF之间的数学相似性，加速DWF狄拉克方程的求解。
• 实现并测试一种新的力梯度积分器，以在保持辛几何性和可逆性的同时提升稳定性与步长。
• 通过联合使用FGI与Hasenbusch质量分裂及商作用量，在生产规模模拟中实现显著加速。

提出的方法

• 采用两步更新方式实现对称、辛几何的力梯度积分器（FGI）：首先使用力梯度的小指数对规范场进行扰动，然后使用更新后的场重新计算力。
• 使用泰勒展开近似计算二阶力梯度项，以 $ U' = e^{-rac{ au^2}{24}F^j T_j}U $ 的规范场位移替代，实现高效计算。
• 通过重用第一步的解作为第二步的初始猜测，实现“力梯度预测”，从而降低力评估中狄拉克求解的成本。
• 引入一种Möbius加速的DWF求解器，通过在Möbius费米子作用量中使用更小的第五维尺寸 $L'$ 来近似完整DWF解，从而减少4维威尔逊狄拉克算符的应用次数。
• 采用缺陷校正方法，通过迭代使用Möbius算符求解残差误差，提高Möbius近似的精度。
• 将FGI与Hasenbusch质量分裂及商作用量结合，允许对与力相关的狄拉克方程进行精度较低的求解，从而降低HMC的计算成本。

实验结果

研究问题

• RQ1 在大规模DWF模拟中，力梯度积分器（FGI）是否能相对于Omelian积分器实现可测量的性能提升？
• RQ2 Möbius费米子形式是否能通过减小有效第五维尺寸，实现DWF狄拉克方程更快的求解？
• RQ3 Möbius近似中，何种参数组合可在DWF求解器中实现精度与计算成本的最佳平衡？
• RQ4 FGI与Hasenbusch质量分裂及商作用量的结合是否能在生产级HMC模拟中实现显著加速？
• RQ5 力梯度预测技术如何在不破坏可逆性或辛几何性的情况下降低HMC中力评估的成本？

主要发现

• 在4.5\textrm{fm}^3 格点、π介子质量为180 MeV时，FGI实现了约20%的适度加速，且在更大体积下具有更大的提升潜力。
• 在 $32^3 \times 64$ 格点上，Möbius加速的DWF求解器将4维威尔逊狄拉克算符的应用次数从3.6×10⁵减少至2.0×10⁵，降幅达1.8倍。
• 与直接求解相比，Möbius加速的DWF求解器在墙钟时间上实现了2倍加速，且其求解仅占总时间的42%。
• 缺陷校正方法显著提升了Möbius近似的精度，同时仅带来极小的额外计算开销。
• Hasenbusch质量分裂结合商作用量使与力相关的狄拉克方程求解精度降低，对有限温度系综贡献了2–3倍的加速。
• Möbius参数仅需为每个系综调优一次，且在不同配置和源向量间保持有效，可在生产模拟中广泛复用。