GaN/AlGaN量子井戸：8バンド六方晶系k·pモデルチュートリアル
1. はじめに
このチュートリアルでは、OghmaNanoが提供する八バンド六方晶系k·pモデルを用いて、2つのAlGaN障壁層に挟まれたGaN量子井戸の電子状態および光学的応答を計算します。本例では、厚さ10nmのGaN量子井戸と、厚さ20nmのAlGaN層が2層存在します。まずは提供されている例題を実行し、その後で伝導帯および価電子帯のサブバンド分散状態と、それに伴う光スペクトルを調査していきます。
六方晶系構造と閃亜鉛鉱構造の違いはここで重要となる。両者とも半導体結晶構造であるが、六方晶系構造には優先的な六方晶軸が存在する。この結晶対称性の違いにより、価電子帯の分裂状態、電子状態間の結合の強さ、さらには異なる方向に偏光した光の強度も変化する。したがって、閃亜鉛鉱構造用の8バンド計算モデルを単純にGaN材料パラメータに適用することは不可能である。
2. 六方晶系八バンドモデルの例を開く
ファイルリボンから「新規シミュレーション」をクリックします。表示されたウィンドウで「GaAsデモ」(??)をダブルクリックします。カテゴリ名はGaAsとなっていますが、実際には窒化物半導体やシリコンの量子井戸モデル例も含まれています。次に表示されるウィンドウで「GaN 六方晶系 8バンド k·p」(??)をダブルクリックします。これにより事前に設定済みの材料層構造が読み込まれ、該当するデモ例が選択されます。
3. GaN/AlGaN量子井戸の確認
メインウィンドウには三層構造が表示されている（??）。実際の量子井戸計算においては、閉じ込め効果はこの積層構造の厚み方向に生じる。電子のエンベロープ関数は各層にわたって変化し、一方で分散プロットはそれらに平行な方向の運動を表している。
左側のツールバーから「Layer Editor」をクリックします。??に示すように、この例では以下のような構造が用いられています：
| 層 | 厚み | 種類 |
|---|---|---|
| アルガン | 20ナノメートル | 活性層／バリア層 |
| GaN（量子井戸） | 10ナノメートル | 量子井戸 |
| アルガン | 20ナノメートル | 活性層／バリア層 |
初回のシミュレーション実行時は、提供されている材料設定をそのまま使用してください。スクリーンショットでは層の厚さが示されていますが、Alのモル分率、ひずみ構成、および分極パラメータは表示されていません。特定の実験用ウェーハを再現する際には、例題の材料設定にてこれらを確認する必要があります。
4. なぜ8バンド六方晶系ハミルトニアンを使用するのか？
単純なパラボリックモデルでは、各キャリアに有効質量を割り当て、そのエネルギーが面内波数ベクトルの二乗に比例して変化すると仮定する。このモデルはバンド極小付近においては有用だが、波数ベクトルが増大するにつれて各種価電子帯状態が互いにどのように混成するかを記述することはできない。六方晶系結晶においては、結晶構造の対称性の違いからc軸方向の運動や光の偏光特性が面内方向のそれらとは異なる点にも注意が必要である。
この8バンドモデルでは、伝導帯に類似した軌道1つと、それぞれに2つのスピン成分を持つ価電子帯に類似した軌道が3つ存在する。概念的には、固有状態は次のような形をとる。
\[ \Psi_{n,\mathbf{k}_{\parallel}}(\mathbf{r})= e^{i\mathbf{k}_{\parallel}\cdot\mathbf{r}_{\parallel}} \sum_{\alpha=1}^{8}F_{n\alpha}(z)u_{\alpha}(\mathbf{r}), \]
ここで、関数\(F_{n\alpha}(z)\)は量子井戸内における電子の閉じ込め状態を表し、\(u_\alpha\)は8つのブロッホ基底状態を意味する。したがって、ある閉じ込め状態は複数の価電子帯軌道からの寄与を含み得る。このため、バンド端から離れた領域においては、各色の曲線を単一の重ホールまたは軽ホールというラベルで表すことは誤解を招く可能性がある。
六方晶系構造における価電子帯の配列は、結晶場相互作用およびスピン軌道相互作用の影響を受ける。簡潔な記法を用いると、この計算は以下の式を解くことに相当する
\[ H_{8\times8}\!\left(k_x,k_y,-i\frac{\partial}{\partial z}\right) \mathbf{F}_n(z)=E_n(\mathbf{k}_{\parallel})\,\mathbf{F}_n(z). \]
面内波数ベクトルは明示的に保持されている一方で、成長方向の波数ベクトルは微分演算子に置き換えられている点に注目してください。異なる\(k_{\parallel}\)値に対してこの固有値問題を繰り返すことで、下図のようなサブバンド分散のプロットが得られます。基底関数、結晶場項、および材料係数に関するより詳細な説明は、六方晶系におけるk·p理論の章に記載されています。
5. シミュレーションを実行する
ファイルリボンに戻り、「シミュレーション実行」をクリックしてください。その後「ターミナル」タブを開いて処理状況を確認します。??内では、ターミナルがWZ KP8ソルバーを認識し、疎行列の構築状況や収束情報を出力するとともに、サンプリングされた各波数ベクトルに対する計算を進めていきます。各点において電子状態および価電子帯状態が個別に解かれるのです。
複数成分からなるエンベロープ関数問題であるため、各バンドごとに独立したスカラー型のパラボリック井戸問題を解く場合よりも処理が複雑になる。結果として生成されるファイルを確認する際には、固有値計算およびそれに続く光学的計算が完了するまで待つ必要がある。
5.1 光遷移強度の確認方法
シミュレーション実行後半において、端末には光学行列要素の強度値が出力される（??）。通常のc面六方晶系構造においては、TEとはc軸に対して垂直な方向の電場成分を指し、TMは同軸に平行な成分を指す。これら二つの偏光は、特定の電子遷移に対して必ずしも同等の強度で結合するわけではない。
光学的遷移強度を考える際に有用な指標は、遷移行列要素の二乗値である。
\[ M_{cv}^{(\mathbf{e})}(\mathbf{k}_{\parallel})= \left\langle\Psi_{c,\mathbf{k}_{\parallel}}\middle| \mathbf{e}\!\cdot\!\mathbf{p}\middle| \Psi_{v,\mathbf{k}_{\parallel}}\right\rangle, \qquad S_{cv}^{(\mathbf{e})}\propto\left|M_{cv}^{(\mathbf{e})}\right|^2. \]
ここで \(\mathbf{e}\) は光の偏光状態を表し、\(\mathbf{p}\) は運動量演算子である。ある状態が電子バンド構造中に存在していても、その行列要素が小さい場合、特定の光偏光への寄与はごくわずかとなる。したがって、最終的な結果におけるTE/TMの相対値は、単にエネルギー準位の位置情報以上の物理的知見を与えるのである。
6. 計算された電子バンド構造の理解
各バンドプロットにおける横軸は面内波数ベクトル \(k_{\parallel}\) であり、量子井戸内の位置を示すものではない。縦軸は閉じ込められたサブバンドのエネルギー値である。\(k_{\parallel}=0\) において各曲線はそれぞれ離散的な閉じ込めエネルギーの値から始まる。右方向へ進むにつれて、キャリアが層に平行に移動する際に同一の電子状態がどのように分散していくかが示される。
6.1 伝導帯状態
伝導帯の出力結果を開き、??内の分散特性を確認してください。GaNのバンド端付近の異なるエネルギー値から始まる複数の量子閉じ込め準位が確認できます。波数ベクトルが小さい領域においては、伝導帯サブバンドに対して次の近似式が有用となります
\[ E_{c,n}(k_{\parallel})\simeq E_{c,n}(0) +\frac{\hbar^2 k_{\parallel}^2}{2m_{e,\parallel}^{*}}. \]
記号\(E_{c,n}(0)\)は量子閉じ込めエネルギーを表し、\(m_{e,\parallel}^{*}\)は平面内運動に対するキャリアの応答を記述している。この関係式は上昇傾向を示す曲線の説明に役立つが、完全な8バンド計算においては固定された放物線形状を仮定するのではなく、バンド間結合が考慮されている。
6.2 価電子帯状態
価電子帯のプロット（??）はより複雑な形状を示している。その上部の曲線はまず\(k_{\parallel}=0\)から上向きに立ち上がった後、下向きに曲がる。また、いくつかの下部の曲線は明らかに異なる曲率を持っている。このような形状は、結晶場分裂、スピン軌道相互作用、量子閉じ込め効果、およびバンド混合が相互に作用する六方晶系の価電子帯結合モデルにおいてのみ実現可能である。
この図は、すべての色付き曲線が波数ベクトルの全範囲にわたって常に特定のA帯・B帯・C帯に属しているという証拠と解釈すべきではない。確定的な帰属を行うためには、各曲線を構成する固有ベクトル成分の寄与度を調べる必要がある。同様に、定量解析目的で本モデルを利用する場合、中心からずれた特異な極値点についても、選定した材料やひずみパラメータと照らし合わせて検討することが望ましい。
7. 光遷移、誘導放出および吸収
電子バンド構造によって可能な光遷移のエネルギーが決定されるが、遷移強度は電子および正孔の波動関数、偏光状態、さらにはキャリアの占有状態にも依存する。同一の面内波数ベクトルにおいて、垂直遷移における光子エネルギーは
\[ \hbar\omega=E_{c,n}(\mathbf{k}_{\parallel}) -E_{v,m}(\mathbf{k}_{\parallel}). \]
フェルミの黄金律によれば、端末に表示される光学的行列要素が重要となる理由が説明される。模式的に表現すると、遷移の確率は以下の通りである
\[ W_{v\rightarrow c}^{(\mathbf{e})}= \frac{2\pi}{\hbar}\left|H_{cv}^{\prime(\mathbf{e})}\right|^2 \delta\!\left(E_c-E_v-\hbar\omega\right). \]
摂動行列要素 \(H'_{cv}\) は光の偏光状態および関与する電子状態に依存する。デルタ関数はエネルギー保存則を表しており、実際のスペクトルを構築する際には許容される遷移の寄与分をサンプリングされた状態全体にわたって合計し、さらに広がりを与える。その後、キャリア占有率が特定の波長において吸収が支配的か誘導放出が支配的かを決定する。便利な計算用係数として \(f_c-f_v\) が存在し、ここで \(f_c\) および \(f_v\) はそれぞれ伝導帯および価電子帯における電子の占有数である。
GaN量子井戸は電子ギャップが広いため、算出されたスペクトルは近紫外領域に現れ、提供された例では主に330～370 nm付近に分布している。波長と光子エネルギーの関係は\(E_{\gamma}=hc/\lambda\)で与えられる。
7.1 誘導放出と正味の光増幅
誘導放出の計算結果図 (??) には、さまざまなキャリア状態におけるスペクトル群が示されている。正味の利得の算出 (??) では、誘導放出と吸収の競合が考慮されている。正の値を示す正味利得は、この電子モデルにおいて光増幅が生じていることを示唆するが、それだけではレーザー共振器が閾値に達しているかどうかは判断できない。閉じ込め係数や共振器の損失もまた考慮する必要がある。
7.2 吸収
吸収特性の結果は ?? に示されている。そのスペクトル端や特徴は、占有されている価電子帯状態と利用可能な伝導帯状態間のエネルギー差、ならびにそれらに対応する光学的行列要素の重み付けによって決定される。六方晶系の量子井戸においては、各状態の軌道成分がそれぞれの偏光モードとの結合強度を決定するため、TEおよびTMモードに対応するスペクトルも互いに異なる場合がある。
これらのスペクトルに現れる特徴は、設定されたモデルおよび広がり効果の結果である。それらを単純に特定の励起子ピークと同一視すべきではない。なぜなら、単一粒子のk·p理論に基づくバンド間計算だけでは、電子と正孔がクーロン力で束縛された励起子の挙動を完全に記述することはできないからである。
8. 自然放出
誘導放出とは異なり、自然放出は遷移を引き起こすための入射光子が存在しなくても発生する。この現象は伝導帯内に電子が存在し、価電子帯内に空の状態（正孔）があることに依存する。したがって、その電子占有率は以下に比例することになる
\[ R_{\mathrm{sp}}(\omega)\ \propto\ \sum_{c,v,\mathbf{k}_{\parallel}} \left|M_{cv}(\mathbf{k}_{\parallel})\right|^2 f_c(1-f_v)\,L\!\left(E_c-E_v-\hbar\omega\right), \]
ここで\(L\)は選択されたスペクトル広がりを表し、比例関係においては光子に関する前因子は省略されている。その結果は??に図示されている。対数スケールを用いることで、何桁にもわたる範囲のスペクトルを比較することが可能となる。スペクトルの短波長側に現れる光学的特性や、各条件における全体的な強度の変化にも注目されたい。
9. 自分で試してみよう
💡 参考になる計算例が2つある：
- GaN量子井戸の幅を変更してみましょう。中央の10nm厚の層を7nmに減らし、障壁層の厚さはそのままにしてみてください。その後、固有状態を再計算し、最も低エネルギーの伝導帯準位および光遷移エネルギーを比較してください。量子閉じ込め効果が強くなると一般に準位間のエネルギー差も増大しますが、価電子帯の応答は複数のバンドに依存するため詳細な挙動は複雑になります。
- TEおよびTMの遷移強度を比較せよ。ターミナルから出力される光学強度の値を用いて、特定の偏光方向を強く選好する遷移を特定すること。また、すべての遷移がスペクトルに均等に寄与すると仮定するのではなく、この強度の違いが価電子帯状態の特性とどのように関連しているかを検討すること。
✅ 何に注目すべきか？
ウェル幅が狭くなると、電子の閉じ込めエネルギーが上昇するとともに正孔のエネルギー準位も変化し、可能なインターバンド遷移エネルギーがシフトする。六方晶系結晶における波動関数の軌道成分の違いを反映して、TE/TM比は状態の組み合わせごとに大きく変動しうる。定量的な材料解析を行う際には、合金組成やひずみ、分極場の影響も併せて検討する必要がある。