Erzeugungs- und Vernichtungsoperator

aus Wikipedia, der freien Enzyklopädie
Wechseln zu: Navigation, Suche

Die Erzeugungs- und Vernichtungsoperatoren sind der Kern einer möglichen Lösung der Schrödingergleichung des harmonischen Oszillators. Diese Operatoren können auch dazu benutzt werden, Probleme mit quantenmechanischem Drehimpuls einfacher zu lösen. Ferner finden die Erzeugungs- und Vernichtungsoperatoren Verwendung bei der Quantisierung von Feldern (der sogenannten zweiten Quantisierung oder Besetzungszahl-Darstellung).

Es gibt eine Vielzahl alternativer Bezeichnungen wie Leiteroperatoren, Kletteroperatoren, Aufsteige- und Absteigeoperatoren sowie Hebe- und Senkoperatoren. Statt „Erzeugungsoperator“ wird manchmal auch Erschaffungsoperator verwendet. Im deutschsprachigen Raum werden darüber hinaus auch die Operatoren \sigma_+ und \sigma_-, die die Zustände eines Atoms ändern, als Erzeugungs- bzw. Vernichtungsoperatoren bezeichnet.

Das Problem des harmonischen Oszillators in der Quantenmechanik lässt sich mithilfe der Methode der Erzeugungs- und Vernichtungsoperatoren behandeln, die auch algebraische Methode genannt wird. Sie wurde hauptsächlich von Paul Dirac entwickelt. Für diesen Lösungsweg definiert man zwei Operatoren \hat a und \hat a^\dagger, die einem Oszillator jeweils ein Energiequant \hbar\omega entziehen oder hinzufügen. Deswegen werden sie Vernichtungs- und Erzeugungsoperator genannt.

Das Zirkumflex (\hat{ }-Symbol) über dem a symbolisiert, dass es sich dabei um einen Operator handelt. Damit gelten nicht dieselben Rechenregeln wie für Skalare, denn die Reihenfolge von Operatoren lässt sich beispielsweise im Allgemeinen nicht vertauschen. Im Folgenden wird auf das Zirkumflex-Symbol zugunsten der Übersichtlichkeit verzichtet, wenn nichts Anderes gesagt ist. Alle lateinischen Großbuchstaben, mit Ausnahme des E, sind Operatoren.

Definition[Bearbeiten]

Man definiert den Erzeugungsoperator a^{\dagger\,} und den dazu adjungierten Vernichtungsoperator a^{\,} über folgende Vertauschungsrelationen mit dem Besetzungszahloperator N:=a^\dagger a:

[N,a]=-a\ ,\quad[N,a^{\dagger}]=a^{\dagger}

Der Besetzungszahloperator N ist ein hermitescher Operator und hat daher reelle Eigenwerte n. Die zugehörige Eigenwertgleichung lautet, wobei \left| n \right\rangle Fock-Zustände sind:

N\left| n \right\rangle =n\left| n \right\rangle

Die Besetzungszahl n ist eine nichtnegative ganze Zahl, also n\in\mathbb{N}_{0}. Bei Fermionen ergibt sich hier noch eine Einschränkung auf die Werte 0 und 1.

Durch Anwendung von a^\dagger bzw. a^{\,} auf den Zustand \left| n \right\rangle erhält man den darüber- bzw. den darunterliegenden Zustand:

a^{\dagger}\left|n\right\rangle =c_{n}^{+}\left|n+1\right\rangle
a\,\left|n\right\rangle =c_{n}^{-}\left|n-1\right\rangle

Die Konstanten c_{n}^{+} und c_{n}^{-} sind davon abhängig, ob a^{\,} und a^\dagger die Kommutator- oder Antikommutator-Vertauschungsrelation erfüllen.

Details[Bearbeiten]

Im Folgenden werden verschiedene Eigenschaften von N abgeleitet. Die Eigenzustände \left| n \right\rangle seien normiert.

  • Der Besetzungszahloperator ist hermitesch, also selbstadjungiert:
N^{\dagger}=(a^{\dagger}a)^{\dagger}=a^{\dagger}a^{\dagger\dagger}=a^{\dagger}a=N
Somit hat N reelle Eigenwerte, die Besetzungszahlen n.
  • Die Eigenwerte sind nicht negativ: n\ge 0
n=\left\langle n|N|n \right\rangle =\left\langle n|a^{\dagger }a|n \right\rangle =\left\langle an|an \right\rangle =\left\| a \left| n \right\rangle  \right\|^{2}\ge 0
Die Ungleichung folgt aus der Tatsache, dass die Norm eines Vektors nicht-negativ ist.
  • Der kleinste Eigenwert ist 0
Der Zustand \left| 0 \right\rangle ist ein Vektor im Hilbertraum und darf nicht mit dem Nullvektor verwechselt werden, sondern wird lediglich mit der Zahl 0 nummeriert.
a^{\dagger }a\left| 0 \right\rangle =0\left| 0 \right\rangle = 0   und   \left\langle  0 | 0 \right\rangle =1
Wegen n\ge 0 muss gelten: a\left| 0 \right\rangle =0. Wendet man also auf den niedrigsten Zustand den Absteigeoperator an, so erhält man den Nullvektor. Dies lässt sich aber nicht umkehren: Durch Anwendung von a^{\dagger } auf den Nullvektor erhält man nicht den Grundzustand sondern wieder den Nullvektor.
  • Die Eigenwerte sind ganzzahlig: n\in \mathbb{N}_{0}
Angenommen die Eigenwerte wären nicht ganzzahlig, so ließen sich ausgehend von einem Eigenzustand durch mehrmalige Anwendung des Absteigeoperators Eigenzustände finden, die negative Eigenwerte besitzen. Dies ist aber ein Widerspruch zur Bedingung n\ge 0. Bei ganzzahligen Eigenwerten erreicht man irgendwann den Grundzustand und durch nochmaliges Anwenden den Nullvektor; ab hier bricht automatisch die Leiter ab.
  • Ist n Eigenwert, dann auch n+1
Na^{\dagger}\left|n\right\rangle =\left(a^{\dagger}N+[N,a^{\dagger}]\right)\left|n\right\rangle =\left(a^{\dagger}N+a^{\dagger}\right)\left|n\right\rangle =a^{\dagger}\left(N+1\right)\left|n\right\rangle =a^{\dagger}\left(n+1\right)\left|n\right\rangle =\left(n+1\right)a^{\dagger}\left|n\right\rangle
Wenn a^{\dagger }\left| n \right\rangle ungleich dem Nullvektor ist, erhält man somit einen neuen Eigenwert (n+1).
a^{\dagger }\left| n \right\rangle ist also Eigenzustand zu N mit Eigenwert (n+1) und somit proportional zu \left| n+1 \right\rangle:
a^{\dagger}\left|n\right\rangle =c_{n}^{+}\left|n+1\right\rangle
  • Ist n>0 Eigenwert, dann auch n-1
Na\left|n\right\rangle =\left(aN+\left[N,a\right]\right)\left|n\right\rangle =\left(aN-a\right)\left|n\right\rangle =a\left(N-1\right)\left|n\right\rangle =a\left(n-1\right)\left|n\right\rangle =\left(n-1\right)a\left|n\right\rangle
Wenn a\left| n \right\rangle ungleich dem Nullvektor ist, erhält man somit einen neuen Eigenwert (n-1).
a^{\,}\left| n \right\rangle ist also Eigenzustand zu N mit Eigenwert (n-1) und somit proportional zu \left| n-1 \right\rangle:
a\,\left|n\right\rangle =c_{n}^{-}\left|n-1\right\rangle

Bosonische Kletteroperatoren[Bearbeiten]

Im bosonischen Fall erfüllen a und a^\dagger die Kommutator-Vertauschungsrelationen:

[a,a^{\dagger}]=1\ ,\quad[a,a]=[a^{\dagger},a^{\dagger}]=0

Somit

a^{\dagger}\left|n\right\rangle =\sqrt{n+1}\left|n+1\right\rangle
a\,\left|n\right\rangle =\sqrt{n}\left|n-1\right\rangle

Im bosonischen Fall können die Besetzungszahlen n beliebig groß werden: n\in\mathbb{N}_{0}.

Details[Bearbeiten]

  • Zunächst ist zu prüfen, ob die obigen Voraussetzungen erfüllt werden:
[N,a]=[a^{\dagger}a,a]=\underbrace{[a^{\dagger},a]}_{-1}a+a^{\dagger}\underbrace{[a,a]}_{0}=-a
[N,a^{\dagger}]=[a^{\dagger}a,a^{\dagger}]=\underbrace{[a^{\dagger},a^{\dagger}]}_{0}a+a^{\dagger}\underbrace{[a,a^{\dagger}]}_{1}=a^{\dagger}
  • Mit a^{\dagger} lässt sich der nächste über \left| n \right\rangle liegende Zustand konstruieren a^{\dagger }\left| n \right\rangle =c\left| n+1 \right\rangle . Der Faktor c ergibt sich aus folgender Rechnung mit dem Kommutator aa^{\dagger}-a^{\dagger}a=1:
\left|c\right|^{2}=\left\Vert c\left|n+1\right\rangle \right\Vert ^{2}=\left\Vert a^{\dagger}\left|n\right\rangle \right\Vert ^{2}=\left\langle a^{\dagger}n|a^{\dagger}n\right\rangle =\left\langle n|aa^{\dagger}|n\right\rangle =\left\langle n|a^{\dagger}a+1|n\right\rangle =\left\langle n|N+1|n\right\rangle =n+1
c=\sqrt{n+1} e^{i\varphi}, die Phase \varphi kann aber vernachlässigt werden, sodass c=\sqrt{n+1}.
  • Mit a lässt sich der unter \left| n \right\rangle liegende Zustand konstruieren a\left| n \right\rangle =c\left| n-1 \right\rangle . Der Faktor c ergibt sich aus folgender Rechnung:
\left|c\right|^{2}=\left\Vert c\left|n-1\right\rangle \right\Vert ^{2}=\left\Vert a\left|n\right\rangle \right\Vert ^{2}=\left\langle an|an\right\rangle =\left\langle n|a^{\dagger}a|n\right\rangle =\left\langle n|N|n\right\rangle =n\quad\Rightarrow\quad c=\sqrt{n}
  • Alle Eigenzustände lassen sich vom Grundzustand ausgehend konstruieren:
\left|n\right\rangle =\frac{1}{\sqrt{n}}a^{\dagger}\left|n-1\right\rangle =\frac{1}{\sqrt{n!}}\left(a^{\dagger}\right)^{n}\left|0\right\rangle \ ,\quad n\in\mathbb{N}_{0}
Auf diese Weise erhält man einen vollständigen diskreten Satz von Eigenzuständen

Fermionische Kletteroperatoren[Bearbeiten]

Im fermionischen Fall erfüllen \,a\, und a^\dagger die Anti-Kommutator-Vertauschungsrelationen:

\{a,a^{\dagger}\}=1\ ,\quad\{a,a\}=\{a^{\dagger},a^{\dagger}\}=0\quad\Rightarrow\quad a^{2}=a^{\dagger\, 2}=0

Somit

a^{\dagger}\left|n\right\rangle =\left(1-n\right)\left|n+1\right\rangle
a\,\left|n\right\rangle =n\left|n-1\right\rangle

Im fermionischen Fall können die Besetzungszahlen n nur die Werte 0 oder 1 annehmen.

Details[Bearbeiten]

  • Mit a\, a=0 und a^{\dagger}a^{\dagger}=0 ist N^{2}=N:
N^{2}=a^{\dagger}aa^{\dagger}a=a^{\dagger}(\underbrace{\{a,a^{\dagger}\}}_{1}-a^{\dagger}a)a=a^{\dagger}a-\underbrace{a^{\dagger}a^{\dagger}}_{0}\underbrace{aa}_{0}=a^{\dagger}a=N
Der Besetzungszahloperator hat also nur die Eigenwerte 0 und 1 und die Eigenzustände |0\rangle und |1\rangle:
n^{2}\left|n\right\rangle =N^{2}\left|n\right\rangle =N\left|n\right\rangle =n\left|n\right\rangle \quad\Rightarrow\quad n^{2}=n\quad\Rightarrow\quad n\in\{0,1\}
  • Zunächst ist zu prüfen, ob die obigen Voraussetzungen erfüllt werden:
[N,a]=[a^{\dagger}a,a]=\underbrace{[a^{\dagger},a]}_{2a^{\dagger}a-\{a,a^{\dagger}\}}a+a^{\dagger}\underbrace{[a,a]}_{0}=2a^{\dagger}\underbrace{aa}_{0}-\underbrace{\{a,a^{\dagger}\}}_{1}a=-a
[N,a^{\dagger}]=[a^{\dagger}a,a^{\dagger}]=\underbrace{[a^{\dagger},a^{\dagger}]}_{0}a+a^{\dagger}\underbrace{[a,a^{\dagger}]}_{\{a,a^{\dagger}\}-2a^{\dagger}a}=a^{\dagger}\underbrace{\{a,a^{\dagger}\}}_{1}-2\underbrace{a^{\dagger}a^{\dagger}}_{0}a=a^{\dagger}
  • Mit a^{\dagger} lässt sich der nächste über \left| n \right\rangle liegende Zustand konstruieren a^{\dagger }\left| n \right\rangle =c\left| n+1 \right\rangle . Der Faktor c ergibt sich aus folgender Rechnung mit dem Anti-Kommutator aa^{\dagger}+a^{\dagger}a=1:
\left|c\right|^{2}=\left\Vert c\left|n+1\right\rangle \right\Vert ^{2}=\left\Vert a^{\dagger}\left|n\right\rangle \right\Vert ^{2}=\left\langle a^{\dagger}n|a^{\dagger}n\right\rangle =\left\langle n|aa^{\dagger}|n\right\rangle =\left\langle n|1-a^{\dagger}a|n\right\rangle =\left\langle n|1-N|n\right\rangle =1-n\quad\Rightarrow\quad c=\sqrt{1-n}
Da n nur 0 oder 1 sein kann, ist c=\sqrt{1-n}=1-n=\delta_{0,n} (dabei ist \delta_{i,j} das Kronecker-Delta).
  • Mit a lässt sich der unter \left| n \right\rangle liegende Zustand konstruieren a\left| n \right\rangle =c\left| n-1 \right\rangle . Der Faktor c ergibt sich aus folgender Rechnung:
\left|c\right|^{2}=\left\Vert c\left|n-1\right\rangle \right\Vert ^{2}=\left\Vert a\left|n\right\rangle \right\Vert ^{2}=\left\langle an|an\right\rangle =\left\langle n|a^{\dagger}a|n\right\rangle =\left\langle n|N|n\right\rangle =n\quad\Rightarrow\quad c=\sqrt{n}
Da n nur 0 oder 1 sein kann, ist c=\sqrt{n}=n=\delta_{1,n}.
  • Alle Eigenzustände lassen sich vom Grundzustand ausgehend konstruieren:
\left|n\right\rangle =\frac{1}{\sqrt{n}}a^{\dagger}\left|n-1\right\rangle =\frac{1}{\sqrt{n!}}\left(a^{\dagger}\right)^{n}\left|0\right\rangle \ ,\quad n\in\{0,1\}
Auf diese Weise erhält man einen vollständigen diskreten Satz von Eigenzuständen

Beispiel für bosonische Kletteroperatoren: Harmonischer Oszillator[Bearbeiten]

Der Hamiltonoperator H des harmonischen Oszillators lautet

H = \frac{P^{2}}{2m}+\frac{m\omega ^{2}Q^{2}}{2}

P Impulsoperator, Q Ortsoperator, m Masse, \omega Eigenfrequenz

Im Folgenden ist die stationäre Schrödingergleichung zu lösen:

H\left| \Psi _{n} \right\rangle =E_{n}\left| \Psi _{n} \right\rangle

E_n Energieeigenwert, \left| \Psi _{n} \right\rangle Energieeigenzustand

Hamiltonoperator umformen[Bearbeiten]

Der Hamiltonoperator lässt sich umformen:

H=\frac{P^{2}}{2m}+\frac{m\omega ^{2}Q^{2}}{2}=\hbar \omega \left( \frac{P^{2}}{2\hbar m\omega }+\frac{m\omega Q^{2}}{2\hbar } \right)

Es werden zwei neue Operatoren definiert:

\tilde{P}:=\frac{P}{\sqrt{2\hbar m\omega }}   und   \tilde{Q}:=\sqrt{\frac{m\omega }{2\hbar }}Q

Der Hamiltonoperator ausgedrückt mit den neuen Operatoren:

H=\hbar \omega \left( \tilde{P}^{2}+\tilde{Q}^{2} \right)

Man versucht nun, den Inhalt der Klammer als Produkt zu schreiben, also (i ist die imaginäre Einheit)

\left( u-iv \right)\left( u+iv \right)=u^{2}+v^{2}+iuv-ivu=u^{2}+v^{2}

Da aber u und v Operatoren sind, die nicht vertauschen, gilt hier das letzte Gleichheitszeichen nicht. Um zwei Operatoren miteinander zu vertauschen, ist der Kommutator vonnöten: \tilde{Q}\tilde{P}=\tilde{P}\tilde{Q}-\left[ \tilde{P},\tilde{Q} \right]

\begin{array}{rcl}
H & = & \hbar\omega\left(\tilde{Q}^{2}+\tilde{P}^{2}\right)\\
 & = & \hbar\omega\left(\tilde{Q}^{2}+\tilde{P}^{2}+i\left[\tilde{Q},\tilde{P}\right]-i\left[\tilde{Q},\tilde{P}\right]\right)\\
 & = & \hbar\omega\left(\tilde{Q}^{2}+\tilde{P}^{2}+i\tilde{Q}\tilde{P}-i\tilde{P}\tilde{Q}-i\left[\tilde{Q},\tilde{P}\right]\right)\\
 & = & \hbar\omega\left((\tilde{Q}-i\tilde{P})(\tilde{Q}+i\tilde{P})-i\left[\tilde{Q},\tilde{P}\right]\right)\end{array}

Der Kommutator \left[ \tilde{Q},\tilde{P} \right] kann auf den Kommutator der ursprünglichen Operatoren Q und P zurückgeführt werden:

\left[\tilde{Q},\tilde{P}\right]=\sqrt{\frac{m\omega}{2\hbar}}\frac{1}{\sqrt{2\hbar m\omega}}\underbrace{\left[Q,P\right]}_{i\hbar}=\frac{i}{2}

Der Hamiltonoperator sieht nun folgendermaßen aus:

H=\hbar \omega \underbrace{\left( \tilde{Q}-i\tilde{P} \right)}_{a^{\dagger }}\underbrace{\left( \tilde{Q}+i\tilde{P} \right)}_{a}+\frac{1}{2}\hbar \omega

Jetzt werden die beiden Leiteroperatoren definiert:

a^{\dagger }:=\tilde{Q}-i\tilde{P}   Erzeugungsoperator
a:=\tilde{Q}+i\tilde{P}   Vernichtungsoperator

Häufig werden sie auch als a+ und a- geschrieben. Man beachte, dass die Leiteroperatoren nicht hermitesch sind, da a\ne a^{\dagger }.

Die Leiteroperatoren ausgedrückt durch Ortsoperator Q und Impulsoperator P:

a^{\dagger} =\sqrt{\frac{m\omega }{2\hbar }}\left( Q-\frac{i}{m\omega }P \right)
a=\sqrt{\frac{m\omega }{2\hbar }}\left( Q+\frac{i}{m\omega }P \right)

Umgekehrt ergibt sich für Q und P:

Q = \sqrt{\frac{\hbar}{2m\omega}}(a + a^{\dagger})
P = \frac{1}{i} \sqrt{\frac{m\omega\hbar}{2}}(a - a^{\dagger})

Mit den Leiteroperatoren schreibt sich der Hamiltonoperator:

H=\hbar \omega \left( a^{\dagger }a+\frac{1}{2} \right)

Eigenschaften der Erzeuger und Vernichter[Bearbeiten]

Zu bestimmen ist noch der Kommutator aus den beiden Leiteroperatoren:

\left[ a,a^{\dagger } \right]=aa^{\dagger }-a^{\dagger }a=\left( \tilde{Q}+i\tilde{P} \right)\left( \tilde{Q}-i\tilde{P} \right)-\left( \tilde{Q}-i\tilde{P} \right)\left( \tilde{Q}+i\tilde{P} \right)=2i\left( \tilde{P}\tilde{Q}-\tilde{Q}\tilde{P} \right)=2i\underbrace{\left[ \tilde{P},\tilde{Q} \right]}_{-i/2}=1

Da außerdem [a,a]=[a^{\dagger},a^{\dagger}]=0 gilt, handelt es sich bei den Kletteroperatoren des harmonischen Oszillators um bosonische Kletteroperatoren. Somit gelten alle obigen Eigenschaften für bosonische Kletteroperatoren.

Das Produkt a^{\dagger }a definiert den Besetzungszahloperator:

N=a^{\dagger }a

Lösung des Eigenwertproblems[Bearbeiten]

Der Hamiltonoperator lässt sich durch den Besetzungszahloperator ausdrücken:

H=\hbar \omega \left( N+\frac{1}{2} \right)

Das Eigenwertproblem H\left| \Psi _{n} \right\rangle =E_{n}\left| \Psi _{n} \right\rangle lässt sich auf die Eigenwertgleichung des Besetzungszahloperators N\left| n \right\rangle =n\left| n \right\rangle zurückführen.

\hbar\omega\left(N+\frac{1}{2}\right)\left|\Psi_{n}\right\rangle =E_{n}\left|\Psi_{n}\right\rangle \quad,\quad\hbar\omega\left(N+\frac{1}{2}\right)\left|n\right\rangle =\hbar\omega\left(n+\frac{1}{2}\right)\left|n\right\rangle

Die Eigenzustände von N sind auch Eigenzustände von H, da \left[ H,N \right]=0. Die Eigenwerte des Hamiltonoperators ergeben sich aus den Eigenwerten des Besetzungszahloperators N:

E_n=\hbar \omega \left( n+\frac{1}{2} \right)   und   \left| \Psi_n \right\rangle = \left| n \right\rangle

Eine besonders wichtige Eigenschaft der Kletteroperatoren ist diese:

a^{\dagger}\left|\Psi_n\right\rangle =\sqrt{n+1}\left|\Psi_{n+1}\right\rangle
a\,\left|\Psi_n\right\rangle =\sqrt{n}\left|\Psi_{n-1}\right\rangle

Ist \left| \Psi_n \right\rangle eine Lösung der Schrödingergleichung für die Energie E_n, so ist a^{\dagger}\left| \Psi_n \right\rangle eine Lösung für die Energie E_n +\hbar \omega und a\left| \Psi_n \right\rangle eine Lösung für die Energie E_n -\hbar \omega. Das bedeutet, dass man aus einer Lösung alle Lösungen erhalten kann, indem man einfach den Erzeugungs- oder Vernichtungsoperator auf diese Lösung anwendet. Dadurch wird eine neue Lösung für das benachbarte Energieniveau erzeugt, das um die Energie \hbar\omega verschoben ist.

Da der Besetzungszahloperator keine negativen Eigenwerte hat, können auch keine negativen Energieeigenwerte existieren. Es gibt also für die minimale Besetzungszahl n=0 eine Lösung \left| \Psi_0 \right\rangle, die auf einem minimalen Energieniveau sitzt (Nullpunktenergie):

E_0 = \frac{1}{2}\hbar\omega

Im Zustand \left| \Psi_n \right\rangle setzt sich die Energie E_n = \hbar\omega (n+\tfrac{1}{2}) zusammen aus der Nullpunktenergie \hbar\omega /2 und n Energiequanten der Größe \hbar\omega. Die Wirkung von a^{\dagger} überführt das System in einen Zustand mit der um \hbar\omega erhöhten Energie. Dies kann man als Erzeugung eines zusätzlichen Energiequants interpretieren, was den Namen Erzeugungsoperator verständlich macht. Analog überführt der Operator a das System in einen um ein Energiequant reduzierten Zustand. Es wird also ein Energiequant vernichtet, deswegen Vernichtungsoperator. Die Eigenwerte des Operators N geben an, wie viele Energiequanten in einem Eigenzustand angeregt sind. Die Besetzung eines Zustandes mit n Energiequanten erklärt den Namen Besetzungszahloperator.

Eigenfunktionen in Ortsdarstellung[Bearbeiten]

Wendet man also auf den niedrigsten Zustand den Absteigeoperator an, so erhält man den Nullvektor a\left| \Psi_0 \right\rangle = 0. Dies lässt sich aber nicht umkehren: Durch Anwendung von a^{\dagger } auf den Nullvektor erhält man nicht den Grundzustand sondern wieder den Nullvektor a^{\dagger }a\left| \Psi_0 \right\rangle =0. Dies liefert eine Gleichung für den Grundzustand:

0=a\left| 0 \right\rangle =\left( \tilde{Q}+i\tilde{P} \right)\left| \Psi _{0} \right\rangle =\left( \sqrt{\frac{m\omega }{2\hbar }}Q+i\frac{P}{\sqrt{2\hbar m\omega }}\right) \left| \Psi _{0} \right\rangle

In der Ortsdarstellung kann man obige Operatorgleichung als Differentialgleichung darstellen und lösen: P=-i\hbar \frac{\mathrm{d}}{\mathrm{d}x} und Q=x

\left( \sqrt{\frac{m\omega }{2\hbar }}x+\sqrt{\frac{\hbar }{2m\omega }}\frac{\mathrm{d}}{\mathrm{d}x} \right)\Psi _{0}\left( x \right)=0   liefert normiert   \Psi _{0}\left( x \right)=\left( \frac{m\omega }{\pi \hbar }\right)^{\frac{1}{4}}\exp \left( -\frac{m\omega }{2\hbar}x^{2} \right)

Durch Anwendung des Aufsteigeoperators auf die Lösung des Grundzustands erhält man alle höheren Eigenfunktionen:

\left|\Psi_{n}\right\rangle =\frac{1}{\sqrt{n}}a^{\dagger}\left|\Psi_{n-1}\right\rangle =\frac{1}{\sqrt{n!}}(a^{\dagger})^{n}\left|\Psi_{0}\right\rangle

In Ortsdarstellung erhält man somit:

\Psi_{n}(x)=\frac{1}{\sqrt{n!}}\left(\sqrt{\frac{m\omega}{2\hbar}}x-\sqrt{\frac{\hbar}{2m\omega}}\frac{\mathrm{d}}{\mathrm{d}x}\right)^{n}\Psi_{0}\left(x\right)

Matrixdarstellung bosonischer Kletteroperatoren[Bearbeiten]

Die Eigenzustände des Besetzungszahloperators \left| n \right\rangle bilden ein vollständiges Orthonormalsystem. Mit Hilfe dieser Hilbertraumbasis soll nun eine Matrixdarstellung der Leiteroperatoren ermittelt werden. Man beachte, dass hier alle Indizes von 0 (nicht von 1) bis unendlich laufen. Die Eigenzustände lassen sich als Vektoren darstellen:

\left| 0 \right\rangle =\left( \begin{matrix}
   1  \\
   0  \\
   0  \\
   0  \\
   \vdots   \\
\end{matrix} \right),\quad \left| 1 \right\rangle =\left( \begin{matrix}
   0  \\
   1  \\
   0  \\
   0  \\
   \vdots   \\
\end{matrix} \right),\quad \left| 2 \right\rangle =\left( \begin{matrix}
   0  \\
   0  \\
   1  \\
   0  \\
   \vdots   \\
\end{matrix} \right)    usw.

Die Vollständigkeit dieser Basis liefert eine Darstellung des Einheitsoperators:

\sum\limits_{n=0}^{\infty }{\left| n \right\rangle \left\langle  n \right|}=1

Erzeugungsoperator[Bearbeiten]

Vor und nach dem Erzeugungsoperator wird eine 1 (Einheitsoperator) eingeschoben:

a^{\dagger }=\sum\limits_{m,n=0}^{\infty }{\left| m \right\rangle \underbrace{\left\langle  m \right|a^{\dagger }\left| n \right\rangle }_{a_{mn}^{\dagger }}\left\langle  n \right|}

Das Matrixelement berechnet sich zu

a_{mn}^{\dagger }=\left\langle  m \right|a^{\dagger }\left| n \right\rangle =\sqrt{n+1}\left\langle  m | n+1 \right\rangle =\sqrt{n+1}\ \delta _{m,n+1}

Der Erzeugungsoperator dargestellt durch die Basisvektoren

a^{\dagger }=\sum\limits_{n=0}^{\infty }{\left| n+1 \right\rangle \sqrt{n+1}\left\langle  n \right|}

Somit ergibt sich die Matrixdarstellung des Erzeugungsoperators bzgl. der Besetzungseigenbasis (alle nicht angegebenen Elemente sind gleich 0):

a^{\dagger }=\left( \begin{matrix}
   0 & {} & {} & {} & {}  \\
   \sqrt{1} & 0 & {} & {} & {}  \\
   {} & \sqrt{2} & 0 & {} & {}  \\
   {} & {} & \sqrt{3} & 0 & {}  \\
   {} & {} & {} & \ddots  & \ddots   \\
\end{matrix} \right)

Vernichtungsoperator[Bearbeiten]

Durch analoge Rechnung erhält man für den Vernichtungsoperator:

a=\sum\limits_{m,n=0}^{\infty }{\left| m \right\rangle \underbrace{\left\langle  m \right|a\left| n \right\rangle }_{a_{mn}}\left\langle  n \right|}=\sum\limits_{m,n=0}^{\infty }{\left| m \right\rangle \sqrt{n}\ \delta _{m+1,n}\left\langle  n \right|}=\sum\limits_{n=0}^{\infty }{\left| n \right\rangle \sqrt{n+1}\ \left\langle  n+1 \right|}

Dabei wurde das Matrixelement schon eingesetzt:

a_{mn}=\left\langle  m \right|a\left| n \right\rangle =\sqrt{n}\left\langle  m | n-1 \right\rangle =\sqrt{n}\ \delta _{m,n-1}=\sqrt{n}\ \delta _{m+1,n}

Matrixdarstellung des Vernichtungsoperators bzgl. der Besetzungseigenbasis:

a=\left( \begin{matrix}
   0 & \sqrt{1} & {} & {} & {}  \\
   {} & 0 & \sqrt{2} & {} & {}  \\
   {} & {} & 0 & \sqrt{3} & {}  \\
   {} & {} & {} & 0 & \ddots   \\
   {} & {} & {} & {} & \ddots   \\
\end{matrix} \right)

Man erkennt, dass die Matrix a^{\dagger} genau die transponierte von a ist. Dies ist verständlich, da die beiden Operatoren zueinander adjungiert (= transponiert + komplex konjugiert) sind.

Besetzungszahloperator[Bearbeiten]

Matrixelement des Besetzungszahloperators bzgl. der Besetzungseigenbasis:

N_{mn}=\left\langle  m \right|N\left| n \right\rangle =n\left\langle  m | n \right\rangle =n\ \delta _{m,n}

alternativ mit den Leiteroperatoren:

N_{mn}=\sum\limits_{k=0}^{\infty }{a_{mk}^{\dagger }a_{kn}}=\sum\limits_{k=0}^{\infty }{\underbrace{\sqrt{k+1}\ \delta _{m,k+1}}_{a_{mk}^{\dagger }}\underbrace{\sqrt{n}\ \delta _{k+1,n}}_{a_{kn}}}=n\ \delta _{m,n}

Matrixdarstellung des Besetzungszahloperators bzgl. der Besetzungseigenbasis:

N=\left( \begin{matrix}
   0 & {} & {} & {}  \\
   {} & 1 & {} & {}  \\
   {} & {} & 2 & {}  \\
   {} & {} & {} & \ddots   \\
\end{matrix} \right)

Hamiltonoperator des harmonischen Oszillators[Bearbeiten]

Matrixelement des Hamiltonoperators für den harmonischen Oszillator bzgl. der Besetzungseigenbasis bzw. der Energieeigenbasis:

H_{mn}=\left\langle  m \right|H\left| n \right\rangle =\left\langle  m \right|\hbar \omega \left( N+\frac{1}{2} \right)\left| n \right\rangle =\hbar \omega \left( n+\frac{1}{2} \right)\left\langle  m | n \right\rangle =\hbar \omega \left( n+\frac{1}{2} \right)\delta _{m,n}

Matrixdarstellung des Hamiltonoperators für den harmonischen Oszillator bzgl. der Besetzungseigenbasis bzw. der Energieeigenbasis:

H=\hbar \omega \left( \begin{matrix}
   \frac{1}{2} & {} & {} & {}  \\
   {} & 1+\frac{1}{2} & {} & {}  \\
   {} & {} & 2+\frac{1}{2} & {}  \\
   {} & {} & {} & \ddots   \\
\end{matrix} \right)

Da die Operatoren N und H hermitesch sind, folgt dass die zugehörigen Matrizen bzgl. der Eigenbasen symmetrisch sind.

Eigenzustände bosonischer Kletteroperatoren („kohärente Zustände“)[Bearbeiten]

Die Eigenzustände des Vernichtungsoperators sind die kohärenten Zustände \left| \alpha \right\rangle. Der Vernichtungsoperator \hat a  (zur Verdeutlichung sind die \hat{ }-Symbole für die Operatoren hier explizit wieder eingeführt) erfüllt folgende Eigenwertgleichung:

\hat a\left| \alpha  \right\rangle =\alpha \left| \alpha  \right\rangle

Für den Erzeugungsoperator ergibt sich daraus, mit einem Linkseigenzustand (Bra-Eigenzustand):

\left\langle  \alpha  \right|\hat a^{\dagger }=\alpha ^{*}\left\langle  \alpha  \right|

Der Vernichtungsoperator \hat a kann also - im Gegensatz zum Erzeugungsoperator \hat a^\dagger - Rechtseigenzustände (Ket-Eigenzustände) besitzen. Der Erzeugungsoperator erhöht die minimale Teilchenzahl eines Zustandes im Fockraum um eins; der damit entstandene Zustand kann also nicht der ursprüngliche sein. Dagegen verringert der Vernichtungsoperator die maximale Teilchenzahl um eins; da ein Zustand im Fockraum aber Komponenten aller Teilchenzahlen (auch beliebig hoher Teilchenzahlen) beinhalten kann, ist damit nicht verboten, dass \hat a Eigenzustände besitzt. Dies sind die kohärenten Zustände:

Der „kohärente Zustand“ |\alpha\rangle ergibt sich als bestimmte Linearkombination aller Zustände fester Teilchenzahl |n\rangle\,, und zwar nach der Formel:

\left| \alpha  \right\rangle :=e^{-\frac{\left| \alpha  \right|^{2}}{2}}\sum\limits_{n=0}^{\infty }{\frac{\alpha ^{n}}{\sqrt{n!}}}\left| n \right\rangle =e^{-\frac{\left| \alpha  \right|^{2}}{2}}\sum\limits_{n=0}^{\infty }{\frac{\alpha ^{n}}{\sqrt{n!}}}\frac{\left(\hat a^{\dagger } \right)^{n}}{\sqrt{n!}}\left| 0 \right\rangle =e^{-\frac{\left| \alpha  \right|^{2}}{2}}e^{\alpha \hat a^{\dagger }}\left| 0 \right\rangle

Dieser Zustand ist also Eigenzustand des Vernichtungsoperators, und zwar zum Eigenwert \alpha\,, während der zugehörige Erzeugungsoperator nur Links-Eigenzustände besitzt. Dabei ist \alpha eine nichtverschwindende komplexe Zahl, die den kohärenten Zustand vollständig definiert und auch explizit von der Zeit abhängen darf.  \left|\alpha \right|^2 ist der Erwartungswert der Besetzungszahl des kohärenten Zustandes.

Kohärente Zustände haben (wie der Grundzustand des Harmonischen Oszillators) minimale Unschärfe und bleiben bei Zeitentwicklung kohärent. Mit ihnen lässt sich die – im Allgemeinen explizit zeitabhängige – elektromagnetische Welle einer Laser-Mode am besten beschreiben (sog. Glauber-Zustände).

Erzeugungs- und Vernichtungsoperatoren in Quantenfeldtheorien[Bearbeiten]

In Quantenfeldtheorie und Vielteilchenphysik verwendet man Ausdrücke der Form \gamma_ia_i^\pm\,, wobei die \gamma_i komplexe Zahlen sind, während die a_i^\pm Erzeugungs- bzw. Vernichtungsoperatoren darstellen: Diese erhöhen bzw. vermindern die Eigenwerte des Anzahloperators \sum a_i^+a_i^- um 1, analog zum harmonischen Oszillator. Die Indizes i berücksichtigen die Freiheitsgrade der Raumzeit und haben i.a. mehrere Komponenten. Wenn die Erzeuger und Vernichter von einer kontinuierlichen Variable abhängen, statt von diskreten Quantenzahlen, schreibt man sie auch als Feldoperatoren \phi(\vec{x}), \phi^+(\vec{x}). Die Anzahloperatoren n_i=a_i^+a_i^- sind selbstadjungiert („hermitesch“) und nehmen alle nicht-negativen ganzzahligen Werte an: n_i\in\{0,\,\,1,\,\,2,\,\,\dots \}\,.
Die nichttrivialen Vertauschungsrelationen sind schließlich, wie beim harmonischen Oszillator:  [a_i^-,a_j^+]:=a_i^-a_j^+-a_j^+a_i^-=\delta_{ij}\,, wobei [.,.] die sog. Kommutatorklammer darstellt, während \delta_{ij} das Kroneckersymbol ist.

Das oben gesagte gilt für Bosonen, wogegen man für Fermionen den Kommutator durch den Antikommutator ersetzen muss, \{a_i^-,a_j^+\}:=a_i^-a_j^++a_j^+a_i^- =\delta_{ij}\,. Als Konsequenz gilt im fermionischen Fall, dass die Anzahloperatoren n_i=a_i^+a_i^- nur die Eigenwerte 0 und 1 haben.

Bezug zu Diagrammtechniken[Bearbeiten]

Konkrete Rechnungen unter Verwendung der Erzeugungs- bzw. Vernichtungsoperatoren kann man in der Regel durch Diagrammtechniken unterstützen (→ Feynman-Diagramme). So kann man z. B. Dreiteilchen-Wechselwirkungen der Form \gamma_{1,2;3}\,a_1^-a_2^-a_3^+ durch drei Linien veranschaulichen, von denen die ersten zwei in einen Vertex einlaufen und dort „vernichtet“ werden, während eine dritte Linie an diesem Vertex „erzeugt“ wird und von ihm ausläuft. Dabei sind in den zugehörigen Regeln Energie- und Impulssatz explizit zu berücksichtigen.

Der angegebene Term, der einen sog. „Konfluenzprozess“ beschreibt, hat bei tiefen Temperaturen i. A. geringere Wahrscheinlichkeit, \propto |\gamma_{1,\,2;\,3}|^2\cdot\{\langle n_1\rangle\,\langle n_2\rangle\,(1+\langle n_3\rangle \}\,, als der inverse sog. „Splitting-Prozess“, der zum adjungierten Term, \gamma_{1,2;3}^*\,a_3^-a_2^+a_1^+\,, gehört. Denn zu jedem Erzeugungsoperator a_i^+ korrespondiert, analog zum harmonischen Oszillator, die Übergangsrate[1] \propto |\gamma^*_{\,i,\dots} |^2\cdot (1+\langle n_i\rangle ), während beim zugehörigen Vernichtungsoperator  der Term \propto 1 fehlt. Auf diese Weise sind bei tiefen Temperaturen die letztgenannten Terme in der Regel wichtiger als die erstgenannten.[2]

Literatur[Bearbeiten]

  • Cohen-Tannoudji, Diu, Laloë: Quantenmechanik 1/2. de Gruyter, Berlin
  • Nolting: Grundkurs theoretische Physik. Bd.5/1 : Quantenmechanik. Springer, Berlin

Siehe auch[Bearbeiten]

Belege[Bearbeiten]

  1. Siehe das Kapitel über zeitabhängige Störungsrechnung in den meisten Standardlehrbüchern der Quantenmechanik II
  2. Siehe z. B. Hermann Haken: Quantenfeldtheorie des Festkörpers, Teubner 1973, ISBN 3-519-03025-X