过渡金属氮化物中单原子固溶和力学性能的第一性原理研究.pdfVIP

第30卷第3期 计 算 物 理 V01．30，No．3 CHINESE OFCOMPUTATIONALPHYSICS 2013年5月 JOURNAL May，2013 文章编号：1001—246X(2013)03-0433-08 过渡金属氮化物中单原子固溶和力学性能的 第一性原理研究 刘学杰， 任 元 (内蒙古科技大学机械工程学院，内蒙古包头014010) 晶体中固溶的稳定结构，讨论置换型和间隙型固溶的低能量稳定结构与晶体间距的关系，研究金属氮化物和固溶 况为，单原子不进入对应过渡金属氮化物晶体中形成间隙固溶或置换固溶，随着晶体间距离变化单原子可以在晶 体之间形成间隙固溶或置换固溶；Si原子在NbN以及B原子在TiN晶体中可以实现间隙固溶，而不能形成置换固 溶；si原子在VN和c原子在TiN晶体中固溶结构形式均为置换固溶．单原子固溶形成低能量置换型固溶体和间 隙型固溶体的弹性常数、体模量和剪切模量均低于原过渡金属氮化物的对应值． 关键词：单原子；过渡金属氮化物；第一性原理；固溶结构；力学性能 中图分类号：0484．2 文献标识码：A O 引言 纳米复合结构超硬性能的研究集中在两个方向：纳米多层膜超硬结构和纳米复合超硬结构．在纳米多 层膜超硬结构的研究中，过渡金属氮化物多层膜超硬性质的研究已经成为研究的热点．例如TiN、TaN、NbN、 GPa～105 GPa的超高硬度的Ti- 他元素形成纳米复合结构的超硬表面．在2000年，Veprek等¨。报道获得80 GPa～110 si-N纳米复合薄膜，硬度值达到了金刚石的硬度(70 GPa)，这一结果激发了该方面的研究热 情‘4—5I． 薄膜抵抗弹性和塑性变形的能力直接决定其硬度的高低，薄膜的变形形式为位错和边界滑移。61．晶粒 尺寸减小使得晶界数量增多从而阻止晶粒位错开动，导致薄膜的硬度升高．但是，随着薄膜中晶粒尺寸的减 小至小于10nm时。7o，薄膜的变形形式由位错转化为边界滑移和晶粒转动，晶粒尺寸继续减小薄膜硬度开始 膜的超高硬度除了取决于TiN纳米晶粒尺寸，由Si形成的界面相结构及其强度同样是重要的影响因素．本 文考察低能量的固溶结构，且比较晶粒相(即过渡金属氮化物)和固溶相的力学性能． 采用实验方法暂时无法提供准确的掺杂元素Si在Ti族、V族氮化物晶体中的固溶形式及力学性能等信 息，而第一性原理计算在这方面的研究中显示出优势，广泛应用于计算多原子系统的电子结构和预测物质结 构的稳定性∽o．本文采用基于密度泛函的第一性原理方法，主要研究6种过渡金属氮化物中，固溶单原子位 于4个不同位置的最低能量与晶体间距的关系，考察可能存在的固溶结构形式以及低能量结构．进而，对过 渡金属氮化物与低能量固溶结构的力学性能进行研究．本文提出一种计算纳米复合薄膜中界面结构预判方 法． 收稿El期：2012—09—06；修回日期：2012—12—10 作者简介：刘学杰(1953一)，男，教授，主要从事纳米材料计算方面研究，E·mail：xuejieliu2000@yahoo．corn 万方数据 434 计 算 物 理 第30卷 1 计算方法 基于密度泛函理论(DFT)第一性原理赝势平面波方法广泛地被用来计算晶体电子结构¨…，本文所有的 计算通过VASP软件包实现，采用周期边界条件进行计算．真实原子有效势(包括原子核对价电子的库仑势 和芯电子对价电子的等效排斥势)采用赝势代替…1，电子波函数采用平面波作基组展开，交换关联能使用 GGA的PW91进行处理．选取5×5×5的Monkorst．parkK点对全布里渊区求和．迭代计算时电子和离子弛