Verlet積分

Verletアルゴリズムとその速度形式は、分子動力学における標準的な積分法であり、時間反転可能でシンプレクティックであり、シミュレーションに必要な数百万ステップにわたってエネルギーを良好に保存するため、高く評価されています。

Definition

Verlet積分は、ニュートン運動方程式を積分するための時間反転可能でシンプレクティックな方法であり、現在の位置と以前の位置、および加速度を使用して粒子の位置を更新し、分子動力学の安定した軌道をもたらします。

Scope

このトピックでは、Verlet積分法のファミリー、すなわち、元の位置Verletスキーム、同等のリープフロッグおよび速度Verlet定式化、それらの時間反転可能性とシンプレクティック構造、および結果として生じる長期的なエネルギー保存について扱います。これらの方法は、ハミルトン系のシンプレクティック積分に関するより広範な理論の中に位置づけられます。

Core questions

• Verletスキームは、力から位置と速度をどのように進めるのでしょうか？
• Verletアルゴリズムはなぜ時間反転可能でシンプレクティックなのでしょうか？
• Verlet積分は、非常に長いシミュレーションにわたってエネルギーをなぜ良好に保存するのでしょうか？
• 位置Verlet、リープフロッグ、および速度Verletの定式化は、互いにどのように関連しているのでしょうか？

Key theories

シンプレクティック構造と時間反転可能性

Verlet積分は、相空間のシンプレクティック幾何学を保存し、時間反転に対して不変であるため、保存系の非シンプレクティック積分器を悩ませる系統的なエネルギードリフトを防ぎます。

シャドウハミルトニアンの保存

離散的なVerlet軌道は真のエネルギーを厳密には保存しませんが、密接に関連するシャドウハミルトニアンをほぼ保存し、エネルギー誤差が成長するのではなく、有界で振動するように保たれます。

同等の定式化

位置Verlet、リープフロッグ、および速度Verletスキームは同じ軌道を生成しますが、速度が利用可能になる方法とタイミングが異なります。同期された位置と速度が必要な場合は、速度Verletが好まれます。

Clinical relevance

Verlet積分は、単純なレナード・ジョーンズ流体から大規模な生体分子シミュレーションまで、実質的にすべての分子動力学コードにおけるデフォルトの時間刻みエンジンであり、同じシンプレクティック原理は天文学における長期軌道積分にも使用されています。

History

このスキームは、20世紀初頭に天文学者カール・ストーマーによって使用され、1967年のレナード・ジョーンズ流体の研究でルー・ヴェルレによって分子シミュレーションに普及しました。その後の分析により、それがシンプレクティック積分器であることが示され、その優れた長期安定性が説明されました。

Key figures

• Loup Verlet
• Carl Stormer
• Ernst Hairer

Seminal works

• verlet1967
• hairer1993

Frequently asked questions

分子動力学において、Verletがより高次のルンゲ＝クッタ法よりも好まれるのはなぜですか？

ルンゲ＝クッタは1ステップあたりの精度が高い場合がありますが、シンプレクティックではなく、長時間の実行ではエネルギーがゆっくりとドリフトします。Verletのシンプレクティックで時間反転可能な構造は、数百万ステップにわたってエネルギーを有界に保ちます。これは、平衡シミュレーションにとって、1ステップあたりの精度よりもはるかに重要です。

Verletはエネルギーを厳密に保存しますか？

いいえ。厳密なエネルギーではなく、近くのシャドウハミルトニアンを保存するため、測定されたエネルギーはドリフトするのではなく、有界な範囲内で振動します。これは、安定した熱力学的平均を計算するには十分です。