tecktom.com

この章のまとめ

2026年10月15日14分で読めますAI生成教科書『虎の巻』第4部 第2章 つなぎ:Boltzmann 方程式と緩和時間近似の電気伝導

AI 生成・木下が未勉強

ハミルトニアンの選び方で固体を考える 第4部 第2章 §4.2.7 ★★★ バンドの電子の Boltzmann 方程式を弾性散乱で線形にすると、輸送緩和時間・Fermi 面の積分・正孔・交流の伝導度・和則が出る。削ったのは1バンドの半古典・独立な Born 散乱・フォノンの緩和時間・電子間の相互作用で、抵抗率の温度変化で確かめる

この章では、周期ポテンシャル・不純物・外場からなる1電子の H^\hat H (4.2.1) を選び、解析力学から様々な理論分野を組み立てる 第4部第10章の半古典運動と解析力学から様々な理論分野を組み立てる 第3部第8章の Liouville+衝突の考え方を借りて、占有率 f(r,k,t)f(\bm{r},\bm{k},t) の Boltzmann 方程式 (4.2.4) を書いた。

当てた道具は、不純物の弾性散乱の衝突の項 (4.2.6) と、電場の1次まで線形にした方程式 (4.2.7) である。Pauli の因子は弾性散乱では打ち消し合い、衝突の項は線形になる。等方なバンドでは、電場の1次のずれに対して衝突の項がちょうど −δf/τtr-\delta f/\tau_{\mathrm{tr}} になり(命題 4.2.1)、緩和時間近似が厳密な帰結になる。τtr\tau_{\mathrm{tr}} は散乱の角度に 1−cos⁡θ1-\cos\theta の重みを付けた率で決まり、前方への散乱は電流をほとんど減らさない。

出てきたものは3つある。第一に、伝導度は −∂f0/∂ε-\partial f_0/\partial\varepsilon の重みで Fermi 面の積分 (4.2.12) になり、等方なら σ=e2D(εF)D\sigma=e^2D(\varepsilon _{\mathrm F})D となる。満ちたバンドは電流を運ばず、満ちかけたバンドの伝導度は正孔の Drude 式と一致する。第二に、交流では σ(ω)=σ0/(1−iωτ)\sigma(\omega)=\sigma_0/(1-\mathrm{i}\omega\tau) で、Re⁡σ\operatorname{Re}\sigma の面積は τ\tau によらずプラズマ振動数で決まる(式 (4.2.16))。第三に、散乱の原因ごとに率を足すと抵抗率は残留抵抗とフォノンの部分の和になり、フォノンの部分は高温で温度に比例し、低温で T5T^5 に従う。T5T^5 の T2T^2 は 1−cos⁡θ1-\cos\theta の重みが温度の冪に現れたものである。

削ったものは、1つのバンドに留まる半古典の波束・独立な1回ずつの Born 散乱(干渉と弱局在を捨てた)・非弾性のフォノン散乱の緩和時間による置き換え・電子間の相互作用の4つで、前の2つは B 型の極限として条件付きで証明があり、3つめは C 型である。次の第4部第3章では、温度の勾配がある場合へ同じ線形の解を広げ、熱電効果と Wiedemann–Franz 則を出す。

総合問題 不純物のある原子の鎖の抵抗 ★★★

オリジナル第2部第3章の原子の鎖(間隔 aa、バンド ε(k)=−2tcos⁡(ka)\varepsilon (k)=-2t\cos(ka)、t>0t>0)に、原子の割合 c≪1c\ll1 で不純物がでたらめに入っているとする。不純物のある原子では軌道のエネルギーが U0U_0 だけずれる(∣U0∣≪t\lvert U_0\rvert \ll t)。Fermi 波数を kFk_{\mathrm F}(0<kFa<π0<k_{\mathrm F}a<\pi)、Fermi 速度を vF=(2ta/ℏ)sin⁡(kFa)v_{\mathrm F}=(2ta/\hbar)\sin(k_{\mathrm F}a) とし、スピンの上下を数える。

  1. (1) τ\tau を一定として、T=0T=0 の単位長さあたりの伝導度が σ=2e2τvF/(πℏ)\sigma=2e^2\tau v_{\mathrm F}/(\pi\hbar) になることを示せ。(使う道具:4.2.3節の式 (4.2.11)・第2部第3章のバンド)
  2. (2) 同じ σ\sigma を、解析力学から様々な理論分野を組み立てる 第3部第8章の形 σ=e2τn⟨1/m∗⟩\sigma=e^2\tau n\langle 1/m^* \rangle(1/m∗(k)=ℏ−2d2ε/dk21/m^*(k)=\hbar^{-2}\mathrm{d}^2\varepsilon /\mathrm{d} k^2 の、占められた状態での平均)で計算し、(1) と一致することを確かめよ。(使う道具:解析力学から様々な理論分野を組み立てる 第3部第8章・4.2.3節)
  3. (3) NN 原子の輪(長さ L=NaL=Na)で、Bloch 状態 ∣k⟩=N−1/2∑jeikja∣j⟩| k \rangle =N^{-1/2}\sum_j\mathrm{e}^{\mathrm{i} kja}| j \rangle の間の不純物の行列要素を書き、その絶対値の2乗の不純物の位置についての平均が、k′≠kk'\neq k で U02c/NU_0^2c/N になることを示せ(c≪1c\ll1)。黄金律から、1次元では弾性散乱の行き先が k′=kk'=k と k′=−kk'=-k の2つしかないことを使い、輸送緩和時間 1/τtr=2U02ca/(ℏ2vF)1/\tau_{\mathrm{tr}}=2U_0^2ca/(\hbar^2v_{\mathrm F}) を導け。(使う道具:4.2.2節の式 (4.2.6)・命題 4.2.1)
  4. (4) (1) と (3) から σ\sigma と平均自由行程 ℓ=vFτtr\ell=v_{\mathrm F}\tau_{\mathrm{tr}} を求め、σ=(4e2/h) ℓ\sigma=(4e^2/h)\,\ell と書けることを示せ。t=1 eVt=1\,\mathrm{eV}、U0=0.5 eVU_0=0.5\,\mathrm{eV}、c=0.01c=0.01、a=0.25 nma=0.25\,\mathrm{nm} の半分詰め(kFa=π/2k_{\mathrm F}a=\pi/2)で、ℓ\ell、kFℓk_{\mathrm F}\ell、長さ 1 μm1\,\mathrm{\mu m} の鎖の抵抗を求めよ。(使う道具:4.2.3節・4.2.5節)
  5. (5) 半分詰めの鎖で、τtr\tau_{\mathrm{tr}} のエネルギー依存まで入れた「エネルギー ε\varepsilon での伝導度」σ(ε)\sigma(\varepsilon ) を求め、第4部第1章の Sommerfeld 展開で温度の補正σ(T)/σ(0)=1−π212(kBT/t)2\sigma(T)/\sigma(0)=1-\frac{\pi^2}{12}\bigl( k_{\mathrm B} T/t\bigr)^2 を導け。t=1 eVt=1\,\mathrm{eV}、T=300 KT=300\,\mathrm K で補正の大きさを求めよ。(使う道具:第4部第1章・4.2.3節)
解答

(1) 例題 4.2.3 の (1) と同じで、−∂f0/∂ε=δ(ε−εF)-\partial f_0/\partial\varepsilon =\delta(\varepsilon -\varepsilon _{\mathrm F}) が2つの Fermi 点で δ(k∓kF)/(ℏvF)\delta(k\mp k_{\mathrm F})/(\hbar v_{\mathrm F}) を与えるのでσ=2e2τ⋅12π⋅2⋅vF2ℏvF=2e2τvFπℏ\sigma=2e^2\tau\cdot\frac{1}{2\pi}\cdot2\cdot\frac{v_{\mathrm F}^2}{\hbar v_{\mathrm F}}=\frac{2e^2\tau v_{\mathrm F}}{\pi\hbar}。

(2) 電子密度は n=2⋅2kF/(2π)=2kF/πn=2\cdot2k_{\mathrm F}/(2\pi)=2k_{\mathrm F}/\pi。1/m∗(k)=2ta2cos⁡(ka)/ℏ21/m^*(k)=2ta^2\cos(ka)/\hbar^2 を ∣k∣<kF\lvert k\rvert <k_{\mathrm F} で平均すると⟨1/m∗⟩=2ta2ℏ2sin⁡(kFa)kFa\langle 1/m^* \rangle =\frac{2ta^2}{\hbar^2}\frac{\sin(k_{\mathrm F}a)}{k_{\mathrm F}a}。したがってe2τn⟨1/m∗⟩=2kFπe2τ2ta2sin⁡(kFa)ℏ2kFa=4tae2τsin⁡(kFa)πℏ2e^2\tau n\langle 1/m^* \rangle =\frac{2k_{\mathrm F}}{\pi}e^2\tau\frac{2ta^2\sin(k_{\mathrm F}a)}{\hbar^2k_{\mathrm F}a}=\frac{4tae^2\tau\sin(k_{\mathrm F}a)}{\pi\hbar^2}で、vF=(2ta/ℏ)sin⁡(kFa)v_{\mathrm F}=(2ta/\hbar)\sin(k_{\mathrm F}a) を入れた (1) の 2e2τvF/(πℏ)2e^2\tau v_{\mathrm F}/(\pi\hbar) と一致する。Fermi 点の電子だけの数え方と、占められた全電子の質量の平均とは、ここでも同じ量である。

(3) 不純物の Hamiltonian は V^=U0∑j∈imp∣j⟩⟨j∣\hat V=U_0\sum_{j\in\mathrm{imp}}| j \rangle \langle j | で、

⟨k′∣V^∣k⟩=U0N∑j∈impei(k−k′)ja.\langle k' | \hat V | k \rangle =\frac{U_0}{N}\sum_{j\in\mathrm{imp}}\mathrm{e}^{\mathrm{i}(k-k')ja}.

k′≠kk'\neq k のとき、絶対値の2乗の和 ∑j,l∈impei(k−k′)(j−l)a\sum_{j,l\in\mathrm{imp}}\mathrm{e}^{\mathrm{i}(k-k')(j-l)a} を不純物の配置について平均すると、j=lj=l の項からは不純物の数 cNcN が残り、j≠lj\neq l の項は位相がでたらめなので(c≪1c\ll1 の範囲で)消える。したがって平均は U02cN/N2=U02c/NU_0^2cN/N^2=U_0^2c/N である(数値で N=400N=400、c=0.05c=0.05 の配置を 400 通り平均すると c(1−c)c(1-c) に合う。c≪1c\ll1 では cc)。黄金律で kk から k′k' への率は 2πℏU02cNδ(εk−εk′)\frac{2\pi}{\hbar}\frac{U_0^2c}{N}\delta(\varepsilon _k-\varepsilon _{k'}) で、∑k′=L2π∫dk′\sum_{k'}=\frac{L}{2\pi}\int\mathrm{d} k' と δ(εk−εk′)=∑±δ(k′∓k)/(ℏvF)\delta(\varepsilon _k-\varepsilon _{k'})=\sum_\pm\delta(k'\mp k)/(\hbar v_{\mathrm F}) を使うと、k′=−kk'=-k への後方散乱の率は

1τb=2πℏU02cN⋅Na2π⋅1ℏvF=U02caℏ2vF\frac{1}{\tau_{\mathrm b}}=\frac{2\pi}{\hbar}\frac{U_0^2c}{N}\cdot\frac{Na}{2\pi}\cdot\frac{1}{\hbar v_{\mathrm F}}=\frac{U_0^2ca}{\hbar^2v_{\mathrm F}}

である。衝突の項 (4.2.6) では、k′=kk'=k への散乱は δfk′−δfk=0\delta f_{k'}-\delta f_k=0 で寄与せず、k′=−kk'=-k への散乱は電場の1次のずれが奇関数(δf−k=−δfk\delta f_{-k}=-\delta f_k)なので −2δfk-2\delta f_k を与える。1次元の cos⁡θ\cos\theta は ±1\pm1 だけで、命題 4.2.1 の 1−cos⁡θ1-\cos\theta が後方散乱で 22 になることにあたる。したがって

1τtr=2τb=2U02caℏ2vF=U02cℏtsin⁡(kFa).\frac{1}{\tau_{\mathrm{tr}}}=\frac{2}{\tau_{\mathrm b}}=\frac{2U_0^2ca}{\hbar^2v_{\mathrm F}}=\frac{U_0^2c}{\hbar t\sin(k_{\mathrm F}a)}.

Fermi 速度が小さいバンドの端ほど状態密度が大きく、散乱の行き先が多いので、τtr\tau_{\mathrm{tr}} は短くなる。

(4) ℓ=vFτtr=ℏ2vF22U02ca=2t2asin⁡2(kFa)U02c\ell=v_{\mathrm F}\tau_{\mathrm{tr}}=\frac{\hbar^2v_{\mathrm F}^2}{2U_0^2ca}=\frac{2t^2a\sin^2(k_{\mathrm F}a)}{U_0^2c}、σ=2e2πℏℓ=4e2t2asin⁡2(kFa)πℏU02c\sigma=\frac{2e^2}{\pi\hbar}\ell=\frac{4e^2t^2a\sin^2(k_{\mathrm F}a)}{\pi\hbar U_0^2c} で、h=2πℏh=2\pi\hbar を使えば σ=(4e2/h)ℓ\sigma=(4e^2/h)\ell である。長さ LL の鎖の抵抗は R=L/σ=h4e2LℓR=L/\sigma=\frac{h}{4e^2}\frac{L}{\ell} で、量子化された抵抗 h/(4e2)h/(4e^2)(スピン2つ・向き2つ)の L/ℓL/\ell 倍になる。この h/e2h/e^2 の単位が1次元の伝導の自然な目盛りであることは、Landauer の公式として第4部第7章で扱う。数値では ℓ=2×12×0.25 nm/(0.52×0.01)=200 nm\ell=2\times1^2\times0.25\,\mathrm{nm}/(0.5^2\times0.01)=200\,\mathrm{nm}、kFℓ=(π/2a)ℓ=1257k_{\mathrm F}\ell=(\pi/2a)\ell=1257 で、Born 散乱の数え方が使える範囲にある。1 μm1\,\mathrm{\mu m} の鎖の抵抗は R=(h/4e2)×5=32.3 kΩR=(h/4e^2)\times5=32.3\,\mathrm k\Omega である。ただし1次元では、どんなに弱い乱れでも鎖が ℓ\ell の数倍より長いと干渉で電子が局在するので、この結果は LL が ℓ\ell の数倍までの目安である(4.2.5節で削った干渉)。

(5) (1) で εF\varepsilon _{\mathrm F} を一般の ε\varepsilon にし、τtr(ε)\tau_{\mathrm{tr}}(\varepsilon ) を (3) から入れるとσ(ε)=2e2πℏv(ε)τtr(ε)=e2ℏ v(ε)2πU02ca\sigma(\varepsilon )=\frac{2e^2}{\pi\hbar}v(\varepsilon )\tau_{\mathrm{tr}}(\varepsilon )=\frac{e^2\hbar\,v(\varepsilon )^2}{\pi U_0^2ca}。ε=−2tcos⁡(ka)\varepsilon =-2t\cos(ka) から v2=(2ta/ℏ)2sin⁡2(ka)=(2ta/ℏ)2(1−ε2/4t2)v^2=(2ta/\hbar)^2\sin^2(ka)=(2ta/\hbar)^2(1-\varepsilon ^2/4t^2) なので

σ(ε)=σ(0)(1−ε24t2).\sigma(\varepsilon )=\sigma(0)\Bigl(1-\frac{\varepsilon ^2}{4t^2}\Bigr).

半分詰めでは状態密度が ε=0\varepsilon =0 について対称なので、化学ポテンシャルは温度によらず μ=0\mu=0 に留まる。Sommerfeld 展開から

σ(T)=∫σ(ε)(−∂f0∂ε)dε=σ(0)+π26(kBT)2σ′′(0)=σ(0)[1−π212(kBTt)2]\sigma(T)=\int\sigma(\varepsilon )\Bigl(-\frac{\partial f_0}{\partial \varepsilon }\Bigr)\mathrm{d}\varepsilon =\sigma(0)+\frac{\pi^2}{6}( k_{\mathrm B} T)^2\sigma''(0)=\sigma(0)\Bigl[1-\frac{\pi^2}{12}\Bigl(\frac{ k_{\mathrm B} T}{t}\Bigr)^2\Bigr]

で、σ(ε)\sigma(\varepsilon ) が2次式なのでこの展開は厳密である(kBT=0.05tk_{\mathrm B} T=0.05t の数値積分と 10−510^{-5} より良く一致)。t=1 eVt=1\,\mathrm{eV}、T=300 KT=300\,\mathrm K では補正は 5.5×10−45.5\times10^{-4} で、不純物だけの散乱では抵抗率はほとんど温度によらない。これが Matthiessen 則 (4.2.18) の残留抵抗 ρ0\rho_0 が一定とみなせる理由である。ここで使った σ(ε)\sigma(\varepsilon ) の ε\varepsilon についての傾きは、第4部第3章で熱電能(Mott の公式)を決める量になる。この半分詰めの鎖では傾きが 0 なので、熱電能も最低次で 0 になる。

木下の勉強ステータス

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

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

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