ハミルトニアンの選び方で固体を考える 第4部 第1章 §4.1.6 ★★ 統計力学の授業の Fermi 積分 f ν ( z ) f_\nu(z) f ν ( z ) と熱的 de Broglie 波長で書いた理想 Fermi 気体は、この章の D ( ε ) D(\varepsilon ) D ( ε ) と Sommerfeld 展開と同じものを別の変数で書いたもの。固体物理の授業の C = γ T + β L T 3 C=\gamma T+\beta_{\mathrm L}T^3 C = γ T + β L T 3 は4.1.5 節の式
統計力学の理想 Fermi 気体 。統計力学の授業では、理想 Fermi 気体の大分配関数を
ln Ξ = ∑ k , σ ln ( 1 + z e − β ε k ) , z = e β μ \ln\Xi=\sum_{\bm{k},\sigma}\ln\bigl(1+z\,\mathrm{e}^{-\beta\varepsilon _{\bm{k}}}\bigr),\qquad z=\mathrm{e}^{\beta\mu} ln Ξ = k , σ ∑ ln ( 1 + z e − β ε k ) , z = e β μ
と書き、k \bm{k} k の和を積分に直して
n λ T 3 g s = f 3 / 2 ( z ) , P λ T 3 g s k B T = f 5 / 2 ( z ) , f ν ( z ) = 1 Γ ( ν ) ∫ 0 ∞ x ν − 1 d x z − 1 e x + 1 , λ T = h 2 π m k B T (4.1.20) \frac{n\lambda_T^3}{g_s}=f_{3/2}(z),\qquad
\frac{P\lambda_T^3}{g_s k_{\mathrm B} T}=f_{5/2}(z),\qquad
f_\nu(z)=\frac{1}{\Gamma(\nu)}\int_0^\infty\frac{x^{\nu-1}\,\mathrm{d} x}{z^{-1}\mathrm{e}^x+1},\qquad
\lambda_T=\frac{h}{\sqrt{2\pi m k_{\mathrm B} T}}\tag{4.1.20} g s n λ T 3 = f 3/2 ( z ) , g s k B T P λ T 3 = f 5/2 ( z ) , f ν ( z ) = Γ ( ν ) 1 ∫ 0 ∞ z − 1 e x + 1 x ν − 1 d x , λ T = 2 π m k B T h ( 4.1.20 )
と習うことが多い(g s = 2 g_s=2 g s = 2 はスピンの縮重度、λ T \lambda_T λ T は熱的 de Broglie 波長、f ν ( z ) = − L i ν ( − z ) f_\nu(z)=-\mathrm{Li}_\nu(-z) f ν ( z ) = − Li ν ( − z ) )。この章の (4.1.8) の第1式に D ( ε ) D(\varepsilon ) D ( ε ) (4.1.5) を入れ、x = β ε x=\beta\varepsilon x = β ε と置いたものが第1式である。実際、n = k F 3 / ( 3 π 2 ) n=k_{\mathrm F}^3/(3\pi^2) n = k F 3 / ( 3 π 2 ) と λ T \lambda_T λ T の定義からn λ T 3 / 2 = 4 3 π ( T F / T ) 3 / 2 n\lambda_T^3/2=\frac{4}{3\sqrt\pi}(T_{\mathrm F}/T)^{3/2} n λ T 3 /2 = 3 π 4 ( T F / T ) 3/2 で、図 4.1.2 の厳密な値はこの式を z z z について解いたものである。
授業の式の読み方は「n λ T 3 n\lambda_T^3 n λ T 3 が気体の縮退の度合いを表す」である。n λ T 3 ≪ 1 n\lambda_T^3\ll1 n λ T 3 ≪ 1 なら z ≪ 1 z\ll1 z ≪ 1 で f ν ( z ) ≈ z f_\nu(z)\approx z f ν ( z ) ≈ z となり、古典的な理想気体(Boltzmann 分布)に戻る。n λ T 3 ≫ 1 n\lambda_T^3\gg1 n λ T 3 ≫ 1 なら z ≫ 1 z\gg1 z ≫ 1 で強く縮退している。この章の言葉では、n λ T 3 ∝ ( T F / T ) 3 / 2 n\lambda_T^3\propto(T_{\mathrm F}/T)^{3/2} n λ T 3 ∝ ( T F / T ) 3/2 なので、「縮退している」は「T ≪ T F T\ll T_{\mathrm F} T ≪ T F 」と同じことである。
Fermi 積分の低温展開 。授業では、z ≫ 1 z\gg1 z ≫ 1 での漸近形
f ν ( z ) = ( ln z ) ν Γ ( ν + 1 ) [ 1 + ν ( ν − 1 ) π 2 6 ( ln z ) − 2 + ν ( ν − 1 ) ( ν − 2 ) ( ν − 3 ) 7 π 4 360 ( ln z ) − 4 + ⋯ ] (4.1.21) f_\nu(z)=\frac{(\ln z)^\nu}{\Gamma(\nu+1)}\Bigl[1+\nu(\nu-1)\frac{\pi^2}{6}(\ln z)^{-2}+\nu(\nu-1)(\nu-2)(\nu-3)\frac{7\pi^4}{360}(\ln z)^{-4}+\cdots\Bigr]\tag{4.1.21} f ν ( z ) = Γ ( ν + 1 ) ( ln z ) ν [ 1 + ν ( ν − 1 ) 6 π 2 ( ln z ) − 2 + ν ( ν − 1 ) ( ν − 2 ) ( ν − 3 ) 360 7 π 4 ( ln z ) − 4 + ⋯ ] ( 4.1.21 )
を使って低温の性質を出すことが多い。これは (4.1.9) を H ( x ) = x ν − 1 / Γ ( ν ) H(x)=x^{\nu-1}/\Gamma(\nu) H ( x ) = x ν − 1 /Γ ( ν ) 、μ → ln z \mu\to\ln z μ → ln z (k B T k_{\mathrm B} T k B T を単位にした化学ポテンシャル)に当てたものそのものである。実際 G ( ln z ) = ( ln z ) ν / Γ ( ν + 1 ) G(\ln z)=(\ln z)^\nu/\Gamma(\nu+1) G ( ln z ) = ( ln z ) ν /Γ ( ν + 1 ) 、H ′ ( ln z ) = ( ν − 1 ) ( ln z ) ν − 2 / Γ ( ν ) H'(\ln z)=(\nu-1)(\ln z)^{\nu-2}/\Gamma(\nu) H ′ ( ln z ) = ( ν − 1 ) ( ln z ) ν − 2 /Γ ( ν ) で、π 2 6 H ′ / G = ν ( ν − 1 ) π 2 6 ( ln z ) − 2 \frac{\pi^2}{6}H'/G=\nu(\nu-1)\frac{\pi^2}{6}(\ln z)^{-2} 6 π 2 H ′ / G = ν ( ν − 1 ) 6 π 2 ( ln z ) − 2 になる。授業の「Fermi 積分の低温展開」と固体物理の「Sommerfeld 展開」は、同じ定理を違う変数で書いたものである。
状態密度の数え方 。固体物理の教科書には、状態密度を1スピンあたりで書くもの、試料全体(体積 V V V を掛けたもの)で書くもの、1原子あたりで書くものがある。この章は単位体積・両スピンで数えた。1スピンあたりで書いた教科書では Pauli の磁化率が 2 μ 0 μ B 2 D 1 スピン ( ε F ) 2\mu_0\mu_{\mathrm B}^2D_{1\text{スピン}}(\varepsilon _{\mathrm F}) 2 μ 0 μ B 2 D 1 スピン ( ε F ) と係数 2 2 2 を付けて現れるので、式を比べるときは数え方を先に確かめる。
比熱の式 。固体物理の授業では、比熱を
C = π 2 2 N k B T T F または C = π 2 3 D ( ε F ) V k B 2 T C=\frac{\pi^2}{2}N k_{\mathrm B}\frac{T}{T_{\mathrm F}}\quad\text{または}\quad C=\frac{\pi^2}{3}D(\varepsilon _{\mathrm F})V k_{\mathrm B}^2T C = 2 π 2 N k B T F T または C = 3 π 2 D ( ε F ) V k B 2 T
と書き、低温の測定を C / T = γ + β L T 2 C/T=\gamma+\beta_{\mathrm L}T^2 C / T = γ + β L T 2 のグラフで読むことを習う。前者は自由電子に限った形、後者は状態密度だけで書いた形で、この章の (4.1.14) と (4.1.13) に当たる。後者が帯の構造を戻しても成り立つ形であることは4.1.5 節で述べた。
例題 4.1.6 Fermi 積分の展開から化学ポテンシャルを出す ★★
オリジナル
(1) (4.1.21) の ν = 3 / 2 \nu=3/2 ν = 3/2 の2項までと (4.1.20) の第1式から、ln z = β μ \ln z=\beta\mu ln z = β μ を T / T F T/T_{\mathrm F} T / T F の2次まで求め、(4.1.11) と一致することを示せ。
(2) 室温 300 K 300\,\mathrm K 300 K の電子の λ T \lambda_T λ T と、銅の n λ T 3 / 2 n\lambda_T^3/2 n λ T 3 /2 を求めよ(h = 6.626 × 10 − 34 J s h=6.626\times10^{-34}\,\mathrm{J\,s} h = 6.626 × 1 0 − 34 J s )。
(3) ln z = 20 \ln z=20 ln z = 20 と 50 50 50 で、(1) に使った2項の展開と f 3 / 2 ( z ) f_{3/2}(z) f 3/2 ( z ) の数値を比べよ。
解答 (1) L = ln z L=\ln z L = ln z 、s = T / T F s=T/T_{\mathrm F} s = T / T F と置く。左辺は 4 3 π s − 3 / 2 \frac{4}{3\sqrt\pi}s^{-3/2} 3 π 4 s − 3/2 、右辺は Γ ( 5 / 2 ) = 3 π 4 \Gamma(5/2)=\frac{3\sqrt\pi}{4} Γ ( 5/2 ) = 4 3 π から 4 3 π L 3 / 2 [ 1 + π 2 8 L − 2 ] \frac{4}{3\sqrt\pi}L^{3/2}\bigl[1+\frac{\pi^2}{8}L^{-2}\bigr] 3 π 4 L 3/2 [ 1 + 8 π 2 L − 2 ] (ν ( ν − 1 ) = 3 4 \nu(\nu-1)=\frac34 ν ( ν − 1 ) = 4 3 )。両辺を比べて L 3 / 2 ( 1 + π 2 8 L − 2 ) = s − 3 / 2 L^{3/2}(1+\frac{\pi^2}{8}L^{-2})=s^{-3/2} L 3/2 ( 1 + 8 π 2 L − 2 ) = s − 3/2 。L = s − 1 ( 1 + c s 2 ) L=s^{-1}(1+c\,s^2) L = s − 1 ( 1 + c s 2 ) と置くと、左辺は s − 3 / 2 ( 1 + 3 2 c s 2 + π 2 8 s 2 + ⋯ ) s^{-3/2}(1+\frac32cs^2+\frac{\pi^2}{8}s^2+\cdots) s − 3/2 ( 1 + 2 3 c s 2 + 8 π 2 s 2 + ⋯ ) なので c = − π 2 12 c=-\frac{\pi^2}{12} c = − 12 π 2 。μ = k B T L = ε F s L = ε F ( 1 − π 2 12 s 2 ) \mu= k_{\mathrm B} TL=\varepsilon _{\mathrm F}\,sL=\varepsilon _{\mathrm F}(1-\frac{\pi^2}{12}s^2) μ = k B T L = ε F s L = ε F ( 1 − 12 π 2 s 2 ) で、(4.1.11) と一致する。
(2) λ T = h / 2 π m k B T = 4.31 n m \lambda_T=h/\sqrt{2\pi m k_{\mathrm B} T}=4.31\,\mathrm{nm} λ T = h / 2 π m k B T = 4.31 nm 。n λ T 3 / 2 = 8.47 × 10 28 × ( 4.31 × 10 − 9 ) 3 / 2 = 3.39 × 10 3 n\lambda_T^3/2=8.47\times10^{28}\times(4.31\times10^{-9})^3/2=3.39\times10^3 n λ T 3 /2 = 8.47 × 1 0 28 × ( 4.31 × 1 0 − 9 ) 3 /2 = 3.39 × 1 0 3 で、1 1 1 よりはるかに大きい。室温の銅の伝導電子は強く縮退した量子気体である。同じ温度・1気圧の窒素の気体(分子の質量が電子の約 5 × 10 4 5\times10^4 5 × 1 0 4 倍、密度が約 3 × 10 − 4 3\times10^{-4} 3 × 1 0 − 4 倍)では n λ T 3 n\lambda_T^3 n λ T 3 は 10 − 7 10^{-7} 1 0 − 7 程度で、古典気体として扱ってよい。
(3) 数値で比べると、ln z = 20 \ln z=20 ln z = 20 でも 50 50 50 でも2項の展開と f 3 / 2 ( z ) f_{3/2}(z) f 3/2 ( z ) の差は相対的に 10 − 5 10^{-5} 1 0 − 5 より小さい(verify_ch4-01.py)。次の項は ( ln z ) − 4 (\ln z)^{-4} ( ln z ) − 4 の大きさで、ln z = T F / T \ln z=T_{\mathrm F}/T ln z = T F / T の金属では2項目すら 10 − 6 10^{-6} 1 0 − 6 程度にとどまる。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com