第三章_经典分子动力学方法汇总.pptVIP

Green－Kubo formula 扩散系数D也可用速度自相关函数通过下面的Green－Kubo 公式表示 （3－3－11） Green－Kubo公式和Einstein 关系是等价的。 Partial diffusion 在多种原子组成的体系中（例如AB2），可以对某一种原子求和，得到Partial diffusion，从而得到原子A或原子B的扩散系数DA、DB。 非均匀体系中的扩散 对一非均匀体系，例如一块材料，有体和晶界部分。可分别对某一类原子做平均，如对晶界原子，或对体中的原子，得到不同部分原子不同的扩散。 与温度有关的扩散 在一个系统里，若原子的扩散是通过跳跃进行的，可将扩散用Arrhenius关系 来表达，式中是pre－factor，是原子迁移率或激发能(activation energy)。对（3－3－12）两边求对数有 (3-3-12) (3-3-13) 与温度有关的扩散 作lnD～1/T图(图3.8)，从图中直线的斜率可以得到Ea，从截距可以得到D0。同理也可以得到表面原子、晶界原子的扩散系数(图3.9)。 图3.8 图3.9 MD模拟材料的溶解 固态到液态的相变（溶解）的过程发生在几ps时间，在MD模拟中可通过以下手段知道是否发生了固-液相变： 观察不同时间段系统的原子位置情况（snapshots） 能量从动能向势能转移 总能量和势能随温度的变化有跳跃 扩散系数有变化 径向分布函数有变化 在做有固定体积的模拟时，压力有突变 MD模拟材料的溶解 溶点对应于固-液共存状态，这时在computational cell 中应含有固态、液态及固-液界面。由图3－10可以得到溶点温度。 图3－10 体膨胀系数?可以由V－T曲线（图3－11）得到 图3－11 对不同的压强P做模拟，可以得到 ?＝?（T，P）。 人有了知识，就会具备各种分析能力， 明辨是非的能力。 所以我们要勤恳读书，广泛阅读， 古人说“书中自有黄金屋。 ”通过阅读科技书籍，我们能丰富知识， 培养逻辑思维能力； 通过阅读文学作品，我们能提高文学鉴赏水平， 培养文学情趣； 通过阅读报刊，我们能增长见识，扩大自己的知识面。 有许多书籍还能培养我们的道德情操， 给我们巨大的精神力量， 鼓舞我们前进。 * 第三章 经典分子动力学方法 3.1 引言 分子动力学 (Molecular Dynamics, 简写为MD) 方法是确定性模拟方法，这方法是按该体系内部的内禀动力学规律来计算确定位形的转变。 首先需要建立一组分子的运动方程，然后通过直接对系统中的每一个原子/分子运动方程进行数值求解，得到每个时刻每个原子/分子的坐标与动量(速度)，即在相空间的运动轨迹，再利用统计方法得到多体系统的静态和动态特性，从而得到系统的宏观性质。 在MD方法的处理过程中，方程组的建立是通过对物理体系的微观数学描述给出的，在这个微观的物理体系中，每个原子/分子都各自服从经典的牛顿力学定律。 MD方法是实现玻尔兹曼的统计力学途径，可以处理与时间有关的过程，因而可以处理非平衡态问题，但是该方法的计算机程序较复杂，计算量大，占内存也多。 MD方法的发展史 MD方法是20世纪50年代后期由B.J Alder和T.E. Wainwright创造发展的。他们在1957年利用MD方法，发现了早在1939年根据统计力学预言的“刚性球组成的集合系统会发生由其液相到结晶相的相转变”。 20世纪70年代，产生了刚性体系的动力学方法被应用于水和氮等分子性溶液体系的处理，取得了成功。1972年，A.W. Less和S.F. Edwards等人发展了该方法，并扩展到了存在速度梯度(即处于非平衡状态)的系统。 之后，此方法被M.J. Gillan等人推广到了具有温度梯度的非平衡系统，从而构造并形成了非平衡MD方法体系。 MD方法的发展史 MD方法真正作为材料科学领域的一个重要研究方法，开始于恒压MD方法(1980)和恒温MD方法(1984)的建立及在应用方面的成功。 1985年人们又提出了将电子论和分子动力学方法有机统一起来的所谓Car－Parrinello方法，即第一性原理MD方法。它不仅可以处理半导体和金属的问题，同时还可应用于处理有机物和化学反应。 1991年有人进一步提出了巨正则系综MD方法，从而又可适用于吸附问题的处理等，该方法还在进一步发展之中。 分子动力学方法的主要发展可见表3.1。 年代 创立者 创造内容工作（MD分类名称） 1957 B. J Alder T. E.