ハミルトニアンの選び方で固体を考える 第2部 第3章 §2.3.4 ★★★ k は隣の原子との位相の差であり、2π/a だけ違う k は同じ状態を表す。帯の底では k は質量 m∗ の自由粒子の波数になり、dε/dk は波束の速さを決める。端のある鎖も N→∞ で同じ帯を与える
隣の原子との位相の差
固有状態 (2.3.7) の原子 j での振幅は eikja/N で、大きさはどの原子でも同じ 1/N である。電子は輪全体に一様に広がり、原子ごとに違うのは位相だけで、隣へ進むたびに位相が ka だけ進む。k の第一の意味は、この「原子1つあたりの位相の進み」を長さ a で割ったものである。
この見方で ε(k) の形が読める。k=0 ではすべての原子が同じ位相で、隣どうしの軌道が強め合う。これは2原子の結合軌道を全体に広げた状態で、エネルギーはいちばん低い。ka=π では隣どうしの符号が逆で、すべての結合が反結合になり、エネルギーはいちばん高い。その間の k では、隣への行列要素 −t が右隣と左隣の位相の因子 e±ika と組になり、和の −2tcoska がエネルギーになる。coska は、隣どうしの軌道がどれだけ同じ符号で重なっているかの割合と読める。帯の底は結合的、天井は反結合的で、帯の中の状態はその中間である。
k と k+2π/a は同じ状態
k を 2π/a だけずらすと、ei(k+2π/a)ja=eikjae2πij=eikja なので、どの原子での振幅も変わらない。波動関数は原子の上でしか値を持たないので、2つの k は区別できない(図 2.3.3)。k は 2π/a を法としてだけ意味を持ち、Brillouin 域の中だけを考えれば足りる。連続の空間を進む波なら波数の違う波は必ず区別できるが、とびとびの点で測る限り、点の間隔より細かく振動する波は粗い波と見分けがつかない。信号処理で標本化の折り返し(エイリアシング)と呼ぶ現象と同じものである。

図 2.3.3 ka=π/3 の波(実線)と ka=π/3+2π の波(点線)。原子の位置(黒丸)では値が完全に一致するので、原子の上にしか値を持たない tight-binding の状態としては同じものである。
帯の底の k:有効質量を持つ自由粒子の波数
帯の底 k≈0 のまわりで ε(k) を展開すると、
ε(k)≈εa−2t+ta2k2=εa−2t+2m∗ℏ2k2,m∗:=2ta2ℏ2(2.3.14)
になる。これは長さ Na の輪の上を動く質量 m∗ の自由粒子のエネルギーと同じ形で、自由粒子の波数もちょうど k=2πν/(Na) に量子化される。帯の底の状態は隣どうしの位相の差が小さく、波長が原子の間隔よりずっと長い。そのような波から見ると原子の並びはほとんど連続体で、電子は「隣へ跳ぶ」のでなく「なめらかに進む」ように振る舞う。t が大きいほど隣へ移りやすく、m∗ は軽くなる。m∗ は帯の曲がり方 d2ε/dk2 の逆数で、記法の正本の有効質量 (m∗−1)=ℏ−2∂2ε/∂k2 の1次元の場合である。自由電子の質量 m とは別の量で、t と a だけで決まる。
dε/dk と波束の速さ
k の近い状態を重ね合わせて波束を作ると、波束は群速度
v(k)=ℏ1dkdε=ℏ2tasin(ka)(2.3.15)
で進む(波の群速度は dω/dk で、量子力学では ω=ε/ℏ)。帯の底 k=0 と天井 k=π/a では v=0 で、電子は進まない。これらの k では ∣k⟩ と ∣−k⟩ が同じ準位に縮退し、実数の定在波を作れる。ka=π では、右へ進む波が原子の列で反射されて左へ進む波と同じ強さで重なっている。これは原子の列による Bragg 反射の1次元版で、ほぼ自由な電子の側から見た同じ現象を第2部第5章で扱う。速さは帯の中央 ka=π/2 で最大の 2ta/ℏ になる。
ℏk を運動量のように扱ってよいか、外から力を加えたとき ℏk がどう変わるかは、ε(k) をハミルトニアンとみなす半古典運動の話で、解析力学から様々な理論分野を組み立てる 第4部第10章で扱う。ℏk が結晶運動量と呼ばれ、2πℏ/a の分だけずれてよい形で保存することは第2部第4章で述べる。
例題 2.3.4 有効質量と最大の速さを数で見る ★★
オリジナルt=1eV、a=0.25nm の輪(または十分長い鎖)を考える。ℏ=1.0546×10−34Js、m=9.109×10−31kg、1eV=1.602×10−19J とする。
- (1) 帯の底の有効質量 m∗ を電子の質量 m との比で求めよ。
- (2) 群速度の最大値を求めよ。
- (3) 帯の幅を Wb=4t に保ったまま a を半分にすると、m∗ と最大の速さはどう変わるか。
解答
(1) m∗/m=ℏ2/(2ta2m)。ℏ2/m=1.221×10−38Jm2=7.62×10−20eVm2 なので、m∗/m=7.62×10−20/(2×1×(0.25×10−9)2)=0.610。この帯の底の電子は、真空の電子より少し軽い粒子として振る舞う。
(2) vmax=2ta/ℏ=2×1.602×10−19×0.25×10−9/1.0546×10−34=7.60×105ms−1。光速の 0.25% 程度である。
(3) m∗∝a−2 なので4倍重くなり、vmax∝a なので半分になる。帯の幅が同じでも、k の範囲 2π/a が2倍に広がるので帯は k の方向に引き延ばされ、曲がりが緩く傾きも小さくなるからである。
端のある鎖
輪を1か所で切って、原子 N 個の直線の鎖にする。行列 (2.3.2) の隅の2つの −t を 0 にしたもので、ブタジエン(N=4)やヘキサトリエン(N=6)の π電子の Hückel 模型である。並進対称性は無くなるが、固有ベクトルは
cj(m)∝sin(N+1mπj),εm=εa−2tcos(N+1mπ)(m=1,…,N)(2.3.16)
と求まる。確かめるには、sin の加法定理 sin(θ(j−1))+sin(θ(j+1))=2cosθsin(θj) を使う。これで内部の原子の式 −t(cj−1+cj+1)=(εm−εa)cj が満たされる。両端の原子 j=1 と j=N では隣が1個しかないが、式の上では「原子 0 と原子 N+1 の振幅が 0」と読めば内部と同じ式になり、sin(0)=sin(mπ)=0 がこれを満たす。
鎖の状態は、km:=mπ/((N+1)a) の右へ進む波と左へ進む波を重ねた定在波である。km は 0<k<π/a に間隔 π/((N+1)a) で N 個並ぶ。輪では −π/a<k≤π/a に間隔 2π/(Na) で N 個並んでいた。輪の k と −k が1組の定在波にまとまると考えると、正の k だけ数えて間隔が半分になったことと辻褄が合う。どちらも ε=εa−2tcos(ka) の曲線の上の点であり、k の並ぶ密度は N→∞ で一致する。したがって状態密度 (2.3.12) も、半分詰めた1原子あたりのエネルギー (2.3.13) も、N→∞ では境界条件によらない。N=1000 で比べると、∣ε−εa∣<t の準位はどちらも 334 個で、1原子あたりのエネルギーは輪で −1.27324t、鎖で −1.27251t と 10−3t 以内で一致する。
境界条件の違いが現れるのは、両端の原子の割合 2/N に比例する量だけである。固体の内部の性質(単位体積あたりのエネルギーや状態密度)を計算するときに、扱いやすい輪(周期境界条件)を使ってよい理由がこれである。逆に、端に特有の状態(表面準位・端の状態)は周期境界条件では見えない。1次元の鎖に端だけに局在した状態が現れる例は第2部第8章の SSH 模型で扱う。
例題 2.3.5 半分詰めの鎖の端のエネルギー ★★
オリジナルεa=0 とし、原子 N 個(N は偶数)の鎖 (2.3.16) に電子 N 個を詰めた基底状態を考える。電子間の相互作用は無視する。
- (1) 全エネルギーがEchain=−4t∑m=1N/2cosN+1mπ=2t−sin(π/(2(N+1)))2tになることを示せ。
- (2) Echain+4tN/π を N=10,50,250 で求め、N→∞ でどの値に近づくかを答えよ。
- (3) 同じ N の輪(N=10,50,250 はどれも 4n+2 の形)の Ering+4tN/π と比べ、輪と鎖の全エネルギーの差が何で決まるかを述べよ。
解答
(1) 下から m=1,…,N/2 の準位に電子が2個ずつ入るので、Echain=2∑m=1N/2(−2tcos(mθ))、θ:=π/(N+1) である。例題 2.3.3 の和の公式 ∑m=−MMcos(mθ)=sin((2M+1)θ/2)/sin(θ/2) から∑m=1Mcos(mθ)=21[sin((2M+1)θ/2)/sin(θ/2)−1] で、M=N/2 では (2M+1)θ/2=π/2 なので分子は 1 になる。よって Echain=−2t(1/sin(θ/2)−1)=2t−2t/sin(π/(2(N+1))) である。
(2) N=10:Echain=−12.053t で、Echain+4tN/π=0.679t。N=50:0.716t。N=250:0.725t。x:=π/(2(N+1)) と置き、1/sinx=1/x+x/6+⋯ を使うと
Echain=−π4t(N+1)+2t−6(N+1)πt+⋯,Echain+π4tN=(2−π4)t−6(N+1)πt+⋯
で、N→∞ で (2−4/π)t≈0.727t に近づく。N=10 の値が 0.048t 低いのは最後の項 πt/66 の分である。
(3) 輪の全エネルギーは例題 2.3.3 (1) の Ering=−4t/sin(π/N) で、Ering+4tN/π=−0.212t, −0.042t, −0.008t と−2πt/(3N) のように 0 へ近づく。輪の全エネルギーは N に比例する −4tN/π に 1/N の補正が付くだけである。鎖の方は、N に比例する部分は輪と同じ −4tN/π で、そこに N によらない定数 (2−4/π)t が残る。輪と鎖の差 Echain−Ering は 0.891t, 0.758t, 0.733t と (2−4/π)t に近づき、N を大きくしても消えない。この差は端の数で決まり、鎖の端は2つなので、端1つあたり (1−2/π)t≈0.363t だけエネルギーが高い。1原子あたりに直すと差は (2−4/π)t/N で、N→∞ で消える。1原子あたりの量が境界条件によらないことと、全エネルギーに端の分が残ることは両立している。
輪では結合1本あたりのエネルギーが −4t/π≈−1.273t なので、鎖が結合を1本失った分をそのまま数えると 1.273t 高くなるはずである。実際の損は 0.727t と小さい。端の近くで波動関数の形が変わり、端に近い結合が内部の結合より強くなって、失った結合の一部を取り戻すからである。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com