Teilchenzahloperator

Die Erzeugungs- und Vernichtungsoperatoren sind der Kern einer eleganten Lösung der Schrödingergleichung des harmonischen Oszillators. Diese Operatoren können auch dazu benutzt werden, gewisse 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 σ ₊ und σ ₋ , 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. Man nennt sie deswegen Vernichtungs- und Erzeugungsoperator

Das Zirkumflex ü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 zugunsten der Übersichtlichkeit verzichtet. Alle lateinischen Großbuchstaben, mit Ausnahme des E, sind Operatoren.

Definition

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 Antikommutatior-Vertauschungsrelation erfüllen.

Details

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 lassen sich ausgehend von einem Eigenzustand durch mehrmaliges 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

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

• 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 φ 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

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

• Mit $a\, a=0$ und $a^{\dagger}a^{\dagger}=0$ ist N² = 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 δ_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

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, ω 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

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-i\frac{1}{m\omega }P \right)$

$a=\sqrt{\frac{m\omega }{2\hbar }}\left( Q+i\frac{1}{m\omega }P \right)$

Mit den Leiteroperatoren schreibt sich der Hamiltonoperator:

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

Eigenschaften der Erzeuger und Vernichter

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$

Der Besetzungszahloperator N ist ein hermitescher Operator und hat daher reelle Eigenwerte, die Besetzungszahlen n. Die Eigenzustände sind Fock-Zustände oder Besetzungszahlzustände. Die zugehörige Eigenwertgleichung lautet:

$N\left| n \right\rangle =n\left| n \right\rangle \ ,\quad n\in\mathbb{N}_{0}$

Lösung des Eigenwertproblems

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 (Nullpunktsenergie):

$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 Nullpunktsenergie $\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, wieviele Energiequanten in einem Eigenzustand angeregt sind. Die Besetzung eines Zustandes mit n Energiequanten erklärt den Namen Besetzungszahloperator.

Eigenfunktionen in Ortsdarstellung

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

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

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

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

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

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

Die Eigenzustände des Vernichtungsoperators sind die kohärenten Zustände $\left| \alpha \right\rangle$. Der Vernichtungsoperator a erfüllt folgende Eigenwertgleichung:

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

Der Erzeugungsoperator erfüllt folgende Eigenwertgleichung, mit einem Linkseigenzustand (Bra-Eigenzustand):

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

Der Vernichtungsoperator a kann im Gegensatz zum Erzeugungsoperator $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 a Eigenzustände besitzt. Dies sind die kohärenten Zustände.

Der kohärente Zustand $|\alpha\rangle$ ergibt sich als Linearkombination von allen Zuständen fester Teilchenzahl $|n\rangle$ nach:

$\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( a^{\dagger } \right)^{n}}{\sqrt{n!}}\left| 0 \right\rangle =e^{-\frac{\left| \alpha \right|^{2}}{2}}e^{\alpha a^{\dagger }}\left| 0 \right\rangle$

Dabei ist α eine nichtverschwindende komplexe Zahl, die den kohärenten Zustand vollständig definiert. $\left|\alpha \right|^2$ ist der Erwartungswert der Besetzungszahl des kohärenten Zustandes. Kohärente Zustände haben minimale Unschärfe und bleiben bei Zeitentwicklung kohärent. Mit ihnen lässt sich die elektromagnetische Welle einer Laser-Mode am besten beschreiben.

Literatur

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

Siehe auch

• Fockraum