Davier B, Larroque J, Dollfus P, Chaput L, Volz S, Lacroix D, Saint-Martin J
C2N UMR 9001, Univ. Paris-Sud-CNRS, Université Paris-Saclay, 91405 Orsay, France.
J Phys Condens Matter. 2018 Dec 12;30(49):495902. doi: 10.1088/1361-648X/aaea4f. Epub 2018 Nov 15.
The Boltzmann transport equation is one of the most relevant frameworks to study the heat transport at the nanoscale, beyond the diffusive regime and up to the micrometer-scale. In the general case of 3D devices, the particle Monte Carlo approach of phonon transport is particularly powerful and convenient, and requires reasonable computational resources. In this work, we propose an original and versatile particle Monte Carlo approach parametrized by using ab initio data. Both the phonon dispersion and the phonon-phonon scattering rates have been computed by DFT calculation in the entire 3D Brillouin zone. To treat the phonon transport at rough interfaces, a combination of specular and diffuse reflections has been implemented in phase space. Thermal transport has been investigated in nanowires and thin films made of cubic and hexagonal Silicon, including edge roughness, in terms of effective thermal conductivity, phonon band contributions and heat flux orientation. It is shown that the effective thermal conductivity in quasi-ballistic regime obtained from our Monte Carlo simulation cannot be accurately fitted by simple semi-analytical Matthiessen-like models and that spectral approaches are mandatory to get good results. Our full band approach shows that some phonon branches exhibiting a negative group velocity in some parts of the Brillouin zone may contribute negatively to the total thermal flux. Besides, the thermal flux clearly appears to be oriented along directions of high density of states. The resulting anisotropy of the heat flux is discussed together with the influence of rough interfaces.
玻尔兹曼输运方程是研究纳米尺度热输运最相关的框架之一,适用于扩散区域之外直至微米尺度的情况。在三维器件的一般情况下,声子输运的粒子蒙特卡罗方法特别强大且方便,并且需要合理的计算资源。在这项工作中,我们提出了一种原始且通用的粒子蒙特卡罗方法,该方法通过使用从头算数据进行参数化。声子色散和声子 - 声子散射率均已通过在整个三维布里渊区的密度泛函理论(DFT)计算得出。为了处理粗糙界面处的声子输运,在相空间中实现了镜面反射和漫反射的组合。我们研究了由立方和六方硅制成的纳米线和薄膜中的热输运,包括边缘粗糙度,涉及有效热导率、声子能带贡献和热流方向。结果表明,通过我们的蒙特卡罗模拟获得的准弹道区域中的有效热导率不能通过简单的类似马蒂森的半解析模型准确拟合,并且光谱方法对于获得良好结果是必不可少的。我们的全带方法表明,在布里渊区某些部分表现出负群速度的一些声子分支可能对总热流有负贡献。此外,热流明显沿着态密度高的方向取向。我们讨论了由此产生的热流各向异性以及粗糙界面的影响。