Kramers-Kronig-Beziehungen

Die Kramers-Kronig-Beziehungen, auch Kramers-Kronig-Relation (nach ihren Entdeckern Hendrik Anthony Kramers und Ralph Kronig), setzen Real- und Imaginärteil bestimmter meromorpher Funktionen in Form einer Integralgleichung miteinander in Beziehung. Sie stellen damit einen Spezialfall der Hilbert-Transformation dar.

Eine wichtige Anwendung ist der Zusammenhang zwischen der Absorption und der Dispersion eines sich in einem Medium ausbreitenden „Lichtstrahls“. Weitere Anwendungen gibt es in der Hochenergiephysik und in den Ingenieurwissenschaften.

Mathematische Formulierung

Sei $F:{\mathbb {C}}\rightarrow {\mathbb {C}}$ eine meromorphe Funktion, deren Polstellen in der unteren Halbebene liegen. Dieser Forderung an die Lage der Polstellen entspricht physikalisch das Kausalitätspostulat. Ferner seien $\operatorname {Re}\,F|_{{\mathbb {R}}}$ bzw. $\operatorname {Im}\,F|_{{\mathbb {R}}}$ Real- und Imaginärteil der Funktion . Es sei vorausgesetzt, dass diese beiden Funktionen gerade bzw. ungerade sind. Das bedeutet, dass durch Fourierintegration nicht aus einer beliebigen komplexen, sondern aus einer reellen Funktion gebildet werden kann.

In der Physik betrachtet man oft statt die Funktion $F/\mathrm {i}$ , wodurch sich die Voraussetzungen bezüglich gerade und ungerade vertauschen.

Schließlich sei $\lim _{{|z|\rightarrow \infty }}|F(z)|=0$ . Dann gelten für $x\in \mathbb {R}$ die folgenden als Kramers-Kronig-Beziehungen bezeichnete Gleichungen:

$\operatorname {Im}\,F(x)=-{\frac {2}{\pi }}\cdot \;{\mathrm {CH}}\,\int _{{0}}^{{+\infty }}{\frac {x\cdot \operatorname {Re}\,F(t)}{t^{2}-x^{2}}}{\mathrm {d}}t$

$\operatorname {Re}\,F(x)={\frac {2}{\pi }}\cdot \;{\mathrm {CH}}\,\int _{{0}}^{{+\infty }}{\frac {t\cdot \operatorname {Im}\,F(t)}{t^{2}-x^{2}}}{\mathrm {d}}t$

${\mathrm {CH}}$ bezeichnet den cauchyschen Hauptwert des auftretenden Integrals.

Real- und Imaginärteil der Funktion bedingen sich also gegenseitig durch Integration. Dies findet Anwendungen in der Optik und in der Systemtheorie wenn die Suszeptibilität eines Systems angibt, siehe Kausalität. Anwendungen finden sich auch in der Hochenergie-Physik bei den Dispersionsrelationen der S-Matrix.

Motivation (Ein Randwertproblem)

Auf der reellen Achse $\mathbb {R}$ sei eine stetige reelle Funktion $\,f$ vorgegeben, die analog zu $\operatorname {Re}\,F$ als gerade vorausgesetzt werden soll. Dazu soll eine in der ganzen oberen Halbebene holomorphe komplexe(!) Funktion $\,F$ so konstruiert werden, dass $\operatorname {Re}\,F|_{{\mathbb {R}}}{\stackrel {!}{=}}f$ gilt.

Es soll also ein Randwertproblem gelöst werden, wobei im Innern des betrachteten Gebietes $\,G$ , d.h. oberhalb von $\mathbb {R}$ , wegen der Holomorphie-Bedingung die cauchy-riemannschen Differentialgleichungen erfüllt werden müssen und auf dem Rand, $\partial G={\mathbb R}\,,$ eine stetige reelle Funktion, , vorgegeben ist, die dort angenommen werden soll.

Eine holomorphe Funktion kann nach dem Residuensatz dargestellt werden als:

$F(z)={\frac {1}{2\pi \mathrm {i} }}\left(\int _{HK_{r}}{\frac {F(t)}{t-z}}\mathrm {d} t+\int _{-r}^{r}{\frac {F(t)}{t-z}}\mathrm {d} t\right)$ ,

wobei $HK_{r}(0)$ den (positiv orientierten) Halbkreis in der oberen Halbebene mit Zentrum $0$ und Radius r>0 bezeichnet. Fällt nun im Unendlichen schnell genug ab, so reduziert sich im Grenzübergang $r \rightarrow \infty$ die Darstellung zu einem Integral über der reellen Achse, also:

$F(z)={\frac {1}{2\pi \mathrm {i} }}\int _{-\infty }^{+\infty }{\frac {F(t)}{t-z}}\mathrm {d} t$

Im Falle $\operatorname {Im}\,z=0$ , und weil $\,f$ bzw. $\operatorname {Re}\,F(t)$ eine gerade Funktion sein soll, ergibt sich schließlich

$\operatorname {Im}\,F(z)=-{\frac {1}{\pi }}\int _{{-\infty }}^{{+\infty }}{\frac {\operatorname {Re}\,F(t)}{t-z}}{\mathrm {d}}t=-{\frac {2}{\pi }}\int _{{0}}^{{+\infty }}{\frac {z\cdot \operatorname {Re}\,F(t)}{t^{2}-z^{2}}}{\mathrm {d}}t$ ,

wobei das auftretende Integral als cauchyscher Hauptwert zu interpretieren ist (Singularität für t=z ) und mit der Hilbert-Transformation von übereinstimmt. Der Residuensatz wird hierbei auf den Integrationsweg $[-r,z-\varepsilon ]\cdot HK_{\varepsilon }(z)\cdot [z+\varepsilon ,r]\cdot HK_{r}(0)$ angewendet. Diese Gleichung entspricht der einen Kramers-Kronig-Beziehung.

Man braucht jetzt zur Lösung des Randwertproblems nur die Beziehung $\operatorname {Re}\,F|_{{\mathbb {R}}}{=}f$ einzusetzen.

Für ungerade Funktionen verfährt man analog und erhält die andere Kramers-Kronig-Beziehung. Eine beliebige Funktion kann immer durch die Vorschrift $\,f=f_{+}+f_{-}$ , mit $f_{\pm }(t)={\frac {1}{2}}\left(f(t)\pm f(-t)\right)$ , in einen geraden bzw. ungeraden Anteil zerlegt werden. Der einfachste Fall einer meromorphen Funktion F(z) mit den vorausgesetzten Eigenschaften ist die lineare Antwortfunktion des gedämpften harmonischen Oszillators, $F(z):=1/(z^{2}-\omega _{0}^{2}+\mathrm {i} \gamma z),$ mit positiver Dämpfungskonstante $\gamma$ und positiver charakteristischer Oszillator-Kreisfrequenz $\omega _{0}.$

Anwendungen

Die Kramers-Kronig-Beziehungen werden dort angewendet, wo eine reelle gerade Funktion $F(\omega )$ – evtl. ungerade gemacht durch einen Zusatzfaktor $1/\mathrm {i}$ – zu einer holomorphen Funktion F(z) ergänzt werden soll. Dies dient meistens der Vereinfachung der auftretenden Rechnungen, insbesondere bei Wellenfunktionen, also hauptsächlich in der Signalverarbeitung und in der Optik, aber auch in der Statistischen Physik im Zusammenhang mit dem Fluktuations-Dissipations-Theorem. Auf diese Weise hängt die Absorption elektromagnetischer Wellen in einem Medium mit dem Brechungsindex zusammen. Es reicht also, die Abhängigkeit einer der beiden Größen von der Frequenz zu kennen, um die andere berechnen zu können.

Die von der Kreisfrequenz $\omega$ abhängige Permittivität $\varepsilon (\omega )$ lässt sich ausdrücken als Integral der von der Kreisfrequenz abhängigen Absorption:

$\operatorname {Re} (\varepsilon (\omega ))=1+{\frac {2}{\pi }}\cdot \;\mathrm {CH} \,\int \limits _{0}^{\infty }{\frac {\Omega \cdot \operatorname {Im} (\varepsilon (\Omega ))}{\Omega ^{2}-\omega ^{2}}}\,\mathrm {d} \Omega$

wobei

die reelle Kreisfrequenz $\Omega$ die Integrationsvariable ist
die ebenfalls reelle Variable $\omega$ eine charakteristische System-Kreisfrequenz darstellt
die Abkürzung ${\mathrm {CH}}\int _{0}^{\infty }\dots$ für den cauchyschen Hauptwert (engl. Cauchy principal value) des Integrals steht ( $\operatorname {Re} \varepsilon (\omega )$ ist gerade, $\operatorname {Im} \varepsilon (\omega )$ ungerade).

Eine alternative Betrachtungsweise ergibt sich mit dem Absorptionskoeffizienten $\alpha$ , dem Brechungsindex und der Lichtgeschwindigkeit :

$n(\omega )=1+{\frac {c}{\pi }}\cdot \;\mathrm {CH} \,\int \limits _{0}^{\infty }{{\alpha (\Omega )} \over {\Omega ^{2}-\omega ^{2}}}\,\mathrm {d} \Omega$

Dadurch lässt sich vor allem in der nichtlinearen Optik aus einer einfachen Absorptionsmessung die komplexe Form des Brechungsindex ableiten. Auch der Name der Dispersionsrelationen der Hochenergiephysik bezieht sich auf dieses Beispiel.

In den Ingenieurswissenschaften kommen die Kramers-Kronig-Beziehungen vor allem im Rahmen von impedanzspektroskopischen Messungen zum Einsatz, wo aus ihrer Nichterfüllung auf eine fehlerhafte Messung des Frequenzganges geschlossen wird.

Die Einschränkung der allgemeiner gültigen Kramers-Kronig-Beziehungen auf Zweipol-Systeme führt zum ZHIT-Algorithmus, der zur Validierung von Impedanzspektren elektrochemischer Systeme (Elektrochemische Impedanzspektroskopie) angewandt werden kann.

Basierend auf einem Artikel in:

Wikipedia.de