《物理问题的计算机模拟方法(1)—分子动力学.doc》由会员分享，可在线阅读，更多相关《物理问题的计算机模拟方法(1)—分子动力学.doc（22页珍藏版）》请在金锄头文库上搜索。
1、硕士研究生课程物理问题的计算机模拟方法讲义适用专业: 凝聚态物理、材料物理与化学、理论物理、光学工程学 时：3040学时参考教材： 1德D.W.Heermann 著，秦克诚译，理论物理中的计算机模拟方法，北京大学出版社，1996。 2荷 Frenkel & Smit 著，汪文川 等译，分子模拟从算法到应用，化学工业出版社，2002。 3M.P.Allen and D.J.Tildesley, Computer Simulation of Liquids, Clarendon Press, Oxford, 1989. 4. A.R.Leach, Molecular Modelling: Prin
2、ciples and Applications, Addison Wesley Longman, England, 1996. 5. 德 D.罗伯 著，计算材料学，化学工业出版社，2002。 6. 英 B. Chopard & Michel Droz 著，物理系统的元胞自动机模拟，祝玉学，赵学龙 译，清华大学出版社，2003。目 录第一章 计算机模拟方法概论1.1 序言1.2 热力学系统物理量的统计平均1.3 分子动力学方法模拟的基本思想1.4 蒙特卡罗方法模拟的基本思想1.5 元胞自动机模拟的基本思想1.5.1 简要的发展历程1.5.2 简单元胞自动机：奇偶规则1.5.3 元胞自动机的一般定
3、义第二章 确定性模拟方法分子动力学方法（MD）2.1 分子动力学方法2.2 微正则系综分子动力学方法2.3 正则系综分子动力学方法2.4 等温等压系综分子动力学方法第三章 随机性模拟方法蒙特卡罗方法（MC）3.1 预备知识3.2 布朗动力学（BD）3.3 蒙特卡罗方法3.4 微正则系综蒙特卡罗方法3.5 正则系综蒙特卡罗方法3.6 等温等压系综蒙特卡罗方法3.7 巨正则系综蒙特卡罗方法第四章 离散性模拟方法原胞自动机（CA）4.1 引言4.2 元胞自动机模拟*4.3元胞自动机模拟的应用第一章 计算机模拟方法概论 1.1 序言1 什么是计算机模拟？ Simulation Modelling2为什
4、么要进行计算机模拟？3常用的计算机模拟方法确定性模拟方法：MD模拟 Molecular Dynamics随机性模拟方法：MC模拟 Monte Carlo离散性模拟方法：CA模拟 Cellular Automata 1.2 热力学系统物理量的统计平均描述系统的坐标（自由度）：x(t)=x1(t),x2(t),xN(t)系统的物理量：A(x(t)1时间平均 分子动力学（MD）模拟 (1-1)2系综平均 蒙特卡罗(MC)模拟 (1-2) 分布函数（几率密度函数） (1-3) 配分函数 (1-4) 相空间 H(x)系统的哈密顿函数对于处于平衡态的系统，可以证明： 对于实际的有限时间内的平均，则有 实际
5、模拟的系统大小也是有限的：有限的粒子数N或有限的系统限度L对统计平均结果有影响。 1.3 分子动力学(MD)方法模拟的基本思想1 基本原理系统：N个粒子，体积V，粒子质量为m描述一个粒子运动状态的自由度：（ri, pi） （pi=mvi）相空间：6N维，相空间中的一点的坐标 XN=rN, (mvN) rN=(r1, r2, , rN), vN=(v1, v2, , vN)粒子间的相互作用势：U(rN)=U(r1, r2, , rN)=决定系统相轨迹XN（t）的运动方程： (1-5)物理量A的宏观值，由A(XN) 的时间平均获得，即 （离散情况：）对于平衡态： 实际模拟时间总是有限的，模拟时间的
6、长短可通过判断时间的增加对平均值的影响来确定，当继续增加时间带来的平均值得变化在允许的误差范围之内时，即可认为模拟足够长了。2 计算步骤运动方程： 即 (1-6) 或 (1-7)数值求解：用差分近似表示微分 （采用不同的差分格式，可得到不同的算法）。用显示中心差分格式，将（7）式写为 (1-8)由（7）和（8）式可得： (1-9)第一步：由（9）式计算第i个粒子在t+t时刻的位置坐标。 要启动计算，我们必须要知道最初两点ri(0)和ri(t)第二步：对不同时刻 t = t , 2t, 3t, , Lt (t0 = 0) 计算物理量A(r1(lt), r2(lt), , rN(lt) (l =
7、1, 2, , N)第三步：计算物理量A的平均值 L的大小由继续增大L 而不变（或变化在误差范围内）来确定。 1.4 蒙特卡罗（MC）方法模拟的基本思想1 基本原理以正则系综（T, V, N）为例正则分布：正则配分函数：系统能量：物理量：A(rN) = A(r1, r1, , rN )系综平均： （1-10） （位形积分） （1-11）用MC方法计算上述多维积分。2 计算步骤（1） 划分原胞 N个粒子3N个空间自由度，3N维空间划分成s个相等的原胞（s1）注意：由于积分中不含动量，所以我们只需要在位置空间积分，而不需要在相空间中积分。当系统的代表点落入第i个原胞时，则认为系统处在状态i，因此，
8、s为系统可能的微观状态数目。于是，积分（10）和（11）可近似表述为 (1-12) (1-13)（2） 建立马尔可夫(Mako)过程（链）将s个状态看作一组随机事件马尔可夫链：从状态ij状态j(eiej)的概率pij ，只与ei和ej 有关。 ，i=1, 2, , s若ei 经历n步到达ej ，其概率表示为，存在极限概率 ( j=1, 2, , s)uj为系统处在状态j的概率。于是，沿无限长的马尔可夫链，物理量A的平均值可写为 (1-14) 选取 ，则（14）式为A的正则系综平均值。（3）抽样方法 采用怎样的抽样方法所构成的马尔可夫链能得到上述平均值？ 粒子位置坐标：粒子编号：g =1, 2,
9、 N坐标的三个方向：a = 1, 2, 3 系统状态：i = 1, 2, , s给定粒子位置坐标的变化量d (小于系统体积的限度)给定系统的初态i，随机选定4个随机数，其中三个xa （a =1， 2，3），且 1 xa 1，一个表示粒子编号 g =1, 2, N，由此随机确定粒子位置的变化： （ 确保 ）若，则 g 运动到新的位置，即系统由状态i过渡到状态j；若，则再选一个随机数x 4（0 x 4 1），若，则粒子保留在原位置，不发生ij的跃迁；若，则发生i j的跃迁。由此进行下去，则形成一个马尔可夫链（或过程），此链的长度L(即粒子行走的步数，远大于s)，由所计算的物理量的平均值 (1-15
10、)不再随链的加长而改变来确定。由此得到的平均值即正则平均值。一般来说，L与N, V, T 有关，比如，N=32108, L=30005000。归纳起来，计算系统物理量的正则系综平均值的具体步骤如下：第一步：给定系统的初始状态(粒子的初始位置)ri和每一步的改变量d；第二步：选择四个随机数，其中一个代表粒子的编号i ( 1 i N )；另外三个表示粒子空间坐标的改变dx, dy, dz ( -d da d , a=1, 2, 3);第三步：计算粒子i的新位置 第四步：计算粒子在新旧两个位置系统的能量之差 第五步：由的大小判断粒子i是否从ri运动到: 若，则ri； 若，则再选一个随机数R(0 R
11、1), 如果 ，则ri； 如果 ，则ri 不变，返回第二步。第六步：计算第七步：重复上述各个步骤，直到完成L步为止，最后利用公式（15）计算A的平均值。3 粒子间相互作用势模型的选取最简单的两种模型： （1）硬球模型 (s 为硬球的直径) （2）LJ 势 1.5 元胞自动机(CA)模拟的基本思想元胞自动机：时间和空间都离散、物理参量只取有限数值集的物理系统的理想化模型cellular automata 或 cellular automaton CA1.5.1 简要的发展历程1自繁殖系统 20世纪40年代，Von Neumann, 构造能解决非常复杂问题的计算机，设想模仿人脑的行为寻求与生物过程
12、无关的情况下自繁殖机理的逻辑抽象。 根据S. Ulam的建议，Von Neumann 在由元胞构成的完全离域的框架下处理这个问题，构造了一个完全离散的动力学系统元胞自动机。第一个自复制元胞自动机由二维方形网络组成，有数千个基本元胞构成的自繁殖结构。（1）一般机器只能构造比自己简单的客体，而采用自复制元胞机，可获得一种能产生新的、具有同样复杂性和功能的“机器”；（2）Von Neumann 的元胞自动机规则具有所谓通用计算的性质，这意味着，存在一种元胞自动机的初始构型，该元胞自动机能产生任何计算机算法的解。通用计算的性质指：用元胞自动机演化规则能够模拟任何计算机流程（逻辑选择器开关）。2生命游戏机1970年，数学家John Conway 生命游戏机的概念，寻找能导致复杂行为的简单规则。设想一个类似于棋盘的二维方形网格，每个元胞可能的状态是活（状态1）或死（状态0），其更新规则是：有三个活元胞包围的一个死元胞恢复为活元胞；由两个以下或三个以上活元胞包围的活元胞因孤立或拥挤而死亡。结果表明，生命游戏机有出乎意料的丰富行为，从原“汤”中显示出
