[論文レビュー] Canonical sampling through velocity-rescaling

Giovanni Bussi, Davide Donadio+1 more|arXiv (Cornell University)|Mar 28, 2008

Spectroscopy and Quantum Chemical Studies被引用数 4

ひとこと要約

本稿では、ガウス分布に従う確率変数から導かれるランダム要因を用いて全粒子の速度を確率的・スケーリングすることで、正準集団の正確なサンプリングを保証する、新しい分子動力学の温度制御法を提案する。この手法は、誤差の監視に役立つ保存量「有効エネルギー」を維持し、Lennard-JonesおよびTIP4P水モデルにおいて優れた性能を示し、動的性質が温度制御の緩和時間にほとんど依存しないことを示している。

ABSTRACT

We present a new molecular dynamics algorithm for sampling the canonical distribution. In this approach the velocities of all the particles are rescaled by a properly chosen random factor. The algorithm is formally justified and it is shown that, in spite of its stochastic nature, a quantity can still be defined that remains constant during the evolution. In numerical applications this quantity can be used to measure the accuracy of the sampling. We illustrate the properties of this new method on Lennard-Jones and TIP4P water models in the solid and liquid phases. Its performance is excellent and largely independent on the thermostat parameter also with regard to the dynamic properties.

研究の動機と目的

• 正準集団を正確かつ効率的にサンプリングできる、シンプルで堅牢な分子動力学用温度制御法の開発。
• Berendsen や Langevin などの人気のある確率的温度制御法に保存量がないという問題を解決し、誤差監視やコード検証を可能にする。
• 適切に選ばれたランダム要因を組み込んだ速度スケーリングの一般化により、正しい統計的サンプリングを保ちつつ数値的安定性を維持する。
• 正準集団違反や時間刻み誤差の信頼性のある評価を可能にする、新たな保存量「有効エネルギー」を導入する。
• Lennard-JonesおよびTIP4P水を含む実際の系（固体および液体相）において、本手法の有効性と遍歴性を実証する。

提案手法

• アルゴリズムは、正準集団の目標温度に対応する運動エネルギーを保証するように、すべての粒子の速度を確率的要因 α でスケーリングする。
• スケーリング要因 α は、緩和時間 τ がその強度を制御するオーナイズ型確率過程に類似したプロセスから導かれる。
• 新しい保存量「有効エネルギー」が動的過程で一定のままであることを示すことで、本手法の形式的妥当性を裏付ける。これは、NVEシミュレーションにおけるエネルギー保存に類似している。
• 有効エネルギーは運動エネルギーと確率的変数の関数として定義され、集団違反や時間刻み誤差の直接的監視が可能になる。
• Berendsen の温度制御法を一般化し、正しくサンプリングできるようにランダム成分を追加しながら、安定性とシンプルさを維持する。
• α² の閉形式表現を用いてアルゴリズムを実装し、ガウス分布に従う確率変数とガンマ分布の和を含むことで、効率的な数値評価が可能になる。

実験結果

リサーチクエスチョン

• RQ1 確率的スケーリングを組み込んだ速度スケーリング法は、正準集団を正しくサンプリングできるように一般化できるか？
• RQ2 確率的スケーリング要因の導入により、誤差監視が可能な保存量が得られるか？
• RQ3 本手法の性能は、Nosé-Hoover や Berendsen といった既存手法と比較して、遍歴性および動的性質の観点で優れているか？
• RQ4 自己拡散係数 D などの系の動的挙動は、温度制御の緩和時間 τ にどの程度依存するか？
• RQ5 本手法で導入された「有効エネルギー」は、有限時間刻み誤差を信頼性高く検出でき、NVT集団の正確なサンプリングを保証できるか？

主な発見

• 提示された温度制御法は、遍歴性の仮定のもとで、正準集団を正しくサンプリングすることが解析的に証明された。
• 運動エネルギーと確率的変数の関数として定義された新たな保存量「有効エネルギー」が同定され、シミュレーション中も一定に保たれることを示した。これにより、誤差監視や時間刻みの選定が可能になった。
• Lennard-JonesおよびTIP4P水モデルの固体および液体相において、本手法は優れた性能を示し、自己拡散係数 D などの動的性質が温度制御の緩和時間 τ にほとんど依存しないことが確認された。
• 300 K における TIP4P 水の自己拡散係数 D は、文献値および NVE シミュレーションの結果と一致しており、τ = 0.002 ps 時のずれは僅かにしか生じなかった。
• 数値的テストでは、Nosé-Hoover 温度制御法よりも本手法の方が遍歴性に優れており、位相空間のより良いサンプリングが得られた。
• アルゴリズムは実装が簡単であり、Berendsen、Nosé、および本研究の温度制御法を特別な場合として含む統一的なフレームワークを提供する。