GaAs チュートリアル（パート A）: 垂直欠陥を持つ 3D ドリフト–拡散シミュレーションを実行する

1. はじめに

このチュートリアルでは、意図的に明瞭な 垂直欠陥（スタックを上から下へ貫通するシャント）を含む GaAs pn ダイオード の 3D ドリフト–拡散 シミュレーションを実行します。これは 1D モデルにとって意図的に「不公平な」問題です。 欠陥によって、純粋な膜厚方向シミュレーションでは存在しえない 横方向電流、電流集中、 および 局所的再結合 が強制されるためです。なぜこれが重要か: 実際のデバイスでは、 性能は理想的な層構造そのものではなく、実際に電流がどこを流れるか によって制限されることがよくあります。 局所欠陥、ピンホール、端部リーク、有限サイズのコンタクト、不均一な材料は、測定された JV を支配することがあります。 3D ドリフト–拡散モデルを用いると、JV 曲線の特徴を デバイス内の特定領域 に直接対応付けることができ、 1D では答えられない問いに答えられます。

自分の研究でこれをどう使うか: ここでの欠陥を、任意の 局在した非理想性 のテンプレートとして扱ってください。 たとえば、シャントフィラメント、ブロッキング層を貫通するピンホール、粒界、損傷したプローブ領域、有限面積の電極などです。 ワークフローは同じです。最小限の 3D 幾何構造を作成し、横方向メッシュを可能な限り小さく保ち、ソルバーを実行し、 その後 コンタクト分解 JV 曲線 と スナップショット （jn.csv、\(\phi\)、再結合）を用いて、電気的挙動を物理的位置へ対応付けます。

目的は二つあります。(1) OghmaNano における 真に 3D な電気モデリング への実用的な入り口を提供すること、 そして (2) モデリング判断、すなわち いつ 3D が本質的に必要か（欠陥、有限コンタクト、横方向変動） およびいつそれが 1D に縮退するかを教えることです。パート A では欠陥を含むデバイスを実行し、 パート B では欠陥を除去して、対称性が回復すると 3D/2D/1D が一致することを示します。 パート C では照明を有効にし、低電流において 数値的感度 がどこに現れるかを確認します。

2. 開始方法

OghmaNano のメインウィンドウから、New simulation をクリックしてデバイスライブラリを開きます （?? を参照）。 GaAs demos をダブルクリックします。これによりガリウムヒ素のサンプルフォルダが開きます （?? を参照）。 最後に、GaAs - 3D defect をダブルクリックし、書き込み権限のあるフォルダへシミュレーションを保存してください。

サンプルを開くと、メインウィンドウには デバイスの 3D 表示 が示されます （?? を参照）。 この構造は、意図的な 垂直欠陥領域（たとえば不純物フィラメントや短絡様の損失経路）を表す 二つの可視ブロックを含む GaAs pn ダイオード です。 欠陥はデバイスを 上から下まで 貫通しており、これにより本質的に 三次元 となり、 純粋な 1D モデルでは存在しない 電流集中 および 空間的に局在した再結合 を生じさせることができます。

材料および欠陥パラメータを確認するには、Device structure タブへ移動し、 Electrical parameters をクリックしてパラメータエディタを開きます （?? を参照）。 そこには次の個別のパラメータセットがあります:

• GaAs p+ および GaAs n+（ダイオード層）、そして
• defect0 および defect1（欠陥領域）。

この例では、物理的に妥当な GaAs に類似した移動度 および キャリア密度 を使用しています。 再結合には単純な二体 自由–自由 チャネル （有機デバイスの文脈ではしばしば n³–p³ と呼ばれます） が含まれており、その代表的な係数は 1 × 10⁻¹⁵ m³s⁻¹ です。

また、平衡 Shockley–Read–Hall（SRH）トラップ が有効になっていることにも気付くでしょう。 多くの有機デバイスシミュレーションでは、トラップ占有が 大きな電荷を蓄積して静電場へフィードバックするため、動的 SRH トラップ が必要です。 しかし、GaAs ではトラップ密度が通常十分低いため、トラップはしばしば純粋に 再結合損失機構 として扱うことができます。 平衡 SRH を使用することで、蓄積トラップ電荷の追加状態変数を導入せずに この物理を捉えることができ、モデルを安定かつ効率的に保つことができます。

3. シミュレーションを調べる: ドーピングとメッシュ

Electrical リボンは、シミュレーションを定義および診断するために用いられる中核的な電気エディタへのアクセスを提供します （?? を参照）。 このセクションでは二つのツールに注目します。 Doping/Ions は pn 接合を定義し、Electrical mesh は ドリフト–拡散計算の空間分解能を制御します。

Doping/Ions をクリックしてドーピングエディタを開きます （?? を参照）。 この例は意図的に単純な pn ダイオード であり、p 型領域と n 型領域が明確に分かれています。 ここでドーピングプロファイルを単純に保つことで、後で現れる異常な挙動を 電気静力学ではなく 幾何形状（欠陥）に帰属しやすくなります。

次に、Electrical mesh をクリックしてメッシュエディタを開きます （?? を参照）。 x、y、および z がすべて有効になっていることがわかります。つまり、このシミュレーションは真に 三次元 です。 この座標系では:

• y が主たる輸送方向（ダイオードの膜厚方向）であり、
• x および z は欠陥と横方向電流を分解するための横方向です。

y 方向 のメッシュは、n 型 層と p 型 層にまたがって意図的に分割されています。 これにより各領域で適切な分解能が維持されます。これは後で層厚を変更したり、 非対称性を導入したりする場合に重要になります。 一方、横方向メッシュ（x と z）では比較的少ない点数 （たとえば、x で約 10 点、z で約 5 点）を使用しています。 これは意図的です。3D では計算コストが急速に増大するためです。 各次元でメッシュ密度を係数 \(k\) だけ増加させると、未知数総数は概ね \(k^3\) に比例して増加します。 実際には、3D シミュレーションはメッシュを細かくするほど三乗的に遅くなります。 したがって、関心のある物理を解像するのに十分な範囲で、横方向分解能は可能な限り低く保つべきです。

4. 3D シミュレーションの実行と結果の表示

シミュレーションを実行するには、Run simulation（青い三角形）をクリックします。これは真に 三次元ドリフト–拡散 問題であるため、1D 計算よりも時間がかかります。 それなりに新しいノート PC であれば、このサンプルはおよそ 30 秒 で完了するはずです。 実行時間がこれより大幅に長い場合は、まずシミュレーションが ローカルドライブ に保存されていることを確認してください。ネットワークマウントやクラウド同期フォルダは、 大量のファイル I/O によって実行を大幅に遅くすることがあります。計算が終了したら、Output タブを開きます （?? を参照）。 このディレクトリには、ドリフト–拡散計算によって生成される標準出力ファイルが含まれており、 内部変数のスナップショットや電流–電圧データなどが含まれます。 この時点で、集約出力とコンタクト分解出力を区別することが重要です。

このデバイスは 有限面積コンタクトを持つ 3D シミュレーションであるため、一般には グローバルな jv.csv ファイルは無視し、代わりに コンタクト分解 JV 曲線 を使用すべきです。 jv_contact0.csv および jv_contact1.csv は、それぞれ 上側および下側コンタクトに関連する電流を報告します。 多次元シミュレーションでは、これらの曲線が電気的挙動を解釈し、 コンタクト固有の損失を診断するための最も物理的に意味のある方法です。jv_contact0.csv を開いて、 上側コンタクトの JV 曲線を表示してください （?? を参照）。 デバイスには垂直欠陥が含まれているため、この曲線は理想的な 1D pn ダイオード特性にはならないはずです。 次のセクション（パート B）では、欠陥を無効にし、JV 応答がどのように変化するか、 そして内在的な三次元非対称性が除去されると同じ構造がどのように 2D および 1D に縮退するかを確認します。

5. スナップショットビューアで内部変数を調べる

ソルバー内部の動作を確認するには、出力フォルダ内の snapshots/ ディレクトリを開きます。 これにより、シミュレーションステップ （通常は印加バイアス）に応じた空間分解されたシミュレータ変数へアクセスできます。 ディレクトリをダブルクリックすると、 ?? および ?? に示されているスナップショットビューアが開きます。

スナップショットウィンドウで、青い + ボタンをクリックしてプロットを追加し、 File to plot ドロップダウンメニューから jn.csv を選択します。 jn.csv は 垂直電子電流密度 に対応します。 メインスライダーを使用すると、印加電圧に応じて電流分布がどのように変化するかを確認でき、 一方で y および z スライダーを使うと 3D 体積のスライスを取ることができます。

この例は 暗状態のダイオード であるため、すべての電流は光電流ではなく電気的に駆動されています。 ここでの重要な観察点は、欠陥が 高度に局在した伝導経路 を生み出しており、 これが 3D で直接可視化され、非理想的な JV 挙動と明確に対応付けられることです。スナップショットビューアは、Ec.csv および Ev.csv（バンド端）、 準フェルミ準位、再結合率、キャリア密度など、他の多くの内部量を調べるためにも使用できます。 JV 曲線と並行してこれらの場を調べることは、 多次元デバイスにおいてドリフト–拡散ソルバーがどのように振る舞うかの直感を養う最も効果的な方法の一つです。

静電ポテンシャル \(\phi\) は印加バイアスに応じて変化しますが、 ドーピング分布（Nad.csv）はドーピングエディタによって固定されており、 電圧によって変化しない ことに注意してください。 それでも、3D でドーピングを可視化することは有用です。とくに後続のチュートリアルで 横方向ドーピング変動や空間的に局在した欠陥を扱う場合に重要です。