ハミルトニアンの選び方で固体を考える 第5部 第1章 §5.1.6 ★★ RPA は誘起された電荷の交換と相関を、ジェリウムはイオンを削った。イオンを戻すと、遮蔽されたイオンの振動が音響フォノンになる。遮蔽・プラズモン・Fermi 液体はそれぞれ電子線の損失・中性子・STM・比熱と磁化率で確かめられる
削ったもの
1つめは誘起された電荷の交換と相関 である。RPA の仮定 (5.1.7) では、電子は平均の場 ϕ \phi ϕ に独立に応える。実際には、ある電子のまわりには交換孔と相関孔があり、その電子が感じる場は平均の場より弱い。この効果を入れるには ε r \varepsilon _{\mathrm r} ε r の v q v_{\bm{q}} v q に補正因子(局所場の補正)を掛ける。Gell-Mann と Brueckner は、高い密度の極限r s / a 0 → 0 r_{\mathrm s}/a_0\to0 r s / a 0 → 0 で、RPA で基底エネルギーの相関の部分の主要項が正しく得られることを示した(M. Gell-Mann and K. A. Brueckner (1957), Phys. Rev. 106, 364)。これは摂動級数の最も発散の強い項を全部の次数で足し上げた結果で、金属の密度 r s / a 0 = 2 r_{\mathrm s}/a_0=2 r s / a 0 = 2 〜6 6 6 はこの極限から遠い。
近似の裏付け
型 B(極限) | 小さいもの r s / a 0 → 0 r_{\mathrm s}/a_0\to0 r s / a 0 → 0 (高い密度。Coulomb エネルギーが運動エネルギーに比べて小さい) | 証明 要確認(Gell-Mann–Brueckner は摂動級数の主要項の再和としての結果。誤差の大きさを定理として評価した原論文は確かめていない) | 行き先 近似の数理 第6章 ・近似の数理 第13章
近似の裏付け
型 B(極限) | 小さいもの q / k F → 0 q/k_{\mathrm F}\to0 q / k F → 0 (ポテンシャルの変化が Fermi 波長に比べてゆっくり。Thomas–Fermi) | 証明 証明済み(RPA の中では Lindhard 関数の極限 L ( 0 ) = 1 L(0)=1 L ( 0 ) = 1 として本文で示した。RPA そのものの正しさは上の札) | 行き先 近似の数理 第6章
2つめはイオンの周期ポテンシャルと運動 である。周期ポテンシャルを戻すと、ε r \varepsilon _{\mathrm r} ε r は q \bm{q} q と q + G \bm{q}+\bm{G} q + G を混ぜる行列になり、バンドの間の遷移が ε r ( ω ) \varepsilon _{\mathrm r}(\omega) ε r ( ω ) に加わる(第4部第6章 )。イオンの運動を戻すと、イオンもまた電荷を持つ粒子の気体で、それ自身のプラズマ振動数 Ω p 2 = n i o n ( Z e ) 2 / ( ε 0 M ) \Omega_{\mathrm p}^2=n_{\mathrm{ion}}(Ze)^2/(\varepsilon _0M) Ω p 2 = n ion ( Z e ) 2 / ( ε 0 M ) を持つ(Z Z Z はイオンの価数、M M M はその質量)。ただしイオンはゆっくり動くので、電子はその動きに追いついて静的に遮蔽する。イオンの振動数は Ω p 2 / ε r ( q , 0 ) \Omega_{\mathrm p}^2/\varepsilon _{\mathrm r}(q,0) Ω p 2 / ε r ( q , 0 ) になり、Thomas–Fermi の ε r ≈ q T F 2 / q 2 \varepsilon _{\mathrm r}\approx q_{\mathrm{TF}}^2/q^2 ε r ≈ q TF 2 / q 2 (q ≪ q T F q\ll q_{\mathrm{TF}} q ≪ q TF )を入れると ω = ( Ω p / q T F ) q \omega=(\Omega_{\mathrm p}/q_{\mathrm{TF}})q ω = ( Ω p / q TF ) q と、波数に比例する音波になる。その速さは
v s = Ω p q T F = Z m 3 M v F (5.1.22) v_{\mathrm s}=\frac{\Omega_{\mathrm p}}{q_{\mathrm{TF}}}=\sqrt{\frac{Zm}{3M}}\,v_{\mathrm F}\tag{5.1.22} v s = q TF Ω p = 3 M Z m v F ( 5.1.22 )
である(Bohm–Staver の関係)。裸のイオンの気体のプラズマ振動が、電子の遮蔽によって音響フォノン(第3部第1章 )に変わる。同じ遮蔽されたイオンの振動を介して、電子どうしには遮蔽された Coulomb の斥力と、フォノンを交換する引力が同時に働く。後者が勝つ低いエネルギーの領域で電子は対を作る(第5部第4章 )。
3つめはFermi 液体であること自体 である。断熱接続は、途中で相転移が起きないことを仮定した。1次元の電子は、どれほど弱い相互作用でも Fermi 液体にならない(準粒子の代わりに電荷とスピンの集団運動が残る。Tomonaga–Luttinger 液体)。2次元・3次元でも、超伝導や磁性への不安定性(F 0 a → − 1 F^{\mathrm a}_0\to-1 F 0 a → − 1 など)の近くでは準粒子の描像が破れる。Fermi 液体が数学として証明されているのは限られた場合だけである。
近似の裏付け
型 B(極限) | 小さいもの Fermi 面の近くの励起だけを残すくりこみ群(ξ / ε F → 0 \xi/\varepsilon _{\mathrm F}\to0 ξ / ε F → 0 ・k B T / ε F → 0 k_{\mathrm B} T/\varepsilon _{\mathrm F}\to0 k B T / ε F → 0 ) | 証明 条件付き(2次元で、Fermi 面が強く凸で k → − k \bm{k}\to-\bm{k} k → − k の反転について対称でなく、相互作用が短距離で弱い場合に、絶対零度の摂動展開が収束することが証明されている〔J. Feldman, H. Knörrer, E. Trubowitz の連作 A Two Dimensional Fermi Liquid。第1部は Commun. Math. Phys. 247 , 1 (2004)〕。長距離の Coulomb 相互作用・3次元・一般の Fermi 面では未解決) | 行き先 近似の数理 第13章
何の実験の量で確かめるか
プラズモン は、薄い金属の膜を通り抜けた電子のエネルギーの損失として見える。例題 5.1.4 のアルミニウムの 15.3 e V 15.3\,\mathrm{eV} 15.3 eV がその例で、同じ測定では表面のプラズモン(ω p / 2 \omega_{\mathrm p}/\sqrt2 ω p / 2 に近い)も見える。損失の角度を変えると q \bm{q} q が変わり、分散 (5.1.15) と Landau 減衰で鋭い損失が消えていく様子を測ることができる。
Kohn 異常 は、中性子の非弾性散乱で測ったフォノンの分散に、q = 2 k F q=2k_{\mathrm F} q = 2 k F に当たる位置の小さな折れ曲がりとして見える。鉛のフォノンの分散でこの折れ曲がりが見つかり、その位置がほぼ自由な電子の Fermi 面から予想される位置とよく合うことが報告されている(B. N. Brockhouse, T. Arase, G. Caglioti, K. R. Rao, A. D. B. Woods (1962), Phys. Rev. 128, 1099)。
Friedel 振動 は、走査トンネル顕微鏡(第4部第7章 )で表面の電子の局所的な状態密度を写すと、表面の段差や吸着原子のまわりの同心円の縞として直接見える。銅の (111) 面の表面に束縛された2次元の電子の定在波がその例である(M. F. Crommie, C. P. Lutz, D. M. Eigler (1993), Nature 363, 524)。
Fermi 液体の数 は、低温の比熱の γ \gamma γ から m ∗ / m m^*/m m ∗ / m を、スピンの磁化率(Knight シフトや伝導電子のスピン共鳴)から 1 + F 0 a 1+F^{\mathrm a}_0 1 + F 0 a を読む。2つを組み合わせた Wilson 比が(5.1.21) の量である。準粒子どうしの散乱による電気抵抗は式 (5.1.17) の T 2 T^2 T 2 に従う(運動量を格子に渡す Umklapp の散乱が要る)。液体の 3 H e {}^3\mathrm{He} 3 He は格子もイオンも無い中性の Fermi 粒子の系で、比熱・磁化率・音速から Landau パラメータが一通り決められ、Landau の理論が最初に確かめられた対象である。
例題 5.1.6 銅の電子どうしの散乱の時間 ★★
オリジナル 銅(ε F = 7.03 e V \varepsilon _{\mathrm F}=7.03\,\mathrm{eV} ε F = 7.03 eV )で、式 (5.1.17) の係数を 1 1 1 として準粒子の寿命を見積もる。
(1) 300 K 300\,\mathrm K 300 K で Fermi 面の上の準粒子(ξ = 0 \xi=0 ξ = 0 )の ℏ / τ \hbar/\tau ℏ/ τ と τ \tau τ を求めよ。
(2) 第4部第2章 の室温の銅の緩和時間 τ = 2.49 × 10 − 14 s \tau=2.49\times10^{-14}\,\mathrm s τ = 2.49 × 1 0 − 14 s と比べ、室温の電気抵抗に電子どうしの散乱が重要かどうかを述べよ。
(3) 光電子分光で ξ = 1 e V \xi=1\,\mathrm{eV} ξ = 1 eV の状態を見るとき、ℏ / τ \hbar/\tau ℏ/ τ と ξ \xi ξ の比を求め、準粒子の描像が成り立つかを述べよ。
解答 (1) π k B T = π × 0.02585 e V = 0.0812 e V \pi k_{\mathrm B} T=\pi\times0.02585\,\mathrm{eV}=0.0812\,\mathrm{eV} π k B T = π × 0.02585 eV = 0.0812 eV 、ℏ / τ = 0.0812 2 / 7.03 = 9.4 × 10 − 4 e V \hbar/\tau=0.0812^2/7.03=9.4\times10^{-4}\,\mathrm{eV} ℏ/ τ = 0.081 2 2 /7.03 = 9.4 × 1 0 − 4 eV 、τ = ℏ / ( 9.4 × 10 − 4 e V ) = 7.0 × 10 − 13 s \tau=\hbar/(9.4\times10^{-4}\,\mathrm{eV})=7.0\times10^{-13}\,\mathrm s τ = ℏ/ ( 9.4 × 1 0 − 4 eV ) = 7.0 × 1 0 − 13 s 。
(2) 7.0 × 10 − 13 / 2.49 × 10 − 14 = 28 7.0\times10^{-13}/2.49\times10^{-14}=28 7.0 × 1 0 − 13 /2.49 × 1 0 − 14 = 28 で、電子どうしの散乱はフォノンによる散乱より30倍ほど遅い。しかも電子どうしの散乱は全体の運動量を保存するので、Umklapp でなければ電流を減らさない。室温の抵抗はフォノンで決まり、電子どうしの散乱は T 2 T^2 T 2 の小さな補正として極低温で現れるだけである。
(3) ℏ / τ ≈ ξ 2 / ε F = 1 / 7.03 = 0.14 e V \hbar/\tau\approx\xi^2/\varepsilon _{\mathrm F}=1/7.03=0.14\,\mathrm{eV} ℏ/ τ ≈ ξ 2 / ε F = 1/7.03 = 0.14 eV (( π k B T ) 2 (\pi k_{\mathrm B} T)^2 ( π k B T ) 2 は ξ 2 \xi^2 ξ 2 の 0.7 % 0.7\% 0.7% で無視できる)。ξ \xi ξ の 1 / 7 1/7 1/7 で、スペクトルの山の幅は位置より小さく、準粒子を見分けることができる。ξ \xi ξ を Fermi 準位に近づけるほど比 ξ / ε F \xi/\varepsilon _{\mathrm F} ξ / ε F は小さくなり、準粒子の山は鋭くなる。
授業で習った形との対応
学部の固体物理の授業では、Thomas–Fermi の遮蔽を次のように導くことが多い。静電ポテンシャル ϕ ( r ) \phi(\bm{r}) ϕ ( r ) の中で、電子の電気化学ポテンシャルが一定に保たれるとすると、各点の局所的な化学ポテンシャルは μ + e ϕ ( r ) \mu+e\phi(\bm{r}) μ + e ϕ ( r ) で、電子の密度は n 0 + D ( ε F ) e ϕ n_0+D(\varepsilon _{\mathrm F})e\phi n 0 + D ( ε F ) e ϕ になる。これを Poisson 方程式に入れると
∇ 2 ϕ = q T F 2 ϕ − ρ e x t ε 0 , q T F 2 = e 2 D ( ε F ) ε 0 (5.1.23) \nabla^2\phi=q_{\mathrm{TF}}^2\phi-\frac{\rho_{\mathrm{ext}}}{\varepsilon _0},\qquad q_{\mathrm{TF}}^2=\frac{e^2D(\varepsilon _{\mathrm F})}{\varepsilon _0}\tag{5.1.23} ∇ 2 ϕ = q TF 2 ϕ − ε 0 ρ ext , q TF 2 = ε 0 e 2 D ( ε F ) ( 5.1.23 )
を得る。この章の式 (5.1.9) の (5.1.6) を実空間で書いたものである。「各点で局所的に平衡」という仮定が q ≪ k F q\ll k_{\mathrm F} q ≪ k F の近似に当たり、Lindhard 関数はその仮定を外したものである。
同じ形は、電解質溶液やプラズマ・半導体の Debye–Hückel の遮蔽にも出る。違いは ∂ n / ∂ μ \partial n/\partial\mu ∂ n / ∂ μ だけで、縮退していない(古典の)電子では n ∝ e μ / k B T n\propto\mathrm{e}^{\mu/ k_{\mathrm B} T} n ∝ e μ / k B T から ∂ n / ∂ μ = n / ( k B T ) \partial n/\partial\mu=n/( k_{\mathrm B} T) ∂ n / ∂ μ = n / ( k B T ) 、遮蔽長は λ D = ε 0 k B T / ( n e 2 ) \lambda_{\mathrm D}=\sqrt{\varepsilon _0 k_{\mathrm B} T/(ne^2)} λ D = ε 0 k B T / ( n e 2 ) になる。縮退した電子では ∂ n / ∂ μ = D ( ε F ) \partial n/\partial\mu=D(\varepsilon _{\mathrm F}) ∂ n / ∂ μ = D ( ε F ) で、k B T k_{\mathrm B} T k B T の代わりに ε F \varepsilon _{\mathrm F} ε F 程度のエネルギーが入る。どちらも「q 2 = ( e 2 / ε 0 ) ∂ n / ∂ μ q^2=(e^2/\varepsilon _0)\,\partial n/\partial\mu q 2 = ( e 2 / ε 0 ) ∂ n / ∂ μ 」の1本の式で、第4部第2章 の伝導度と拡散係数の関係 σ = e 2 ( ∂ n / ∂ μ ) D \sigma=e^2(\partial n/\partial\mu)D σ = e 2 ( ∂ n / ∂ μ ) D と同じ ∂ n / ∂ μ \partial n/\partial\mu ∂ n / ∂ μ が現れる。
プラズマ振動は、授業では電子全体を変位 x x x だけずらした板を考え、表面に現れる電荷 ± n e x \pm nex ± n e x の電場 n e x / ε 0 nex/\varepsilon _0 n e x / ε 0 が復元力になるとして m x ¨ = − ( n e 2 / ε 0 ) x m\ddot x=-(ne^2/\varepsilon _0)x m x ¨ = − ( n e 2 / ε 0 ) x から ω p \omega_{\mathrm p} ω p を導く。この章では ε r ( q , ω ) \varepsilon _{\mathrm r}(q,\omega) ε r ( q , ω ) の零点として同じものが出て、さらに q q q を変えたときの分散と、粒子–正孔の連続帯へ入って減衰することまで分かる。
例題 5.1.7 半導体と金属の遮蔽長を比べる ★
オリジナル
(1) 第4部第5章の例題 で、Si(比誘電率 11.7 11.7 11.7 )の中の縮退していない電子(n = 10 22 m − 3 n=10^{22}\,\mathrm{m^{-3}} n = 1 0 22 m − 3 、300 K 300\,\mathrm K 300 K )の Debye 長を λ D = 41 n m \lambda_{\mathrm D}=41\,\mathrm{nm} λ D = 41 nm と求めた。これを例題 5.1.3 の銅の 1 / q T F 1/q_{\mathrm{TF}} 1/ q TF と比べよ。
(2) 銅の伝導電子を仮に古典の気体とみなした Debye 長を 300 K 300\,\mathrm K 300 K で求め、例題 5.1.3 の 1 / q T F 1/q_{\mathrm{TF}} 1/ q TF と比べよ。比を ε F \varepsilon _{\mathrm F} ε F と k B T k_{\mathrm B} T k B T で書け。
解答 (1) 第4部第5章の値 λ D = 4.09 × 10 − 8 m \lambda_{\mathrm D}=4.09\times10^{-8}\,\mathrm m λ D = 4.09 × 1 0 − 8 m を銅の 1 / q T F = 5.52 × 10 − 11 m 1/q_{\mathrm{TF}}=5.52\times10^{-11}\,\mathrm m 1/ q TF = 5.52 × 1 0 − 11 m で割ると 7.4 × 10 2 7.4\times10^2 7.4 × 1 0 2 で、Si の遮蔽長は銅の約 740 倍長い。ドープした半導体の遮蔽長は格子定数(0.54 n m 0.54\,\mathrm{nm} 0.54 nm )の数十倍で、不純物の Coulomb ポテンシャルは長距離に及ぶ。金属では同じポテンシャルが格子間隔より短い距離で消える。
(2) λ D = ε 0 k B T / ( n e 2 ) = 4.1 × 10 − 12 m \lambda_{\mathrm D}=\sqrt{\varepsilon _0 k_{\mathrm B} T/(ne^2)}=4.1\times10^{-12}\,\mathrm m λ D = ε 0 k B T / ( n e 2 ) = 4.1 × 1 0 − 12 m で、1 / q T F = 5.52 × 10 − 11 m 1/q_{\mathrm{TF}}=5.52\times10^{-11}\,\mathrm m 1/ q TF = 5.52 × 1 0 − 11 m の 1 / 13.5 1/13.5 1/13.5 である。比は 1 / q T F λ D = n / ( k B T ) 3 n / ( 2 ε F ) = 2 ε F 3 k B T = 13.5 \frac{1/q_{\mathrm{TF}}}{\lambda_{\mathrm D}}=\sqrt{\frac{n/( k_{\mathrm B} T)}{3n/(2\varepsilon _{\mathrm F})}}=\sqrt{\frac{2\varepsilon _{\mathrm F}}{3 k_{\mathrm B} T}}=13.5 λ D 1/ q TF = 3 n / ( 2 ε F ) n / ( k B T ) = 3 k B T 2 ε F = 13.5 。縮退した電子は Pauli の原理のために Fermi 準位の近くの D ( ε F ) k B T D(\varepsilon _{\mathrm F}) k_{\mathrm B} T D ( ε F ) k B T 程度しか動けず、古典の気体より遮蔽が弱い。それでも金属の遮蔽長は格子間隔より短い。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com