[论文解读] Implementation of the HMC algorithm on the tempered Lefschetz thimble method

Masafumi Fukuma, Nobuyuki Matsumoto+1 more|arXiv (Cornell University)|Dec 31, 2019

Theoretical and Computational Physics参考文献 37被引用 24

一句话总结

该论文在热化Lefschetz指环方法（TLTM）上实现了混合蒙特卡洛（HMC）算法，以解决格点场论中，特别是费米子系统中的数值符号问题。通过将RATTLE算法适配于流动表面的分子动力学，并开发出对费米子行列式零点的稳健处理方法，作者在2×2 Hubbard模型测试中，将计算成本降低了70%，且自相关时间显著低于Metropolis算法。

ABSTRACT

The tempered Lefschetz thimble method (TLTM) is a parallel-tempering algorithm towards solving the numerical sign problem, where the system is tempered by the antiholomorphic gradient flow to tame both the sign and ergodicity problems simultaneously. In this paper, we implement the hybrid Monte Carlo (HMC) algorithm for transitions on each flowed surface, expecting that this implementation on TLTM will give a useful framework for future computations of large-scale systems including fermions. Although the use of HMC in Lefschetz thimble methods has been proposed so far, our crucial achievement here is that HMC is implemented on TLTM so as to work within the parallel-tempering algorithm in TLTM, especially by developing an algorithm to handle zeros of fermion determinants in the course of the molecular-dynamics process. We confirm that the algorithm works correctly by applying it to the sign problem of the Hubbard model on a small lattice, for which the TLTM is known to work with the Metropolis algorithm. We show that the use of HMC significantly reduces the autocorrelation times with less computational times compared to the Metropolis algorithm.

研究动机与目标

• 为克服Metropolis算法在TLTM中因高自相关性和混合缓慢而带来的局限性，尤其在费米子系统中。
• 在TLTM上实现HMC，以提升涉及费米子的大规模模拟的效率。
• 开发一种处理在流动指环表面进行分子动力学时出现的费米子行列式零点的方法。
• 使用小尺寸格点的Hubbard模型作为基准，验证HMC在TLTM上的实现。
• 证明与Metropolis算法相比，HMC在TLTM上能降低计算成本和自相关时间。

提出的方法

• 将RATTLE算法适配于在流动表面Σₜ上执行分子动力学，通过速度-Verlet积分方法确保约束的保持。
• 引入一种带自适应步长控制和缩放的牛顿-拉夫森迭代方法，以求解Σₜ上的约束动力学方程。
• 采用一种判据来检测并避免进入Dₜ区域（即流向费米子行列式零点的区域），以接近Dₜ作为停止条件。
• 当牛顿迭代失败或进入Dₜ时，使用动量翻转（Ψ）作为备用方案，以在HMC过程中保持细致平衡。
• 将HMC动力学整合进TLTM框架，实现在不同流动时间t的平行热化。
• 将该算法应用于Nₛ=2×2和Nₜ=5的Hubbard模型，与精确解及基于Metropolis的TLTM结果进行比较。

实验结果

研究问题

• RQ1 HMC能否在TLTM框架中成功实现，以提升费米子系统采样效率？
• RQ2 在流动指环表面进行分子动力学时，如何处理费米子行列式零点而不破坏可逆性？
• RQ3 与Metropolis算法相比，HMC在TLTM上是否能降低符号问题严重系统中的自相关时间与计算成本？
• RQ4 在远离半满填充时，HMC扩展的TLTM方法在Hubbard模型上是否保持稳定与准确？
• RQ5 在接近费米子行列式零点的区域，该算法能否在保持细致平衡与遍历性的同时运行？

主要发现

• 在2×2 Hubbard模型中，HMC在TLTM上的实现将每组独立配置的计算成本降低至Metropolis算法的约30%。
• HMC显著降低了自相关时间，表明混合速度更快，采样效率更高。
• 该算法正确复现了Hubbard模型中数密度的精确结果，证实了其准确性。
• 通过迭代牛顿求解结合缩放与区域检测处理费米子行列式零点，确保了方法的鲁棒性与稳定性。
• 当解发散或进入禁区时，通过使用动量翻转，方法保持了细致平衡与可逆性。
• 该框架具有可扩展性，预计在自由度更高的大系统中将获得更大的效率提升。