Mori-Zwanzig-Formalismus

Der Mori-Zwanzig-Formalismus, benannt nach den Physikern Hajime Mori und Robert Zwanzig, ist eine Methode der statistischen Physik. Er ermöglicht, die Dynamik eines Systems mithilfe von Projektionsoperatoren in einen relevanten und einen irrelevanten Teil zu zerlegen, wodurch sich geschlossene Bewegungsgleichungen für den relevanten Teil finden lassen.

Idee

Makroskopische Systeme mit einer großen Zahl mikroskopischer Freiheitsgrade werden typischerweise sehr gut durch eine kleine Zahl relevanter Variablen beschrieben, etwa der Gesamtmagnetisierung bei Spinsystemen. Der Mori-Zwanzig-Formalismus ermöglicht, basierend auf bekannten mikroskopischen Bewegungsgleichungen eines Systems – typischerweise basierend auf der Hamilton-Funktion der klassischen Mechanik bzw. dem quantenmechanischen Hamiltonoperator – makroskopische Gleichungen zu erhalten, welche nur von den relevanten Variablen abhängen. Der irrelevante Teil tritt in diesen Gleichungen als Rauschen auf. Der Formalismus ermöglicht keine Aussagen darüber, was die relevanten Variablen sind, diese ergeben sich typischerweise aus den Eigenschaften des untersuchten Systems.

Die Observablen, welche das System beschreiben, bilden einen Hilbertraum. Der Projektionsoperator entspricht dann einer Abbildung auf den Unterraum, welcher nur von den relevanten Variablen aufgespannt wird. Der irrelevante Teil der Dynamik hängt von zu den relevanten Variablen orthogonalen Observablen ab. Als Skalarprodukt im Raum der dynamischen Variablen dient dabei eine Korrelationsfunktion. Infolgedessen kann der Mori-Zwanzig-Formalismus auch zur Behandlung von Korrelationsfunktionen verwendet werden.

Herleitung

Eine nicht explizit zeitabhängige Observable $A$ gehorcht im Heisenbergbild der hamiltonischen Bewegungsgleichung

${\frac {d}{dt}}A=\mathrm {i} LA,$

wobei der Liouville-Operator $L$ im quantenmechanischen Fall durch den Kommutator $L(\cdot )={\frac {1}{\hbar }}[H,\cdot ]$ und im klassischen Fall durch die Poissonklammer $L(\cdot )=-\mathrm {i} \{H,\cdot \}$ definiert ist (es wird hier angenommen, dass der Hamiltonoperator $H$ nicht von der Zeit abhängt, Verallgemeinerungen auf den zeitabhängigen Fall existieren auch). Diese Gleichung wird formal gelöst durch

$A(t)=e^{\mathrm {i} Lt}A.$

Der Projektionsoperator, der auf eine Observable $X$ wirkt, ist definiert durch

$PX=(A,A)^{-1}(X,A)A,$

wobei $A$ die relevante Observable ist (dies kann auch ein Vektor aus mehreren relevanten Größen sein). Als Skalarprodukt – hier mit runden Klammern notiert – wird typischerweise das Mori-Produkt, eine Verallgemeinerung der üblichen Korrelationsfunktion, verwendet. Dieses ist für Observablen $X,Y$ definiert durch

$(X,Y)={\frac {1}{\beta }}\int _{0}^{\beta }\mathrm {d} \alpha \,\operatorname {Tr} ({\bar {\rho }}Xe^{-\alpha H}Ye^{\alpha H}).$

Hierbei ist $\beta =(k_{B}T)^{-1}$ die inverse Temperatur, $\operatorname {Tr}$ ist die Spur (im klassischen Fall das Phasenraumintegral) und $H$ der Hamiltonoperator bzw. die Hamiltonfunktion. ${\bar {\rho }}$ ist die relevante Wahrscheinlichkeitsdichte (bzw. der Dichteoperator in der Quantenmechanik). Diese wird so gewählt, dass sie lediglich als Funktion der relevanten Observablen geschrieben werden kann, gleichzeitig aber die (i. d. R. unbekannte) tatsächliche Wahrscheinlichkeitsdichte möglichst gut approximiert, d. h. insbesondere die gleichen Erwartungswerte für die relevanten Variablen liefert wie diese.

Nun wendet man die Operatoridentität

$e^{\mathrm {i} Lt}=e^{\mathrm {i} (1-P)Lt}+\int _{0}^{t}\mathrm {d} s\,e^{\mathrm {i} L(t-s)}PLe^{\mathrm {i} (1-P)Ls}$

auf die Größe

$(1-P)\mathrm {i} LA$

an. Unter Ausnutzung der obigen Definition des Projektionsoperators und der Definitionen

$\Omega =(\mathrm {i} LA,A)(A,A)^{-1}$

(Frequenzmatrix),

$F(t)=e^{t(1-P)L}(1-P)\mathrm {i} LA$

(stochastische Kraft) und

$K(t)=(\mathrm {i} LF(t),A)(A,A)^{-1}$

(Gedächtnisfunktion) lässt sich dies schreiben als

${\dot {A}}(t)=\Omega A(t)+\int _{0}^{t}\mathrm {d} s\,K(s)A(t-s)+F(t).$

Dies ist eine Bewegungsgleichung für die relevante Observable $A(t)$, welche von dem Wert der Observable zum Zeitpunkt $t$, dem Wert zu früheren Zeitpunkten (Gedächtnisterm) und der stochastischen Kraft (Rauschen, entspricht dem Einfluss der zu $A(t)$ orthogonalen Dynamik) abhängt.

Markov-Näherung

Aufgrund des Faltungsterms ist die obigen Gleichung meist nur schwer zu lösbar. Typischerweise ist man an langsamen (makroskopischen) Variablen interessiert, deren Änderung auf größeren Zeitskalen stattfindet als das (mikroskopische) Rauschen. Entwickelt man die Gleichung bis zur zweiten Ordnung in $\mathrm {i} LA(t)$, so erhält man

${\dot {A}}(t)\approx \Omega A(t)+\int _{0}^{\infty }\mathrm {d} s\,K(s)A(t)+F(t)$,

wobei

$K(t)=-(e^{\mathrm {i} Lt}(1-P)\mathrm {i} LA,(1-P)\mathrm {i} LA)(A,A)^{-1}$

ist.

Verallgemeinerungen

Für den Fall beliebiger Abweichungen vom thermodynamischen Gleichgewicht wird die allgemeinere Form des Mori-Zwanzig-Formalismus verwendet, aus welcher sich obige Gleichungen als Linearisierung ergeben. In diesem Fall hängt der Projektionsoperator explizit von der Zeit ab. In diesem Fall kann die Gleichung für eine relevante Observable

$A(t)=a(t)-\delta A(t)$

(wobei $a(t)$ der Mittelwert und $\delta A(t)$ die Fluktuation ist) geschrieben werden als (verwende Indexnotation mit Summenkonvention)

${\dot {A}}_{i}(t)=v_{i}(t)+\Omega _{ij}(t)\delta A_{j}(t)+\int _{0}^{t}\mathrm {d} s\,K_{i}(t,s)+\phi _{ij}(t,s)\delta A_{j}(t)+F_{i}(t,0)$,

wobei

$v_{i}(t)=\operatorname {Tr} \left({\bar {\rho }}(t)A_{i}\right)$,

$\Omega _{ij}(t)=\operatorname {Tr} \left({\frac {\partial {\bar {\rho }}(t)}{\partial a_{j}(t)}}{\dot {A}}_{i}\right)$,

$K_{i}(t,s)=\operatorname {Tr} \left({\bar {\rho }}(s)\mathrm {i} L(1-P(s))G(s,t){\dot {A}}_{i}\right),$

und

$\phi _{ij}(t,s)=\operatorname {Tr} \left({\frac {\partial {\bar {\rho }}(t)}{\partial a_{j}(t)}}\mathrm {i} L(1-P(s))G(s,t){\dot {A}}_{i}\right)-{\dot {a}}_{k}(t)\operatorname {Tr} \left({\frac {\partial ^{2}{\bar {\rho }}(t)}{\partial a_{k}(t)\partial a_{j}(t)}}G(s,t){\dot {A}}_{i}\right)$.

Hierbei wurden das zeitgeordnete Exponential

$G(s,t)=T_{-}\exp \left(\int _{s}^{t}\mathrm {d} u\,\mathrm {i} L(1-P(u))\right)$

sowie der zeitabhängige Projektionsoperator

$P(t)X=\operatorname {Tr} \left({\bar {\rho }}(t)X)+(A_{i}-a_{i}(t)\right)\operatorname {Tr} \left({\frac {\partial {\bar {\rho }}(t)}{\partial a_{i}(t)}}X\right)$

verwendet. Auch diese Gleichungen können mittels einer Verallgemeinerung des Mori-Produkts durch Korrelationsfunktionen dargestellt werden. Weitere Varianten des Mori-Zwanzig-Formalismus dienen der Beschreibung von Systemen mit zeitabhängigen Hamiltonoperatoren bzw. von allgemeinen dynamischen Systemen.

Siehe auch

• Nakajima-Zwanzig-Gleichung

Literatur

• Hermann Grabert Projection operator techniques in nonequilibrium statistical mechanics, Springer Tracts in Modern Physics, Band 95, 1982
• Robert Zwanzig Nonequilibrium Statistical Mechanics 3rd ed., Oxford University Press, New York, 2001