[論文レビュー] Common Molecular Dynamics Algorithms Revisited

Alexey K. Mazur|arXiv (Cornell University)|Sep 15, 1997

Scientific Research and Discoveries参考文献 23被引用数 31

ひとこと要約

この論文は分子動力学におけるルンゲ・クッタ積分法の再評価を行い、そのエネルギー保存性および切り捨て誤差の性質が従来の認識を上回ることを示している。標準的なエネルギー揺らぎ測定における補間誤差が見かけの不安定さの原因であると特定し、そのような誤差を含まない修正手法を提案する。これにより、現在のステップサイズが正確なサンプリングにあまりに小さく設定されている可能性が示唆される。

ABSTRACT

The Stormer?Verlet?leapfrog group of integrators commonly used in molecular dynamics simulations has long become a textbook subject and seems to have been studied exhaustively. There are, however, a few striking effects in performance of algorithms which are well known but have not received adequate attention in the literature. A closer view of these unclear observations results in unexpected conclusions. It is shown here that contrary to the conventional point of view, the leapfrog scheme is distinguished in this group both in terms of the order of truncation errors and the conservation of the total energy. In this case the characteristic square growth of fluctuations of the total energy with the step size, commonly measured in numerical tests, results from additional interpolation errors with no relation to the accuracy of the computed trajectory. An alternative procedure is described for checking energy conservation of leapfrog-like algorithms which is free from interpolation errors. Preliminary tests on a representative model system suggest that standard step size values used at present are lower than necessary for accurate sampling.

研究の動機と目的

• 分子動力学シミュレーションにおけるルンゲ・クッタ積分法の性能を再表現・再評価すること。
• 標準的な数値的テストにおける誤ったエネルギー揺らぎ測定の原因を特定・解消すること。
• 補間誤差を回避するエネルギー保存性の評価手法を提案すること。
• 現在のステップサイズ選択に関する一般的な認識に疑問を呈し、正確なサンプリングに必要以上に小さなステップサイズが使われている可能性を示すこと。

提案手法

• 理論的および数値的分析を用いて、ルンゲ・クッタ積分法の切り捨て誤差およびエネルギー保存性の性質を分析する。
• 見かけのエネルギー揺らぎがステップサイズに比例して二次的に増大するのは、本質的な不正確さではなく、補間誤差に起因することを特定する。
• 一貫した時刻点でのサンプリングにより補間を回避する代替エネルギー評価手順を提案する。
• 代表的なモデル系を用いて新しい手法をテストし、標準的手法と比較してエネルギー保存性を評価する。
• 数値実験を用いて、異なるステップサイズにおけるエネルギー揺らぎおよび軌道の正確性を比較する。

実験結果

リサーチクエスチョン

• RQ1 なぜ標準的な数値的テストにおいて、ルンゲ・クッタ積分法のエネルギー揺らぎがステップサイズに比例して二次的に増大するのか？
• RQ2 理論的エネルギー保存性と、ルンゲ・クッタシミュレーションでの観測揺らぎの間に生じる乖離の原因は何か？
• RQ3 補間誤差を導入せずにエネルギー保存性を正確に測定するにはどうすればよいか？
• RQ4 現在の分子動力学シミュレーションにおける標準的なステップサイズは、正確なサンプリングに必要以上に小さいのだろうか？
• RQ5 ルンゲ・クッタ型アルゴリズムにおいて、ステップサイズとエネルギー保存性の真の関係は何か？

主な発見

• ルンゲ・クッタ積分法は、ヴェルレ族の他のメンバーと比較して、優れた切り捨て誤差特性およびエネルギー保存性を示す。
• エネルギー揺らぎがステップサイズに比例して二次的に増大するという観測は、標準的な測定法における補間誤差に起因する偽の現象であり、精度の低さの兆候ではない。
• 提案された代替エネルギー評価手法は補間誤差を排除し、エネルギー保存性のより正確な評価を可能にする。
• 予備的なテストでは、現在の分子動力学における標準的なステップサイズが、正確なサンプリングに必要な最小値を下回っている可能性がある。
• エネルギー保存性が正しく測定された場合、ルンゲ・クッタ法の性能は従来の認識を上回っていることが示され、特に正確な測定が行われた場合には顕著である。