ハミルトニアンの選び方で固体を考える 第4部 第3章 §4.3.3 ★★★ k B T ≪ ε F k_{\mathrm B} T\ll\varepsilon _{\mathrm F} k B T ≪ ε F で3つの積分を展開すると、κ / ( σ T ) \kappa/(\sigma T) κ / ( σ T ) は物質によらない数になり、S S S は σ ( ε ) \sigma(\varepsilon ) σ ( ε ) の対数微分で決まる
3つの積分を展開する
金属では Fermi 準位の近くで σ ( ε ) \sigma(\varepsilon ) σ ( ε ) がゆっくり変わり、− ∂ f 0 / ∂ ε -\partial f_0/\partial\varepsilon − ∂ f 0 / ∂ ε の幅 k B T k_{\mathrm B} T k B T は ε F \varepsilon _{\mathrm F} ε F に比べて小さい。そこで 第4部第1章 の Sommerfeld 展開
∫ − ∞ ∞ ( − ∂ f 0 ∂ ε ) G ( ε ) d ε = G ( μ ) + π 2 6 ( k B T ) 2 G ′ ′ ( μ ) + O ( ( k B T ) 4 ) \int_{-\infty}^{\infty}\left(-\frac{\partial f_0}{\partial \varepsilon }\right)G(\varepsilon )\,\mathrm{d}\varepsilon =G(\mu)+\frac{\pi^2}{6}( k_{\mathrm B} T)^2G''(\mu)+\mathcal{O}\bigl(( k_{\mathrm B} T)^4\bigr) ∫ − ∞ ∞ ( − ∂ ε ∂ f 0 ) G ( ε ) d ε = G ( μ ) + 6 π 2 ( k B T ) 2 G ′′ ( μ ) + O ( ( k B T ) 4 )
を G ( ε ) = ( ε − μ ) n σ ( ε ) G(\varepsilon )=(\varepsilon -\mu)^n\sigma(\varepsilon ) G ( ε ) = ( ε − μ ) n σ ( ε ) に当てる。n = 0 n=0 n = 0 では G ( μ ) = σ ( μ ) G(\mu)=\sigma(\mu) G ( μ ) = σ ( μ ) 、n = 1 n=1 n = 1 では G ( μ ) = 0 G(\mu)=0 G ( μ ) = 0 ・G ′ ′ ( μ ) = 2 σ ′ ( μ ) G''(\mu)=2\sigma'(\mu) G ′′ ( μ ) = 2 σ ′ ( μ ) 、n = 2 n=2 n = 2 では G ( μ ) = 0 G(\mu)=0 G ( μ ) = 0 ・G ′ ′ ( μ ) = 2 σ ( μ ) G''(\mu)=2\sigma(\mu) G ′′ ( μ ) = 2 σ ( μ ) なので
K 0 = σ ( μ ) , K 1 = π 2 3 ( k B T ) 2 σ ′ ( μ ) , K 2 = π 2 3 ( k B T ) 2 σ ( μ ) (4.3.14) \mathcal K_0=\sigma(\mu),\qquad \mathcal K_1=\frac{\pi^2}{3}( k_{\mathrm B} T)^2\sigma'(\mu),\qquad \mathcal K_2=\frac{\pi^2}{3}( k_{\mathrm B} T)^2\sigma(\mu)\tag{4.3.14} K 0 = σ ( μ ) , K 1 = 3 π 2 ( k B T ) 2 σ ′ ( μ ) , K 2 = 3 π 2 ( k B T ) 2 σ ( μ ) ( 4.3.14 )
を得る(それぞれ O ( ( k B T / ε F ) 2 ) \mathcal{O}(( k_{\mathrm B} T/\varepsilon _{\mathrm F})^2) O (( k B T / ε F ) 2 ) の相対的な補正を除く)。
Wiedemann–Franz 則
(4.3.14) を (4.3.11) へ入れる。K 1 2 / K 0 \mathcal K_1^2/\mathcal K_0 K 1 2 / K 0 は ( k B T ) 4 ( k_{\mathrm B} T)^4 ( k B T ) 4 の大きさなので主要項では落ち、
κ = π 2 3 k B 2 T e 2 σ , κ σ T = L 0 ≔ π 2 3 ( k B e ) 2 = 2.443 × 10 − 8 W Ω K − 2 (4.3.15) \kappa=\frac{\pi^2}{3}\frac{ k_{\mathrm B}^2T}{e^2}\sigma,\qquad \frac{\kappa}{\sigma T}=L_0\coloneqq\frac{\pi^2}{3}\left(\frac{ k_{\mathrm B}}{e}\right)^2=2.443\times10^{-8}\,\mathrm{W}\,\Omega\,\mathrm{K^{-2}}\tag{4.3.15} κ = 3 π 2 e 2 k B 2 T σ , σ T κ = L 0 : = 3 π 2 ( e k B ) 2 = 2.443 × 1 0 − 8 W Ω K − 2 ( 4.3.15 )
となる。これが Wiedemann–Franz 則 で、L 0 L_0 L 0 を Lorenz 数 と呼ぶ。金属の熱伝導率と電気伝導度の比が金属の種類によらずほぼ同じであることはG. Wiedemann, R. Franz (1853), Ann. Phys. 165 , 497 が測定から見つけ、L 0 L_0 L 0 の値 π 2 / 3 \pi^2/3 π 2 /3 は A. Sommerfeld (1928), Z. Phys. 47 , 1 が Fermi 統計から導いた。
(4.3.15) では、σ ( ε ) \sigma(\varepsilon ) σ ( ε ) の形も τ \tau τ の値も消えている。理由は例題 4.3.1 で見たとおりで、電荷を運ぶ電子と熱を運ぶ電子は同じ電子であり、Fermi 準位の上下 k B T k_{\mathrm B} T k B T の幅の電子が同じ緩和時間で散乱される。電荷 e e e の代わりにエネルギー ε − μ \varepsilon -\mu ε − μ を運ばせたときの比は、その幅の電子のエネルギーの揺らぎ ⟨ ( ε − μ ) 2 ⟩ = π 2 3 ( k B T ) 2 \langle(\varepsilon -\mu)^2\rangle=\frac{\pi^2}{3}( k_{\mathrm B} T)^2 ⟨( ε − μ ) 2 ⟩ = 3 π 2 ( k B T ) 2 だけで決まる。この理由から、Wiedemann–Franz 則の前提は「散乱が弾性的で、電荷の流れと熱の流れが同じ時間で緩和すること」である。4.3.5 節で、この前提が破れる場合を扱う。
Mott 公式
(4.3.14) を (4.3.10) へ入れると
S = − π 2 3 k B 2 T e d ln σ ( ε ) d ε ∣ ε = μ (4.3.16) S=-\frac{\pi^2}{3}\frac{ k_{\mathrm B}^2T}{e}\left.\frac{\mathrm{d}\ln\sigma(\varepsilon )}{\mathrm{d}\varepsilon }\right|_{\varepsilon =\mu}\tag{4.3.16} S = − 3 π 2 e k B 2 T d ε d ln σ ( ε ) ε = μ ( 4.3.16 )
を得る。これを Mott 公式 と呼ぶ(N. F. Mott, H. Jones (1936), The Theory of the Properties of Metals and Alloys , Oxford)。Seebeck 係数は、輸送関数が Fermi 準位で上り坂か下り坂かで符号が決まる。上り坂(ε > μ \varepsilon >\mu ε > μ の電子の方が電気を多く運ぶ)なら上の電子の流れが勝って S < 0 S<0 S < 0 、下り坂なら S > 0 S>0 S > 0 である。後者はバンドの上端の近く、正孔の描像が自然な所で起きる。
自由電子で、緩和時間がエネルギーのべき τ ( ε ) ∝ ε r \tau(\varepsilon )\propto\varepsilon ^r τ ( ε ) ∝ ε r で変わるとする。n ( ε ) ∝ ε 3 / 2 n(\varepsilon )\propto\varepsilon ^{3/2} n ( ε ) ∝ ε 3/2 なので σ ( ε ) ∝ ε 3 / 2 + r \sigma(\varepsilon )\propto\varepsilon ^{3/2+r} σ ( ε ) ∝ ε 3/2 + r で、d ln σ / d ε = ( 3 / 2 + r ) / ε F \mathrm{d}\ln\sigma/\mathrm{d}\varepsilon =(3/2+r)/\varepsilon _{\mathrm F} d ln σ / d ε = ( 3/2 + r ) / ε F となり
S = − π 2 3 k B e k B T ε F ( 3 2 + r ) = − π 2 2 k B e T T F ( 1 + 2 r 3 ) (4.3.17) S=-\frac{\pi^2}{3}\frac{ k_{\mathrm B}}{e}\frac{ k_{\mathrm B} T}{\varepsilon _{\mathrm F}}\left(\frac32+r\right)=-\frac{\pi^2}{2}\frac{ k_{\mathrm B}}{e}\frac{T}{T_{\mathrm F}}\left(1+\frac{2r}{3}\right)\tag{4.3.17} S = − 3 π 2 e k B ε F k B T ( 2 3 + r ) = − 2 π 2 e k B T F T ( 1 + 3 2 r ) ( 4.3.17 )
を得る。k B / e = 86.17 μ V K − 1 k_{\mathrm B}/e=86.17\,\mu\mathrm{V\,K^{-1}} k B / e = 86.17 μ V K − 1 が熱電能の自然な単位で、金属の S S S はその T / T F ∼ 10 − 3 T/T_{\mathrm F}\sim10^{-3} T / T F ∼ 1 0 − 3 倍、つまり数 μ V K − 1 \mu\mathrm{V\,K^{-1}} μ V K − 1 になる。電子のうち Fermi 準位の近く k B T k_{\mathrm B} T k B T の幅にいるものだけが動けるので、熱電能も比熱と同じく T / T F T/T_{\mathrm F} T / T F の小ささを持つ。
図 4.3.1 と図 4.3.2 に、自由電子で τ \tau τ が一定(r = 0 r=0 r = 0 )の場合の S S S と κ / ( σ T ) \kappa/(\sigma T) κ / ( σ T ) を、(4.3.8) の積分を数値で行った厳密な値と、Sommerfeld 展開の値とで比べた。μ \mu μ は電子の数を一定に保つように温度ごとに決めた(第4部第1章 )。T = 0.1 T F T=0.1\,T_{\mathrm F} T = 0.1 T F でも S S S の差は 1% ほど、L L L の差は 3% ほどである。高温では電子は古典的な気体に近づき、κ / ( σ T ) \kappa/(\sigma T) κ / ( σ T ) は非縮退の値 5 2 ( k B / e ) 2 = 0.760 L 0 \frac52( k_{\mathrm B}/e)^2=0.760\,L_0 2 5 ( k B / e ) 2 = 0.760 L 0 へ向かう。
図 4.3.1 自由電子・緩和時間一定の Seebeck 係数。実線は Mott 公式 (4.3.17) の − π 2 2 T / T F -\frac{\pi^2}{2}T/T_{\mathrm F} − 2 π 2 T / T F 、点線は (4.3.8) の数値積分(μ ( T ) \mu(T) μ ( T ) は電子数一定から決めた)。金属の室温は T / T F ∼ 0.004 T/T_{\mathrm F}\sim0.004 T / T F ∼ 0.004 で、図の左端のごく近くにある。
図 4.3.2 自由電子・緩和時間一定の Lorenz 比 κ / ( σ T ) \kappa/(\sigma T) κ / ( σ T ) を L 0 L_0 L 0 で割ったもの。実線は Wiedemann–Franz 則、点線は数値積分、灰色の破線は古典的な気体(非縮退)の極限 15 / ( 2 π 2 ) 15/(2\pi^2) 15/ ( 2 π 2 ) 。
例題 4.3.3 自由電子の銅の熱電能 ★★★
オリジナル Cu を自由電子の金属とみなし、T F = 8.16 × 10 4 K T_{\mathrm F}=8.16\times10^4\,\mathrm K T F = 8.16 × 1 0 4 K (第1部第1章 )を使う。T = 300 K T=300\,\mathrm K T = 300 K とする。
(1) 緩和時間が一定(r = 0 r=0 r = 0 )のとき、Seebeck 係数 S S S と Peltier 係数 Π \Pi Π を求めよ。
(2) 緩和時間が τ ∝ ε 3 / 2 \tau\propto\varepsilon ^{3/2} τ ∝ ε 3/2 で変わるとき、S S S は何倍になるか。
(3) 熱伝導がすべて電子によるとし、Wiedemann–Franz 則を使って Z T = S 2 σ T / κ ZT=S^2\sigma T/\kappa Z T = S 2 σ T / κ を (1) の場合に求めよ。
解答 (1) T / T F = 300 / ( 8.16 × 10 4 ) = 3.68 × 10 − 3 T/T_{\mathrm F}=300/(8.16\times10^4)=3.68\times10^{-3} T / T F = 300/ ( 8.16 × 1 0 4 ) = 3.68 × 1 0 − 3 。(4.3.17) で r = 0 r=0 r = 0 とするとS = − π 2 2 × 86.17 μ V K − 1 × 3.68 × 10 − 3 = − 1.56 μ V K − 1 S=-\frac{\pi^2}{2}\times86.17\,\mu\mathrm{V\,K^{-1}}\times3.68\times10^{-3}=-1.56\,\mu\mathrm{V\,K^{-1}} S = − 2 π 2 × 86.17 μ V K − 1 × 3.68 × 1 0 − 3 = − 1.56 μ V K − 1 。Π = T S = 300 K × ( − 1.56 μ V K − 1 ) = − 0.469 m V \Pi=TS=300\,\mathrm K\times(-1.56\,\mu\mathrm{V\,K^{-1}})=-0.469\,\mathrm{mV} Π = T S = 300 K × ( − 1.56 μ V K − 1 ) = − 0.469 mV 。1 A の電流は、電流と逆向きに 0.469 m W 0.469\,\mathrm{mW} 0.469 mW の熱を運ぶ(電子は電流と逆向きに動き、平均として μ \mu μ より上のエネルギーを運ぶから)。
(2) 因子 1 + 2 r / 3 1+2r/3 1 + 2 r /3 が 1 1 1 から 2 2 2 になるので、2倍の − 3.13 μ V K − 1 -3.13\,\mu\mathrm{V\,K^{-1}} − 3.13 μ V K − 1 になる。S S S の大きさは緩和時間のエネルギー依存性に直接よる。σ \sigma σ や κ \kappa κ が τ ( ε F ) \tau(\varepsilon _{\mathrm F}) τ ( ε F ) の値だけで決まるのと対照的である。
(3) Wiedemann–Franz 則 κ = L 0 σ T \kappa=L_0\sigma T κ = L 0 σ T を使うと Z T = S 2 / L 0 = ( 1.56 × 10 − 6 ) 2 / ( 2.443 × 10 − 8 ) = 1.0 × 10 − 4 ZT=S^2/L_0=(1.56\times10^{-6})^2/(2.443\times10^{-8})=1.0\times10^{-4} Z T = S 2 / L 0 = ( 1.56 × 1 0 − 6 ) 2 / ( 2.443 × 1 0 − 8 ) = 1.0 × 1 0 − 4 。S S S が k B / e k_{\mathrm B}/e k B / e の T / T F T/T_{\mathrm F} T / T F 倍しかないので、金属は熱電変換に向かない。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com