[論文レビュー] Event-driven Monte Carlo algorithm for general potentials

Etienne P. Bernard, Werner Krauth|arXiv (Cornell University)|Nov 29, 2011

Theoretical and Computational Physics参考文献 8被引用数 15

ひとこと要約

この論文は、微正準集団における一般のポテンシャルへイベントチェーン・モンテカルロアルゴリズムを拡張し、拒否なしで非局所的な更新が可能となり、詳細つり合いを破る。この手法はエネルギーの量子化を用い、離散化の精度が向上するにつれて漸近的に一定の性能を示し、局所的モンテカルロ法を上回り、連続極限でも効率的なシミュレーションを可能にする。

ABSTRACT

We extend the event-chain Monte Carlo algorithm from hard-sphere interactions to the micro-canonical ensemble (constant potential energy) for general potentials. This event-driven Monte Carlo algorithm is non-local, rejection-free, and allows for the breaking of detailed balance. The algorithm uses a discretized potential, but its running speed is asymptotically independent of the discretization. We implement the algorithm for the cut-off linear potential, and discuss its possible implementation directly in the continuum limit.

研究の動機と目的

• 硬球系に特化した元来のイベントチェーン・モンテカルロアルゴリズムを、微正準集団における一般の二体ポテンシャルへ一般化すること。
• 詳細つり合いを破る非局所的で拒否なしのアルゴリズムを構築し、局所的モンテカルロ法よりも高速なサンプリングを可能にすること。
• エネルギーの離散化ステップ ∆E が 0 に近づくにつれて性能が低下しないように保証すること。これにより、イベント駆動分子動力学の制限を克服すること。
• 粒子の軌道に沿ったエネルギープロファイルを分析することで、連続極限における直接的実装を可能にすること。

提案手法

• ポテンシャルエネルギーを ∆E のステップで離散化し、V(r) と全エネルギー E がすべて ∆E の整数倍となるようにすることで、離散的なイベント追跡を可能にする。
• ランダムに選ばれた粒子を、全移動距離 ℓ に達するまでまたはポテンシャルエネルギーの増加が E を超えると予想されるまで直線的に移動させる。
• エネルギー E を超えると予想される場合、チェーン反応が発生する：経路に続く次の粒子が残りの距離だけ移動し、全 ℓ が消費されるまで繰り返される。
• 軌道に沿ったエネルギー変化を、離散化されたポテンシャルにおける不連続性を特定することで追跡し、ソート済みの交差点を用いて E(x) を段階的に計算する。
• 連続極限では、Vcont(r) = max(0, 1−r) (r≤1) としてポテンシャルをモデル化し、E(x) が連続的かつ区分的 C∞ で、導関数が単調に減少するようにする。
• C∞ 区間上の意思決定木を用いて、Econt(xroot) = E を満たす唯一の xroot を特定し、∆E に依存しない正確なイベント検出を可能にする。

実験結果

リサーチクエスチョン

• RQ1 イベントチェーン・アルゴリズムは、硬球ポテンシャルに限らず、微正準集団における任意の二体ポテンシャルへ一般化可能か？
• RQ2 エネルギーの離散化 ∆E → 0 の際にも、イベント駆動分子動力学とは異なり性能が一定を保つのか？
• RQ3 離散的エネルギーステップに依存せずに、連続極限において直接実装可能か？
• RQ4 詳細つり合いを破りながらも、微小反転性を保ち、正しいサンプリングが可能か？

主な発見

• イベント駆動モンテカルロアルゴリズムは、詳細つり合いが破られても拒否なしで微小反転可能であり、微正準集団における効率的なサンプリングを可能にする。
• 3 GHz のワークステーションで単純な実装でも、1時間あたり約 5×10⁹ 個の衝突を達成し、局所的モンテカルロ法を著しく上回る。
• 実行時間は ∆E に漸近的に依存せず、任意に小さな離散化ステップに対しても効率的なシミュレーションが可能である。
• N=1282 個の粒子に対して、イベント駆動モンテカルロで得られた対相関関数は、局所的モンテカルロと完全に一致し、正しさが確認された。
• 連続極限において、エネルギープロファイル E(x) は連続的かつ区分的 C∞ であり、導関数は単調に減少するため、意思決定木による根の特定が安定に可能である。
• 衝突イベントが Econt(x) = E の一意な解として扱えることから、離散的ステップに依存せず、連続極限への直接実装が可能である。