tecktom.com

当てる道具:線形応答と感受率

2026年10月15日11分で読めますAI生成教科書『虎の巻』第10章 物質中の Maxwell 方程式と光学定数
目次
  1. 応答関数と感受率
  2. 感受率の3つの性質
  3. 自由電荷の電流を誘電関数に入れる
  4. 木下の層

AI 生成・木下が未勉強

光学と量子エレクトロニクス 第10章 §10.2 ★★★ 線形・時間並進不変・因果的・局所的な応答は1つの関数 χ(ω)\chi(\omega) で書け、χ\chi は実数性・上半面での解析性・受動性の3つを満たす

応答関数と感受率

物質の応答に次の4つを仮定する。線形:電場を2倍にすれば分極も2倍になる。時間並進不変:物質の性質が時間によらない。因果的:時刻 tt の分極は時刻 tt より前の電場だけで決まる。局所的:点 r\bm{r} の分極は同じ点 r\bm{r} の電場だけで決まる。さらに媒質が等方的であるとすると(異方性は第15章)、最も一般の形は

P(r,t)=ε0∫0∞χ(τ) E(r,t−τ) dτ(10.6)\bm{P}(\bm{r},t)=\varepsilon_0\int_0^\infty\chi(\tau)\,\bm{E}(\bm{r},t-\tau)\,\mathrm{d}\tau\tag{10.6}

である。χ(τ)\chi(\tau) は応答関数で、時刻 00 に一瞬だけ加えた電場に対する分極の時間変化を表す。積分を τ≥0\tau\ge0 に限ったことが因果性で、τ<0\tau<0 では χ(τ)=0\chi(\tau)=0 と言っても同じである。電場も分極も実数の量なので、χ(τ)\chi(\tau) は実数である。

時間の Fourier 変換(記法の正本の流儀:f(t)=∫dω2πf~(ω)e−iωtf(t)=\int\frac{\mathrm{d}\omega}{2\pi}\tilde f(\omega)\mathrm{e}^{-\mathrm{i}\omega t})を取ると、畳み込みは積になる:

P~(r,ω)=ε0χ(ω)E~(r,ω),χ(ω)=∫0∞χ(τ) eiωτ dτ.(10.7)\tilde{\bm{P}}(\bm{r},\omega)=\varepsilon_0\chi(\omega)\tilde{\bm{E}}(\bm{r},\omega),\qquad \chi(\omega)=\int_0^\infty\chi(\tau)\,\mathrm{e}^{\mathrm{i}\omega\tau}\,\mathrm{d}\tau .\tag{10.7}

この χ(ω)\chi(\omega) を感受率と呼ぶ。応答関数 χ(τ)\chi(\tau) と感受率 χ(ω)\chi(\omega) は同じ字で書き、引数で区別する。角周波数 ω\omega の単色の電場 E=Re⁡[E0e−iωt]\bm{E}=\operatorname{Re}[\bm{E}_0\mathrm{e}^{-\mathrm{i}\omega t}] に対しては、式 (10.6) に入れて積分すれば

P=Re⁡[ε0χ(ω)E0e−iωt](10.8)\bm{P}=\operatorname{Re}[\varepsilon_0\chi(\omega)\bm{E}_0\mathrm{e}^{-\mathrm{i}\omega t}]\tag{10.8}

となる。χ(ω)\chi(\omega) の絶対値が分極の振幅を、偏角が電場に対する分極の位相の遅れを表す。第1章の分極率 α(ω)\alpha(\omega) は1原子の双極子 d=αE\bm{d}=\alpha\bm{E} の係数で、原子が希薄で互いの場を感じないなら P=NαE\bm{P}=N\alpha\bm{E}、つまり χ(ω)=Nα(ω)/ε0\chi(\omega)=N\alpha(\omega)/\varepsilon_0 である。密な媒質で原子が感じる場が E\bm{E} からずれる補正(局所場)は第11章(第11章)で扱う。

感受率の3つの性質

命題 10.1(感受率の性質 ★★★)

式 (10.6) の χ(τ)\chi(\tau) が実数で、τ≥0\tau\ge0 で積分可能なとき、感受率 χ(ω)\chi(\omega) は次の3つを満たす。

  1. (i) 実数性:実数の ω\omega について χ(−ω)=χ(ω)∗\chi(-\omega)=\chi(\omega)^*。したがって Re⁡χ\operatorname{Re}\chi は ω\omega の偶関数、Im⁡χ\operatorname{Im}\chi は奇関数である。
  2. (ii) 解析性:χ(ω)\chi(\omega) は複素数の ω\omega に拡張でき、上半面 Im⁡ω>0\operatorname{Im}\omega>0 で正則である。
  3. (iii) 受動性:物質がエネルギーを外へ出さない(熱平衡にある)なら、ω>0\omega>0 で Im⁡χ(ω)≥0\operatorname{Im}\chi(\omega)\ge0 である。
証明

(i) χ(τ)\chi(\tau) が実数なので χ(−ω)=∫0∞χ(τ)e−iωτdτ=(∫0∞χ(τ)eiωτdτ)∗\chi(-\omega)=\int_0^\infty\chi(\tau)\mathrm{e}^{-\mathrm{i}\omega\tau}\mathrm{d}\tau=\bigl(\int_0^\infty\chi(\tau)\mathrm{e}^{\mathrm{i}\omega\tau}\mathrm{d}\tau\bigr)^*。

(ii) ω=ω1+iω2\omega=\omega_1+\mathrm{i}\omega_2(ω2>0\omega_2>0)とすると ∣eiωτ∣=e−ω2τ≤1\lvert \mathrm{e}^{\mathrm{i}\omega\tau}\rvert =\mathrm{e}^{-\omega_2\tau}\le1 が τ≥0\tau\ge0 で成り立つので、積分 (10.7) は収束し、ω\omega について微分しても収束する。積分範囲が τ≥0\tau\ge0 に限られていること、つまり因果性が、上半面での収束を保証している。

(iii) 単色の電場が媒質にする仕事の時間平均を計算する。単位体積あたりの仕事率は E⋅∂tP\bm{E}\cdot\partial_t\bm{P} で、複素表示の時間平均の公式 ⟨Re⁡[Ae−iωt]Re⁡[Be−iωt]⟩=12Re⁡(AB∗)\langle \operatorname{Re}[A\mathrm{e}^{-\mathrm{i}\omega t}]\operatorname{Re}[B\mathrm{e}^{-\mathrm{i}\omega t}] \rangle =\frac12\operatorname{Re}(AB^*) をA=−iωε0χE0A=-\mathrm{i}\omega\varepsilon_0\chi\bm{E}_0, B=E0B=\bm{E}_0 に使うと

⟨E⋅∂P∂t⟩=12Re⁡(−iωε0χ(ω))∣E0∣2=12ωε0Im⁡χ(ω) ∣E0∣2(10.9)\langle \bm{E}\cdot\frac{\partial \bm{P}}{\partial t} \rangle =\frac12\operatorname{Re}\bigl(-\mathrm{i}\omega\varepsilon_0\chi(\omega)\bigr)\lvert \bm{E}_0\rvert ^2=\frac12\omega\varepsilon_0\operatorname{Im}\chi(\omega)\,\lvert \bm{E}_0\rvert ^2\tag{10.9}

となる。熱平衡にある物質は、周期的な外力からエネルギーを取り出されることはない(取り出せるなら第二種の永久機関になる)ので、この量は 00 以上である。ω>0\omega>0 なら Im⁡χ≥0\operatorname{Im}\chi\ge0。

(ii) は第12章(第12章)で Kramers–Kronig の関係を出す出発点になる。Re⁡χ\operatorname{Re}\chi と Im⁡χ\operatorname{Im}\chi は独立でなく、一方を全周波数で知れば他方が決まる。(iii) は、反転分布を作った物質(第21章)では成り立たない。そのとき Im⁡χ<0\operatorname{Im}\chi<0 になり、媒質は光を増幅する。レーザーの利得は、受動性を外部からの励起で破った状態である。

自由電荷の電流を誘電関数に入れる

金属のように自由に動ける電荷があるとき、その電流にも同じ形の線形応答を仮定し、J~f(ω)=σ(ω)E~(ω)\tilde{\bm{J}}_{\mathrm f}(\omega)=\sigma(\omega)\tilde{\bm{E}}(\omega) と書く(σ(ω)\sigma(\omega) は複素数の伝導度)。単色の場では Ampère–Maxwell の式の右辺は

Jf+∂D∂t → σE0−iωε0(1+χ)E0=−iωε0(1+χ(ω)+iσ(ω)ε0ω)E0(10.10)\bm{J}_{\mathrm f}+\frac{\partial \bm{D}}{\partial t}\ \to\ \sigma\bm{E}_0-\mathrm{i}\omega\varepsilon_0(1+\chi)\bm{E}_0=-\mathrm{i}\omega\varepsilon_0\Bigl(1+\chi(\omega)+\frac{\mathrm{i}\sigma(\omega)}{\varepsilon_0\omega}\Bigr)\bm{E}_0\tag{10.10}

となり、伝導電流は感受率に iσ/(ε0ω)\mathrm{i}\sigma/(\varepsilon_0\omega) を足したのと区別がつかない。そこで束縛電荷と自由電荷をまとめた誘電関数

εr(ω)=1+χ(ω)+iσ(ω)ε0ω(10.11)\varepsilon_{\mathrm r}(\omega)=1+\chi(\omega)+\frac{\mathrm{i}\sigma(\omega)}{\varepsilon_0\omega}\tag{10.11}

を定義し、εr=ε1+iε2\varepsilon_{\mathrm r}=\varepsilon_1+\mathrm{i}\varepsilon_2 と実部と虚部に分けて書く(第12・13章とハミルトニアンの選び方で固体を考える 第4部第6章と同じ書き方)。ω≠0\omega\neq0 では、ある電荷を「束縛」と呼ぶか「自由」と呼ぶかは帳簿の付け方の違いにすぎない。違いが見えるのは ω→0\omega\to0 で、自由電荷があると ε2∝1/ω\varepsilon_2\propto1/\omega と発散する(直流の電流が流れ続ける)。以下では、自由電荷がある場合も εr(ω)\varepsilon_{\mathrm r}(\omega) に含め、ρf=0\rho_{\mathrm f}=0, Jf=0\bm{J}_{\mathrm f}=0 と書き直した方程式を使う。命題 10.1 の3つの性質は εr−1\varepsilon_{\mathrm r}-1 についても成り立つ(σ(τ)\sigma(\tau) も実数で因果的なため)。ただし ω=0\omega=0 の極は除く。

例題 10.2 Debye の緩和で3つの性質を確かめる ★★

オリジナル分子の永久双極子が電場の向きにそろう速さが有限で、応答関数が χ(τ)=χ0τDe−τ/τD\chi(\tau)=\frac{\chi_0}{\tau_{\mathrm D}}\mathrm{e}^{-\tau/\tau_{\mathrm D}}(τ≥0\tau\ge0、χ0>0\chi_0>0)で与えられる媒質を考える(Debye の緩和)。

  1. (1) 感受率 χ(ω)\chi(\omega) を求めよ。
  2. (2) 命題 10.1 の (i)(ii)(iii) を確かめよ。
  3. (3) 時刻 00 に大きさ EsE_{\mathrm s} の一定の電場を入れたとき(t<0t<0 では 00)、分極 P(t)P(t) を求めよ。
解答

(1) χ(ω)=χ0τD∫0∞e−τ/τD+iωτdτ=χ0τD⋅11/τD−iω=χ01−iωτD\chi(\omega)=\frac{\chi_0}{\tau_{\mathrm D}}\int_0^\infty\mathrm{e}^{-\tau/\tau_{\mathrm D}+\mathrm{i}\omega\tau}\mathrm{d}\tau=\frac{\chi_0}{\tau_{\mathrm D}}\cdot\frac{1}{1/\tau_{\mathrm D}-\mathrm{i}\omega}=\dfrac{\chi_0}{1-\mathrm{i}\omega\tau_{\mathrm D}}。

(2) (i) χ(−ω)=χ0/(1+iωτD)=χ(ω)∗\chi(-\omega)=\chi_0/(1+\mathrm{i}\omega\tau_{\mathrm D})=\chi(\omega)^*。(ii) χ(ω)\chi(\omega) の極は 1−iωτD=01-\mathrm{i}\omega\tau_{\mathrm D}=0、つまり ω=−i/τD\omega=-\mathrm{i}/\tau_{\mathrm D} で、下半面にある。上半面では正則である。(iii) Im⁡χ=χ0ωτD/(1+ω2τD2)\operatorname{Im}\chi=\chi_0\omega\tau_{\mathrm D}/(1+\omega^2\tau_{\mathrm D}^2) は ω>0\omega>0 で正。緩和の時間スケール ωτD=1\omega\tau_{\mathrm D}=1 で最大になり、ω→0\omega\to0(双極子が追いつく)と ω→∞\omega\to\infty(全く追いつかない)では 00 に近づく。分極が電場に遅れて追いつくときにエネルギーが失われる。

(3) 式 (10.6) で E(t−τ)=EsE(t-\tau)=E_{\mathrm s}(τ<t\tau<t)、00(τ>t\tau>t)なので

P(t)=ε0Es∫0tχ0τDe−τ/τDdτ=ε0χ0Es(1−e−t/τD)(t≥0)P(t)=\varepsilon_0E_{\mathrm s}\int_0^t\frac{\chi_0}{\tau_{\mathrm D}}\mathrm{e}^{-\tau/\tau_{\mathrm D}}\mathrm{d}\tau=\varepsilon_0\chi_0E_{\mathrm s}\bigl(1-\mathrm{e}^{-t/\tau_{\mathrm D}}\bigr)\qquad(t\ge0)

で、t<0t<0 では 00 である。分極は電場を入れる前には現れず(因果性)、時間 τD\tau_{\mathrm D} をかけて静的な値 ε0χ0Es=ε0χ(0)Es\varepsilon_0\chi_0E_{\mathrm s}=\varepsilon_0\chi(0)E_{\mathrm s} に近づく。極が下半面の −i/τD-\mathrm{i}/\tau_{\mathrm D} にあることと、時間の応答が e−t/τD\mathrm{e}^{-t/\tau_{\mathrm D}} で減衰することは、同じ1つの事実の2つの表現である。

ほかの巻では

木下の勉強ステータス

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

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

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