[論文レビュー] 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₂のような系において、弱い結合と高い界面品質が特徴である場合の熱伝達係数の定量的推定を提供すること。
• 界面スプリング定数とフォノンダイナミクスがナノスケール界面における熱伝導度に与える影響を明確にすること。
• 理論的予測とナノ材料、特に2次元膜/基板系における実験的測定結果を一致させること。

提案手法

• 界面応力を表面変位差の線形関数として定式化する：$\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_0$と$u_1$の連立式を$u_{\text{0f}}$で表現する。
• 揺動散逸定理を適用し、熱揺動のパワースペクトル密度を動的行列の虚部と関連付ける：$\langle |u_{\text{0f}}(\mathbf{q},\omega)|^2 \rangle \propto \text{Im} M_1(\mathbf{q},\omega)$。
• 時間平均された応力と速度の積を用いてエネルギーフラックスを導出し、フーリエ空間に変換し、波数および周波数で統合する。
• 最終的に、界面結合$K$、動的行列$M_0$および$M_1$の虚部、および系応答の全周波数および波数依存性を含む熱伝達係数$\alpha$の式を導出する。

実験結果

リサーチクエスチョン

• RQ1 ファンデルワールス力のような弱い界面結合は、固体界面を貫通する熱伝導度にどのように影響するか？
• RQ2 音響的・拡散的不一致モデルを超える一般化された解析的モデルを、平坦で弱い結合性を示す界面を介した熱伝導に適用できるか？
• RQ3 2次元膜/基板系の熱伝導度を決定づける界面スプリング定数とフォノンダイナミクスの役割は何か？
• 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₂上にグラフェンを配置した系の熱伝達係数を、実験データと良好に一致する定量的予測を提供し、理論的手法の妥当性を裏付けた。
• 理論は、粗さが複数の長尺度にまたがる場合に、熱伝導が実接触面積から分離可能であることを示しており、これは直感に反するが、特定の条件下で成り立つ。
• 本モデルは、拡散抵抗と界面温度ジャンプ抵抗の両方を考慮しており、全熱抵抗がこれら2つの寄与の直列和として表されることを示している：$1/\alpha \approx 1/\alpha_{\text{spred}} + 1/\alpha_{\text{c}}$。
• この枠組みは、2次元材料/基板のような任意の平坦で弱い結合性界面に一般化可能であり、波数カットオフを調整することで低温領域への応用も可能である。