tecktom.com

この章で選ぶもの:バンドの電子と不純物、占有率 f(r,k,t)f(\bm{r},\bm{k},t)

2026年10月15日12分で読めますAI生成教科書『虎の巻』第4部 第2章 つなぎ:Boltzmann 方程式と緩和時間近似の電気伝導
目次
  1. 選ぶ
  2. 借りるもの:半古典運動と Liouville の定理
  3. Boltzmann 方程式
  4. 木下の層

AI 生成・木下が未勉強

ハミルトニアンの選び方で固体を考える 第4部 第2章 §4.2.1 ★★★ 選ぶのは「周期ポテンシャルのバンド+ばらまかれた不純物+ゆっくり変わる外場」の1電子の H^\hat H。半古典運動と Liouville の定理を借りると、占有率 f(r,k,t)f(\bm{r},\bm{k},t) の Boltzmann 方程式が書ける

前の第4部第1章では、金属の伝導電子を自由電子気体とみなし、Fermi 分布の下の積分を Sommerfeld 展開で kBT/εFk_{\mathrm B} T/\varepsilon _{\mathrm F} の冪に展開した。そこで扱ったのは平衡の性質(比熱・化学ポテンシャル)である。この章では電場をかけて電流を流す。電流が流れている金属は平衡にないので、電子の分布が平衡の Fermi 分布からどれだけずれ、そのずれが散乱でどう引き戻されるかを決める方程式が要る。それが Boltzmann 方程式である。

Boltzmann 方程式そのものの導出は、解析力学から様々な理論分野を組み立てる 第3部第8章へ移した。そこでは NN 粒子の Liouville 方程式を粒子1つ分に縮め、希薄な気体で分子的混沌を仮定して衝突の項を書き、衝突の項を −(f−f0)/τ-(f-f_0)/\tau に置き換える緩和時間近似で Drude の伝導度 ne2τ/mne^2\tau/m まで進んだ。バンドのエネルギー εn(k)\varepsilon _n(\bm{k}) を Hamiltonian とみなす半古典運動は、解析力学から様々な理論分野を組み立てる 第4部第10章へ移した。この章の最初の節は、その2つから何を借りるかを1頁にまとめる「つなぎ」である(4.2.1節)。

残りの節で、この巻の主題である「どの H^\hat H を選び、何を削ったか」に沿って金属の電気伝導を計算する。当てる道具は、不純物による弾性散乱の衝突の項を Fermi の黄金律で書き、電場の1次まで線形にした Boltzmann 方程式である(4.2.2節)。ここで、散乱の角度に 1−cos⁡θ1-\cos\theta の重みを付けた輸送緩和時間 τtr\tau_{\mathrm{tr}} が現れ、等方なバンドでは緩和時間近似が近似でなく厳密な帰結になる。出てくるものの第一は、伝導度が Fermi 面の上の積分で書けること、満ちたバンドは電流を運ばないこと、縮退していない電子では緩和時間の平均の取り方が変わることである(4.2.3節)。第二は、交流の伝導度 σ(ω)=σ0/(1−iωτ)\sigma(\omega)=\sigma_0/(1-\mathrm{i}\omega\tau) と、τ\tau によらない和則、プラズマ振動数である(4.2.4節)。そのあと、削ったものを並べて近似の裏付けの札を貼り、抵抗率の温度変化で確かめる(4.2.5節)。最後に授業の Drude 模型・移動度の式と並べ(4.2.6節)、まとめと総合問題で終える。熱の流れを含む輸送(熱電効果・Wiedemann–Franz 則)は第4部第3章、磁場の中の輸送(Hall 効果)は第4部第4章で扱う。

この章では、電子の電荷を −e-e(e>0e>0)、抵抗率を ρ\rho と書く(この章に相空間の密度は出ない)。1電子のエネルギーは ε\varepsilon、化学ポテンシャルは μ\mu で、移動度には必ず添字を付けて μe\mu_{\mathrm e}・μh\mu_{\mathrm h} と書く。状態密度 D(ε)D(\varepsilon ) は記法の正本どおり単位体積・両スピンあたりで、拡散係数は引数なしの DD と書く。

選ぶ H^\hat H

この章で選ぶのは、電子1個の Hamiltonian

H^=H^cr+∑ju(r^−Rj)−eϕ(r^,t)(4.2.1)\hat H=\hat H_{\mathrm{cr}}+\sum_{j}u(\hat{\bm{r}}-\bm{R}_j)-e\phi(\hat{\bm{r}},t)\tag{4.2.1}

である。第1項 H^cr\hat H_{\mathrm{cr}} は完全な結晶の周期ポテンシャルの中の電子の Hamiltonian で、第2部で見たとおりバンド εn(k)\varepsilon _n(\bm{k}) と Bloch 関数 ψnk\psi_{n\bm{k}} で対角化される。第2項は位置 Rj\bm{R}_j にばらまかれた不純物(または欠陥)のポテンシャルで、不純物の数密度を nin_{\mathrm i} と書く。不純物は動かないとし、その位置はでたらめとする。第3項は外から加えた電場のポテンシャル ϕ\phi によるエネルギー(電子の電荷 −e-e)である。磁場を入れるときは運動量をゲージ場の分だけずらすが、この章の計算では磁場を 0 にする。電子どうしの相互作用と、原子核の振動(フォノン)は、この H^\hat H には入れていない。フォノンによる散乱は4.2.5節で衝突の項にだけ戻す。

借りるもの:半古典運動と Liouville の定理

この H^\hat H の Schrödinger 方程式を解く代わりに、電子を Bloch 関数を重ねた波束とみなし、その中心の位置 r\bm{r} と波数 k\bm{k} の運動を追う。外場が格子間隔に比べてゆっくり変わり、1つのバンドの中に留まるとき、波束の中心は

r˙=vn(k)=1ℏ∇kεn(k),ℏk˙=−e(E+vn×B)(4.2.2)\dot{\bm{r}}=\bm{v}_n(\bm{k})=\frac{1}{\hbar}\nabla_{\bm{k}}\varepsilon _n(\bm{k}),\qquad \hbar\dot{\bm{k}}=-e\bigl(\bm{E}+\bm{v}_n\times\bm{B}\bigr)\tag{4.2.2}

に従う。これは εn(k)−eϕ(r)\varepsilon _n(\bm{k})-e\phi(\bm{r}) を Hamiltonian、(r,ℏk)(\bm{r},\hbar\bm{k}) を正準変数とした Hamilton 方程式で、導出は解析力学から様々な理論分野を組み立てる 第4部第10章にある。Hamilton の流れなので、(r,k)(\bm{r},\bm{k}) の空間の体積は運動で保たれる(Liouville の定理。磁場があっても ∇k⋅(v×B)=0\nabla_{\bm{k}}\cdot(\bm{v}\times\bm{B})=0 なので保たれることは解析力学から様々な理論分野を組み立てる 第3部第8章で確かめた)。

相空間の体積 d3r d3k\mathrm{d}^3r\,\mathrm{d}^3k の中にある1電子状態の数は、スピン1つあたり d3r d3k/(2π)3\mathrm{d}^3r\,\mathrm{d}^3k/(2\pi)^3 である(体積 VV の結晶では k\bm{k} が (2π)3/V(2\pi)^3/V ごとに1つ許されるため。第1部第2章)。そこで、この章では分布関数 f(r,k,t)f(\bm{r},\bm{k},t) を「その状態が電子で占められている確率」(占有率、0≤f≤10\le f\le1)として定める。電子密度と電流密度は

n(r,t)=2∫d3k(2π)3 f(r,k,t),J(r,t)=−2e∫d3k(2π)3 v(k) f(r,k,t)(4.2.3)n(\bm{r},t)=2\int\frac{\mathrm{d}^3k}{(2\pi)^3}\,f(\bm{r},\bm{k},t),\qquad \bm{J}(\bm{r},t)=-2e\int\frac{\mathrm{d}^3k}{(2\pi)^3}\,\bm{v}(\bm{k})\,f(\bm{r},\bm{k},t)\tag{4.2.3}

で、先頭の 22 はスピンの向き、積分は第1 Brillouin 域にわたる。解析力学から様々な理論分野を組み立てる 第3部第8章の ff は相空間の密度(2/h32/h^3 を含む)だったので、この章の ff はそれを 2/h32/h^3 で割った量にあたる。占有率で書くと、平衡の分布は Fermi 分布 f0(ε)=1/(eβ(ε−μ)+1)f_0(\varepsilon )=1/(\mathrm{e}^{\beta(\varepsilon -\mu)}+1) になり、Pauli の排他律が「f≤1f\le1」として見える。

Boltzmann 方程式

散乱が無ければ、各電子は式 (4.2.2) の流れに沿って動き、Liouville の定理から ff は流れに沿って一定である。散乱はこれに飛び移りの項を足す。したがって

∂f∂t+v⋅∇rf−eℏ(E+v×B)⋅∇kf=(∂f∂t)coll(4.2.4)\frac{\partial f}{\partial t}+\bm{v}\cdot\nabla_{\bm{r}}f-\frac{e}{\hbar}\bigl(\bm{E}+\bm{v}\times\bm{B}\bigr)\cdot\nabla_{\bm{k}}f=\Bigl(\frac{\partial f}{\partial t}\Bigr)_{\mathrm{coll}}\tag{4.2.4}

が成り立つ。これが固体の電子の Boltzmann 方程式である。左辺は ff を式 (4.2.2) の軌道に沿って時間で微分したもの(df/dt=∂f/∂t+r˙⋅∇rf+k˙⋅∇kf\mathrm{d} f/\mathrm{d} t=\partial f/\partial t+\dot{\bm{r}}\cdot\nabla_{\bm{r}}f+\dot{\bm{k}}\cdot\nabla_{\bm{k}}f)で、右辺の衝突の項が不純物などによる散乱を表す。解析力学から様々な理論分野を組み立てる 第3部第8章の気体との違いは、粒子どうしの衝突でなく、散乱の相手が主に動かない不純物であることである。不純物は電子の運動量を受け取れるので、電子の全運動量は保存せず、そのために電流が減衰して抵抗が生じる。

ここで借りたものと、この章で計算するもの

解析力学から様々な理論分野を組み立てる 第3部第8章からは、Boltzmann 方程式の形(Liouville+衝突)、緩和時間近似 (∂f/∂t)coll≈−(f−f0)/τ(\partial f/\partial t)_{\mathrm{coll}}\approx-(f-f_0)/\tau、それで解いた伝導度 σij=e2τn⟨(m∗−1)ij⟩\sigma_{ij}=e^2\tau n\langle (m^{*-1})_{ij} \rangle と Drude の伝導度 ne2τ/mne^2\tau/m を借りる。解析力学から様々な理論分野を組み立てる 第4部第10章からは式 (4.2.2) を借りる。この章で新しく計算するのは、不純物の衝突の項から τ\tau が何で決まるか、Fermi 面の積分の形、満ちたバンドと正孔、交流の伝導度、抵抗率の温度変化である。

例題 4.2.1 銅の伝導電子の長さとエネルギーの順序 ★★

オリジナル銅の伝導電子を自由電子とみなし、電子密度 n=8.47×1028 m−3n=8.47\times10^{28}\,\mathrm{m^{-3}}、室温(293 K293\,\mathrm K)の抵抗率 ρ=1.68×10−8 Ω m\rho=1.68\times10^{-8}\,\Omega\,\mathrm m を与える。

  1. (1) Drude の式 ρ=m/(ne2τ)\rho=m/(ne^2\tau) から τ\tau を求め、Fermi 波長 λF=2π/kF\lambda_{\mathrm F}=2\pi/k_{\mathrm F}、Fermi 速度 vFv_{\mathrm F}、平均自由行程 ℓ=vFτ\ell=v_{\mathrm F}\tau を求めよ。
  2. (2) 散乱による準位の幅 ℏ/τ\hbar/\tau を εF\varepsilon _{\mathrm F} と比べ、その比が 2/(kFℓ)2/(k_{\mathrm F}\ell) に等しいことを示せ。
  3. (3) 電流密度 1 A mm−21\,\mathrm{A\,mm^{-2}} を流すときの電場 EE、ドリフト速度 vd=J/(ne)v_{\mathrm d}=J/(ne)、平均自由行程の間に電場がする仕事 eEℓeE\ell を求め、vFv_{\mathrm F}・kBTk_{\mathrm B} T と比べよ。
解答

(1) τ=m/(ne2ρ)=2.49×10−14 s\tau=m/(ne^2\rho)=2.49\times10^{-14}\,\mathrm s。kF=(3π2n)1/3=1.36×1010 m−1k_{\mathrm F}=(3\pi^2n)^{1/3}=1.36\times10^{10}\,\mathrm{m^{-1}} で λF=0.462 nm\lambda_{\mathrm F}=0.462\,\mathrm{nm}、vF=ℏkF/m=1.57×106 m s−1v_{\mathrm F}=\hbar k_{\mathrm F}/m=1.57\times10^6\,\mathrm{m\,s^{-1}}、ℓ=3.9×10−8 m=39 nm\ell=3.9\times10^{-8}\,\mathrm m=39\,\mathrm{nm}。

(2) ℏ/τ=0.0264 eV\hbar/\tau=0.0264\,\mathrm{eV}、εF=ℏ2kF2/(2m)=7.03 eV\varepsilon _{\mathrm F}=\hbar^2k_{\mathrm F}^2/(2m)=7.03\,\mathrm{eV} で、比は 3.75×10−33.75\times10^{-3}。εF=12ℏkFvF\varepsilon _{\mathrm F}=\frac12\hbar k_{\mathrm F}v_{\mathrm F} を使うと ℏ/(εFτ)=2/(kFvFτ)=2/(kFℓ)\hbar/(\varepsilon _{\mathrm F}\tau)=2/(k_{\mathrm F}v_{\mathrm F}\tau)=2/(k_{\mathrm F}\ell) で、kFℓ=533k_{\mathrm F}\ell=533 である。

(3) E=ρJ=1.68×10−2 V m−1E=\rho J=1.68\times10^{-2}\,\mathrm{V\,m^{-1}}、vd=7.37×10−5 m s−1v_{\mathrm d}=7.37\times10^{-5}\,\mathrm{m\,s^{-1}} で vd/vF=4.7×10−11v_{\mathrm d}/v_{\mathrm F}=4.7\times10^{-11}。eEℓ=6.6×10−10 eVeE\ell=6.6\times10^{-10}\,\mathrm{eV} で、kBT=0.0252 eVk_{\mathrm B} T=0.0252\,\mathrm{eV} の約 4×1074\times10^7 分の1である。

この3つの比が、この章の計算の前提をそれぞれ数で示している。λF≪ℓ\lambda_{\mathrm F}\ll\ell なので、波束の広がりを ℓ\ell より小さく取って位置と波数を同時に指定でき、式 (4.2.4) の半古典の描像を使うことができる。ℏ/τ≪εF\hbar/\tau\ll\varepsilon _{\mathrm F} は散乱がまれで、散乱と散乱の間は Bloch 状態がよく定まっていることを表す。eEℓ≪kBTeE\ell\ll k_{\mathrm B} T は、分布が平衡からわずかしかずれないので、電場の1次で十分なことを表す。

木下の勉強ステータス

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

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

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