[论文解读] Heat Transfer between Weakly Coupled Systems

B. N. J. Persson, A. I. Volokitin+1 more|arXiv (Cornell University)|Aug 26, 2010

Adhesion, Friction, and Surface Interactions参考文献 1被引用 4

一句话总结

本文提出了一种针对弱耦合平面界面（如石墨烯在非晶SiO₂上）的通用解析框架，结合线性响应理论与涨落-耗散定理，推导出考虑界面弹簧常数与声子动力学的热传导系数的闭式表达式，其结果与石墨烯/SiO₂体系的实验数据高度一致。

ABSTRACT

We study the heat transfer between weakly coupled systems with flat interface. We present simple analytical results which can be used to estimate the heat transfer coefficient. As applications we consider the heat transfer across solid-solid contacts, and between a membrane (graphene) and a solid substrate (amorphous SiO2). For the latter system the calculated value of the heat transfer coefficient is in good agreement with the value deduced from experimental data.

研究动机与目标

• 开发一种适用于平面弱耦合固态界面的通用解析理论，尤其关注界面相互作用主要由范德华力主导的情形。
• 解决传统声学与扩散失配模型的局限性，后者假设强界面散射且未考虑弱耦合效应。
• 为如石墨烯在非晶SiO₂上的体系提供热传导系数的定量估算，此类体系中弱耦合与高界面质量是关键因素。
• 阐明界面弹簧常数与声子动力学在决定纳米尺度界面热导率中的作用。
• 调和理论预测与纳米材料实验测量结果的一致性，尤其针对二维膜在基底上的体系。

提出的方法

• 将界面应力表述为表面位移差的线性函数：$\sigma(\mathbf{q},\omega) = K[u_0(\mathbf{q},\omega) - u_1(\mathbf{q},\omega)]$，其中$K$为界面弹簧常数。
• 利用线弹性理论，通过动力学矩阵的逆将表面位移与应力关联：$u_1(\mathbf{q},\omega) = M_1(\mathbf{q},\omega)\sigma(\mathbf{q},\omega)$。
• 引入一个描述非相互作用固体中热噪声的随机涨落位移$u_{\text{0f}}$，并推导出以$u_{\text{0f}}$表示的$u_0$与$u_1$的耦合方程。
• 应用涨落-耗散定理，将热涨落的功率谱密度与动力学矩阵的虚部关联：$\langle |u_{\text{0f}}(\mathbf{q},\omega)|^2 \rangle \propto \text{Im} M_1(\mathbf{q},\omega)$。
• 通过时间平均应力与速度的乘积推导净能量通量，转换至傅里叶空间，并对波矢与频率进行积分。
• 最终得到热传导系数$\alpha$的表达式，其依赖于界面耦合$K$、动力学矩阵$M_0$与$M_1$的虚部，以及系统响应的完整频率与波矢依赖性。

实验结果

研究问题

• RQ1 弱界面耦合（如范德华力）如何影响固态界面间的热导率？
• RQ2 能否为平面弱耦合界面的热传导建立一个超越声学与扩散失配模型的通用解析模型？
• RQ3 界面弹簧常数与声子动力学在决定二维膜/基底体系热导率中的作用是什么？
• RQ4 该理论模型对石墨烯在非晶SiO₂上体系的实验观测热传导系数预测效果如何？
• RQ5 在具有多尺度粗糙度的体系中，热传导在何种条件下会与实际接触面积解耦？

主要发现

• 推导出的热传导系数$\alpha$由一个包含界面耦合$K$、动力学矩阵$M_0$与$M_1$的虚部以及完整响应函数的频率与波矢积分表达式给出：$\alpha = \frac{4A^{*}}{(2\pi)^3} \int_0^\infty d\omega \frac{\partial \Pi(\omega)}{\partial T} \frac{\text{Im}(KM_0)\text{Im}(KM_1)}{|1 + K(M_0 + M_1)|^2}$。
• 在高温与弱耦合条件下，表达式简化为$\alpha \approx \frac{4k_B A^{*}}{(2\pi)^3} \int_0^\infty d\omega \, \text{Im}(KM_0)\text{Im}(KM_1)$，该结果在爱因斯坦振子极限下与已知理论一致。
• 该模型对石墨烯在非晶SiO₂上的热传导系数预测与实验数据具有良好定量一致性，验证了理论方法的有效性。
• 理论表明，当粗糙度跨越多个长度尺度时，热传导可与实际接触面积解耦，这一反直觉结果在特定条件下成立。
• 该模型同时考虑了扩展电阻与界面温度跳跃电阻，表明总热阻为两者串联之和：$1/\alpha \approx 1/\alpha_{\text{spred}} + 1/\alpha_{\text{c}}$。
• 该框架具有普适性，适用于任何平面弱耦合界面，包括二维材料与基底的界面，且可通过调整波矢截断扩展至低温区域。