tecktom.com

この章のまとめ

2026年10月15日10分で読めますAI生成教科書『虎の巻』第2部 第5章 ほぼ自由な電子:ブラッグ反射でギャップが開く

AI 生成・木下が未勉強

ハミルトニアンの選び方で固体を考える 第2部 第5章 §2.5.7 ★★★ 弱い周期ポテンシャルは境界で縮退した2つの平面波だけを強く混ぜ、Bragg 反射としてギャップ 2∣VG∣2\lvert V_G\rvert を開く。強い端は tight-binding につながる

この章では、周期ポテンシャルの中の1電子の HH を、tight-binding と反対の端、自由電子の側から解いた。平面波を基底にすると、ポテンシャルの Fourier 成分 VGV_G は波数が GG だけ違う平面波どうしだけを結び、固有値問題は結晶運動量 kk ごとの中心方程式に分かれる(第1節)。VGV_G が小さいとき、エネルギーの離れた平面波どうしは ∣VG∣2\lvert V_G\rvert ^2 の大きさでしか混ざらないが、Brillouin 域の境界 ∣k∣=∣k−G∣\lvert k\rvert =\lvert k-G\rvert では2つの平面波が縮退していて、VGV_G はそれらを対等に混ぜて幅 2∣VG∣2\lvert V_G\rvert のギャップを開く。境界の条件は Bragg の条件で、境界の状態は定在波、群速度は 00、ギャップの端の有効質量は ∣VG∣\lvert V_G\rvert が小さいほど小さい(第2節)。厳密に解ける δ 関数の列では、弱い端でギャップが 2VG2V_G に一致し、強い端では帯が細くなって tight-binding の帯に移る(第3節)。2次元以上では境界ごとのギャップがエネルギーの上でずれて重なり、2価の結晶でも金属になりうる(第4節)。この HH は電子間の相互作用と内殻を弱い擬ポテンシャルにまとめ、平面波を数本で打ち切っていて、その3つに裏付けの札を置いた。ギャップの大きさは光の吸収で、バンドの形は ARPES で確かめる(第5節)。

次の第2部第6章では、ここで得たバンドの上に状態がどれだけあるか(状態密度)を数え、ギャップの端や境界で群速度が 00 になる所に van Hove 特異点が現れることを見る。どのバンドが埋まるかで金属・絶縁体が決まる話は第2部第7章で扱う。

総合問題 引力の δ 関数の列:NFE から tight-binding まで ★★★

オリジナル引力の δ 関数の列 V(x)=α∑nδ(x−na)V(x)=\alpha\sum_n\delta(x-na)(α<0\alpha<0、g=mαa/ℏ2<0g=m\alpha a/\hbar^2<0)を考える。エネルギーの単位は εX=ℏ2π2/(2ma2)\varepsilon _{\mathrm X}=\hbar^2\pi^2/(2ma^2) とする。

  1. (1) ∣g∣\lvert g\rvert が小さいとき、第1の境界 k=π/ak=\pi/a のギャップの大きさを NFE で求め、下の端と上の端の状態がそれぞれ cos⁡(πx/a)\cos(\pi x/a) 型と sin⁡(πx/a)\sin(\pi x/a) 型のどちらかを答えよ。g=−0.5g=-0.5 の厳密なギャップは 0.213εX0.213\varepsilon _{\mathrm X} である。斥力 g=0.5g=0.5 の厳密な値 0.193εX0.193\varepsilon _{\mathrm X}(例題 2.5.4)と比べよ。(使う道具:第2節・第3節)
  2. (2) 負のエネルギー ε=−ℏ2p2/(2m)\varepsilon =-\hbar^2p^2/(2m)(p>0p>0)の Bloch 状態が満たす条件を、式 (2.5.16) から求めよ。(使う道具:第3節・第2部第4章の Bloch の定理)
  3. (3) δ 関数1個 αδ(x)\alpha\delta(x) の束縛状態の波動関数とエネルギー εb\varepsilon _{\mathrm b} を求め、εb=−(g2/π2)εX\varepsilon _{\mathrm b}=-(g^2/\pi^2)\varepsilon _{\mathrm X} となることを示せ。(使う道具:1次元の束縛状態。第2部第1章で原子1個の準位に当たるもの)
  4. (4) ∣g∣≫1\lvert g\rvert \gg1 のとき、(2) の条件を pa=∣g∣pa=\lvert g\rvert のまわりで展開して、最も低い帯が ε(k)≈εa−2tcos⁡ka\varepsilon (k)\approx\varepsilon _{\mathrm a}-2t\cos ka の形になることを示し、t=2∣εb∣e−κat=2\lvert \varepsilon _{\mathrm b}\rvert \mathrm{e}^{-\kappa a}(κ=m∣α∣/ℏ2\kappa=m\lvert \alpha\rvert /\hbar^2)を導け。第2部第3章の tight-binding の帯と比べよ。(使う道具:第3節・第2部第3章)
  5. (5) g=−5g=-5 のとき εb\varepsilon _{\mathrm b}・tt・帯の幅 4t4t を εX\varepsilon _{\mathrm X} で求めよ。厳密な条件の数値解では、帯は k=0k=0 で −2.598εX-2.598\varepsilon _{\mathrm X}、k=π/ak=\pi/a で −2.461εX-2.461\varepsilon _{\mathrm X} にある。(使う道具:第3節)
解答

(1) Fourier 成分は VG=α/a=(2g/π2)εX<0V_G=\alpha/a=(2g/\pi^2)\varepsilon _{\mathrm X}<0 で、ギャップは 2∣VG∣=(4∣g∣/π2)εX=0.203εX2\lvert V_G\rvert =(4\lvert g\rvert /\pi^2)\varepsilon _{\mathrm X}=0.203\varepsilon _{\mathrm X} である。VG<0V_G<0 なので、第2節の議論から、δ 関数(引力)の上に密度の山を持つ cos⁡(πx/a)\cos(\pi x/a) 型が下の端、節を持つ sin⁡(πx/a)\sin(\pi x/a) 型が上の端である。sin⁡\sin 型は δ 関数を感じないので、上の端は厳密に εX\varepsilon _{\mathrm X} に留まる(例題 2.5.4 (3) と同じ理由)。厳密なギャップは引力で 0.2130.213、斥力で 0.1930.193 で、NFE の 0.2030.203 をはさんでほぼ対称に ±0.01\pm0.01 ずれる。gg の1次では符号によらず同じで、ずれは gg の2次の補正である。

(2) 式 (2.5.16) は qq について解析的なので、q=ipq=\mathrm{i} p を代入してよい。cos⁡(ipa)=cosh⁡pa\cos(\mathrm{i} pa)=\cosh pa、sin⁡(ipa)/(ipa)=sinh⁡(pa)/(pa)\sin(\mathrm{i} pa)/(\mathrm{i} pa)=\sinh(pa)/(pa) から

cos⁡ka=cosh⁡pa+g sinh⁡papa\cos ka=\cosh pa+g\,\frac{\sinh pa}{pa}

である。g<0g<0 なら右辺の第2項が負になり、右辺が [−1,1][-1,1] に入る pp が存在しうる。

(3) x≠0x\neq0 では ψ′′=p2ψ\psi''=p^2\psi で、遠くで 00 になる解は ψ=e−p∣x∣\psi=\mathrm{e}^{-p\lvert x\rvert } である。導関数の跳び ψ′(0+)−ψ′(0−)=−2p=2mαℏ2ψ(0)=2mαℏ2\psi'(0^+)-\psi'(0^-)=-2p=\frac{2m\alpha}{\hbar^2}\psi(0)=\frac{2m\alpha}{\hbar^2} から p=−mα/ℏ2=m∣α∣/ℏ2=κp=-m\alpha/\hbar^2=m\lvert \alpha\rvert /\hbar^2=\kappa で、

εb=−ℏ2κ22m=−ℏ22ma2(κa)2=−ℏ2g22ma2=−g2π2εX\varepsilon _{\mathrm b}=-\frac{\hbar^2\kappa^2}{2m}=-\frac{\hbar^2}{2ma^2}(\kappa a)^2=-\frac{\hbar^2g^2}{2ma^2}=-\frac{g^2}{\pi^2}\varepsilon _{\mathrm X}

である(κa=∣g∣\kappa a=\lvert g\rvert)。束縛状態の広がりは 1/κ=a/∣g∣1/\kappa=a/\lvert g\rvert で、∣g∣≫1\lvert g\rvert \gg1 なら隣の δ 関数までの距離 aa より十分短い。

(4) x≔pax\coloneqq pa、f(x)≔cosh⁡x−∣g∣ sinh⁡xxf(x)\coloneqq\cosh x-\lvert g\rvert \,\frac{\sinh x}{x} とおくと、条件は f(x)=cos⁡kaf(x)=\cos ka である。x=∣g∣x=\lvert g\rvert では f=cosh⁡∣g∣−sinh⁡∣g∣=e−∣g∣f=\cosh\lvert g\rvert -\sinh\lvert g\rvert =\mathrm{e}^{-\lvert g\rvert } で、ほぼ 00 である。f′(x)=sinh⁡x−∣g∣cosh⁡xx+∣g∣sinh⁡xx2f'(x)=\sinh x-\lvert g\rvert \frac{\cosh x}{x}+\lvert g\rvert \frac{\sinh x}{x^2} は x=∣g∣x=\lvert g\rvert で sinh⁡∣g∣−cosh⁡∣g∣+sinh⁡∣g∣∣g∣≈e∣g∣2∣g∣\sinh\lvert g\rvert -\cosh\lvert g\rvert +\frac{\sinh\lvert g\rvert }{\lvert g\rvert }\approx\frac{\mathrm{e}^{\lvert g\rvert }}{2\lvert g\rvert }(∣g∣≫1\lvert g\rvert \gg1)なので、x=∣g∣+δx=\lvert g\rvert +\delta と書くと

e−∣g∣+e∣g∣2∣g∣δ=cos⁡ka⟹δ=2∣g∣ e−∣g∣(cos⁡ka−e−∣g∣)\mathrm{e}^{-\lvert g\rvert }+\frac{\mathrm{e}^{\lvert g\rvert }}{2\lvert g\rvert }\delta=\cos ka\quad\Longrightarrow\quad\delta=2\lvert g\rvert \,\mathrm{e}^{-\lvert g\rvert }\bigl(\cos ka-\mathrm{e}^{-\lvert g\rvert }\bigr)

である。エネルギーは ε=−ℏ2x22ma2≈−ℏ22ma2(g2+2∣g∣δ)=εb−2∣εb∣∣g∣δ\varepsilon =-\frac{\hbar^2x^2}{2ma^2}\approx-\frac{\hbar^2}{2ma^2}(g^2+2\lvert g\rvert \delta)=\varepsilon _{\mathrm b}-\frac{2\lvert \varepsilon _{\mathrm b}\rvert }{\lvert g\rvert }\delta なので

ε(k)≈εb+4∣εb∣e−2∣g∣⏟εa−2⋅2∣εb∣e−∣g∣⏟tcos⁡ka\varepsilon (k)\approx\underbrace{\varepsilon _{\mathrm b}+4\lvert \varepsilon _{\mathrm b}\rvert \mathrm{e}^{-2\lvert g\rvert }}_{\varepsilon _{\mathrm a}}-2\cdot\underbrace{2\lvert \varepsilon _{\mathrm b}\rvert \mathrm{e}^{-\lvert g\rvert }}_{t}\cos ka

となる。∣g∣=κa\lvert g\rvert =\kappa a なので t=2∣εb∣e−κat=2\lvert \varepsilon _{\mathrm b}\rvert \mathrm{e}^{-\kappa a} である。これは第2部第3章の tight-binding の帯 εa−2tcos⁡ka\varepsilon _{\mathrm a}-2t\cos ka と同じ形で、ホッピング tt は束縛状態の裾 e−κ∣x∣\mathrm{e}^{-\kappa\lvert x\rvert } が隣の原子まで届く分 e−κa\mathrm{e}^{-\kappa a} に比例する。原子1個の準位 εb\varepsilon _{\mathrm b} が、隣へのトンネルで幅 4t4t の帯に広がる、という tight-binding の描像が、同じ δ 関数の列の強い端として出てきた。(1) の弱い端の NFE と合わせて、1本の模型の両端が2つの描像になっている。

(5) εb=−25/π2 εX=−2.533εX\varepsilon _{\mathrm b}=-25/\pi^2\,\varepsilon _{\mathrm X}=-2.533\varepsilon _{\mathrm X}、t=2×2.533×e−5εX=0.0341εXt=2\times2.533\times\mathrm{e}^{-5}\varepsilon _{\mathrm X}=0.0341\varepsilon _{\mathrm X}、4t=0.137εX4t=0.137\varepsilon _{\mathrm X} である。厳密な帯の幅は −2.461−(−2.598)=0.137εX-2.461-(-2.598)=0.137\varepsilon _{\mathrm X}(4桁目まで見ると 0.13720.1372 と 0.13650.1365 で差は 0.5%)で一致し、帯の中心 −2.529εX-2.529\varepsilon _{\mathrm X} は εb\varepsilon _{\mathrm b} のすぐ上にある。k=0k=0 の方が低いのは t>0t>0 だからで、引力の井戸の束縛状態は同じ符号でつながる方(結合性)がエネルギーが低い、という第2部第1章の結合軌道の結果と同じである。

木下の勉強ステータス

木下はまだこの記事を勉強していません。

AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com

← AI生成教科書『虎の巻』の記事一覧へ