纳米金球熔化的分子动力学研究.pdfVIP

第22卷 第 1期 河南教育学院学报 (自然科学版) V01．22 No．1 2013年 3月 JournalofHenanInstituteofEducation(NaturalScienceEdition) Mar．2O13 doi：10．3969／j．issn．1007—0834．2013．O1．003 纳米金球熔化 的分子动力学研究 王三军 ，王雪青2，饶凤飞2，宋友林’，李金铭 ，王 飞2 (1．河南教育学院 理论物理重点学科组，河南 郑州450046；2．郑州大学 物理工程学院，河南 郑州 450001) 摘要 ：利用Johnson的嵌入原子势(EAM)和分子动力学的方法对贵金属金 (Au)的有关性质进行 了模拟计算 ， 得到了Au在不同温度下的晶格常数和平均势能，以及Au在不 同温度下的径向分布函数 (RDF)，在此基础上得到 Au的熔点为 1350K左右，结果与实验值 1337K符合很好．进一步构建 了一个Au纳米球 ，研究给 出了纳米球 的平 均势能、热熔和RDF随温度变化关系。得到其在各个温度下的熔化模拟 图像．得到直径4nm左右含有 2ll2个原 子的Au纳米球的熔点为 900K，和体材料的熔点 l350K有较大的差别，最后分析并讨论 了纳米金球 的熔化过程 ． 关键词 ：分子动力学；晶格常数 ；熔点 中图分类号 ：0562．4 文献标识码 ：A 文章编号 ：1007—0834(2013)O1—0010—04 0 引言 金是位于元素周期表第 79位的金属 ，是一种贵金属 ，在几世纪 以来都被用作货币、保值物及珠宝．由于 金 的化学性质十分稳定 ，在 自然界中，经常出现在岩石中的金块或金粒、地下矿脉及冲积层中．金亦是货币金 属之一．在室温下为固体、密度高、柔软、光亮，其延展性及延性均是已知金属中最高的．纯金的亮黄色在传统 上被认为具有吸引力．而纳米状态的金却具有很高的催化活性 ，在催化剂领域应用极其广泛 ，此外在航天、生 物 医学等领域也有很大应用范围 I2．但是催化应用中的金主要为几十至上百个原子的金 团簇和纳米状态 金 ，这引起 了人们很大的研究兴趣．大量实验验证 ，当材料的尺寸减小到纳米量级时，纳米材料的熔融性质明 显低于其块体材料 ，因此金纳米团簇的熔化过程和熔点也和块体材料 的金有很大的区别 ，在应用上也是 需要考虑的一个重要 内容． 然而由于实验条件 、时间和空间分辨率的限制 ，对于纳米尺度的熔点和熔化过程很难进行细致研究．本 文利用分子动力学计算模拟的方法先对金块体材料的晶格常数和平均势能随温度变化关系及熔点、弹性系 数和切变常数 以及径向分布函数进行 了模拟计算 ，模拟计算得到的结果和实验符合得很好．然后构造一个 Au纳米球，模拟计算了其势能随温度变化关系，利用径向分布函数和热容性质计算给出Au纳米球的熔点， 还利用可视化程序研究了Au纳米球在各个温度下的熔化图像 ，分析讨论 了Au纳米球的熔化过程． 1 计算方法 本文中所有的计算都是采用 Johnson的Au的嵌入原子势(EAM)，这一势 的相关参数可以查询文献 ． 尽管这一势函数是半经验的，它还是很好地给出了Au的一些基本特性 ，比如弹性常数、体模量 、声子色散曲 线 以及单个空位和原子的扩散特性． 模拟采用周期性边界条件 ，计算 中时间步长采取 1．28×10 s．对块体材料的计算 中包含 4000个原 子 ，模拟温度从 100K至 1400K，每 100K取一次结果．100K模拟开始时原子位置采用 Au—FCC晶体结构 设置 ，晶格常数 0=4．08A，以后每一次模拟均 以上一次所产生的构型为基础进行模拟计算．计算采用 Verlet 算法积分原子的牛顿运动方程．计算体系选用 Ⅳ尸 系综(粒子数 Ⅳ、压力P和温度 不变)．采用 ANDERS． EN等压算法 计算了不同温度下Fe的晶格常数以便于在有限温度下模拟时用．在构建 Au纳米球的过程 中，首先建立一个大小为 15 0。×15 D。×15 n。(其中o。为晶格常数)的超原胞，利用命令构建一个半径 为5：．=o。的球 ，去掉 Au纳米球之外 的原子 ，得到该 Au纳米球中的原子数 目(2112个)和坐标．在以下的模拟 收稿 日期 ：2012一l1—12 基金项 目：河南省 自然科学基金资助项 目(112300410023) 作者简介：王三