[論文レビュー] 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）シミュレーションにおけるエネルギー保存の著しい向上を図るために、高次シンプレクティック積分法—特に力勾配型およびヘッセ行列に基づく積分法—を提案する。ポisson括弧を用いてシャロウハミルトニアンを導出し、積分法パラメータを最適化することで、ステップサイズ依存性を低減しつつ安定性を維持し、格子場理論シミュレーションにおける計算コストの低減と高い受容率を実現する。

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）におけるダイナミカルフェルミオンシミュレーションの計算コストを低減すること。
• 高次シンプレクティック積分法を用いることで、積分ステップサイズを小さくせずにエネルギー保存を改善すること。
• シャロウハミルトニアンにおける測定可能なポisson括弧項に基づく積分法最適化フレームワークの構築。
• 格子場理論におけるHMCにおける力勾配型およびヘッセ行列に基づく積分法の性能向上の可能性の解明。
• パラメータ最適化により真のハミルトニアンとシャロウハミルトニアンの差を最小化することで、より効率的なHMCの実現。

提案手法

• Baker–Campbell–Hausdorff（BCH）展開を用いて導出された、真のハミルトニアンに近い「シャロウ」ハミルトニアンを厳密に保存するシンプレクティック積分法を用いる。
• ポisson括弧を用いてシャロウハミルトニアンを計算し、高次積分法ではO(δτ⁴)まで補正を加える。
• 作用の2階微分を組み込むことでエネルギー保存を向上させる力勾配型およびヘッセ行列に基づく積分法を導出する。
• シャロウハミルトニアンにおけるO(δτ²)およびO(δτ⁴)項を最小化することで、積分法パラメータ（例：PQPQPスキームにおけるα）を最適化する。
• 測定可能なポisson括弧平均を用いて積分法パラメータを事前最適化し、δH = H̃ − H を最小化する手法を提案する。
• リー群手法を用いて多様体上での古典力学を実装し、複雑な積分法における力項およびポアソン括弧の構造的導出を可能にする。

実験結果

リサーチクエスチョン

• RQ1 高次シンプレクティック積分法は、ステップサイズを小さくせずにエネルギードリフトを低減できるか？
• RQ2 BCH展開におけるポアソン括弧構造は、シャロウハミルトニアンおよびシミュレーションの安定性にどのように影響するか？
• RQ3 力勾配型およびヘッセ行列に基づく積分法は、ダイナミカルフェルミオンシミュレーションにおけるエネルギー保存をどの程度向上できるか？
• RQ4 ポアソン括弧測定に基づく積分法パラメータの事前最適化は、計算コストの低減に寄与するか？
• RQ5 シャロウハミルトニアンにおける高次項は、受容率およびシミュレーション効率にどのような影響を及えるか？

主な発見

• 最適化されたα = (1 - √3)/6を有するPQPQP積分法は、シャロウハミルトニアンにおける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シミュレーションにおける受容率の向上を実現し、性能向上を達成する。