能量的梯度的负数就是受力,该粒子因为该能量项所受到的作用力.一般而言,体系有多种能量项,非成键作用如范德华作用和静电作用,成键作用如2个原子的伸缩,3个原子的键角弯曲和4个原子的二面角变化(有时还考虑还有4个原子同面的作用).你将这所有和该粒子有关的各能量项求梯度,就可以得到该项作用相应的作用力,再对所有的力进行求和则是该粒子收到的总作用力.有了作用力,a=F/m,有加速度就可以求下一步的速度,再是下一步得到位置,有了新的位置就可以求新的势能再是受力...迭代下去,这就是分子模拟的基本流程.
LZ可以去分子动力学模拟吧讨论讨论分子模拟的东西!