tecktom.com

出てくるもの2:動的な遮蔽、プラズモンと粒子–正孔の連続帯

2026年10月15日10分で読めますAI生成教科書『虎の巻』第5部 第1章 遮蔽と Fermi 液体
目次
  1. 高い振動数の極限
  2. プラズモン
  3. 粒子–正孔の連続帯と Landau 減衰
  4. 木下の層

AI 生成・木下が未勉強

ハミルトニアンの選び方で固体を考える 第5部 第1章 §5.1.4 ★★★ 速く変わる場には電子が追いつけず、遮蔽が起きなくなる。誘電関数の零点は電子気体全体の縦の振動(プラズモン)で、Coulomb 相互作用の長距離に及ぶ部分はこの高いエネルギーの集団運動に取り込まれる

高い振動数の極限

振動数が ω≫qvF\omega\gg qv_{\mathrm F} のとき、式 (5.1.5) を ℏω\hbar\omega で展開する。第2項を k→k−q\bm{k}\to\bm{k}-\bm{q} と付け替えるとχ0=1V∑kσfk[1ℏω−a−1ℏω−b]\chi_0=\frac1V\sum_{\bm{k}\sigma}f_{\bm{k}}\bigl[\frac{1}{\hbar\omega-a}-\frac{1}{\hbar\omega-b}\bigr]、a=εk+q−εka=\varepsilon _{\bm{k}+\bm{q}}-\varepsilon _{\bm{k}}、b=εk−εk−qb=\varepsilon _{\bm{k}}-\varepsilon _{\bm{k}-\bm{q}} と書ける。1ℏω−a=1ℏω+a(ℏω)2+a2(ℏω)3+a3(ℏω)4+⋯\frac{1}{\hbar\omega-a}=\frac{1}{\hbar\omega}+\frac{a}{(\hbar\omega)^2}+\frac{a^2}{(\hbar\omega)^3}+\frac{a^3}{(\hbar\omega)^4}+\cdots を使うと、角かっこの中はa−b(ℏω)2+a2−b2(ℏω)3+a3−b3(ℏω)4+⋯\frac{a-b}{(\hbar\omega)^2}+\frac{a^2-b^2}{(\hbar\omega)^3}+\frac{a^3-b^3}{(\hbar\omega)^4}+\cdots になる。a−b=ℏ2q2/ma-b=\hbar^2q^2/m は k\bm{k} によらず、a2−b2a^2-b^2 は Fermi 球について平均すると 00(k⋅q\bm{k}\cdot\bm{q} の奇関数)、a3−b3a^3-b^3 の平均は ℏ2q2m⋅35ℏ2q2vF2\frac{\hbar^2q^2}{m}\cdot\frac35\hbar^2q^2v_{\mathrm F}^2(q6q^6 の項を除く)である。したがって

χ0(q,ω)≈nq2mω2[1+35q2vF2ω2+⋯ ],εr(q,ω)≈1−ωp2ω2[1+35q2vF2ω2],ωp2=ne2ε0m(5.1.14)\chi_0(q,\omega)\approx\frac{nq^2}{m\omega^2}\Bigl[1+\frac35\frac{q^2v_{\mathrm F}^2}{\omega^2}+\cdots\Bigr],\qquad \varepsilon _{\mathrm r}(q,\omega)\approx1-\frac{\omega_{\mathrm p}^2}{\omega^2}\Bigl[1+\frac35\frac{q^2v_{\mathrm F}^2}{\omega^2}\Bigr],\qquad \omega_{\mathrm p}^2=\frac{ne^2}{\varepsilon _0m}\tag{5.1.14}

を得る。q→0q\to0 の εr=1−ωp2/ω2\varepsilon _{\mathrm r}=1-\omega_{\mathrm p}^2/\omega^2 は、第4部第2章で Drude 模型から出した誘電関数で緩和時間を無限大にしたものと同じである。静的な極限では χ0<0\chi_0<0 で遮蔽が起きたのに、高い振動数では χ0>0\chi_0>0 になり、εr<1\varepsilon _{\mathrm r}<1 と符号が逆になる。電子は質量を持つので、速く変わる力には慣性のために遅れて応え、集まるべき所から逆に離れる位相で動く。

プラズモン

外の電荷が無くても ϕ\phi が 00 でない解がありうるのは、式 (5.1.8) で εr(q,ω)=0\varepsilon _{\mathrm r}(q,\omega)=0 のときである。これは電子気体が外から押されなくても自分で続ける振動、つまり電子の密度の縦波である。式 (5.1.14) を 00 とおいて ω2\omega^2 について解き、q2q^2 の1次まで残すと

ω(q)2≈ωp2+35vF2q2(5.1.15)\omega(q)^2\approx\omega_{\mathrm p}^2+\frac35v_{\mathrm F}^2q^2\tag{5.1.15}

になる。q→0q\to0 の振動数は ωp\omega_{\mathrm p} で有限に残り、音波のように 00 へ行かない。電子全体が正の背景に対してずれると、ずれに比例する電場が全体に一様にかかり、ばねの定数が ne2/ε0ne^2/\varepsilon _0 の振動子になるからである。この集団運動を量子化したものをプラズモンと呼ぶ。

Coulomb 相互作用の長距離に及ぶ部分(小さい qq)の効果は、プラズモンの形で現れる。小さい qq の電荷の揺らぎは、エネルギー ℏωp\hbar\omega_{\mathrm p} の集団運動としてしか存在できず、金属では ℏωp\hbar\omega_{\mathrm p} が εF\varepsilon _{\mathrm F} と同じかそれ以上の大きさなので、低温・低エネルギーの性質には現れない。Fermi 面の近くの電子に残るのは、遮蔽された短距離の相互作用だけである。この役割分担が、次の節で電子を準粒子として扱える理由の1つになる。

粒子–正孔の連続帯と Landau 減衰

χ0\chi_0 の虚部は、分母の ℏω+εk−εk+q\hbar\omega+\varepsilon _{\bm{k}}-\varepsilon _{\bm{k}+\bm{q}} が 00 になる所から出る。これは、エネルギー ℏω\hbar\omega と運動量 ℏq\hbar\bm{q} を受け取って、Fermi 球の中の電子 k\bm{k} を外の k+q\bm{k}+\bm{q} へ移す(粒子と正孔の対を1つ作る)過程である。εk+q−εk=ℏ2m(k⋅q)+ℏ2q22m\varepsilon _{\bm{k}+\bm{q}}-\varepsilon _{\bm{k}}=\frac{\hbar^2}{m}(\bm{k}\cdot\bm{q})+\frac{\hbar^2q^2}{2m} で ∣k∣<kF\lvert \bm{k}\rvert <k_{\mathrm F} だから、絶対零度でこの対が作れるのは

max⁡(0, ℏ2q22m−ℏvFq)≤ℏω≤ℏvFq+ℏ2q22m(5.1.16)\max\Bigl(0,\ \frac{\hbar^2q^2}{2m}-\hbar v_{\mathrm F}q\Bigr)\le\hbar\omega\le\hbar v_{\mathrm F}q+\frac{\hbar^2q^2}{2m}\tag{5.1.16}

の範囲である(下限の 00 は、q<2kFq<2k_{\mathrm F} なら Fermi 面の近くで小さなエネルギーの対が作れることによる)。この範囲を粒子–正孔の連続帯と呼ぶ。連続帯の外では χ0\chi_0 は実数で、プラズモンは減衰しない。qq を大きくするとプラズモンの分散 (5.1.15) は連続帯の上の縁 ℏvFq+ℏ2q2/(2m)\hbar v_{\mathrm F}q+\hbar^2q^2/(2m) に追いつき、その先ではプラズモンが1つの粒子–正孔対に崩れることができて、急に減衰する。これを Landau 減衰と呼ぶ。

図

図 5.1.3 アルミニウムの密度(rs/a0=2.07r_{\mathrm s}/a_0=2.07、ℏωp/εF=1.354\hbar\omega_{\mathrm p}/\varepsilon _{\mathrm F}=1.354)の電子気体の励起。灰色が粒子–正孔の連続帯 (5.1.16)、実線が小さい qq の近似 (5.1.15)、点線が Lindhard 関数で εr=0\varepsilon _{\mathrm r}=0 を数値で解いた RPA のプラズモン。近似は q=0.635kFq=0.635k_{\mathrm F}、RPA の値は q=0.739kFq=0.739k_{\mathrm F} で連続帯に入り、その先は Landau 減衰で鋭い励起でなくなる。

例題 5.1.4 アルミニウムのプラズモン ★★

オリジナルアルミニウムは fcc(格子定数 0.405 nm0.405\,\mathrm{nm}、与えられた値)で、1原子あたり3個の伝導電子を持つ。

  1. (1) nn, kFk_{\mathrm F}, εF\varepsilon _{\mathrm F}, rs/a0r_{\mathrm s}/a_0 を求めよ。
  2. (2) ℏωp\hbar\omega_{\mathrm p} を求めよ。電子線のエネルギー損失の測定では、アルミニウムの体積のプラズモンは 15.3 eV15.3\,\mathrm{eV} に見える(C. J. Powell and J. B. Swan (1959), Phys. Rev. 115, 869)。比を求めよ。
  3. (3) 近似 (5.1.15) を ℏω=ℏωp+Aq2\hbar\omega=\hbar\omega_{\mathrm p}+Aq^2 と書いたときの AA を eV A˚2\mathrm{eV\,\AA^2} で求めよ。
  4. (4) 近似 (5.1.15) が連続帯の上の縁に達する qcq_{\mathrm c} を kFk_{\mathrm F} 単位で求め、そのときの ℏω\hbar\omega を eV で答えよ。
解答

(1) fcc の単位胞には4原子あるので n=4×3/(0.405×10−9)3=1.81×1029 m−3n=4\times3/(0.405\times10^{-9})^3=1.81\times10^{29}\,\mathrm{m^{-3}}。kF=(3π2n)1/3=1.75×1010 m−1k_{\mathrm F}=(3\pi^2n)^{1/3}=1.75\times10^{10}\,\mathrm{m^{-1}}、εF=ℏ2kF2/(2m)=11.65 eV\varepsilon _{\mathrm F}=\hbar^2k_{\mathrm F}^2/(2m)=11.65\,\mathrm{eV}、rs=(3/(4πn))1/3=2.07a0r_{\mathrm s}=(3/(4\pi n))^{1/3}=2.07a_0。

(2) ωp=ne2/(ε0m)\omega_{\mathrm p}=\sqrt{ne^2/(\varepsilon _0m)} から ℏωp=15.8 eV\hbar\omega_{\mathrm p}=15.8\,\mathrm{eV}。実測との比は 15.3/15.8=0.9715.3/15.8=0.97 である。周期ポテンシャルも内殻の電子も無視した模型で、集団運動のエネルギーが 3%3\% の範囲で合う。ℏωp/εF=1.354\hbar\omega_{\mathrm p}/\varepsilon _{\mathrm F}=1.354 で、プラズモンは、Fermi 面の近くの熱的な励起(kBT≈0.026 eVk_{\mathrm B} T\approx0.026\,\mathrm{eV})よりはるかに高い所にある。

(3) ω≈ωp+3vF210ωpq2\omega\approx\omega_{\mathrm p}+\frac{3v_{\mathrm F}^2}{10\omega_{\mathrm p}}q^2 から A=3ℏvF210ωpA=\frac{3\hbar v_{\mathrm F}^2}{10\omega_{\mathrm p}}。vF=ℏkF/m=2.02×106 m/sv_{\mathrm F}=\hbar k_{\mathrm F}/m=2.02\times10^6\,\mathrm{m/s} を入れるとA=3.38 eV A˚2A=3.38\,\mathrm{eV\,\AA^2}。q=0.1 A˚−1q=0.1\,\mathrm{\AA^{-1}} で 0.034 eV0.034\,\mathrm{eV} 上がる程度の、ゆるい分散である。

(4) x=q/kFx=q/k_{\mathrm F}、W=ℏω/εFW=\hbar\omega/\varepsilon _{\mathrm F} とすると、ℏvFq=2xεF\hbar v_{\mathrm F}q=2x\varepsilon _{\mathrm F}、ℏ2q2/(2m)=x2εF\hbar^2q^2/(2m)=x^2\varepsilon _{\mathrm F} で、近似は W2=1.3542+125x2W^2=1.354^2+\frac{12}{5}x^2、上の縁は W=2x+x2W=2x+x^2。等しいとおいて数値で解くと x=0.635x=0.635、W=1.674W=1.674、ℏω=19.5 eV\hbar\omega=19.5\,\mathrm{eV} である。Lindhard 関数で εr=0\varepsilon _{\mathrm r}=0 を解いた RPA の値は x=0.739x=0.739 で少し先になる(図 5.1.3)。波長にして 2π/qc≈0.6 nm2\pi/q_{\mathrm c}\approx0.6\,\mathrm{nm} より短い電荷の揺らぎは、もはや集団運動として振る舞わない。

ほかの巻では

木下の勉強ステータス

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

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

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