tecktom.com

この章のまとめ

2026年10月15日12分で読めますAI生成教科書『虎の巻』第12章 Kramers–Kronig の関係と総和則

AI 生成・木下が未勉強

光学と量子エレクトロニクス 第12章 §12.7 ★★★ 線形・因果的な応答は上半面で正則で、実部と虚部が KK の関係で結ばれる。正常分散・総和則・反射率から位相を出す方法はすべてその帰結

この章では物質のモデルを選ばず、分極が過去の電場にだけ線形に応答するという性質だけを仮定した。応答関数が t<0t<0 で 00 になることは、感受率 χ(ω)\chi(\omega) が上半面で正則であることと同じ内容を持ち、Cauchy の積分から実部と虚部が Hilbert 変換で結ばれる(KK の関係)。実数性を使うと正の振動数だけの形になり、細い吸収線からは Sellmeier の式が、透明な領域では正常分散 dε1/dω>0\mathrm{d}\varepsilon_1/\mathrm{d}\omega>0 が出た。金属では原点の極を引けば同じ形が使え、虚部の側に直流の伝導度の項が残る。高い振動数では電子が自由な粒子として振る舞い χ→−ωp2/ω2\chi\to-\omega_{\mathrm p}^2/\omega^2 になるので、吸収を ω\omega で重み付けた面積が電子の数密度で決まる(ff 総和則)。同じ議論は 1/εr1/\varepsilon_{\mathrm r} と n~\tilde n にも当てはまり、屈折率は振動数について平均すると 11 になる。実験では反射率だけを広く測り、ln⁡R\ln R の KK の関係で位相を出して nn と κ\kappa を組み立てる。この操作で近似が入るのは測っていない振動数の外挿で、その減り方の指数が結果を大きく左右する。KK の関係そのものは近似を含まず、置いた仮定は線形性・局所性・時間的な一様性の3つである。次の第13章の Drude 模型と、第4部の2準位系の線形応答は、どちらもこの章の枠に入る具体例になる。

この章の道具

  • 応答関数:P(t)=ε0∫χ(t−t′)E(t′)dt′P(t)=\varepsilon_0\int\chi(t-t')E(t')\mathrm{d} t'、χ(t)=0\chi(t)=0(t<0t<0)、χ(ω)=∫0∞χ(t)eiωtdt\chi(\omega)=\int_0^\infty\chi(t)\mathrm{e}^{\mathrm{i}\omega t}\mathrm{d} t、χ(−ω)=χ(ω)∗\chi(-\omega)=\chi(\omega)^*、吸収の仕事率 ω2ε0ε2∣E0∣2\frac\omega2\varepsilon_0\varepsilon_2\lvert E_0\rvert ^2。
  • 全実軸の KK:Re⁡χ=1πP ⁣∫Im⁡χ(ω′)ω′−ωdω′\operatorname{Re}\chi=\frac1\pi\mathrm P\!\int\frac{\operatorname{Im}\chi(\omega')}{\omega'-\omega}\mathrm{d}\omega'、Im⁡χ=−1πP ⁣∫Re⁡χ(ω′)ω′−ωdω′\operatorname{Im}\chi=-\frac1\pi\mathrm P\!\int\frac{\operatorname{Re}\chi(\omega')}{\omega'-\omega}\mathrm{d}\omega'(因果律だけ)。
  • 正の振動数の KK:ε1−1=2πP ⁣∫0∞ω′ε2ω′2−ω2dω′\varepsilon_1-1=\frac2\pi\mathrm P\!\int_0^\infty\frac{\omega'\varepsilon_2}{\omega'^2-\omega^2}\mathrm{d}\omega'、ε2=−2ωπP ⁣∫0∞ε1−1ω′2−ω2dω′\varepsilon_2=-\frac{2\omega}\pi\mathrm P\!\int_0^\infty\frac{\varepsilon_1-1}{\omega'^2-\omega^2}\mathrm{d}\omega'(+金属では σ0/(ε0ω)\sigma_0/(\varepsilon_0\omega))。
  • 細い線:ε1−1=∑jfjωp2/(ωj2−ω2)\varepsilon_1-1=\sum_jf_j\omega_{\mathrm p}^2/(\omega_j^2-\omega^2)。透明な領域で dε1/dω>0\mathrm{d}\varepsilon_1/\mathrm{d}\omega>0。
  • 総和則:∫0∞ωε2dω=π2ωp2\int_0^\infty\omega\varepsilon_2\mathrm{d}\omega=\frac\pi2\omega_{\mathrm p}^2、ε1(0)−1=2π∫ε2/ω dω\varepsilon_1(0)-1=\frac2\pi\int\varepsilon_2/\omega\,\mathrm{d}\omega、∫ωIm⁡(−1/εr)dω=π2ωp2\int\omega\operatorname{Im}(-1/\varepsilon_{\mathrm r})\mathrm{d}\omega=\frac\pi2\omega_{\mathrm p}^2、∫(n−1)dω=0\int(n-1)\mathrm{d}\omega=0。
  • 反射の位相:θ=−ωπ∫0∞ln⁡R(ω′)−ln⁡R(ω)ω′2−ω2dω′\theta=-\frac\omega\pi\int_0^\infty\frac{\ln R(\omega')-\ln R(\omega)}{\omega'^2-\omega^2}\mathrm{d}\omega'、n=1−R1+R−2Rcos⁡θn=\frac{1-R}{1+R-2\sqrt R\cos\theta}、κ=2Rsin⁡θ1+R−2Rcos⁡θ\kappa=\frac{2\sqrt R\sin\theta}{1+R-2\sqrt R\cos\theta}。

総合問題 回転波近似の感受率を因果律と実数性で直す ★★★

オリジナル数密度 NN の2準位原子の気体(遷移角周波数 ω0\omega_0、双極子の行列要素 dd、線幅 γ≪ω0\gamma\ll\omega_0)が、弱い光に応答する。第1章で、2準位系の分極率は α2(ω)=2ω0d2/[ℏ(ω02−ω2)]\alpha_2(\omega)=2\omega_0d^2/[\hbar(\omega_0^2-\omega^2)] で、回転波近似ではそのうち共鳴の項 d2/[ℏ(ω0−ω)]d^2/[\hbar(\omega_0-\omega)] だけを残すことを見た。減衰を入れた回転波近似の感受率(局所場は無視する)を

χRWA(ω)=Aω0−ω−iγ/2,A=Nd2ε0ℏ\chi_{\mathrm{RWA}}(\omega)=\frac{A}{\omega_0-\omega-\mathrm{i}\gamma/2},\qquad A=\frac{Nd^2}{\varepsilon_0\hbar}

とする。

  1. (1) (使う道具:12.1節・12.2節)χRWA(ω)\chi_{\mathrm{RWA}}(\omega) の極の位置を述べ、対応する応答関数 χRWA(t)\chi_{\mathrm{RWA}}(t) を求めよ。それが因果律を満たすが実数性を満たさないことを示せ。ω>0\omega>0 で Im⁡χRWA>0\operatorname{Im}\chi_{\mathrm{RWA}}>0 を確かめよ。
  2. (2) (使う道具:12.2節)χRWA\chi_{\mathrm{RWA}} が全実軸の KK の関係 (12.8) を満たすことを示せ。
  3. (3) (使う道具:12.3節)測定にかかるのは ω>0\omega>0 の吸収だけである。ε2(ω)=Im⁡χRWA(ω)\varepsilon_2(\omega)=\operatorname{Im}\chi_{\mathrm{RWA}}(\omega)(ω>0\omega>0)を正の振動数の KK (12.10) に入れ、γ→0\gamma\to0 の細い線の極限で ε1−1\varepsilon_1-1 を求めよ。結果を Nα2/ε0N\alpha_2/\varepsilon_0 と比べ、回転波近似で落とした項がどこから戻ったかを述べよ。
  4. (4) (使う道具:12.4節・第1章の振動子強度)この線の重み ∫0∞ωε2dω\int_0^\infty\omega\varepsilon_2\mathrm{d}\omega を ff 総和則と比べて、振動子強度 ff を dd で表せ。f=0.73f=0.73、λ0=780 nm\lambda_0=780\,\mathrm{nm} のとき dd を ea0ea_0 を単位にして求めよ。
  5. (5) (使う道具:(3)・12.4節の静的な和・第1章の総合問題)ω<ω0\omega<\omega_0 の光に対する分散 ε1−1\varepsilon_1-1 の、正しい値と回転波近似の値の比を求めよ。λ0=780 nm\lambda_0=780\,\mathrm{nm} の原子に 1064 nm1064\,\mathrm{nm} の光を当てたときの比と、静的な極限の比を求めよ。比が 1.11.1(回転波近似の誤差が 10%10\%)になる波長を求めよ。
解答

(1) 分母が 00 になるのは ω=ω0−iγ/2\omega=\omega_0-\mathrm{i}\gamma/2 で、極は下半面に1つだけある。逆変換 χ(t)=∫dω2πχ(ω)e−iωt\chi(t)=\int\frac{\mathrm{d}\omega}{2\pi}\chi(\omega)\mathrm{e}^{-\mathrm{i}\omega t} は、t<0t<0 では上半面で閉じると内側に極が無いので 00、t>0t>0 では下半面で時計回りに閉じて留数を拾い

χRWA(t)=iA e−iω0t−γt/2 Θ(t)\chi_{\mathrm{RWA}}(t)=\mathrm{i} A\,\mathrm{e}^{-\mathrm{i}\omega_0t-\gamma t/2}\,\Theta(t)

を得る。確かめに、∫0∞iAei(ω−ω0)t−γt/2dt=iA/[γ/2−i(ω−ω0)]=A/(ω0−ω−iγ/2)\int_0^\infty\mathrm{i} A\mathrm{e}^{\mathrm{i}(\omega-\omega_0)t-\gamma t/2}\mathrm{d} t=\mathrm{i} A/[\gamma/2-\mathrm{i}(\omega-\omega_0)]=A/(\omega_0-\omega-\mathrm{i}\gamma/2) である。t<0t<0 で 00 なので因果律は満たすが、χRWA(t)\chi_{\mathrm{RWA}}(t) は複素数で、実数の電場に複素数の分極を返す。振動数の側では χRWA(−ω)=A/(ω0+ω−iγ/2)\chi_{\mathrm{RWA}}(-\omega)=A/(\omega_0+\omega-\mathrm{i}\gamma/2) が χRWA(ω)∗=A/(ω0−ω+iγ/2)\chi_{\mathrm{RWA}}(\omega)^*=A/(\omega_0-\omega+\mathrm{i}\gamma/2) と一致せず、実数性 (12.5) が破れる。虚部は Im⁡χRWA=A(γ/2)/[(ω0−ω)2+γ2/4]>0\operatorname{Im}\chi_{\mathrm{RWA}}=A(\gamma/2)/[(\omega_0-\omega)^2+\gamma^2/4]>0 で、ω0\omega_0 を中心とする幅 γ\gamma の Lorentz 形の吸収線である。

(2) 式 (12.8) を出すのに使ったのは、上半面で正則であることと ∣ω∣→∞\lvert \omega\rvert \to\infty で 00 に減ることだけである。χRWA\chi_{\mathrm{RWA}} は極が下半面にしか無く、∣ω∣→∞\lvert \omega\rvert \to\infty で A/ωA/\omega のように減るので、どちらも満たし、式 (12.8) は成り立つ。直接確かめるには、例題 12.2 と同じく ω′\omega' について部分分数に分ければよい(χRWA\chi_{\mathrm{RWA}} は Debye 形の ω\omega を ω−ω0\omega-\omega_0 に、τ\tau を 2/γ2/\gamma に置き換え、i\mathrm{i} を掛けたものなので、同じ計算がそのまま使える)。

(3) γ→0\gamma\to0 で Im⁡χRWA→πA δ(ω−ω0)\operatorname{Im}\chi_{\mathrm{RWA}}\to\pi A\,\delta(\omega-\omega_0) である(面積 π\pi の山)。式 (12.10) の第1式に入れると

ε1(ω)−1=2π⋅ω0⋅πAω02−ω2=2Aω0ω02−ω2=A[1ω0−ω+1ω0+ω].\varepsilon_1(\omega)-1=\frac2\pi\cdot\frac{\omega_0\cdot\pi A}{\omega_0^2-\omega^2}=\frac{2A\omega_0}{\omega_0^2-\omega^2}=A\Bigl[\frac{1}{\omega_0-\omega}+\frac{1}{\omega_0+\omega}\Bigr].

A=Nd2/(ε0ℏ)A=Nd^2/(\varepsilon_0\hbar) を入れると Nα2(ω)/ε0N\alpha_2(\omega)/\varepsilon_0 そのもので、回転波近似で落とした反共鳴の項 1/(ω0+ω)1/(\omega_0+\omega) が戻っている。正の振動数の形は、ω>0\omega>0 の吸収線を実数性で ω<0\omega<0 へ折り返した線(−ω0-\omega_0 にある鏡像)を含んで KK の関係を当てたものになる。反共鳴の項は、その鏡像の線が作る分散である。回転波近似は因果律を保つが実数性を壊す近似で、実数性を課し直すと落とした項が戻る。

(4) 細い線の重みは ∫0∞ωε2dω=ω0⋅πA=πω0Nd2/(ε0ℏ)\int_0^\infty\omega\varepsilon_2\mathrm{d}\omega=\omega_0\cdot\pi A=\pi\omega_0Nd^2/(\varepsilon_0\hbar) である。ff 総和則の1本の線の重み π2f Ne2/(ε0m)\frac\pi2f\,Ne^2/(\varepsilon_0m) と等しいとおくと

f=2mω0d2ℏe2f=\frac{2m\omega_0d^2}{\hbar e^2}

で、第1章の振動子強度の式と一致する。ω0=2πc/(780 nm)=2.415×1015 s−1\omega_0=2\pi c/(780\,\mathrm{nm})=2.415\times10^{15}\,\mathrm{s^{-1}} のとき d2=fℏe2/(2mω0)d^2=f\hbar e^2/(2m\omega_0) からd=0.73×1.0546×10−34/(2×9.109×10−31×2.415×1015)  e=1.32×10−10 m×e=2.50 ea0=2.12×10−29 C md=\sqrt{0.73\times1.0546\times10^{-34}/(2\times9.109\times10^{-31}\times2.415\times10^{15})}\;e=1.32\times10^{-10}\,\mathrm m\times e=2.50\,ea_0=2.12\times10^{-29}\,\mathrm{C\,m} になる。第1章で使った値と同じである。

(5) 正しい値は (3) の 2Aω0/(ω02−ω2)2A\omega_0/(\omega_0^2-\omega^2)、回転波近似の値は Re⁡χRWA≈A/(ω0−ω)\operatorname{Re}\chi_{\mathrm{RWA}}\approx A/(\omega_0-\omega)(線から十分離れた所)なので、比は

2Aω0/(ω02−ω2)A/(ω0−ω)=2ω0ω0+ω\frac{2A\omega_0/(\omega_0^2-\omega^2)}{A/(\omega_0-\omega)}=\frac{2\omega_0}{\omega_0+\omega}

である。ω/ω0=780/1064=0.733\omega/\omega_0=780/1064=0.733 では 2/1.733=1.1542/1.733=1.154 で、回転波近似は分散を 15%15\% 小さく見積もる。第1章の総合問題で光トラップの深さについて出した比 1.1541.154 と同じ数で、双極子力のポテンシャルが分散(分極率の実部)で決まるからである。静的な極限 ω→0\omega\to0 では比は 22 で、回転波近似は静的な分極率をちょうど半分にする。(2/π)∫ε2/ω dω(2/\pi)\int\varepsilon_2/\omega\,\mathrm{d}\omega の静的な和で計算しても 2A/ω02A/\omega_0 になり、正しい値の側と一致する。比が 1.11.1 になるのは ω/ω0=2/1.1−1=0.818\omega/\omega_0=2/1.1-1=0.818、波長で 780/0.818=953 nm780/0.818=953\,\mathrm{nm} で、これより長い波長では回転波近似の誤差が 10%10\% を超える。

木下の勉強ステータス

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

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

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