分子動力学シミュレーション

分子動力学は、相互作用する原子系のニュートン運動方程式を統合し、構造的、動的、熱力学的特性が計算される軌跡を生成する。

Definition

ポテンシャルエネルギー関数からの力の下で原子の位置と速度を時間的に伝播させ、位相空間を通る決定論的な軌跡を生成するシミュレーション技術。

Scope

古典的な運動方程式の数値積分、Verletファミリーのような時間可逆積分器、サーモスタットとバロスタットによる温度と圧力の制御、長距離静電相互作用と周期的境界の処理、および軌跡からの平衡特性と輸送特性の抽出を扱う。

Core questions

• ニュートンの運動方程式は、多数のタイムステップにわたってどのように安定して統合されるのか？
• 選択された熱力学的アンサンブルをサンプリングするために、温度と圧力はどのように制御されるのか？
• 長距離静電相互作用はどのように効率的に処理されるのか？
• 有限の軌跡から観測可能な特性はどのように回復されるのか？

Key theories

Verlet積分

運動方程式を積分するためのシンプルで時間可逆的かつシンプレクティックなスキームであり、長時間のシミュレーションでエネルギーを良好に保存し、ほとんどの分子動力学コードの基礎となっている。

エルゴードサンプリング

エルゴード仮説の下では、十分に長い軌跡に沿った時間平均はアンサンブル平均に等しく、シミュレーションを統計力学的観測量に結びつける。

Mechanisms

各ステップでは、ポテンシャル勾配から力が計算され、積分器によって位置と速度が進められ、サーモスタットまたはバロスタットの補正が適用される。これを繰り返すことで軌跡が構築され、その平均値から熱力学的および動的特性が得られる。

Clinical relevance

分子動力学は、タンパク質、膜、材料におけるコンフォメーション変化、結合イベント、拡散、および反応に隣接するプロセスを明らかにし、生物物理学および創薬における実験を補完するメカニズム的洞察を提供する。

History

AlderとWainwrightの剛体球シミュレーション、およびRahmanの液体アルゴンの連続ポテンシャル研究に続き、Verletの1967年のアルゴリズムと、その後のサーモスタット、バロスタット、およびEwaldベースの静電相互作用の開発により、分子動力学は成熟したシミュレーション分野として確立された。

Key figures

• Loup Verlet
• Berni Alder
• Aneesur Rahman
• Herman Berendsen

Seminal works

• verlet1967
• frenkel2002

Frequently asked questions

分子動力学におけるタイムステップを制限するものは何か？

最も速い運動、典型的には水素を含む結合振動が限界を設定する。これらの結合を拘束することで、数フェムト秒のタイムステップが可能になり、到達可能な総時間はマイクロ秒以上に達する。

分子動力学は決定論的か？

初期条件が与えられれば積分は決定論的であるが、軌跡はカオス的であるため、ごくわずかな違いが急速に増大する。意味のある結果は個々の経路からではなく、統計的平均から得られる。