クラマース・クローニッヒの関係式

クラマース・クローニッヒの関係式(—かんけいしき、Kramers-Kronig relation)とは線形応答における周波数応答関数の実部と虚部の間に成り立つ関係式のことである。 1926年にラルフ・クローニッヒ、1927年にヘンドリック・アントニー・クラマースによって電磁波の分散現象に対して導かれた。

クラマース・クローニッヒの関係式

周波数応答関数H(ω)=H_R(ω)+i H_I(ω)に対して(ただし、H_R はHの実部、H_I はHの虚部である。)

$H_R(\omega) = \frac{1}{\pi} \mathcal{P} \int_{-\infty}^{\infty} {\frac {H_I(\omega')} {\omega' - \omega}} \, d\omega'$
$H_I(\omega) = -\frac{1}{\pi} \mathcal{P} \int_{-\infty}^{\infty} {\frac {H_R(\omega')} {\omega' - \omega}} \, d\omega'$

がクラマース・クローニッヒの関係式である。（$\mathcal{P}$はコーシーの主値をとることを表す。)

後述するインパルス応答h(t) が恒に実数であるという条件を付けると、周波数応答関数の実部は偶関数、虚部は奇関数になる。 これを用いて積分範囲を正の部分にするようにクラマース・クローニッヒの関係式を変形すると

$H_R(\omega) = \frac{2}{\pi} \mathcal{P} \int_{0}^{\infty} {\frac {\omega' H_I(\omega')} {{\omega'}^2 - \omega^2}} \, d\omega'$
$H_I(\omega) = -\frac{2}{\pi} \mathcal{P} \int_{0}^{\infty} {\frac {\omega H_R(\omega')} {{\omega'}^2 - \omega^2}} \, d\omega'$

となる。

因果律からの導出

クラマース・クローニッヒの関係式は、刺激よりも前に応答は起こりえないという因果律から導かれる。

線形応答においてはt=0におけるインパルスδ(t)に対する応答h(t)が決まれば、任意の刺激に対する応答が決定される。h(t)を偶関数h_e(t)と奇関数h_o(t)の和の形

$h(t) = \frac{h(t)+h(-t)}{2} + \frac{h(t)-h(-t)}{2} \equiv h_\mathrm{e}(t) + h_\mathrm{o}(t)$

に分解すると、因果律よりt<0でh(t)=0なのでh_o(t)=h_e(t)·sgn(t)、h_e(t)=h_o(t)·sgn(t)となる。

インパルス応答をフーリエ変換して周波数応答関数を求めると、

$\begin{align} H(\omega) &= \int_{-\infty}^{\infty} dt\, h_\mathrm{e}(t) e^{i\omega t} + \int_{-\infty}^{\infty} dt\, h_\mathrm{o}(t) e^{i\omega t} \\ &= \int_{-\infty}^{\infty} dt\, h_\mathrm{e}(t) \cos(\omega t) + i \int_{-\infty}^{\infty} dt\, h_\mathrm{o}(t) \sin(\omega t) \end{align}$

となり、偶関数部h_e(t)のフーリエ変換は周波数応答関数の実部、奇関数部h_o(t)のフーリエ変換は周波数応答関数の虚部にあたることが分かる。

それぞれに対して積関数のフーリエ変換が畳み込みになることを使えば、クラマース・クローニッヒの関係式が導かれる。ここで、$\widehat{\sgn}(\omega)$は符号関数のフーリエ変換を表す。

$\begin{align} H_\mathrm{R}(\omega) &= \int_{-\infty}^{\infty} dt\, h_\mathrm{e}(t) e^{i\omega t} = \int_{-\infty}^{\infty} dt\, h_\mathrm{o}(t) \sgn(t) e^{i\omega t} \\ &= \frac{1}{2\pi} \int_{-\infty}^{\infty} d\omega'\, i H_\mathrm{I}(\omega') \widehat{\sgn}(\omega-\omega') \\ &= \frac{1}{\pi} {\mathcal P} \int_{-\infty}^{\infty} d\omega'\, \frac{H_\mathrm{I}(\omega')}{\omega'-\omega} \\ i H_\mathrm{I}(\omega) &= \int_{-\infty}^{\infty} dt\, h_\mathrm{e}(t) \sgn(t) e^{i\omega t} \\ &= \frac{1}{2\pi} \int_{-\infty}^{\infty} d\omega'\, H_\mathrm{R}(\omega') \widehat{\sgn}(\omega-\omega') \\ &= \frac{1}{\pi i} {\mathcal P} \int_{-\infty}^{\infty} d\omega'\, \frac{H_\mathrm{R}(\omega')}{\omega'-\omega} \end{align}$

複素関数を用いた導出

クラマース・クローニッヒの関係式の導出における複素積分

またH(ω) を複素平面に解析接続した複素関数H(z) が、実軸より上側で正則かつ|z|→∞ で一様にH(z)→0 であるときにはH(ω) がクラマース・クローニッヒの関係式を満たすことを示すことができる。

H(z)/(z-ω) を複素平面上で、以下の4つの区間からなる閉曲線上で複素積分する。

1. 実軸上の(-R,0)→(ω-r,0)
2. (ω,0)を中心とする半径r の半円(ω-r,0)→(ω+r,0)
3. 実軸上の(ω+r,0)→(R,0)
4. 原点を中心とする半径Rの半円(R,0)→(-R,0)

実軸より上側で正則であるという条件から、コーシーの積分定理によりこの閉曲線上の積分は0になる。 ここでR→∞、r→0の極限をとると区間4の積分は|z|→∞で一様にH(z)→0の条件より0となる。 区間2の積分はr→0で-iπH(ω)となる。 したがって区間1と3の積分の和はR→∞、r→0の極限で

$\lim_{R \to \infty , r \to 0} \left \{ \int_{-R}^{\omega-r} {\frac {H(z)}{z-\omega}} \, dz + \int_{\omega+r}^{R} {\frac {H(z)}{z-\omega}} \, dz \right \} = \mathcal{P} \int_{-\infty}^{\infty} \frac{H(z)}{z-\omega} \, dz = i\pi H(\omega)$

この式の実部と虚部を比較することでクラマース・クローニッヒの関係式が導出される。

応用

クラマース・クローニッヒの関係式を用いることで周波数応答関数の実部か虚部の一方からもう一方を計算で求めることが可能になる。これをクラマース・クローニッヒ解析という。これは複素屈折率の測定などで用いられる。

フーリエ変換核磁気共鳴分光法で行なわれるように、インパルス応答であるFIDのn 点のサンプリングデータから周波数応答関数(スペクトル)を離散フーリエ変換で求めた場合、それぞれ独立に実部と虚部のスペクトルデータがn 点得られる。しかし実際にはクラマース・クローニッヒの関係式により実部と虚部は独立ではないので、虚部の自由度を実部に移してスペクトルの分解能を2倍にすることが可能である。（ただしSN比はその分低下する。）このテクニックは時間的制約の大きい二次元NMRなどにおいてデータ点数を二倍にするゼロフィリングなどを用いて実行される。

関連項目

• ヒルベルト変換