氩晶体薄膜法向热导率的分子动力学模拟!-物理学报.PDF

第 卷 第 期 年 月 物 理 学 报 ## %$$( ， ， ， NC8 )## /C ) O61A6EP %$$( （ ） $$$5%3$B%$$(B## $ B$$$5$# GHIG JFKLMHG LM/MHG !%$$( H971 ) J9PQ ) LC+ ) 氩晶体薄膜法向热导率的分子动力学模拟! ! 吴国强 孔宪仁 孙兆伟 王亚辉 （哈尔滨工业大学卫星技术研究所，哈尔滨 #$$$） （ 年 月 日收到； 年 月 日收到修改稿） %$$# %’ %$$# ( ( 结合卫星“微型核”的特点，研究电介质薄膜中的导热机理以及薄膜厚度对导热系数的影响 以结构较为简单、 ) 具有可靠势能函数，实验数据较为丰富和可靠的氩的（ ）晶体为模型，采用平衡分子动力学方法（ ）和各向异 *++ ,-. 性非平衡分子动力学方法（ ）计算了氩晶体及其法向薄膜的热导率，并与实验结果进行比较 模拟结果表明， /,-. ) 氩晶体纳米薄膜的热导率显著小于对应大体积晶体的实验值，具有明显的尺寸效应 在氩薄膜厚度为 — ) % 0 %’ # 0$12 的模拟范围内，薄膜的法向热导率随着薄膜厚度的增加而呈近似线性增加) 关键词：热导率，纳米薄膜，尺寸效应，平衡分子动力学，非平衡分子动力学 ： !## $%$$ 算机技术的高速发展，在国内使用分子动力学的方 0 引 言 法从事微观尺度下材料的物理特性的研究逐渐增 多，并取得了一些采用分子动力学研究材料微尺度 ［— ］ 分子动力学计算机模拟是研究微尺度热传输特 特性的最新成果% 4 ) 从统计物理的角度可以将分 性的有力工具，特别是对于一些理论上难以说明或 子动力学分为平衡分子动力学模拟（ ）和非平衡 ,-. 实验中难以观察的现象给出基于物质微观结构及分 分子动力学模拟（ ）两种模拟方法 本文采用 /,-. ) 子动力学关系的解释，尤其适合研究薄膜材料的微