[論文レビュー] Determination of Forces from a Potential in Molecular Dynamics

Bernard Monasse, Frédéric Boussinot|arXiv (Cornell University)|Jan 6, 2014

Origins and Evolution of Life参考文献 4被引用数 5

ひとこと要約

この論文は、一般的な分子力学ポテンシャル—結合、Valence角、二面角、Lennard-Jones相互作用—から原子力のきめ細やかなベクトル形式の導出を提示する。力とトルクのバランスを保つことで、系のエネルギーを維持する。導出された力がニュートンの第三法則を満たし、力学的平衡を保存することを証明しており、人工的なエネルギーのずれが生じない安定な分子動力学的シミュレーションを可能にする。

ABSTRACT

In Molecular Dynamics (MD), the forces applied to atoms derive from potentials which describe the energy of bonds, valence angles, torsion angles, and Lennard-Jones interactions of which molecules are made. These de finitions are classic; on the contrary, their implementation in a MD system which respects local equilibrium of mechanical conditions is usually not described. The precise derivation of the forces from the potential and the proof that their application preserves energy is the object of this note. This work is part of the building of a multi-scale MD system, presently under development.

研究の動機と目的

• 分子動力学（MD）における標準的な分子力学ポテンシャルから、数学的に正確な力の導出を提供すること。
• 導出された力がニュートンの第三法則を満たし、力学的平衡を保存することを保証すること。
• これらの力の適用が、MD シミュレーション中に人工的なエネルギー変化を引き起こさないことを証明すること。
• 一貫性のある力場実装フレームワークを確立することで、マルチスケール MD システムの開発を支援すること。

提案手法

• 力はスカラー型ポテンシャルエネルギー関数の負の勾配として導出する：$ \vec{F} = -\nabla \mathcal{U} $。
• ベクトル微積分を用いて、ポテンシャルの原子位置に関する偏微分を計算し、正しい方向成分を保証する。
• 作用・反作用の法則を適用して力のバランスを確保する：結合した原子対に対して $ \vec{f}_a + \vec{f}_b = 0 $。
• 角度および二面角の場合は、ベクトルの射影および正規化された外積を用いて、分子面内の力の方向を定義する。
• Lennard-Jonesポテンシャルの場合は、径方向の力を $ \vec{f}_a = \frac{24\epsilon}{r} \left(2(\sigma/r)^{12} - (\sigma/r)^6 \right) \hat{u} $ として導出する。ここで $ \hat{u} $ は a から b への単位ベクトルである。
• すべての原子における力およびトルクの合計がゼロであることを検証し、力学的平衡およびエネルギー保存を保証する。

実験結果

リサーチクエスチョン

• RQ1 どのようにしてベクトル微積分を用いて、標準的な分子力学ポテンシャルから力 systematically に導出できるか？
• RQ2 導出された力が MD シミュレーションにおける力学的平衡を保つために満たすべき条件は何か？
• RQ3 結合、角度、二面角、Lennard-Jones 相互作用の各タイプについて、力とトルクのバランスを数学的にどのように証明できるか？
• RQ4 マルチスケール MD システムに力の実装を行うために、一貫性があり座標に依存しないベクトル形式を用いることができるか？
• RQ5 これらの力の適用が、シミュレーション中に人工的なエネルギー変化を引き起こすか？

主な発見

• 結合、価角、ねじれ角、Lennard-Jones相互作用のすべてのタイプについて、力はポテンシャルエネルギー関数の負の勾配として導出される。
• 導出された力はニュートンの第三法則を満たす：$ \vec{f}_a + \vec{f}_b = 0 $ であり、結合原子対における力のバランスが保たれる。
• どの相互作用タイプに対しても、力が同一直線上にあるか、または回転効果を相殺するように対称的に分布しているため、合成トルクはゼロとなる。
• 系に含まれるすべての原子における力の合計がゼロであるため、力学的平衡がシミュレーション中に維持されていることが確認された。
• 力の導出プロセスによりエネルギー保存が保証され、力の適用時に人工的なエネルギーが追加されない。
• この手法は座標に依存しないベクトル形式で実装されており、堅牢でポータブルな MD コード開発に適している。