クローニッヒ・ペニーのモデル

クローニッヒ・ペニーのモデルは結晶内での電子の挙動を近似的に記述する量子力学的なモデルの1つである。周期的な井戸型ポテンシャル型の一次元のモデルであり、狭義には周期的にデルタ関数型のポテンシャルを持つモデルを指すこともある。1931年にラルフ・クローニッヒとウィリアム・ペニーによって提出された。バンド理論の基本的な枠組みをこのモデルで説明することができる。

1. クローニッヒ・ペニー・ポテンシャル
--------------------

クローニッヒ・ペニーのモデルのポテンシャル V}} は n}} を任意の整数として以下のように表される。

このポテンシャルは周期 //a//}} を持っている。

特に重要なのは //b//→0}} かつ //V//0→∞}} の極限を取ったモデルでこれはディラックのデルタ関数を用いて以下のようなくし型関数で表される。

これは間隔 a}} で一次元に配列している原子によるポテンシャルを荒く近似したものと考えることができる。

2. ブロッホの定理・周期的境界条件
------------------

ポテンシャルが周期的な場合、ブロッホの定理よりシュレーディンガー方程式の固有関数は次を満たさなければならない。

ここで//u//}}は、//u// =}} //u//}}を満たす周期関数である。数学においてはフロケ指数と呼ばれる。

格子の両端付近では、境界条件が問題となる。ここでボルン=フォン・カルマン境界条件を課す。

ただし格子の長さL}}は//L// ≫ //a//}}であるとする。格子中のイオンの数をN}}とすると、//aN// =}} //L//}}である。

ブロッホの定理を適用すると、k}}が量子化される。

3. シュレーディンガー方程式の解
-----------------

ブロッホの定理を用いると、1周期での解だけを見つければ良いことになる。 ポテンシャルの1周期の中には2つの領域があり、それぞれを独立に解く。

本来シュレーディンガー方程式はエネルギーについての固有値方程式であるが、ここでは一先ずエネルギー固有値Eは求めるものではないと見なす。 するとシュレーディンガー方程式は微分方程式となる。 そして微分方程式の解を固有値問題に代入してEを求め、解としての妥当性を検証する。

まず//E//が井戸の高さより高いとして、2つの領域の解を求める。

でのシュレーディンガー方程式は、

この微分方程式の解は、ある//α//}}を用いて次のように表される。

ブロッホの定理より、

固有値方程式に代入することで、エネルギーは//α//}}を用いて次のように求められる。

同様に、 でのシュレーディンガー方程式は、

この微分方程式の解は、ある//β//}}を用いて次のように表される。

ブロッホの定理より、

4. 解の存在条件
---------

以下、//α//}}と//β//}}と//k//}}が満たすべき条件について考える。

クローニッヒ・ペニーのモデルのシュレーディンガー方程式の解の存在条件は、以下の2つの条件から導出される永年方程式を解くことで導出される。

波動関数 ψ}} とその一次微分が 1=//x// = 0}} および 1=//x// = //a//}} で連続でなくてはならない。
周期的ポテンシャルに対する波動関数がブロッホの定理を満たさなければならない。これらの条件により、次の行列が得られる。

自明でない解を得るためには、この行列の行列式は0でなければならない。よって//α//}}と//β//}}と//k//}}は次式を満たさなければならない。

ここで簡単のため次の近似を行い、ポテンシャルをデルタ関数型にして考える。

すると、//α//}}と//k//}}は次式を満たさなければならない。

次に//E//が井戸の高さより低い場合を考える。この場合、//α//}}と//β//}}と//k//}}は次式を満たさなければならない。

先ほどと同じ近似により、//α//}}と//k//}}は次式を満たさなければならない。

5. バンドギャップ
----------

これまでの議論により、エネルギー固有値//E//}}とブロッホ関数 1=//ψ////k// = //u//exp}} で状態を指定する波数//k//}}が満たさなければならない条件が得られた。 ある//E//}}の値を選べば//α//}}と//β//}}が求まり、cos}}を計算することができる。そして両辺のをとることで//k//}}を計算でき、//E//}}と//k//}}の関係が得られる。

ただし電子が束縛されている場合、cos}}が1以上または-1以下になる//E//}}が存在し、その時この方程式を満たす//k//}}は存在しない。 逆に、1=//k// = //nπ/////a//}} においてエネルギー値が不連続に変化し、シュレーディンガー方程式の解が存在しない E}} が現れる。 このことは、ポテンシャルが周期的になったことである特別な波数//k//}}ではシュレーディンガー方程式の固有関数が存在しない//E//}}が存在することを意味している。

すなわち、以下の2つの区間が存在することになる。

E}} について解が存在する = その E}} の値をとることが許容された区間。これをエネルギーバンドと呼ぶ。
E}} について解が存在しない = その E}} の値をとることが禁止された区間。これをバンドギャップと呼ぶ。また k}} に対して E}} が連続な一つの区間はブリュアン領域に当たる。

クローニッヒ・ペニーモデルは、バンドギャップを示す最も単純な周期的ポテンシャルの1つである。

//b//→0}} かつ //U//0→∞}} の極限を取ったモデルにおける分散関係は、1=//k// = //nπ/////a//}} 以外の点では連続であり、k}} の絶対値の増加につれて E}} も増加する関数となる。

6. バンドギャップの生じる理由
----------------

ポテンシャルの無い自由電子モデルにおいては波動関数は 1=//ψk// = //u//exp}} の形を持つ。一方、周期 a}} のポテンシャルを持つモデルにおいては、これに対応する波動関数はブロッホの定理より

の形を持つ。各項の係数 cm}} の絶対値は 1=//m// = 0}} が最大である。

クローニッヒ・ペニーのデルタ関数型のポテンシャルでは係数 cm}} は大雑把には 2−//k//2}} の絶対値が小さいほど大きくなる。もっとも大きい係数 //c//0}} の項と二番目に大きい絶対値を持つ項 cm}} の2項を用いて波動関数を

と近似できる。

//k//>0, //U//0 > 0}} の条件を前提とすると、0 < //k// < //π/////a//}} においては、//c//0}} と 1=//c////m// }} は反符号であり、cm}} の絶対値は 0}} から k}} が増加するにつれて増加し、//π/////a//}} で //c//0}} と等しくなる。//π/////a// < //k// < //π/////a//}} においては、//c//0}} と 1=//c////m// }} は同符号であり、1=//c////m// }} の絶対値は 2//π/////a//}} において //c//0}} と等しく、k}} が増加するにつれて減少する。1=//k// = //π/////a//}} において 2−//k//2}} の絶対値が 1=//m//=1}} と 1=//m//=2}} で等しくなり、これより k}} が大きくなると 1=//m//=2}} の項の寄与の方が大きくなる。 //π/////a// < //k// < 2//π/////a//}} においては //c//0}} と 1=//c////m// }} は反符号であり、1=//c////m// }} の絶対値は k}} が増加するにつれて増加し、2//π/////a//}} で //c//0}} と等しくなる。

2//π/////a// < //k// < //π/////a//}} においては //c//0}} と 1=//c////m// }} は同符号であり、1=//c////m// }} の絶対値は 2//π/////a//}} において //c//0}} と等しく、k}} が増加するにつれて減少する。1=//k// = //π/////a//}}において 2−//k//2}} の絶対値が 1=//m//=2}} と 1=//m//=3}} で等しくなり、これより k}} が大きくなると 1=//m//=3}} の項の寄与の方が大きくなる。//π/////a// < //k// < 3//π/////a//}} においては //c//0}} と 1=//c////m// }} は反符号であり、1=//c////m// }} の絶対値は k}} が増加するにつれて増加し、3//π/////a//}} で //c//0}} と等しくなる。3//π/////a// < //k// < //π/////a//}} においては //c//0}} と 1=//c////m// }} は同符号であり、1=//c////m// }} の絶対値は 2//π/////a//}} において //c//0}} と等しく、k}} が増加するにつれて減少する。

以上のように波動関数は変化していくが、1=//k// = //nπ/////a//}} においては2つの波動関数が解となっている。すなわち k}} を小さい側から //k// → //nπ/////a//}} に近づけた場合の解

と k}} を大きい側から //k// → //nπ/////a//}} に近づけた場合の解

がある。差の形式の解においては波動関数はポテンシャルが値を持つ 1=//x// = //na//}} の位置で 0}} となりポテンシャルの影響を受けず、自由電子モデルの場合と同じエネルギー固有値を持つ。一方、和の形式の解においては、1=//x// = //na//}} の位置で波動関数は値を持つのでポテンシャルの影響を受けた分だけ高いエネルギー固有値を持つ。これにより 1=//k// = //nπ/////a//}} においてエネルギーが不連続にジャンプすることになり、バンドギャップが生じることになる。

7. クローニッヒ・ペニーモデル: 別解
--------------------

ここでデルタ型の周期ポテンシャルを考える。

A}}はある定数で、a}}は格子定数。 このポテンシャルは周期的であるため、これをフーリエ級数として展開できる。

ここで

.ブロッホの定理によると波動関数はと表せ、は格子の周期性を持つ関数である。 このことは、もフーリエ級数として展開できることを意味する。

よって波動関数は、

これをシュレーディンガー方程式に代入すると、

または、

ここで新しい関数を定義する。

これをシュレーディンガー方程式に代入すると、

これをについて解くと、

全てのK}}についてこの式を足し合わせると、

または、

都合がよいことに、は打ち消しあい、

または、

ここで式を簡単にするため、新しい変数を定義する。

これを用いると、

ここでK}}は逆格子ベクトルである。つまりK}}についての和は、の整数倍にわたる和である。よって、

ここで部分分数分解を用いて式を変形すると、

cot関数の和の良い恒等式 

を代入すると、

cot}}の和を用い、sin}}の積

この式は、α}}を通じたエネルギーと波数ベクトルk}}の関係を示す。見てわかるように、この式の右辺は-1}}から1}}の範囲のみであるため、これらの式の解が存在しないα}}がある。つまり系がとることができないある範囲のエネルギーがある。これがいわゆるエネルギーギャップで、デルタ型または長方形型の障壁だけでなく全ての形の周期ポテンシャルで存在する。

バンド間のギャップについての別の詳細な計算について、また1次元シュレーディンガー方程式の固有値の準位分裂については文献 を参照。コサイン型ポテンシャルでの結果についても、文献で詳細に与えられている。

8. 参考文献
-------