[論文レビュー] GiftBTE: An efficient deterministic solver for non-gray phonon Boltzmann transport equation

Yue Hu, Ru Jia+6 more|arXiv (Cornell University)|Jun 25, 2023

Thermal properties of materialsMaterials Science被引用数 3

ひとこと要約

GiftBTE は、結晶性材料におけるサブミクロンスケールの熱輸送を効率的かつ高精度にシミュレートできる、オープンソースで決定論的な非灰色フォノンBoltzmann輸送方程式（BTE）の数値解法である。定常状態では2次精度の陰解法と合成反復スキームを採用し、過渡問題では陽解法を用いることで、優れた計算効率を実現し、ナノ構造、トランジスタ、レーザ加熱プロセスの現実的な3次元シミュレーションを可能にした。

ABSTRACT

Advances in nanotechnology have facilitated the exploration of submicron thermal transport. At this scale, Fourier's law is no longer applicable, and the governing equation for thermal transport is the phonon Boltzmann transport equation (BTE). However, the availability of open-source solvers for the phonon BTE is limited, impeding progress in this field. This study introduces an open-source package, GiftBTE, for numerically solving the non-gray phonon BTE. GiftBTE employs deterministic solutions and provides both steady-state and transient solvers. For the steady-state solver, GiftBTE employs the implicit discrete ordinates method (DOM) with second-order spatial accuracy and the synthetic iterative scheme. For the transient solver, GiftBTE employs the explicit DOM with second-order spatial accuracy. This package demonstrates excellent computational efficiency, enabling realistic three-dimensional simulations of devices and materials. By interfacing with first-principles calculations, this solver enables parameter-free computation of submicron thermal transport. The application of GiftBTE includes, but is not limited to, computing the thermal conductivity of nanostructures, predicting temperature rises in transistors, and simulating laser heating processes.

研究の動機と目的

• 複雑な幾何形状における非灰色フォノンBTEの一般用途で利用可能なオープンソースのソルバーが不足しているという問題に対処すること。
• バルクまたは単純な幾何形状に限定され、過渡解析やパラメータフリーな能力に欠ける既存のソルバーの限界を克服すること。
• 第一原理フォノン特性を用いて、現実的な3次元材料およびデバイスにおけるサブミクロンスケールの熱輸送を、高精度で効率的かつスケーラブルにシミュレートすること。
• ナノエレクトロニクスおよび熱電材料における熱輸送解析のためのCAE風のツールを提供し、熱伝導率の予測やトランジスタにおける温度上昇の推定といった応用を支援すること。
• 原子スケールの手法（例：MD, NEMD）とマクロスケールの熱拡散方程式の間のギャップを埋め、サブミクロンスケールで決定論的かつ非フォーリエ的シミュレーションを可能にすること。

提案手法

• 非灰色フォノンBTEを解く決定論的アプローチとして、離散座標法（DOM）を採用し、モンテカルロ法に内在する統計的ノイズを回避する。
• 定常状態および過渡問題の両方のソルバーにおいて、空間的に2次精度を実現し、メッシュ要件を低減するとともに解の精度を向上させる。
• 近拡散領域における収束速度を向上させるために、定常状態ソルバーで合成反復スキーム（SIS）を採用し、標準的な反復スキームが遅く収束する問題を解決する。
• 過渡シミュレーションに、2次精度を有する陽解法のDOM定式化を用い、時間分解能のある熱的ダイナミクスの分析を可能にする。
• phono3py や類似のツールを通じて第一原理フォノン計算とインターフェースを設け、経験的フィッティングなしに温度依存および波数依存のフォノン特性を取得する。
• 任意の結晶性材料および複雑な3次元幾何形状（ナノポーラス媒体、薄膜、デバイススケール構造など）をサポートする。

実験結果

リサーチクエスチョン

• RQ1 決定論的で2次精度のBTEソルバーは、現実的な3次元サブミクロンスケールの熱輸送問題に対して、高い計算効率と精度を達成できるか？
• RQ2 合成反復スキームは、標準的な反復スキームと比較して、近拡散領域における収束速度をどのように向上させるか？
• RQ3 非フォーリエ効果は、時間領域熱反射率測定（TDTR）実験のような過渡的熱応答において、どの程度顕著に現れるか？
• RQ4 第一原理入力に基づいて、ナノスケールのトランジスタやレーザ加熱材料における温度上昇を高精度に予測できるか？
• RQ5 GiftBTEフレームワークは、単純な幾何形状を越えた複雑なナノ構造において、パラメータフリーな熱輸送シミュレーションを可能にするか？

主な発見

• 2次空間精度と合成反復スキームを有する定常状態ソルバーは、1次精度手法と比較してメッシュ数を顕著に削減し、計算効率を向上させた。
• 過渡ソルバーは非フォーリエ的熱的ダイナミクスを適切に捉えており、TDTRシミュレーションにおいてフォーリエ法と比較して温度の低下が遅い傾向を示し、顕著なボールティック輸送効果を示した。
• TDTRの事例研究において、表面温度の低下がフォノンBTEの下では熱拡散方程式よりも遅く、64コアクラスタ上で16時間のシミュレーションを要した。
• 薄膜や多孔質材料のようなナノ構造における熱伝導率の計算を、第一原理フォノンデータからの完全なパラメータ化により高精度に実行可能である。
• 変調されたパルスを用いたレーザ加熱プロセスを、実験的挙動を再現する温度応答と現実的なビームプロファイルおよび浸透深さを伴ってシミュレート可能である。
• トランジスタや熱電的ヘテロ構造を含む複雑なデバイスの3次元シミュレーションにおいて、スケーラビリティとロバスト性を示し、フォノンモード分解された熱フラックスおよび温度場への完全なアクセスを提供した。