Revisiting Theoretical Modeling of Seebeck Coefficient of Semiconductors

技術要約：半導体のゼーベック係数の理論的モデリングの再検討

問題提起
本論文は、ボルツマン輸送方程式（BTE）や比熱流の定義に依存する既存の手法において、半導体のゼーベック係数（$\alpha_S$）をモデル化する際に生じる根本的な熱力学的欠陥を特定している。著者は、従来のアプローチにおける3つの主要な問題を指摘している：

1. 潜在的な混同： 定義において、電気化学ポテンシャル（フェルミ準位、$E_F$）と静電ポテンシャル（$\varphi$）が頻繁に混同されており、それが誤った導出を招いている。
2. 恣意的な無視： BTEによる導出において、分布関数の勾配を計算する際、伝導帯底（$E_C$）およびフェルミ準位（$E_F$）の空間勾配が、開回路（OC）条件下であってもサンプル内で変化しているにもかかわらず、それらが恣意的に無視されている。
3. 不適切な熱流： 局所熱流（$J_Q = T J_S = J_E - \mu J_N$）の定義は、エントロピー生成がゼロである局所的な準静的（可逆）プロセスの強い仮定に依存している。これは基本的な非平衡熱力学の原理に矛盾するものであるが、BTE、ランダウアー・ビュリカー、およびクボ公式の適用において広く用いられている。

手法
これらの欠陥に対処するため、著者は開回路条件下におけるソレ効果を組み込んだドリフト拡散方程式に基づく代替手法を提案する。このアプローチは、明示的な熱流の提示や、電気化学ポテンシャルと静電ポテンシャルの間の曖昧さを回避するものである。

導出は以下の3段階で進められる：

1. 近平衡状態（バンド項）： 局所平衡（LEA）および局所的な電荷中性を仮定し、$\alpha_S$ のバンド項に関する閉じた形式のモデルを導出する。正味の電荷束をゼロに設定し、一般化されたアインシュタイン関係式を利用することで、「可逆モデル」が回収される。このモデルは、電子速度や緩和時間の項を明示的に排除し、代わりにフェルミ準位の空間勾配に依存する。
2. フォノン・ドラッグ項： 非平衡フォノンとの運動量交換から生じる電荷束（$J_d$）を加えることで、手法をフォノン・ドラッグ効果へと拡張する。その結果得られるフォノン・ドラッグ項の式は、従来のBTEベースのモデルと同一の形式であることが示される。
3. 弾道的ケース： 弾道的領域（LEAが成立しない領域）では、ランダウアー定式化を用いて高温端と低温端からの電荷束を計算し、従来の弾道的ゼーベック係数の式を導出する。

本論文は、これら（BTE、可逆、および弾道的モデル）を、エネルギー指数 $r$ に依存する単一の積分形式へと統合する。ここで $r$ は、散乱メカニズムや材料の次元性に基づいて変化する。

主な貢献

• 熱力学的修正： 標準的なBTEの導出および熱流の定義における熱力学的な不整合を厳密に批判している。
• 可逆モデルの導出： 著者の先行研究で提示された「可逆モデル」が、速度や緩和時間の項を落としているにもかかわらず、近平衡状態におけるゼーベック係数のバンド項として妥当であることを示している。
• 理論的境界： 本研究は、可逆モデルがゼーベック係数のバンド項に対する理論的な下限（3D材料の場合は $r=1/2$ に相当）を提供することを確立している。一方で、弱スクリーニングされたイオン化不純物散乱を用いたBTEモデルは、上限（$r=3$）を提供する。
• 統一された枠組み： 状態密度、速度、および緩和時間のエネルギー依存性に基づき、可逆、弾道的、およびBTEモデルを比較するための統一された数学的枠組みを提示している。

結果
高濃度n型シリコン（$n_d > 1 \times 10^{19} \text{ cm}^{-3}$）を用いたケーススタディを行い、GeballeらおよびSadhuらの実験データに対して3つのモデルを比較した：

• 可逆モデル： 低い計算コストで「まずまず良好な」予測を提供した。このモデルは一貫して実験データの低限として機能し、ドーピング濃度が $2.7 \times 10^{19}$、$6.0 \times 10^{19}$、$7.0 \times 10^{19} \text{ cm}^{-3}$ の場合、それぞれ 13%、18%、5% の平均相対偏差を示した。
• 弾道的モデル： 混合した性能を示し、偏差は 3%、6%、16% であった。ナノワイヤのケース（$7.0 \times 10^{19} \text{ cm}^{-3}$）ではデータを過大評価した。
• BTEモデル： 弱スクリーニングされたイオン化不純物散乱（上限の仮定）を用いた場合、このモデルは実験値を大幅に過大評価し、平均相対偏差は 60% から 90% に達した。これは、BTEの結果が緩和時間定式化の選択に対して極めて敏感であることを浮き彫りにしている。

意義および主張
本論文は、提案された可逆モデルが、特に近平衡および可逆的なケースにおいて、熱電特性を評価するための簡潔で、理論的に健全かつ計算効率の高い代替手段を提供すると主張している。複雑な散乱時間の計算を排除し、熱力学的パラドックスを回避することで、このモデルはハイスループット計算や材料スクリーニングのための実行可能なツールとして提示されている。著者は、このモデルは不可逆的なメカニズム（散乱）を無視しているものの、ゼーベック係数の下限を正確に捉えており、熱電特性評価のための堅牢なベースラインとなることを強調している。