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

1、分子动力学在材料科学中的应用摘要：本文综述了几种常见条件下的分子动力学模拟方法以及分子动力学模拟的最新发展趋势.介绍用分子动力学模拟方法研究周休的休相结构，外表问题，界面问题以及薄膜形成过程等方面的研究成果。关键词：分子动力学；电脑模拟；材料科学1引言分子动力学MolecularDyanmica,简称MD用于计算以固体、液体、气体为模型的单个分子运动，它是探索各种现象本质和某些新规律的一种强有力的电脑模拟方法，具有沟通宏观特性与微观结构的作用，对于许多在理论分析和实验观察上难以理解的现象可以做出一定的解释：1：。MD方法不要求模型过分简化，可以基于分子原子、离子的排列和运动的模拟结果直接计算求

2、和以实现宏观现象中的数值估算。可以直接模拟许多宏观现象，取得和实验相符合或可以比较的结果，还可以提供微观结构、运动以及它们和体系宏观性质之间关系的极其明确的图象。MD以其不带近似、跟踪粒子轨迹、模拟结果准确，而倍受研究者的关注，在物理、化学、材料、摩擦学等学科及纳米机械加工中得到广泛而成功的应用。本文主要评述MD方法在材料科学中的应用.目前在材料微观结构的研究中，由于实验条件的限制，使得许多重要的微观结构的信息难以得到，如，对于由液态金属快速凝固的非晶转变过程，其微观结构的瞬时变化根本无法用实验仪器去测量。理论分析、实验测定及模拟计算已成为现代材料科学研究的3种主要方法。20世纪90年代以来，

3、由于电脑科学和技术的飞速发展，模拟计算的地位日渐突显。电脑模拟可以提供实验上尚无法获得或很难获得的信息。虽然电脑模拟不能完全取代实验，但可以用来指导实验的进行，从而促进理论和实践的发展，所以有必要对这一领域进行介绍。2分子动力学基本原理分子动力学将连续介质看成由N个原子或分子组成的粒子系统，各粒子之间的作用力可以通过量子力学势能函数求导得出，忽略量子效应后，运用经典牛顿力学建立系统粒子运动数学模型，通过数值求解得到粒子在相空间的运动轨迹，然后由统计物理学原理得出该系统相应的宏观动态、静态特性。图1所示是MD学习文档仅供参考模拟过程。MD具体的做法是电脑上求运动方程的数值解。通过适当的格式对方程

4、进行近似，使之适于在电脑求数值解。从使用连续变数和微分算符的描述过渡到使用离散变数和有限差分算符的描述，显然会有误差，误差的阶数取决于具体的近似机制，即所用的算法。模拟首先是规定初始条件。为了确定起见，可令初始位置在格子的格点上，而初始速度则从波尔兹曼分布得出。一个按上述方法建立的系统不会具有所要的能量，而且，很可能这个状态并不对应于一个平衡态。II为了推动系统到达平衡，需要一个趋衡阶段。可以通过增加或从系统中移走能量，对运动方程向前积分假设干时间步，使系统弛豫到平衡态。接着是物理量的计算阶段，沿着系统在相空间中的轨道计算一切令人感兴趣的量。初始化参数初始化原子性性分子动力学模拟图1分子动力学

5、模拟过程图Fig.1Processschemeofmoleculardynamicssimulation模拟中，MC»用周期边界条件和最小镜像原理，可以大幅度减少计算工作量1。周期边界条件是将一定数量的粒子N集中在一定的容积V中，这个容积V称为原胞，原胞周围的部分可以看作是原胞的复制，它们称作镜像细胞。这些镜像细胞的尺寸和形状与原胞完全相同，并且每个镜像细胞所包含的N个粒子是原胞中粒子的镜像，原胞在各个方向周期复制便形成了宏观物质样本。这样只需根据原胞周围的边界条件计算原胞内粒子的运动，幅度减少了工作量。原子间作用势能模型的构造对于；MD法的应用至关重要。最简单的偶势模型只考虑两体作

6、用，而与其它原子无关，在模拟中运算量小。20世纪80年代以来，各种经验或半经验的多体势模型迅速发展，特别是镶嵌原子法EAM）力既克服了偶势的缺陷，又不会使计算量太大。学习文档仅供参考3分子模拟主要技术细节3.1分子间作用势分子动力学模拟的首要条件就是要知道分子间的相互作用势，分子间作用势函数确定后，通过势函数对“求导即可得出分子间的作用力。所以必须对分：子间的作用势进行研究，目前存在多种势能模型。分子问势函数的发展经历了对势一多体势的过程。对势认为原子之间的相互作用是两两之间的作用，与其他原子的位置无关；而实际上，在多原子体系中，一个原子的位置不同，将影响空间一定范围内的电子云分布，从而影响其

7、他原子之间的有效相互作用，故多原子体系的势函数需更准确地用多体势表示。在分子动力学模拟的初期，人们经常采用的是对势。常用的对势模型包括LennardJones势、Morse势、BornLande势及Johnson势。在对势模型中，系统能量表示为原子对相互作用能量的加和。对于简单的、高对称性的晶体，对势能够描述原子间的相互作用。对势模型的主要缺点是忽略了多体原子间的相互作用，更具体地说，它是忽略了由于原子局域环境的变化引起的原子间相互作用的变化，因此需要建立新的模型。基于EAM&的势函数还有很多种。这些多体势大都用于金属的微观模拟。为了将EAM势推广到共价键材料，需要考虑电子云的非球形对

8、称。于是，Baskes等提出了修正型嵌入原子核法（MEAM）经过修正，Baskes理论已能解决上述问题，但其确定参数的过程相当复杂，应用上仍受到很大的限制。Pasianot等试图在总能量中加入一修正项，以修正原EAMK型中2条基本假设所引起的能量差。这一模型成功地描述了Nb>Fe、Cr等元素，但应用起来很繁杂。张邦维等综合考虑以上EAM模型的优缺点后，提出了分析型EAM理论，成为一个普适分析型EAM真型。胡望宇等在张邦维分析型EAM真型的基础上进行了改良，提出了改良的分析型EAMg型（MAEAM）此外，还有许多形式的多体势函数形式，如Jacobsen等在等效介质原理（EMT）的基础上提出

9、的另一种函数形式，由于其简单、有效，因此也得到了广泛的应用。势函数确定后，通过对势函数求导即可得出分子间的作用力。3.2周期性边界条件原子和分子体系的分子模拟的目的是提供一个宏观样本的物性信息。在具有自学习文档仅供参考由边界的三维N个粒子的体系中处于界面的分子数正比于N1脂。在使用有限的原子数来模拟实际体系中原子的运动时，必须考虑外表对体系中原子运动的影响。为防止这种影响，可以通过周期性边界条件来实现。将含有N个粒子的体积当作具有与其相同单元无限周期点阵的原始单元（如图1是一个二维的示意图，在这个二维图象中每个单胞被其他的8个单胞所包围；在三维方向上每个单胞就会被26个单胞所包围），一给定粒子

10、则与在此无限周期体系中的全部其他粒子相互作用。对有些模拟，在所有方向都用周期性边界条件是不合适的。比方在研究外表的分子吸附时，在与外表垂直方向上不能用周期性边界条件，而仅在平行于外表的2个方向需要应用周期性边界条件。4分子动力学在材料科学中的应用4.1分子动力学的适用范围材料科学中计算模拟研究范围极为广泛，从埃量级的量子力学计算到连续介质层次的有限元或有限差分模型，可分为4个层次：电子、原子、显微组织和宏观层次见图2。；MD主要是原子尺度上研究体系中与时间和温度有关的性质的模拟方法。学习文档仅供参考s-观层次Ip.I微观层次Ins-|原子层汨ps-电子层次|nm出mmmm图2材料模型的层次划分

11、最早将分子动力学方法用于材料研究中的是Vineyard于1960年探讨材料辐射损伤的动力学规律。模拟结果给出了原子轨迹，这一工作使得过去对热力学性能的定性估计迈向对微观过程的定量研究.1964年RahmarfflMDJ法模拟液体氮，同时加进了周期性边界条件，结果他惊奇地发现可以用很少的粒子864个来反映真实系统的热力学性质。自此，凝聚态物质的分子动力学模拟成为可能，许多研究者纷纷投入这一研究工作。最初应用是基于偶函数，如Lendard-Jones势函数和Morse势函数，模型简单，运算量小，而得以在材料科学中广泛应用。但由于其未考虑到体积相关项，在研究材料的弹性系数性质和预言金属的结合能及空位

12、形成能时，难以获得准确的结果2。EAM体势主要用于fcc型金属及其合金中，处理其结构、热力学、外表、缺陷及液态金属等问题，也应用于hcp及bcc型金属及合金，以及半导体Si。一般，MET法在中型机或微机上进行时，由于其内存和运算速度的限制，模拟研究只能限于500-1000个原子的小系统。因而模拟结果虽然也能揭示一些微观结构的特征和规律，但与实际的大系统情况有较大差异。在并行处理系统上对更大量的原子系统进行模拟研究网，其结果必然会接近于实际，从而对生产实践将会更有实际指导意义。4.2分子动力学的应用4.2.1 金属的液态结构在目前实验条件下，液态金属的结构及其变化尚很难精确测定。王鲁红等人学习文

13、档仅供参考7采用F-S型多体势描述了8种hcp型金属的液态微观结构并与实验相比较，模拟结果说明，MgCo和Zn的势函数可以较好的描述其液态结构，Ti和Zr的势函数则不能；由Be和Ru的势函数描述的液态结构较为合理，Hf则与一般液态结构特点不一致。李辉等人采用EAM多体势模型，很好地描述了液态过渡金属Ni的结构变化特性。4.2.2 薄膜形成过程薄膜研究是当今科学研究的热点之一。目前在很多薄膜制备方法中，都应用了低能离子轰击技术，在这些方法中，低能离子/外表相互作用在控制薄膜的微观结构方面起着重要作用。由于离子/外表相互作用发生在时间间隔小于10-12s内，因而特别适合于用M防法对这一过程进行描述

14、。薄腊研究是当今科学研究的热点之一目前在很多薄摸制备方法中，都应用了低能离子轰击技术，如离子束增强沉积，等离子休辅助化学汽相沉积、溅射沉积、离子辅助分子束处延等。Garrison9,Kitabtake网和郊正明四等人分别用低能粒子轰击Si(ool)-2xl外表，由此可用分子动力学方法研究低能粒子对外表原子行为的影响。郊等人的研究说明，10ev,100ev粒子的轰击，一方面增强了外表原子形成二聚体(dimer)的能力，使外表二聚体键数增加，另一方面也使外表原子的排列更趋无序。Ethier和Lewis10模拟了纯Si、SiGe和纯Ge在si(loo)-2xl再构外表上用分子束处延(MBE)t生长膜

15、的过程，其结果给出了薄膜质量与衬底温度之间的关系，即衬底温度较低时，形成的结构有序性较差；在高的衬底温度下，发生外延生长。对于纯Ge的外延生长，只有最初的三层结晶，以后便出现岛状结构，这在定性上与实验和理论结果相一致。文献中大量报道了近十儿年来关于品界的一些分子动力学模拟结果.目前有关界面的分子动力学模拟研究开展的不多.金属一金属界面的分子动力学模拟研究还有一些报道.李明研究了Ag/Ni界面处裂纹扩展行为.Ag/Ni和Cu/Ni界面在弯曲状态下的力学行为作了分子动力学模拟研究，结果给出的力学性能曲线与宏观规律相符合，并给出了进行界面模拟时计算元胞的原子数为1000个左右.同时证明了界面的存在对

16、复合材料性能影响很大，界面结构不同复合材料的性能也不相同。Yang等人研究了Ni(100)面涂弋层的结构和动力学行为，并对嵌人原子法学习文档仅供参考(EAM)和对势结果作了比较.Rahrnan等人对上述问题作了进一步研究，给出了单层Ag在Ni(loo)面上的结构与温度之间的关系。即室温下,Ag在Ni衬底上前后滑动，距Ni上外表平均距离为2.15埃，温度为1200K时,Ag在Ni上形成孔泡;1300K时,Ag在Ni上形成单晶.这些工作足以说明用分子动力学方法研究界面问题是可行的.采用分子动力学模拟方法进行界面研究是十分有意义的工作，尤其是目前正急于解决的金属一陶瓷界面问题。Parrinell和R

17、ahmanll最先采用分子动力学方法研究固体性能的，以此代替Milstein和Faber的静力学计算.其对象是Ni单晶在单轴压力下由面心立方结构向密排六方结构转变的过程，这是等温一等压分子动力学问题，所用势函数是Morse对势，其结果给出了应力一应变曲线，并与Milstein和Faber的汁算结果符合得相当好。Zhond815和RahmanParrincllo分子动力学形式。结果说明氢脆的微观起因是氢饱和使某些区域的韧性和塑性增强，这与假定的氢增强局域塑性的机制相一致。5结论从国内外分子动力学的模拟研究可以看出，随着电脑的普及及计算能力的提高，分子模拟方法的推广应用日益受到高度重视。模拟已逐渐

18、从分子结构简单的体系扩展到分子结构复杂的体系，模拟粒子数已从一二百扩展到几千，并从单相扩展到多相、从均相扩展到非均相，位能函数也从简单扩展到复杂乃至要考虑量子效应。模拟方法也已经从Mc转至MD以模拟真正的动力学过程。MD真拟方法也从平衡态模拟发展到外加场的非平衡态模拟。可以预料，分子模拟与材料的关系将更加紧密，它对材料科学工程的基础研究、工艺过程以及新产品的开发都会发挥更为明显的、不可替代的作用。参考文献学习文档仅供参考1罗熙淳，梁迎春，董中.分子动力学在纳米机械加工技术中的应用J中国机械工程,1999,10(60:6922胡壮麒，王鲁红，刘轶电子和原子层次材料行为的计算机模拟J材料研究学报，1099,12(1):13 NoseS.J.Chem.Phys.1984.85:5114 RiceBD,Chcn