量子井戸とキャリアの閉じ込め

量子井戸とは何か？

量子井戸とは、その周囲よりもポテンシャルエネルギーが低くなっている空間領域であり、このエネルギー障壁によって電子のような粒子が閉じ込められる。井戸の幅が電子の波長と同程度になると、量子力学的な効果が特に顕著に現れる。

??に示すように、典型的な量子井戸の構造が描かれている。半導体デバイスにおいては、電子のバンドエッジが異なる2種類の材料の間に薄い層を成長させることでこのような構造を作製することができる。よく見られる例としては、AlGaAsで形成されたポテンシャル障壁に囲まれたGaAs量子井戸が挙げられる。

無限のポテンシャル障壁を持つ量子井戸をシュレーディンガー方程式で記述する

ポテンシャル井戸の幅が電子の波長と同程度になると、電子を明確な位置を持つ古典的な粒子として記述することはもはや不可能となる。その代わりに、電子の波動的性質や井戸内の各位置で電子が見つかる確率を記述する必要がある。

このため、時間に依存しないシュレーディンガー方程式を解く必要がある。

\[ \hat H\psi(x)=E\psi(x). \]

ここで、\(\psi(x)\)は電子の波動関数、\(E\)はそのエネルギー、\(\hat H\)は電子の運動エネルギーおよびポテンシャルエネルギーを記述するハミルトニアンです。ある位置に電子が存在する確率密度は\(|\psi(x)|^2\)によって与えられます。

単純な一次元量子井戸においては、シュレーディンガー方程式は次のように表される：

\[ \boxed{ -\frac{\hbar^2}{2m^*} \frac{d^2\psi(x)}{dx^2} +V(x)\psi(x) =E\psi(x). } \]

この方程式は一見複雑に見えるかもしれませんが、各項はそれぞれ明確な物理的意味を持っています。

• \(V(x)\)は、量子井戸およびその周囲のポテンシャル障壁のポテンシャルエネルギー分布を表している。
• \(\psi(x)\)は電子の波動関数であり、その二乗値は位置\(x\)において電子が存在する確率密度を表します。
• \(E\) これがその特定の量子状態に伴うエネルギーである。
• \(m^*\) これが電子の有効質量であり、半導体内で電子が力にどのように反応するかを示す指標となる。最初の項は電子の運動エネルギーを表している。

ポテンシャル分布 \(V(x)\) を与え、この方程式を解くことで、許容される電子のエネルギー準位およびそれに対応する波動関数が得られる。これらがまさに ?? に示されている結果、すなわち水平なエネルギー準位および各状態に付随する色付きの波動関数である。次の節では、単純な量子井戸に対してこれらの解をどのように算出するかを説明する。

この方程式は現実において何を意味するのか？

シュレーディンガー方程式の解を直感的に理解するためには、両端が固定されたバイオリン弦を例に考えるのが有効です。このような弦は、離散的な調和振動モードで振動することができます。基本調和振動モードにはピークが1つだけ存在しますが、より高次の調和振動モードではピークや節点が複数現れます。各調和振動モードはそれぞれ固有の振動周波数を持ちます。

量子井戸内に閉じ込められた電子にも同様の原理が適用される。量子井戸の寸法によって電子は特定の許容される波動関数群に制限され、それぞれが固有のエネルギー準位に対応する。ヴァイオリン弦が基本調和振動およびそれより高次の調和振動を持つのと同様に、量子井戸にも最低エネルギー状態とそれより高いエネルギーを持つ一連の状態が存在するのである。

これらの状態のうち最初の3つが??に示されている。点線の水平線は許容されるエネルギー値\(E_1\)、\(E_2\)、\(E_3\)を示しており、色付きの曲線はそれぞれに対応する波動関数を表している。最もエネルギーの低い状態は最も単純な形状をしており、エネルギーが高くなるにつれて波のパターンも次第に複雑になっていく。

では、これらの波動関数は実際に何を示しているのでしょうか？振動するバイオリン弦の変位とは異なり、電子の波動関数は確率的な意味を持ちます。その絶対値の二乗は、特定の位置に電子が存在する確率密度を表すのです。

\[ \boxed{ P(x)=|\psi(x)|^2 } \]

例えば、最もエネルギーの低い状態にある電子は、量子井戸の中心付近に存在する確率が最も高くなる。より高いエネルギー状態では確率分布が複雑になり、確率がゼロとなる節点によって隔てられた複数のピークが現れる。また、波動関数は周囲のポテンシャル障壁内にもわずかに広がっているため、電子が量子井戸の外側に存在する可能性もごくわずかながら存在する。

したがって、シュレーディンガー方程式を解くことで、電子が取り得るエネルギー準位と、それぞれのエネルギー準位に伴う確率密度分布の両方が明らかになるのである。

単純な解法：無限高さの量子井戸

では、量子井戸内の電子の波動関数を計算してみましょう。計算を簡略化するため、まずは井戸を取り囲むポテンシャル障壁が無限に高いと仮定します。この場合、電子は井戸の外側には存在できず、波動関数は井戸の両端でゼロにならなければなりません。

幅\(L\)の量子井戸を考え、その両端をそれぞれ\(x=0\)および\(x=L\)とする。井戸内部ではポテンシャルエネルギーをゼロと設定するため\(V(x)=0\)となる。まずシュレーディンガー方程式から始める。

\[ -\frac{\hbar^2}{2m^*} \frac{d^2\psi(x)}{dx^2} +V(x)\psi(x) =E\psi(x), \]

ポテンシャルエネルギー項を取り除くと、次式が得られる：

\[ -\frac{\hbar^2}{2m^*} \frac{d^2\psi(x)}{dx^2} =E\psi(x). \]

さて、この方程式を満たす波動関数を見つける必要がある。適切な出発点としては、正弦波解を仮定することだ。

\[ \psi(x)=A\sin(kx). \]

ここで\(A\)は波動関数の振幅を決定し、\(k\)は位置に応じて波動関数がどの程度速く振動するかを示します。\(k\)の値が大きいほど波長は短くなり、同じ距離内により多くの振動が含まれることになります。

この関数が正しい解であるかどうかを確認するため、2回微分してみる。

\[ \frac{d^2\psi(x)}{dx^2} =-k^2A\sin(kx) =-k^2\psi(x). \]

推測した関数を二度微分すると、元の関数に定数を掛けた形に戻ることがわかる。この性質こそが、量子井戸内部におけるシュレーディンガー方程式が求める数学的挙動なのである。

仮定した関数を方程式に代入することで、対応するエネルギー値が以下の条件を満たす場合にのみ、正弦波状の波動関数が妥当な解となることが確認できる。

\[ E=\frac{\hbar^2k^2}{2m^*}. \]

まずはエネルギー値の話に戻りますが、その前により興味深い問題があります。すなわち、これらの正弦波解のうちどれが実際に量子井戸内に存在し得るのか、という点です。

どのような波長の波が量子井戸内に存在し得るのか？

シュレーディンガー方程式の正弦波解が得られたものの、すべての正弦波が量子井戸内で物理的に存在し得るわけではない。ここで注意すべき点は、ポテンシャル障壁の高さが無限大であるため、波動関数は井戸の両端において必ずゼロにならなければならないということである。

これはまさにバイオリンの弦の振る舞いと同様である。??に示すように、両端が固定された弦は離散的な調和振動数で振動することができる。基本調和振動ではピークが1つ、第2調和振動では2つ、第3調和振動では3つと、それぞれの調和振動において弦の変位が両端でゼロになるという条件が満たされるのである。

無限に広がる量子井戸の場合も、まさに同様の数学的制約が適用される。許容される波動関数は、\(x=0\)と\(x=L\)の間に収まり、両端の境界においてその値がゼロでなければならない。

我々が仮定した正弦波型の波動関数は、\(\sin(0)=0\)であるため最初の境界条件を既に満たしている。二番目の境界条件を満たすためには、以下が必要となる：

\[ \psi(L)=A\sin(kL)=0. \]

正弦関数は、その引数が\(\pi\)の整数倍である場合にのみゼロとなります。したがって、許容される値は以下の通りです：

\[ kL=n\pi, \qquad n=1,2,3,\ldots \]

これを整理すると次のようになる：

\[ \boxed{ k_n=\frac{n\pi}{L}, \qquad n=1,2,3,\ldots } \]

したがって我々は重要な事実を明らかにした。シュレーディンガー方程式は正弦波解を許容するものの、量子井戸の境界条件により許される波長は特定の値に限定されるのである。これらはまさにバイオリン弦の倍音と類似している。

対応する波動関数は次のように表せる：

\[ \boxed{ \psi_n(x)= \sqrt{\frac{2}{L}} \sin\left(\frac{n\pi x}{L}\right), \qquad 0 \lt x \lt L } \]

正弦関数の前に付く係数は単に波動関数を規格化するためのものであり、これにより井戸内部のどこかに電子が存在する確率の総和が1となる。\(n\)の各値に対して、それぞれ異なる許容される確率分布が得られる。

エネルギー準位の計算

量子井戸内に存在しうる波動関数がどのようなものかを理解したところで、次はそれらが持つエネルギーについて考察していきましょう。

バイオリンの例えと関連性が見出せる。バイオリン弦の高次倍音は周波数が次第に高くなるが、同様に量子井戸における高次の状態も電子のエネルギーが次第に高くなることを意味する。

先ほど正弦波解をシュレーディンガー方程式に代入した際に得られた関係式を用いることで、これらのエネルギー値を計算することができる。

\[ E=\frac{\hbar^2k^2}{2m^*}. \]

既に説明した通り、\(k\)として許容される値は特定のもののみである。

\[ k_n=\frac{n\pi}{L}. \]

これをエネルギー式に代入すると次式が得られる：

\[ E_n= \frac{\hbar^2}{2m^*} \left(\frac{n\pi}{L}\right)^2. \]

式を整理すると、無限量子井戸における許容エネルギー値に関するよく知られた結果が得られる。

\[ \boxed{ E_n=\frac{n^2\pi^2\hbar^2}{2m^*L^2} } \]

この結果から、エネルギーは状態番号の二乗に比例して増加することがわかる。 \(n\)つまり、2番目の状態のエネルギーは1番目の状態の4倍であり、3番目の状態のエネルギーは9倍になるのである。

この結果から半導体工学において特に有用な知見が得られる。すなわち、エネルギー値は\(1/L^2\)に依存するため、量子井戸の幅を狭めるだけで閉じ込めエネルギーを増大させることが可能となる。また電子の有効質量も重要な要素であり、他の条件が同一であれば有効質量が小さいほどエネルギー準位は高くなる。

例えば、幅が\(12\,\mathrm{nm}\)で電子の有効質量が\(m^*=0.067m_0\)の無限量子井戸において、最も低いエネルギー準位はおよそ\(39\,\mathrm{meV}\)となる。したがって、許容される最初の3つのエネルギー準位は以下の通りである：

\[ E_1\approx39\,\mathrm{meV}, \qquad E_2=4E_1, \qquad E_3=9E_1. \]

これまでに、許容される波動関数およびそれらに対応するエネルギー準位を算出してきました。ただし、ポテンシャル障壁の高さが無限大であるという仮定は数学的な単純化に過ぎません。図??に示すより現実的な量子井戸を記述するためには、次に障壁の高さが有限である場合にどのような現象が生じるかを考察する必要があります。

無限高さのポテンシャル障壁ではない場合、どのような現象が生じるのか？

実際の半導体量子井戸においては、周囲のポテンシャル障壁の高さは有限である。このことにより重要な違いが生じる。すなわち、電子の波動関数はもはや井戸の端でゼロにならなくてもよく、たとえ電子のエネルギーが障壁の高さよりも低い場合であっても、周囲の物質中へわずかに広がることが可能となるのである。

この現象も同じシュレーディンガー方程式を用いて記述できるが、ここでは二つの異なる領域を考慮する必要がある。すなわち、量子井戸自体の領域である\(V(x)=0\)と、周囲のポテンシャル障壁領域である\(V(x)=V_0\)である。量子井戸の中心が\(x=0\)にある場合、ポテンシャルは次のように表される：

\[ V(x)= \begin{cases} 0, & \lvert x\rvert \lt \dfrac{L}{2}, \\[6pt] V_0, & \lvert x\rvert \geq \dfrac{L}{2}. \end{cases} \]

Let us consider an electron with an energy below the barrier height, so that \(0

重要な違いは、波動関数が二つの領域において異なる振る舞いを見せる点です。量子井戸内部では、無限井戸モデルの場合と同様に波動関数は振動を続けます。一方、ポテンシャル障壁領域では波動関数は直ちにゼロにはならず、指数関数的に減衰していきます。

これは、ポテンシャルエネルギーが電子のエネルギーを上回る場合にシュレーディンガー方程式が変化するために生じる現象です。波動関数がポテンシャル障壁内で減衰していく速度は次式で表されます：

\[ \boxed{ \kappa = \frac{\sqrt{2m^{*}(V_0-E)}}{\hbar} } \]

許容されるエネルギー準位を算出するためには、量子井戸内部における振動する波動関数が障壁領域において指数関数的に減衰する波動関数へと滑らかに接続するような解を見出す必要がある。この作業は通常、数値計算によって行われる。

その結果は??に示されている。色付きの波動関数が量子井戸の端を超えて広がっている点に注目されたい。波動関数の二乗値が確率密度を表すため、これは周囲のポテンシャル障壁内に電子が存在する確率がごくわずかであることを意味している。

電子が実質的にやや広い領域に閉じ込められるため、そのエネルギー準位は無限井戸近似によって予測される値よりも一般に低くなる。さらに、有限井戸では存在し得る束縛状態の数も有限に限定される。

OghmaNanoにおける量子井戸の計算方法

単純な無限井戸モデルの計算式はキャリア閉じ込めの仕組みを理解する上で有用だが、実際のデバイスを扱う際には材料層の構造、バンドオフセット、有効質量、電界の影響も計算に含める必要がある。OghmaNanoでは、考慮すべき物理的効果に応じて複数の量子井戸モデリングレベルが用意されている：

• A 放物線型バンドソルバー 単純な電子および正孔の閉じ込め計算を行うために用いられる。
• 多バンドk·pモデルとは、異なる電子バンド間の結合や混合を記述するためのモデルであり、8バンド亜鉛鉱型構造、10バンド亜鉛鉱型構造、8バンドウルツ鉱型構造といったバリエーションが存在する。
• 適切な半導体系向けに、専門的なシリコンバレー処理法およびルッティンガー・コーン処理法が存在する。

得られた固有値および波動関数は、その後光学的遷移や利得、キャリア散乱の計算に利用可能である。GaAs/AlGaAs抛物線型量子井戸チュートリアルが実践的な入門資料となっている。さらに多バンドk·pモデル、シリコン/SiGe量子井戸、光スペクトル、キャリア移動度などを扱う高度な計算用にはk·pバンド構造・量子井戸・移動度計算ソルバーが用意されている。

なぜ量子井戸が注目されるのか？

量子井戸が持つ最も有用な特性の一つは、井戸の幅を変えるだけで電子および正孔のエネルギー準位を制御できる点である。井戸幅を狭くすると閉じ込められたエネルギー準位間の隔たりが大きくなり、逆に幅を広くするとその隔たりが小さくなる。このことにより、材料の化学組成を変更することなくその電子的・光学的特性を制御することが可能となる。

量子井戸は半導体デバイスにおいて極めて重要であり、これにより電荷キャリアの存在位置や利用可能なエネルギー遷移を制御することが可能となる。その応用例としては以下が挙げられる：

• 半導体レーザーおよびLEDにおいては、量子井戸が電子および正孔を薄い活性層内に閉じ込めることで、井戸の幅や材料組成を調整することにより発光波長を制御することが可能となる。
• 光検出器：量子井戸内におけるエネルギー準位間の隔たりを制御することで、赤外線を含む特定の波長の光を吸収させることが可能となる。
• 高電子移動度トランジスタ（HEMT）：電子は薄い導電性チャネル内に閉じ込められると同時に、界面に平行な方向へ自由に移動できるため、二次元電子ガスが形成される。
• 量子カスケードレーザー：複数の量子井戸を結合させることで、電子が移動するための人工的に設計されたエネルギー準位を形成する。電子がこれらの準位間を遷移する際に光子が放出されるのである。

重要な点は、量子井戸構造の幾何形状を調整することで電子のエネルギー準位を自在に設計できるということです。こうしたエネルギー準位がどのように形成されるかを理解するためには、井戸内に閉じ込められた粒子に対してシュレーディンガー方程式を解く必要があります。