---
title: Kadanoff–Baym Equations for Nonequilibrium Dynamics
url: https://www.emergentmind.com/topics/kadanoff-baym-equations
type: topic
---

# Kadanoff–Baym Equations for Nonequilibrium Dynamics

Kadanoff–Baym equations are the nonequilibrium Dyson equations for contour-ordered single-particle Green’s functions on the Keldysh contour. They provide a two-time, non-Markovian formulation of interacting quantum dynamics in which the one-body propagator evolves under a one-body Hamiltonian and a time-nonlocal self-energy, while initial correlations are incorporated through the imaginary branch of the contour. In this form, they are a central tool of nonequilibrium Green’s-function theory for correlated dynamics, transport, spectroscopy, and open-system evolution [2112.05670, 1006.2912].

## 1. Contour formulation and real-time structure

The basic object is the contour-ordered Green’s function
\[
G(1,2)=-i\langle T_\gamma[\psi(1)\psi^\dagger(2)]\rangle,
\]
with \(1\equiv(r_1,\sigma_1,t_1)\), \(2\equiv(r_2,\sigma_2,t_2)\), and \(T_\gamma\) denoting ordering on the Keldysh contour. In compact notation, the Kadanoff–Baym equation can be written as
\[
\left(i\partial_{t_1}-h(t_1)\right)G(t_1,t_2)=\delta(t_1,t_2)+\int_\gamma dt\,\Sigma(t_1,t)\,G(t,t_2),
\]
together with the adjoint equation differentiated with respect to \(t_2\). The same content may be expressed as the contour Dyson equation
\[
G=G_0+G_0\Sigma G,
\]
with \(\Sigma=\Sigma[G]\) the many-body self-energy functional [1009.2917, 1312.0214].

For nonequilibrium initial-value problems, the contour is typically L-shaped: forward real time, backward real time, and an imaginary branch for the initial thermal state. Projection onto contour sectors yields the standard set of components \(G^{<,>}\), \(G^{R/A}\), \(G^{\rceil,\lceil}\), and \(G^{\mathrm M}\). The mixed and Matsubara components encode initial correlations, and KMS boundary conditions enforce consistency with the initial equilibrium density matrix. This decomposition converts the contour equation into coupled integro-differential equations for real-time propagation, but it preserves the causal memory structure: the value at time \(t\) depends on the entire earlier history through self-energy convolutions [1006.2912, 2112.05670].

This two-time formulation is the defining structural feature of the Kadanoff–Baym framework. It differs from time-local kinetic equations by retaining spectral evolution, memory, and initial-correlation effects at the level of the one-particle propagator itself. That structure is also what makes the equations computationally demanding.

## 2. Self-energies, conserving structure, and excitation spectra

In practical calculations the Kadanoff–Baym equations are closed by a self-consistent approximation for \(\Sigma[G]\). The papers considered here repeatedly emphasize the role of \(\Phi\)-derivable approximations, for which
\[
\Sigma=\delta\Phi/\delta G.
\]
When solved self-consistently, Hartree–Fock, second Born, \(GW\), and \(T\)-matrix approximations are conserving in the Baym–Kadanoff sense: they enforce exact conservation of total particle number, energy, and momentum during the time evolution [1006.2912, 1103.4054].

The diagrammatic content of these approximations matters. Hartree–Fock is time-local and describes only mean-field renormalization. Second Born adds the second-order bubble and exchange diagrams and is the simplest time-nonlocal approximation. \(GW\) replaces bare interactions by a screened interaction \(W\) obtained from bubble resummation. The \(T\)-matrix approximation resums particle-particle ladder diagrams and is therefore particularly sensitive to repeated two-particle scattering processes [1009.2917].

A major consequence is visible in excitation spectra. Time-local Hartree–Fock can give incomplete neutral spectra, while second Born already reproduces many additional excitations characterized as double excitations [1108.4106]. In finite quantum wells, real-time propagation of the full two-time Kadanoff–Baym equations in second Born captures correlation-induced double-excitation features absent at the Hartree–Fock level [1203.4765]. In the language used there, propagating the Kadanoff–Baym equations with a conserving self-energy is equivalent to solving a Bethe–Salpeter equation with a correspondingly dressed kernel while fulfilling a frequency sum rule [1108.4106].

This establishes a central methodological point: the Kadanoff–Baym equations are not only transport equations for one-particle observables. Through the choice of \(\Sigma[G]\), they define a self-consistent approximation to neutral and charged excitation spectra, damping, and correlation satellites.

## 3. Higher-order observables and positivity constraints

Although the Kadanoff–Baym equations evolve a one-particle Green’s function, equations of motion permit extraction of certain two-particle observables. For a local Hubbard interaction, the local double occupancy
\[
d_R=\langle \hat n_{R\uparrow}\hat n_{R\downarrow}\rangle
\]
can be written as
\[
d_R=-\frac{i}{U_R}\left[\int_\gamma \Sigma(1,3)G(3,1^+)\,d3\right]_{RR}.
\]
This relation is exact at the formal level; the approximation enters through the chosen self-energy \(\Sigma[G]\) [1009.2917, 1103.4054].

That identity exposes an important limitation of conserving approximations. In a Hubbard dimer driven out of equilibrium, second Born and \(GW\) can yield negative double occupancy at finite times, even though both approximations are \(\Phi\)-derivable and conserve particle number and energy. Since \(d_R\) is the expectation value of a positive operator, negative values are unphysical. Among the tested schemes, only the \(T\)-matrix approximation yields double occupancies that remain non-negative in all studied cases [1009.2917].

The same issue appears in entanglement calculations. For a non-magnetic spin-\(1/2\) system, the local single-site entanglement entropy can be reconstructed from the density \(n_R\) and double occupancy \(d_R\). In a transport setup with an interacting impurity, \(d_R\) obtained from Kadanoff–Baym propagation within second Born or \(GW\) can again become negative, which makes the entropy formula ill-defined. In the \(T\)-matrix approximation, by contrast, the pair correlation function is proven to remain non-negative; this provides a controlled route from Kadanoff–Baym dynamics to entanglement entropy in the low-density regime [1103.4054].

These examples sharply delimit what conservation does and does not guarantee. Conserving approximations preserve macroscopic balance laws, but they do not automatically preserve microscopic positivity constraints on pair densities or related local observables. The Kadanoff–Baym framework therefore requires approximation diagnostics beyond conservation alone.

## 4. Finite, open, and dissipative systems

The formalism applies to both isolated and embedded systems. In quantum transport, the central interacting region is coupled to noninteracting leads that are integrated out into an embedding self-energy,
\[
\Sigma_{\rm em}(z,z')=\sum_\alpha H_{C\alpha}\,g_{\alpha\alpha}(z,z')\,H_{\alpha C}.
\]
The resulting embedded Kadanoff–Baym equation contains \(\Sigma^{\rm MB}+\Sigma_{\rm em}\) and consistently incorporates initial correlations, initial embedding, and arbitrary AC or DC biases. In this form, the framework yields time-dependent densities, currents, dipole moments, and nonequilibrium spectral functions for open correlated conductors [1006.2912].

The same structure extends to explicitly open quantum systems coupled to thermal environments. For fermions in an attractive one-dimensional potential interacting with a bosonic heat bath through elastic \(2\to2\) collisions, the nonequilibrium Kadanoff–Baym equations describe equilibration, thermalization with the bath, and decoherence of off-diagonal density-matrix elements. The retarded propagator simultaneously determines the evolving spectral properties of bound and scattering states [2308.07659]. A later extension generalizes this construction to arbitrary fermionic or bosonic systems and baths, including three spatial dimensions, different traps, and different interaction potentials, and uses it to study bound-state formation, decay, regeneration, and thermodynamics in open systems [2506.13302].

A different extension reformulates Keldysh theory for Lindbladian dynamics. In that setting, a contour-ordered NEGF can still be defined, but the contour Hamiltonian acquires branch-dependent non-Hermitian and cross-branch terms. The result is a generalized Kadanoff–Baym framework for interacting systems with dissipative Lindblad evolution, preserving a diagrammatic many-body structure while accommodating loss and gain processes [2402.10824].

Together these formulations show that the Kadanoff–Baym equations are not restricted to isolated Hamiltonian dynamics. They also organize embedding, reservoirs, heat baths, and even Lindbladian dissipation within closely related two-time structures.

## 5. Memory kernels, computational complexity, and algorithmic reformulations

The main numerical obstacle is the memory kernel. On a time grid with \(N_t\) points, the direct two-time formulation stores \(G(t_n,t_{n'})\) on an \(N_t\times N_t\) mesh, with memory scaling as \(O(N_t^2)\), while evaluation of the history integrals leads to a computational cost that typically scales as \(O(N_t^3)\) [2112.05670, 1312.0214]. This is the canonical cost barrier of full Kadanoff–Baym propagation.

Several algorithmic responses have been developed. One is an auxiliary-Hamiltonian representation for spatially local self-energies. There the nonequilibrium Dyson equation is mapped to a noninteracting auxiliary Hamiltonian with explicit bath orbitals chosen so that the bath hybridization reproduces the original self-energy kernel. The original non-Markovian dynamics is then recast as a time-local evolution in an enlarged Hilbert space, avoiding explicit memory-integral evaluation and substantially reducing storage demands in nonequilibrium DMFT applications [1312.0214].

A second strategy is direct memory truncation. If the self-energy decays beyond a finite correlation time \(t_c\), one may truncate the memory integrals to \(|t-t'|\le t_c\). In the truncation scheme analyzed for the Hubbard model, the cost drops from cubic to linear in \(t_{\max}\) for fixed \(t_c\), and simulation times up to two orders of magnitude longer become accessible in DMFT calculations. The truncation is especially well controlled when the self-energy is local or momentum independent [2112.05670].

A third direction explores global-in-time iterative solvers. Fixed-point, Jacobian-free, and Newton–Krylov variants have been tested in nonequilibrium DMFT. Several global methods remain stable at large propagation times, but a standard forward fixed-point iteration does not. Even for the stable methods, the number of iterations required for fixed time step and target accuracy scales roughly linearly with the number of time steps, a behavior associated with a propagating front in the residual error. This identifies a concrete obstacle to making global solvers competitive with causal time-stepping schemes [2512.11371].

These developments do not alter the formal content of the equations. They address how the same two-time physics can be represented or approximated more efficiently.

## 6. Reduced descriptions, practical scope, and known artifacts

The generalized Kadanoff–Baym ansatz replaces the full two-time evolution of \(G^{<,>}(t_1,t_2)\) by a reconstruction from equal-time densities and retarded/advanced propagators. In the scenarios tested by Reeves and collaborators, the terms neglected when deriving the GKBA are explicitly computed and shown to be orders of magnitude smaller than the terms retained, so that they provide only a small correction in the full Kadanoff–Baym equations. In the same study, both GKBA and Kadanoff–Baym calculations capture the dynamics of interacting systems with moderate and even strong interactions well [2306.06225].

More drastic reductions lead toward kinetic theory. In first-principles carrier dynamics for bulk silicon, the generalized Baym–Kadanoff ansatz combined with the complete collision approximation yields a scattering term similar to that of the semiclassical Boltzmann equation [1503.00866]. At a more formal level, analytical WKB solutions of the Kadanoff–Baym equations show that a generalized Boltzmann equation emerges whenever the WKB approximation holds, including far from equilibrium and in time-dependent backgrounds [1202.1301].

At the same time, the full two-time framework has known pathologies when combined with approximate self-energies. In finite Hubbard clusters, self-consistent many-body approximations can generate correlation-induced damping and steady states under strong driving, even for isolated finite systems where such damping is purely artificial. The same study also finds that, for isolated clusters, the steady state reached is not unique but depends on how the external field is switched on. This behavior is tied to approximate self-energies based on infinite partial summations [1004.3662].

A plausible implication is that the practical value of the Kadanoff–Baym equations depends on two distinct choices: whether the two-time structure must be retained, and whether the chosen self-energy respects the microscopic constraints relevant to the observable of interest. The formalism is exact at the contour level; its successes and failures in applications are governed by the approximation used for \(\Sigma[G]\) and by how the memory structure is represented numerically.

Source: https://www.emergentmind.com/topics/kadanoff-baym-equations