Quantenmechanik/ Harmonischer Oszillator

Der Harmonische Oszillator

Betrachten wir den Hamiltonoperator in einer Dimension:

$\hat{H}=\frac{p^2}{2m}+\frac{m{\omega }^2}{2}{x}^2$

und suchen wir die dazugehörenden Eigenfunktionen und Eigenwerte.

Es gibt zwei grundsätzliche Lösungswege. Zum einen könnten wir den Operator in die Ortsdarstellung der Schrödingergleichung einsetzen und die resultierende Differentialgleichung mit den entsprechenden Randbedingungen lösen. Der zweite Weg zeigt eine grundsätzliche Herangehensweise an Probleme in der Quantenmechanik auf und nutzt die Operatordarstellung. Dieser Lösungsweg wird hier gerechnet.

Zerlegung des Hamiltonoperators

Zunächst führen wir drei Operatoren ein:

$a=\sqrt{\frac{m\omega }{2\hslash }}(x+\frac{i}{m\omega }p)$

$a^{†}=\sqrt{\frac{m\omega }{2\hslash }}(x-\frac{i}{m\omega }p)$

$N=a{a}^{†}$

Zu beachten ist, dass $a,a^{†}$ keine hermiteschen Operatoren sind, $N$ allerdings schon. Man nennt Operatoren, die in ähnlicher Weise wirken wie $a,a^{†}$ Vernichtungs- und Erzeugungsoperatoren, bzw. auch Ab- und Aufsteigeoperatoren. Warum und worauf auf- und abgestiegen wird, werden wir im nächsten Kapitel sehen.

Zuvor müssen wir uns aber einige Eigenschaften klarmachen

1. $[a,a^{†}]=1$

   Beweis: $[a,a^{†}]=\frac{m\omega }{2\hslash }[x+\frac{i}{m\omega }p,x-\frac{i}{m\omega }p]=\frac{m\omega }{2\hslash }(-\frac{i}{m\omega }[x,p]+\frac{i}{m\omega }[p,x])=\frac{m\omega }{2\hslash }\frac{-i}{m\omega }2i\hslash =1$

Der Hamiltonoperator lässt sich nun darstellen als $\hat{H}=\hslash \omega (a{a}^{†}-\frac{1}{2})$. Wir können jetzt schon das Energiespektrum dieses Systems bestimmen! Nehmen wir an den Grundzustand $|0⟩$ mit Energie $E_0$ zu kennen, das heißt

$\hat{H}|0⟩=\hslash \omega (a{a}^{†}-\frac{1}{2})|0⟩=E_0|0⟩$.

Wenden wir nun den Operator $a^{†}$ auf diese Gleichung an und benutzen den Kommutator von oben, so erhalten wir

$a^{†}\hat{H}|0⟩=\hslash \omega (a^{†}a{a}^{†}-\frac{1}{2}{a}^{†})|0⟩=\hslash \omega (a{a}^{†}{a}^{†}-a^{†}-\frac{1}{2}{a}^{†})|0⟩=(\hat{H}-\hslash \omega )a^{†}|0⟩=E_0{a}^{†}|0⟩$.

Wir definieren nun den Zustand $|1⟩=a^{†}|0⟩$, und bringen ein $\hslash \omega$ auf die andere Seite. Damit erhalten wir

$\hat{H}|1⟩=(E_0+\hslash \omega )|1⟩$,

was die Schrödingergleichung für den ersten angeregten Zustand ist. Dieser hat also die Energie $E_0+\hslash \omega$. Da das gleiche Argument wieder auf den ersten Zustand angewandt werden kann, ist die Energie für den n-ten Zustand:

$E_n=E_0+n\cdot \hslash \omega$

Wir haben also gesehen, dass $a^{†}$ der Operator ist, der einen Zustand in den nächst höheren übergehen lässt, also $a^{†}|n⟩=|n+1⟩$. Der Operator $a$ tut das gleiche in umgekehrter Richtung, er senkt die Energie ab. Daher heißen diese Operatoren Erzeuger/Vernichter, sie erzeugen/vernichten ein Energiequant.

Nun können wir noch die Energie $E_0$ des Grundzustands $|0⟩$ bestimmen. Da es per Definition keinen Zustand unterhalb des Grundzustands gibt, gilt $a|0⟩=0$. Nutzen wir wieder die Kommutatorbeziehung von oben aus, erhalten wir

$H|0⟩=\hslash \omega (a{a}^{†}-\frac{1}{2})|0⟩=\hslash \omega (a^{†}a+\frac{1}{2})|0⟩=\frac{1}{2}\hslash \omega |0⟩$.

Die Wellenfunktion in Ortsdarstellung

Die ursprüngliche Differentialgleichung lässt sich zwar nicht ohne weiteres lösen, mit dem oben hergeleiteten Formalismus lässt sich aber relativ leicht ein Konstruktionsschema für die Wellenfunktion in Ortsdarstellung angeben.

Dazu nutzen wir die Eigenschaft des Grundzustands aus, dass sich kein Energiequant vernichten lässt:

$a|0⟩=0.$

Der Einfachheit zuliebe betrachten wir nur eine Raumdimension, der Impulsoperator wird also im Ortsraum $p=\frac{\hslash }{i}\frac{d}{dx}$. Setzen wir dies also in die Definition von $a$ ein ergibt sich die Differentialgleichung

$(x+\frac{\hslash }{m\omega }\frac{d}{dx}){\psi }_0(x)=0.$

Diese Gleichung lässt sich durch Separation der Variablen lösen:

$\frac{d{\psi }_0}{{\psi }_0}=-\frac{m\omega }{\hslash }xdx,$

was nach Integration die Wellenfunktion des Grundzustand liefert:

${\psi }_0(x)=A\exp (-\frac{m\omega }{2\hslash }{x}^2).$

Der Faktor A entsteht als Integrationskonstante und dient der Normierung der Wellenfunktion. Um ihn zu bestimmen wird das uneigentliche Integral $\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}dx{e}^{-a{x}^2}=\sqrt{\frac{\pi }{a}}$ benötigt:

$\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}dx|{\psi }_0(x){|}^2=|A{|}^2\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}dx\exp (-\frac{m\omega }{\hslash }{x}^2)=|A{|}^2\sqrt{\frac{\pi \hslash }{m\omega }}\stackrel{!}{=}1,$

womit, bis auf eine komplexe Phase, A bestimmt ist als $A={\left(\frac{m\omega }{\pi \hslash }\right)}^{1/4}$.

Vom Grundzustand aus kann man nun schrittweise die Wellenfunktionen der angeregten Zustände berechnen. Dazu muss nur der Operator $a^{†}$ n-mal auf ${\psi }_0$ angewandt werden, um den n-ten Zustand zu erreichen. Zum Beispiel entsteht so der erste angeregte Zustand:

${\psi }_1(x)=a^{†}{\psi }_0(x)=\sqrt{\frac{m\omega }{2\hslash }}(x-\frac{\hslash }{m\omega }\frac{d}{dx}){\psi }_0(x)=2{\left(\frac{m\omega }{\pi \hslash }\right)}^{1/4}\sqrt{\frac{m\omega }{2\hslash }}x\exp (-\frac{m\omega }{2\hslash }{x}^2)$

Eigenwerte des quantenmechanischen harmonischen Oszillators (Rechnung)

In diesem Artikel wird gezeigt, wie die Energie-Eigenwerte des quantenmechanischen harmonischen Oszillators in einer Dimension berechnet werden können. Die Rechnung versteht sich auch als ein repräsentatives Beispiel für die Lösung quantenmechanischer Probleme mit Hilfe der Methode der zweiten Quantisierung, etwa die Bestimmung von Drehimpulseigenwerten oder in der Festkörperphysik.

Problemstellung und Vorgehensweise

Allgemeines

Wir bestimmen die Energie-Eigenwerte des eindimensionalen quantenmechanischen harmonischen Oszillators, also sämtliche Eigenwerte $ϵ$ des Eigenwertproblems

$H|\psi ⟩=ϵ|\psi ⟩,$

wobei der Hamilton-Operator $H$ gegeben ist durch $H=\frac{1}{2m}{P}^2+\frac{1}{2}m{\omega }^2{X}^2,$ $P$ ist der Impulsoperator, $X$ der Ortsoperator. Dabei werden wir wie folgt vorgehen:

Einführung von Auf- und Absteigern

Zunächst werden wir zueinander adjungierte lineare Operatoren $a$ und $a^{+}$ einführen, mit denen sich der Hamilton-Operator $H$ schreiben lässt als

$H=\hslash \omega (a^{+}a+\frac{1}{2}).$

Untersuchung des Operators $N:=a^{+}a$

Anstelle der Eigenzustände von $H$ untersuchen wir die Eigenzustände von $N:=a^{+}a$.Es lässt sich leicht zeigen, dass wir alle Eigenzustände von $H$ aus denen von $N$ konstruieren können, ebenso können wir die Eigenwerte von $H$ leicht aus denen von $N$ gewinnen. Wir konzentrieren uns daher auf das Eigenwertproblem

$N|{\psi }_n^j⟩=n|{\psi }_n^j⟩.$

Wir haben dabei die Eigenzustände von $N$ mit zwei Indizes "nummeriert": Der untere Index gibt den zu $|{\psi }_n^j⟩$ gehörigen Eigenwert an und durchläuft die (noch unbekannte) Menge aller Eigenwerte. Auch wenn die Wahl des Buchstabens bereits die spätere Erkenntnis andeutet, dass es sich bei den Eigenwerten von $N$ um natürliche Zahlen handelt, so ist an dieser Stelle noch keinerlei Einschränkung vorgenommen; die Menge der Eigenwerte und damit der Index $n$ kann durchaus eine kontinuierliche Größe sein. Gleiches gilt für den oberen Index $j$, der dazu dient, Eigenfunktionen zum gleichen Eigenwert zu unterscheiden, ohne ein neues Symbol verwenden zu müssen. Da der Eigenzustand $|{\psi }_n^j⟩$ durch die Angabe der beiden Indizes $n$ und $j$ bereits vollständig bestimmt ist, ist es in der Literatur vielfach üblich, das Symbol $\psi$ komplett wegzulassen und den Zustand einfach nur durch

$|{\psi }_n^j⟩=:|n,j⟩$

zu kennzeichnen. Dieser Konvention werden wir uns aber nicht anschließen.

Bestimmung der Eigenwerte von $N$

Um die Eigenwerte von $N$ zu bestimmen, werden wir die Zustände $a|{\psi }_n^j⟩$ und $a^{+}|{\psi }_n^j⟩$ genauer untersuchen: wir werden sehen, dass $a|{\psi }_n^j⟩$ unter bestimmten Bedingungen ein Eigenzustand zum Eigenwert $n-1$ und $a^{+}|{\psi }_n^j⟩$ ein Eigenzustand zum Eigenwert $n+1$ ist, darüber hinaus werden wir zeigen, dass $n=0$ ein Eigenwert von $N$ ist, indem wir in der Ortsdarstellung eine konkrete Wellenfunktion des zu $0$ gehörigen Eigenvektors $|{\psi }_0⟩$ ermitteln; Da, wie wir sehen werden, alle Eigenwerte von $N$ nicht-negativ sind, ist dann schon bewiesen, dass alle $n\in {\mathbb{N}}_0$ zum Eigenwertspektrum von $N$ gehören:

1. 0 ist Eigenwert, wie wir durch explizite Angabe eines Eigenzustandes $|{\psi }_0⟩$ zu diesem Eigenwert sehen werden.
2. $|{\psi }_1⟩:=a^{+}|{\psi }_0⟩$ ist ein Eigenzustand zum Eigenwert $0+1=1$, also ist auch die Zahl 1 aus dem Eigenwertspektrum von $N$. Ebenso ist $a^{+}|{\psi }_1⟩$ ein Eigenzustand zum Eigenwert $1+1=2$, also ist auch 2 Eigenwert von $N$.

Sukzessive folgt so, dass alle natürlichen Zahlen zum Eigenwertspektrum von $N$ gehören (formal vollziehen wir den Beweis mit Hilfe der vollständigen Induktion). Induktiv werden wir abschließend zeigen, dass $N$ keine weiteren Eigenwerte besitzt. Damit ist das Eigenwertspektrum von $N$ und somit auch das von $H$ vollständig bestimmt. Es sei angemerkt, dass wir die Rechnung in demjenigen Hilbertraum $\mathscr{H}$ durchführen, der der Bewegung eines spinlosen Teilchens in einer Dimension zugeordnet ist. Wir werden daher die Begriffe Eigenzustand und Eigenvektor synonym verwenden.

Die Rechnung im Detail

Schritt 1: Umschreiben des Hamiltonoperators

Zunächst führen wir die Operatoren $a$ und $a^{+}$ ein. Aus Gründen, die später klar werden, nennen wir $a^{+}$ den Aufsteigeoperator, $a$ den Absteigeoperator (übliche Bezeichnungen sind auch Erzeuger und Vernichter.). Sie sind definiert durch:

$\begin{array}{ccc}a& :=& \sqrt{\frac{m\omega }{2\hslash }}X+i\frac{1}{\sqrt{2m\hslash \omega }}P\\ a^{+}& :=& \sqrt{\frac{m\omega }{2\hslash }}X-i\frac{1}{\sqrt{2m\hslash \omega }}P.\end{array}$

Sowohl $a$ als auch $a^{+}$ sind lineare Operatoren, da $P$ und $X$ linear sind. Man beachte in allen nachfolgenden Rechnungen, dass der Einsoperator, definiert durch

$\mathrm{𝟏}|\psi ⟩=|\psi ⟩\mathrm{\forall }|\psi ⟩\in \mathscr{H},$

stets unterdrückt wird: so schreiben wir z. B. stets $[X,P]=i\hslash$ anstatt $[X,P]=i\hslash \mathrm{𝟏}$. Diese Unterdrückung gestaltet Rechnungen übersichtlicher und ist an gängige Konventionen in üblicher Quantenmechaniklehrbücher angepasst. Die Operatornatur des Kommutators $[\cdot ,\cdot ]$ sollte dabei aber nicht aus den Augen verloren werden! Widmen wir uns aber nun der Bestimmung der Eigenwerte von $H$. Wie bereits in der Einleitung angedeutet, zeigen wir zunächst:

Satz 1: Der Hamiltonoperator ist gegeben durch $H=\hslash \omega (a^{+}a+\frac{1}{2})$.

Beweis: Den Beweis führen wir direkt durch Einsetzen, wobei wir ausnutzen, dass die Menge der linearen Operatoren auf $H$ einen (nichtkommutativen) Ring bildet:

$\begin{array}{ccc}{a}^{+}a& =& (\sqrt{\frac{m\omega }{2\hslash }}X-i\frac{1}{\sqrt{2m\hslash \omega }}P)\cdot (\sqrt{\frac{m\omega }{2\hslash }}X+i\frac{1}{\sqrt{2m\hslash \omega }}P)\\ & =& \frac{m\omega }{2\hslash }{X}^2+i\frac{1}{2\hslash }XP-i\frac{1}{2\hslash }PX+\frac{1}{2m\omega \hslash }{P}^2\\ & =& \frac{m\omega }{2\hslash }{X}^2+\frac{i}{2\hslash }[X,P]+\frac{1}{2m\hslash \omega }{P}^2\end{array}$

Bekanntlich ist $[X,P]=i\hslash$, so dass insgesamt folgt:

$a^{+}a=\frac{m\omega }{2\hslash }{X}^2-\frac{1}{2}+\frac{1}{2m\hslash \omega }{P}^2$

oder

$\hslash \omega a^{+}a=\underset{=H}{\underbrace{\frac{1}{2}m{\omega }^2{X}^2+\frac{1}{2m}{P}^2}}-\frac{1}{2}\hslash \omega$

und somit ist alles gezeigt.

Wir führen nun das Eigenwertspektrum des Hamiltonoperators auf das Eigenwertspektrum des Operators $N:=a^{+}a$ zurück:

Satz 2: Jeder Eigenvektor von $H$ ist auch Eigenvektor von $N$ und umgekehrt. Weiterhin gilt: Alle Eigenwerte $ϵ$ von $H$ sind durch $ϵ=(n+\frac{1}{2})\hslash \omega$ gegeben, wenn $n$ die Eigenwerte von $N$ durchläuft.

Beweis:: Sei $|\psi ⟩$ ein Eigenvektor von $N$ zum Eigenwert $n$. Dann gilt:

$\begin{array}{ccc}H|\psi ⟩& =& \hslash \omega (N+\frac{1}{2})|\psi ⟩\\ & =& \hslash \omega N|\psi ⟩+\frac{\hslash \omega }{2}|\psi ⟩\\ & =& \hslash \omega n|\psi ⟩+\frac{\hslash \omega }{2}|\psi ⟩\\ & =& \hslash \omega (n+\frac{1}{2})|\psi ⟩,\end{array}$

$|\psi ⟩$ ist also ein Eigenvektor von $H$, wie behauptet. Aus der Rechnung folgt weiterhin direkt, dass jede Zahl der Form $(n+\frac{1}{2})\hslash \omega$ mit einem Eigenwert $n$ von $N$ ein Eigenwert des Hamiltonoperators $H$ ist. Sei nun umgekehrt $|\varphi ⟩$ ein Eigenvektor von $H$ zum Eigenwert $ϵ$. Wegen $N=\frac{1}{\hslash \omega }H-\frac{1}{2}$ folgt mit einer analogen Rechnung, dass $|\varphi ⟩$ auch Eigenvektor von $N$ ist. $H$ und $N$ besitzen also die gleichen Eigenzustände, da jeder Eigenzustand von $N$ auch Eigenzustand von $H$ ist und umgekehrt. Bleibt noch zu zeigen, dass jeder Eigenwert von $H$ von der Form $\hslash \omega (n+\frac{1}{2})$ mit einem Eigenwert $n$ von $N$ ist. Sei dazu $ϵ$ ein Eigenwert von $H$. Dann existiert dazu eine Eigenfunktion $|{\psi }_{ϵ}⟩$. Nun ist $|{\psi }_{ϵ}⟩$ aber auch Eigenfunktion von $N$, wie weiter oben gezeigt. Damit gilt:

$ϵ=H|{\psi }_{ϵ}⟩=\hslash \omega (N+\frac{1}{2})|{\psi }_{ϵ}⟩=\hslash \omega (n+\frac{1}{2})|{\psi }_{ϵ}⟩,$

da $|{\psi }_{ϵ}⟩\ne 0$ folgt $ϵ=\hslash \omega (n+\frac{1}{2})$, wie behauptet.

Der letzte Satz garantiert uns, dass wir sämtliche Eigenvektoren sowie das gesamte Eigenwertspektrum von $H$ aus den Eigenvektoren bzw. Eigenwerten von $N$ rekonstruieren können. Wenn wir alle Eigenvektoren und -Werte von $N$ kennen, ist das Problem der Bestimmung des Eigenwertspektrums des eindimensionalen harmonischen Oszilators im Prinzip gelöst. Daher werden wir nun das Eigenwertproblem

$N|{\psi }_n^j⟩=n|{\psi }_n^j⟩$

genauer untersuchen, wobei wir konsequent die bereits in der Einführung beschriebene Notation verwenden werden: der Zustand $|{\psi }_n^j⟩$ steht in allen nachfolgenden Rechnungen stets für einen Eigenzustand von $N$ zum Eigenwert $n$, auch wenn dies keine explizite Erwähnung findet! Per Definition eines Eigenzustandes ist damit auch $|{\psi }_n^j⟩\ne 0$. Wir zeigen zunächst die folgende wichtige Tatsache:

Satz 3: $a^{+}$ ist adjungiert zu $a$:$a^{+}=a^{†}$

Beweis: Dies erkennt man sofort aus den Regeln zur Bildung des adjungierten Operators zu einem Operator, der Linearkombination anderer linearer Operatoren ist:

$\begin{array}{ccc}{a}^{†}& =& {(\sqrt{\frac{m\omega }{2\hslash }}X+i\frac{1}{\sqrt{2m\hslash \omega }}P)}^{†}\\ & =& \sqrt{\frac{m\omega }{2\hslash }}{X}^{†}-i\frac{1}{\sqrt{2m\hslash \omega }}{P}^{†}\end{array}$

Da $X$ und $P$ selbstadjungiert sind (d.h. $X^{†}=X,P^{†}=P$), ist die Behauptung gezeigt.

Schritt 2: Bestimmung der Eigenwerte von $N$

Aus Satz 3 folgt sofort:

Satz 4: Alle Eigenwerte von $N$ sind reell und nicht negativ.

Beweis: $n$ ist reell, da der Operator $N=a^{+}a$ selbstadjungiert ist:

$N^{†}=(a^{+}a{)}^{†}=a^{†}(a^{+}{)}^{†}=a^{†}(a^{†}{)}^{†}=a^{†}a=a^{+}a=N$

Nun berechnen wir die Norm des Zustandes $a|{\psi }_n^j⟩$

$\begin{array}{ccc}\Vert a|{\psi }_n^j⟩{\Vert }^2& =& ⟨a|{\psi }_n^j⟩,a|{\psi }_n^j⟩⟩\\ & =& ⟨|{\psi }_n^j⟩,a^{+}a|{\psi }_n^j⟩⟩\\ & =& ⟨{\psi }_n^j|a^{+}a|{\psi }_n^j⟩\\ & =& ⟨{\psi }_n^j|N|{\psi }_n^j⟩\\ & =& n⟨{\psi }_n^j|{\psi }_n^j⟩\end{array}$

Man beachte, dass wir im zweiten Schritt benutzt haben, dass $a^{+}=a^{†}$, im vorletzten Schritt haben wir $N|{\psi }_n^j⟩=n|{\psi }_n^j⟩$ verwendet. Wegen $⟨{\psi }_n^j|{\psi }_n^j⟩\ge 0$ und $\Vert a|{\psi }_n^j⟩{\Vert }^2\ge 0$ folgt sofort $n\ge 0$.

Wir untersuchen nun den Zustand $a|{\psi }_n^j⟩$ etwas genauer. Da wir in den nachfolgenden Rechnungen häufiger die Kommutatoren $[a,N]$, $[a^{+},N]$ und $[a,a^{+}]$ benötigen werden, wollen wir diese Ausdrücke an dieser Stelle berechnen. Zunächst rechnet man mit Hilfe der Definition von $a$ und $a^{+}$ unter Benutzung der Linearität des Kommutators und dem bekannten Postulat $[X,P]=i\hslash$ sofort nach, dass gilt:

$[a,a^{+}]=1.$

Damit lässt sich der Kommutator von $a$ und $N$ direkt ermitteln:

$\begin{array}{ccc}[a,N]& =& [a,a^{+}a]=a{a}^{+}a-a^{+}aa\\ & =& (a{a}^{+}-a^{+}a)a\\ & =& [a,a^{+}]a\end{array}$

Es folgt $[a,N]=a.$ Ganz analog lässt sich beweisen, dass

$[a^{+},N]=-a^{+}$

gilt. Wir haben nun alle nötigen Vorbereitungen für die nächsten beiden Sätze und für einen großen Schritt in Richtung der Lösung unseres Problems getroffen.

Satz 5: Existiert zu $n=0$ ein Eigenzustand $|{\varphi }_0⟩$, so ist $a|{\varphi }_0⟩=0$. Ist $n\ne 0$, so ist der Zustand $a|{\psi }_n^j⟩$ ein Eigenzustand zum Eigenwert $n-1$.

Beweis: Die erste Behauptung folgt direkt aus dem Beweis zum vorangegangenen Satz. Ist nämlich $n=0$, so gilt: $\Vert a|{\psi }_n^j⟩{\Vert }^2=n⟨{\psi }_n^j|{\psi }_n^j⟩=0,$ der Nullvektor ist aber der einzige Vektor mit einer Norm von 0, also ist $a|{\psi }_n^j⟩=0$. Sei nun $n\ne 0$. Um zu zeigen, dass $a|{\psi }_n^j⟩$ Eigenzustand zum Eigenwert $n-1$ ist, müssen wir beweisen:

• $a|{\psi }_n^j⟩$ genügt der Eigenwertgleichung $N(a|{\psi }_n^j⟩)=(n-1)(a|{\psi }_n^j⟩)$
• Es ist $a|{\psi }_n^j⟩\ne 0$.

Zum Beweis der ersten Behauptung untersuchen wir den Ausdruck $N(a|{\psi }_n^j⟩)$. Da $|{\psi }_n^j⟩$ Eigenzustand von $N$ ist, wäre es günstig, $a$ mit $N$ vertauschen zu können. Dies ist zwar nicht möglich, aber wir können die zuvor berechneten Kommutatoren ausnutzen, denn es gilt:

$[a,N]=aN-Na⟹Na=aN-[a,N].$

Der Kommutator $[a,N]$ ist aber bekannt: wie wir zuvor gezeigt haben, ist $[a,N]=a$, und es folgt:

$\begin{array}{ccc}N(a|{\psi }_n^j⟩)& =& (Na)|{\psi }_n^j⟩=(aN-[a,N])|{\psi }_n^j⟩\\ & =& (aN-a)|{\psi }_n^j⟩\\ & =& (aN)|{\psi }_n^j⟩-a|{\psi }_n^j⟩\\ & =& aN|{\psi }_n^j⟩-a|{\psi }_n^j⟩\\ & =& n\cdot a|{\psi }_n^j⟩-a|{\psi }_n^j⟩\\ & =& (n-1)(a|{\psi }_n^j⟩),\end{array}$

wie behauptet. $\Vert a|{\psi }_n^j⟩\Vert \ne 0$ folgt aus dem Beweis von Satz 4. Dort wurde gezeigt, dass die Norm von $a|{\psi }_n^j⟩$ durch

$\Vert a|{\psi }_n^j⟩\Vert =n⟨{\psi }_n^j|{\psi }_n^j⟩$

gegeben ist. Da $|{\psi }_n^j⟩$ als Eigenvektor per Definition nicht der Nullvektor ist und damit auch eine von 0 verschiedene Norm aufweist, ist auch $n⟨{\psi }_n^j|{\psi }_n^j⟩$ und damit $\Vert |{\psi }_n^j⟩\Vert$ von 0 verschieden, es gilt also $|{\psi }_n^j⟩\ne 0$.

Der letzte Satz gibt Aufschluss darüber, warum der Operator $a$ auch Absteigeoperator genannt wird: $a$ erzeugt aus einem Eigenvektor zum Eigenwert $n$ einen Eigenvektor zu dem um 1 verminderten Eigenwert $n-1$. Eine ganz analoge Eigenschaft besitzt auch der Aufsteigeoperator $a^{+}$: Dieser erzeugt aus einem Eigenvektor zum Eigenwert $n$ einen Eigenvektor zu dem um 1 erhöhten Eigenwert $n+1$, wie der nachfolgende Satz und sein Beweis zeigen! Man beachte, dass wir allerdings im Falle des Aufsteigeoperators die Einschränkung $n\ne 0$ fallen lassen können.

Satz 6: Der Zustand $a^{+}|{\psi }_n^j⟩$ ist Eigenzustand zum Eigenwert $n+1$.

Beweis: Der Beweis verläuft vollkommen analog zum Beweis von Satz 5. Wir bestimmen zunächst die Norm von $a^{+}|{\psi }_n^j⟩$. Wenn wir berücksichtigen, dass gilt:

$[a,a^{+}]=1⟹a{a}^{+}-a^{+}a=1⟹a{a}^{+}=1+a^{+}a,$

so finden wir:

$\begin{array}{ccc}\Vert a^{+}|{\psi }_n^j⟩{\Vert }^2& =& ⟨{\psi }_n^j|a{a}^{+}|{\psi }_n^j⟩\\ & =& ⟨{\psi }_n^j|1+a^{+}a||{\psi }_n^j⟩\\ & =& ⟨{\psi }_n^j|{\psi }_n^j⟩+⟨{\psi }_n^{j}N|{\psi }_n^j⟩\\ & =& ⟨{\psi }_n^j|{\psi }_n^j⟩+n⟨{\psi }_n^j|{\psi }_n^j⟩\\ & =& (n+1)⟨{\psi }_n^j|{\psi }_n^j⟩.\end{array}$

Da $n\ge 0$ gilt, ist $n+1>0$ für alle Eigenwerte $n$ von $N$. Da $|{\psi }_n^j⟩$ ein Eigenvektor ist, ist seine Norm stets von 0 verschieden, daher ist auch $\Vert a^{+}|{\psi }_n^j⟩\Vert \ne 0$, folglich ist $a^{+}|{\psi }_n^j⟩$ nicht der Nullvektor. Nun zeigen wir, dass $a^{+}|{\psi }_n^j⟩$ die Eigenwertgleichung

$N(a^{+}|{\psi }_n^j⟩)=(n+1)(a^{+}|{\psi }_n^j⟩)$

erfüllt:

$\begin{array}{ccc}N(a^{+}|{\psi }_n^j⟩)& =& (N{a}^{+})|{\psi }_n^j⟩=(a^{+}+a^{+}N)|{\psi }_n^j⟩\\ & =& a^{+}|{\psi }_n^j⟩+a^{+}(N|{\psi }_n^j⟩)=a^{+}|{\psi }_n^j⟩+n{a}^{+}|{\psi }_n^j⟩\\ & =& (n+1)(a^{+}|{\psi }_n^j⟩).\end{array}$

Bei der Rechnung haben wir verwendet, dass

$[a^{+},N]=-a^{+}⟹a^{+}N-N{a}^{+}=-a^{+}⟹N{a}^{+}=a^{+}N+a^{+}.$

Damit ist alles gezeigt.

Die immense Bedeutung des letzten Satzes ergibt sich anhand folgender Überlegung: können wir zeigen, dass die Zahl 0 ein Eigenwert des Operators $N$ darstellt, so können wir mit Hilfe des letzten Satzes beweisen, dass alle natürlichen Zahlen im Eigenwertspektrum von $N$ enthalten sind. Ist nämlich 0 ein Eigenwert, dann auch 1, schließlich ist $|{\psi }_1⟩:=a^{+}|{\psi }_0⟩$ nach Satz 6 ein Eigenzustand zum Eigenwert 1. Mit der gleichen Argumentation ist dann auch $|{\psi }_2⟩:=a^{+}|{\psi }_1⟩$ ein Eigenzzustand zum Eigenwert 2 usw., oder allgemeiner: ist irgendeine beliebige natürliche Zahl $n$ Eigenwert von $N$, so ist auch $n+1$ Eigenwert von $N$, wir müssen ja lediglich $a^{+}$ auf $|{\psi }_n^j⟩$ anwenden, um einen Eigenzustand zu $n+1$ zu erhalten. Der letzte Satz liefert uns also letztlich den Induktionsschluss für einen Induktionsbeweis, den wir im folgenden führen werden.

Theorem 1: Sei $n\in {\mathbb{N}}_0$ beliebig. Dann gehört $n$ zum Eigenwertspektrum von $N$.

Beweis: Wie angekündigt, führen wir den Beweis mittels der vollständigen Induktion. Induktionsanfang: $n=0$. Wir müssen beweisen, dass 0 im Eigenwertspektrum von $N$ enthalten ist. Dazu betrachten wir das gesamte Problem im Ortsraum. Mit Hilfe von Satz 5 finden wir leicht eine notwendige Bedingung an einen möglichen Eigenzustand $|{\psi }_0⟩$ zum Eigenwert 0: Er muss der Bedingung

$a|{\psi }_0⟩=0$

genügen. Der Operator $a$ ist aber in der Ortsdarstellung durch

$a=\sqrt{\frac{m\omega }{2\hslash }}X+i\frac{1}{\sqrt{2m\hslash \omega }}P$

gegeben, wobei $X=x$ der Ortsoperator und $P=\frac{\hslash }{i}\frac{d}{dx}$ der Impulsoperator ist. Die Bedingung $a|\psi ⟩=0$wird dann nach wenigen leichten Umformungen zu der linearen Differenzialgleichung

$m\omega x\psi +\hslash \frac{d\psi }{dx}=0$

oder

$\frac{d\psi }{dx}=-\frac{m\omega x\psi }{\hslash }.$

Wie aus der Mathematik bekannt (durch "Trennung der Variablen" ), besitzt diese einzig die Lösungsschar

${\psi }_0(x)=A{e}^{-\frac{1}{2}\frac{m\omega }{\hslash }{x}^2}$

Der Einfachheit halber können wir sogar $A=1$ wählen (auch wenn $\psi$ dann i.A. nicht normiert ist), dann ist ${\psi }_0$ mit Sicherheit nicht die Nullfunktion, und wir müssen nun nur noch zeigen, dass diese Funktion tatsächlich Eigenfunktion zum Eigenwert 0 des Operators $N$ ist. Wegen

$N=\frac{H}{\hslash \omega }-\frac{1}{2}$

und

$H=-\frac{{\hslash }^2}{2m}\frac{d^2}{d{x}^2}+\frac{1}{2}m{\omega }^2{x}^2,$

in der Ortsdarstellung lässt sich dies sofort durch Einsetzen von ${\psi }_0$ in die Gleichung

$N{\psi }_0=0{\psi }_0=0,$

verifizieren. Damit ist der Induktionsanfang getan. Induktionsschluss:$n\to n+1$. Sei nun $n\in {\mathbb{N}}_0$ ein Eigenwert des Operators $N$, $|{\psi }_n^j⟩$ ein zugehöriger Eigenzustand. Nach Satz 6 ist $a^{+}|{\psi }_n^j⟩$ Eigenzustand zu $n+1$, also ist mit $n$ auch $n+1$ ein Eigenwert von $N$ und es ist alles gezeigt.

Schritt 3: Ausschluss weiterer Eigenwerte

Wir wissen nun, dass alle natürlichen Zahlen zum Eigenwertspektrum von $N$ gehören. Nun müssen wir uns die Frage stellen, ob es noch weitere Eigenwerte geben kann. Wegen Satz 4 kommen nur noch positive reelle Zahlen in Frage, zu untersuchen wären als mögliche weitere Eigenwerte also noch positive, nicht-natürliche Zahlen. Der nächste Satz zeigt, dass keiner dieser Zahlen zum Eigenwertspektrum von $N$ gehört.

Theorem 2: $N$ besitzt keine nicht-natürlichen Eigenwerte.

Beweis:: Sei $\nu >0$ irgendeine reelle Zahl. Ist $\nu$ nicht-natürlich, so gibt es eine natürliche Zahl $n$ mit

$n<\nu <n+1,$

sowohl $n$ also auch $n+1$ gehören aber nach Theorem 1 zum Eigenwertspektrum von $N$. Wir zeigen nun induktiv, dass zwischen einer natürlichen Zahl $n$ und ihrem Nachfolger $n+1$ kein Eigenwert liegen kann. Induktionsanfang: $n=0$. Angenommen, es gäbe einen Eigenwert $\nu$ von $N$ mit $0<\nu <1$. Dann gäbe es zu $\nu$ einen Eigenzustand $|{\psi }_{\nu }^j⟩$ und $a|{\psi }_{\nu }^j⟩$ wäre wegen $\nu \ne 0$ nach Satz 6 Eigenzustand zum Eigenwert $\nu -1$. Wegen $\nu <1$ ist aber $\nu -1<0$, Widerspruch zu Satz 4. Induktionsschritt: $n\to n+1$. Sei der Satz für ein $n\in \mathbb{N}$ bereits bewiesen. Nach Induktionsvoraussetzung gibt es dann keinen Eigenwert $\nu$ von $N$ mit $\nu \in ]n,n+1[$. Angenommen, es gäbe einen Eigenwert $\nu$ von $N$ mit $n+1<\nu <n+2$. Dann gäbe es zu $\nu$ einen Eigenzustand $|{\psi }_{\nu }^j⟩$, wegen $\nu \ne 0$ wäre nach Satz 5 $a|{\psi }_{\nu }^j⟩$ aber ein Eigenzustand zum Eigenwert $\nu -1\in ]n,n+1[$, was aber einen Widerspruch zur Induktionsvoraussetzung darstellt.

Ergebnis

Zusammen mit Satz 2 folgt aus den letzten beiden Theoremen abschließend:

Eigenwertspektrum des harmonischen Oszillators. Der Hamilton-Operator $H$ des eindimensionalen harmonischen Oszillators, $H=\frac{1}{2m}{P}^2+\frac{1}{2}m{\omega }^2{X}^2$ besitzt das Eigenwertspektrum ${ϵ}_n=(n+\frac{1}{2})\hslash \omega ,n\in {\mathbb{N}}_0=\{0,1,2,...\}.$

Animierte Schrödingerwellen des Oszillators

Die stationären Wellenfunktionen des Oszillators werden andeutungsweise in Quantenmechanik/ Darstellungen gezeigt. Nicht-stationäre kohärente Gauss-Wellenpakete kommen vor in Quantenmechanik/ Unschärferelation.

In den Animationen auf der Unterseite Quantenmechanik/ Oszillatorbilder wird die zeitveränderliche Amplitude von Oszillator-Wellenfunktion gezeigt und in Farbe ist die Phase hinzugefügt. Vorsichtig anzuklicken wegen Netzwerk-, Speicher- und Prozessorlast?

Der Grundzustand ist eine stehende Gausskurve, deren Phase gleichmäßig um den Einheitskreis läuft.

Ein kohärenter Zustand schwingt bei gleichbleibender Gauss-Form wie ein klassischer Oszillator. An der Phase, die mehr oder weniger kurze Wellen vorweist, liest man die Geschwindigkeit ab.

Ein Katzen-Zustand ist die Superposition von zwei kohärenten Zuständen. Ein zweigeteiltes Pendel, so unstabil wie Schrödingers gleichzeitig tote und lebendige Katze.

Ein gequetschter Zustand ist eine Verformung des kohärenten Zustands. Die Amplitude oder die Phase einer solchen 'atmenden' Welle streut wesentlich weniger als bei den formstabilen kohärenten Zuständen.

Gequetschte elektromagnetische Zustände werden mit Erfolg benutzt, um das Quantenrauschen entweder der Amplitude oder der Phase von Lichtwellen zu vermindern. Beides zugleich verbietet die Unschärferelation. Bei hochsensiblen Experimenten wie den LIGO-Detektoren für Gravitationswellen werden Interferometer mit phasengequetschtem Licht eingesetzt.

Literatur

• w:Claude Cohen-Tannoudji u.a.: Quantenmechanik. Band 1.