它等于常数-固体微结构研究室.pdfVIP

第四章 分子动力学方法 §4.1 分子动力学方法 N 2 1 2 Δr t ⎡r t =−r 0 ⎤ (-3) () ∑ () ( ) ⎣ i i ⎦ N i 1 对于固体，它等于常数（图-2 ）。而在液态时，它是与时间成正比的，与 （-7 ）式相同，这是因为液态中的原子可以流动扩散。 第三个因子是对关联分布函数g r （-2 ）式，晶态时它有特定的结构。 ( ) 4.1.4 运动方程积分算法 Verlet 型算法 计算了粒子间的相互作用力后，就可以对Newton运动方程积分。最简单通 常也是最好的方法是用Verlet算法，将（-3 ）式写为 x t +h 2x t =−x t −h +h2 F t m +O h4 ， (-1) ( ) ( ) ( ) ( ) ( ) 在对位置的大量积分步之后，其坐标的累计误差将达到O h2 。该算法在时间反 ( ) 演（h →−h ）下是不变的，意味着将运动方程积分一段时间后再用当前的两个 点起始做时间反演的话，相空间中的轨迹将完全沿原路返回。 注意Verlet算法中并不使用速度来计算新的位置。但我们需要知道速度以计 算动能K ，以此来监测总能E K +V 是否守恒，它是分子动力学中验证计算结 果是否正确的最重要判据。根据空间轨迹的知识可以导得速度，它为中心差分 x t +h −x t −h ( ) ( ) 2 v t =+O h ， (-2) () 2h ( ) 即用该算法得到的速度是一阶的。Verlet算法中的速度和位置是不对称的，最早 的方法中速度远没位置重要，这是与相空间中轨迹的速度坐标和位置坐标等价的 概念不符的，更为现代的方法已经改正了这点。 利用半步长的速度值， x t +h −x t ( ) ( ) 2 v t +h 2 =+O h ， (-3) ( ) h ( ) 可以将Verlet算法改写成对精度误差敏感性较小的蛙跳法， v t +h 2 v t =−h 2 +hF t m +O h2 ( ) ( ) ( ) ( ) 。 (-4) x t +h x t =+hv t +h 2 +O h4 ( ) () ( ) ( ) 该方法是从Verlet算法导出的，产生相同的空间轨迹。需注意的是，速度和位置 并未在同一时间定义，结果是动能和势能也并未同时定义。因此蛙跳法中不能直 接计算总能量。 还可以将Verlet算法重写成一种同时计算位置和速度的形式－速度Verlet 法。 坐标用Taylor展开， x t +h x t =+v t h +h2 F t 2m ， (-5) ( ) ( ) ( )