Al-Zr二元合金的第一原理计算.pdfVIP

第 32卷 第 1期 广西大学学报 (自然科学版) Vo1．32．No．1 2007年 3月 JournalofGuangxiUniversity(NatSciEd) M ar．．2007 文章编号 ：1001-7445(2007)01-0035—04 Al—Zr二元合金的第一原理计算 欧阳义芳 ，曾凡江 ，王陆阳 ，陶小马 (1．广西大学 物理科学与工程技术学院，广西 南宁 530004；2．中南大学 材料学院，湖南 长沙 410083) 摘要：采用基于密度泛函理论的第一原理方法，计算了Al—Zr二元合金系金属间化合物的晶格常数、体弹性模 量和形成焓，计算结果与已有 的实验值和其他理论结果符合得很好． 关键词 ：Al～Zr合金 ；第一原理 ；形成焓 中图分类号 ：TU512．4；0648．12+2 文献标识码 ：A 由于铝合金在工程和基础研究中的重要性，特别是近年来对高 比强度铝基大块非晶的研究，使得对 铝基合金的研究成为材料研究领域的前沿课题．在对铝基非晶的研究中，合金热力学性质对研究非晶形 成、相图计算等所涉及的热力学参数有着非常重要的作用，因此精确确定合金体系的热力学性质是铝基 合金研究中的重点之一． 确定合金的热力学性质的方法较多，其中基于密度函数理论 (DFT)[】 的第一原理计算方法是 比 较精确的方法之一．近年来，基于密度函数理论的第一原理计算与动力学相结合，在材料设计、合成和材 料模拟等多方面有许多突破性的进展，已经成为计算材料科学的基本方法之一． 本文拟对Al—Zr合金系统的热力学性质进行研究．用第一原理方法计算Al，Zr及AI—Zr二元合金系 统的总能量、晶体结构、体弹性模量、体弹性模量对压强的微分和金属间化合物的形成焓． 1 计算方法 对纯元素和金属间化合物，其总能量采用基于密度函数理论的VASP程序[s来计算，计算 中布里 渊区积分采用Monkhorst—Pack方案Ⅲ，赝势采用投影缀加波 (PAW)方法 ，交换关联势采用Perdew— Wang广义梯度近似(GGA—PW91)[9]，能量的计算精度为10 eV，力的精度为10～eV／A．所有计算 中均 考虑 自旋，能量计算的误差小于10meV／atom．金属间化合物的形成焓为金属间化合物的总能量与纯元 素的总能量之差． 在丁一0K时，Al．Zr合金的形成焓 由下面的公式得到 △H (A11Zr)一E(A11．Zr)一 (1一 )E(A1)一xE(Zr)， (1) 式中E(AI)，E(Zr)，E(Al．Zr)分别为Al，Zr和相应的金属间化合物在平衡状态下总能量．z为合金中 Zr的原子成分． 2 结果和讨论 ．‘ Al和Zr在室温和低温下分别为fcc和hcp结构 ，在 已确定的Al—Zr相图中，稳定相有AI。Zr(tI16)， Al2Zr(hP12)，A13Zr2(oF40)，AlZr(oC8)，AI3Zr4(hP7)，A12Zr3(tP20)，A13Zr5(hP16)，A13Zr5(ti32)， A1Zr(hP6)，A1Zr。(cP4)．我们首先对纯元素及金属间化合物进行结构优化，然后采用计算其在不同体 积下所对应的总能量 ，得到总能量随体积的变化关系，通过拟合Rose状态方程n ，得到相应结构在平衡 收稿 日期 ：2006—12—14；修订 日期 ：2007—01—09 基金项 目：国家 自然科学基金 和广西 回国基金(0448007) 作者简介 ：欧阳义芳 (1965一)，男 ，湖南衡东人 。博士，广西大学教授． 广西大学学报 (自然科学舨) 第 32卷 状态下对应的平衡体积、体弹性模量及体弹性模量对压强的微分，所有 E— 关系与Rose状态方程均符 台得非常好 ，结果如表 1所示．从表 1可见，所有计算的晶格常数均 比实验值[1。稍大，这与采用GGA近 似计算得到的结果一致，计算结果与实验值的差别很小，除了Al，Zr(hP16)(误差为2．1 )外，其余误差 均小于2 ，计算值与Gb,oshn目的计算值一致．计算得到的Al