铜纳米线熔化过程的分子动力学模拟.pdfVIP

2007年 赣南师范学院学报 №．6 第六期 Journal of Gannan Normal University Dec．2007 铜纳米线熔化过程的分子动力学模拟’ 陈远铭，曾瑞淮 (赣南师范学院物质设备处，江西赣州 341000) 摘 要：本文采用分子动力学方法对111轴铜纳米线的熔化过程进行了计算机模拟研究，结果表明：铜纳 米线的熔点随体系的大小而变化，体系越小熔点越低，并与其横截面积的平方根的倒数近似成线性关系． 关键词：铜纳米线；熔化；分子动力学模拟 中图分类号：TP391．9 文献标识码：A 文章编号：1004—8332(2007)06—0036—04 近几年随着碳纳米管和硅纳米线等一维纳米材料的研究进一步深人，人们发现金属纳米线集合了一维 纳米材料及金属本身的特性，也具备优异的物理和化学性能，在宇航、电子、冶金、化工、生物和医学领域展现 出广阔的应用前景．因此，对金属纳米线的结构和性能的研究具有重要的理论和现实意义．一维纳米材料是 指两个维度处于纳米尺度(0．1 nm到100 ntn之间)的材料，如纳米线、纳米管等，金属纳米线是由金属原子 构成的实心的一维纳米材料．在目前条件下，通过实验来研究纳米线的微观变化过程比较困难，而近几年发 展起来的分子动力学模拟方法已成为研究微观物质性质的重要手段，现广泛用于凝聚态物理、纳米力学、材 料科学、核技术、化学反应动力学、生物学、药物设计等领域． 本文采用分子动力学方法，利用已应用于许多结构研究中的Ackland【1 镶嵌原子势(EAM)，对 111 轴铜纳米线的熔化过程、熔点等进行了模拟，并通过势能一温度曲线和热容一温度曲线对其熔化过程、熔点 及与体系尺寸的关系进行分析． 1 模拟方法与势函数 分子动力学模拟中原子间的相互作用以势函数描述，其准确性对计算结果的影响较大，在对所选定的系 统模拟之前，选定合适的势函数至关重要．本文采用的镶嵌原子势(EAM)，Ackland等LlI2 在此基础上通过 拟合金属的弹性常数、点阵常数、空位形成能、聚合能及压强和体积的关系，给出了Cu、Au、Ag、Ni的多体势 及对势函数，它们分别为： (r )=∑Ak(Rk—r ) ( 一r ) (1) (r )=∑ ( 一r ) 日( 一r0) (2) )={ 。0 (3) A ，R ， ， 为拟合的常数，且有R。R：，r。r：… r ，它们的取值随金属物质不同而有所不同． 本文模拟的纳米线的初始模型在 、y、z三个方向的晶胞数分别为28 18 52，三维尺寸分别为7．22，I 8．04 32．83 nm，晶体结构为金属cu的面心立方晶格，晶格常数a。=0．3645 nm，坐标系的x方向沿晶体的 方向，y方向沿[1l2]，z方向沿方向[111]．模拟时在 、y轴方向运用周期边界条件，这使得体系相当于一 个无穷长的二维纳米线阵列结构． 模拟中采用Gear 算法求解体系运动方程，时间步长为10 s．分析计算时，我们选取了六种不同尺寸 的纳米线，分别取名为Cul、Cu2、Cu3、Cu4、Cu5和Cu6，它们的晶胞数、原子数及尺寸见表1． · 收稿日期：2007—07—16 修回日期：2007一lO—l5 作者简介：陈远铭(1965一)．男．江鹳 l：犹人．赣南师范学院物质设备处高级实验师、硕士，主要从事材料物理实验教学和研究 第6期 陈远铭，曾瑞淮 铜纳米线熔化过程的分子动力学模拟 37 表1 原子数及其尺寸大小 2 模拟结果与分析 模拟时系统驰