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.

The equation belongs to the Mori-Zwanzig formalism within the statistical mechanics of irreversible processes (named after Hazime Mori).

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).

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

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

In the Nakajima-Zwanzig formulation, a key step involves defining a projection operator

that projects the total density 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

), 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

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.

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.

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.

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.

presents difficulties due to its time-dependent nature.

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