この記事は検証可能 な参考文献や出典 が全く示されていないか、不十分です。 出典を追加 して記事の信頼性向上にご協力ください。(このテンプレートの使い方 ) 出典検索? : "クローニッヒ・ペニーのモデル" – ニュース · 書籍 · スカラー · CiNii · J-STAGE · NDL · dlib.jp · ジャパンサーチ · TWL (2015年9月 )
クローニッヒ・ペニーのモデル (英 : Kronig–Penney model )は結晶 内での電子 の挙動を近似 的に記述する量子力学 的なモデルの1つである。周期的な井戸型ポテンシャル 型の一次元のモデルであり、狭義には周期的にデルタ関数 型のポテンシャルを持つモデルを指すこともある。1931年にラルフ・クローニッヒ とウィリアム・ペニーによって提出された。バンド理論 の基本的な枠組みをこのモデルで説明することができる。
クローニッヒ・ペニーのモデルのポテンシャル V は n を任意の整数 として以下のように表される。
V ( x ) = { 0 ( n a ≤ x < ( n + 1 ) a − b ) V 0 ( ( n + 1 ) a − b ≤ x < ( n + 1 ) a ) {\displaystyle V(x)={\begin{cases}0&(na\leq x<(n+1)a-b)\\V_{0}&((n+1)a-b\leq x<(n+1)a)\end{cases}}} このポテンシャルは周期 a を持っている。
特に重要なのは b →0 かつ V 0 →∞ の極限を取ったモデルでこれはディラックのデルタ関数を用いて以下のようなくし型関数 (comb関数)で表される。
V ( x ) = ∑ n δ ( x − n a ) {\displaystyle V(x)=\sum _{n}\delta (x-na)} これは間隔 a で一次元に配列している原子によるポテンシャルを荒く近似したものと考えることができる。
ポテンシャルが周期的な場合、ブロッホの定理 [1] よりシュレーディンガー方程式の固有関数は次を満たさなければならない。
ψ ( x ) = e i k x u ( x ) {\displaystyle \psi (x)=e^{ikx}u(x)} ここでu (x ) は、u (x + a ) = u (x ) を満たす周期関数 である。数学において k {\displaystyle k} はフロケ指数と呼ばれる。
格子の両端付近では、境界条件が問題となる。ここでボルン=フォン・カルマン境界条件 を課す。
ψ ( 0 ) = ψ ( L ) {\displaystyle \psi (0)=\psi (L)} ただし格子の長さL はL ≫ a であるとする。格子中のイオン(つまりポテンシャル井戸)の数をN とすると、aN = L である。
ブロッホの定理を適用すると、k が量子化される。
ψ ( 0 ) = e i k ⋅ 0 u ( 0 ) = e i k L u ( L ) = ψ ( L ) {\displaystyle \psi (0)=e^{ik\cdot 0}u(0)=e^{ikL}u(L)=\psi (L)} ⇒ u ( 0 ) = e i k L u ( L ) = e i k L u ( N a ) {\displaystyle \Rightarrow u(0)=e^{ikL}u(L)=e^{ikL}u(Na)} ⇒ e i k L = 1 {\displaystyle \Rightarrow e^{ikL}=1} ⇒ k L = 2 π n {\displaystyle \Rightarrow kL=2\pi n} ⇒ k = 2 π L n ( n = 0 , ± 1 , … , ± N 2 ) {\displaystyle \Rightarrow k={2\pi \over L}n\qquad \left(n=0,\pm 1,\dots ,\pm {N \over 2}\right)} ブロッホの定理を用いると、1周期での解だけを見つければ良いことになる。 ポテンシャルの1周期の中には2つの領域があり、それぞれを独立に解く。
本来シュレーディンガー方程式はエネルギーについての固有値方程式 であるが、ここでは一先ずエネルギー固有値Eは求めるものではないと見なす。 するとシュレーディンガー方程式は微分方程式となる。 そして微分方程式の解を固有値問題に代入してEを求め、解としての妥当性を検証する。
まずE が井戸の高さより高い(E>0)として、2つの領域の解を求める。
0 < x < a − b {\displaystyle 0<x<a-b} でのシュレーディンガー方程式は、
− ℏ 2 2 m ψ x x = E ψ {\displaystyle {-\hbar ^{2} \over 2m}\psi _{xx}=E\psi } この微分方程式の解は、あるα を用いて次のように表される。
ψ = A e i α x + A ′ e − i α x = e i k x ( A e i ( α − k ) x + A ′ e − i ( α + k ) x ) {\displaystyle \psi =Ae^{i\alpha x}+A'e^{-i\alpha x}=e^{ikx}\left(Ae^{i(\alpha -k)x}+A'e^{-i(\alpha +k)x}\right)} ブロッホの定理より、
u ( x ) = A e i ( α − k ) x + A ′ e − i ( α + k ) x {\displaystyle u(x)=Ae^{i(\alpha -k)x}+A'e^{-i(\alpha +k)x}\,\!} 固有値方程式に代入することで、エネルギーはα を用いて次のように求められる。
α 2 ≡ 2 m E ℏ 2 {\displaystyle \alpha ^{2}\equiv {2mE \over \hbar ^{2}}} 同様に、 − b < x < 0 {\displaystyle -b<x<0} でのシュレーディンガー方程式は、
− ℏ 2 2 m ψ x x = ( E + V 0 ) ψ {\displaystyle {-\hbar ^{2} \over 2m}\psi _{xx}=(E+V_{0})\psi } この微分方程式の解は、あるβ を用いて次のように表される。
ψ = B e i β x + B ′ e − i β x ( β 2 ≡ 2 m ( E + V 0 ) ℏ 2 ) . {\displaystyle \psi =Be^{i\beta x}+B'e^{-i\beta x}\quad \left(\beta ^{2}\equiv {2m(E+V_{0}) \over \hbar ^{2}}\right).} ブロッホの定理より、
u ( x ) = B e i ( β − k ) x + B ′ e − i ( β + k ) x {\displaystyle u(x)=Be^{i(\beta -k)x}+B'e^{-i(\beta +k)x}} 以下、α とβ (またはE )とk が満たすべき条件について考える。
クローニッヒ・ペニーのモデルのシュレーディンガー方程式 の解の存在条件は、以下の2つの条件から導出される永年方程式 を解くことで導出される。
波動関数 ψ とその一次微分が x = 0 および x = a で連続でなくてはならない(接続条件)。 ψ ( 0 − ) = ψ ( 0 + ) ψ ′ ( 0 − ) = ψ ′ ( 0 + ) {\displaystyle \psi (0^{-})=\psi (0^{+})\qquad \psi '(0^{-})=\psi '(0^{+})} 周期的ポテンシャルに対する波動関数がブロッホの定理 を満たさなければならない。 u ( − b ) = u ( a − b ) u ′ ( − b ) = u ′ ( a − b ) {\displaystyle u(-b)=u(a-b)\qquad u'(-b)=u'(a-b)} これらの条件により、次の行列が得られる。
( 1 1 − 1 − 1 α − α − β β e i ( α − k ) ( a − b ) e − i ( α + k ) ( a − b ) − e − i ( β − k ) b − e i ( β + k ) b ( α − k ) e i ( α − k ) ( a − b ) − ( α + k ) e − i ( α + k ) ( a − b ) − ( β − k ) e − i ( β − k ) b ( β + k ) e i ( β + k ) b ) ( A A ′ B B ′ ) = ( 0 0 0 0 ) . {\displaystyle {\begin{pmatrix}1&1&-1&-1\\\alpha &-\alpha &-\beta &\beta \\e^{i(\alpha -k)(a-b)}&e^{-i(\alpha +k)(a-b)}&-e^{-i(\beta -k)b}&-e^{i(\beta +k)b}\\(\alpha -k)e^{i(\alpha -k)(a-b)}&-(\alpha +k)e^{-i(\alpha +k)(a-b)}&-(\beta -k)e^{-i(\beta -k)b}&(\beta +k)e^{i(\beta +k)b}\end{pmatrix}}{\begin{pmatrix}A\\A'\\B\\B'\end{pmatrix}}={\begin{pmatrix}0\\0\\0\\0\end{pmatrix}}.} 自明でない解を得るためには、この行列の行列式は0でなければならない。よってα とβ (つまりE )とk は次式を満たさなければならない。
cos ( k a ) = cos ( β b ) cos [ α ( a − b ) ] − α 2 + β 2 2 α β sin ( β b ) sin [ α ( a − b ) ] {\displaystyle \cos(ka)=\cos(\beta b)\cos[\alpha (a-b)]-{\alpha ^{2}+\beta ^{2} \over 2\alpha \beta }\sin(\beta b)\sin[\alpha (a-b)]} ここで簡単のため次の近似を行い、ポテンシャルをデルタ関数型にして考える。
b → 0 , V 0 → ∞ , V 0 b = c o n s t a n t {\displaystyle b\to 0,\quad V_{0}\to \infty ,\quad V_{0}b=\mathrm {constant} } ⇒ β 2 b = c o n s t a n t , α 2 b → 0 {\displaystyle \Rightarrow \beta ^{2}b=\mathrm {constant} ,\quad \alpha ^{2}b\to 0} ⇒ β b → 0 , ; sin ( β b ) → β b , cos ( β b ) → 1 {\displaystyle \Rightarrow \beta b\to 0,;\quad \sin(\beta b)\to \beta b,\quad \cos(\beta b)\to 1} すると、α (つまりE )とk は次式を満たさなければならない。
cos ( k a ) = cos ( α a ) − P sin ( α a ) α a ( P ≡ m V 0 b a ℏ 2 ) {\displaystyle \cos(ka)=\cos(\alpha a)-P{\frac {\sin(\alpha a)}{\alpha a}}\qquad \left(P\equiv {\frac {mV_{0}ba}{\hbar ^{2}}}\right)} 次にE が井戸の高さより低い場合(E >0)を考える。この場合、α とβ とk は次式を満たさなければならない。
cos ( k a ) = cos ( β b ) cosh [ α ( a − b ) ] + β 2 − α 2 2 α β sin ( β b ) sinh [ α ( a − b ) ] ( α 2 ≡ 2 m | E | ℏ 2 , β 2 ≡ 2 m ( V 0 − | E | ) ℏ 2 ) {\displaystyle \cos(ka)=\cos(\beta b)\cosh[\alpha (a-b)]+{\beta ^{2}-\alpha ^{2} \over 2\alpha \beta }\sin(\beta b)\sinh[\alpha (a-b)]\quad \left(\alpha ^{2}\equiv {2m|E| \over \hbar ^{2}},\quad \beta ^{2}\equiv {2m(V_{0}-|E|) \over \hbar ^{2}}\right)} 先ほどと同じ近似( b → 0 , V 0 → ∞ , V 0 b = c o n s t a n t {\displaystyle b\to 0,\quad V_{0}\to \infty ,\quad V_{0}b=\mathrm {constant} } )により、α とk は次式を満たさなければならない。
cos ( k a ) = cos ( α a ) + P sin ( α a ) α a {\displaystyle \cos(ka)=\cos(\alpha a)+P{\frac {\sin(\alpha a)}{\alpha a}}} P = 1.5をもつ、分散関係におけるcos(k a)に等しい式の値。黒線は k が計算できる α a {\displaystyle \alpha a} の領域。 P = 1.5をもつクローニッヒペニー模型の分散関係。 これまでの議論により、エネルギー固有値E とブロッホ関数 ψ k (x ) = u (x )exp(ikx ) で状態を指定する波数 (結晶波数 )k が満たさなければならない条件が得られた。 あるE の値を選べばα とβ が求まり、cos(ka) を計算することができる。そして両辺の arccos {\displaystyle \arccos } をとることでk を計算でき、E とk の関係(分散関係 )が得られる。
ただし電子が束縛されている場合(E < 0 )、cos(ka) が1以上または-1以下になるE が存在し、その時この方程式を満たすk は存在しない。 逆に、k = nπ /a においてエネルギー値が不連続に変化し、シュレーディンガー方程式の解が存在しない E が現れる。 このことは、ポテンシャルが周期的になったことである特別な波数(結晶波数)k ではシュレーディンガー方程式の固有関数が存在しないE が存在することを意味している。
すなわち、以下の2つの区間が存在することになる。
E について解が存在する = その E の値をとることが許容された区間。これをエネルギーバンド と呼ぶ。 E について解が存在しない = その E の値をとることが禁止された区間。これをバンドギャップ と呼ぶ。 また k に対して E が連続な一つの区間はブリュアン領域 に当たる。
クローニッヒ・ペニーモデルは、バンドギャップを示す最も単純な周期的ポテンシャルの1つである。
b →0 かつ U 0 →∞ の極限を取ったモデルにおける分散関係は、k = nπ /a (n は整数)以外の点では連続であり、k の絶対値の増加につれて E も増加する関数となる。
ポテンシャルの無い自由電子モデルにおいては波動関数は ψk (x ) = u (x )exp(ikx ) の形を持つ。一方、周期 a のポテンシャルを持つモデルにおいては、これに対応する波動関数はブロッホの定理より
ψ k ( x ) = ∑ m c m e i ( k − 2 π m a ) x {\displaystyle \psi _{k}(x)=\sum _{m}c_{m}e^{i(k-{\frac {2\pi m}{a}})x}} の形を持つ。各項の係数 cm の絶対値(その2乗が波動関数への寄与と考えられる)は m = 0 が最大である。
クローニッヒ・ペニーのデルタ関数型のポテンシャルでは係数 cm は大雑把には (k −2πm /a )2 −k 2 の絶対値が小さいほど大きくなる。もっとも大きい係数 c 0 の項と二番目に大きい絶対値を持つ項 cm の2項を用いて波動関数を
ψ k ( x ) = c 0 e i k x + c m e i ( k − 2 π m a ) x {\displaystyle \psi _{k}(x)=c_{0}e^{ikx}+c_{m}e^{i(k-{\frac {2\pi m}{a}})x}} と近似できる。
k >0, U 0 > 0 の条件を前提とすると、0 < k < π /a においては、c 0 と c m (m =1) は反符号であり、cm の絶対値は 0 から k が増加するにつれて増加し、π /a で c 0 と等しくなる。π /a < k < (5/3)π /a においては、c 0 と c m (m =1) は同符号であり、c m (m =1) の絶対値は 2(n +1)π /a において c 0 と等しく、k が増加するにつれて減少する。k = (5/3)π /a において (k −2πm /a )2 −k 2 の絶対値が m =1 と m =2 で等しくなり、これより k が大きくなると m =2 の項の寄与の方が大きくなる。 (5/3)π /a < k < 2π /a においては c 0 と c m (m =2) は反符号であり、c m (m =2) の絶対値は k が増加するにつれて増加し、2π /a で c 0 と等しくなる。
2π /a < k < (13/5)π /a においては c 0 と c m (m =2) は同符号であり、c m (m =2) の絶対値は 2(n +1)π /a において c 0 と等しく、k が増加するにつれて減少する。k = (13/5)π /a において (k −2πm /a )2 −k 2 の絶対値が m =2 と m =3 で等しくなり、これより k が大きくなると m =3 の項の寄与の方が大きくなる。(13/5)π /a < k < 3π /a においては c 0 と c m (m =3) は反符号であり、c m (m =2) の絶対値は k が増加するにつれて増加し、3π /a で c 0 と等しくなる。3π /a < k < (25/7)π /a においては c 0 と c m (m =3) は同符号であり、c m (m =1) の絶対値は 2(n +1)π /a において c 0 と等しく、k が増加するにつれて減少する。
以上のように波動関数は変化していくが、k = nπ /a においては2つの波動関数が解となっている。すなわち k を小さい側から k → nπ /a に近づけた場合の解
ψ k ( x ) = c 0 e i k x − c 0 e − i k x {\displaystyle \psi _{k}(x)=c_{0}e^{ikx}-c_{0}e^{-ikx}} と k を大きい側から k → nπ /a に近づけた場合の解
ψ k ( x ) = c 0 e i k x + c 0 e − i k x {\displaystyle \psi _{k}(x)=c_{0}e^{ikx}+c_{0}e^{-ikx}} がある。差の形式の解においては波動関数はポテンシャルが値を持つ x = na の位置で 0 となりポテンシャルの影響を受けず、自由電子モデルの場合と同じエネルギー固有値を持つ。一方、和の形式の解においては、x = na の位置で波動関数は値を持つのでポテンシャルの影響を受けた分だけ高いエネルギー固有値を持つ。これにより k = nπ /a においてエネルギーが不連続にジャンプすることになり、バンドギャップが生じることになる。
ここでデルタ型の周期ポテンシャルを考える。
V ( x ) = A ⋅ ∑ n = − ∞ ∞ δ ( x − n ⋅ a ) . {\displaystyle V(x)=A\cdot \sum _{n=-\infty }^{\infty }\delta (x-n\cdot a).} A はある定数で、a は格子定数(各サイト間の間隔)。 このポテンシャルは周期的であるため、これをフーリエ級数として展開できる。
V ( x ) = ∑ K V ~ ( K ) ⋅ e i ⋅ K ⋅ x , {\displaystyle V(x)=\sum _{K}{\tilde {V}}(K)\cdot e^{i\cdot K\cdot x},} ここで
V ~ ( K ) = 1 a ∫ − a / 2 a / 2 d x V ( x ) e − i K x = 1 a ∫ − a / 2 a / 2 d x ∑ n = − ∞ ∞ A δ ( x − n a ) e − i K x = A a {\displaystyle {\tilde {V}}(K)={\frac {1}{a}}\int _{-a/2}^{a/2}dx\,V(x)\,e^{-iKx}={\frac {1}{a}}\int _{-a/2}^{a/2}dx\sum _{n=-\infty }^{\infty }A\delta (x-na)\,e^{-iKx}={\frac {A}{a}}} . ブロッホの定理によると波動関数は ψ k ( x ) = e i k x u k ( x ) {\displaystyle \psi _{k}(x)=e^{ikx}u_{k}(x)} と表せ、 u k ( x ) {\displaystyle u_{k}(x)} は格子の周期性を持つ関数である。 このことは、 u k ( x ) {\displaystyle u_{k}(x)} もフーリエ級数として展開できることを意味する。
u k ( x ) = ∑ K u ~ k ( K ) e i K x . {\displaystyle u_{k}(x)=\sum _{K}{\tilde {u}}_{k}(K)e^{iKx}.} よって波動関数は、
ψ k ( x ) = ∑ K u ~ k ( K ) e i ( k + K ) x . {\displaystyle \psi _{k}(x)=\sum _{K}{\tilde {u}}_{k}(K)\,e^{i(k+K)x}.} これをシュレーディンガー方程式に代入すると、
[ ℏ 2 ( k + K ) 2 2 m − E k ] ⋅ u ~ k ( K ) + ∑ K ′ V ~ ( K − K ′ ) u ~ k ( K ′ ) = 0 {\displaystyle \left[{\frac {\hbar ^{2}(k+K)^{2}}{2m}}-E_{k}\right]\cdot {\tilde {u}}_{k}(K)+\sum _{K'}{\tilde {V}}(K-K')\,{\tilde {u}}_{k}(K')=0} または、
[ ℏ 2 ( k + K ) 2 2 m − E k ] ⋅ u ~ k ( K ) + A a ∑ K ′ u ~ k ( K ′ ) = 0 {\displaystyle \left[{\frac {\hbar ^{2}(k+K)^{2}}{2m}}-E_{k}\right]\cdot {\tilde {u}}_{k}(K)+{\frac {A}{a}}\sum _{K'}{\tilde {u}}_{k}(K')=0} ここで新しい関数を定義する。
f ( k ) := ∑ K ′ u ~ k ( K ′ ) {\displaystyle f(k):=\sum _{K'}{\tilde {u}}_{k}(K')} これをシュレーディンガー方程式に代入すると、
[ ℏ 2 ( k + K ) 2 2 m − E k ] ⋅ u ~ k ( K ) + A a f ( k ) = 0 {\displaystyle \left[{\frac {\hbar ^{2}(k+K)^{2}}{2m}}-E_{k}\right]\cdot {\tilde {u}}_{k}(K)+{\frac {A}{a}}f(k)=0} これを u ~ k ( K ) {\displaystyle {\tilde {u}}_{k}(K)} について解くと、
u ~ k ( K ) = 2 m ℏ 2 A a f ( k ) 2 m E k ℏ 2 − ( k + K ) 2 = 2 m ℏ 2 A a 2 m E k ℏ 2 − ( k + K ) 2 f ( k ) {\displaystyle {\tilde {u}}_{k}(K)={\frac {{\frac {2m}{\hbar ^{2}}}{\frac {A}{a}}f(k)}{{\frac {2mE_{k}}{\hbar ^{2}}}-(k+K)^{2}}}={\frac {{\frac {2m}{\hbar ^{2}}}{\frac {A}{a}}}{{\frac {2mE_{k}}{\hbar ^{2}}}-(k+K)^{2}}}\,f(k)} 全てのK についてこの式を足し合わせると、
∑ K u ~ k ( K ) = ∑ K 2 m ℏ 2 A a 2 m E k ℏ 2 − ( k + K ) 2 f ( k ) {\displaystyle \sum _{K}{\tilde {u}}_{k}(K)=\sum _{K}{\frac {{\frac {2m}{\hbar ^{2}}}{\frac {A}{a}}}{{\frac {2mE_{k}}{\hbar ^{2}}}-(k+K)^{2}}}\,f(k)} または、
f ( k ) = ∑ K 2 m ℏ 2 A a 2 m E k ℏ 2 − ( k + K ) 2 f ( k ) {\displaystyle f(k)=\sum _{K}{\frac {{\frac {2m}{\hbar ^{2}}}{\frac {A}{a}}}{{\frac {2mE_{k}}{\hbar ^{2}}}-(k+K)^{2}}}\,f(k)} 都合がよいことに、 f ( k ) {\displaystyle f(k)} は打ち消しあい、
1 = ∑ K 2 m ℏ 2 A a 2 m E k ℏ 2 − ( k + K ) 2 {\displaystyle 1=\sum _{K}{\frac {{\frac {2m}{\hbar ^{2}}}{\frac {A}{a}}}{{\frac {2mE_{k}}{\hbar ^{2}}}-(k+K)^{2}}}} または、
ℏ 2 2 m a A = ∑ K 1 2 m E k ℏ 2 − ( k + K ) 2 {\displaystyle {\frac {\hbar ^{2}}{2m}}{\frac {a}{A}}=\sum _{K}{\frac {1}{{\frac {2mE_{k}}{\hbar ^{2}}}-(k+K)^{2}}}} ここで式を簡単にするため、新しい変数を定義する。
α 2 := 2 m E k ℏ 2 {\displaystyle \alpha ^{2}:={\frac {2mE_{k}}{\hbar ^{2}}}} これを用いると、
ℏ 2 2 m a A = ∑ K 1 α 2 − ( k + K ) 2 {\displaystyle {\frac {\hbar ^{2}}{2m}}{\frac {a}{A}}=\sum _{K}{\frac {1}{\alpha ^{2}-(k+K)^{2}}}} ここでK は逆格子ベクトルである。つまりK についての和は、 2 π a {\displaystyle {\frac {2\pi }{a}}} の整数倍にわたる和である。よって、
ℏ 2 2 m a A = ∑ n = − ∞ ∞ 1 α 2 − ( k + 2 π n a ) 2 {\displaystyle {\frac {\hbar ^{2}}{2m}}{\frac {a}{A}}=\sum _{n=-\infty }^{\infty }{\frac {1}{\alpha ^{2}-(k+{\frac {2\pi n}{a}})^{2}}}} ここで部分分数分解 を用いて式を変形すると、
ℏ 2 2 m a A = ∑ n = − ∞ ∞ 1 α 2 − ( k + 2 π n a ) 2 = − 1 2 α ∑ n = − ∞ ∞ [ 1 ( k + 2 π n a ) − α − 1 ( k + 2 π n a ) + α ] = − a 4 α ∑ n = − ∞ ∞ [ 1 π n + k a 2 − α a 2 − 1 π n + k a 2 + α a 2 ] = − a 4 α [ ∑ n = − ∞ ∞ 1 π n + k a 2 − α a 2 − ∑ n = − ∞ ∞ 1 π n + k a 2 + α a 2 ] {\displaystyle {\begin{aligned}{\frac {\hbar ^{2}}{2m}}{\frac {a}{A}}&=\sum _{n=-\infty }^{\infty }{\frac {1}{\alpha ^{2}-(k+{\frac {2\pi n}{a}})^{2}}}\\&=-{\frac {1}{2\alpha }}\sum _{n=-\infty }^{\infty }\left[{\frac {1}{(k+{\frac {2\pi n}{a}})-\alpha }}-{\frac {1}{(k+{\frac {2\pi n}{a}})+\alpha }}\right]\\&=-{\frac {a}{4\alpha }}\sum _{n=-\infty }^{\infty }\left[{\frac {1}{\pi n+{\frac {ka}{2}}-{\frac {\alpha a}{2}}}}-{\frac {1}{\pi n+{\frac {ka}{2}}+{\frac {\alpha a}{2}}}}\right]\\&=-{\frac {a}{4\alpha }}\left[\sum _{n=-\infty }^{\infty }{\frac {1}{\pi n+{\frac {ka}{2}}-{\frac {\alpha a}{2}}}}-\sum _{n=-\infty }^{\infty }{\frac {1}{\pi n+{\frac {ka}{2}}+{\frac {\alpha a}{2}}}}\right]\end{aligned}}} cot関数の和の良い恒等式 (Equation 18 )
cot ( x ) = ∑ n = − ∞ ∞ 1 n π + x {\displaystyle \cot(x)=\sum _{n=-\infty }^{\infty }{\frac {1}{n\pi +x}}} を代入すると、
ℏ 2 2 m a A = − a 4 α [ cot ( k a 2 − α a 2 ) − cot ( k a 2 + α a 2 ) ] {\displaystyle {\frac {\hbar ^{2}}{2m}}{\frac {a}{A}}=-{\frac {a}{4\alpha }}\left[\cot \left({\tfrac {ka}{2}}-{\tfrac {\alpha a}{2}}\right)-\cot \left({\tfrac {ka}{2}}+{\tfrac {\alpha a}{2}}\right)\right]} cot の和を用い、sin の積(cot の和の公式の一部)
cos ( k a ) = cos ( α a ) + m A ℏ 2 α sin ( α a ) {\displaystyle \cos(ka)=\cos(\alpha a)+{\frac {mA}{\hbar ^{2}\alpha }}\sin(\alpha a)} この式は、α を通じたエネルギーと波数ベクトルk の関係を示す。見てわかるように、この式の右辺は-1 から1 の範囲のみであるため、これらの式の解が存在しないα がある。つまり系がとることができないある範囲のエネルギーがある(エネルギーギャップ)。これがいわゆるエネルギーギャップで、デルタ型または長方形型の障壁だけでなく全ての形の周期ポテンシャルで存在する。
バンド間のギャップについての別の詳細な計算について、また1次元シュレーディンガー方程式の固有値の準位分裂については文献[2] を参照。コサイン型ポテンシャル(マシュー方程式)での結果についても、文献で詳細に与えられている。
^ F. Bloch, Z. Physik 52 (1928) 555 ^ Harald J. W. Muller-Kirsten, Introduction to Quantum Mechanics: Schrodinger Equation and Path Integral, 2nd ed., World Scientific (Singapore, 2012), 325?329, 458?477.