[论文解读] Comparison of the mass preconditioned HMC and the DD-HMC algorithm for two-flavour QCD

Marina Krstić Marinković, Stefan Schaefer|arXiv (Cornell University)|Nov 3, 2010

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

一句话总结

本论文通过使用共享的局部降阶求解器，对两味格点QCD中的质量预条件HMC（MP-HMC）与DD-HMC进行了公平的性能评估比较。在mπ ≈ 420 MeV的48×24³格点上，两种算法表现出相近的效率，具有相似的积分自相关时间与接受率，表明采用降阶的MP-HMC是真实模拟中DD-HMC的可行替代方案。

ABSTRACT

Mass preconditioned HMC and DD-HMC are among the most popular algorithms to simulate Wilson fermions. We present a comparison of the performance of the two algorithms for realistic quark masses and lattice sizes. In particular, we use the locally deflated solver of the DD-HMC environment also for the mass preconditioned simulations.

研究动机与目标

• 评估并比较质量预条件HMC（MP-HMC）与DD-HMC在真实夸克质量下的两味QCD性能。
• 评估在两种算法中使用同一套本地降阶求解器对实现公平、与实现无关的比较的影响。
• 研究在轨迹成本与自相关时间方面，采用降阶的MP-HMC是否能与DD-HMC达到相当的效率。
• 探讨MP-HMC在参数调优与可扩展性方面的潜在优势，例如未来扩展至非 degenerate 夸克或SF边界条件。

提出的方法

• 研究采用威尔逊规范作用项与非微扰O(a)-改进的威尔逊费米子，κ_sea = 0.13625，在β = 5.3下于48×24³格点上得到mπ ≈ 420 MeV。
• MP-HMC通过对称的偶奇预条件实现，其中舒尔补项通过正质量Δm ≈ 0.09进行平移，以预条件费米子行列式。
• DD-HMC采用块大小为6²×12²的区域分解，导致R = 36%的活跃链接，并使用相同的本地降阶GCR求解器于两种算法。
• 力被划分为块力（F₁）与相互作用力（F₂），其中F₂在更大的时间尺度上积分，以提高效率。
• 通过舒尔补分解构建有效作用量，实现对低动量与高动量模式的独立处理。
• 性能通过每条轨迹的时钟时间、接受率、积分自相关时间（τ_int）及平面子期望值进行评估。

实验结果

研究问题

• RQ1 在真实夸克质量下，采用降阶求解器的MP-HMC是否在两味QCD模拟中实现与DD-HMC相当的性能？
• RQ2 在相同的求解器与模拟条件下，MP-HMC与DD-HMC的积分自相关时间与接受率如何比较？
• RQ3 MP-HMC中预条件质量的连续调节是否相比DD-HMC中离散的块大小选择具有优势？
• RQ4 使用共享的降阶求解器在多大程度上消除了算法比较中的实现偏差？

主要发现

• MP-HMC与DD-HMC算法表现出相近的性能，相同可观测量统计量下，其每条轨迹的时钟时间分别为1530秒与2010秒。
• 将平面子的积分自相关时间按活跃链接比例（R）归一化后，DD-HMC为10(5)，MP-HMC为10(4)，表明两者的去相关效率相当。
• DD-HMC的接受率为90(1)%，MP-HMC为85(2)%，表明尽管力结构不同，两种算法均保持了高接受率。
• 使用共享的本地降阶求解器显著加速了两种算法，缩小了性能差距，实现了公平比较。
• 能量偏差（ΔH的尖峰）与F₂力的不规则性相关，而MP-HMC由于预条件质量的连续调节，能更有效地控制此类不规则性。
• MP-HMC展现出更大的可扩展性潜力，尤其在引入非简并或重夸克方面，且更易于扩展代码功能，如支持Schrödinger函数边界条件。