Two-fluid nature of phonon heat conduction in a monatomic lattice

The thermal resistance of a crystal lattice with a monatomic unit cell due to three-phonon scattering processes is investigated in detail theoretically. A general expression for the lattice thermal conductivity is derived from a combined analysis based on: (i) the Boltzmann equation and (ii) data on the heat current autocorrelation function obtained via molecular dynamics simulations in conjunction with the Green–Kubo formalism. It is argued that the phonon gas in a monatomic lattice conducts heat as if it consisted of two distinct parts (two ‘thermal fluids’), so that the lattice thermal conductivity can be decomposed into contributions from these two parts. The origin of the behaviour of the phonon gas, which is explored in the present work, is due to an intrinsic interplay between Umklapp and normal three-phonon scattering processes. New insight into the nature of the lattice thermal conductivity is demonstrated and the results of the present work are in agreement with previous studies in this area.     

Phonon relaxation and heat conduction in one-dimensional Fermiben Pastaben Ulam β lattices by molecular dynamics simulations

The phonon relaxation and heat conduction in one-dimensional Fermi-Pasta-Ulam (FPU) β lattices are studied by using molecular dynamics simulations. The phonon relaxation rate, which dominates the length dependence of the FPU β lattice, is first calculated from the energy autocorrelation function for different modes at various temperatures through equilibrium molecular dynamics simulations. We find that the relaxation rate as a function of wave number k is proportional to k^1.688, which leads to a N^0.41 divergence of the thermal conductivity in the framework of Green-Kubo relation. This is also in good agreement with the data obtained by non-equilibrium molecular dynamics simulations which estimate the length dependence exponent of the thermal conductivity as 0.415. Our results confirm the N^2/5 divergence in one-dimensional FPU β lattices. The effects of the heat flux on the thermal conductivity are also studied by imposing different temperature differences on the two ends of the lattices. We find that the thermal conductivity is insensitive to the heat flux under our simulation conditions. It implies that the linear response theory is applicable towards the heat conduction in one-dimensional FPU β lattices.     

Phonon-mediated heat dissipation in a monatomic lattice: case study on Ni

The recently introduced analytical model for the heat current autocorrelation function of a crystal with a monatomic lattice [Evteev et al., Phil. Mag. 94 (2014) p. 731 and 94 (2014) p. 3992] is employed in conjunction with the Green–Kubo formalism to investigate in detail the results of an equilibrium molecular dynamics calculations of the temperature dependence of the lattice thermal conductivity and phonon dynamics in f.c.c. Ni. Only the contribution to the lattice thermal conductivity determined by the phonon–phonon scattering processes is considered, while the contribution due to phonon–electron scattering processes is intentionally ignored. Nonetheless, during comparison of our data with experiment an estimation of the second contribution is made. Furthermore, by comparing the results obtained for f.c.c. Ni model to those for other models of elemental crystals with the f.c.c. lattice, we give an estimation of the scaling relations of the lattice thermal conductivity with other lattice properties such as the coefficient of thermal expansion and the bulk modulus. Moreover, within the framework of linear response theory and the fluctuation-dissipation theorem, we extend our analysis in this paper into the frequency domain to predict the power spectra of equilibrium fluctuations associated with the phonon-mediated heat dissipation in a monatomic lattice. The practical importance of the analytical treatment lies in the fact that it has the potential to be used in the future to efficiently decode the generic information on the lattice thermal conductivity and phonon dynamics from a power spectrum of the acoustic excitations in a monatomic crystal measured by a spectroscopic technique in the frequency range of about 1–20 THz.     

Insight into lattice thermal impedance via equilibrium molecular dynamics: case study on Al

Using results of equilibrium molecular dynamics simulation in conjunction with the Green–Kubo formalism, we present a general treatment of thermal impedance of a crystal lattice with a monatomic unit cell. The treatment is based on an analytical expression for the heat current autocorrelation function which reveals, in a monatomic lattice, an energy gap between the origin of the phonon states and the beginning of the energy spectrum of the so-called acoustic short-range phonon modes. Although, we consider here the f.c.c. Al model as a case example, the analytical expression is shown to be consistent for different models of elemental f.c.c. crystals over a wide temperature range. Furthermore, we predict a frequency ‘window’ where the thermal waves can be generated in a monatomic lattice by an external periodic temperature perturbation.     

石墨烯的声子热学性质研究

声子是石墨烯导热过程中的主要载体,而声子的弛豫时间又是其中最基本、最重要的物理量.本文采用简正模式分解法研究了石墨烯声子的弛豫时间,并且借此分析了不同声子在导热过程中的贡献.该方法通过平衡分子动力学模拟实现,首先通过模拟得到单个声子的能量自相关函数衰减曲线,并进一步采用拟合和积分两种方法得到单个声子的弛豫时间.然后,研究了弛豫时间与波矢、频率和温度的关系.结果发现,弛豫时间随波矢的变化与对应的色散关系相近,弛豫时间与频率和温度的关系符合理论模型:1/τ=νnTm,其中声学支的n为1.56,而光学支结果较为发散,指数m对于不同声子支结果略有不同.最后,还研究了不同频率声子对导热的贡献,发现低频声子在态密度上占有绝对优势,并且其弛豫时间整体高于高频声子,所以低频声子对导热的贡献占据主导地位.     

基于量子修正的石墨烯纳米带热导率分子动力学表征方法

本文提出了基于量子修正的非平衡态分子动力学模型,可用于石墨烯纳米带热导率的表征.利用该模型对不同温度下,不同手性及宽度的石墨烯纳米带热导率进行了研究,结果发现:相较于经典分子动力学模型给出的热导率随温度升高而单调下降的结论,在低于Debye温度的情况下,量子修正模型的计算结果出现了反常现象.本文研究还发现,石墨烯纳米带的热导率呈现出明显的边缘效应及尺度效应:锯齿型石墨烯纳米带的热导率明显高于扶手椅型石墨烯纳米带;全温段的热导率及热导率在低温段随温度变化的斜率均随宽度的增加而增大.最后,文章用Boltzmann声子散射理论对低温段的温度效应及尺度效应进行了阐释,其理论分析结果说明文章所建模型适合在全温段范围内对不同宽度和不同手性的热导率进行精确计算,可为石墨烯纳米带在传热散热领域的应用提供理论计算和分析依据.     

How to improve the accuracy of equilibrium molecular dynamics for computation of thermal conductivity?

Equilibrium molecular dynamics (EMD) simulations through Green-Kubo formula (GKF) have been widely used in the study of thermal conductivity of various materials. However, there exist controversial simulation results which have huge discrepancies with experimental ones in literatures. In this Letter, we demonstrate that the fluctuation in calculated thermal conductivity is due to the uncertainty in determination of the truncation time, which is related to the ensemble and size dependent phonon relaxation time. We thus propose a new scheme in the direct integration of heat current autocorrelation function (HCACF) and a nonzero correction in the double-exponential-fitting of HCACF to describe correctly the contribution to thermal conductivity from low frequency phonons. By using crystalline Silicon (Si) and Germanium (Ge) as examples, we demonstrate that our method can give rise to the values of thermal conductivity in an excellent agreement with experimental ones.     

Effects of fullerene coalescence on the thermal conductivity of carbon nanopeapods

The heat conduction and its dependence on fullerene coalescence in carbon nanopeapods (CNPs) have been investigated by equilibrium molecular dynamics simulations. The effects of fullerene coalescence on the thermal conductivity of CNPs were discussed under different temperatures. It is shown that the thermal conductivity of the CNPs decreases with the coalescence of encapsulated fullerene molecules. The thermal transmission mechanism of the effect of fullerene coalescence was analysed by the mass transfer contribution, the relative contributions of phonon oscillation frequencies to total heat current and the phonon vibrational density of states (VDOS). The mass transfer in CNPs is mainly attributed to the motion of encapsulated fullerene molecule and it gets more restricted with the coalescence of the fullerene. It shows that the low-frequency phonon modes below 20 THz contribute mostly to thermal conductivity in CNPs. The analysis of VDOS demonstrates that the dominating contribution to heat transfer is from the inner fullerene chain. With the coalescence of fullerene, the interfacial heat transfer between the CNT and fullerene chain is strengthened; however, the heat conduction of the fullerene chain decreases more rapidly at the same time.     

Effect of electrostatics techniques on the estimation of thermal conductivity via equilibrium molecular dynamics simulation: application to methane hydrate

Equilibrium molecular dynamics (MD) simulations for three system sizes of fully occupied methane hydrate have been performed at around 265 K to estimate the thermal conductivity using the Ewald, Lekner, reaction field, shifted-force and undamped Fennell–Gezelter methods. The TIP4P water model was used in conjunction with a fully atomistic methane potential with which it had been parameterized from quantum simulation. The thermal conductivity was evaluated by integration of the heat flux autocorrelation function (ACF) derived from the Green–Kubo formalism; this approach vas validated by estimation of the average phonon mean free path. The thermal conductivities predicted by non-periodic techniques were in reasonable agreement with the experimental results of 0.62 and 0.68 W/m K, although it was found that the estimates by the non-periodic techniques were up to 25% larger than those of Lekner and Ewald estimates, particularly for larger systems. The results for the Lekner method exhibited the least variation with respect to system size. A decomposition of the heat flux vector into its respective contributions revealed the importance of electrostatic interactions, and how different electrostatic treatments affect the contribution to the thermal conductivity.     

离子辐照对磷烯热导率的影响及其机制分析

通过离子辐照产生缺陷,可以非常有效地调控磷烯诸多物理性质.本文应用分子动力学方法模拟离子辐照磷烯的过程,给出了缺陷的形成概率与入射离子能量、离子种类以及离子入射角度之间的关系,并且应用非平衡态分子动力学计算辐照后磷烯热导率的变化.以缺陷形成概率为切入点,系统地研究了辐照离子的能量、辐照剂量、离子的种类以及离子的入射角度对磷烯热导率的影响.应用晶格动力学方法研究了空位缺陷对磷烯声子参与率的影响,并计算了声子局域模式的空间分布.基于量子微扰和键弛豫理论,指出空位缺陷明显降低磷烯热导率的最重要物理机制是空位缺陷附近的低配位原子对声子强烈散射.本文研究可为缺陷工程调控磷烯的热输运性质提供理论参考.     

Thermal conductivity of carbon nanoring linked graphene sheets:A molecular dynamics investigation

Improving the thermal conduction across graphene sheets is of great importance for their applications in thermal management. In this paper, thermal transport across a hybrid structure formed by two graphene nanoribbons and carbon nanorings(CNRs) was investigated by molecular dynamics simulations. The effects of linker diameter, number, and height on thermal conductivity of the CNRs–graphene hybrid structures were studied respectively, and the CNRs were found effective in transmitting the phonon modes of GNRs. The hybrid structure with 2 linkers showed the highest thermal conductivity of 68.8 W·m~(-1)·K~(-1). Our work presents important insight into fundamental principles governing the thermal conduction across CNR junctions and provides useful guideline for designing CNR–graphene structure with superior thermal conductivity.     

三氨基三硝基苯基高聚物粘结炸药热力学性质的理论计算研究

高聚物粘结炸药(PBX)的热力学性质是用于炸药结构响应、安全性评估、数值模拟分析等的重要参数.由于PBX结构的多尺度特性,完全采取实验方法精细表征这些参数存在巨大的挑战.本文运用第一性原理和分子动力学计算的方法,系统研究了三氨基三硝基苯(TATB)基高聚物粘结炸药的热力学参数和界面热传导性质.利用散射失配模型研究了TATB与聚偏二氟乙烯(PVDF)界面的热传导过程,发现热导率随温度升高而上升,并且在高温情况下接近于定值.基于分子动力学获得的TATB热导率并结合界面热导率,分析了PBX炸药的热导与颗粒尺寸的关系,当颗粒尺寸大于100 nm时,界面热阻对于PBX热导率的影响有限.     

A molecular dynamics simulation: the effect of finite size on the thermal conductivity in a single crystal silicon

Non-equilibrium molecular dynamics (NEMD) simulations are performed to calculate thermal conductivity. The environment-dependent interatomic potential (EDIP) potential on crystal silicon is adopted as a model system. The issues are related to nonlinear response, local thermal equilibrium and statistical averaging. The simulation results by non-equilibrium molecular dynamics show that the calculated thermal conductivity decreases almost linearly as the film thickness reduced at the nanometre scale. The effect of size on the thermal conductivity is also obtained by a theoretic analysis of the kinetic theory and formulas of the heat capacity. The analysis reveals that the contributions of phonon mean free path (MFP) and phonon number in a finite cell to thermal conductivity are very important.     

导热系数的分子动力学模拟研究及相关问题的探讨

对于采用分子动力学方法研究导热系数的背景、研究现状及存在的问题进行了综述和分析。总结了通过分子动力学模拟方法求得导热系数的物理模型和基本算法。讨论了目前在微尺度导热问题的研究中引入分子动力学模拟方法需要考察的影响因素和几个重要问题。     

带孔硅纳米薄膜热整流及声子散射特性研究

本文基于非平衡的分子动力学模拟方法计算了带有三角形孔的硅纳米薄膜的界面热阻特性,结果表明300-1100 K范围内随着热流方向的改变,在含有三角形孔的硅纳米薄膜中存在热整流效应,热整流系数达28％.同时借助于声子波包动力学模拟方法,获得了不同频率下的纵波声子在三角形孔处的散射特性,结果表明纵波声子在散射过程中产生了横波声子,并且从三角形底部向顶部入射的声子能量透射系数比反向时平均低22％.不对称结构引起的声子透射率的差异是引起热整流效应的主要因素.     

多约束纳米结构的声子热导率模型研究

纳米技术的快速发展使得对微纳尺度导热机理的深入研究变得至关重要. 理论和实验都表明, 在纳米尺度下声子热导率将表现出尺寸效应. 基于声子玻尔兹曼方程和修正声子平均自由程的方法得到了多约束纳米结构的声子热导率模型, 可以描述多个几何约束共同作用下热导率的尺寸效应. 不同几何约束对声子输运的限制作用可以分开计算, 总体影响则通过马西森定则进行耦合. 对于热流方向的约束, 采用扩散近似的方法求解声子玻尔兹曼方程; 对于侧面边界约束, 采用修正平均自由程的方法计算边界散射对热导率的影响. 得到的模型能够预测纳米薄膜(法向和面向)及有限长度方形纳米线的热导率随相应特征尺寸的变化. 与蒙特卡罗模拟及硅纳米结构热导率实验值的对比验证了模型的正确性.     

硅纳米薄膜法向热导率的进一步分析

应用非平衡分子动力学方法进一步研究了平均温度为300K、厚度为2.715nm-43.44nm的单晶硅薄膜的法向热导率,模拟结果表明,薄膜热导率低于同温度下单晶硅的实验值,存在显著的尺寸效应,当膜厚度在20nm以下时,法向热导率随尺度减小而线性减小,当膜厚度大于20nm时法向热导率随尺度呈现二阶多项式变化。法向热导率的变化规律与面向热导率的变化规律类似,表明薄膜厚度和表面晶格结构对声子传热影响的复杂性。     

孔洞缺陷对二维石墨烯/氮化硼横向异质结热传导的调控

本文采用孔洞缺陷来实现对二维石墨烯/氮化硼横向异质结热导率的调控.平衡态分子动力学(EMD)计算结果表明,界面孔洞的引入会降低二维石墨烯/氮化硼横向异质结的热导率.相较于有序的孔洞分布,无序的孔洞分布能够更有效地降低异质结的热导率,这一现象可通过声子安德森局域化来解释.孔洞缺陷的存在导致声子的频率和波失发生变化,从而使声子散射变得更加频繁,孔洞随机分布时,则导致声子波在材料中发生多次反射和散射,最终形成局域振动模式.本研究揭示了孔洞缺陷降低二维石墨烯/氮化硼横向异质结热导率的物理机制,对二维热电材料的结构设计有一定的指导意义.     

碳纳米管的热导率:从弹道到扩散输运

采用非平衡分子动力学方法研究了300 K和1000 K时（5,5）碳纳米管热导率随长度的变化.在室温下,碳纳米管长度小于40 nm时热导率与长度呈线性关系,此时导热处于弹道输运阶段,单位面积弹道热导为5.88×10⁹ Wm^-2K^-1．随着碳纳米管长度的增加,其热导率逐渐增加,但增加速度随长度逐渐减小,此时导热处于弹道—扩散输运阶段,并随长度的增加从以弹道输运为主向以扩散输运为主转变．长度大于10 μm时由于弹道输运可以忽略,导热近似达到完全 关键词： 碳纳米管 热导率 弹道输运 低维导热     

Ballistic–diffusive phonon transport and size induced anisotropy of thermal conductivity of silicon nanofilms

The effective thermal conductivity of nanofilms is size dependent due to the diffusive–ballistic transport of phonons. In this paper, we investigate the cross-plane phonon transport from the viewpoint of the phonon Boltzmann equation. A predictive model for the size dependent thermal conductivity is proposed and agrees well with the results of molecular dynamics simulation for silicon nanofilms. The ballistic transport has different effects on the heat conduction in the in-plane or cross-plane directions, which causes the anisotropy of thermal conductivity of nanofilms. Such anisotropy is also size dependent and vanishes with the increase of film thickness.