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.

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.[1] 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.[2] Infolgedessen kann der Mori-Zwanzig-Formalismus auch zur Behandlung von Korrelationsfunktionen verwendet werden.[3]

Herleitung

Bearbeiten

Eine nicht explizit zeitabhängige Observable[Anm. 1]  gehorcht im Heisenbergbild der hamiltonischen Bewegungsgleichung

 

wobei der Liouville-Operator   im quantenmechanischen Fall durch den Kommutator   und im klassischen Fall durch die Poissonklammer   definiert ist (es wird hier angenommen, dass der Hamiltonoperator   nicht von der Zeit abhängt, Verallgemeinerungen auf den zeitabhängigen Fall existieren auch[4]). Diese Gleichung wird formal gelöst durch

 

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

 

wobei   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   definiert durch[5]

 

Hierbei ist   die inverse Temperatur,   ist die Spur (im klassischen Fall das Phasenraumintegral) und   der Hamiltonoperator bzw. die Hamiltonfunktion.   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.[6]

Nun wendet man die Operatoridentität

 

auf die Größe

 

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

 

(Frequenzmatrix),

 

(stochastische Kraft) und

 

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

 

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

Markov-Näherung

Bearbeiten

Aufgrund des Faltungsterms ist die obigen Gleichung meist nur schwer zu lösen. 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  , so erhält man[7]

 ,

wobei

 

ist.

Verallgemeinerungen

Bearbeiten

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.[8] In diesem Fall hängt der Projektionsoperator explizit von der Zeit ab.[Anm. 2] In diesem Fall kann die Gleichung für eine relevante Observable

 

(wobei   der Mittelwert und   die Fluktuation ist) geschrieben werden als (verwende Indexnotation mit Summenkonvention)[9]

 ,

wobei

 ,
 ,
 

und

 .

Hierbei wurden das zeitgeordnete Exponential

 

sowie der zeitabhängige Projektionsoperator

 

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

Siehe auch

Bearbeiten

Anmerkungen

Bearbeiten
  1. Eine analoge Herleitung findet sich, z. B., in Robert Zwanzig Nonequilibrium Statistical Mechanics 3rd ed., Oxford University Press, New York, 2001, S. 149 ff.
  2. Für eine detailliertere Beschreibung der Herleitung der verallgemeinerten Gleichungen siehe Hermann Grabert Nonlinear Transport and Dynamics of Fluctuations Journal of Statistical Physics, Vol. 19, No. 5, 1978 bzw. Hermann Grabert Projection operator techniques in nonequilibrium statistical mechanics, Springer Tracts in Modern Physics, Band 95, 1982

Literatur

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

Einzelnachweise

Bearbeiten
  1. Robert Zwanzig Nonequilibrium Statistical Mechanics 3rd ed., Oxford University Press, New York, 2001, S. 144 ff.
  2. a b Hermann Grabert Nonlinear Transport and Dynamics of Fluctuations Journal of Statistical Physics, Vol. 19, No. 5, 1978
  3. Jean-Pierre Hansen und Ian R. McDonald, Theory of Simple Liquids: with Applications to Soft Matter 4th ed. (Elsevier Academic Press, Oxford, 2009), S. 363 ff.
  4. a b M. te Vrugt and R. Wittkowski Mori-Zwanzig projection operator formalism for far-from-equilibrium systems with time-dependent Hamiltonians Physical Review E 99, 062118 (2019)
  5. Hermann Grabert Projection operator techniques in nonequilibrium statistical mechanics, Springer Tracts in Modern Physics, Band 95, 1982, S. 37
  6. Hermann Grabert Projection operator techniques in nonequilibrium statistical mechanics, Springer Tracts in Modern Physics, Band 95, 1982, S. 13
  7. Robert Zwanzig Nonequilibrium Statistical Mechanics 3rd ed., Oxford University Press, New York, 2001, S. 165 ff.
  8. Hermann Grabert Projection operator techniques in nonequilibrium statistical mechanics, Springer Tracts in Modern Physics, Band 95, 1982, S. 36
  9. Hermann Grabert Projection operator techniques in nonequilibrium statistical mechanics, Springer Tracts in Modern Physics, Band 95, 1982, S. 18
  10. Hugues Meyer, Thomas Voigtmann und Tanja Schilling On the dynamics of reaction coordinates in classical, time-dependent, many-body processes J. Chem. Phys. 150, 174118 (2019)
  11. A. J. Chorin, O. H. Hald und R. Kupferman Optimal prediction with memory Physica D: Nonlinear Phenomena 166, 239{257 (2002)