tecktom.com

出てくるもの(2):総和則

2026年10月15日14分で読めますAI生成教科書『虎の巻』第12章 Kramers–Kronig の関係と総和則
目次
  1. 高い振動数での振る舞い
  2. 総和則
  3. 有効電子数
  4. 木下の層

AI 生成・木下が未勉強

光学と量子エレクトロニクス 第12章 §12.4 ★★★ 高い振動数で電子は自由な粒子として振る舞い χ→−ωp2/ω2\chi\to-\omega_{\mathrm p}^2/\omega^2 になる。これを KK の関係と合わせると、吸収スペクトルの重み付きの面積が電子の数密度だけで決まる

高い振動数での振る舞い

光の振動が速すぎて、電子が束縛の力を感じる前に向きを変える振動数を考える。例題 12.1 で見たように、撃力の直後の応答は χ(0+)=0\chi(0^+)=0、χ˙(0+)=ωp2\dot\chi(0^+)=\omega_{\mathrm p}^2 である。この2つは束縛にも減衰にもよらない一般の性質である。撃力の一瞬に電子が受け取るのは運動量 −eE1-eE_1 だけで、位置はまだ動いていないから、物質の中のすべての電子(数密度 NN、内殻の電子も含める)について P˙(0+)=Ne2E1/m\dot P(0^+)=Ne^2E_1/m になる。原子核も同じように動くが、その寄与は質量の比 m/Mm/M だけ小さいので無視する。式 (12.3) を部分積分で 1/ω1/\omega について展開すると

χ(ω)=iχ(0+)ω−χ˙(0+)ω2+O(ω−3)=−ωp2ω2+O(ω−3),ωp2=Ne2ε0m(12.14)\chi(\omega)=\frac{\mathrm{i}\chi(0^+)}{\omega}-\frac{\dot\chi(0^+)}{\omega^2}+\mathcal{O}(\omega^{-3})=-\frac{\omega_{\mathrm p}^2}{\omega^2}+\mathcal{O}(\omega^{-3}),\qquad \omega_{\mathrm p}^2=\frac{Ne^2}{\varepsilon_0m}\tag{12.14}

になる。十分に高い振動数(X 線の領域)では、どの物質も自由電子の気体と同じ誘電関数 1−ωp2/ω21-\omega_{\mathrm p}^2/\omega^2 に近づく。

ff 総和則

一方、式 (12.10) の第1式で ω\omega を大きくすると、ε2\varepsilon_2 が ω\omega より十分低い振動数で減ってしまう限り、分母を −ω2-\omega^2 で置き換えることができ

ε1(ω)−1≈−2πω2∫0∞ω′ε2(ω′) dω′\varepsilon_1(\omega)-1\approx-\frac{2}{\pi\omega^2}\int_0^\infty\omega'\varepsilon_2(\omega')\,\mathrm{d}\omega'

となる。これを式 (12.14) と比べると

総和則

∫0∞ω ε2(ω) dω=π2 ωp2=π2 Ne2ε0m(12.15)\int_0^\infty\omega\,\varepsilon_2(\omega)\,\mathrm{d}\omega=\frac\pi2\,\omega_{\mathrm p}^2=\frac\pi2\,\frac{Ne^2}{\varepsilon_0m}\tag{12.15}

ε1(0)−1=2π∫0∞ε2(ω)ω dω(絶縁体)(12.16)\varepsilon_1(0)-1=\frac2\pi\int_0^\infty\frac{\varepsilon_2(\omega)}{\omega}\,\mathrm{d}\omega\qquad(\text{絶縁体})\tag{12.16}

∫0∞ω Im⁡[−1εr(ω)]dω=π2 ωp2,∫0∞ω κ(ω) dω=π4 ωp2,∫0∞[n(ω)−1]dω=0(12.17)\int_0^\infty\omega\,\operatorname{Im}\Bigl[-\frac{1}{\varepsilon_{\mathrm r}(\omega)}\Bigr]\mathrm{d}\omega=\frac\pi2\,\omega_{\mathrm p}^2,\qquad \int_0^\infty\omega\,\kappa(\omega)\,\mathrm{d}\omega=\frac\pi4\,\omega_{\mathrm p}^2,\qquad \int_0^\infty\bigl[n(\omega)-1\bigr]\mathrm{d}\omega=0\tag{12.17}

を得る。式 (12.15) を ff 総和則と呼ぶ。吸収スペクトルを ω\omega で重み付けた面積は、吸収の形や物質の種類によらず、電子の数密度だけで決まる。吸収の山を幅の狭いものから広いものに変えても、山を別の振動数へ動かしても、その合計の面積を変えることはできない。

同じ式は、上半面で正則な ωχ(ω)\omega\chi(\omega) の、実軸の [−R,R][-R,R] と半径 RR の大きな半円からなる閉曲線の積分からも出る。閉曲線の積分は 00 で、大きな半円の上では式 (12.14) から ωχ≈−ωp2/ω\omega\chi\approx-\omega_{\mathrm p}^2/\omega なので、その積分は −ωp2⋅iπ-\omega_{\mathrm p}^2\cdot\mathrm{i}\pi である。よって ∫−RRωχ dω=iπωp2\int_{-R}^R\omega\chi\,\mathrm{d}\omega=\mathrm{i}\pi\omega_{\mathrm p}^2 で、虚部をとり、被積分関数 ωIm⁡χ\omega\operatorname{Im}\chi が偶関数であることを使うと式 (12.15) になる。

名前の ff は振動子強度に由来する。式 (12.11) の細い線の重みを入れると、左辺は π2ωp2∑jfj\frac\pi2\omega_{\mathrm p}^2\sum_jf_j になるので、式 (12.15) は ∑jfj=1\sum_jf_j=1(電子1個あたり)、原子1個あたりでは振動子強度の和が電子の数 ZZ に等しいことを言っている。これは第1章で引いた Thomas–Reiche–Kuhn の総和則(解析力学から様々な理論分野を組み立てる 第4部第12章)と同じ内容で、量子力学の側では交換関係 [x^, [H^, x^]]=−ℏ2/m\left[\hat x,\,\left[\hat H,\,\hat x\right]\right]=-\hbar^2/m から、光学の側では因果律と高い振動数の自由電子の振る舞いから出る。

伝導度で書くと、εr=1+iσ/(ε0ω)\varepsilon_{\mathrm r}=1+\mathrm{i}\sigma/(\varepsilon_0\omega) から ε2=Re⁡σ/(ε0ω)\varepsilon_2=\operatorname{Re}\sigma/(\varepsilon_0\omega) なので、式 (12.15) は∫0∞Re⁡σ(ω) dω=π2ε0ωp2\int_0^\infty\operatorname{Re}\sigma(\omega)\,\mathrm{d}\omega=\frac\pi2\varepsilon_0\omega_{\mathrm p}^2 になる。ハミルトニアンの選び方で固体を考える 第4部第2章で Drude の伝導度について確かめた和則は、この一般の式の一例である。

式 (12.16) は式 (12.10) の第1式で ω=0\omega=0 とおいたものである。分母の ω′\omega' のために低い振動数の吸収ほど重く数えられ、赤外のわずかな吸収(イオンの振動)が静的な比誘電率を大きく持ち上げることがある。式 (12.17) は、εr\varepsilon_{\mathrm r} の代わりに 1/εr1/\varepsilon_{\mathrm r} と n~\tilde n に同じ議論を当てたものである。どちらも上半面で正則で(εr\varepsilon_{\mathrm r} は上半面で 00 にならない)、高い振動数で 1/εr≈1+ωp2/ω21/\varepsilon_{\mathrm r}\approx1+\omega_{\mathrm p}^2/\omega^2、n~−1≈−ωp2/(2ω2)\tilde n-1\approx-\omega_{\mathrm p}^2/(2\omega^2) になる。最後の ∫(n−1)dω=0\int(n-1)\mathrm{d}\omega=0 は、κ\kappa が ω−1\omega^{-1} より速く減ることから出る(KK の第2式で ω\omega を大きくすると、κ≈2πω∫(n−1)dω′\kappa\approx\frac{2}{\pi\omega}\int(n-1)\mathrm{d}\omega' が残るため)。屈折率は平均すると 11 で、n>1n>1 の領域があれば必ず n<1n<1 の領域がある1屈折率の和則の整理は M. Altarelli, D. L. Dexter, H. M. Nussenzveig, D. Y. Smith, Phys. Rev. B 6, 4502 (1972)。下の脚注へ。n<1n<1 は位相速度が cc を超えることを表すが、信号の先頭が伝わる速さは cc を超えない。n~\tilde n が上半面で正則で、∣ω∣→∞\lvert \omega\rvert \to\infty で 11 に近づくことから、この結論が出る(A. Sommerfeld と L. Brillouin の議論)。

有効電子数

測定した ε2\varepsilon_2 を ω\omega まで積分した量を、式 (12.15) の右辺で割って

Neff(ω)=2πωp2∫0ωω′ε2(ω′) dω′(12.18)N_{\mathrm{eff}}(\omega)=\frac{2}{\pi\omega_{\mathrm p}^2}\int_0^{\omega}\omega'\varepsilon_2(\omega')\,\mathrm{d}\omega'\tag{12.18}

と書く。ここで ωp\omega_{\mathrm p} は全電子の数密度で作る。Neff(ω)N_{\mathrm{eff}}(\omega) は「振動数 ω\omega までの吸収に参加した電子の割合」で、ω→∞\omega\to\infty で 11 になる。図 12.3 は、電子の 1/41/4 が角周波数 ωv\omega_{\mathrm v} の価電子の山、残りの 3/43/4 が 8 ωv8\,\omega_{\mathrm v} の内殻の山に吸収の重みを持つ模型の NeffN_{\mathrm{eff}} である。価電子の山を過ぎたところで NeffN_{\mathrm{eff}} は 0.250.25 の付近で平らになり、内殻の山を過ぎてから 11 に近づく。平らな段の高さから、光学で測っただけで「この振動数の範囲で応答している電子が原子あたり何個か」を数えることができる。

図

図 12.3 2つの Lorentz の山(価電子:ωv\omega_{\mathrm v}・幅 0.15 ωv0.15\,\omega_{\mathrm v}・重み 1/41/4、内殻:8 ωv8\,\omega_{\mathrm v}・幅 1.5 ωv1.5\,\omega_{\mathrm v}・重み 3/43/4)を持つ模型の有効電子数 (12.18)。この模型の中では厳密な値で、線は点をつないだだけである。ω→∞\omega\to\infty で 11 に近づく

例題 12.4 Drude 金属の3つの和 ★★

オリジナルDrude 模型の比誘電率 εr=1−ωp2/[ω(ω+iγ)]\varepsilon_{\mathrm r}=1-\omega_{\mathrm p}^2/[\omega(\omega+\mathrm{i}\gamma)] を考える(σ0=ε0ωp2/γ\sigma_0=\varepsilon_0\omega_{\mathrm p}^2/\gamma)。

  1. (1) ε2(ω)\varepsilon_2(\omega) を求め、ff 総和則 (12.15) を直接の積分で確かめよ。結果が γ\gamma によらない理由を述べよ。
  2. (2) Im⁡[−1/εr]\operatorname{Im}[-1/\varepsilon_{\mathrm r}] を求め、それが ωp\omega_{\mathrm p} に共鳴を持つ Lorentz 形であることを示せ。式 (12.17) の第1式を確かめよ。
  3. (3) 金属では式 (12.16) の右辺が発散する。極の項を引いた ε1(0)−1=2π∫0∞[ε2−σ0/(ε0ω)]ω−1dω\varepsilon_1(0)-1=\frac2\pi\int_0^\infty[\varepsilon_2-\sigma_0/(\varepsilon_0\omega)]\omega^{-1}\mathrm{d}\omega を計算し、ε1(0)\varepsilon_1(0) と一致することを確かめよ。
解答

(1) 1ω(ω+iγ)=ω−iγω(ω2+γ2)\frac{1}{\omega(\omega+\mathrm{i}\gamma)}=\frac{\omega-\mathrm{i}\gamma}{\omega(\omega^2+\gamma^2)} なので ε1=1−ωp2/(ω2+γ2)\varepsilon_1=1-\omega_{\mathrm p}^2/(\omega^2+\gamma^2)、ε2=ωp2γ/[ω(ω2+γ2)]\varepsilon_2=\omega_{\mathrm p}^2\gamma/[\omega(\omega^2+\gamma^2)]。

∫0∞ωε2 dω=ωp2γ∫0∞dωω2+γ2=ωp2γ⋅π2γ=π2ωp2.\int_0^\infty\omega\varepsilon_2\,\mathrm{d}\omega=\omega_{\mathrm p}^2\gamma\int_0^\infty\frac{\mathrm{d}\omega}{\omega^2+\gamma^2}=\omega_{\mathrm p}^2\gamma\cdot\frac{\pi}{2\gamma}=\frac\pi2\omega_{\mathrm p}^2.

γ\gamma を大きくすると吸収の山は低く広がるが、γ\gamma は山の高さ(∝1/γ\propto1/\gamma)と幅(∝γ\propto\gamma)に逆向きに入るので面積は変わらない。物理的には、散乱が起きるより前の短い時間に電子が J˙=(Ne2/m)E\dot J=(Ne^2/m)E で加速されることが、γ\gamma と無関係だからである。

(2) 1/εr=ω(ω+iγ)/[ω(ω+iγ)−ωp2]1/\varepsilon_{\mathrm r}=\omega(\omega+\mathrm{i}\gamma)/[\omega(\omega+\mathrm{i}\gamma)-\omega_{\mathrm p}^2] なので

−1εr=−1−ωp2ω2−ωp2+iγω=−1+ωp2ωp2−ω2−iγω,Im⁡[−1εr]=ωp2γω(ωp2−ω2)2+γ2ω2.-\frac{1}{\varepsilon_{\mathrm r}}=-1-\frac{\omega_{\mathrm p}^2}{\omega^2-\omega_{\mathrm p}^2+\mathrm{i}\gamma\omega}=-1+\frac{\omega_{\mathrm p}^2}{\omega_{\mathrm p}^2-\omega^2-\mathrm{i}\gamma\omega},\qquad \operatorname{Im}\Bigl[-\frac{1}{\varepsilon_{\mathrm r}}\Bigr]=\frac{\omega_{\mathrm p}^2\gamma\omega}{(\omega_{\mathrm p}^2-\omega^2)^2+\gamma^2\omega^2}.

これは固有角周波数 ωp\omega_{\mathrm p}・プラズマ周波数 ωp\omega_{\mathrm p} の Lorentz 振動子の Im⁡χ\operatorname{Im}\chi と同じ形で、ω=ωp\omega=\omega_{\mathrm p} に山(プラズモン)を持つ。Lorentz 振動子の ff 総和則(例えば (1) と同じ閉曲線の議論)から、∫0∞ωIm⁡[−1/εr]dω=π2ωp2\int_0^\infty\omega\operatorname{Im}[-1/\varepsilon_{\mathrm r}]\mathrm{d}\omega=\frac\pi2\omega_{\mathrm p}^2 である。速い電子が薄膜を通るときに失うエネルギーの分布はこの量に比例する(ハミルトニアンの選び方で固体を考える 第4部第6章の電子エネルギー損失分光)。

(3) ε2−σ0/(ε0ω)=ωp2ω[γω2+γ2−1γ]=−ωp2ωγ(ω2+γ2)\varepsilon_2-\sigma_0/(\varepsilon_0\omega)=\frac{\omega_{\mathrm p}^2}{\omega}\Bigl[\frac{\gamma}{\omega^2+\gamma^2}-\frac1\gamma\Bigr]=-\frac{\omega_{\mathrm p}^2\omega}{\gamma(\omega^2+\gamma^2)} で、

2π∫0∞(−ωp2γ(ω2+γ2))dω=−2π⋅ωp2γ⋅π2γ=−ωp2γ2.\frac2\pi\int_0^\infty\Bigl(-\frac{\omega_{\mathrm p}^2}{\gamma(\omega^2+\gamma^2)}\Bigr)\mathrm{d}\omega=-\frac2\pi\cdot\frac{\omega_{\mathrm p}^2}{\gamma}\cdot\frac{\pi}{2\gamma}=-\frac{\omega_{\mathrm p}^2}{\gamma^2}.

一方 (1) から ε1(0)=1−ωp2/γ2\varepsilon_1(0)=1-\omega_{\mathrm p}^2/\gamma^2 で、一致する。極の項を引けば、金属でも静的な和を使うことができる。

1 屈折率の和則の整理は M. Altarelli, D. L. Dexter, H. M. Nussenzveig, D. Y. Smith, Phys. Rev. B 6, 4502 (1972)。

木下の勉強ステータス

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

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

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