[论文解读] Common Molecular Dynamics Algorithms Revisited

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

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

一句话总结

本文重新評估分子動力學中的Leapfrog積分器，證明其能量守恆與截斷誤差特性優於傳統認知。研究識別出標準能量波動測量中的插值誤差為表觀不穩定的根源，並提出一種無此誤差的修正方法，建議目前的步長可能過於細小，不足以實現精確取樣。

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.

研究动机与目标

• 重新表達並重新評估分子動力學模擬中Leapfrog積分器的性能。
• 識別並解決標準數值測試中誤導性能量波動測量的來源。
• 提出一種消除插值誤差的新型能量守恆評估方法。
• 透過證明目前的步長值可能過於細小，挑戰分子動力學中關於步長選擇的傳統觀點。

提出的方法

• 使用理論與數值分析，分析Leapfrog積分器的截斷誤差與能量守恆特性。
• 識別出表觀能量波動隨步長呈二次增長，是因插值誤差所致，而非本質上的不準確。
• 提出一種替代的能量評估程序，透過使用一致的時間點取樣來避免插值。
• 使用代表性模型系統測試新方法，並比較標準與修正程序下的能量守恆表現。
• 透過數值實驗，比較不同步長下能量波動與軌跡準確性的差異。

实验结果

研究问题

• RQ1 為何在標準數值測試中，Leapfrog積分器的能量波動會隨步長呈二次增長？
• RQ2 理論上的能量守恆與Leapfrog模擬中觀察到的波動之間的差異是由何原因造成？
• RQ3 如何在不引入插值誤差的情況下，精確測量能量守恆？
• RQ4 目前分子動力學模擬中的標準步長是否過於細小，導致無法實現精確取樣？
• RQ5 Leapfrog類演算法中，步長與能量守恆之間的真實關係為何？

主要发现

• Leapfrog積分器在截斷誤差特性與能量守恆方面，優於Verlet家族的其他成員。
• 觀察到的能量波動隨步長呈二次增長，是標準測量方法中插值誤差的產物，並非準確度不佳的徵兆。
• 所提出的替代能量評估方法消除了插值誤差，提供了更精確的能量守恆評估。
• 初步測試顯示，目前分子動力學中使用的標準步長低於精確取樣所需的水準。
• 當能量守恆以正確方式測量時，Leapfrog方法的表現優於傳統認知，尤其在高準確度要求下。