tecktom.com

当てる道具:光学定数・反射率・総和則

2026年10月15日8分で読めますAI生成教科書『虎の巻』第13章 Drude と金属
目次
  1. 屈折率・表皮の深さ・反射率
  2. 因果律と Kramers–Kronig の関係
  3. 総和則
  4. 木下の層

AI 生成・木下が未勉強

光学と量子エレクトロニクス 第13章 §13.2 ★★ εr\varepsilon_{\mathrm r} から屈折率・消衰係数・垂直入射の反射率を出す。Drude の式は因果律を満たし、γ\gamma によらず総和則 ∫ωε2 dω=πωp2/2\int\omega\varepsilon_2\,\mathrm{d}\omega=\pi\omega_{\mathrm p}^2/2 を満たす

屈折率・表皮の深さ・反射率

第10章(第10章)で定めたとおり、複素屈折率は n~=n+iκ=εr\tilde n=n+\mathrm{i}\kappa=\sqrt{\varepsilon_{\mathrm r}}(κ≥0\kappa\ge0 の枝)で、実部と虚部は

n=∣εr∣+ε12,κ=∣εr∣−ε12(13.6)n=\sqrt{\frac{\lvert \varepsilon_{\mathrm r}\rvert +\varepsilon_1}{2}},\qquad \kappa=\sqrt{\frac{\lvert \varepsilon_{\mathrm r}\rvert -\varepsilon_1}{2}}\tag{13.6}

である。zz 方向に進む平面波は ein~ωz/c=einωz/ce−κωz/c\mathrm{e}^{\mathrm{i}\tilde n\omega z/c}=\mathrm{e}^{\mathrm{i} n\omega z/c}\mathrm{e}^{-\kappa\omega z/c} なので、強度の吸収係数は第10章の αabs=2κω/c\alpha_{\mathrm{abs}}=2\kappa\omega/c、振幅が 1/e1/\mathrm{e} になる深さは

δs=cκω=2αabs(13.7)\delta_{\mathrm s}=\frac{c}{\kappa\omega}=\frac{2}{\alpha_{\mathrm{abs}}}\tag{13.7}

である。金属ではこれを表皮の深さと呼ぶ(変分の δ\delta と区別するため添字 s を付ける)。ε1<0\varepsilon_1<0 で ε2\varepsilon_2 が小さいときは n≪κn\ll\kappa で、波は進まずに減衰する。

真空から垂直に入射した光の反射率は、第14章(第14章)で境界条件から導く Fresnel の式の垂直入射の場合

R=∣n~−1n~+1∣2=(n−1)2+κ2(n+1)2+κ2(13.8)R=\lvert \frac{\tilde n-1}{\tilde n+1}\rvert ^2=\frac{(n-1)^2+\kappa^2}{(n+1)^2+\kappa^2}\tag{13.8}

で与えられる。この章ではこの式を結果として使う。κ≫1\kappa\gg1 なら nn によらず R≈1R\approx1 で、1−R1-R は nn に比例する。金属の反射の強さを決めるのは、κ\kappa が大きいことと nn が小さいことの2つである。

因果律と Kramers–Kronig の関係

第12章(第12章)では、応答関数が上半面で正則で無限遠で 00 に近づくことから Kramers–Kronig の関係を出した。(13.3) の εr−εb\varepsilon_{\mathrm r}-\varepsilon_{\mathrm b} は ω=0\omega=0 と ω=−iγ\omega=-\mathrm{i}\gamma に極を持つ。下半面の極 −iγ-\mathrm{i}\gamma は、電流が時間 τ\tau で減衰することを表し、因果律と矛盾しない。実軸上の極 ω=0\omega=0 は、一定の電場で電流が流れ続けることを表す。この極があるために、ε2\varepsilon_2 から ε1\varepsilon_1 を出す向きの関係

ε1(ω)−εb=2π P ⁣∫0∞ω′ε2(ω′)ω′2−ω2 dω′(13.9)\varepsilon_1(\omega)-\varepsilon_{\mathrm b}=\frac{2}{\pi}\,\mathrm{P}\!\int_0^\infty\frac{\omega'\varepsilon_2(\omega')}{\omega'^2-\omega^2}\,\mathrm{d}\omega'\tag{13.9}

はそのまま成り立つが(分子の ω′\omega' が ε2\varepsilon_2 の 1/ω′1/\omega' を打ち消す)、逆向きの関係には σ0/(ε0ω)\sigma_0/(\varepsilon_0\omega) の項を足す必要がある。

金属の Kramers–Kronig

逆向きの関係は ε2(ω)=−2ωπP ⁣∫0∞ε1(ω′)−εbω′2−ω2dω′+σ0ε0ω\varepsilon_2(\omega)=-\frac{2\omega}{\pi}\mathrm{P}\!\int_0^\infty\frac{\varepsilon_1(\omega')-\varepsilon_{\mathrm b}}{\omega'^2-\omega^2}\mathrm{d}\omega'+\frac{\sigma_0}{\varepsilon_0\omega} になる。最後の項は実軸上の極の留数から出る。σ(ω)\sigma(\omega) 自身は実軸上に極を持たないので、金属では誘電関数の代わりに伝導度について Kramers–Kronig の関係を書くほうが扱いやすい。

総和則

ω\omega が大きいとき (13.4) は ε1→εb−ωp2/ω2\varepsilon_1\to\varepsilon_{\mathrm b}-\omega_{\mathrm p}^2/\omega^2 になる。一方 (13.9) の右辺は、ω\omega が ε2\varepsilon_2 の広がりより十分大きいとき −2πω2∫0∞ω′ε2 dω′-\frac{2}{\pi\omega^2}\int_0^\infty\omega'\varepsilon_2\,\mathrm{d}\omega' に近づく。2つを比べると

∫0∞ω ε2(ω) dω=π2 ωp2=π2nee2ε0m(13.10)\int_0^\infty\omega\,\varepsilon_2(\omega)\,\mathrm{d}\omega=\frac{\pi}{2}\,\omega_{\mathrm p}^2=\frac{\pi}{2}\frac{n_{\mathrm e}e^2}{\varepsilon_0m}\tag{13.10}

を得る。第12章の総和則を伝導電子に当てはめたものである。Drude の式では直接 ∫0∞γωp2 dω/(ω2+γ2)=π2ωp2\int_0^\infty\gamma\omega_{\mathrm p}^2\,\mathrm{d}\omega/(\omega^2+\gamma^2)=\frac{\pi}{2}\omega_{\mathrm p}^2 と確かめられ、結果は γ\gamma によらない。衝突は吸収の山の形と幅を変えるが、山の面積は電子の数だけで決まる。面積が電子の数だけで決まるのは、高い周波数では束縛があっても無くても電子が同じに振る舞うからである(例題 13.1 の (2))。

例題 13.2 Drude の ε2\varepsilon_2 から ε1\varepsilon_1 を作る ★★

オリジナルDrude の ε2\varepsilon_2 (13.4) を (13.9) に代入し、ε1\varepsilon_1 が (13.4) に戻ることを示せ。主値積分 P ⁣∫0∞dω′/(ω′2−ω2)=0\mathrm{P}\!\int_0^\infty\mathrm{d}\omega'/(\omega'^2-\omega^2)=0(ω>0\omega>0)を使ってよい。

解答

被積分関数は ω′ε2(ω′)=γωp2/(ω′2+γ2)\omega'\varepsilon_2(\omega')=\gamma\omega_{\mathrm p}^2/(\omega'^2+\gamma^2) なので

γωp2(ω′2+γ2)(ω′2−ω2)=γωp2ω2+γ2(1ω′2−ω2−1ω′2+γ2)\frac{\gamma\omega_{\mathrm p}^2}{(\omega'^2+\gamma^2)(\omega'^2-\omega^2)} =\frac{\gamma\omega_{\mathrm p}^2}{\omega^2+\gamma^2}\Bigl(\frac{1}{\omega'^2-\omega^2}-\frac{1}{\omega'^2+\gamma^2}\Bigr)

と部分分数に分ける。第1項の主値積分は 00、第2項は ∫0∞dω′/(ω′2+γ2)=π/(2γ)\int_0^\infty\mathrm{d}\omega'/(\omega'^2+\gamma^2)=\pi/(2\gamma) なので

ε1(ω)−εb=2π⋅γωp2ω2+γ2(0−π2γ)=−ωp2ω2+γ2\varepsilon_1(\omega)-\varepsilon_{\mathrm b}=\frac{2}{\pi}\cdot\frac{\gamma\omega_{\mathrm p}^2}{\omega^2+\gamma^2}\Bigl(0-\frac{\pi}{2\gamma}\Bigr)=-\frac{\omega_{\mathrm p}^2}{\omega^2+\gamma^2}

で、(13.4) の ε1\varepsilon_1 に一致する。吸収 ε2\varepsilon_2 を全周波数で測れば、分散 ε1\varepsilon_1 は計算で決まる。ω≫γ\omega\gg\gamma で右辺は −ωp2/ω2-\omega_{\mathrm p}^2/\omega^2 になり、総和則 (13.10) の面積 πωp2/2\pi\omega_{\mathrm p}^2/2 を 2/(πω2)2/(\pi\omega^2) 倍したものになっている。

木下の勉強ステータス

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

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

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