中介尺度Au纳米团簇熔化的分子动力学模拟.pdfVIP

物理化学学报 (WuliHuaxueXuebao) ActaPhys．一Chim．Sin．，2003。19(1)：35～39 35 中介尺度Au纳米团簇熔化的分子动力学模拟 张妍宁 王 丽 边秀房 (山东大学材料液态结构及其遗传性教育部重点实验室，济南 250061) 摘要 采用分子动力学模拟技术，研究了原子个数为 16—8628的Au纳米团簇的熔化过程 ．采用Johnson的 EAM(embeddedatommethod)模型，模拟结果表明，金属纳米团簇存在一中介尺度区域 ．对Au纳米团簇而言， 当原子个数 N456时，团簇的热力学性质与团簇尺寸呈线性关系，熔化首先从表面开始，逐步向中心区域推 进，且满足 一 (N)=nⅣ 的关系 ．另外，计算了中介区域的团簇的尺寸、熔化温度、表面能、熵、焓等热 力学量以及均方根位移(RMSD)等动力学量 ，为研究纳米团簇提供定量数据 ． 关键词： 分子动力学模拟， Au纳米团簇， 熔化过程 中图分类号： O641 纳米材料的出现引起了世界各国的广泛关注， 模拟分两部分进行，首先是对纳米团簇 的模 并相继开展了对这种新材料的研究工作 ．随着物 拟．我们选定不同半径的纳米球，其构型是完整的 质的超微化，纳米材料表面电子结构和晶体结构发 FCC结构，采用在无周期性边界条件下的等温一等 生变化，产生了宏观物体所不具有的四大效应，使 容(NVT)模拟方法 ． 得其具有传统材料所不具备的一系列优异的宏观 其次是对大块体系的模拟 ．我们以一个施加 特性 ．这也给材料科学带来了更多的课题，对纳米 周期性边界条件的立方盒子中的500个原子为研 材料微观结构的研究即是一个方面，因为纳米材料 究对象，采用等温 一等压 (NPT)模拟方法 ．同样 ，从 的晶粒尺寸变化和界面结构以及 晶粒 内部的变化 完整的FCC几何结构开始 ．模拟的时间步长取为 对其性能有直接的影响…． 1．28fS，首先在 300K下运行 5000个时间步长，使 中介尺度的纳米团簇通常指直径在 1～10nn]， 不同尺寸的原子团簇及大块体系处于平衡状态，然 介于分子和晶体之间的团簇 ．其性质可以认为是 后以1．56×10 K ·S 的速度持续升温，每隔2OK 随微粒数 目的增多而逐步由分子向宏观晶体过渡， 记录体系的能量和构型，直至2000K．进而获得体 但从实验上得出中介尺度的范围目前还比较困难， 系的热力学参数和动力学参数 ． 这种体系的热力学的数据也比较少 ．近年来发展 的分子动力学模拟无疑成为了解纳米团簇性质的 2 热力学计算 理想工具 。．另外，纳米团簇的熔化温度随原子尺 2．1 热容的计算 寸的增加而升高已得到人们的普遍公认，但我们对 通过对加热过程 中得到的能量 曲线求导来计 这一 问题只是定性的分析而缺乏定量的数据 ．本 算体系的热容，其计算公式为 文采用分子动力学模拟方法，利用Johnson的EAM Cp(T)=d(Ep)／d丁+3R／2 (1) (embeddedatommethod)模型 ，计算了中介尺度 其中R为气体常数 ， 为体系的能量 ．我们定义 Au纳米团簇熔化中的热力学量，如熔点、表面能、 熔点为热容取得最大值时所对应的温度 ． 熵及焓等的变化，并通过计算均方根位移来描述体 2．2 团簇尺寸的计算 系熔化过程中团簇 的变化规律，并将原子位置可视 原子团簇的半径 R。可通过下式计算： 化 ． R