[论文解读] Speeding up HMC with better integrators

M. A. Clark, A.D. Kennedy|ArXiv.org|Oct 18, 2007

Numerical methods for differential equations参考文献 7被引用 10

一句话总结

本文提出使用高阶辛积分器——特别是力梯度和海森矩阵基积分器——显著提升哈密顿蒙特卡洛（HMC）模拟中动力费米子计算的能量守恒。通过利用泊松括号推导影子哈密顿量并优化积分器参数，该方法在不减小步长的情况下降低了对步长的依赖性，同时保持了稳定性，从而在格点场论模拟中实现了更低的计算成本和更高的接受率。

ABSTRACT

We discuss how dynamical fermion computations may be made yet cheaper by using symplectic integrators that conserve energy much more accurately without decreasing the integration step size. We first explain why symplectic integrators exactly conserve a ``shadow'' Hamiltonian close to the desired one, and how this Hamiltonian may be computed in terms of Poisson brackets. We then discuss how classical mechanics may be implemented on Lie groups and derive the form of the Poisson brackets and force terms for some interesting integrators such as those making use of second derivatives of the action (Hessian or force gradient integrators). We hope that these will be seen to greatly improve energy conservation for only a small additional cost and that their use will significantly reduce the cost of dynamical fermion computations.

研究动机与目标

• 降低格子量子色动力学（QCD）中动力费米子模拟的计算成本。
• 通过使用高阶辛积分器，在不减小积分步长的情况下改善HMC中的能量守恒。
• 基于影子哈密顿量中的可观测量泊松括号项，开发积分器优化框架。
• 探索力梯度和海森矩阵基积分器在格场理论HMC中的性能提升。
• 通过参数优化最小化真实哈密顿量与影子哈密顿量之间的差异，实现更高效的HMC。

提出的方法

• 使用精确守恒一个接近真实哈密顿量的‘影子’哈密顿量的辛积分器，该影子哈密顿量通过Baker–Campbell–Hausdorff（BCH）展开推导得出。
• 应用泊松括号计算影子哈密顿量，对高阶积分器的修正项可达O(δτ⁴)。
• 推导出包含作用量二阶导数的力梯度和海森矩阵基积分器，以改善能量守恒。
• 通过最小化影子哈密顿量中的O(δτ²)和O(δτ⁴)项，优化积分器参数（如PQPQP方案中的α）。
• 提出基于可观测量泊松括号平均值的离线积分器参数优化方法，以最小化δH = H̃ − H。
• 采用李群方法在流形上实现经典力学，从而系统化推导复杂积分器的力项和泊松括号。

实验结果

研究问题

• RQ1 高阶辛积分器是否能在不减小步长的情况下减少HMC中的能量漂移？
• RQ2 BCH展开中的泊松括号结构如何影响影子哈密顿量和模拟稳定性？
• RQ3 力梯度和海森矩阵基积分器在动力费米子模拟中能在多大程度上改善能量守恒？
• RQ4 基于泊松括号测量的积分器参数离线优化是否能降低计算成本？
• RQ5 影子哈密顿量中的高阶项对接受率和模拟效率有何影响？

主要发现

• 经过优化的PQPQP积分器（α = (1 - √3)/6）可最小化影子哈密顿量中的O(δτ²)项，从而减少能量漂移。
• 带有O(δτ⁴)修正项的力梯度积分器在几乎不增加额外成本的情况下，展现出显著改善能量守恒的潜力。
• 基于泊松括号推导系数的海森矩阵基积分器实现了O(δτ⁴)阶的影子哈密顿量修正，与标准积分器相比，使δH降低数个数量级。
• 力梯度积分器的影子哈密顿量包含复杂的嵌套泊松括号组合，其系数由BCH展开至O(δτ⁴)推导得出。
• 形式为e^{(1/6)δτ T}e^{(3/8)δτ S}e^{(1/3)δτ T}...的积分器，其影子哈密顿量修正为δH ∝ δτ⁴ × (2259{S,{S,{S,{S,T}}}} + ...)，系数由BCH展开得出。
• 该方法通过直接测量泊松括号平均值来优化积分器参数，从而减小δH，并在HMC模拟中提升接受率。