量子井戸モデルの選定方法
学習目標：このページを読み終えた後、以下のことができるようになっている必要があります。
- 各OghmaNanoモデルが適用可能な材料系を対応付けてください。
- 各モデルにおける伝導帯および価電子帯の取り扱い方を比較してみましょう。
- 精度と計算コストのバランスを考慮する必要がある。
- 対象となる物理的問題に応じて適切なモデルを選択すること。
- 最適なモデルは一般性だけで決まるのではなく、対象となる系に依存することを理解しておく必要がある。
前提条件：各モデルに関する説明ページを一読すること。本ページではそれらをまとめて解説する。
想定読了時間：20分。
OghmaNanoには複数の量子井戸モデルが用意されていますが、あらゆるシステムに最適なモデルは存在しません。スカラー型の放物線モデルは迅速な見積もりに適しており、正確な価電子帯の物理特性や利得の算出には8バンドのk·pモデルが必要です。また、間接遷移型バンドギャップを持つシリコン材料の場合にはシリコン専用モデルが求められます。本ページではこれらのモデルを比較することで、単に最も複雑なモデルを選択するのではなく、目的に応じて適切なモデルを選べるようにしています。
1. 比較の概要
| モデル | 材料 | 伝導帯 | 価電子帯 | 非放物線性 | ひずみ | バンド混合 | 計算コスト |
|---|---|---|---|---|---|---|---|
| 放物線型 | 任意の値（推定値） | 単一質量 | 単一質量 | なし | オフセットのみ* | 該当するテキストなし | 極めて低い |
| 8バンド亜鉛鉱型モデル | ガリウムヒ素、インジウムガリウムヒ素、リン化インジウム、ヒ化インジウム、ヒ素化ガリウム… | 非放物線的 | HH、LH、SOが結合している | はい | 完全なBir–Pikusモデル | 伝導帯・価電子帯および価電子帯モデル | 高い |
| 10バンド亜鉛鉱型モデル | 希薄窒化物半導体（GaInNAs） | BACおよび非放物線性 | HH、LH、SOの結合 | はい | 完全なBir–Pikusモデル | 伝導帯・価電子帯＋N＋価電子帯 | 高い |
| 8バンドウルツァイトモデル | GaN、AlN、InN、およびそれらの合金 | 異方性 | HH、LH、CHの結合状態 | はい | 完了 (\(a_i,D_i\)) | 伝導帯・価電子帯および価電子帯モデル | 高い |
| シリコンのΔバレー | シリコンおよびシリコンゲルマニウム（伝導帯） | 6つの異方性を持つバレー | — | バレー質量 | \(\Xi_d,\Xi_u\) 分裂 | ひずみによるバレー形成 | 低い |
| ルッティンガー・コーン6バンドモデル | シリコン、ゲルマニウム、SiGe（価電子帯） | — | HH、LH、SOの結合 | はい | 完全なビル・ピクスモデル | 価電子帯 | 中程度 |
この放物線型固有値ソルバーでは、バンドオフセットや（任意で）バンド端の剛体的なひずみシフトを考慮可能だが、ひずみによって生じる価電子帯の分裂は考慮されない。
2. 伝導帯の取り扱い
直接遷移型のIII–V族半導体材料においては、伝導帯はケイン結合を明示的に用いることで非放物線性を表現できる8バンドモデルによって適切に記述される。希薄窒化物半導体の場合、伝導帯はさらにバンド反交差現象によって分裂するため、10バンドモデルが必要となる。シリコンやSiGeにおいては、伝導帯の極小点は6つのΔバレーであり、これらは間接遷移型かつ異方性を持つため、いかなるゾーン中心用のk·pモデルでも記述できない。このためバレーソルバーの使用が不可欠である。放物線近似ソルバーでは伝導帯を単一の固定質量として扱うため、大まかなエネルギー準位の推定には適しているが、より高次のサブバンドや狭いバンドギャップの解析には不十分である。
3. 価電子帯の取り扱い
価電子帯においては、重い正孔、軽い正孔、およびスプリットオフ正孔が常に強く結合しているため、単一の有効質量では記述が極めて困難となる。したがって、正孔のサブバンド間隔、利得、あるいはTE/TM比といった定量的な計算を行うには、多バンドモデル、すなわち8バンドモデルが不可欠である。 亜鉛鉱型構造 または ウルツァイト III族化合物半導体用のモデル、あるいは 6バンドのラッティンガー・コーン模型 IV族材料用のモデルです。放物線型ホールモデルは、あくまで初期見積もり用としてのみ使用すべきです。
4. 計算コストと実用性
計算コストは基底関数の規模および\(k_\parallel\)点の数に比例して増加する。スカラーモデル（パラボリックモデルやシリコン・バレーモデル）では三対角行列を解くだけで済み、計算負荷は低い。一方、マルチバンドモデルでは各\(k_\parallel\)点において次元数が\(6N\)から\(10N\)に及ぶ複雑なエルミート行列を解く必要があるため、完全な分散関係を求めるにはかなりの計算コストがかかる。ただし、シフト反転固有値解法を用いれば重要な状態のみを算出できるため、実用的な範囲内に抑えられる。実務的には、まずスカラーモデルを用いてエネルギー準位を特定し構造を設定した上で、最終的な定量的計算には適切なマルチバンドモデルを採用するのが望ましい。
5. モデルの選択
- GaAs、InGaAs、InP、InAs、GaSbの量子井戸およびレーザー用モデルとしては8バンド亜鉛鉱型モデルが適しています。
- GaInNAsおよびその他の希薄窒化物半導体：10バンド亜鉛鉱型モデル。
- GaN／InGaN／AlGaN（窒化物系LEDおよびレーザー）：8バンドウルツァイトモデルを用い、TEモードおよびTMモードにおける利得値を算出する。
- シリコン／SiGeの伝導帯（輸送現象および量子井戸）：Δバレー固有値計算機。
- シリコン／ゲルマニウム／SiGeの価電子帯（正孔）：6バンドのラッティンガー・コーン模型。
- 簡易的なレベル推定や教育目的には、パラボリック型固有値ソルバーが適しています。
適切なモデルとは、対象となる物理現象を記述するのに必要最小限の要素のみを含む最も単純なモデルである。より複雑なモデルであっても、不適切なパラメータを用いたり想定範囲外の系に適用したりすれば必ずしも正確性が向上するわけではない。逆に、その適用範囲内で用いられる単純なモデルであれば、正確かつ高速な計算が可能となる。
要点
- どのモデルも普遍的に最適というわけではなく、各モデルは材料系や解析対象に応じて選定される。
- 多バンドモデルは、価電子帯や利得に関する定量的な計算を行う際に不可欠である。
- シリコンの場合、専用の伝導バンド用ソルバーおよびLuttinger–Kohn型価電子帯用ソルバーが必要となる。
- まずはスカラーモデルを用いて設定と概算を行い、最終的な結果を得る際にはマルチバンドモデルを使用してください。