tecktom.com

この章のまとめ

2026年10月15日13分で読めますAI生成教科書『虎の巻』第4部 第1章 自由電子気体と Sommerfeld 展開

AI 生成・木下が未勉強

ハミルトニアンの選び方で固体を考える 第4部 第1章 §4.1.7 ★★★ 自由電子気体は密度だけで決まる Fermi 球を持ち、低温の熱と応答は Sommerfeld 展開で D(εF)D(\varepsilon _{\mathrm F}) に帰着する。比熱 γT\gamma T・Pauli 常磁性・圧縮率が並び、Wilson 比は 11。γ\gamma の実測は貴金属で自由電子の値の 1.01.0〜1.41.4 倍

この章では、金属の伝導電子からイオンの周期ポテンシャルと電子間の相互作用を削った自由電子気体 (4.1.2) を選んだ。絶対零度の基底状態は半径 kF=(3π2n)1/3k_{\mathrm F}=(3\pi^2n)^{1/3} の Fermi 球で、Fermi エネルギーは銅で 7.03 eV7.03\,\mathrm{eV}、温度に直して 8.16×104 K8.16\times10^4\,\mathrm K である。排他律のために電子は絶対零度でも1個あたり 35εF\frac35\varepsilon _{\mathrm F} の運動エネルギーと、数十 GPa の縮退圧を持つ。

有限温度の量はどれも ∫H(ε)f(ε) dε\int H(\varepsilon )f(\varepsilon )\,\mathrm{d}\varepsilon の形をとり、−∂f/∂ε-\partial f/\partial\varepsilon が幅 kBTk_{\mathrm B} T の対称な山であることから、Sommerfeld 展開 (4.1.9) が得られる。最初の補正は π26(kBT)2H′(μ)\frac{\pi^2}{6}( k_{\mathrm B} T)^2H'(\mu) で、級数は漸近級数である。これを当てると、化学ポテンシャルは D′/DD'/D の向きに (kBT)2( k_{\mathrm B} T)^2 だけ動き (4.1.10)、内部エネルギーは π26D(εF)(kBT)2\frac{\pi^2}{6}D(\varepsilon _{\mathrm F})( k_{\mathrm B} T)^2 だけ増え (4.1.12)、比熱は cV=π23D(εF)kB2Tc_V=\frac{\pi^2}{3}D(\varepsilon _{\mathrm F}) k_{\mathrm B}^2T (4.1.13) になる。動ける電子が Fermi 準位のまわり幅 kBTk_{\mathrm B} T の中だけにいることが、比熱が古典値の T/TFT/T_{\mathrm F} 倍しかない理由である。

同じ理由で、Pauli 常磁性 (4.1.15) と圧縮率 (4.1.17) も D(εF)D(\varepsilon _{\mathrm F}) で決まり、温度にほとんどよらない。比熱と磁化率の比である Wilson 比は自由電子で 11 になり (4.1.18)、11 からのずれが相互作用の目印になる。

削ったものは、周期ポテンシャル・電子間の相互作用・格子振動と電子–格子の結合・有限の大きさの4つで、札は4枚置いた。4.1.3節以降の結果は D(εF)D(\varepsilon _{\mathrm F}) だけで書いたので、帯の構造を戻しても実際の状態密度を入れれば成り立つ。実験では、低温の比熱 C=γT+βLT3C=\gamma T+\beta_{\mathrm L}T^3 から γ\gamma を読み、銅・銀・金で自由電子の値の 1.371.37・0.9970.997・1.081.08 倍を得た。

次の第4部第2章では、この自由電子気体を電場の中に置き、Boltzmann 方程式と緩和時間近似で電気伝導を求める。そこでも Sommerfeld 展開の考え方(Fermi 準位のまわり幅 kBTk_{\mathrm B} T の電子だけが流れに関わる)が、伝導度・熱伝導度・熱電係数のすべてに現れる(第4部第3章)。

総合問題 1次元の tight-binding 鎖の電子に Sommerfeld 展開を当てる ★★★

オリジナル原子を間隔 aa で一列に並べ、隣どうしをホッピング t>0t>0 でつないだ tight-binding 鎖を考える(εa=0\varepsilon _{\mathrm a}=0)。1電子の帯は ε(k)=−2tcos⁡(ka)\varepsilon (k)=-2t\cos(ka) で、電子間の相互作用は無視する。状態密度は1原子・1スピンあたりの状態密度 Dsite(ε)D_{\mathrm{site}}(\varepsilon )(第2部第3章・第2部第6章)で書き、原子1個あたり・両スピンで数えた状態密度を 2Dsite(ε)2D_{\mathrm{site}}(\varepsilon ) とする。原子1個あたりの電子の数(詰まり方)を ν\nu(0<ν<20<\nu<2)とする。

  1. (1) 2Dsite(ε)=2π4t2−ε22D_{\mathrm{site}}(\varepsilon )=\dfrac{2}{\pi\sqrt{4t^2-\varepsilon ^2}}(∣ε∣<2t\lvert \varepsilon \rvert <2t)を示し、絶対零度の Fermi エネルギーが εF=−2tcos⁡(πν/2)\varepsilon _{\mathrm F}=-2t\cos(\pi\nu/2) になることを示せ。(使う道具:第2部第3章・第2部第6章・第2部第7章)
  2. (2) 4分の1詰め(ν=1/2\nu=1/2)で、μ(T)−εF\mu(T)-\varepsilon _{\mathrm F} を (kBT)2( k_{\mathrm B} T)^2 の次数で求めよ。kBT=0.05tk_{\mathrm B} T=0.05t での μ\mu を、数値で求めた厳密な値 −1.41121t-1.41121t と比べよ。(使う道具:4.1.2節・4.1.3節)
  3. (3) 半分詰め(ν=1\nu=1)では、μ(T)=0\mu(T)=0 が温度によらず厳密に成り立つことを示せ。(使う道具:4.1.2節)
  4. (4) ν=1/2\nu=1/2 の電子比熱 C=γatTC=\gamma_{\mathrm{at}}T(原子1個あたり)を求め、帯の底の有効質量 m∗=ℏ2/(2ta2)m^*=\hbar^2/(2ta^2) を持つ1次元の自由電子気体(同じ ν\nu)の値との比を求めよ。(使う道具:4.1.3節・第2部第3章)
  5. (5) ほとんど空の帯(ν=0.05\nu=0.05)で kBT=0.05tk_{\mathrm B} T=0.05t と 0.005t0.005t のとき、(2) と同じ式で求めた μ\mu を厳密な値 −2.031t-2.031t・−1.9925t-1.9925t と比べ、Sommerfeld 展開が使える条件を述べよ。(使う道具:4.1.2節・4.1.5節)
解答

(1) NN 原子の輪で kk は 2π/(Na)2\pi/(Na) おきに NN 個あり、kk と −k-k が同じエネルギーを持つ。[k,k+dk][k,k+\mathrm{d} k] と [−k−dk,−k][-k-\mathrm{d} k,-k] の2か所に Na2πdk\frac{Na}{2\pi}\mathrm{d} k 個ずつ、スピン2つでdN=4Na2πdk\mathrm{d} N=4\frac{Na}{2\pi}\mathrm{d} k の状態が dε=∣ε′(k)∣dk=2ta∣sin⁡ka∣dk\mathrm{d}\varepsilon =\lvert \varepsilon '(k)\rvert \mathrm{d} k=2ta\lvert \sin ka\rvert \mathrm{d} k に入る。原子1個あたりにすると

2Dsite(ε)=1NdNdε=2π12t∣sin⁡ka∣=2π4t2−ε22D_{\mathrm{site}}(\varepsilon )=\frac{1}{N}\frac{\mathrm{d} N}{\mathrm{d}\varepsilon }=\frac{2}{\pi}\frac{1}{2t\lvert \sin ka\rvert }=\frac{2}{\pi\sqrt{4t^2-\varepsilon ^2}}

(4t2sin⁡2ka=4t2−ε24t^2\sin^2ka=4t^2-\varepsilon ^2)。帯の端 ε=±2t\varepsilon =\pm2t で発散するのは1次元の van Hove 特異点である。∫−2tεF2Dsite dε=2π[arcsin⁡εF2t+π2]=ν\int_{-2t}^{\varepsilon _{\mathrm F}}2D_{\mathrm{site}}\,\mathrm{d}\varepsilon =\frac2\pi\bigl[\arcsin\frac{\varepsilon _{\mathrm F}}{2t}+\frac\pi2\bigr]=\nu から εF=2tsin⁡π(ν−1)2=−2tcos⁡πν2\varepsilon _{\mathrm F}=2t\sin\frac{\pi(\nu-1)}{2}=-2t\cos\frac{\pi\nu}{2}。kk で言えば、両スピンで詰めるので kF=πν/(2a)k_{\mathrm F}=\pi\nu/(2a) まで埋まり、εF=−2tcos⁡(kFa)\varepsilon _{\mathrm F}=-2t\cos(k_{\mathrm F}a) と同じである。

(2) (4.1.10) は DD を 2Dsite2D_{\mathrm{site}} に読み替えてそのまま使える(原子1個あたり・1スピンあたりで数えても比 Dsite′/DsiteD_{\mathrm{site}}'/D_{\mathrm{site}} は変わらない)。Dsite′/Dsite=ε/(4t2−ε2)D_{\mathrm{site}}'/D_{\mathrm{site}}=\varepsilon /(4t^2-\varepsilon ^2) なので、εF=−2t\varepsilon _{\mathrm F}=-\sqrt2t では Dsite′/Dsite=−22tD_{\mathrm{site}}'/D_{\mathrm{site}}=-\frac{\sqrt2}{2t} で

μ−εF=−π26(kBT)2Dsite′Dsite=π262(kBT)2t=1.163(kBT)2t\mu-\varepsilon _{\mathrm F}=-\frac{\pi^2}{6}( k_{\mathrm B} T)^2\frac{D_{\mathrm{site}}'}{D_{\mathrm{site}}}=\frac{\pi^2}{6\sqrt2}\frac{( k_{\mathrm B} T)^2}{t}=1.163\frac{( k_{\mathrm B} T)^2}{t}

である。3次元の自由電子と逆に μ\mu は上がる。εF\varepsilon _{\mathrm F} は帯の下半分にあり、状態密度は下(帯の底の特異点)へ向かって増えるので、下に席が多いからである。kBT=0.05tk_{\mathrm B} T=0.05t では μ=−1.41421t+0.00291t=−1.41131t\mu=-1.41421t+0.00291t=-1.41131t で、厳密な −1.41121t-1.41121t との差は 1.0×10−4t1.0\times10^{-4}t、補正 0.0029t0.0029t の 3%3\% である。

(3) 帯は ε→−ε\varepsilon \to-\varepsilon で対称(Dsite(−ε)=Dsite(ε)D_{\mathrm{site}}(-\varepsilon )=D_{\mathrm{site}}(\varepsilon ))で、ν=1\nu=1 では εF=0\varepsilon _{\mathrm F}=0 である。μ=0\mu=0 と置くと、f(−ε)=1−f(ε)f(-\varepsilon )=1-f(\varepsilon ) から∫2Dsite(ε)f(ε) dε=∫−2t02Dsite (f+1−f) dε\int2D_{\mathrm{site}}(\varepsilon )f(\varepsilon )\,\mathrm{d}\varepsilon =\int_{-2t}^{0}2D_{\mathrm{site}}\,(f+1-f)\,\mathrm{d}\varepsilon(ε>0\varepsilon >0 の部分を −ε-\varepsilon に移した)=∫−2t02Dsite dε=1=\int_{-2t}^02D_{\mathrm{site}}\,\mathrm{d}\varepsilon =1 で、どの温度でも ν=1\nu=1 を満たす。したがって μ(T)=0\mu(T)=0 は厳密であり、Sommerfeld 展開の各項も Dsite′(0)=Dsite′′′(0)=⋯=0D_{\mathrm{site}}'(0)=D_{\mathrm{site}}'''(0)=\cdots=0 によってすべて 00 になる(例題 4.1.2 と違い、この場合はべきの外の項も無い)。

(4) (4.1.13) を原子1個あたりに書き直して C=π23 2Dsite(εF)kB2TC=\frac{\pi^2}{3}\,2D_{\mathrm{site}}(\varepsilon _{\mathrm F}) k_{\mathrm B}^2T、2Dsite(−2t)=2π2 t=2πt=0.450t2D_{\mathrm{site}}(-\sqrt2t)=\frac{2}{\pi\sqrt2\,t}=\frac{\sqrt2}{\pi t}=\frac{0.450}{t}。自由電子(ε=−2t+ta2k2\varepsilon =-2t+ta^2k^2)では、同じ kF=πν/(2a)k_{\mathrm F}=\pi\nu/(2a) で 2Dsite自由=2aπ1∣dε/dk∣=2aπ12ta2kF=2π2tν2D_{\mathrm{site}}^{\text{自由}}=\frac{2a}{\pi}\frac{1}{\lvert \mathrm{d}\varepsilon /\mathrm{d} k\rvert }=\frac{2a}{\pi}\frac{1}{2ta^2k_{\mathrm F}}=\frac{2}{\pi^2t\nu}、ν=12\nu=\frac12 で 4π2t=0.405t\frac{4}{\pi^2t}=\frac{0.405}{t}。比は 2/π4/π2=π22=1.111\frac{\sqrt2/\pi}{4/\pi^2}=\frac{\pi}{2\sqrt2}=1.111 である。帯は kFk_{\mathrm F} で放物線より平らになり(群速度が小さい)、同じ詰まり方でも Fermi 準位の席が多くなる。帯の底の有効質量で見積もると γ\gamma を 10%10\% 小さく見積もることになる。数値で厳密に求めた比熱(kBT=0.02tk_{\mathrm B} T=0.02t で C/kB=0.02976C/ k_{\mathrm B}=0.02976)は、Sommerfeld の 0.029620.02962 と 0.5%0.5\% で一致する。

(5) ν=0.05\nu=0.05 では εF=−1.99383t\varepsilon _{\mathrm F}=-1.99383t で、帯の底から 0.00617t0.00617t しか離れていない。係数は −π26Dsite′/Dsite=133.2/t-\frac{\pi^2}{6}D_{\mathrm{site}}'/D_{\mathrm{site}}=133.2/t と大きくなる。kBT=0.05tk_{\mathrm B} T=0.05t では μ=−1.99383t+133.2×0.0025t=−1.661t\mu=-1.99383t+133.2\times0.0025t=-1.661t となり、厳密な −2.031t-2.031t(帯の底より下)とまったく合わない。kBT=0.005tk_{\mathrm B} T=0.005t でも −1.9905t-1.9905t で、厳密な −1.9925t-1.9925t と 0.002t0.002t ずれ、補正の大きさ 0.0033t0.0033t と同じ桁の誤差がある。Sommerfeld 展開は DsiteD_{\mathrm{site}} を εF\varepsilon _{\mathrm F} のまわり幅 kBTk_{\mathrm B} T で Taylor 展開したので、kBTk_{\mathrm B} T が εF\varepsilon _{\mathrm F} から帯の端(特異点)までの距離0.00617t0.00617t より十分小さくなければ使えない。金属で kBT≪εFk_{\mathrm B} T\ll\varepsilon _{\mathrm F} と言うときの εF\varepsilon _{\mathrm F} は、正確には「Fermi 準位から状態密度が急に変わる所までの距離」である。kBT=0.05tk_{\mathrm B} T=0.05t では帯の底の状態の占有は 0.350.35 しかなく、電子は縮退した気体と古典気体の中間にある。これは半導体のキャリアの統計(第4部第5章)の状況である。

木下の勉強ステータス

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

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

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