Markow-Näherung
# Markow-Näherung

Als Markow-Näherung (nach Andrei Andrejewitsch Markow) bezeichnet man eine Näherungsmethode der Quantenoptik. In der Markow-Näherung geht man davon aus, dass ein Reservoir, an das ein quantenmechanisches System gekoppelt ist, nach der Kopplung instantan in seinen Ausgangszustand zurückkehrt. Anschaulich hat das Reservoir damit kein „Gedächtnis“ und verhält sich sehr schnell wieder so, als hätte die Wechselwirkung mit dem angekoppelten System überhaupt nicht stattgefunden.
Die Markow-Näherung wird beispielsweise bei der Bestimmung der Mastergleichung von gedämpften harmonischen Oszillatoren verwendet. Außerdem ist sie eine nützliche Annahme für die numerische Berechnung von quantenphysikalischen Vorgängen in diskreten Zeitschritten (engl. coarse graining).

## Herleitung
### Grundlagen

Dieser Abschnitt behandelt ein quantenmechanisches System (Index S) und ein Reservoir (Index R) sowie die Wechselwirkung zwischen beiden (Index SR). Die Dynamik des Gesamtsystems wird mit einem Hamilton-Operator $H$ beschrieben, der sich aus den einzelnen Hamiltonians des Systems und des Reservoirs sowie einem Wechselwirkungs-Hamiltonian zusammensetzt:
    $H=H_{S}+H_{R}+H_{SR}$
Der Dichteoperator des Gesamtsystems ist $\chi$. Der Dichteoperator $\rho$ des einzelnen Systems und der des Reservoirs $R$ sind durch Bildung der Partialspur über die Freiheitsgrade des Reservoirs bzw. des angekoppelten Systems bestimmt:
    $\rho =\operatorname {Spur} _{R}(\chi )$
    $R=\operatorname {Spur} _{S}(\chi )$
Die Zeitentwicklung von $\chi$ wird durch die Von-Neumann-Gleichung beschrieben:
    ${\dot {\chi }}={\frac {\partial \chi }{\partial t}}={\frac {1}{\mathrm {i} \hbar }}[H,\chi ]$

### Übergang ins Wechselwirkungsbild

Nun transformiert man die Von-Neumann-Gleichung ins Wechselwirkungsbild (gekennzeichnet durch Tilden ${\tilde}$), indem man den unitären Zeitentwicklungsoperator $U$ anwendet.
    ${\dot {\tilde {\chi }}}=\overbrace {e^{\mathrm {i} /\hbar \left(H_{S}+H_{R}\right)t}} ^{U^{\dagger }}\;{\dot {\chi }}\;\overbrace {e^{-\mathrm {i} /\hbar \left(H_{S}+H_{R}\right)t}} ^{U}$
Die unitäre Transformation beim Übergang ins Wechselwirkungsbild absorbiert quasi die schnelle Zeitentwicklung aufgrund der ungestörten Terme $H_{S}$ und $H_{R}$, sodass lediglich die Kopplung von System und Reservoir für die zeitliche Entwicklung von ${\tilde {\chi }}$ verantwortlich ist. Weil $U$ eine Funktion von $H_{S}+H_{R}$ ist, kommutiert $U$ mit $H_{S}+H_{R}$; jedoch nicht mit $H_{SR}$. Daraus folgt:
    ${\dot {\tilde {\chi }}}={\frac {1}{\mathrm {i} \hbar }}\left[{\tilde {H}}_{SR},{\tilde {\chi }}\right]$
Diese Differentialgleichung lässt sich formal durch
    ${\tilde {\chi }}(t)=\chi (0)+{\frac {1}{\mathrm {i} \hbar }}\int _{0}^{t}\left[{\tilde {H}}_{SR}(\tau ),{\tilde {\chi }}(\tau )\right]\mathrm {d} \tau$
lösen und wird wieder in die Von-Neumann-Gleichung eingesetzt:
    ${\dot {\tilde {\chi }}}={\frac {1}{\mathrm {i} \hbar }}\left[{\tilde {H}}_{SR}(t),{\tilde {\chi }}(0)\right]-{\frac {1}{\hbar ^{2}}}\int _{0}^{t}\left[{\tilde {H}}_{SR}(t),\left[{\tilde {H}}_{SR}(\tau ),{\tilde {\chi }}(\tau )\right]\right]\mathrm {d} \tau$
Mit einem weiteren Iterationsschritt erhält man folgende immer noch exakte Lösung für ${\tilde {\chi }}(t)$.
    ${\tilde {\chi }}(t)={\tilde {\chi }}(0)+{\frac {1}{\mathrm {i} \hbar }}\int _{0}^{t}\left[{\tilde {H}}_{SR}(\tau ),{\tilde {\chi }}(0)\right]\mathrm {d} \tau -{\frac {1}{\hbar ^{2}}}\int _{0}^{t}\int _{0}^{\tau }\left[{\tilde {H}}_{SR}(\tau ),\left[{\tilde {H}}_{SR}(\tau '),{\tilde {\chi }}(\tau ')\right]\right]\mathrm {d} \tau '\mathrm {d} \tau$
Aus ${\tilde {\chi }}(t)$ lässt sich durch Ausspuren des Reservoirs auch ${\tilde {\rho }}(t)$ bestimmen:
    ${\begin{aligned}{\tilde {\rho }}(t)&=\operatorname {Spur} _{R}\left({\tilde {\chi }}(t)\right)\\&=\operatorname {Spur} _{R}\left({\tilde {\chi }}(0)\right)+{\frac {1}{\mathrm {i} \hbar }}\int _{0}^{t}\underbrace {\operatorname {Spur} _{R}\left(\left[{\tilde {H}}_{SR}(\tau ),{\tilde {\chi }}(0)\right]\right)} _{=0}\mathrm {d} \tau -{\frac {1}{\hbar ^{2}}}\int _{0}^{t}\int _{0}^{\tau }\operatorname {Spur} _{R}\left(\left[{\tilde {H}}_{SR}(\tau ),\left[{\tilde {H}}_{SR}(\tau '),{\tilde {\chi }}(\tau ')\right]\right]\right)\mathrm {d} \tau '\mathrm {d} \tau \end{aligned}}$

### Annahmen an das Reservoir / Bornsche Näherung

Zur Vereinfachung nimmt man an, dass die Auswirkungen der an das System koppelnden Operatoren des Reservoirs zum Zeitpunkt $t=0$ im Mittel keinen Einfluss auf das Reservoir haben. Der Erwartungswert von ${\tilde {H}}_{SR}$ verschwindet also.
    $\operatorname {Spur} _{R}\left({\tilde {H}}_{SR}R_{0}\right)=0$
Dies lässt sich zeigen, indem man ${\tilde {H}}_{SR}$ wie im folgenden Abschnitt in ein Produkt aus zwei Operatoren zerlegt, von denen einer aus der Spur herausgezogen werden kann. Daraus folgt direkt, dass
    $\operatorname {Spur} _{R}\left(\left[{\tilde {H}}_{SR},{\tilde {\chi }}(0)\right]\right)=0$
Des Weiteren nimmt man an, dass sich die Dichtematrix des Gesamtsystems zum Zeitpunkt $\chi (0)=\rho (0)R(0)=\rho (0)R_{0}$ faktorisieren lässt. Nun wendet man die Bornsche Näherung an, d. h., man vernachlässigt Beiträge höherer Ordnungen in ${\tilde {H}}_{SR}$:
    ${\tilde {\chi }}(t)\approx {\tilde {\rho }}(t)R_{0}+{\mathcal {O}}({\tilde {H}}_{SR})$
So ergibt sich schließlich für die Dichtematrix des Systems:
    ${\tilde {\rho }}(t)={\tilde {\rho }}(0)-{\frac {1}{\hbar ^{2}}}\int _{0}^{t}\int _{0}^{\tau }\operatorname {Spur} _{R}\left(\left[{\tilde {H}}_{SR}(\tau ),\left[{\tilde {H}}_{SR}(\tau '),{\tilde {\rho }}(\tau ')R_{0}\right]\right]\right)\mathrm {d} \tau '\mathrm {d} \tau$

### Reservoir-interne Korrelationsfunktionen

Man schreibt den Wechselwirkungs-Hamiltonian $H_{SR}$ mit Hilfe von Operatoren $s_{i}$, die nur auf das angekoppelte System wirken und Operatoren $r_{i}$, die nur auf das Reservoir wirken. Somit ist:
    ${\tilde {H}}_{SR}=\hbar \sum _{i}{\tilde {s}}_{i}{\tilde {r}}_{i}$
Diesen Ausdruck für ${\tilde {H}}_{SR}$ setzt man nun in die letzte Gleichung des vorherigen Abschnittes ein. Dabei ist zu beachten, dass Operatoren, die nur auf dem angekoppelten System wirken, aus der Spur über die Freiheitsgrade des Reservoirs herausgezogen werden können. Des Weiteren nutzt man die Invarianz der Spur unter zyklischen Permutationen aus und erhält somit
    ${\begin{aligned}{\tilde {\rho }}={\tilde {\rho }}(0)-\sum _{i,j}\int _{0}^{t}\int _{0}^{\tau }{\Big&\left({\tilde {s}}_{i}\left(\tau \right){\tilde {s}}_{j}\left(\tau '\right){\tilde {\rho }}\left(\tau '\right)-{\tilde {s}}_{j}\left(\tau '\right){\tilde {\rho }}\left(\tau '\right){\tilde {s}}_{i}\left(\tau \right)\right){\color {red}\operatorname {Spur} _{R}\left(R_{0}{\tilde {r}}_{i}\left(\tau \right){\tilde {r}}_{j}\left(\tau '\right)\right)}\\+&\left({\tilde {\rho }}\left(\tau '\right){\tilde {s}}_{j}\left(\tau '\right){\tilde {s}}_{i}\left(\tau \right)-{\tilde {s}}_{i}\left(\tau \right){\tilde {\rho }}\left(\tau '\right){\tilde {s}}_{j}\left(\tau '\right)\right){\color {red}\operatorname {Spur} _{R}\left(R_{0}{\tilde {r}}_{j}\left(\tau '\right){\tilde {r}}_{i}\left(\tau \right)\right)}{\Big )}\mathrm {d} \tau '\mathrm {d} \tau \end{aligned}}$
Die rot markierten Ausdrücke sind Korrelationsfunktionen des Reservoirs. Sie hängen nur von der Zeitdifferenz $\Delta \tau =\tau -\tau '$ ab und lassen sich in Exponentialfunktionen entwickeln, die mit verschiedenen Frequenzen $\omega _{\mu \nu }$ oszillieren.
    ${\begin{aligned}\operatorname {Spur} _{R}\left(R_{0}{\tilde {r}}_{i}\left(\tau \right){\tilde {r}}_{j}\left(\tau '\right)\right)&=\operatorname {Spur} _{R}\left(R_{0}{\tilde {r}}_{i}\left(\Delta \tau \right){\tilde {r}}_{j}\left(0\right)\right)\\&=\operatorname {Spur} _{R}\sum _{\mu }\left(p_{\mu }|\mu \rangle \langle \mu |{\tilde {r}}_{i}\left(\Delta \tau \right){\tilde {r}}_{j}\left(0\right)\right)\\&=\sum _{\mu }p_{\mu }\langle \mu |{\tilde {r}}_{i}\left(\Delta \tau \right){\tilde {r}}_{j}\left(0\right)|\mu \rangle \\&=\sum _{\mu }\sum _{\nu }p_{\mu }\left|R_{\mu \nu }\right|^{2}e^{i\omega _{\mu \nu }t}\end{aligned}}$
In der zweiten Zeile wurde hierbei angenommen, dass das Reservoir in einem Eigenzustand von $H_{R}$ ist und $R_{0}$ somit nur Diagonal-Einträge $p_{\mu }$ hat. Dabei ist $R_{\mu \nu }=\langle \mu |r|\nu \rangle$ und $\hbar \omega _{\mu \nu }=\hbar \left(\omega _{\mu }-\omega _{\nu }\right)=E_{\mu }-E_{\nu }$ die Energiedifferenz zweier Eigenzustände des Reservoirs.

### Durchführen der Markow-Näherung

Die Markow-Näherung besteht nun darin, anzunehmen, dass im Reservoir eine Vielzahl von Energieniveaus vorhanden sind, deren zugeordnete Korrelationsfunktionen destruktiv zu null interferieren, sobald $\Delta \tau$ größer ist als die internen Zeitskalen des Reservoirs. Man nimmt an, dass die Korrelationsfunktionen als Deltafunktionen genähert werden können.
    $\operatorname {Spur} _{R}\left(R_{0}{\tilde {r}}_{i}(\tau ){\tilde {r}}_{j}(\tau ')\right)\propto \delta \left(\tau -\tau '\right)$
Damit wird das Integral zu einer Faltung mit einer Deltafunktion. Dies hat zur Folge, dass in der letzten Gleichung des vorherigen Abschnittes $\rho (\tau ')$ durch $\rho (\tau )$ ersetzt werden kann. Anschaulich bedeutet die Näherung mit einer Deltafunktion, dass das Reservoir nach der Ankopplung des Systems instantan in seinen Ausgangszustand zurückkehrt. Dies geschieht auf einer Zeitskala, die wesentlich kürzer ist als die, auf der die relevanten Prozesse des angekoppelten Systems ablaufen. Im Frequenzraum ist die Korrelationsfunktion des Reservoirs somit viel breiter als die des angekoppelten Systems. Diese deutliche Separation zweier Zeitskalen ist die wesentliche Bedingung für die Gültigkeit der Markow-Näherung.
Soll also das Verhalten von $\rho$ untersucht werden, so kann man in der Born-Markow-Näherung annehmen, dass das Reservoir nach der Ankopplung des Systems in seinen Ausgangszustand zurückkehrt (Bornsche Näherung) und dies quasi sofort geschieht (Markow-Näherung).