[論文レビュー] Direct Calculation of Modal Contributions to Thermal Conductivity via Green-Kubo Modal Analysis: Crystalline and Amorphous Silicon

Wei Lv, Asegun Henry|arXiv (Cornell University)|Mar 27, 2015

Thermal properties of materials参考文献 53被引用数 155

ひとこと要約

本稿では、フォノン速度やフォノンガスモデルに依存せずに、熱伝導率へのモード寄与を直接計算する第一原理的手法であるGreen-Kuboモーダル解析（GKMA）を紹介する。格子力学とGreen-Kubo形式を組み合わせることで、完全な非調和性およびモード間相関を捉え、結晶性および非晶性ケイ素における温度依存熱伝導率の高精度な予測が可能となり、量子補正を施した場合、非晶性ケイ素に関してはこれまでで最も実験データと一致する結果を得た。

ABSTRACT

In this letter we derive a new method for direct calculation of the modal contributions to thermal conductivity, which is termed Green-Kubo modal analysis (GKMA). The GKMA method combines the lattice dynamics formalism with the Green-Kubo formula for thermal conductivity, such that the thermal conductivity becomes a direct summation of modal contributions, where one need not define the phonon velocity. As a result the GKMA method can be applied to any material/group of atoms where the atoms vibrate around stable equilibrium positions, which includes not only crystalline line compounds, but also random alloys, amorphous materials and even molecules. The GKMA method provides new insight into the nature of phonon transport, as it casts the problem in terms of mode-mode correlation instead of scattering, and provides a general unified formalism that can be used to understand phonon-phonon interactions in essentially any class of materials or structures where the atoms vibrate around stable equilibrium sites.

研究の動機と目的

• フォノン速度やフォノンガスモデルに依存しない、モード寄与を計算する一般化された形式を構築すること。
• 周期的でない材料（非晶性ケイ素など）における熱輸送を記述する既存モデルの限界を克服すること。フォノンモードが明確な平面波として定義されない状況を想定する。
• 熱伝導率計算に完全な非調和性およびモード間相関を組み込み、従来の手法が失敗する系に対しても有効な手法を提供すること。
• 古典的分子動力学とモード固有の比熱に対する量子補正を組み合わせることで、温度依存熱伝導率の高精度な予測を可能にすること。

提案手法

• 格子力学を用いて、原子の変位および速度を変換行列を介した正規モード成分の重ね合わせとして表現する。
• モード速度成分に基づき、ハードイの熱フラックス演算子を用いて、各モードの時間依存的熱フラックス寄与を導出する。
• モード分解熱フラックスにGreen-Kubo公式を適用し、各モードが熱伝導率に寄与する部分を直接計算可能にする。
• 時間相関関数を用いてモード間クロス相関を計算し、異なる周波数のモード間の相互作用を明らかにする。
• 各モードの寄与を、量子比熱と古典比熱の比を用いてスケーリングすることで、古典的MD結果に対する量子補正を可能にする。
• 結晶性および非晶性ケイ素への適用を通じて、実験データおよび先行理論モデルと比較することで、手法の妥当性を検証する。

実験結果

リサーチクエスチョン

• RQ1 フォノン速度や群速度を仮定せずに、熱伝導率へのモード寄与を直接計算できるか？
• RQ2 非晶性材料（ケイ素など）におけるモード間相関、特に非対角項が熱伝導率にどのように寄与するか？
• RQ3 分子動力学で捉えた完全な非調和性は、調和近似と比較して熱伝導率の予測精度をどの程度向上させるか？
• RQ4 GKMAフレームワークにおいて古典的MD結果に量子補正を施すことで、非晶性ケイ素における温度依存熱伝導率を高精度に予測できるか？
• RQ5 GKMA形式は、局在モード（locon）と拡散モード（propagon/diffuson）の熱輸送における役割に関する新たな物理的洞察を明らかにできるか？

主な発見

• GKMA手法はフォノン速度を必要とせず、非周期的系（非晶性材料など）への応用が可能である。
• 非晶性ケイ素において、量子補正を施したGKMA手法は、これまでで最も実験データと一致する結果を得ており、先行モデルを上回っている。
• 室温では、モード間クロス相関が全熱伝導率の約30%を占め、非調和性の顕著な寄与が示された。
• クロス相関行列は、約16 THzでモード相互作用に顕著な遷移を示しており、事前にモードの性質を知らなくても、局在モード（locon）の発生と一致している。
• loconモードの自己相関は無視できるほど小さく、熱伝導率への寄与が極めて小さいことと整合的である。
• 低温域におけるモード寄与のシフトを正しく捉えており、GKMAで得た周波数分解寄与に量子補正を施すことで、温度依存性の高精度な再現が可能であることが示された。