[论文解读] 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

一句话总结

本文提供了从常见的分子力学势能函数（键合、键角、二面角及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系统的开发。

提出的方法

• 将受力表示为标量势能函数的负梯度：$ \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 如何利用矢量微积分系统地从标准分子力学势能中推导受力？
• RQ2 为确保所推导受力在分子动力学模拟中保持机械平衡，需满足哪些条件？
• RQ3 如何对键合、键角、二面角及Lennard-Jones相互作用的受力与力矩平衡进行数学证明？
• RQ4 能否使用一致的、与坐标无关的矢量形式，实现多尺度MD系统中的受力计算？
• RQ5 在模拟过程中应用这些受力是否会引入非物理的能量变化？

主要发现

• 对于所有相互作用类型——键合、键角、二面角及Lennard-Jones势能——受力均被推导为势能函数的负梯度。
• 所推导的受力满足牛顿第三定律：$ \vec{f}_a + \vec{f}_b = 0 $，确保成键对之间的受力平衡。
• 所有相互作用类型均不产生净力矩，因为受力要么共线，要么对称分布以抵消旋转效应。
• 系统中所有原子受力的总和为零，证实机械平衡在模拟过程中得以保持。
• 受力推导确保能量守恒，因为受力施加过程中未引入非物理能量。
• 该方法采用与坐标无关的矢量形式实现，适用于稳健且可移植的分子动力学代码开发。