光学と量子エレクトロニクス 第12章 §12.2 ★★★ χ ( t ) = 0 \chi(t)=0 χ ( t ) = 0 (t < 0 t<0 t < 0 )は「χ ( ω ) \chi(\omega) χ ( ω ) が上半面で正則」と同じ内容を持ち、Cauchy の積分から実部と虚部を結ぶ KK の関係が出る
複素数の振動数へ延長する
式 (12.3) の ω \omega ω を複素数 ω = ω R + i ω I \omega=\omega_{\mathrm R}+\mathrm{i}\omega_{\mathrm I} ω = ω R + i ω I に置き換えると、被積分関数は χ ( t ) e i ω R t e − ω I t \chi(t)\mathrm{e}^{\mathrm{i}\omega_{\mathrm R}t}\mathrm{e}^{-\omega_{\mathrm I}t} χ ( t ) e i ω R t e − ω I t になる。積分は t > 0 t>0 t > 0 だけにわたるので、ω I > 0 \omega_{\mathrm I}>0 ω I > 0 (上半面)なら因子 e − ω I t \mathrm{e}^{-\omega_{\mathrm I}t} e − ω I t が t → ∞ t\to\infty t → ∞ で減衰し、χ ( t ) \chi(t) χ ( t ) が有界である限り積分は収束する。ω \omega ω で微分しても t e − ω I t t\,\mathrm{e}^{-\omega_{\mathrm I}t} t e − ω I t が掛かるだけで収束は保たれるので、χ ( ω ) \chi(\omega) χ ( ω ) は上半面で正則(複素微分可能)である。さらに ∣ ω ∣ → ∞ \lvert \omega\rvert \to\infty ∣ ω ∣ → ∞ で上半面のどの向きに進んでも χ ( ω ) → 0 \chi(\omega)\to0 χ ( ω ) → 0 になる。実軸の上では振動の打ち消し合い(Riemann–Lebesgue の補題)により、上半面ではそれに減衰が加わるからである。
もし χ ( t ) \chi(t) χ ( t ) が t < 0 t<0 t < 0 にも値を持てば、積分には t → − ∞ t\to-\infty t → − ∞ の部分が加わり、そこでは e − ω I t \mathrm{e}^{-\omega_{\mathrm I}t} e − ω I t が上半面で発散する。したがって、上半面で正則であることは因果律と結びついている。この対応は逆向きにも成り立つ。
定理 12.1(Titchmarsh の定理 ★★★)
χ ( ω ) \chi(\omega) χ ( ω ) が実軸の上で2乗可積分なとき、次の3つは同値である。
(i) 逆変換した χ ( t ) \chi(t) χ ( t ) が t < 0 t<0 t < 0 で 0 0 0 になる。
(ii) χ ( ω ) \chi(\omega) χ ( ω ) は、上半面で正則な関数 χ ( ω R + i ω I ) \chi(\omega_{\mathrm R}+\mathrm{i}\omega_{\mathrm I}) χ ( ω R + i ω I ) の ω I → 0 + \omega_{\mathrm I}\to0^+ ω I → 0 + の極限で、∫ ∣ χ ( ω R + i ω I ) ∣ 2 d ω R \int\lvert \chi(\omega_{\mathrm R}+\mathrm{i}\omega_{\mathrm I})\rvert ^2\mathrm{d}\omega_{\mathrm R} ∫ ∣ χ ( ω R + i ω I ) ∣ 2 d ω R が ω I > 0 \omega_{\mathrm I}>0 ω I > 0 について有界である。
(iii) χ ( ω ) \chi(\omega) χ ( ω ) の実部と虚部が、次の式 (12.8) で互いに決まる。
証明は E. C. Titchmarsh の Fourier 積分の教科書(1937年)にある1 1 E. C. Titchmarsh, Introduction to the Theory of Fourier Integrals (Clarendon Press, Oxford, 1937)。物理の因果律との関係を整理したのは J. S. Toll, Phys. Rev. 104 , 1760 (1956)。 。この章では (i) から (iii) への向きだけを使い、その道筋を次に示す。誘電体の χ ( ω ) \chi(\omega) χ ( ω ) は Lorentz 振動子の例で見たように ∣ ω ∣ → ∞ \lvert \omega\rvert \to\infty ∣ ω ∣ → ∞ で ω − 2 \omega^{-2} ω − 2 で減るので、2乗可積分の条件は満たされる。金属の χ ( ω ) \chi(\omega) χ ( ω ) は ω = 0 \omega=0 ω = 0 に極を持ち、この条件から外れる。その扱いは 12.3 節で述べる。
Cauchy の積分から KK の関係へ
実軸の上の1点 ω \omega ω をとり、上半面で正則な関数 χ ( ω ′ ) / ( ω ′ − ω ) \chi(\omega')/(\omega'-\omega) χ ( ω ′ ) / ( ω ′ − ω ) を、図 12.1 の閉曲線に沿って ω ′ \omega' ω ′ で積分する。閉曲線は実軸の上を − R -R − R から R R R まで進み、ω ′ = ω \omega'=\omega ω ′ = ω の点だけを半径 δ \delta δ の小さな半円で上へ避け、半径 R R R の大きな半円で上半面を回って戻る。内側に極が無いので、Cauchy の積分定理から積分は 0 0 0 になる。R → ∞ R\to\infty R → ∞ で大きな半円の寄与は χ → 0 \chi\to0 χ → 0 により消える。小さな半円を時計回りに半周する寄与は − i π χ ( ω ) -\mathrm{i}\pi\chi(\omega) − i π χ ( ω ) で、実軸の残りの部分は δ → 0 \delta\to0 δ → 0 で主値積分 P ∫ \mathrm P\!\int P ∫ (ω \omega ω の前後を対称に δ \delta δ だけ切り取ってから δ → 0 \delta\to0 δ → 0 とする積分)になる。よって
P ∫ − ∞ ∞ χ ( ω ′ ) ω ′ − ω d ω ′ − i π χ ( ω ) = 0 , χ ( ω ) = 1 i π P ∫ − ∞ ∞ χ ( ω ′ ) ω ′ − ω d ω ′ (12.7) \mathrm P\!\int_{-\infty}^{\infty}\frac{\chi(\omega')}{\omega'-\omega}\,\mathrm{d}\omega'-\mathrm{i}\pi\chi(\omega)=0,\qquad
\chi(\omega)=\frac{1}{\mathrm{i}\pi}\,\mathrm P\!\int_{-\infty}^{\infty}\frac{\chi(\omega')}{\omega'-\omega}\,\mathrm{d}\omega'\tag{12.7} P ∫ − ∞ ∞ ω ′ − ω χ ( ω ′ ) d ω ′ − i π χ ( ω ) = 0 , χ ( ω ) = i π 1 P ∫ − ∞ ∞ ω ′ − ω χ ( ω ′ ) d ω ′ ( 12.7 )
を得る。右辺の 1 / i 1/\mathrm{i} 1/ i が実部と虚部を入れ替えるので、両辺の実部と虚部をとって
Kramers–Kronig の関係(全実軸の形)
Re χ ( ω ) = 1 π P ∫ − ∞ ∞ Im χ ( ω ′ ) ω ′ − ω d ω ′ , Im χ ( ω ) = − 1 π P ∫ − ∞ ∞ Re χ ( ω ′ ) ω ′ − ω d ω ′ (12.8) \operatorname{Re}\chi(\omega)=\frac1\pi\,\mathrm P\!\int_{-\infty}^{\infty}\frac{\operatorname{Im}\chi(\omega')}{\omega'-\omega}\,\mathrm{d}\omega',\qquad
\operatorname{Im}\chi(\omega)=-\frac1\pi\,\mathrm P\!\int_{-\infty}^{\infty}\frac{\operatorname{Re}\chi(\omega')}{\omega'-\omega}\,\mathrm{d}\omega'\tag{12.8} Re χ ( ω ) = π 1 P ∫ − ∞ ∞ ω ′ − ω Im χ ( ω ′ ) d ω ′ , Im χ ( ω ) = − π 1 P ∫ − ∞ ∞ ω ′ − ω Re χ ( ω ′ ) d ω ′ ( 12.8 )
が出る。虚部(吸収)がすべての振動数で分かれば実部(分散)が決まり、逆も成り立つ。この形が使ったのは因果律と、χ → 0 \chi\to0 χ → 0 (∣ ω ∣ → ∞ \lvert \omega\rvert \to\infty ∣ ω ∣ → ∞ )だけである。関係は H. A. Kramers と R. de L. Kronig が独立に見いだした2 2 R. de L. Kronig, J. Opt. Soc. Am. 12 , 547 (1926);H. A. Kramers, Atti Congr. Int. Fisici (Como) 2 , 545 (1927)。 。数学では、式 (12.8) の右辺の積分変換を Hilbert 変換 と呼ぶ。
図 12.1 KK の関係を出す積分路。χ ( ω ′ ) \chi(\omega') χ ( ω ′ ) は上半面で正則で、極(Lorentz 振動子なら ω ′ = − i γ / 2 ± ω ˉ \omega'=-\mathrm{i}\gamma/2\pm\bar\omega ω ′ = − i γ /2 ± ω ˉ )はすべて下半面にある。実軸の上の点 ω \omega ω を小さな半円で避け、大きな半円で上半面を閉じる
同じ内容は、実軸から少しだけ上へずらした分母の公式
1 ω ′ − ω − i 0 = P 1 ω ′ − ω + i π δ ( ω ′ − ω ) (12.9) \frac{1}{\omega'-\omega-\mathrm{i}0}=\mathrm P\frac{1}{\omega'-\omega}+\mathrm{i}\pi\,\delta(\omega'-\omega)\tag{12.9} ω ′ − ω − i 0 1 = P ω ′ − ω 1 + i π δ ( ω ′ − ω ) ( 12.9 )
(Sokhotski–Plemelj の公式)でも書くことができる。左辺の虚部 η / [ ( ω ′ − ω ) 2 + η 2 ] \eta/[(\omega'-\omega)^2+\eta^2] η / [( ω ′ − ω ) 2 + η 2 ] は面積 π \pi π の山で、η → 0 + \eta\to0^+ η → 0 + でデルタ関数の π \pi π 倍になる。式 (12.7) の小さな半円の寄与 − i π χ ( ω ) -\mathrm{i}\pi\chi(\omega) − i π χ ( ω ) は、この山の面積に当たる。
なぜ因果律が実部と虚部を結ぶのか
階段関数 Θ ( t ) \Theta(t) Θ ( t ) を掛けても χ ( t ) \chi(t) χ ( t ) が変わらない、という式 χ ( t ) = Θ ( t ) χ ( t ) \chi(t)=\Theta(t)\chi(t) χ ( t ) = Θ ( t ) χ ( t ) を Fourier 変換すると、積は畳み込みになる。Θ ( t ) \Theta(t) Θ ( t ) の Fourier 変換は i / ( ω + i 0 ) \mathrm{i}/(\omega+\mathrm{i}0) i / ( ω + i 0 ) で、その実部はデルタ関数、虚部は P ( 1 / ω ) \mathrm P(1/\omega) P ( 1/ ω ) である。畳み込みの結果がもとの χ ( ω ) \chi(\omega) χ ( ω ) に戻るという条件が、実部と虚部を P ( 1 / ω ) \mathrm P(1/\omega) P ( 1/ ω ) との畳み込みで入れ替える式、つまり式 (12.8) になる。時間の半直線に応答が閉じ込められていることが、振動数の側では「実部と虚部が独立でない」ことに対応する。
例題 12.2 Debye 緩和で KK の関係を確かめる ★★
オリジナル 極性分子の液体の向きの緩和を、応答関数 χ ( t ) = ( Δ χ / τ ) e − t / τ Θ ( t ) \chi(t)=(\Delta\chi/\tau)\,\mathrm{e}^{-t/\tau}\Theta(t) χ ( t ) = ( Δ χ / τ ) e − t / τ Θ ( t ) で表す(Debye 緩和。Δ χ > 0 \Delta\chi>0 Δ χ > 0 、τ > 0 \tau>0 τ > 0 )。
(1) χ ( ω ) \chi(\omega) χ ( ω ) を求め、実部と虚部を書け。極はどこにあるか。
(2) 式 (12.8) の第1式の右辺を部分分数で計算し、実部が再現されることを確かめよ。
(3) 実部が同じで虚部の符号だけが逆の関数 χ × ( ω ) = Δ χ / ( 1 + i ω τ ) \chi_{\times}(\omega)=\Delta\chi/(1+\mathrm{i}\omega\tau) χ × ( ω ) = Δ χ / ( 1 + i ω τ ) を考える。対応する χ × ( t ) \chi_{\times}(t) χ × ( t ) を求め、式 (12.8) が成り立たないことを確かめよ。
解答 (1) χ ( ω ) = ( Δ χ / τ ) ∫ 0 ∞ e ( i ω − 1 / τ ) t d t = Δ χ / ( 1 − i ω τ ) \chi(\omega)=(\Delta\chi/\tau)\int_0^\infty\mathrm{e}^{(\mathrm{i}\omega-1/\tau)t}\mathrm{d} t=\Delta\chi/(1-\mathrm{i}\omega\tau) χ ( ω ) = ( Δ χ / τ ) ∫ 0 ∞ e ( i ω − 1/ τ ) t d t = Δ χ / ( 1 − i ω τ ) で、
Re χ = Δ χ 1 + ω 2 τ 2 , Im χ = Δ χ ω τ 1 + ω 2 τ 2 . \operatorname{Re}\chi=\frac{\Delta\chi}{1+\omega^2\tau^2},\qquad \operatorname{Im}\chi=\frac{\Delta\chi\,\omega\tau}{1+\omega^2\tau^2}. Re χ = 1 + ω 2 τ 2 Δ χ , Im χ = 1 + ω 2 τ 2 Δ χ ω τ .
極は ω = − i / τ \omega=-\mathrm{i}/\tau ω = − i / τ で、下半面の虚軸の上に1つある。虚部は ω τ = 1 \omega\tau=1 ω τ = 1 で最大値 Δ χ / 2 \Delta\chi/2 Δ χ /2 をとり、実部はそこで静的な値 Δ χ \Delta\chi Δ χ の半分に下がっている。
(2) ω ′ \omega' ω ′ の関数として
ω ′ τ ( 1 + ω ′ 2 τ 2 ) ( ω ′ − ω ) = a ω ′ − ω + b ω ′ + c 1 + ω ′ 2 τ 2 \frac{\omega'\tau}{(1+\omega'^2\tau^2)(\omega'-\omega)}=\frac{a}{\omega'-\omega}+\frac{b\,\omega'+c}{1+\omega'^2\tau^2} ( 1 + ω ′2 τ 2 ) ( ω ′ − ω ) ω ′ τ = ω ′ − ω a + 1 + ω ′2 τ 2 b ω ′ + c
と分ける。分母を払って ω ′ 2 \omega'^2 ω ′2 ・ω ′ \omega' ω ′ ・定数の係数を比べると a τ 2 + b = 0 a\tau^2+b=0 a τ 2 + b = 0 、c − b ω = τ c-b\omega=\tau c − bω = τ 、a − c ω = 0 a-c\omega=0 a − c ω = 0 で、解はa = ω τ / ( 1 + ω 2 τ 2 ) a=\omega\tau/(1+\omega^2\tau^2) a = ω τ / ( 1 + ω 2 τ 2 ) 、b = − a τ 2 b=-a\tau^2 b = − a τ 2 、c = τ / ( 1 + ω 2 τ 2 ) c=\tau/(1+\omega^2\tau^2) c = τ / ( 1 + ω 2 τ 2 ) である。全実軸の主値積分では、a / ( ω ′ − ω ) a/(\omega'-\omega) a / ( ω ′ − ω ) の項は ω \omega ω の前後で打ち消し合って 0 0 0 、b ω ′ / ( 1 + ω ′ 2 τ 2 ) b\omega'/(1+\omega'^2\tau^2) b ω ′ / ( 1 + ω ′2 τ 2 ) の項は奇関数なので 0 0 0 になり、残る c ∫ d ω ′ / ( 1 + ω ′ 2 τ 2 ) = c π / τ c\int\mathrm{d}\omega'/(1+\omega'^2\tau^2)=c\pi/\tau c ∫ d ω ′ / ( 1 + ω ′2 τ 2 ) = c π / τ だけが寄与する。よって
1 π P ∫ − ∞ ∞ Im χ ( ω ′ ) ω ′ − ω d ω ′ = Δ χ π ⋅ π τ ⋅ τ 1 + ω 2 τ 2 = Δ χ 1 + ω 2 τ 2 = Re χ ( ω ) \frac1\pi\,\mathrm P\!\int_{-\infty}^\infty\frac{\operatorname{Im}\chi(\omega')}{\omega'-\omega}\mathrm{d}\omega'=\frac{\Delta\chi}{\pi}\cdot\frac{\pi}{\tau}\cdot\frac{\tau}{1+\omega^2\tau^2}=\frac{\Delta\chi}{1+\omega^2\tau^2}=\operatorname{Re}\chi(\omega) π 1 P ∫ − ∞ ∞ ω ′ − ω Im χ ( ω ′ ) d ω ′ = π Δ χ ⋅ τ π ⋅ 1 + ω 2 τ 2 τ = 1 + ω 2 τ 2 Δ χ = Re χ ( ω )
で、実部が再現される。
(3) χ × ( ω ) \chi_\times(\omega) χ × ( ω ) は χ ( ω ) \chi(\omega) χ ( ω ) の複素共役で、極は上半面の ω = + i / τ \omega=+\mathrm{i}/\tau ω = + i / τ にある。χ × ( t ) = ( Δ χ / τ ) e t / τ Θ ( − t ) \chi_\times(t)=(\Delta\chi/\tau)\,\mathrm{e}^{t/\tau}\Theta(-t) χ × ( t ) = ( Δ χ / τ ) e t / τ Θ ( − t ) とおくと∫ − ∞ 0 ( Δ χ / τ ) e ( i ω + 1 / τ ) t d t = Δ χ / ( 1 + i ω τ ) \int_{-\infty}^0(\Delta\chi/\tau)\mathrm{e}^{(\mathrm{i}\omega+1/\tau)t}\mathrm{d} t=\Delta\chi/(1+\mathrm{i}\omega\tau) ∫ − ∞ 0 ( Δ χ / τ ) e ( i ω + 1/ τ ) t d t = Δ χ / ( 1 + i ω τ ) になるので、これが対応する応答関数である。この応答は t < 0 t<0 t < 0 にだけ値を持ち、電場が来る前に分極が立ち上がる。(2) と同じ計算で虚部の符号だけが変わるので、式 (12.8) の右辺は − Re χ × ( ω ) -\operatorname{Re}\chi_\times(\omega) − Re χ × ( ω ) になり、左辺と符号が合わない。実部の形だけを見ても応答が因果的かどうかは分からず、実部と虚部の組で初めて判定できる。
1 E. C. Titchmarsh, Introduction to the Theory of Fourier Integrals (Clarendon Press, Oxford, 1937)。物理の因果律との関係を整理したのは J. S. Toll, Phys. Rev. 104 , 1760 (1956)。
2 R. de L. Kronig, J. Opt. Soc. Am. 12 , 547 (1926);H. A. Kramers, Atti Congr. Int. Fisici (Como) 2 , 545 (1927)。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com