分子动力学作业.doc

《分子动力学作业.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、分子动力学(MD)1分子动力学（MD）基础1.1MD分类1.2MD简介1.3MD适用范围2分子动力学运动方程数值求解2.1基础知识2.1.1运动方程2.1.2空间描述2.1.3最小作用量原理2.1.4拉格朗日(Lagrange)方程2.1.5哈密顿(Hamilton)方程2.2粒子运动方程的数值解法2.2.1Verlet算法2.2.2欧拉(Euler)预测—矫正公式2.2.3Gear预测—矫正方法3分子动力学原胞与边界条件3.1分子动力学原胞3.2边界条件3.2.1自由表面边界3.2.2固定边界3.2.3柔性边界3.2.

2、4周期性边界4势函数与分子力场4.1势函数4.1.1两体势4.1.2多体势4.2分子力场4.2.1分子力场函数的构成4.2.2常用力场函数和分类5分子动力学模拟的基本步骤5.1设定模拟所采用的模型5.2给定初始条件5.3趋于平衡计算5.4宏观物理量的计算6平衡态分子动力学模拟6.1系综6.2微正则系综的分子动力学模拟6.3正则系综的分子动力学模拟1分子动力学（MD）基础1.1MD分类微正则系综（VNE）正则系综（VNP）平衡态MD等温等压系综（NPT）经典MD等焓等压系综（NPH）巨正则系综（ＶＴμ）非平衡态ＭＤ量子MD

3、1.2分子动力学(MD)简介分子动力学是在原子、分子水平上求解多体问题的重要的计算机模拟方法。分子动力学方法为确定性模拟方法，广泛地用于研究经典的多粒子体系的研究中,是按该体系内部的内禀动力学规律来计算并确定位形的转变。分子动力学方法是通过建立一组分子的运动方程，并通过直接对系统中的一个个分子运动方程进行数值求解，得到每个时刻各个分子的坐标与动量，即在相空间的运动轨迹，再利用统计计算方法得到多体系统的静态和动态特性,从而得到系统的宏观性质。在分子动力学中，粒子的运动行为是通过经典的Newton运动方程所描述。系统的所有粒

4、子服从经典力学的运动规律，它的动力学方程就是从经典力学的运动方程——拉格朗日(lagrange)方程和哈密顿(Hamilton)方程导出。1.3适用范围原则上，分子动力学方法所适用的微观物理体系并无什么限制。这个方法适用的体系既可以是少体系统，也可以是多体系统；既可以是点粒子体系，也可以是具有内部结构的体系；处理的微观客体既可以是分子，也可以是其它的微观粒子。实际上，分子动力学模拟方法和随机模拟方法一样都面临着两个基本限制：一个是有限观测时间的限制；另一个是有限系统大小的限制。通常人们感兴趣的是体系在热力学极限下（即粒子

5、数目趋于无穷时）的性质。但是计算机模拟允许的体系大小要比热力学极限小得多，因此可能会出现有限尺寸效应。为了减小有限尺寸效应，人们往往引入周期性、全反射、漫反射等边界条件。当然边界条件的引入显然会影响体系的某些性质。2分子动力学运动方程数值求解2.1基础知识2.1.1运动方程系统的动力学机制决定运动方程的形式。在分子动力学方法处理过程中，方程组的建立是通过对物理体系的微观数学描述给出的。在这个微观的物理体系中，每个分子都各自服从经典的牛顿力学。每个分子运动的内禀动力学是用理论力学上的哈密顿量或者拉格朗日量来描述，也可以直接

6、用牛顿运动方程来描述。采用分子动力学方法时，必须对一组分子运动微分方程做数值求解。从计算数学的角度来看，这是个求一个初值问题的微分方程的解。实际上计算数学为了求解这种问题已经发展了许多的算法。但是并不是所有的这些算法都可以用来解决物理问题。2.1.2空间描述在空间描述如何物体的运动，如果其本身的大小可以忽略时，就可以将其看作是粒子（或质点）。粒子描述：空间位置：r速度：v=dr/dt加速度：若一个系统由N个粒子组成，则粒子描述：空间位置：r1,r2,r3…,rN笛卡尔坐标系，粒子有3N个自由度设系统有s个自由度广义坐标：

7、q1,q2,q3,…,qN广义速度：q1,q2,q3,…,qN2.1.3最小作用量原理莫培督1744年提出最小作用量原理：保守的、完整的力学系统，由某一初位形转变到另一位形的一切具有相同能量的可能运动中，真实的运动是其作用量具有极小值的那种运动。力学系统中，构造能量函数L及其作用量S作用量的积分式叫做泛函(functional)，作用量取极值的方法就是求其变分δS=0。2.1.4拉格朗日(Lagrange)方程由最小作用量原理可导出拉格朗日方程对于孤立的保守系统，每个粒子在势场U中运动，则系统整体的Lagrange函数是

8、得到第i个粒子的牛顿运动方程（α指每个粒子的自由度）2.1.5哈密顿(Hamilton)方程哈密顿(Hamilton)原理:保守的、完整的力学系统在相同时间内，由某一初位形转移到另一已知位形的一切可能运动中，真实运动的作用函数具有极值，即作用函数的变分等于零。哈密顿(Hamilton)方程，Lagrange函数全微分