基于dft 计算的alnnim ( n + m = 13 ,19) 团簇的 - 河北师范大学学报.pdf

第 32 卷/ 第 5 期/ 河北师范大学学报/ 自然科学版/ Vol . 32 No . 5 2008 年 9 月 J OURNAL OF HEBEI NORMAL UNIVERSITY/ Natural Science Edition/ Sep . 2008 基于 D F T 计算的Al n Ni m ( n + m = 13 ,19) 团簇的稳态结构 李文娜 , 　魏国强 , 　刘 　英 (河北师范大学 物理科学与信息工程学院 ,河北 石家庄 　0500 16) ( ) 摘 要 :基于密度函数理论 DF T 和 ChenMbius 晶格反演方法计算了Al - Al ,Ni - Ni ,Al - Ni 原子间的相互 作用势. 根据获得的原子间相互作用势 ,通过采用最陡下降法和共轭梯度法实现对势能面全局最小点的搜索 ,获得 n + m = 13 ,19 的Ni - Al 团簇的稳态结构. 关键词 :ChenMbius 晶格反演 ; 原子间作用势 ; 稳态结构 中图分类号 :O 562 . 4 　　　文献标识码 :A 　　　文章编号 (2008) 团簇几何结构是研究团簇的各种特性和团簇生长模式的基础. 探索团簇结构的常用方法有从头计算 法[ 1 ] 、基于密度泛函的紧束缚近似[2 ] 、分子动力学[3 ,4 ] 、模拟退火方法 、基于各种经验势[5 ,6 ] 的共轭梯度法[7 ] 、 遗传算法[ 8 ] 等. Ni - Al 金属间化合物及其合金由于具有熔点高 、密度低 、抗氧化性好和热导率高等特点 ,受 到了人们广泛关注[9～11 ] . 本文中 ,笔者基于对 Ni - Al 化合物总能的密度泛函理论的计算结果 ,采用 Chen Mbius 反演[ 12 ] ,获得了 Ni - Al 的原子间的相互作用势 ;采用最陡下降法 、共轭梯度法和模拟退火法相结 合 ,对不同组分Ni - Al 的稳态结构进行了计算 ,为研究二元 混合团簇的稳态结构提供了一条可行的途径. 1 　研究方法 采用 GGA (gradient corrected LDA) 计算了 Fcc - Ni ,Fcc - Al 和 AlNi3 晶体的总能. 在计算中 ,使用第一原理正则守 ( ) 恒赝势 normconserving potential , 晶体波函数 由平面波展 开 , 自洽迭代过程中 ,使用 Monkhor stPack 方法在简约布里 渊区产生适当的 k 点以满足自洽迭代的精度要求 , 自洽迭代 的收敛判据为 5 ×10 - 7 eV/ atom . 首先用 CA STEP (cambridge sequential tot al energy p ackage) 计算得到总能 ,通过平移总能 曲线适应结合能的实验值得到结合能曲线 ,计算结果见图 1 , 表 1. 　图 1 　晶体结合能随晶格常数的变化 同种原子构成的晶体的结合能可表示为原子间势函数的 　　　表 1 　理论计算值与实验值 eV/ atom 求和