令/()=x),利用牛顿第二定律公式F()=mx,公式6.2.1)写为坐标的递推公式 x(+b)=2()-x-b)+h2 (6.2.13) 公式(6.2.12)写为计算动量的公式得到 p()=m()=m()=[x(+b)-x-h) (6.2.14) 这样我们就推导出了一个比(6.2.6)和(6.2.⑦)更精确的递推公式。这是二步法的一种,称为 Verlet 方法。 当然我们还可以建立更高阶的多步算法,然而大部分更高阶的方法所需要的内存比一步法和二步 法所需要的大得多,并且有些更高阶的方法还需要用迭代来解出隐式给定的变量,内存的需求量就更 大。并且当今的计算机都仅仅只有有限的内存,因而并不是所有的高阶算法都适用于物理系统的计算 机计算。 Verlet算法是分子动力学模拟中求解常微分方程最通用的方法令 ( ) = (txtf )，利用牛顿第二定律公式 ( ) 2 2 dt xd = mtF ，公式(6.2.11)写为坐标的递推公式 ( ) () ( ) ( ) m tF hhtxtxhtx 2 2 +−−=+ . (6.2.13) 公式(6.2.12)写为计算动量的公式得到 () () () [ ] ( )( ) htxhtx hm tmvtxmtp −−+=== 2 & . (6.2.14) 这样我们就推导出了一个比(6.2.6)和(6.2.7)更精确的递推公式。这是二步法的一种, 称为 Verlet 方法。 当然我们还可以建立更高阶的多步算法，然而大部分更高阶的方法所需要的内存比一步法和二步 法所需要的大得多，并且有些更高阶的方法还需要用迭代来解出隐式给定的变量，内存的需求量就更 大。并且当今的计算机都仅仅只有有限的内存，因而并不是所有的高阶算法都适用于物理系统的计算 机计算。Verlet 算法是分子动力学模拟中求解常微分方程最通用的方法