GaAs/AlGaAsパラボリック量子井戸チュートリアル
1. はじめに
このチュートリアルでは、OghmaNanoにおける放物線型バンドモデルを用いて、単純なGaAs/AlGaAs量子井戸の電子的・光学的特性を計算する。幅の狭いGaAs層が、バンドギャップの広いAlGaAs障壁層の間に配置されている。このようなバンド端エネルギーの差異により、成長方向において電子および正孔が閉じ込められる一方で、層面に平行な方向には自由に移動可能となる。その結果、離散的な閉じ込めエネルギーが生じ、これらが二次元の電子サブバンドを形成するのである。
まずは既存のシミュレーションを開き、層構造を確認した上でモデルを実行し、計算された伝導帯および価電子帯のサブバンドを調べます。次に、これらの状態が吸収現象、誘導放出、正味の光利得、自然放出にどのように寄与するかを検討します。この分野の物理に不慣れな方は、放物線型バンド理論ページでより詳細な解説をご覧いただけます。
2. 量子井戸のサンプルを開く
OghmaNanoを起動し、ファイルリボンから「新規シミュレーション」を選択します。デバイス選択画面で「GaAsデモ」カテゴリをダブルクリックします（??）。これにより半導体および量子井戸に関する各種サンプルが表示されます。続いて「GaAs/AlGaAs QW parabolic」をダブルクリックします（??）。そうすると事前に設定済みのサンプルがメインのシミュレーションウィンドウに表示されます。
3. デバイス構造の理解
メインウィンドウには3層構造が表示されている（??）。GaAs層が井戸層であり、その両側にある2つのAlGaAs層が障壁層となっている。寸法を確認するにはメインウィンドウ左側の「レイヤー編集器」をクリックする。提供されている例では以下の構造が用いられている：
| 層 | 厚み | 目的 |
|---|---|---|
| アルミニウムガリウム砒素 | 40ナノメートル | 最初の障壁層 |
| ガリウム砒素（量子井戸） | 10ナノメートル | 量子井戸 |
| アルミニウムガリウム砒素 | 40ナノメートル | 第二の障壁層 |
中央の層には「アクティブ層」ではなく「量子井戸層」というタイプが明示的に設定されている点にご注意ください。各設定内容および光学材料の選択状況は??に記載されています。バンドギャップの大きい障壁層がキャリアの閉じ込めポテンシャルを形成しており、具体的なバンドオフセットの値は本例で用いられている材料パラメータに依存します。
なぜ量子井戸の幅を変更することが重要なのでしょうか？
成長方向におけるキャリアの閉じ込め作用により、電子および正孔の波動関数は特定の形態に限定される。幅が\(L\)である一次元無限井戸という最も単純なモデルにおいて、井戸底部から測定した電子の閉じ込めエネルギーは以下の通りである
\[ E_n=\frac{n^2\pi^2\hbar^2}{2m^{*}L^2},\qquad n=1,2,3,\ldots \]
ここで\(m^{*}\)はキャリアの有効質量を表し、\(\hbar\)は換算プランク定数です。実際の例ではアルミニウムガリウム砒素からなる障壁層が有限の幅を持つため、OghmaNanoは無限井戸モデルを用いるのではなく実際の材料プロファイルに基づいて閉じ込め状態を計算します。それでもなお、この式は有用な予測を与えてくれます。すなわち、ガリウム砒素の井戸層の幅を狭めると電子の閉じ込めエネルギーが一般に増加するのです。正孔にもそれぞれ固有の有効質量および閉じ込めエネルギーが存在します。
電気的パラメータボタンを使用することで、基礎となる材料関数を確認または編集することが可能です。スクリーンショットにおいて、AlGaAs材料の定義には組成および温度に依存するバンドギャップの計算処理が含まれています（??）。この入門チュートリアルを完了するにあたり、材料コードを編集する必要はありません。
4. シミュレーションの実行
前の画面に戻る ファイル リボンを引き、押す シミュレーションを実行 （青色の三角形ボタン）。 ターミナル このタブには、ソルバーが実行中に出力するメッセージが表示される。
提供されている例においては、ターミナルに固有状態のエネルギー値や繰り返し計算結果が報告される。
block converged after ... iterations メッセージ??）。計算が完了するまで待ってからプロットを開いてください。
ソルバーは各層における量子化された状態を計算し、次にそれらが面内波数ベクトルによってどのように変化するかを調べる。 \(k_{\parallel}\)光学的計算が設定されている箇所では、例示されている各種キャリア状態に対応するスペクトルも生成される。
5. 電子サブバンド：バンド構造図の解釈方法
計算されたバンド分散プロット（??）を開いてください。そこには明らかに二つの異なる曲線群が存在します。正のエネルギー値を持つ伝導帯サブバンド群と、価電子帯基準値付近またはそれ以下のエネルギー値を持つ価電子帯サブバンド群です。これらの曲線群間の垂直方向の隔たりは、半導体のバンドギャップならびに閉じ込められた電子・正孔のエネルギー値を反映しています。各グループ内でエネルギー値が異なるのは、量子井戸内における量子化状態に起因するものです。
井戸層内での運動は束縛されているものの、キャリアは層に沿って移動することも可能である。このため、束縛エネルギー値を平行方向の波数ベクトルに対してプロットすると、放物線形状となるのである。
\[ E_{e,n}(k_{\parallel})=E_{e,n}(0) +\frac{\hbar^2 k_{\parallel}^2}{2m_e^{*}}. \]
ここで\(E_{e,n}(0)\)は電子サブバンド\(n\)の底部におけるエネルギー値である。対応する価電子サブバンドにおいては、単純な放物線形状が下向きに湾曲する。
\[ E_{v,m}(k_{\parallel})=E_{v,m}(0) -\frac{\hbar^2 k_{\parallel}^2}{2m_h^{*}}. \]
これが本チュートリアルにおける基本的な近似条件である。すなわち、量子閉じ込め方向においては状態のエネルギーが量子化されるが、面内方向においては有効質量に基づく分散関係が維持されるのである。より高度な多バンドk·pモデルを用いれば、各バンド間の混成が生じ非放物線性が現れることも可能となる。 8バンド亜鉛ブレンデ理論のページ なぜそれが重要となるのかを説明している。
6. サブバンドから光スペクトルへ
電子は価電子帯および伝導帯のサブバンド間で光遷移を起こすことができる。通常の光入射条件下では、光子の運動量は電子の波数ベクトルに比べて非常に小さいため、面内波数ベクトルを保存する遷移のみを考慮する単純なモデルが用いられる。分散曲線の中心付近において、閉じ込められた電子状態\(n\)と正孔状態\(m\)間の遷移に要する光子エネルギーの近似値は以下の通りである。
\[ \hbar\omega\simeq E_g+E_{e,n}^{\mathrm{conf}}+E_{h,m}^{\mathrm{conf}}. \]
これら2つの閉じ込めエネルギーは、それぞれの井戸層におけるバンド端から上方へ測定される。利用可能な遷移過程、その始状態および終状態のキャリア占有状態、さらにはスペクトルの広がり具合が、計算された光応答の形状を決定する。これらの事象については、さらに詳しく調査することが可能である。 光遷移理論のページ.
6.1 誘導放出と正味利得
誘導放出の結果プロット（??）には、例で用いられた各種キャリア状態における波長ごとの発光寄与度が示されている。スペクトル上の段差や肩の部分は、各サブバンド間遷移が寄与するエネルギー値に対応している。また、別途出力される正味利得の結果（??）には誘導吸収の影響も含まれるため、その値は正にも負にもなり得る。
概略的には、確率\(f_c\)で占有されている伝導帯状態および確率\(f_v\)で占有されている価電子帯状態を考えると、正味の誘導放出寄与量にはこれらの占有因子が含まれる。
\[ \underbrace{f_c(1-f_v)}_{\text{stimulated emission}} -\underbrace{f_v(1-f_c)}_{\text{absorption}} =f_c-f_v. \]
正の正味誘導放出寄与を得るためには、遷移に関与する状態間で十分なキャリア反転が必要である。プロットされた全スペクトルは、光遷移の強度および利用可能な状態の数にも依存する。
6.2 吸収と自然放出
吸収特性のプロット（??）は、許容される光学遷移による吸収の寄与を示している。吸収が生じるためには、初期状態として価電子帯に電子が存在し、かつ伝導帯に遷移先となる空の状態が存在する必要がある。一方、自然放出のプロット（??）は、外部からの光子の誘導作用を必要とせず、伝導帯に存在する電子が空の価電子帯へと放射再結合する過程を表している。この自然放出における占有状態の計算には\(f_c(1-f_v)\)が用いられる。
これらのプロットにおいて、顕著なスペクトル特性が同様の波長領域に現れている点に注目してください。これは当然のことで、同じ電子・正孔サブバンド構造が可能な光遷移エネルギーを決定しており、キャリア占有状態の違いによって吸収現象と各種発光現象が区別されるからです。自然放出のプロットでは縦軸が対数スケールとなっているため、曲線の高さを比較する際には注意が必要です。
7. 計算結果ファイルの確認方法
シミュレーション終了後、メインウィンドウで出力タブを開いてください（??）。この例では、各種光学的計算結果が個別のCSVファイルとして出力される。ファイル名はgain_00.csvというパターンに従い、吸収、正味利得、自然放出、誘導放出といった接尾辞が付加されており、さらに別の計算条件に対応する同様のファイルも存在する。結果を確認するには該当ファイルをダブルクリックしてプロットウィンドウで表示すればよい。出力結果がCSV形式で保存されるため、OghmaNano外部でも解析や再プロットが可能である。
番号付きのファイル群を参照することで、単一の光学計算結果だけでなく、さまざまなキャリア状態における計算スペクトルの変化の様子を追うことができる。
8. 試してみよう：井戸層の幅を変更する
推奨実験：レイヤーエディターに戻り、GaAs（量子井戸）の厚みを10 nmから8 nmに変更し、AlGaAs障壁層はそのままにしておく。その後シミュレーションを再実行し、得られたサブバンドのエネルギー値および光スペクトルを元の例と比較してみる。
- 最も低エネルギーの電子閉じ込め準位は、GaAsの伝導帯端に対して上方へ移動するのでしょうか？
- 電子および正孔の遷移エネルギーはどのように変化するのか？
- 光スペクトルの開始波長は、より短波長側へ移動するのでしょうか、それとも長波長側へ移動するのでしょうか？
✅ 予想される物理的傾向
井戸層の幅が狭くなるほど、キャリアの閉じ込めエネルギーは増大する傾向にある。そのため最も低エネルギーのバンド間遷移はより高い光子エネルギーへと移動し、結果として波長が短くなる。ただしこのシフトの正確な度合いは、障壁層の高さや電子・正孔の有効質量に依存する。\(\lambda=hc/E_{\mathrm{photon}}\)によれば、遷移エネルギーが高くなるほど発光または吸収の波長は短くなる。
この例に慣れたら、次に比較すべきはGaAs/AlGaAs量子井戸の8バンド亜鉛スピネル型k·pモデルです。このモデルでは量子閉じ込め効果は維持されるものの、独立した放物線型バンドにとどまらず、伝導帯と価電子帯間の相互作用や混成も記述可能となります。