纳米线热导率的分子动力学研讨.pdfVIP

中国工程热物理学会 传热传质学 2006年学术会议论文 编号：063339 纳米线热导率的分子动力学研究 顾骁坤，王威扬，刘靳，曹炳阳 (清华大学航天航空学院，传热强化与过程节能教育部重点实验室。北京100084) 摘要：采用非平衡分子动力学方法研究了氩纳米线的热导率，模拟得到了相同截面积的圆形、方形以 及三角形截面纳米线的热导率。‘模拟发现纳米线的热导率小于体材料的值，体现}【{明显的尺寸效应， 并且在较低温时尺寸效应更显著。纳米线截面积相等时截面形状对纳米线的热导率有非常重要的影响， 揭示出面体比足表征纳米线热导率的关键因素，纳米线的热导率随着面体比的增大而减小。 关键词； 热导率纳米线尺寸效应分子动力学 1引言 近几年来，微电子、光电子、生物芯片实验室和微纳电子机械系统等纳米科学和技 术获得了飞速发展，人们已经能够在纳米尺度构建工程器件。纳米线是构建微纳器件的 基本构件之一，其导热性质影响着器件性能的可靠性和稳定性，了解纳米线的导热性质 对纳米技术中的热设计和热管理具有重要的意义[I,21。通过实验测定纳米尺度材料的热导 率，不仅对设备要求高，同时测量技术上也有很大的难度，目前实验能够测量的最小尺 度还不到100nm。分子动力学模拟是研究微尺度热传导特性的有力工具，特别是对于一 些理论上难以说明或实验中难以观察的现象给出原子水平的微观解释。目前的研究发现 纳米线的热导率比相应体材料的值要小，其尺寸效应归结为由于尺寸减小而引起的尺寸 效应，本质上源自纳米线表面对载热粒子的散射作用增强p制。 在芯片业，不同的工艺往往可以制备得到不同截面形状的纳米线，但是，目前纳米 线截面形状对其热导率的影响还鲜有报道。为了研究形状等囚素对纳米线热导率的影 响，本文采用氩作为研究对象，针对不同形状和不同温度下的氩纳米线采用非平衡分子 动力学算法方案(RNEMD)研究其热传导率并对影响其的因素做探讨。 2分子动力学模拟 本文采用拟扰动非平衡分子动力学算 J：将 法(RNEMD)，具体的实现如图1所示【7 系统沿X方向划分为N个断面(Slab)，将 中间的断面定义为热源，最左边的定义为 冷源。通过交换热源中动能最小的原子与 冷源中动能最大的原子的速度向量，来产 生热流。当这个过程进行了一段时间后系 统达到平衡，此后通过速度向苗交换产生 图1分子动力学模拟方法示意图 的热流会被材料的热传导所平衡。根据经典的Fourier定律可以得到热导率： ∑詈(砰一曙) 五=一—— 2tSaT}ax 其中m为原子质量；下标h，c分别指交换速度的原子中动能大的和动能小的原子：f指 的是统计热流的时间；S是纳米线的截面积；aT／igx是x方向上的平均温度梯度。 选用有可靠实验数据和原子间作用势的氩晶体为研究对象。原子依照氩的真实晶体 结构按面心立方(FCC)品格排列。本文所计算的系统的尺度，由于计算不同形状的纳 。 ， 44000个氩原子。为了讨论形状对热导率的影响，模拟中将纳米线的截面积取为定值 【，(，．)=4“(=)”一(二)6】 (1) ，． ，． 其中，占为能量参数，盯为原子直径。厂为原子中心之间的距离。 表1氪的LJ势能的参数 系统的温度通过阻尼力约束方法控制在确定值‘81。．计算分子间作用势能