平衡态分子动力学Green-Kubo方法计算氮化硼单层结构热导率的模型尺寸效应研究  

Effect of Domain Size on Thermal Conductivity of Boron Nitride Monolayer in Equilibrium Molecular Dynamics with Green-Kubo Method

作　　者：鲁济豹 杨楠楠 孙蓉[1] 汪正平[1,3] LU Jibao1 YANG Nannan1,2 SUN Rong1 WONG ChingPing1,3 1(Guangdong Provincial Key Laboratory of Materials for High Density Electronic Packaging, Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences, Shenzhen 518055, China ） 2（ College of Nanoscience and Technology of University of Science and Technology of China, Suzhou 215123, China ） 3（ Department of Electronics Engineering, The Chinese University of Hong Kong, Hong Kong 999077, China)

机构地区：[1]中国科学院深圳先进技术研究院广东省高密度电子封装关键材料重点实验室,深圳518055 [2]中国科学技术大学纳米科学技术学院,苏州215123 [3]香港中文大学电子工程系,香港999077

出　　处：《集成技术》2018年第2期12-21,共10页Journal of Integration Technology

基　　金：国家科技重点研发计划(2017YFB0406000);先进电子封装材料国地联合实验室(深圳市发展和改革委员会[2017-934]);中国科学院前沿科学重点研究项目(QYZDY-SSW-JSC010);广东省重点实验室(2014B030301014);深圳市基础研究学科布局(JCYJ20160331191741738)

摘　　要：分子动力学模拟可以直接表征体系原子的行为,因此成为研究氮化硼相关材料微观导热机理的重要工具,但目前尚没有关于氮化硼材料模型尺寸对其热传导相关性质影响规律的研究。该文采用平衡态分子动力学并结合Green-Kubo方法,研究了纯净氮化硼单层结构热导率、声子色散关系以及态密度随模拟尺寸的变化规律,并解释了其内部机理。实验发现,氮化硼单层材料热导率随着模拟尺寸的增大而减小,并在单层面积约4.1 nm×4.1 nm时收敛于(349±22)W/(m·K),此收敛值远小于平衡态分子动力学计算中石墨烯热导率的收敛尺寸(10 nm×10 nm),这说明氮化硼单层中声子之间的散射大于石墨烯。此外,不同于热导率,氮化硼单层结构的声子色散曲线、态密度几乎不受模拟尺寸的影响。该研究结果可为采用平衡态分子动力学研究氮化硼相关材料的微观导热机理提供重要参考。Molecular dynamics simulation can directly model the molecular behavior, making it convenient to investigate the microscopic mechanism of thermal conductance of boron nitride （BN） nanomaterials. However, there still no explicitly investigation to the size effect of the BN monolayer on its thermal properties, in this paper, the equilibrium molecular dynamics （EMD） combined with the Green-Kubo method was used to unravel and explain the relations between the system size and the thermal conductivity, phonon dispersion, and phonon density of states （DOS） of pristine BN monolayer. It was found that the thermal conductivity of the BN monolayer decreased with increasing the size of the sturcture, reaching a converged value （349±22） W/（m.K） at 4.1 nm×4.1 nm. This value was much smaller than the converging size （10 nm×10 nm） of graphene in calculating its thermal conductivity using EMD simulations, which implied the phonon-phonon scattering in BN monolayer was larger than that in graphene. Different from the thermal conductivity, the phonon dispersion and phonon DOS of the BN monolayer did not depend on the size of the structure. Our findings provide important reference for investigating the microscopic mechanisms of BN related materials by using the equilibrium molecular dynamics modeling.

关 键 词：氮化硼 平衡态分子动力学 Green-Kubo方法 热导率 

分 类 号：O736[理学—晶体学]