分子动力学在材料科学中的应用.docVIP

分子动力学在材料科学中的应用 摘要： 本文综述了种常见条件下的分子动力学模拟方法以及分子动力学模拟的最新发展趋势.介绍用分子动力学模拟方法研究固休的休相结构,表面问题,界面问题以及薄膜形成过程等方面的研究成果 关键词分子动力学; 计算机模拟; 材料 模拟中，MD采用周期边界条件和最小镜像原理，可以大幅度减少计算工作量［1］。周期边界条件是将一定数量的粒子N集中在一定的容积V中，这个容积V称为原胞，原胞周围的部分可以看作是原胞的复制，它们称作镜像细胞。这些镜像细胞的尺寸和形状与原胞完全相同，并且每个镜像细胞所包含的N个粒子是原胞中粒子的镜像，原胞在各个方向周期复制便形成了宏观物质样本。这样只需根据原胞周围的边界条件计算原胞内粒子的运动，幅度减少了工作量。 原子间作用势能模型的构造对于;MD法的应用至关重要。最简单的偶势模型只考虑两体作用，而与其它原子无关，在模拟中运算量小。20世纪80年代以来，各种经验或半经验的多体势模型迅速发展，特别是镶嵌原子法（EAM)［4］既克服了偶势的缺陷，又不会使计算量太大。 3分子模拟主要技术细节 3．1 分子间作用势 分子动力学模拟的首要条件就是要知道分子间的相互作用势，分子间作用势函数确定后，通过势函数对“求导即可得出分子间的作用力。所以必须对分：子间的作用势进行研究，目前存在多种势能模型。分子间势函数的发展经历了对势一多体势的过程。对势认为原子之间的相互作用是两两之间的作用，与其他原子的位置无关；而实际上，在多原子体系中，一个原子的位置不同，将影响空间一定范围内的电子云分布，从而影响其他原子之间的有效相互作用，故多原子体系的势函数需更准确地用多体势表示。 在分子动力学模拟的初期，人们经常采用的是对势。常用的对势模型包括Lennard—Jones势、Morse势、Born—Lande势及Johnson势。在对势模型中，系统能量表示为原子对相互作用能量的加和。对于简单的、高对称性的晶体，对势能够描述原子间的相互作用。对势模型的主要缺点是忽略了多体原子间的相互作用，更具体地说，它是忽略了由于原子局域环境的变化引起的原子间相互作用的变化，因此需要建立新的模型。 基于EAM势的势函数还有很多种。这些多体势大都用于金属的微观模拟。为了将EAM势推广到共价键材料，需要考虑电子云的非球形对称。于是，Baskes等提出了修正型嵌入原子核法(MEAM)。经过修正，Baskes理论已能解决上述问题，但其确定参数的过程相当复杂，应用上仍受到很大的限制。Pasianot等试图在总能量中加入一修正项，以修正原EAM模型中2条基本假设所引起的能量差。这一模型成功地描述了Nb、Fe、Cr等元素，但应用起来很繁杂。张邦维等综合考虑以上EAM模型的优缺点后，提出了分析型EAM理论，成为一个普适分析型EAM模型。胡望宇等在张邦维分析型EAM模型的基础上进行了改进，提出了改进的分析型EAM模型(MAEAM)。此外，还有许多形式的多体势函数形式，如Jacobsen等在等效介质原理(EMT)的基础上提出的另一种函数形式，由于其简单、有效，因此也得到了广泛的应用。势函数确定后，通过对势函数求导即可得出分子间的作用力。 3．2周期性边界条件 原子和分子体系的分子模拟的目的是提供一个宏观样本的物性信息。在具有自由边界的三维N个粒子的体系中处于界面的分子数正比于N1脂。在使用有限的原子数来模拟实际体系中原子的运动时，必须考虑表面对体系中原子运动的影响。为避免这种影响，可以通过周期性边界条件来实现。将含有N个粒子的体积当作具有与其相同单元无限周期点阵的原始单元(如图1是一个二维的示意图，在这个二维图象中每个单胞被其他的8个单胞所包围；在三维方向上每个单胞就会被26个单胞所包围)，一给定粒子则与在此无限周期体系中的全部其他粒子相互作用。对有些模拟，在所有方向都用周期性边界条件是不合适的。比如在研究表面的分子吸附时，在与表面垂直方向上不能用周期性边界条件，而仅在平行于表面的2个方向需要应用周期性边界条件。 4 分子动力学在材料科学中的应用 4.1 分子动力学的适用范围 材料科学中计算模拟研究范围极为广泛，从埃量级的量子力学计算到连续介质层次的有限元或有限差分模型，可分为4个层次：电子、原子、显微组织和宏观层次（见图2）。;MD主要是原子尺度上研究体系中与时间和温度有关的性质的模拟方法。 最早将分子动力学方法用于材料研究中的是Vineyard 于1960年探讨材料辐射损伤的动力学规律。模拟结果给出了原子轨迹，这一工作使得过去对热力学性能的定性估计迈向对微观过程的定量研究.1964年Rahman用MD方法模拟液体氩，同时加进了周期性边界条件，结果他惊奇地发现可以用很少的粒子（864个）来反映真实系统的热力学性质。自此，凝聚态物质