GaAs/AlGaAs量子井戸：8バンド亜鉛鉱型結晶構造におけるk・pモデルチュートリアル
1. はじめに
このチュートリアルでは、OghmaNanoの8バンド亜鉛鉱型k・pモデルを用いてGaAs/AlGaAs量子井戸の電子的・光学的特性を計算します。まず既存のサンプルを確認し、層構造を調べた後、シミュレーションを実行して得られた伝導帯および価電子帯のサブバンドを検討します。最後に、これらの状態から算出された光スペクトルについても見ていきます。
このパラボリック量子井戸チュートリアルでは、比較的単純な有効質量近似に基づくバンドモデルを用いて電子および正孔の振る舞いを扱っている。ここで用いられるハミルトニアンには、伝導帯の基底関数に加え、重い正孔（HH）、軽い正孔（LH）、およびスピットオフ（SO）からなる価電子帯の基底関数が含まれ、各基底関数には2つのスピン成分が存在するため、合計で8つの基底状態が存在する。これら状態間の結合作用により、本モデルは個別のパラボリックモデルでは記述不可能な現象、特に価電子帯混成現象や面内波数ベクトルに応じたバンド曲率の変化をも記述することができる。
2. 8バンドモデルの例を開く
OghmaNanoにおいて、ファイルリボン上の「新規シミュレーション」をクリックします。??フォルダ内の「GaAsデモ」をダブルクリックします。次に、量子井戸関連のサンプル群から「GaAs/AlGaAs QW 8 band k.p」をダブルクリックします。??これにより事前に設定済みのシミュレーションが起動され、8×8のハミルトニアンを手動で入力する必要はありません。
3. 量子井戸構造および材料パラメータ
メインウィンドウには、3層構造のAlGaAs/GaAs/AlGaAs構造が表示されている（??）。左側のボタンを使用してレイヤーエディターを開く。提供されている例では、厚さ40nmのAlGaAs障壁層2枚が厚さ10nmのGaAs井戸層を挟み込む構造となっている。
| 層 | 厚み | 層の種類 |
|---|---|---|
| アルガサス | 40ナノメートル | 活性障壁層 |
| GaAs（量子井戸層） | 10ナノメートル | 量子井戸 |
| AlGaAs | 40ナノメートル | 活性障壁層 |
GaAsのバンド端は周囲のAlGaAsのバンドギャップ内に位置しており、層成長方向における量子閉じ込め効果をもたらす。レイヤエディタにおいて、中央のGaAs層は特に量子井戸として明示的にマークされている（??）。本例では、光学材料の設定および層の種類も既に構成済みである。
8バンドソルバーでは、スカラー型の放物線バンドモデルよりも多くの材料情報が必要となります。すなわち、バンド端や有効質量に加えて、バンド間結合を記述するパラメータ、スプリットオフエネルギー、そして適用可能な場合はひずみの情報もハミルトニアンに用いられます。GaAsに関する関数情報は、電気的パラメータ機能および材料用コードエディター（??）から確認可能です。スクリーンショットには価電子帯の変形ポテンシャルであるbおよびdが示されていますが、この入門的なシミュレーションにおいてはこれらを変更する必要はありません。
8バンド計算において何が変化するのか？
放物線近似モデルでは各バンドごとに個別の曲率が与えられる。一方、k・p法では対角要素以外の行列要素が基底状態間の結合を担う形で行列が構築される。模式的に表すと、閉じ込め状態の問題は以下のようになる。
\[ H_{8\times8}(k_x,k_y,-i\partial_z)\,\Psi_n(z) =E_n(k_x,k_y)\,\Psi_n(z). \]
ここで\(z\)は閉じ込め方向を表し、\(k_x\)および\(k_y\)は層に平行な運動を記述する。エンベロープ波動関数\(\Psi_n\)は8つの成分から成り、これらが各基底状態の寄与度を表している。本チュートリアルでは、完全な行列演算ではなく計算された固有エネルギーに焦点を当てる。詳細な説明は8バンド閃亜鉛鉱構造理論の章に記載されている。
4. シミュレーションの実行
ファイルリボンに戻り、青色の「シミュレーション実行」ボタンをクリックしてください。8バンドハミルトニアンが構築され、一連の面内波数ベクトルにおいて固有状態が計算されます。計算の進行状況はターミナルタブ（??）で確認可能です。
例示されているターミナル画面には、スパース行列の構築状況や固有エネルギーの値、そして block converged after 12 iterations のようなメッセージが表示されている。伝導帯および価電子帯の状態は、連続する波数ベクトルごとに個別に計算される。計算が完了し、ボタンが再び「Run simulation」に戻るまで待ってから、出力結果を確認するようにしてください。
5. 計算されたバンド構造の理解
このソルバーは、各面内波数ベクトルにおいて、閉じ込められた固有状態の離散的な集合を計算する。 \(k_{\parallel}\)この図はしたがって、以下の内容を示している。 サブバンド分散層間での空間的なバンド端分布ではなく、
OghmaNanoを起動してください。 出力結果
タブを開いて中身を確認してください bands フォルダを開くか、表示されているバンド概要ファイルを開いてください ??.
5.1 伝導帯サブバンド
伝導帯のプロット（??）には、上方に湾曲する2つの分枝が示されている。\(k_{\parallel}=0\)付近において、垂直方向の隔たりはそれぞれ異なる閉じ込め状態に起因する。また\(k_{\parallel}\)に伴うエネルギーの増加は、キャリアが閉じ込められていない井戸層に沿った運動を反映している。
伝導サブバンドの極小点付近において、その分散関係は次式で近似できる
\[ E_n(k_{\parallel})\simeq E_n(0) +\frac{\hbar^2 k_{\parallel}^{2}}{2m_{n,\parallel}^{*}}. \]
多バンドモデルにおいてこれは局所的な近似に過ぎず、曲率、ひいては見かけ上の有効質量はプロットされた波数ベクトル範囲全体にわたって一定である必要はない。伝導帯と価電子帯間の結合がこの挙動に寄与している。
5.2 価電子帯サブバンド：混成現象と曲率
価電子帯のプロットの方がより興味深い（??）。表示されている2つの分枝はそれぞれ異なるエネルギー値から始まり、下向きに湾曲しており、\(k_{\parallel}\)の値において互いに交わるように見える。亜鉛ブレンデ型結晶構造における価電子帯では、重ホール状態と軽ホール状態はまったく独立した2つの放物線のように振る舞わない。有限な波数ベクトルにおいては、ハミルトニアン内の項がこれらの状態の性質を混成させるのである。
これを表現するための有用な方法として、行列要素によって結合している2つの近接した状態を考慮することが挙げられる \(V\)。
\[ H_{\mathrm{illustrative}}= \begin{pmatrix}E_{\mathrm{HH}}(k)&V(k)\\ V^{*}(k)&E_{\mathrm{LH}}(k)\end{pmatrix}. \]
もしこれらの状態間で結合が許される場合、互いに近接するときにそれぞれのエネルギー値が反発し合い、HH状態とLH状態の性質が各分岐間で入れ替わる可能性がある。また、対称性により結合がゼロとなる場合には実際の交差点が生じることもある。ただし、現在のスクリーンショットだけではこれらのケースを区別することはできない。状態の性質を特定するためには対応する固有ベクトルを調べ、必要に応じて見かけ上の交差点付近でさらに多くの点を計算する必要がある。上記の2×2行列は単なる例示用であり、完全な8バンドハミルトニアンではない。
5.3 伝導帯および価電子帯状態を同時に表示する
合成されたプロット図（??）では、電子状態および価電子状態が共通のエネルギースケール上に示されている。この例においては、価電子状態はエネルギー基準値付近に位置し、閉じ込められた伝導帯状態は約1.5 eV以上の領域に存在する。任意の2つの状態間のエネルギー差は、対応する光学行列要素やキャリア占有状態に応じて、垂直方向の光遷移に伴う光子エネルギーとなる。
光子の運動量は比較的小さいため、光遷移においては面内波数ベクトルが保存されるという近似が有効となる。
\[ \hbar\omega=E_{c,n}(k_{\parallel})-E_{v,m}(k_{\parallel}), \qquad \lambda=\frac{hc}{\hbar\omega}. \]
したがって、異なる電子状態と正孔状態の組み合わせがそれぞれ異なる波長において寄与することになる。また、多バンド波動関数は、単にエネルギー値からだけでは判断できない各遷移の強度や偏光特性も決定するのである。
6. 量子井戸の光学特性
8バンド計算により、光学モデル用の電子エネルギーおよび混成状態の波動関数が得られる。この例では、誘導放出、吸収、正味利得、自発放出それぞれに対するスペクトルが生成される。各プロット内の複数の曲線は、シミュレーションで用いられたキャリア状態の順序に対応するものであり、8つの個別なハミルトニアン帯域を表すものではない。
6.1 誘導放出と吸収
誘導放出とは、入射した光子によって電子が満たされている伝導帯から空いている価電子帯へと遷移する現象である。吸収はその逆の過程であり、価電子帯の電子が光子を吸収して利用可能な伝導帯へ移動する。これらの過程における電子の占有状態は、概念的には以下のように表される。
\[ \begin{aligned} \text{stimulated emission:}&\quad f_c(1-f_v),\\ \text{absorption:}&\quad f_v(1-f_c). \end{aligned} \]
誘導放出の曲線（??）および吸収の曲線（??）は、共に同じバンド間遷移エネルギーに基づいているが、キャリア占有状態の変化に伴いその強度は異なる。得られた結果では、スペクトル構造の大部分が790～850 nm付近に現れており、具体的な特徴は設定されたサブバンド、行列要素、キャリア密度および線幅広がりに依存する。
6.2 正の純利得：発光が吸収を上回る場合
正味の利得のグラフ（??）は、競合する誘導放出と吸収の寄与を合わせたものである。これらの占有確率を差し引くと、極めて単純な結果が得られる：
\[ f_c(1-f_v)-f_v(1-f_c)=f_c-f_v. \]
したがって、関与する伝導帯状態の電子占有数が対応する価電子帯状態よりも多い場合、この遷移は利得を生み出すことになる。負の値を示す純利得係数は吸収が支配的な状態を表し、正の値は刺激放出が吸収を上回る状態に対応する。ここで示されている純利得係数は量子井戸における光学的な結果であり、閉じ込め効果やその他の光学的損失にも依存する実際のレーザー共振器におけるモード利得と同一視すべきではない。
6.3 自然放出
自然放出現象においては、遷移を引き起こすための入射光子は必要とされない。これは、電子が占有している状態から利用可能な価電子帯状態へ放射的に再結合する際に生じる現象である。このときの占有因子もまた\(f_c(1-f_v)\)となるが、自然放出と誘導放出ではスペクトルの前係数や物理的意味合いが異なる。??に示されているスペクトルでは縦軸が対数スケールとなっているため、縦方向の等間隔が必ずしも放出される光子数の絶対的増加量が等しいことを意味するわけではない。
遷移強度、キャリア占有率、および線幅の拡がりについてより詳細な説明は、光学的遷移および利得理論のページをご参照ください。
7. 出力ファイルの確認方法
bandsおよびgainディレクトリにはその他の計算結果が格納されています。シミュレーション終了後、メインウィンドウで「Output」を選択してください（??）。表示される結果の中には、バンド構造図用のEc.csv、Ev.csv、Ec_Ev.csvのほか、gain.csv、absorption.csv、net_gain.csv、spontaneous.csvも含まれます。bandsおよびgainフォルダーには、さらに詳細なソルバー出力データが格納されています。CSVファイルをダブルクリックするとOghmaNano上でプロットが開きますし、数値データを他の解析用にエクスポートすることも可能です。
8. 自分で試してみよう：放物線近似モデルと比較してみる
推奨される実験として、パラボリック型GaAs/AlGaAs例をこの8バンド計算と並行して実行してください。この際、GaAs井戸層の厚みは10 nmに保ち、使用材料やその他の設定が同等であることを確認してください。その後、得られたバンド構造図および光学スペクトルを比較してください。
- \(k_{\parallel}=0\)付近において、最も低エネルギーの伝導帯サブバンドはほぼ放物線形状を呈しているでしょうか？
- 価電子帯の分散特性は、個々のパラボリックなバンドから成るモデルとどのように異なるのでしょうか？
- 計算された光学的閾値や相対的なスペクトル特性はどのように変化するだろうか？
次の実験として、GaAs井戸層の厚みを10nmから8nmに変更し、40nm厚のAlGaAs障壁層の厚みはそのままにしておきます。その後モデルを再実行し、閉じ込められた状態の固有エネルギーおよび光遷移が生じる波長範囲を比較してください。
✅ 確認すべき点
バンド極小点付近では、8バンドモデルによる計算結果が有効質量近似に基づく放物線状の分布と類似する場合がある。極小点から離れるにつれて、特に価電子帯においては混成効果や非放物線性がより明確に現れる。井戸層の幅が狭くなると一般に閉じ込めエネルギーが増大し、低エネルギー側のバンド間遷移がより高い光子エネルギー（より短い波長）へとシフトする。この変化の度合いは材料パラメータや結合した固有状態に依存する。