Ni合金结构分子动力学模拟.PDFVIP

第 15卷 第4期 化 学 物 理 学 报 V.1.15. No.4 2002年8月 CHINESEJOURNALOFCHEMICALPHYSICS Aug.2002 1003-7713/2002/04-288-7 Cu-Ni合金结构的分子动力学模拟* 丛红日”，‘边秀房，‘李 辉b，王 丽 “ (a.山东大学材料液态结构及其遗传性教育部重点实验室，济南 2500( b南京大学物理系固体微结构国家重点实验室，南京 210093) 摘 要:米用EAM作用势对Cu-N，合金的结构特性进行了MD模拟研究通过Fz结构因子可发现，Cu含舌的 变化对结构因子的波动影响很小，健取向序参数和健对也表现出相似的变化规律，这表明液态Cu-Ni合金对成份 变化不敏感，体系中的化学序较弱一将Cu,Ni、合金熔体的FZ结构因子与W-d.的实验结果进行对比，发现二 者吻合得较好，表明EAM势可以很好地描绘Cu-Ni合金的结构特性.在快速冷却过程中，除了Cu,,Ni‘合金外，其 他合金成份的双体分布函数的第二峰都发生了臂裂，标志着体系最终形成了非晶结构，而Cu,,Ni。合金的双体分 布函数却表现出晶体峰的特征 通过对健取向序参数、健型指数以及钥镶原子的有效扩散系数的分析表明，在快 速冷却过程中，Cu,,Ni。合金最终形成了hep晶体结构 关键词:Cu-Ni合金;MD模拟;结构因子;中程有序 中图价敬昙 .()56卜4 女献坛识码 .A J 引 言 2 镶嵌原子势 分子动力学模拟的结果主要取决于原子间相互 按照EAM势[’」的基本表述，对任意原子构型， 作用势的选择17-3]，而过渡金属中强烈的多体效应 系统的总能量可描述为: 使得偶势描述往往不能给出合理的结果.80年代中 E,.,=},F,p(})·奋菩，‘r。，‘ ，，‘ 期，Daw和Raskes根据准原子概念和密度泛函理 了俩少 论、，‘·，提出了一个能较好描述过渡金属各种性质的 式中，F,(p,)为镶嵌一个原子到电子密度为p的局 新方法— 镶嵌原子法(EmbeddedAtomMethod，简 域环境中的能量;p为i处由周围贡献的电r密度 称EAM).后经Johnson的改进16]，EAM法在固态材 的总和，可表示为: 料的诸多领域已经取得了很大成功[-1.Fodes计 ，二艺石(r,,) (2) 算了Cu-Ni合金表面的偏析il「l，结果表明，在近表 少尹 ‘ 原子的电子密度f(r)由数值求解Hartree-Fock波函 面处，合金的成分与表面层是一种非线性的关系.李 数得到.1988年，Johnson开发了用于fee结构的过 辉等研究了快速冷却过程中Cu,,Ni,。合金微观结构 渡金属EAM势161，其结构采用衰减指数的形式为: 的变化 ，」‘，指出体系的微观构型受温度和冷却速率 的影响.而对于整个Cu-Ni合金体系的微观结构研 “，·‘=一“·{，一(育一)]‘ 究却很少. 本工作选择Cu-Ni合金作为研究对象，分析了 xexp卜a(r一1)] (3) 不同成分的Cu-Ni合金在快速冷却过程中的结构特 其中，a二3(从叹/E,)`2,E为内聚能;E.,几分 证，发现除了Cu2oNi,。合金外，其他成分的Cu-Ni合 别为体积弹性模量和原子体积.Mei等根据