[论文解读] Tuning HMC using Poisson brackets

M. A. Clark, A.D. Kennedy+1 more|arXiv (Cornell University)|Oct 7, 2008

Markov Chains and Monte Carlo Methods被引用 12

一句话总结

本文提出将真实哈密顿量与影子哈密顿量之间差异ΔH的方差最小化，作为优化混合蒙特卡洛（HMC）积分器的更优策略，而非直接最小化|ΔH|本身。通过将ΔH表示为泊松括号的级数展开，并测量其平衡态平均值，可仅通过一次模拟即确定最优积分器参数（如步长和STSTS方案中的系数），显著提升格点场论中HMC的效率。

ABSTRACT

We discuss how the integrators used for the Hybrid Monte Carlo (HMC) algorithm not only approximately conserve some Hamiltonian $H$ but exactly conserve a nearby shadow Hamiltonian ( ilde H), and how the difference $ΔH \equiv ilde H - H $ may be expressed as an expansion in Poisson brackets. By measuring average values of these Poisson brackets over the equilibrium distribution $\propto e^{-H}$ generated by HMC we can find the optimal integrator parameters from a single simulation. We show that a good way of doing this in practice is to minimize the variance of $ΔH$ rather than its magnitude, as has been previously suggested. Some details of how to compute Poisson brackets for gauge and fermion fields, and for nested and force gradient integrators are also presented.

研究动机与目标

• 通过优化积分器参数，提升格格点场论中混合蒙特卡洛（HMC）模拟的效率。
• 解决在动力费米子模拟中为实现高接受率和低自相关性而调优HMC积分器的挑战。
• 开发一种仅需单次模拟运行即可确定最优积分器参数的实用方法。
• 证明最小化ΔH方差比最小化|ΔH|在积分器调优中更为有效。
• 为规范理论及含动力费米子系统的泊松括号计算提供一个框架。

提出的方法

• 将真实哈密顿量与影子哈密顿量之间的差异ΔH = H̃ - H 表示为动能与势能项泊松括号的级数展开。
• 利用Baker–Campbell–Hausdorff（BCH）公式，推导出对称辛积分器（如STS和STSTS）的影子哈密顿量H̃。
• 通过单条已平衡的HMC轨迹，测量泊松括号项⟨{S,{S,T}}⟩和⟨{T,{S,T}}⟩的平衡态平均值，以预测积分器性能。
• 通过在平衡分布exp(−H)上最小化ΔH的方差，优化积分器参数（如STSTS中的α）。
• 将该方法应用于嵌套积分器和力梯度积分器，并在格点规范理论中通过数值结果加以验证。
• 利用单组参数设置下测得的泊松括号，预测参数空间中ΔH及其方差的行为。

实验结果

研究问题

• RQ1 在平衡分布上最小化ΔH的方差，是否能比最小化|ΔH|带来更优的HMC积分器调优效果？
• RQ2 如何高效计算并测量规范系统与费米子系统中的泊松括号，以实现积分器优化？
• RQ3 通过测量泊松括号平均值，能否仅用一次模拟运行即实现积分器参数的充分优化？
• RQ4 当采用ΔH方差最小化策略调优时，力梯度积分器与标准Omelian积分器的性能表现有何差异？
• RQ5 所提出的方法能否推广至格点场论中的嵌套积分器与多尺度模拟？

主要发现

• 最小化ΔH的方差可获得比最小化|ΔH|更优的HMC性能，因为ΔH的波动远小于其均值。
• 通过最小化Var(ΔH)确定的STSTS积分器最优α值约为α ≈ 0.1826，与其它性能指标（如⟨|δH|⟩和⟨1−min(1,e−δH)⟩）的最小值高度一致。
• ΔH方差最小值的实测位置与预测位置高度吻合，验证了仅在单一α设置下测量泊松括号即可准确预测最优参数的可行性。
• 当采用方差最小化策略调优时，力梯度积分器相比二阶Omelian积分器展现出更优的标度行为。
• 该方法可实现积分器参数的离线优化，包括步长、阶数及赝费米子场数量，且仅需一次已平衡的HMC运行。
• 规范场与费米子场的泊松括号可系统性地计算并用于表达ΔH，从而在真实的格点QCD模拟中实现实用化的参数调优。