tecktom.com

出てくるもの(1):正の振動数の KK の関係と分散の形

2026年10月15日11分で読めますAI生成教科書『虎の巻』第12章 Kramers–Kronig の関係と総和則
目次
  1. 正の振動数だけで書く
  2. 細い吸収線と Sellmeier の式
  3. 透明な領域では屈折率が必ず増える
  4. 金属:原点の極
  5. 木下の層

AI 生成・木下が未勉強

光学と量子エレクトロニクス 第12章 §12.3 ★★★ 実数性で負の振動数を消すと、測定できる ω>0\omega>0 だけの式になる。吸収の山はその外側の屈折率を持ち上げ、透明な領域では屈折率が必ず振動数とともに増える

正の振動数だけで書く

実験で測れるのは ω>0\omega>0 の量だけである。式 (12.8) の積分を ω′<0\omega'<0 と ω′>0\omega'>0 に分け、負の側で ω′→−ω′\omega'\to-\omega' と置き換えて実数性 (12.5)(虚部は奇、実部は偶)を使う。第1式では Im⁡χ(ω′)[1/(ω′−ω)+1/(ω′+ω)]=2ω′Im⁡χ(ω′)/(ω′2−ω2)\operatorname{Im}\chi(\omega')[1/(\omega'-\omega)+1/(\omega'+\omega)]=2\omega'\operatorname{Im}\chi(\omega')/(\omega'^2-\omega^2)、第2式では Re⁡χ(ω′)[1/(ω′−ω)−1/(ω′+ω)]=2ωRe⁡χ(ω′)/(ω′2−ω2)\operatorname{Re}\chi(\omega')[1/(\omega'-\omega)-1/(\omega'+\omega)]=2\omega\operatorname{Re}\chi(\omega')/(\omega'^2-\omega^2) がまとまり、

Kramers–Kronig の関係(正の振動数の形)

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

を得る。この形は因果律と実数性の両方を使っている。数値で計算するときは、P ⁣∫0∞dω′/(ω′2−ω2)=0\mathrm P\!\int_0^\infty\mathrm{d}\omega'/(\omega'^2-\omega^2)=0 を使って被積分関数の分子から ω=ω′\omega=\omega' での値を引いておくと、特異点の無い普通の積分になる。

図 12.2 は、Lorentz 振動子(γ=0.1 ω0\gamma=0.1\,\omega_0、ωp=1.5 ω0\omega_{\mathrm p}=1.5\,\omega_0)の ε1\varepsilon_1 と ε2\varepsilon_2 を、式 (12.10) で ε2\varepsilon_2 から作った ε1\varepsilon_1 と比べたものである。ただし実際の測定をまねて、ε2\varepsilon_2 を ω≤1.6 ω0\omega\le1.6\,\omega_0 の範囲だけで与え、その外を 00 として積分した。厳密な ε1\varepsilon_1 との差は低い振動数で 0.020.02、ω=2 ω0\omega=2\,\omega_0 で 0.020.02 の程度で、共鳴の形はほぼ完全に再現される。打ち切った先の吸収は小さいが、それを捨てたことで低い振動数の ε1\varepsilon_1 が一様に少し下がった。式 (12.10) の分母が ω′2\omega'^2 の程度なので、高い振動数の吸収は低い振動数の屈折率に正の定数を足す形で入るからである。

図

図 12.2 Lorentz 振動子(γ=0.1 ω0\gamma=0.1\,\omega_0、ωp=1.5 ω0\omega_{\mathrm p}=1.5\,\omega_0)の誘電関数。破線は与えた ε2\varepsilon_2、点線は厳密な ε1\varepsilon_1、実線は ω≤1.6 ω0\omega\le1.6\,\omega_0 の ε2\varepsilon_2 だけから式 (12.10) で作った ε1\varepsilon_1(打ち切りの近似)。共鳴の手前で ε1\varepsilon_1 が増え(正常分散)、山の幅の中でだけ減る(異常分散)

細い吸収線と Sellmeier の式

吸収が角周波数 ωj\omega_j の細い線に集まっている場合を考える。線の重みを振動子強度 fjf_j で表し、ω>0\omega>0 の側でε2(ω)=π2fjωp2ωj δ(ω−ωj)\varepsilon_2(\omega)=\frac\pi2\frac{f_j\omega_{\mathrm p}^2}{\omega_j}\,\delta(\omega-\omega_j) とおく(係数の選び方は 12.4節の総和則で意味が分かる)。式 (12.10) の第1式に入れると

ε1(ω)−1=∑jfj ωp2ωj2−ω2(12.11)\varepsilon_1(\omega)-1=\sum_j\frac{f_j\,\omega_{\mathrm p}^2}{\omega_j^2-\omega^2}\tag{12.11}

となる。減衰の無い Lorentz 振動子の和と同じ形で、波長で書けば光学ガラスの屈折率を表す Sellmeier の式になる。第11章ではこの形を振動子の運動方程式から出したが、ここでは「吸収が ωj\omega_j に集まっている」ことだけから出た。物質の中身が何であっても、吸収線の位置と重みが分かれば、その外側の屈折率が決まる。

透明な領域では屈折率が必ず増える

ω\omega の近くに吸収が無い透明な領域では、式 (12.10) の第1式の被積分関数は ω′=ω\omega'=\omega で特異にならず、ω\omega で微分して積分の中へ入れることができる。∂ω[1/(ω′2−ω2)]=2ω/(ω′2−ω2)2>0\partial_\omega[1/(\omega'^2-\omega^2)]=2\omega/(\omega'^2-\omega^2)^2>0 なので

dε1dω=4ωπ∫0∞ω′ε2(ω′)(ω′2−ω2)2 dω′>0(透明な領域・ε2≥0)(12.12)\frac{\mathrm{d} \varepsilon_1}{\mathrm{d} \omega}=\frac{4\omega}{\pi}\int_0^\infty\frac{\omega'\varepsilon_2(\omega')}{(\omega'^2-\omega^2)^2}\,\mathrm{d}\omega'>0\qquad(\text{透明な領域・}\varepsilon_2\ge0)\tag{12.12}

である。受動的な物質の透明な領域では、ε1\varepsilon_1 は振動数とともに必ず増える。これが正常分散で、ガラスが青い光を赤い光より強く曲げる理由である。ε1\varepsilon_1 が振動数とともに減る異常分散は、ε2\varepsilon_2 がゼロでない所、つまり吸収のある所でしか起きない(図 12.2 の山の幅の中)。第11章で Lorentz 振動子について見た正常分散と異常分散の配置は、この不等式によってどの物質でも同じになる。反転分布のある利得媒質では ε2<0\varepsilon_2<0 の線があり、その外側で ε1\varepsilon_1 が振動数とともに減ることがある。これは第11章の遅い光と速い光の話につながる。

金属:原点の極

金属では直流の伝導度 σ0\sigma_0 のために、ω→0\omega\to0 で εr≈iσ0/(ε0ω)\varepsilon_{\mathrm r}\approx\mathrm{i}\sigma_0/(\varepsilon_0\omega) と発散する。この極は実軸の上にあるので、図 12.1 の積分路で原点も小さな半円で避ける必要がある。極の部分を引いた χ(ω)−iσ0/(ε0ω)\chi(\omega)-\mathrm{i}\sigma_0/(\varepsilon_0\omega) は原点で正則で、これに KK の関係を当てればよい。引いた項は実軸の上で純虚数なので実部は変わらず、虚部だけが σ0/(ε0ω)\sigma_0/(\varepsilon_0\omega) だけ変わる。正の振動数の形に直すと

ε1(ω)−1=2π P ⁣∫0∞ω′ε2(ω′)ω′2−ω2 dω′,ε2(ω)=σ0ε0ω−2ωπ P ⁣∫0∞ε1(ω′)−1ω′2−ω2 dω′(12.13)\varepsilon_1(\omega)-1=\frac2\pi\,\mathrm P\!\int_0^\infty\frac{\omega'\varepsilon_2(\omega')}{\omega'^2-\omega^2}\,\mathrm{d}\omega',\qquad \varepsilon_2(\omega)=\frac{\sigma_0}{\varepsilon_0\omega}-\frac{2\omega}\pi\,\mathrm P\!\int_0^\infty\frac{\varepsilon_1(\omega')-1}{\omega'^2-\omega^2}\,\mathrm{d}\omega'\tag{12.13}

である。第1式が変わらないのは、引いた項 σ0/(ε0ω′)\sigma_0/(\varepsilon_0\omega') を第1式に入れると ω′\omega' と打ち消して定数になり、P ⁣∫0∞dω′/(ω′2−ω2)=0\mathrm P\!\int_0^\infty\mathrm{d}\omega'/(\omega'^2-\omega^2)=0 で消えるからである。第2式の σ0/(ε0ω)\sigma_0/(\varepsilon_0\omega) は、ε1\varepsilon_1 を知っても直流の伝導度が決まらないことを表す。ε1\varepsilon_1 から作れるのは、直流の極を除いた吸収だけである。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)で式 (12.13) を数値で確かめると、どちらの式も成り立つ(検算に入れた)。

例題 12.3 1本の紫外の吸収線で決まるガラスの屈折率 ★★

オリジナル可視光で透明なガラスの吸収が、波長 λ1=0.100 μm\lambda_1=0.100\,\mathrm{\mu m} の紫外の1本の細い線に集まっているとみなす(赤外の吸収は無視する)。波長 589 nm589\,\mathrm{nm} で屈折率が 1.5001.500 だった。

  1. (1) 式 (12.11) を波長 λ=2πc/ω\lambda=2\pi c/\omega で書き直し、n2−1=Bλ2/(λ2−λ12)n^2-1=B\lambda^2/(\lambda^2-\lambda_1^2) の形にせよ。BB は何を表すか。
  2. (2) BB を求め、波長 400 nm400\,\mathrm{nm} と 700 nm700\,\mathrm{nm} の屈折率、および静的な極限 λ→∞\lambda\to\infty の屈折率を求めよ。可視域の屈折率の順序が式 (12.12) と合うことを確かめよ。
  3. (3) 吸収が少しでもある実際のガラスでも、(2) の順序は変わらないと言えるか。式 (12.12) の仮定に照らして述べよ。
解答

(1) 透明な領域では κ=0\kappa=0 で ε1=n2\varepsilon_1=n^2 である。ω=2πc/λ\omega=2\pi c/\lambda、ω1=2πc/λ1\omega_1=2\pi c/\lambda_1 を入れると

n2−1=f1ωp2ω12−ω2=f1ωp2ω12⋅λ2λ2−λ12,B=f1ωp2ω12.n^2-1=\frac{f_1\omega_{\mathrm p}^2}{\omega_1^2-\omega^2}=\frac{f_1\omega_{\mathrm p}^2}{\omega_1^2}\cdot\frac{\lambda^2}{\lambda^2-\lambda_1^2},\qquad B=\frac{f_1\omega_{\mathrm p}^2}{\omega_1^2}.

BB は λ→∞\lambda\to\infty での n2−1n^2-1、つまり紫外の線が静的な比誘電率に寄与する分である。

(2) λ=0.589 μm\lambda=0.589\,\mathrm{\mu m} で λ2/(λ2−λ12)=0.346921/0.336921=1.02968\lambda^2/(\lambda^2-\lambda_1^2)=0.346921/0.336921=1.02968 なので、B=(1.5002−1)/1.02968=1.214B=(1.500^2-1)/1.02968=1.214。400 nm400\,\mathrm{nm} では 0.16/0.15=1.066670.16/0.15=1.06667 で n=1+1.214×1.06667=1.515n=\sqrt{1+1.214\times1.06667}=1.515、700 nm700\,\mathrm{nm} では 0.49/0.48=1.020830.49/0.48=1.02083 で n=1.496n=1.496、λ→∞\lambda\to\infty では n=1+B=1.488n=\sqrt{1+B}=1.488 である。λ\lambda が短い(ω\omega が大きい)ほど屈折率が大きく、式 (12.12) の正常分散と合う。

(3) 式 (12.12) は、注目する振動数の近くで ε2=0\varepsilon_2=0 であることと、どこでも ε2≥0\varepsilon_2\ge0 であることだけを仮定している。吸収線がいくつあっても、形がどうであっても、可視域に吸収が無い限り被積分関数は正なので、可視域で ε1\varepsilon_1 が振動数とともに増えることは変わらない。ただし ε1=n2−κ2\varepsilon_1=n^2-\kappa^2 なので、可視域にわずかな吸収 κ\kappa がある場合に nn の増加まで保証するには、κ\kappa が小さく κ2\kappa^2 の変化が無視できることを別に確かめる必要がある。

木下の勉強ステータス

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

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

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