電気熱シミュレーションチュートリアル（Part B）：熱メッシュ、境界条件、および実験フィット
1. 導入：Part B の範囲と焦点
Part A では、電気熱シミュレーションを実行し、 主な発熱メカニズム （キャリア輸送による発熱、再結合による発熱、および寄生損失） を特定し、印加バイアスの関数として格子温度の出力を抽出する方法を学びました。 最初のパートでは、どのような熱が発生し、それがデバイス内のどこに現れるかという what に重点を置きました。この Part B では、問題のもう一方の側面、 すなわち 熱がどのようにデバイスから輸送されて離れるか に焦点を当てます。 これは熱メッシュによって支配され、特に 熱境界条件 によって決まります。 境界条件は、実際にデバイスがどのように実装され、冷却され、またはヒートシンクに接続されているかを表現し、 その結果生じる温度上昇に一次の影響を与えることがよくあります。 電気的挙動と発熱項が同一の二つのシミュレーションでも、 熱境界の仮定が異なるだけで、温度分布は大きく変わり得ます。
したがって、電気熱シミュレーションでは通常、 電気領域よりも大きい熱領域 を用いる理由、 明示的にメッシュを切らずに巨視的ヒートシンクを表現するために有効境界条件がどのように使われるか、 そして温度依存の材料量がなぜ離散温度グリッド上で事前計算されるのかが分かるようになります。 このチュートリアルでは、有機ダイオードの例を引き続き使用します。これは実験的な裏付けがあり、 測定可能な自己発熱効果を含んでいるためです。しかし、ワークフロー自体は完全に一般的です。同じ概念は、無機半導体、薄膜デバイス、 パワーエレクトロニクス、光検出器、そして電力散逸が温度を介して輸送にフィードバックするあらゆる構造に適用されます。
Part B の終了時には、以下を理解しているはずです：
- なぜ熱メッシュが電気メッシュとは別であるのか、
- 熱境界条件がどのように熱抽出を制御するのか、
- 効率のために内部で温度離散化がどのように使われるのか、そして
- 実験データと一致させるために、なぜ自己発熱を組み込む必要がある場合があるのか。
2. 熱メッシュと境界条件
メインのシミュレーションウィンドウから開始し、Thermal リボンへ移動します。これは Part A で導入したものと同じリボンです： ここには発熱源のトグルと、主要な設定コントロール（熱パラメータ、熱メッシュ、境界条件）が含まれています。 このリボンは ?? に示されています。
2.1 熱メッシュと、それが電気メッシュではない理由
Thermal mesh をクリックして熱メッシュエディタを開きます。一見すると、OghmaNano の他の場所で見たメッシュエディタと似ていますが、 これは異なる物理問題、すなわち 熱拡散 を解いています。重要な点は、電気熱シミュレーションでは 2 つのメッシュ を使うということです。なぜなら、電気問題と熱問題は通常、 根本的に異なる長さスケール上に存在するためです：
- 電気メッシュ は輸送（drift–diffusion + Poisson + 再結合/トラップ）を解きます。通常これは、電場やキャリア密度が急激に変化する 活性層およびその近傍界面など、電気的な変化が最も強い場所に集中しています。 それらの変化は ナノメートルからマイクロメートル スケールで起こり得ます。
- 熱メッシュ は格子の熱方程式を解きます。熱は活性層で止まりません。積層全体、 コンタクト、基板、パッケージング、そしてデバイスから熱を除去するあらゆるものへと広がります。これらの経路は一般に ミリメートルからセンチメートル スケールです。
この不一致こそが、単一の「共有メッシュ」が通常は不適切な抽象化である理由です。電気問題は活性領域における高い分解能を必要とし、 熱問題は熱流経路を表現できる十分に大きな領域を必要とします。 例えば：
- 薄膜ダイオードの活性層厚さは 100–300 nm かもしれませんが、基板は 0.5–1 mm の厚さ（ガラス）を持ち、取り付けブロックやステージは センチメートル スケールになり得ます。
- 金属コンタクトは電気的には数十ナノメートルの厚さしかない場合でも、 有機層や酸化物層に比べて熱伝導率が大きいため、熱的には支配的な横方向熱拡散体となることがあります。
この例では、熱領域の Y 高さ が電気的に活性な領域だけに制限されず、 デバイス積層全体の厚さにわたっていることが分かります。これは典型的です：熱輸送には、電流が流れる層だけでなく、熱を伝導するすべての層が関与します。 次に、境界条件を用いて、この熱領域から どのように熱が出ていくことが許されるか を設定します。メッシュは 熱が拡散できる領域を定義し、境界条件はその領域の端で何が起こるかを定義します。
2.2 境界条件：断熱面と有効ヒートシンク
?? に示す境界条件エディタを開いてください。 この例では、大部分の面が Neumann（しばしば「Neumann (==0)」と表示されます）に設定されています。 物理的には、Neumann 境界条件は境界における 法線方向熱流束 を指定します。 それがゼロに設定されると、次を課します：
\[ -k \nabla T \cdot \hat{n} = 0 \]
すなわち、その境界を通って熱は流れない ということです。これらの面は熱的に絶縁されているものとして扱われます。 これはデバイス全体が熱的に孤立していることを意味するわけではなく、単にそれらの 境界が意図された熱除去経路の一部ではないことを熱ソルバーに伝えているだけです。 ここでの例外は ymax 面で、これは Heatsink に設定されています。 この境界は、ymax におけるデバイス温度（ここでは約 300 K）を指定し、 それに加えて有効な ヒートシンク熱伝導率 と ヒートシンク長 （この例ではミリメートル程度） を指定します。
概念的には、これは 有効熱除去モデル です。これにより、熱メッシュに巨視的物体を含めることなく、 シミュレーションされたデバイス積層から、明示的にはメッシュ化されていないシンクへ熱を流出させることができます。 これは重要です。なぜなら、実際のヒートシンクは通常、 薄膜デバイスより何桁も大きいからです：ミリメートルからセンチメートルスケールのシンクをサブミクロン分解能で明示的にメッシュ化することは、 計算効率の面でも物理的観点からも不要です。 熱問題においては、境界条件がデバイスの温度をどれだけ高くするかを大きく決定する ことを強調しておく価値があります。 熱抽出経路が不十分であれば温度は非常に急速に上昇し得ますが、効率的な経路であれば 相当な電力散逸下でも温度を周囲温度近くに保つことができます。
有用なたとえとして、デバイスを 水で満たされる浴槽 と考えることができます。 発熱項は蛇口であり、系にエネルギーを注ぎ込みます。 熱境界条件は排水口です：
- 排水口が全開なら（優れた熱抽出）、水位は低いままです。
- 排水口が部分的に詰まっているなら（限定的な熱抽出）、水位は上昇します。
- 排水口が閉じているなら（熱抽出なし）、浴槽はいずれ溢れます。
このたとえでは、水位が格子温度です。 したがって熱境界条件の目的は、見かけ上のものでも副次的なものでもありません： それらは、最終的にデバイスが低温で動作するのか、 高温になるのか、あるいは破局的に高温になるのかを支配する物理環境を符号化しています。
2.3 温度「メッシュ点」（なぜ温度が離散化されるのか）
熱メッシュエディタには、温度範囲 と 温度点数 も表示されます。これは空間的な熱メッシュではなく、 事前計算に用いられる温度グリッドです。電気／熱モデルの多くの部分では、温度依存量を結合計算中に何度も評価する必要があり、 それらを中間温度ごとに毎回ゼロから再計算するのは計算コストが高くなります。 例えば、選択した統計モデルや状態密度モデルによっては、ソルバーはキャリア密度に関する温度依存関係、 準熱力学関数、および関連するルックアップ量を必要とする場合があります。結合ループが \(T_L\) を更新するたびにそれらを「その場で」評価するのではなく、OghmaNano は 有限個の温度点でバックグラウンドテーブルを事前計算し、実行中にその間を 補間 します。 この例では、テーブルは 290 K から 350 K の間の 7 点 で生成されます。
したがって、温度範囲はモデル選択の一部となります：自己発熱中に予想される温度を十分に包含する必要があります。デバイス温度が 上限温度を超えると、外挿アーティファクト（あるいは設定によってはクランプ）が生じる危険があり、これは通常望ましくありません。実用的な指針としては：予想される動作温度を 余裕を持ってカバーする範囲を選んでください。
3. 熱パラメータ（熱伝導率と緩和時間）
電気熱シミュレーションでは、境界条件が熱がデバイスから どのように 出ていくかを定義しますが、材料パラメータは 熱が積層内を どのように 移動するかを定義します。実際には、温度上昇とホットスポット形成を制御する 最も重要な量の一つが 熱伝導率 \(k\)（\(\kappa\) と書かれることもあります）です。
これらのパラメータを表示または編集するには、Thermal リボンを開き、 Thermal parameters （\(k\) / \(\kappa\) と表示されたボタン）をクリックします。これにより、 ?? に示す層ごとの熱パラメータエディタが開きます。
このエディタは各層に対して 3 つの主要パラメータを提供します：
- 熱伝導率 \(k\)：その層を通して熱がどれだけ効率よく拡散するかを支配します。 \(k\) の低い層は熱ボトルネックとして働き、予測される格子温度上昇を大きく増加させることがあります。
- 電子緩和時間 と 正孔緩和時間：これらは hydrodynamic / energy-transport モデルを使用する場合にのみ必要で、その場合にはキャリア温度が格子温度から乖離し得ます。 標準の格子熱モデルでは、キャリア温度は \(T_L\) に固定されているため不要です。
ほとんどのワークフローでは、熱伝導率は公表値（または利用可能であれば測定値）から選ばれ、 明確な実験的根拠がある場合にのみ調整されます。境界条件とともに、\(k\) は絶対デバイス温度と ホットスポットの空間構造を決める主要因の一つです。
4. 結合電気熱ソルバーの動作原理（およびなぜ遅くなるのか）
OghmaNano における電気熱シミュレーションは、2 つの結合ソルバー から構成されます：
- 完全結合された電気ソルバー（drift–diffusion と Poisson。必要に応じて 再結合とトラップを含む）、および
- 電気解から計算された発熱項を用いて 格子温度場 \(T_L\) のための熱拡散方程式を解く 熱ソルバー。
ある印加電圧において、ソルバーはこれらの段階を一度だけ順番に実行するのではありません。 代わりに、自己無撞着な外側ループを実行します：
- 現在の温度場を用いて電気輸送問題を解く、
- 発熱項（輸送、再結合、寄生、および任意で光学的なもの）を計算する、
- 熱拡散問題を解いて \(T_L\) を更新する、そして
- 電気残差と熱残差の両方が収束するまで繰り返す。
この反復結合こそが、電気熱シミュレーションが純粋な電気シミュレーションより遅い理由です： 各バイアスポイントに対して、電気 Newton 計算と熱計算の両方が、 結合系が安定する前に複数回実行される可能性があります。その結果得られるのは、 温度を後処理として加えた電気計算ではなく、真に 自己無撞着なマルチフィジックス解 です。
この結合は不可欠です。なぜなら、電気問題と熱問題は通常、 非常に異なる物理的長さスケール 上に存在するからです。 電気輸送は活性領域によって支配され、そこでは電場、キャリア密度、 および再結合率が急激に変化します。 熱輸送は、コンタクト、 基板、封止、ヒートシンクを含む熱流経路全体を考慮しなければなりません。 このため、OghmaNano は 熱メッシュ と 熱境界条件 に対して独立した制御を提供します。
また、電気熱フレームワークには自然に 複数の温度 が関与することを認識することも重要です。 一般に、このモデルは 格子温度 \(T_l\)、 すなわち原子格子の温度、電子温度 \(T_e\)、 および 正孔温度 \(T_h\) を区別します。
電子と正孔は移動可能なキャリア集団であり、原理的には、 それぞれの準フェルミ準位を持つのと同様に、それぞれ独自の局所エネルギーを持ち得ます。 ここで用いている標準的な電気熱モデルでは、電子温度と正孔温度は 格子温度に固定されています。したがって、すべてのサブシステムは単一の温度場を共有します。
極端な条件では、OghmaNano はさらに hydrodynamic transport model もサポートしており、その場合 電子温度と正孔温度が格子温度からずれることが許されます。 この高度なモデルは別の場所で説明されており、特殊な状況でのみ必要です。 このチュートリアルを含む大多数のデバイスシミュレーションに対しては、格子ベースの 電気熱モデルで適切かつ十分です。
5. 自己発熱と実験との比較
この例のデバイスは、電気熱モデルを用いて実験データと比較されています。 このチュートリアルではフィッティング自体は行いません。代わりに、既存の比較を用いて、 自己発熱がデバイスの電気応答をどのように変化させるか、 そして高動作電流においてなぜ結合電気熱解析が必要なのかを示します。
比較ワークフローには、Automation ツールバーからアクセスします。これは ?? に示されています。 fit-to-experiment ツールを開くと、このデバイスの実験 JV データが表示されます。 これは ?? に示されています。 この曲線は、バイアス下で測定されたデバイス挙動を表しています。
1 回の反復 を実行してください。1 回の電気熱計算後、シミュレーションされた JV 曲線が実験データに重ねて表示されます。これは ?? に示されています。 この場合、一致は良好であり、モデルがデバイスの支配的な 輸送および熱物理を捉えていることを示しています。
自己発熱の役割は、同じ比較を 熱モデルを無効にして繰り返すと明確になります：
- Thermal リボンに戻ります。
- thermal model を無効にします。
- 比較を再び 1 回の反復 実行します。
結果として得られる JV 曲線は、もはや実験データと一致しません。これは ?? に示されています。 このずれは、電流密度が大きく、 キャリア輸送による自己発熱 が顕著になる高バイアスで最も顕著です。
この比較は、電気熱モデリングの中心的な結果を示しています： デバイスが高電流密度領域で動作すると、電力散逸が 温度上昇を介して輸送にフィードバックします。 温度固定の電気シミュレーションではこの効果を捉えることができませんが、 結合電気熱モデルでは可能です。