Zeitentwicklungsoperator

Der Zeitentwicklungsoperator ist ein quantenmechanischer Operator, mit dem sich die zeitliche Entwicklung eines physikalischen Systems berechnen lässt. Der quantenmechanische Operator ist eng verwandt mit dem Propagator in der Quantenfeld- oder Vielteilchentheorie. Üblicherweise wird er als $U(t,t_{0})$ geschrieben und bezeichnet die Entwicklung des Systems vom Zeitpunkt $t_{0}$ zum Zeitpunkt $$ t $$ .

Konstruktion

Die Konstruktion des Zeitentwicklungsoperators kann über die formale Exponentiation der Schrödingergleichung oder direkt über seine Eigenschaften erfolgen. Die nachstehende Erläuterung folgt zunächst der ersten Route. Der Zeitentwicklungsoperator $U(t,t_{0})$ wird über die Zeitentwicklung eines beliebigen Zustandes $|\psi \rangle$ zu einem Zeitpunkt $t_{0}$ bis zum Zeitpunkt $$ t $$ definiert:

|\psi (t)\rangle =U(t,t_{0})|\psi (t_{0})\rangle \quad \forall |\psi \rangle

Einsetzen in die Schrödingergleichung liefert einen Satz gewöhnlicher Differentialgleichungen 1. Ordnung:

\mathrm {i} \hbar {\tfrac {\partial }{\partial t}}U(t,t_{0})=H(t)U(t,t_{0})

Diese Gleichungen sind zur Schrödingergleichung insofern äquivalent, als sie die Erweiterung des Zeitentwicklungsoperators um einen infinitesimalen Zeitschritt $\delta t$ beschreiben:

U(t+\delta t,t_{0})=(1-{\tfrac {\mathrm {i} }{\hbar }}H(t)\delta t)U(t,t_{0})+O(\delta t^{2})

Obige Gleichung macht die Rolle des Hamiltonoperators als Erzeuger der Zeitentwicklungen deutlich.

Aus diesen Gleichungen können einige Eigenschaften von $U(t,t_{0})$ abgelesen werden:

Kontinuität: $U(t_{0},t_{0})=1$
Unitarität: $U^{\dagger }(t,t_{0})U(t,t_{0})=1$
Propagatoreigenschaft: $U(t,t_{0})=U(t,t^{\prime })U(t^{\prime },t_{0})\quad \forall t^{\prime }$

Eigenschaften 1 und 3 folgen direkt aus der Definition von $$ U $$ . Die Unitarität folgt aus der Selbstadjungiertheit von $$ H $$ und sichert die Erhaltung der Norm und damit der Gesamtwahrscheinlichkeit:

\langle \psi (t)|\psi (t)\rangle =\langle \psi (t_{0})|U^{\dagger }(t,t_{0})U(t,t_{0})|\psi (t_{0})\rangle =\langle \psi (t_{0})|\psi (t_{0})\rangle

Es kann auch der „umgekehrte“ Weg gegangen und die Zeitentwicklung als Ausgangspunkt genommen werden: Das System wird durch einen Zeitentwicklungsoperator $U(t,t_{0})$ definiert, der den Kriterien 1 bis 3 genügen muss. Es kann dann gezeigt werden, dass dieser durch einen selbstadjungierten Operator $$ H(t) $$ erzeugt wird, womit die Brücke zum Hamiltonoperator und zur Schrödingergleichung geschlagen ist.

Explizite Form

Aus der Schrödingergleichung erhalten wir zunächst nur die infinitesimale Form des Zeitentwicklungsoperators (siehe oben):

U(t+\delta t,t)=1-{\tfrac {\mathrm {i} }{\hbar }}H(t)\delta t+O(\delta t)^{2}

Endliche Zeitentwicklungen erhält man grob gesprochen entweder durch Verknüpfung von unendlich vielen infinitesimalen Zeitschritten oder durch eine Reihenentwicklung.

Zeitunabhängige Systeme

Falls der Hamiltonoperator $$ H $$ nicht explizit zeitabhängig ist, so lässt sich eine analytische Lösung für den Zeitentwicklungsoperator finden. Für die Aneinanderreihung von kleinen Zeitentwicklungen $\delta t$ erhält man im Grenzwert $\delta t$ gegen 0 ein Matrixexponential:

U(t,t_{0})=\lim _{\delta t\to 0}\left(1-{\frac {\mathrm {i} }{\hbar }}H\delta t\right)^{\frac {t-t_{0}}{\delta t}}=\exp \left(-{\frac {\mathrm {i} }{\hbar }}H\cdot (t-t_{0})\right)

Dasselbe Ergebnis folgt auch direkt aus der Schrödingergleichung für $$ U $$ , wenn man den Operator in eine Potenzreihe in $t-t_{0}$ entwickelt und einen Koeffizientenvergleich durchführt. Man erhält dann die Potenzreihenentwicklung der Exponentialfunktion^{[NB 1]}:

U(t,t_{0})=\sum _{n=0}^{\infty }{\frac {1}{n!}}\left(-{\frac {\mathrm {i} }{\hbar }}H\cdot (t-t_{0})\right)^{n}=\exp \left(-{\frac {\mathrm {i} }{\hbar }}H\cdot (t-t_{0})\right)

$U(t,t_{0})$ ist für zeitunabhängige Systeme nur von der Zeitdifferenz $\Delta t=t-t_{0}$ abhängig, was die Unabhängigkeit des Systems gegenüber der Wahl des Zeitursprungs ausdrückt (zeitliche Translationsinvarianz).

Für praktische Rechnungen verwendet man meist die Spektraldarstellung des Zeitentwicklungsoperators, bei der der „unpraktische“ Operator im Exponenten zu einem Phasenfaktor wird:

U(t,t_{0})=\sum _{n}\exp \left(-{\frac {\mathrm {i} }{\hbar }}E_{n}(t-t_{0})\right)|n\rangle \langle n|

,

wobei ${|n\rangle }$ die Energieeigenbasis ( $H|n\rangle =E_{n}|n\rangle$ ) darstellt ^{[NB 2]}. Diese Form passt mit der Lösung der zeitseparierten Schrödingergleichung zusammen.

Explizit zeitabhängige Systeme

Ist $$ H=H(t) $$ zeitabhängig, so sind im Allgemeinen nur mehr numerische Lösungen für $$ U $$ möglich, es sei denn, der Hamiltonoperator kommutiert, zu verschiedenen Zeitpunkten ausgewertet, mit sich selbst. Dann kann wie im eindimensionalen Fall die Differenzialgleichung durch Exponentiation gelöst werden und man erhält:

U(t,t_{0})=\exp \left(-{\frac {\mathrm {i} }{\hbar }}\int _{t_{0}}^{t}H(\tau )d\tau \right)\quad \Leftrightarrow \quad [H(t),H(t^{\prime })]=0

In der Regel ist dies nicht der Fall und man muss eine der beiden obigen Techniken für zeitunabhängige Systeme verallgemeinern, um zu einer Lösung zu kommen.

Betrachtet man die Zeitentwicklung wieder als Aneinanderreihung kleiner (in diesem Fall nicht äquivalenter) Zeitschritte $\delta t$ , so definiert der Grenzwert $\delta t$ gegen 0 formal ein Produktintegral nach Vito Volterra:

U(t,t_{0})=\lim _{\delta t\to 0}\prod _{t_{i}}\left(1-{\frac {\mathrm {i} }{\hbar }}H(t_{i})\delta t\right)=:\prod _{\tau =t_{0}}^{t}\left(1-{\frac {\mathrm {i} }{\hbar }}H(\tau )d\tau \right)

Dieses Bild der aneinandergereihten kleinen Zeitschritte („time slicing“) ist eine wesentliche Zutat für die Definition des Pfadintegrals. Für praktische Rechnungen betreibt man in der Regel Zeitabhängige Störungstheorie, bei der man eine Reihenentwicklung für $$ U $$ in die Schrödingergleichung in Integralform einsetzt. Die Störungsreihe ergibt sich dann zur sogenannten Dyson-Reihe:

U(t,t_{0})=\sum _{n=0}^{\infty }\left(-{\frac {\mathrm {i} }{\hbar }}\right)^{n}\int _{t_{0}}^{t}dt_{1}H(t_{1})\int _{t_{0}}^{t_{1}}dt_{2}H(t_{2})\cdots \int _{t_{0}}^{t_{n-1}}dt_{n}\;H(t_{n})

=\sum _{n=0}^{\infty }{\frac {1}{n!}}\left(-{\frac {\mathrm {i} }{\hbar }}\right)^{n}\int _{t_{0}}^{t}dt_{1}\cdots \int _{t_{0}}^{t}dt_{n}\;T\left[H(t_{1})\cdots H(t_{n})\right]

=:T\left[\exp \left(-{\frac {\mathrm {i} }{\hbar }}\int _{t_{0}}^{t}H(\tau )d\tau \right)\right].

mit dem Zeitordnungs-Operator $$ T $$ . Im eindimensionalen Fall ist dieser Operator trivial und man erhält erneut obige Gleichung. Während man sich für viele Systeme auf wenige Terme der Störungsreihe beschränken kann, gibt es manche Systeme wie etwa Nicht-Fermi-Flüssigkeiten, bei denen die Reihe nicht konvergiert.

Literatur

J. J. Sakurai: Modern Quantum Mechanics. Prentice Hall, 1993. ISBN 978-0201539295

Anmerkungen

↑ Oft wird die Exponentialfunktion für Operatoren auch über ihre Potenzreihe definiert
↑ Die Formel gilt in dieser Form nur für abzählbare Basen, kann aber leicht auch kontinuierliche Spektren verallgemeinert werden (siehe auch Spektralsatz).

[1] Oft wird die Exponentialfunktion für Operatoren auch über ihre Potenzreihe definiert

[2] Die Formel gilt in dieser Form nur für abzählbare Basen, kann aber leicht auch kontinuierliche Spektren verallgemeinert werden (siehe auch Spektralsatz).

[NB 1]

[NB 2]