第六章-分子动力学方法.pdfVIP

分子动力学方法 （ＭＤ） 1.引言 什么是分子动力学方法(MD)?  确定论的模拟方法实际上在本书第一章介绍的的第 一个计算机模拟工作,即费米--巴斯塔--乌拉姆实 验中已经产生了. 但是在统计物理中所用的分子动 力学方法则产生于1957年. B.J.Alder等首先用此 方法研究了硬球气体的状态方程.  经典分子动力学方法就是通过直接数值解N个分子 的运动方程得到这些分子每个时刻的坐标与动量， 即相空间中的轨迹，从而用统计方法计算出多体系 统的宏观性质。 分子：微观客体 （不一定是分子）. 运动方程：经典牛顿方程或哈密顿方程. 数值积分：欧拉法、龙格--库塔法、辛算法等. 相互作用：Lennard--Jones 势等. 初始条件：有序、无序等. 边界条件：周期性、全反射、漫反射等. 经典MD 平衡态MD 微正则系综 （NVE) 正则系综 （NVT） 等温等压系综 （NPT） 等焓等压系综 （NPH） 非平衡态 MD 量子MD （密度泛函分子动力学，第一原理计算） MD与MC的比较 分子动力学方法实现了 Boltzmann 的统计力 学途径,而蒙特卡罗方法实现了 Gibbs 的统计力学 途径. 分子动力学方法可处理与时间有关的过程，因 而可处理非平衡态，是更有力的模拟方法，但计算 量大、占内存多、程序也较复杂; 蒙特卡罗方法则具有程序简单、占内存少的优 点,但难以处理非平衡态. 2.单粒子问题初值问题的数值解法 2.1.一个例子一维简谐振子 我们先通过简谐振子的例子看看如何用计算机数值解单粒子问题, 一维简谐振子的哈密顿量为 p 2 1 2 H  kx 2m 2 初始条件为x(0), p (0) ，哈密顿量 （能量）为守恒量。 运动方程可用拉格朗日形式和哈密顿形式，等价，但计算方法有 区别，哈密顿形式对时间一阶导数： dx H p dt p m dp H  kx dt x 用数值积分计算在相空间中的轨道 (x(t), p (t)) 方法是将微分变成差分，引入时间步子 h 作泰勒展 开 df [f (t h) f (t)]/ h dt dx [x (t h) x (t)] p (t) 于是，哈密顿方程变为  dt h m dp [p (t h) p (t)]  kx(t)