金属钒及其氢化物晶体结构模拟.pdfVIP

　第34卷第5期 原 子 能 科 学 技 术 Vol. 34 ,No. 5 　2000年9月 Atomic Energy Science and Technology Sep. 2000 金属钒及其氢化物的晶体结构模拟 彭述明 ,赵鹏骥 ,杨茂年 ,姚书久 ,罗顺忠 ,徐志磊 ( 中国工程物理研究院 核物理与化学研究所, 四川 成都　610003) ( ) 摘要 :用 Castep 程序的局域密度泛函近似方法 LDA 对金属钒及不同原子比的钒氢化物进行了结 构优化和总能计算。给出了 VHx (0 ≤x ≤2) 体系的各种氢化物的晶体结构的计算结果 ;分析了氢 含量对 VH 体系相组成的影响。与现有的实验结果相比较 ,二者基本相符。 关键词 :钒 ;钒氢化物 ; 晶体结构 中图分类号: TG11531 　　　文献标识码 :A 　　　文章编号 :(2000) 钒基贮氢合金作为一种先进的贮氢材料在国防和民用工业中具有广阔的应用前景。 N Budylkin 等[1 ]研究了氢同位素在 VCrTi 合金中的行为 ,该材料成功地应用于聚变堆中处 理氢同位素 ;VL Arbuzov 等[2 ]研究了 H 在 V 中的“陷井”效应; G. G. Libowitz 等[3 ] 研究了立 ( ) ( ) [4 ] 方体心 bcc 的金属钒及钒基合金 以VTi 合金为主 的吸、放氢行为;ON Stephen 研究了 ( ) [5 ] VNb 、VCr 等 bcc 合金固溶氢的行为 ; HD Carstanjen 等 测定了 VD 、VT 体系的相图 ; A J Maeland 等[6 ,7 ]证明了 V 基 bcc 合金在室温下无需活化便能与氢快速反应。目前 ,准确实 验确定 H/ V 原子比、H 在金属 V 中的位置和存在形式有一定难度 ,对 H/ V 原子比为非整数 氢化物的晶体结构的测定尤为困难。因此 ,对贮氢材料的研究主要集中在经验配方和性能的 实验研究上 ,对贮氢机制和材料设计的研究报道所见甚少。本工作用 Castep 程序的局域密度 ( ) 泛函近似方法 LDA 对金属钒及不同原子比的金属钒氢化物进行结构优化和总能计算。 1 　计算方法 2 ( ) 计算在 INDIGO 和 IND Y 两台工作站上完成。利用局域密度泛函近似 LDA 与赝势技 ( ) 巧相结合的方法 , 采用 BF GS 算法 NewRaphson 算法的 BF GS 版本 , 动能域值设为 1360 eV ,设置 Subspace rotation of wave functions 以加速总能收敛,整个体系作为金属来考 虑。总能收敛指标为 2 ×10 - 4 eV ,位移均方差设为 001 nm ,力场的均方差小于 50 GeV/ m , 对不同含氢量的 V H x 体系的晶体结构进行全优化。在建立 V H x 体系的初始构型时 , 当 x ≤ 10 时 ,金属原子的格点位置与金属钒的原子位置相同 ,H 填在四面体间隙中 ;当 1 ≤x