ハミルトニアンの選び方で固体を考える 第2部 第3章 §2.3.3 ★★★ N → ∞ N\to\infty N → ∞ で準位は幅 4 t 4t 4 t の区間を密に埋め、数え方は状態密度 D s i t e ( ε ) D_{\mathrm{site}}(\varepsilon ) D site ( ε ) になる。電子を半分詰めると、帯の中ほどまで埋まった金属の状態ができる
準位の間隔と帯の幅
式 (2.3.8) から帯の性質を読み取る。第一に、固有値は k = 0 k=0 k = 0 の ε a − 2 t \varepsilon _{\mathrm a}-2t ε a − 2 t から k = π / a k=\pi/a k = π / a の ε a + 2 t \varepsilon _{\mathrm a}+2t ε a + 2 t までにあり、帯の幅は
W b ≔ ε ( π / a ) − ε ( 0 ) = 4 t (2.3.10) W_{\mathrm b}\coloneqq\varepsilon (\pi/a)-\varepsilon (0)=4t\tag{2.3.10} W b : = ε ( π / a ) − ε ( 0 ) = 4 t ( 2.3.10 )
で N N N によらない(W b W_{\mathrm b} W b は記法の正本に無い記号で、仕事 W W W と分けるため添字 b を付けた)。帯の幅は原子の数でなく、隣の原子とのつながりの強さ t t t と隣の数(1次元では2)で決まる。
第二に、k k k は間隔 2 π / ( N a ) 2\pi/(Na) 2 π / ( N a ) で並ぶので、隣り合う k k k の準位の差は ∣ d ε / d k ∣ ⋅ 2 π / ( N a ) = 2 t ∣ sin k a ∣ ⋅ 2 π / N \lvert \mathrm{d}\varepsilon /\mathrm{d} k\rvert \cdot2\pi/(Na)=2t\lvert \sin ka\rvert \cdot2\pi/N ∣ d ε / d k ∣ ⋅ 2 π / ( N a ) = 2 t ∣ sin k a ∣ ⋅ 2 π / N である。± k \pm k ± k の縮退を除いた異なる準位の間隔もこれと同じ大きさで、帯の中央(k a = π / 2 ka=\pi/2 k a = π /2 )で最大の 4 π t / N 4\pi t/N 4 π t / N 、帯の底と天井では sin k a → 0 \sin ka\to0 sin k a → 0 なので小さくなる。帯の底の付近では ε ( k ) − ε ( 0 ) ≈ t a 2 k 2 \varepsilon (k)-\varepsilon (0)\approx ta^2k^2 ε ( k ) − ε ( 0 ) ≈ t a 2 k 2 なので、底から1番目の準位までの間隔は t ( 2 π / N ) 2 = 4 π 2 t / N 2 t(2\pi/N)^2=4\pi^2t/N^2 t ( 2 π / N ) 2 = 4 π 2 t / N 2 で、N − 2 N^{-2} N − 2 で狭くなる。どちらにしても N → ∞ N\to\infty N → ∞ で間隔は 0 0 0 になり、準位は区間 [ ε a − 2 t , ε a + 2 t ] [\varepsilon _{\mathrm a}-2t,\varepsilon _{\mathrm a}+2t] [ ε a − 2 t , ε a + 2 t ] を連続に埋める。
状態密度
準位が十分に詰まると、1つ1つの準位を数える代わりに「エネルギー ε \varepsilon ε のまわりの幅 d ε \mathrm{d}\varepsilon d ε に準位がいくつあるか」で数える方が便利になる。この章では状態密度 を、原子1個あたり・スピンの向き1つあたりで
D s i t e ( ε ) ≔ 1 N ∑ ν = 0 N − 1 δ ( ε − ε ( k ν ) ) (2.3.11) D_{\mathrm{site}}(\varepsilon )\coloneqq\frac1N\sum_{\nu=0}^{N-1}\delta\bigl(\varepsilon -\varepsilon (k_\nu)\bigr)\tag{2.3.11} D site ( ε ) : = N 1 ν = 0 ∑ N − 1 δ ( ε − ε ( k ν ) ) ( 2.3.11 )
と定める。これは1原子・1スピンあたりの状態密度で、正本の D ( ε ) D(\varepsilon ) D ( ε ) (単位体積・両スピン)とは D s i t e = ( v c / 2 ) D D_{\mathrm{site}}=(v_{\mathrm c}/2)D D site = ( v c /2 ) D の関係にある(第2部第6章 )。1次元の輪では原子1個あたりで数える方が N N N を比べやすい。N → ∞ N\to\infty N → ∞ では k k k の和を積分に置き換える。k k k は間隔 2 π / ( N a ) 2\pi/(Na) 2 π / ( N a ) で並ぶので、∑ ν → N a 2 π ∫ − π / a π / a d k \sum_\nu\to\frac{Na}{2\pi}\int_{-\pi/a}^{\pi/a}\mathrm{d} k ∑ ν → 2 π N a ∫ − π / a π / a d k であり、
D s i t e ( ε ) = a 2 π ∫ − π / a π / a δ ( ε − ε ( k ) ) d k . D_{\mathrm{site}}(\varepsilon )=\frac{a}{2\pi}\int_{-\pi/a}^{\pi/a}\delta\bigl(\varepsilon -\varepsilon (k)\bigr)\,\mathrm{d} k . D site ( ε ) = 2 π a ∫ − π / a π / a δ ( ε − ε ( k ) ) d k .
ε ( k ) = ε \varepsilon (k)=\varepsilon ε ( k ) = ε を満たす k k k は ± k ε \pm k_{\varepsilon } ± k ε の2つあり、デルタ関数の公式 δ ( f ( k ) ) = ∑ δ ( k − k i ) / ∣ f ′ ( k i ) ∣ \delta(f(k))=\sum\delta(k-k_i)/\lvert f'(k_i)\rvert δ ( f ( k )) = ∑ δ ( k − k i ) / ∣ f ′ ( k i )∣ で積分すると、∣ d ε / d k ∣ = 2 t a ∣ sin k a ∣ = a 4 t 2 − ( ε − ε a ) 2 \lvert \mathrm{d}\varepsilon /\mathrm{d} k\rvert =2ta\lvert \sin ka\rvert =a\sqrt{4t^2-(\varepsilon -\varepsilon _{\mathrm a})^2} ∣ d ε / d k ∣ = 2 t a ∣ sin k a ∣ = a 4 t 2 − ( ε − ε a ) 2 を使って
D s i t e ( ε ) = a 2 π ⋅ 2 2 t a ∣ sin k ε a ∣ = 1 π 4 t 2 − ( ε − ε a ) 2 ( ∣ ε − ε a ∣ < 2 t ) (2.3.12) D_{\mathrm{site}}(\varepsilon )=\frac{a}{2\pi}\cdot\frac{2}{2ta\lvert \sin k_{\varepsilon } a\rvert }=\frac{1}{\pi\sqrt{4t^2-(\varepsilon -\varepsilon _{\mathrm a})^2}}\qquad(\lvert \varepsilon -\varepsilon _{\mathrm a}\rvert <2t)\tag{2.3.12} D site ( ε ) = 2 π a ⋅ 2 t a ∣ sin k ε a ∣ 2 = π 4 t 2 − ( ε − ε a ) 2 1 (∣ ε − ε a ∣ < 2 t ) ( 2.3.12 )
を得る。帯の外では D s i t e ( ε ) = 0 D_{\mathrm{site}}(\varepsilon )=0 D site ( ε ) = 0 である。全体を積分すると ∫ D s i t e d ε = 1 π [ arcsin ( x / 2 t ) ] − 2 t 2 t = 1 \int D_{\mathrm{site}}\,\mathrm{d}\varepsilon =\frac1\pi[\arcsin(x/2t)]_{-2t}^{2t}=1 ∫ D site d ε = π 1 [ arcsin ( x /2 t ) ] − 2 t 2 t = 1 で、原子1個あたり1つ(スピン1つあたり)の状態があることと合う。
図 2.3.2 輪の状態密度。柱は N = 3000 N=3000 N = 3000 の輪の固有値を幅 0.125 t 0.125t 0.125 t の区間ごとに数え、N × 0.125 t N\times0.125t N × 0.125 t で割ったもの。点線は N → ∞ N\to\infty N → ∞ の式 (2.3.12) (厳密)である。帯の端で D s i t e ( ε ) D_{\mathrm{site}}(\varepsilon ) D site ( ε ) が発散するのは、端で ε ( k ) \varepsilon (k) ε ( k ) が平らになり(d ε / d k = 0 \mathrm{d}\varepsilon /\mathrm{d} k=0 d ε / d k = 0 )、狭いエネルギーに多くの k k k が集まるからである。
帯の端の発散(図 2.3.2 )は、1次元の帯に必ず現れる。ε ( k ) \varepsilon (k) ε ( k ) が極値を取る所では、k k k が変わってもエネルギーがほとんど変わらないので、そのエネルギーのまわりに準位が溜まる。発散は ( ε − ε e d g e ) − 1 / 2 (\varepsilon -\varepsilon _{\mathrm{edge}})^{-1/2} ( ε − ε edge ) − 1/2 の形で、積分すると有限なので状態の数は有限である。この種の特異点は van Hove 特異点と呼ばれ、2次元・3次元での形は第2部第6章 で扱う。
状態密度の使い方の例として、帯の中央の幅 2 t 2t 2 t 、つまり ∣ ε − ε a ∣ < t \lvert \varepsilon -\varepsilon _{\mathrm a}\rvert <t ∣ ε − ε a ∣ < t にある準位の割合を求めると∫ − t t D s i t e d ε = 2 π arcsin 1 2 = 1 3 \int_{-t}^{t}D_{\mathrm{site}}\,\mathrm{d}\varepsilon =\frac2\pi\arcsin\frac12=\frac13 ∫ − t t D site d ε = π 2 arcsin 2 1 = 3 1 である。帯の幅の半分の区間に、準位は3分の1しか無い。N = 1000 N=1000 N = 1000 の輪で実際に数えると、この区間の準位は 334 334 334 個である。
電子を詰める
原子1個が電子1個を出す場合(π電子の炭素や、価電子1個のアルカリ金属の模型)を考える。電子は Pauli の原理に従い、各 k k k の準位にスピンの上下で2個まで入る。電子は N N N 個、準位は N N N 個なので、帯のちょうど半分が埋まる。基底状態では低い方から埋めるので、− π / ( 2 a ) < k ≤ π / ( 2 a ) -\pi/(2a)<k\le\pi/(2a) − π / ( 2 a ) < k ≤ π / ( 2 a ) の準位が埋まり、最後に埋まる準位のエネルギー(Fermi 準位)は ε ( π / 2 a ) = ε a \varepsilon (\pi/2a)=\varepsilon _{\mathrm a} ε ( π /2 a ) = ε a 、帯の真ん中である。
N → ∞ N\to\infty N → ∞ での全エネルギーを原子1個あたりで求めると、
E N − ε a = 2 ⋅ a 2 π ∫ − π / 2 a π / 2 a ( − 2 t cos k a ) d k = − 4 t π ≈ − 1.273 t (2.3.13) \frac{E}{N}-\varepsilon _{\mathrm a}=2\cdot\frac{a}{2\pi}\int_{-\pi/2a}^{\pi/2a}\bigl(-2t\cos ka\bigr)\mathrm{d} k=-\frac{4t}{\pi}\approx-1.273\,t\tag{2.3.13} N E − ε a = 2 ⋅ 2 π a ∫ − π /2 a π /2 a ( − 2 t cos k a ) d k = − π 4 t ≈ − 1.273 t ( 2.3.13 )
である(先頭の 2 2 2 はスピン)。原子を2個ずつ孤立した2量体に組んだ場合は、1原子あたり − t -t − t であった(第2部第1章 )。輪全体に電子を広げると、1原子あたり 0.273 t 0.273t 0.273 t だけさらにエネルギーが下がる。電子が狭い範囲に閉じ込められず全体に広がると運動エネルギーが下がる、という第1部第1章 の6個の輪の例題の結論が、N → ∞ N\to\infty N → ∞ でも成り立っている。
埋まった準位と空の準位の境目に注目すると、いちばん上の埋まった準位(HOMO)といちばん下の空の準位(LUMO)の差は、帯の中央の準位間隔 4 π t / N 4\pi t/N 4 π t / N と同じ程度で、N → ∞ N\to\infty N → ∞ で 0 0 0 になる。電子は無限に小さいエネルギーで空の準位へ移れるので、わずかな電場で電流が流れる。半分埋まった帯は金属 を表す。バンドの埋まり方で金属・絶縁体が分かれることは第2部第7章 で一般に述べ、この半分埋まった1次元の金属が実際には格子を歪めて絶縁体になること(Peierls 転移)は第2部第8章 で扱う。
例題 2.3.3 有限の輪の全エネルギーと 4 n + 2 4n+2 4 n + 2 ★★★
オリジナル ε a = 0 \varepsilon _{\mathrm a}=0 ε a = 0 の N N N 個の輪に電子 N N N 個を詰めた基底状態(N N N は偶数)を考える。
(1) N = 4 n + 2 N=4n+2 N = 4 n + 2 (n = 1 , 2 , … n=1,2,\dots n = 1 , 2 , … )のとき、埋まる準位は ν = − n , … , n \nu=-n,\dots,n ν = − n , … , n である。全エネルギーが E = − 4 t / sin ( π / N ) E=-4t/\sin(\pi/N) E = − 4 t / sin ( π / N ) になることを示せ。
(2) N = 4 n N=4n N = 4 n のとき、ν = ± n \nu=\pm n ν = ± n の2つの準位はエネルギー 0 0 0 で縮退し、そこに電子が2個入る。全エネルギーが E = − 4 t cot ( π / N ) E=-4t\cot(\pi/N) E = − 4 t cot ( π / N ) になることを示せ。
(3) 1原子あたりのエネルギー E / N E/N E / N を N = 4 , 6 , 8 , 10 N=4,6,8,10 N = 4 , 6 , 8 , 10 で求め、N → ∞ N\to\infty N → ∞ の − 4 t / π -4t/\pi − 4 t / π にどちら側から近づくかを、1 / N 2 1/N^2 1/ N 2 の項まで展開して答えよ。
解答 (1) 埋まる準位それぞれに電子が2個入るので、E = 2 ∑ ν = − n n ( − 2 t ) cos ( 2 π ν / N ) E=2\sum_{\nu=-n}^{n}(-2t)\cos(2\pi\nu/N) E = 2 ∑ ν = − n n ( − 2 t ) cos ( 2 π ν / N ) 。和の公式 ∑ ν = − n n cos ( ν θ ) = sin ( ( 2 n + 1 ) θ / 2 ) / sin ( θ / 2 ) \sum_{\nu=-n}^{n}\cos(\nu\theta)=\sin\bigl((2n+1)\theta/2\bigr)/\sin(\theta/2) ∑ ν = − n n cos ( ν θ ) = sin ( ( 2 n + 1 ) θ /2 ) / sin ( θ /2 ) を θ = 2 π / N \theta=2\pi/N θ = 2 π / N で使うと、( 2 n + 1 ) π / N = π / 2 (2n+1)\pi/N=\pi/2 ( 2 n + 1 ) π / N = π /2 なので分子は 1 1 1 で、E = − 4 t / sin ( π / N ) E=-4t/\sin(\pi/N) E = − 4 t / sin ( π / N ) となる。
(2) エネルギー 0 0 0 の2つの準位は全エネルギーに寄与しないので、E = − 4 t ∑ ν = − ( n − 1 ) n − 1 cos ( 2 π ν / N ) = − 4 t sin ( ( 2 n − 1 ) π / N ) / sin ( π / N ) E=-4t\sum_{\nu=-(n-1)}^{n-1}\cos(2\pi\nu/N)=-4t\,\sin\bigl((2n-1)\pi/N\bigr)/\sin(\pi/N) E = − 4 t ∑ ν = − ( n − 1 ) n − 1 cos ( 2 π ν / N ) = − 4 t sin ( ( 2 n − 1 ) π / N ) / sin ( π / N ) 。( 2 n − 1 ) π / N = π / 2 − π / N (2n-1)\pi/N=\pi/2-\pi/N ( 2 n − 1 ) π / N = π /2 − π / N なので分子は cos ( π / N ) \cos(\pi/N) cos ( π / N ) で、E = − 4 t cot ( π / N ) E=-4t\cot(\pi/N) E = − 4 t cot ( π / N ) となる。2個の電子を縮退した2つの準位にどう入れるか(スピンの向き)は、電子間の相互作用を入れないこの模型では決まらない。
(3) N = 4 N=4 N = 4 :E / N = − 4 t cot ( π / 4 ) / 4 = − t E/N=-4t\cot(\pi/4)/4=-t E / N = − 4 t cot ( π /4 ) /4 = − t 。N = 6 N=6 N = 6 :− 4 t / ( 6 sin 30 ∘ ) = − 1.333 t -4t/(6\sin30^\circ)=-1.333t − 4 t / ( 6 sin 3 0 ∘ ) = − 1.333 t 。N = 8 N=8 N = 8 :− 4 t cot ( 22.5 ∘ ) / 8 = − 1.207 t -4t\cot(22.5^\circ)/8=-1.207t − 4 t cot ( 22. 5 ∘ ) /8 = − 1.207 t 。N = 10 N=10 N = 10 :− 4 t / ( 10 sin 18 ∘ ) = − 1.294 t -4t/(10\sin18^\circ)=-1.294t − 4 t / ( 10 sin 1 8 ∘ ) = − 1.294 t 。展開は N sin ( π / N ) = π − π 3 / ( 6 N 2 ) + ⋯ N\sin(\pi/N)=\pi-\pi^3/(6N^2)+\cdots N sin ( π / N ) = π − π 3 / ( 6 N 2 ) + ⋯ 、N tan ( π / N ) = π + π 3 / ( 3 N 2 ) + ⋯ N\tan(\pi/N)=\pi+\pi^3/(3N^2)+\cdots N tan ( π / N ) = π + π 3 / ( 3 N 2 ) + ⋯ から
E N ∣ 4 n + 2 = − 4 t π ( 1 + π 2 6 N 2 + ⋯ ) , E N ∣ 4 n = − 4 t π ( 1 − π 2 3 N 2 + ⋯ ) \frac EN\Big|_{4n+2}=-\frac{4t}{\pi}\Bigl(1+\frac{\pi^2}{6N^2}+\cdots\Bigr),\qquad \frac EN\Big|_{4n}=-\frac{4t}{\pi}\Bigl(1-\frac{\pi^2}{3N^2}+\cdots\Bigr) N E 4 n + 2 = − π 4 t ( 1 + 6 N 2 π 2 + ⋯ ) , N E 4 n = − π 4 t ( 1 − 3 N 2 π 2 + ⋯ )
である。4 n + 2 4n+2 4 n + 2 個の輪は − 4 t / π -4t/\pi − 4 t / π より低い側から、4 n 4n 4 n 個の輪は高い側から近づく。4 n + 2 4n+2 4 n + 2 個の輪ではいちばん上の埋まった準位までがちょうど閉じた殻になり、4 n 4n 4 n 個の輪では半端に埋まった縮退準位が残るためである。化学で「4 n + 2 4n+2 4 n + 2 個の π電子の環は安定(芳香族)」と言う Hückel 則そのものは、第2部第2章 で準位の閉殻として見た(ベンゼン N = 6 N=6 N = 6 と、正方形の形を保てないシクロブタジエン N = 4 N=4 N = 4 )。ここではそれが全エネルギーの 1 / N 2 1/N^2 1/ N 2 の差として現れる。どちらの差も 1 / N 2 1/N^2 1/ N 2 で消えるので、大きな N N N の帯には効果が残らない。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com