Faddeeva-Funktion

Die Faddeeva-Funktion (auch Kramp-Funktion oder relativistische Plasma-Dispersions-Funktion) ist eine skalierte komplexe komplementäre Fehlerfunktion,

w (z) := e^{- z^{2}} erfc (- i z) = e^{- z^{2}} (1 + \frac{2 i}{\sqrt{π}} \int_{0}^{z} e^{t^{2}} d t) .

Sie ist verwandt mit den Fresnel-Integralen, den Dawson-Integralen und dem Voigt-Profil. Die Funktion ist nach Wera Nikolajewna Faddejewa benannt.

Eigenschaften

Real- und Imaginärteil

Für genauere Betrachtungen lässt sich $w (z)$ mit $z = x + i y$ wie folgt zerlegen:

w (z) = V (x, y) + i L (x, y)

,

$V$ und $L$ stellen hierbei die reale und imaginäre Voigt-Funktion dar, da es sich bei $V (x, y)$ bis auf Vorfaktoren um das Voigt-Profil handelt.

Integraldarstellung

Die Faddeeva-Funktion besitzt die Integraldarstellung

w (z) = \frac{i}{π} \int_{- \infty}^{\infty} \frac{e^{- t^{2}}}{z - t} d t ℑ (z) > 0

sprich sie ist die Konvolution einer Gauß-Funktion und einer einfachen Polstelle.
Die reale und imaginäre Voigt-Funktion lassen sich in ähnlicher Weise darstellen:

V (x, y) = \frac{y}{π} \int_{- \infty}^{\infty} \frac{e^{- t^{2}}}{{(x - t)}^{2} + y^{2}} d t

L (x, y) = \frac{1}{π} \int_{- \infty}^{\infty} \frac{(x - t) \cdot e^{- t^{2}}}{{(x - t)}^{2} + y^{2}} d t

Verhalten bei Vorzeichenumkehr

Bei einer Vorzeichenumkehr von $z$ kann bei Berechnungen auf die folgenden Zusammenhänge zurückgegriffen werden:

w (- z) = 2 e^{- z^{2}} - w (z)

sowie

w (- z) = \overline{w (\overline{z})}

$\overline{z}$ ist die Konjugation von $z$ .

Ableitung

In manchen Anwendungen muss nicht nur die Faddeeva-Funktion selbst, sondern auch ihre Ableitungen berechnet werden, beispielsweise bei der Nichtlinearen Regression in der Spektroskopie. Ihre analytische Ableitung lautet:

\frac{d w (z)}{d z} = \frac{2 i}{\sqrt{π}} - 2 \cdot z \cdot w (z)

Dieser Ausdruck kann auch herangezogen werden, um die Änderungen im Real- und Imaginärteil der Faddeeva-Funktion $ℜ (w (z)) = ℜ_{w}$ und $ℑ (w (z)) = ℑ_{w}$ nachzuvollziehen. Im Prinzip muss dafür das Produkt $z \cdot w (z)$ eingehender betrachtet werden. Mit der obigen Definition des Arguments $z = x + i y$ , kann die Ableitung auch in ihre partiellen Ableitungen nach $x$ und $y$ zerlegt werden:

\frac{d ℜ_{w}}{d x} = 2 \cdot (y \cdot ℑ_{w} - x \cdot ℜ_{w}) = \frac{d ℑ_{w}}{d y}

\frac{d ℜ_{w}}{d y} = - 2 \cdot (\frac{1}{\sqrt{π}} - x \cdot ℑ_{w} - y \cdot ℜ_{w}) = - \frac{d ℑ_{w}}{d x}

\frac{d ℑ_{w}}{d x} = 2 \cdot (\frac{1}{\sqrt{π}} - x \cdot ℑ_{w} - y \cdot ℜ_{w}) = - \frac{d ℜ_{w}}{d y}

\frac{d ℑ_{w}}{d y} = 2 \cdot (y \cdot ℑ_{w} - x \cdot ℜ_{w}) = \frac{d ℜ_{w}}{d x}

Beziehungen zu anderen Funktionen

Dawsonsche Funktion

Es gilt folgende Beziehung zur Dawsonschen Funktion $D_{+} (x)$

w (x) = e^{- x^{2}} erfc (- i x) = e^{- x^{2}} + \frac{2 i}{\sqrt{π}} D_{+} (x) .

Komplementäre Fehlerfunktion

Für rein imaginäre Argumente $i y$ entspricht die Faddeeva-Funktion der skalierten Komplementären Fehlerfunktion $e r f c x$

w (i y) = e r f c x (y) = e^{y^{2}} e r f c (y)

,

mit der Komplementären Fehlerfunktion $e r f c$ .

Geschichte

Die Funktion wurde 1954 von Wera Faddejewa und Terentjew tabuliert. Sie erscheint als namenlose Funktion $w (z)$ im Standardwerk von Abramowitz-Stegun (1964), Formel 7.1.3. Der Name Faddeeva function wurde anscheinend 1990 von Poppe und Wijers eingeführt.

Implementierungen

Steven G. Johnson hat eine Implementierung als freie und offene Software veröffentlicht, die auf einer Kombination der Algorithmen 680 und 916 beruht. Sie liegt der Funktion scipy.special.wofz in der Python-Bibliothek SciPy zugrunde, und sie ist auch in Form einer C-Bibliothek libcerf verfügbar.
Für Matlab existiert eine öffentlich einsehbare Implementierung, die auf einer Approximation durch Fourierreihen sowie einer Kettenbruchdarstellung basiert.

Literatur

W. Gautschi ACM Transactions on Mathematical Software (1969?): ACM Algorithmus 363.
W. Gautschi SIAM J. Numer. Anal. 7, 187 (1970).
G. P. M. Poppe, C. M. J. Wijers, ACM Transactions on Mathematical Software 16, 38–46 (1990): ACM Algorithm 680.
J. A. C. Weideman, SIAM J. Numer. Anal. 31, 1497–1518 (1994): Besonders kompakter Algorithmus in 8 Zeilen Matlab.
M. R. Zaghloul and A. N. Ali, ACM Transactions on Mathematical Software 38, 15 (2011): ACM Algorithm 916.
S. M. Abrarov and B. M. Quine, Appl. Math. Comp. 218, 1894–1902 (2011).
S. M. Abrarov and B. M. Quine, Arxiv, Preprint 2012

Quellen

1 2 3 4 V. G. Avetisov: A Least-Squares Fitting Technique for Spectral Analysis of Direct and Frequency-Modulation Lineshapes. In: Fakultät für Physik der Universität Lund (Hrsg.): Lund Reports in Atomic Physics. LRAP-186, 1995 (lu.se [PDF]).
↑ W. Gautschi: Efficient Computation of the Complex Error Function. In: Society for Industrial and Applied Mathematics (Hrsg.): SIAM Journal on Numerical Analysis. Band 7, Nr. 1, 1970, S. 187–198 (nasa.gov [PDF]).
↑ National Institute of Standards and Technology (NIST): 7 Error functions, Dawson's and Fresnel integrals - 7.10 Derivatives, 15. März 2023, abgerufen am 14. Mai 2023
↑ J. H. McCabe: A Continued Fraction Expansion, with a Truncation Error Estimate, for Dawson's Integral. In: American Mathematical Society (Hrsg.): Mathematics of Computation. Band 28, Nr. 127, 1974, S. 811–816.
↑ V. N. Faddeeva, N. N. Terent'ev: Tables of values of the function $w (z) = \exp (- z^{2}) (1 + 2 i / \sqrt{π} \int_{0}^{z} \exp (t^{2}) d t)$ for complex argument. Gosud. Izdat. Teh.-Teor. Lit., Moscow, 1954; English transl., Pergamon Press, New York, 1961.
↑ Google-Scholar-Recherche laut engl. Wikipedia.
↑ Faddeeva Package, unter MIT-Lizenz.
↑ Archivierte Kopie (Memento vom 17. Februar 2013 im Webarchiv archive.today)
↑ Sanjar Abrarov: The Voigt/complex error function (second version), MATLAB Central File Exchange, 10. Juli 2016, abgerufen am 14. Mai 2023

[Avetisov_1995-1] 1 2 3 4 V. G. Avetisov: A Least-Squares Fitting Technique for Spectral Analysis of Direct and Frequency-Modulation Lineshapes. In: Fakultät für Physik der Universität Lund (Hrsg.): Lund Reports in Atomic Physics. LRAP-186, 1995 (lu.se [PDF]).

[2] W. Gautschi: Efficient Computation of the Complex Error Function. In: Society for Industrial and Applied Mathematics (Hrsg.): SIAM Journal on Numerical Analysis. Band 7, Nr. 1, 1970, S. 187–198 (nasa.gov [PDF]).

[3] National Institute of Standards and Technology (NIST): 7 Error functions, Dawson's and Fresnel integrals - 7.10 Derivatives, 15. März 2023, abgerufen am 14. Mai 2023

[4] J. H. McCabe: A Continued Fraction Expansion, with a Truncation Error Estimate, for Dawson's Integral. In: American Mathematical Society (Hrsg.): Mathematics of Computation. Band 28, Nr. 127, 1974, S. 811–816.

[5] V. N. Faddeeva, N. N. Terent'ev: Tables of values of the function $w (z) = \exp (- z^{2}) (1 + 2 i / \sqrt{π} \int_{0}^{z} \exp (t^{2}) d t)$ for complex argument. Gosud. Izdat. Teh.-Teor. Lit., Moscow, 1954; English transl., Pergamon Press, New York, 1961.

[6] Google-Scholar-Recherche laut engl. Wikipedia.

[7] Faddeeva Package, unter MIT-Lizenz.

[8] Archivierte Kopie (Memento vom 17. Februar 2013 im Webarchiv archive.today)

[9] Sanjar Abrarov: The Voigt/complex error function (second version), MATLAB Central File Exchange, 10. Juli 2016, abgerufen am 14. Mai 2023