Wignerfunktion

Die Wignerfunktion (Wigner-Quasi-Wahrscheinlichkeitsverteilung) wurde 1932 von Eugene Wigner eingeführt, um Quantenkorrekturen der klassischen Statistischen Mechanik zu untersuchen. Das Ziel bestand darin, die Wellenfunktion der Schrödingergleichung durch eine Wahrscheinlichkeitsverteilung im Phasenraum zu ersetzen. Eine solche Verteilung wurde unabhängig 1931 von Hermann Weyl als Dichtematrix in der Darstellungstheorie gefunden. Ein weiteres Mal wurde sie durch J. Ville 1948 als quadratische (als Funktion des Signals) Darstellung der örtlichen Zeit-Frequenz Energie eines Signals entdeckt. Diese Verteilung ist auch unter den Namen „Wignerfunktion“, „Wigner-Weyl-Transformation“ oder „Wigner-Ville-Verteilung“ bekannt. Sie findet Anwendung in der Statistischen Mechanik, Quantenchemie, Quantenoptik, klassischen Optik und der Signalanalyse, sowie in einer Reihe von Gebieten der Elektrotechnik, Seismologie, Biologie und Motorendesign.

Ein klassisches Teilchen besitzt eine definierte Lage und einen definierten Impuls und kann daher zu jedem Zeitpunkt durch einen Punkt im Phasenraum dargestellt werden. Für ein Ensemble von Teilchen lässt sich eine Wahrscheinlichkeitsverteilung definieren, die die Wahrscheinlichkeit angibt, mit der sich ein Teilchen an einem bestimmten Ort im Phasenraum befindet. Dies ist jedoch nicht für ein Quantenteilchen möglich, welches der Unschärferelation genügen muss. Stattdessen lässt sich eine Quasi-Wahrscheinlichkeitsverteilung definieren, die notwendigerweise nicht alle Eigenschaften einer gewöhnlichen Wahrscheinlichkeitsverteilung aufweist. Die Wignerverteilung kann zum Beispiel für nicht-klassische Zustände negative Werte annehmen und kann daher verwendet werden, um solche Zustände zu identifizieren.

Bei gegebener Wellenfunktion $\psi$ wird die zugehörige Wignerverteilung $P(x,p)$ definiert als:

P(x,p)={\frac {1}{\pi \hbar }}\int _{-\infty }^{\infty }dy\,\psi ^{*}(x+y)\psi (x-y)e^{2ipy/\hbar }

mit dem Ort $x$ sowie Impuls $p$ . Letztere können aber auch ein beliebiges Paar konjugierter Variablen sein (z. B. Real- und Imaginärteil des elektrischen Feldes oder Frequenz und Dauer eines Signals). Die Verteilung ist symmetrisch in $x$ und $p$ :

P(x,p)={\frac {1}{\pi \hbar }}\int _{-\infty }^{\infty }dq\,\phi ^{*}(p+q)\phi (p-q)e^{-2ixq/\hbar }

wobei $\phi$ die Fouriertransformierte von $\psi$ , also die Impulswellenfunktion ist.

Für ein Zustandsgemisch ist:

P(x,p)={\frac {1}{\pi \hbar }}\int _{-\infty }^{\infty }dy\,\langle x-y|{\hat {\rho }}|x+y\rangle e^{2ipy/\hbar }

wobei ${\hat {\rho }}$ den Dichteoperator bezeichnet.

Mathematische Eigenschaften

1. $P(x,p)$ ist reell.

2. Die Wahrscheinlichkeitsverteilungen von $x$ und $p$ ergeben sich aus:

$\int _{-\infty }^{\infty }dp\,P(x,p)=\langle x|{\hat {\rho }}|x\rangle$ . Wenn sich das System durch einen reinen Zustand beschreiben lässt, folgt $\int _{-\infty }^{\infty }dp\,P(x,p)=|\psi (x)|^{2}$ .
$\int _{-\infty }^{\infty }dx\,P(x,p)=\langle p|{\hat {\rho }}|p\rangle$ . Wenn sich das System durch einen reinen Zustand beschreiben lässt, folgt $\int _{-\infty }^{\infty }dx\,P(x,p)=|\phi (p)|^{2}$ .
$\int _{-\infty }^{\infty }dx\int _{-\infty }^{\infty }dp\,P(x,p)=\operatorname {Spur} ({\hat {\rho }})$ .
Die Spur von ${\hat {\rho }}$ ist gewöhnlich gleich 1.
Aus 1. und 2. folgt, dass $P(x,p)$ an einigen Stellen negativ ist, falls es sich nicht um einen kohärenten Zustand (oder eine Mischung kohärenter Zustände) oder nicht um einen gequetschten Vakuumzustand handelt.

3. $P(x,p)$ besitzt die folgenden Symmetrien:

Zeitumkehr: $\psi (x)\rightarrow \psi (x)^{*}\Rightarrow P(x,p)\rightarrow P(x,-p)$ .
Raumspiegelung: $\psi (x)\rightarrow \psi (-x)\Rightarrow P(x,p)\rightarrow P(-x,-p)$ .

4. $P(x,p)$ ist Galilei-invariant:

$\psi (x)\rightarrow \psi (x+y)\Rightarrow P(x,p)\rightarrow P(x+y,p)$ .
Sie ist nicht invariant unter der Lorentztransformation.

5. Die Bewegungsgleichung eines Punktes im Phasenraum ist ohne Kräfte klassisch:

${\frac {\partial P(x,p)}{\partial t}}={\frac {-p}{m}}{\frac {\partial P(x,p)}{\partial x}}$ .

6. Die Überlappung zweier Zustände berechnet sich als:

$|\langle \psi |\theta \rangle |^{2}=2\pi \hbar \int _{-\infty }^{\infty }dx\,\int _{-\infty }^{\infty }dp\,P_{\psi }(x,p)P_{\theta }(x,p)$ .

7. Operatoren und Erwartungswerte (Mittelwerte) werden folgendermaßen berechnet:

$A(x,p)=\int _{-\infty }^{\infty }dy\,\langle x-y/2|{\hat {A}}|x+y/2\rangle e^{ipy/\hbar }$ .
$\langle \psi |{\hat {A}}|\psi \rangle =Tr({\hat {\rho }}{\hat {A}})=\int _{-\infty }^{\infty }dx\,\int _{-\infty }^{\infty }dpP(x,p)A(x,p)$ .

8. Damit $P(x,p)$ physikalische (positive) Dichtematrizen beschreibt, muss gelten:

$\int _{-\infty }^{\infty }dx\,\int _{-\infty }^{\infty }dp\,P(x,p)P_{\theta }(x,p)\geq 0$

wobei $|\theta \rangle$ ein reiner Zustand ist.

Anwendung außerhalb der Quantenmechanik

In der Modellierung optischer Systeme wie Teleskope oder Glasfasern in Geräten zur Telekommunikation füllt die Wignerfunktion die Lücke zwischen dem einfachen Ray tracing und der vollständigen Wellenanalyse des Systems. Dabei wird in der Näherung kleiner Winkel (paraxial) $p/\hbar$ durch $k=|k|\sin \theta \approx |k|\theta$ ersetzt. In diesem Zusammenhang ist die Wignerfunktion die beste Näherung zu einer Beschreibung des Systems mit Hilfe von Strahlen mit dem Ort $x$ und Winkel $\theta$ unter Einschluss von Interferenzeffekten. Falls diese negative Werte an irgendeinem Punkt annimmt, lässt sich das System nicht mit der einfachen Raytracing-Methode beschreiben.
In der Signalanalyse wird ein zeitabhängiges elektrisches Signal, mechanische Vibrationen oder Schallwellen durch die Wignerfunktion dargestellt. Dabei wird $x$ durch die Zeit und $p/\hbar$ durch die Kreisfrequenz $\omega =2\pi f$ ersetzt. Hierbei bezeichnet $f$ die gewöhnliche Frequenz.
Auf dem Gebiet der ultraschnellen Optik werden kurze Laserpulse durch die Wignerfunktion mittels der gleichen Substitution von Frequenz und Zeit charakterisiert. Bestimmte Pulseigenschaften wie ein Chirp (Änderung der Frequenz mit der Zeit) lassen sich durch die Wignerfunktion darstellen.
In der Quantenoptik werden $x$ und $p/\hbar$ durch $X$ und $P$ Quadraturen ersetzt, welche den Real- und Imaginärteil des elektrischen Feldes bezeichnen (siehe kohärenter Zustand).

Messung

Tomographie
Homodyne Detektion
FROG Frequency-resolved optical gating

Weitere Quasi-Wahrscheinlichkeitsverteilungen

Die Wignerverteilung war die erste Quasi-Wahrscheinlichkeitsverteilung, aber viele weitere mit verschiedenen Vorteilen folgten, darunter:

Historische Anmerkung

Wie in der Einleitung angeführt, wurde die Wignerfunktion viele Male unabhängig in verschiedenen Zusammenhängen gefunden. Tatsächlich scheint es so, dass Wigner nicht wusste, dass diese Funktion selbst innerhalb der Quantentheorie zuvor von Heisenberg und Dirac eingeführt worden war. Diese sahen jedoch nicht deren Bedeutung und glaubten, dass diese Funktion lediglich eine Näherung der exakten quantenmechanischen Beschreibung des Systems war. Im Übrigen wurde Dirac später der Schwager von Wigner (siehe Literatur).

Literatur

E.P. Wigner: On the quantum correction for thermodynamic equilibrium. In: Phys. Rev. Band 40, Juni 1932, S. 749–759.
H. Weyl: Quantenmechanik und Gruppentheorie. In: Z. Phys. Band 46, Juni 1927, S. 1–46.

H. Weyl: Gruppentheorie und Quantenmechanik. S. Hirzel, Leipzig 1928.

H. Weyl: The Theory of Groups and Quantum Mechanics. Dover, New York 1931.
J. Ville: Théorie et Applications de la Notion de Signal Analytique. In: Cables et Transmission. Band 2, Nr. 1, 1948, S. 61–74.
W. Heisenberg: Über die inkohärente Streuung von Röntgenstrahlen. In: Physik. Zeitschr. Band 32, 1931, S. 737–740.

P.A.M. Dirac: Note on exchange phenomena in the Thomas atom. In: Proc. Camb. Phil. Soc. Band 26, 1930, S. 376–395.

C. Zachos, D. Fairlie, T. Curtright: Quantum Mechanics in Phase Space. World Scientific, Singapore 2005, S. 737–740.

Weblinks

A gallery of Wigner functions at Department of Physics and Astronomy, University of Calgary
A Gallery of Quantum States
Wigner-Verteilung