[論文レビュー] FourPhonon: An extension module to ShengBTE for computing four-phonon scattering rates and thermal conductivity

Zherui Han, Xiaolong Yang+3 more|arXiv (Cornell University)|Apr 10, 2021

Thermal properties of materials参考文献 43被引用数 433

ひとこと要約

FourPhonon は、適応的幅広化スキームを用いて第一原理的に四フォノン散乱率および格子熱伝導率を計算可能な、ShengBTE フレームワークへのオープンソース拡張である。既存のワークフローと統合可能で、空間群と並進対称性を考慮したスクリプトにより四階相互作用力定数を効率的に計算し、正確な熱伝導率予測を実現する。Si、BAs、LiCoO2 における検証では、四フォノン散乱が熱輸送に顕著に寄与しており、特に室温における BAs や LiCoO2 ではその影響が顕著であることが示された。

ABSTRACT

FourPhonon is a computational package that can calculate four-phonon scattering rates in crystals. It is built within ShengBTE framework, which is a well-recognized lattice thermal conductivity solver based on Boltzmann transport equation. An adaptive energy broadening scheme is implemented for the calculation of four-phonon scattering rates. In analogy with $thirdorder.py$ in ShengBTE, we also provide a separate python script, $Fourthorder.py$, to calculate fourth-order interatomic force-constants. The extension module preserves all the nice features of the well-recognized lattice thermal conductivity solver ShengBTE, including good parallelism and straightforward workflow. In this paper, we discuss the general theory, program design, and example calculations on Si, BAs and $\mathrm{LiCoO_2}$.

研究の動機と目的

• 結晶における四フォノン散乱率および熱伝導率を計算する初のオープンソースツールを開発すること。
• 広く使用されている ShengBTE フレームワークに四フォノン機能を追加し、並列処理およびワークフローの単純さを維持すること。
• 四フォノン散乱率計算の収束性と精度を向上させるために、適応的エネルギー幅広化スキームを実装すること。
• 点群および並進対称性を活用した、四階相互作用力定数（IFC）を効率的に計算する Python スクリプト（Fourthorder.py）を提供すること。
• 代表的な材料（Si、BAs、LiCoO2）を対象に手法を検証し、熱輸送における四フォノン効果の重要性を示すこと。

提案手法

• Boltzmann輸送方程式（BTE）フレームワーク内での四フォノン散乱率計算を可能にする新モジュール FourPhonon を ShengBTE に拡張する。
• 集団速度の差に基づく適応的エネルギー幅広化スキームを実装し、散乱率の収束性と精度を向上させる。
• 第一原理計算と対称性低減を活用して四階相互作用力定数（IFC）を計算する Python スクリプトである Fourthorder.py を導入する。
• 四フォノン散乱を単一モード緩和時間近似（SMRTA）で計算し、反復的三フォノン BTE 解と統合するように変更されたワークフローを採用する。
• CONTROL ファイルで four_phonon=.TRUE. を指定することで、四フォノン計算を有効化可能にし、既存の ShengBTE 入出力構造を再利用できる。
• 再結合、再分配、分裂の各チャネルへの散乱プロセスの分解を可能とし、Umklapp および通常過程の寄与を別々に出力する。

実験結果

リサーチクエスチョン

• RQ1 Si や BAs、LiCoO2 などの材料において、四フォノン散乱は格子熱伝導率にどの程度寄与しているか？
• RQ2 適応的幅広化スキームは、四フォノン散乱率計算の精度と収束性を効果的に向上させることができるか？
• RQ3 BAs や LiCoO2 のような強い非調和性を示す材料では、四フォノン過程が三フォノン過程をどれほど上回るか？
• RQ4 三フォノンのみを考慮する従来のアプローチと比較して、四フォノン散乱を組み込むことで予測される熱伝導率にどのような影響があるか？
• RQ5 四フォノン計算の計算コストとスケーラビリティはどの程度で、オープンソースツールによって実用的に有効化可能か？

主な発見

• シリコンでは、室温における四フォノン散乱率は無視できるほど小さく、許容可能なプロセス数は約 10^3 個にとどまり、この系ではその寄与が小さいことが確認された。
• ホウ素ヒ素化物（BAs）では、四フォノン散乱が極めて強く、利用可能な散乱プロセス数は約 10^10 個に達し、特定の周波数範囲では三フォノン散乱を上回ることが確認された。これは先行予測および実験結果とも整合的である。
• LiCoO2 では、全周波数スペクトルにわたり四フォノン散乱率が三フォノン率と同等の大きさを示し、特に低周波数の音響モードおよび光学モードで顕著である。
• 再分配チャネル（λ + λ′ → λ′′ + λ′′′）が四フォノン散乱に最も寄与しており、LiCoO2 では Umklapp 过程がこのチャネルを支配している。
• 300 K における LiCoO2 の予測熱伝導率は、κ∥ = 9.35 W/mK および κ⊥ = 1.39 W/mK であり、文献値（9.7 および 1.4 W/mK）と良好に一致している。
• プログラムは詳細な散乱率、状態密度、および異なる散乱チャネルからの寄与を正しく計算・出力できており、四フォノン効果の体系的分析を可能にしている。