Nakajima–Zwanzig equation

The Nakajima–Zwanzig equation (named after the physicists who developed it, Sadao Nakajima[1] and Robert Zwanzig[2]) is an integral equation describing the time evolution of the "relevant" part of a quantum-mechanical system. It is formulated in the density matrix formalism and can be regarded as a generalization of the master equation.

Four approaches to obtaining the memory kernel of the GQME considering the system-bath form of the Hamiltonian.

The equation belongs to the Mori-Zwanzig formalism within the statistical mechanics of irreversible processes (named after Hazime Mori). By means of a projection operator, the dynamics is split into a slow, collective part (relevant part) and a rapidly fluctuating irrelevant part. The goal is to develop dynamical equations for the collective part.

The Nakajima-Zwanzig (NZ) generalized master equation is a formally exact approach for simulating quantum dynamics in condensed phases. This framework is particularly designed to address the dynamics of a reduced system interact with a larger environment, often represented as a system coupled to a bath.  Within the NZ framework, one can choose between time convolution (TC) and time convolution less (TCL) forms of the quantum master equations.

The TC approach involves memory effects, where the future state of the system depends on its entire history (Non-Markovian dynamics). The TCL approach formulates the dynamics where the system's rate of change at any moment depends only on its current state, simplifying calculations by neglecting memory effects (Markovian dynamics).

Derivation

edit

The total Hamiltonian of a system interacting with its environment (or bath) is typically expressed in system-bath form,

 

where   is the system Hamiltonian,   is the bath Hamiltonian, and   describes the coupling between them.

The starting point[note 1] is the quantum mechanical version of the von Neumann equation, also known as the Liouville equation:

 

where the Liouville operator   is defined as  .

In the Nakajima-Zwanzig formulation, a key step involves defining a projection operator   that projects the total density operator   onto the subspace of the system of interest. The complementary operator   projects onto the orthogonal subspace, effectively separating the system from the bath.

The Liouville – von Neumann equation can thus be represented as

 

The dynamics of the projected state  , under any idempotent projection operator (where ), is described by the NZ generalized master equation (GQME). This equation can be used to obtain a closed equation of motion for the reduced system dynamics, focusing solely on the dynamics within the subsystem of interest.  


In practice, the specific form of the projection operator   can be chosen based on the problem at hand. One common choice involves defining   using a reference nuclear density operator   such that  .

 

This ensures that   remains idempotent. Using this projection, tracing over the nuclear Hilbert space leads to a generalized quantum master equation that describes the reduced electronic density operator which accounts for both Markovian dynamics generated by the Hamiltonian and non-Markovian dynamics due to coupling between electronic and nuclear degrees of freedom.

 

This  

describes the dynamics driven by the Hamiltonian,  which are Hamiltonian and Markovian in nature, while the other two terms on the right-hand side represent the non-Hamiltonian and non-Markovian dynamics that arise from the interactions between the electronic and nuclear degrees of freedom.

The memory kernel   captures the effects of the bath on the system over the time interval from (0, t), reflecting non-Markovian dynamics where the system's history influences its future evolution.

 

The inhomogeneous term   represents the influence of the initial state of the bath on the system at time t, which is crucial for accurately describing the system dynamics from an initial condition.

 

The memory kernel is crucial for simulating the dynamics of the electronic degrees of freedom. However, calculating   presents difficulties due to its time-dependent nature. Additionally, the time dependency of   is complex because it is governed by the projection-dependent propagator,  . Therefore, the exact memory kernel is difficult to calculate except for several analytically solvable models proposed by Shi-Geva to remove the projection operator  .


See also

edit

Notes

edit
  1. ^ A derivation analogous to that presented here is found, for instance, in Breuer, Petruccione The theory of open quantum systems, Oxford University Press 2002, S.443ff

References

edit
  1. ^ Nakajima, Sadao (1 December 1958). "On Quantum Theory of Transport Phenomena: Steady Diffusion". Progress of Theoretical Physics. 20 (6): 948–959. Bibcode:1958PThPh..20..948N. doi:10.1143/PTP.20.948. ISSN 0033-068X.
  2. ^ Zwanzig, Robert (1960). "Ensemble Method in the Theory of Irreversibility". The Journal of Chemical Physics. 33 (5): 1338–1341. Bibcode:1960JChPh..33.1338Z. doi:10.1063/1.1731409.
  • E. Fick, G. Sauermann: The Quantum Statistics of Dynamic Processes Springer-Verlag, 1983, ISBN 3-540-50824-4.
  • Heinz-Peter Breuer, Francesco Petruccione: Theory of Open Quantum Systems. Oxford, 2002 ISBN 9780198520634
  • Hermann Grabert Projection operator techniques in nonequilibrium statistical mechanics, Springer Tracts in Modern Physics, Band 95, 1982
  • R. Kühne, P. Reineker: Nakajima-Zwanzig's generalized master equation: Evaluation of the kernel of the integro-differential equation, Zeitschrift für Physik B (Condensed Matter), Band 31, 1978, S. 105–110, doi:10.1007/BF01320131
  • Xu, M.; Yan, Y.; Liu, Y.; Shi, Q. Convergence of High Order Memory Kernels in the Nakajima-Zwanzig Generalized Master Equation and Rate Constants: Case Study of the Spin-Boson Model. Journal of Chemical Physics 2018, 148 (16). https://doi.org/10.1063/1.5022761.
  • Mulvihill, E.; Geva, E. A Road Map to Various Pathways for Calculating the Memory Kernel of the Generalized Quantum Master Equation. Journal of Physical Chemistry B 2021, 125 (34), 9834–9852. https://doi.org/10.1021/acs.jpcb.1c05719.
  • Mulvihill, E.; Schubert, A.; Sun, X.; Dunietz, B. D.; Geva, E. A Modified Approach for Simulating Electronically Nonadiabatic Dynamics via the Generalized Quantum Master Equation. Journal of Chemical Physics 2019, 150 (3). https://doi.org/10.1063/1.5055756.
edit