基于声子输运理解非平衡分子动力学

分子动力学是微纳尺度传热研究中常用的计算方法。因为其使用简单、容易展示，在微纳尺度传热研究中起着至关重要的作用。分子动力学计算热导率的主流方法有两种，一种是平衡态分子动力学：运用涨落耗散的Green-Kubo公式计算，后面简称EMD。第二种是非平衡分子动力学方法：通过引入热源和冷源，在体系中形成温度梯度来计算热导率，后面简称NEMD。这两种方法各有优劣，但是一个核心问题在于，其背后的声子图像一直以来是不甚清晰的。

 

本文主要讨论NEMD。众所周知，如何建立非平衡对于计算结果影响极大。其中的主要问题包括：

1．体系长度、热源长度、热源类型很大程度上会影响计算结果，导致了对于同样的原子模型，同样的相互作用势，算出来的热导率不一样。

2. 提取热导率的方式也会影响计算结果。有时在热源附近会有明显的非线性现象，因此一般用线性区间的温差和温度梯度来提取热导率。因为计算有尺寸效应，因此一般使用线性外推的方式来得到体块材料的热导率。

3. NEMD两侧边界的形貌会明显影响界面热阻的计算结果(Liang and Keblinski Phys. Rev. B 90, 075411,2014)。

 

上述问题直接导致了对NEMD计算结果缺乏真正的理解。那么，NEMD是否有统一的基于声子图像的物理解释呢？

 

最近，上海交通大学密西根学院鲍华课题组与美国橡树岭国家实验室冯天力博士，芬兰阿尔托大学樊哲勇博士，上海交通大学机动学院顾骁坤副教授以及美国普渡大学的Mark Lundstrom教授等人合作，在Phys. Rev. B上发表了一篇题为“Unification of Non-equilibrium