基于微观相场DO22原子间相互作用势的计算.pdfVIP

第39卷 第ll期 稀有金属材料与工程 V01．39，No．11 RARE 2010年 11月 METALMATERIALSANDENGINEERING November2010 基于微观相场D022原子间相互作用势的计算模型 董卫平，陈铮，王永欣，徐聪，卢艳丽 (西北工业大学凝固技术国家重点实验事，陕西西安710072) 的计算模型，其输入参数为相变点的温度和原子浓度。针对Ni3V-D022相的反演计算表明：随温度和浓度增大，第一近 邻原f间相互作用势W1增大。将计算所得随温度和浓度变化着的原F间相互作用势代入微观相场方程中，表明该方法 得到的原了演化结构的温度和浓度相关性更加明显。 关键词：原f间相互作用势；微观相场：D022相；温度；浓度 中图法分类号：0572．2 文献标识码：A 文章编号：1002-185X(2010)11-1970．04 从原子层次上对凝聚态物质的性质和过程进行计 展，并在应用中随温度、压强等参数变化而不断修正。 算机模拟，其结果的准确性和计算速度很大程度上依 近来陈铮课题组已经进行了L12相原子间相互作 赖于原子间相互作用势的精度，因此原子间相互作用 用势的研究且模拟结果与实际情况相符合[12】。本研究 势的研究一直是研究领域的一个重要课题。 就是基于此课题组利用微观相场原理进行多年研究的 原子间相互作用势的研究始于20世纪20年代， 基础上，研究了D022相第一近邻原子间相互作用势的 20世纪40年代以后由于计算机的出现开始了计算机计算以及应用计算出的原子间相互作用势在微观相场 模拟【¨。随着计算机在科学研究中的应用，于20世纪 方法中模拟合金的沉淀过程． 50到60年代初出现了原子层次的计算机模拟。在80 1 公式推导 年代初，大多数经验原子间相互作用势分成两组。一 组是由对势组成，最著名的是Lennard—Jones“6．12”本研究第一近邻原子间相互作用势计算模型推 势和指数形式的Morse势，另一组势函数是用来精确 描述更为复杂的系统，如金刚石结构的半导体基态性 质，最著名的是Keating模型。另一方面，对势采用 另一种形式展开成N．体势的形式[21。近年来，已经建 子(B原子)占位几率和自由能之间关系方程： 立了许多计算原子间相互作用势的模型，其中有： EAM，MEAM，MAEAM，F—S，SMA—TBt孓8】和CW， cI1+∑仉E(r)I= ＼ 5‘1 ／ CPA，KKR．CPA[明以及XEFS，EXAFS[m】等，都得到 ，l、 W Q等人通过一系列的第 了很好的结果。国内Zhang H坐竽掣]+l『¨¨ 一性原理能带计算方法，如LAPW方法、点缀球面波 (ASW)方法、正则守恒赝势及超软赝势法等，得到了式中，C为溶质原子(B原子)的浓度，桃是长程序参 各种结合能曲线，并依据陈氏反演得出SC、BCC、数，疋(力是与晶格对称性有关的函数，∥是A、B原子 FCC、HCP、L12、B2、Diamond、。c—ZnS、∥．ZnS、NaCl 之间的化学势，卜1是超晶格结构非零向量‰的数目， 型等多种晶体结构的反演系数，得到同种及异种原子 ％为玻尔兹曼常数，r为绝对温度。 间相互作用对势，并将它们应用于各种复杂体系[11】。 根据式(1)进一步得出了D022相的占位几率和能 不难看出，人们一直寻找着一种通用的势函数形式， 量之间关系方程： 但没有一种势是完全适用于不同结构的，它将不断发 收到初稿日期