《分子动力学》doc版.doc

袃羄莂蒀羅腿芈葿蚅羂膄薈螇膇肀薇衿羀荿薆蕿膆芅薅螁羈芁薅袃芄膇薄羆肇蒅薃蚅袀莁薂螈肅芇蚁袀袈膃蚀薀肃聿虿蚂袆蒈虿袄肂莄蚈羇羄芀蚇蚆膀膆蚆蝿羃蒅蚅袁膈莀螄羃羁芆螃蚃膆膂莀螅罿肈荿羈膅蒇莈蚇肈莃莇蝿芃艿莇袂肆膅莆羄衿蒄莅蚄肄莀蒄螆袇芆蒃袈肂膂蒂蚈袅膈蒁螀膁蒆蒁袃羄莂蒀羅腿芈葿蚅羂膄薈螇膇肀薇衿羀荿薆蕿膆芅薅螁羈芁薅袃芄膇薄羆肇蒅薃蚅袀莁薂螈肅芇蚁袀袈膃蚀薀肃聿虿蚂袆蒈虿袄肂莄蚈羇羄芀蚇蚆膀膆蚆蝿羃蒅蚅袁膈莀螄羃羁芆螃蚃膆膂莀螅罿肈荿羈膅蒇莈蚇肈莃莇蝿芃艿莇袂肆膅莆羄衿蒄莅蚄肄莀蒄螆袇芆蒃袈肂膂蒂蚈袅膈蒁螀膁蒆蒁袃羄莂蒀羅腿芈葿蚅羂膄薈螇膇肀薇衿羀荿薆蕿膆芅薅螁羈芁薅袃芄膇薄羆肇蒅薃蚅袀莁薂螈肅芇蚁袀袈膃蚀薀肃聿虿蚂袆蒈虿袄肂莄蚈羇羄芀蚇蚆膀膆蚆蝿羃蒅蚅袁膈莀螄羃羁芆螃蚃膆膂莀螅罿肈荿羈膅蒇莈蚇肈莃莇蝿芃艿莇袂肆膅莆羄衿蒄莅蚄肄莀蒄螆袇芆蒃袈肂膂蒂蚈袅膈蒁螀膁蒆蒁袃羄莂蒀羅腿芈葿蚅羂膄薈螇膇肀薇衿羀荿薆蕿膆芅薅螁羈芁薅袃芄膇薄羆肇蒅 第五章 分子动力学第一节 Verlet算法 牛顿方程 记 方程写为三点公式 如果给出初始条件和，可求解方程，但常常给出的初始条件是，那么 (为什么？ 因为，所以，所以，)方法的优点：保持时间反演不变性，即令 , 方程形式不变 （尽管误差会破坏这一对称性）如果问题与无关，计算精度相当高方法的缺点： 必须用到（为什么是缺点？）另一方案 缺点：失去时间反演不变性第二节 多体问题的基本方法 （阅读材料）全同粒子，概率分布为 物理量平均值分子动力学 个粒子处于的分布密度函数来自个粒子中取个的组合数例如：是1 是通常记，称系统的粒子密度定义则证明：这是显然的这里假设了是关于交换和对称的还可证明证明：如如多出一项，来自 的贡献。我们定义粒子对分布函数 如下 当系统的密度比较均匀时，退化为粒子对分布函数包含体系丰富的关于平移对称性的性质l 对固体，粒子对分布函数在晶体格距呈现尖锐峰值l 对液体，分布函数只呈现平坦峰值，而且随距离迅速消失类似地，还可以定义关于对称性的物理量。第三节 分子动力学的简单应用1二维固液相变的磁偶极子模型 Hamiltonian H=K+VK 是动能项，势能项 在实际模拟中，为了节省计算时间，可以切断相互作用的力程。但无论如何，带有相互作用的系统的模拟比硬碟模型困难多了。 我们特别关注对称性空间关联函数时间关联函数 数值模拟结果与实验结果较好吻合 2二维理论的Hamiltonian 动力学假设是孤立系统，Hamiltonian为其中，Hamiltonian 方程为应当指出，这里我们已经把定义在格点上。在连续极限下，这便是Ginsburg-Landau 理论。应用l 场论l 宇宙学l 统计物理学l 凝聚态物理学 .Verlet算法在相变点附近，由于动力学慢化，求解方程到平衡态比较困难。点阵太小，存在有限点阵效应。点阵太大，关联时间长，难以达到平衡态，误差难以控制。如果我们已经非平衡态动力学，这一困难不存在。假设初始状态是高温态，即随机态。我们测量宏观物理量，如磁化等，随时间的演化，可以确定相变点以及相关的临界指数。物理量的测量，例如，磁化强度和它的二次矩 ， k=1,2自关联函数 磁化的标度行为从这式子我们可以测量相变点（即相变能量），指数和 从时间自关联函数和磁化的二次矩可以测量指数z 和 结果可以和Ising模型以及Monte Carlo 动力学比较1.252.165(10).191(1) Ising.95(5).24(3)2.148(20).176(7) Zq 关键是Lorentz不变性被破坏，所以，3 一维热传导的简单模型热传导已经是一个古老的物理问题。现在人们对它又感兴趣，一方面是纳米材料的兴起，另一方面是低维热传导有些不同于高维的特点，如热传导系数发散等。在环境温度差的驱动下，产生能量的定向流动，由能量守恒，我们得到热传导方程 其中f(x,t) 是能量密度分布函数，是能流密度矢量。在稳态时，Fourier定律假设 常数k 称热传导系数。对一维系统，k发散。一个简单模型一根空心管，管内壁设置一些障碍物，最简单情形，是一些半园。管子两端分别射出一些粒子，出射粒子的速度由两端的温度决定。温度高的粒子速度快，温度低的速度慢。用分子动力学方法模拟粒子的运动，可以看到能量从高温端向低温端传递。按照温度是平均动能的概念， 再测量能流密度， 从而计算热传导系数。一般地， 其中L是体系的尺寸，是正数，其数值与体系有关。参考文献：D. Alonso, R. Artuso, G. Casati, I. Guarneri，Phys. Rev. Lett. 82, 1859 (1999)小结：l 分子动力学方法 求解多粒子系统的基本微观运动方程 广泛应用 比较耗时，误差有时不易控制l Monte Carlo 方法 求解多粒子系统的平衡态或非平衡态问题 处于微观或介观层次 较广泛应用 简单实用，比较节省时间l 有限元方法 求解宏观或介观运动方程例如，静电势的Poisson方程把空间分割成许多小块，每块用坐标 标记。设其中定义于 附近的局域函数。 显然，如果n足够大，可以逼近方程的解。如果n有限，记方程的误差为 现在，我们的目标是选取恰当的使极小。例如，引入其中是一个权重函数，然后取使为零。这样，条件便等价于一个n元的线性方程组 Aab a是的列矩阵，而例如，Galerkin方法设，取，这里试题：I （50分）1）设积分 ，试证明 ，其中 。2） 设 ，具体写出上述两个表达式。II （50分）1）设积分 ，假设我们可以按照分布W(x)得到，则，如果用Markov过程产生，转移矩阵应当满足什么条件？2）设，写出相应的Metropolis算法的转移矩阵。 羆肈节螄羅膀蒈蚀肄芃芀薆肃羂蒆蒂肂肅艿袁肁芇薄螇肀荿莇蚃肀聿薃蕿聿膁莅袇肈芄薁螃膇莆莄虿膆肆蕿薅螃膈莂薁螂莀蚇袀螁肀蒀螆螀膂蚆蚂蝿芄蒈薈螈莇芁袆袇肆蒇螂袇腿芀蚈袆莁蒅蚄袅肁芈薀袄膃薃衿袃芅莆螅袂莈薂蚁袁肇莄薇羁膀薀蒃羀节莃螁罿羂薈螇羈膄蒁蚃羇芆蚇蕿羆莈葿袈羆肈节螄羅膀蒈蚀肄芃芀薆肃羂蒆蒂肂肅艿袁肁芇薄螇肀荿莇蚃肀聿薃蕿聿膁莅袇肈芄薁螃膇莆莄虿膆肆蕿薅螃膈莂薁螂莀蚇袀螁肀蒀螆螀膂蚆蚂蝿芄蒈薈螈莇芁袆袇肆蒇螂袇腿芀蚈袆莁蒅蚄袅肁芈薀袄膃薃衿袃芅莆螅袂莈薂蚁袁肇莄薇羁膀薀蒃羀节莃螁罿羂薈螇羈膄蒁蚃羇芆蚇蕿羆莈葿袈羆肈节螄羅膀蒈蚀肄芃芀薆肃羂蒆蒂肂肅艿袁肁芇薄螇