AuCu团簇分子动力学研究.PDF

第 38 卷　第 6 期 . 38　 . 6 ( ) V o l N o 厦门大学学报 自然科学版 　1999 年 11 月 ( ) N ov. 1999　 Jou rn a l o f X iam en U n iver sity N atu ra l Scien ce A uCu 团 簇 分 子 动 力 学 研 究 ① 谭　凯　林梦海　莫亦荣　张乾二 ( 厦门大学化学系, 物化所, 固体表面物理化学国家重点实验室　厦门　36 1005) 摘要　根据 的嵌入原子法( ) , 结合 双体势和 多体势, 对 团簇 D B EAM John son M u r re ll A uCu 结构进行了分子动力学模拟. 研究表明 EAM 势函数可较好地描述非晶态合金体系. 非晶态团簇形成 ( ) 中以多面体方式生长, 无序态 采用二十面体及十面体结构 比有序态更稳定. 关键词　 团簇, 非晶结构, 分子动力学, 模拟 A uCu 中国图书分类号　 64 1. 1 O 1　简　介 原子簇研究是化学中引人关注的领域. 它可以帮助人们认识原子团簇是怎样作为基本单元, 构造起微观原子与宏观材料间的桥梁. 尽管等离子体质谱实验能提供各种原子簇生成的信息, 但 ( 仍然需要理论计算去研究其结构和能量. 早期的理论计算假设原子簇是某种结构 例如fcc 和 b cc) 大块固体的碎片, 进一步研究发现原子簇以多面体方式生长, 非晶态的结构 占较为重要地 ( ) 位 即含有 5 次对称轴 . 我们应用分子动力学( ) 方法研究了金铜合金体系. 比较原子簇增大过程中非晶态与晶 M D 态结构稳定性差异. 在分子动力学模拟中, 描述原子间相互作用关系的势函数是整个模拟的关 键. 传统的两体势在研究材料的弹性系数性质和预言金属的结合能及空位形成能时难以获得准 确 的结果. 80 年代 以来, 各种经验或半经验 的多体模型得到迅速发展, 特别是嵌入原子法 (EAM ) , 即克服了两体势的缺陷, 又不会使计算量太大. 2　模拟计算方法 由D aw 和B a sk e s 提出的嵌入原子法(EAM ) , 把每一个原子看成是镶嵌到其他所有原子组 成的基体中, 嵌入能则依赖于电子密度, 引入一个杂质后, 总势能就是基体和杂质势能的泛函, 也 就是无杂质基体电子密度的泛函及杂质类型和位置的函数. 体系的总能量 　 E = F () + ( r ) 2 ( 1) tot ∑ i i ∑ ij ij 其中F 是嵌入能, 是在r