非調和力定数・格子熱伝導率の算出：ALAMODEとの連携#

Advance/PHASE Logo

半導体デバイスや熱電変換材料の設計において、物質の「熱の伝わりやすさ」を微視的な原子振動（フォノン）の観点から予測することは重要です。本解析では、第一原理計算ソフトウェアAdvance/PHASEを力の評価エンジンとして活用し、フォノン解析パッケージALAMODEと連携させることで、シリコン（Si）のフォノン分散および格子熱伝導率を高精度に算出した事例を紹介します。

Keywords: 第一原理計算 (DFT), ALAMODE, フォノン分散, 非調和力定数, 格子熱伝導率, ボルツマン輸送方程式 (BTE), 有限サイズ効果

理論背景と計算手順#

格子熱伝導率の定式化#

格子の振動による熱伝導率は、フォノンの緩和時間近似（RTA）に基づき、ボルツマン輸送方程式（BTE）を解くことで得られます。フォノン解析パッケージALAMODEは、第一原理計算から抽出した調和・非調和力定数を用いてこのBTEを解くことにより、物質固有の熱伝導率を微視的な観点から定量的に予測します [1, 2]。熱伝導率 $\kappa$ は以下の式で表されます。

$\kappa_{\alpha\beta} = \frac{1}{\Omega N_q} \sum_{q, \nu} C_{q, \nu} v_{q, \nu, \alpha} v_{q, \nu, \beta} \tau_{q, \nu}$

ここで、式中の各変数は以下の通りです。

$\alpha, \beta$ : 空間の方向（例: $x, y, z$ ）
$\Omega$ : ユニットセルの体積
$N_q$ : サンプリングした $q$ 点（波数ベクトル）の総数
$q$ : フォノンの波数ベクトル
$\nu$ : フォノンの分枝（バンド）インデックス
$C_{q, \nu}$ : 波数 $q$ 、バンド $\nu$ のフォノンのモード比熱
$v_{q, \nu, \alpha}$ : フォノンの群速度の $\alpha$ 成分
$\tau_{q, \nu}$ : フォノンの緩和時間（寿命）

寿命 $\tau$ を決定するには、原子間の3次以上の非調和力定数（IFC）を考慮したフォノン-フォノン散乱強度の評価が不可欠です。

Advance/PHASEを用いた解析フロー#

本事例では、以下のステップで解析を実施しました。

ステップ	作業内容	使用ツール
1. 変位構造生成	平衡構造から原子を僅かに動かした複数のスーパーセル構造を生成。	ALAMODE (alm suggest)
2. 第一原理計算	生成された各構造に対し、各原子に働く「力」を算出。	Advance/PHASE
3. 力定数の抽出	変位量と力の関係から、フィッティングにより2次および3次の力定数を決定。	ALAMODE (alm optimize)
4. 物性値の算出	力定数を用いてフォノン分散、熱伝導率、寿命などを計算。	ALAMODE (anphon)

熱伝導率の計算精度の確保
高精度な熱伝導率を得るためには、Advance/PHASEのDFT計算において十分な計算条件（k点メッシュなど）を用いることが必須です。ソフトウェア間の連携にあたっては、専用のスクリプトを利用して入出力の整合性（特にAngstromやBohrといった単位系、原子の座標・並び順）を厳密に保つ必要があります。また、BTEを解く際には、長波長フォノンの寄与を正しく評価するために十分な密度のq点メッシュが必要となります。加えて、非調和力定数を空間的にどこまで考慮するか（カットオフ半径）も重要です。

計算結果と考察#

1. フォノン分散関係（Phonon Dispersion）#

Advance/PHASEで算出した2次の力定数（調和項）に基づき、Siのブリルアンゾーンにおけるフォノン分散曲線を評価しました（図1）。

Si Phonon Band Structure

図1. シリコンのフォノン分散関係。

音響フォノンの立ち上がりや $\Gamma$ 点での縮退挙動など、シリコン特有のフォノン特性がよく再現されています。これは、DFT計算が原子に働く力を高い精度で評価できていることを示しています。

2. 格子熱伝導率の温度依存性#

3次の非調和力定数を用い、10 K から 1000 K までの格子熱伝導率を算出しました。図2に、今回の計算結果（Advance/PHASE + ALAMODE）および実験値の比較を示します。

Thermal Conductivity Comparison

図2. 格子熱伝導率の温度依存性。緑線(calculation)が本ワークフローによる結果、青四角（experiment）が実験値 [3] です。

計算の結果、本ワークフローによる格子熱伝導率は、温度上昇に伴ってフォノン-フォノン散乱（Umklapp散乱）が活発化し、典型的な $1/T$ の減衰挙動を示すことが確認できました。注目すべき点として、室温付近から高温域に至るまで、本計算結果は実験値 [3] と良好に一致しています。800 K以上の高温域で見られるわずかな過大評価は、本計算において4次以上の高次のフォノン散乱プロセスが省略されていることに起因すると考えられます。全体として、非調和性を決める微小な力の変動をDFT計算で高精度に評価したことで、実用上十分な精度を持つ熱伝導シミュレーションが実現できたことを裏付けています。

3. 試料の有限サイズ効果（境界散乱）の考慮#

バルクの熱伝導率計算では、低温域でフォノン散乱が消失するため値が発散しますが、実際の測定試料では端（境界）での散乱が支配的になります。マシーセンの法則に基づき、一般的なバルク単結晶の測定サイズに相当する 1 mm とした場合の有限サイズ効果を評価しました（図3）。

Finite Size Effect on Thermal Conductivity

図3. 有限な試料サイズによる熱伝導率の変化。

計算の結果、有限なサイズ（1 mm）では、100 K以下の極低温域において境界散乱により熱伝導率が著しく抑制され、実験値の傾向に近づくことが確認できます。もし、この試料サイズをナノワイヤや薄膜といった数十〜数百 nm スケールまで微細化した場合、室温付近でも顕著なサイズ効果が現れます。このような解析により、ナノデバイスにおける放熱特性の設計指針を、第一原理計算から定量的に見積もることが可能となります。

まとめ#

本事例では、第一原理計算ソフトウェアAdvance/PHASEとフォノン解析パッケージALAMODEをシームレスに連携させることにより、原子スケールの第一原理計算からマクロな熱伝導物性までを一気通貫で解析するパイプラインを実証しました。特に、シリコンの熱伝導率において実験値とよく一致する結果が得られました。この解析手法は、次世代パワーデバイスの熱マネジメントや熱電材料の開発において、強力な設計ツールとなります。

本解析の詳細や、研究への適用可能性に関するご相談はこちら

お問い合わせ

参考文献#

T. Tadano, Y. Gohda, and S. Tsuneyuki, "Anharmonic force constants extracted from first-principles molecular dynamics: applications to heat transfer simulations", J. Phys.: Condens. Matter 26, 225402 (2014).
A. Togo, L. Chaput, and I. Tanaka, "Distributions of phonon lifetimes in Brillouin zones", Phys. Rev. B 91, 094306 (2015).
C. J. Glassbrenner and Glen A. Slack, "Thermal Conductivity of Silicon and Germanium from 3°K to the Melting Point," Phys. Rev. 134, A1058 (1964).