---
title: 'Time-Dependent NEGF: Real-Time Quantum Transport'
url: https://www.emergentmind.com/topics/time-dependent-negf-td-negf
type: topic
---

# Time-Dependent NEGF: Real-Time Quantum Transport

Time-dependent nonequilibrium Green's function (TD-NEGF) theory is a quantum many-body framework that enables the real-time description of electronic, spin, and excitonic transport in open systems subjected to arbitrary time-dependent perturbations. TD-NEGF generalizes the conventional steady-state NEGF formalism to problems far from equilibrium, allowing for the treatment of transient dynamics, ultrafast response, quantum pumping, and dynamical functional phenomena in nanostructures, molecular junctions, and mesoscale devices. The approach leverages contour-ordered Green's functions, self-energy techniques, and memory-kernel equations of motion to capture the influence of leads, interactions, and external drives, and incorporates both transient and asymptotic steady-state behavior [1505.03567, 2201.02646, 1707.03524].

## 1. Fundamental Formalism: Contour-Ordered and Real-Time Green’s Functions

TD-NEGF theory is fundamentally built on the single-particle (or, more generally, many-body) Green's function defined on the Keldysh contour:
\[
G(\tau, \tau') = -i \langle T_C[\psi_H(\tau)\psi_H^\dagger(\tau')]\rangle
\]
where \( T_C \) orders operators along the Keldysh time contour γ, which runs from the initial time to infinity (forward branch), then back to the initial time (backward branch), closed with a descent along the imaginary axis for finite temperature [2201.02646].

Physical quantities such as particle density, correlation functions, and currents are obtained by projecting \( G(\tau,\tau') \) onto the real time branches, yielding the retarded (\( G^R \)), advanced (\( G^A \)), lesser (\( G^< \)), and greater (\( G^> \)) components. For noninteracting and interacting systems, these have explicit operator definitions (see, e.g., [1707.03524]):
\[
\begin{aligned}
G^R(t, t') &= -i\Theta(t-t')\langle[\psi(t),\psi^\dagger(t')]\rangle \\
G^<(t, t') &= i\langle \psi^\dagger(t')\psi(t) \rangle
\end{aligned}
\]
The local density is obtained as \( n(x, t) = -i G^<(x, t; x, t) \) [1505.03567].

## 2. Kadanoff–Baym and Dyson Equations in the Time Domain

The time evolution of the system is governed by the integro-differential Kadanoff–Baym equations (KBEs), which, for the single-particle sector, read:
\[
\begin{aligned}
\left[i \partial_{t_1} - H(t_1)\right] G^R(t_1, t_2) &= \delta(t_1-t_2) + \int_{t_2}^{t_1} \Sigma^R(t_1, t') G^R(t', t_2) dt' \\
G^<(t_1, t_2) &= \int_{-\infty}^{t_1} dt \int_{-\infty}^{t_2} dt' G^R(t_1, t) \Sigma^<(t, t') [G^R(t_2, t')]^\dagger
\end{aligned}
\]
Here, \( H(t) \) represents the one-particle part of the Hamiltonian (including time-dependent and interaction terms), and the self-energies \( \Sigma^{R,<} \) encode the effect of reservoir coupling (embedding) and many-body interactions [1505.03567, 2201.02646]. For interacting systems, the self-energy splits as \( \Sigma = \Sigma_{em} + \Sigma_{MB} \) with embedding and many-body contributions [2201.02646].

In practical partitioning schemes, only the central region is treated explicitly, with the effect of semi-infinite leads entering through boundary self-energies localized to the edges, ensuring exact transparent boundary conditions and nonreflected wavefunctions [1505.03567, 1307.6419].

## 3. Embedding Self-Energy, Initial Condition, and Observables

The embedding self-energy \( \Sigma_{em} \) provides a nonperturbative account of the reservoirs:
\[
\Sigma_{em,\alpha}(t, t') = h_{C\alpha}(t) g_{\alpha\alpha}(t, t') h_{\alpha C}(t')
\]
where \( h_{C\alpha} \) is the tunneling Hamiltonian, and \( g_{\alpha\alpha} \) is the Green’s function of an isolated lead α [2201.02646].

Transmission, reflection, and local observables are computed from \( G^{<} \) or associated wave functions. For example, the local current between sites \( i \) and \( i+1 \):
\[
J_{i+1/2}(t) = \frac{1}{m a} \mathrm{Im}\left[\Psi_i^*(t)\Psi_{i+1}(t)\right]
\]
with incoming and outgoing currents yielding the time-resolved transmission \( T(t) \) [1505.03567]. Physical charge and energy currents in general open systems follow the Meir–Wingreen formula, which in time domain reads:
\[
I_\alpha(t) = 2\,\mathrm{Re}\,\mathrm{Tr}_C \int_{t_0}^t d\tau \, [\Sigma^<_{em,\alpha}(t, \tau) G^A(\tau, t) + \Sigma^R_{em,\alpha}(t, \tau) G^<(\tau, t)]
\]
[2201.02646, 1707.03524]. The correspondence to the steady-state (Landauer–Büttiker) limit is recovered as \( t_0 \to -\infty \) and drive potentials become stationary [1412.4118, 2201.02646].

## 4. Numerical Implementation and Computational Schemes

TD-NEGF calculations require efficient handling of memory kernels and matrix-valued propagation:
- **Spatial discretization**: The device region is discretized into a finite grid or basis \( x_i \), with kinetic terms represented by finite-difference or tight-binding Hamiltonians [1505.03567].
- **Memory kernels**: The integral terms arising from lead self-energies are history-dependent and can be re-expressed via auxiliary variables or reduced using pole decompositions (Padé, Lorentzian) for spectral densities and Fermi functions, allowing linear-scaling algorithms in the time domain [2211.15635, 1304.0192, 1810.02577, 1907.05252].
- **Wide-band limit (WBL)**: In many applications, the embedding self-energies are approximated as energy-independent (constant broadening \( \Gamma_\alpha \)), greatly simplifying the time-evolution equations to closed forms [1412.4118, 1810.02577].
- **ODE formulations**: For large-scale simulations, the generalized Kadanoff–Baym ansatz (GKBA) allows the equations to be rewritten as a hierarchy of ODEs for single-time quantities (density matrices and embedding correlators), achieving linear scaling in time [2211.15635, 1907.05252].

Hierarchical equations of motion (HEOM) and pole-fitting schemes are leveraged to handle complex and structured spectral densities beyond the WBL [1304.0192].

## 5. Extensions: Interactions, Correlations, and Multiphasic Physics

TD-NEGF accommodates many-body interactions via diagrammatic self-energies (Hartree–Fock, second Born, GW, GD, T-matrix, Fan–Migdal) [2201.02646, 2211.15635]. The collision integral arising from interactions can be closed using small sets of ODEs involving higher-order correlation matrices. The GKBA allows for conserving dynamics at the level of the one-particle density matrix. Electron-electron, electron-phonon, and even electron-photon couplings can be incorporated in this framework.

In multiscale or hybrid schemes, TD-NEGF can be coupled to classical equations of motion for slow variables, e.g., spin dynamics via Landau–Lifshitz–Gilbert equations in spintronics [1802.05682].

Transient charge transfer and Hartree corrections are treated efficiently by coupling the quantum equations to analytical Poisson solvers (e.g., muffin-tin networks), ensuring charge conservation under strong time-dependent driving [1810.02577, 1907.05252].

## 6. Applications: Ultrafast Scattering, Quantum Pumping, and Device Physics

TD-NEGF is employed to simulate time-resolved quantum transport under arbitrary field protocols, such as voltage pulses, AC or multi-harmonic drives, gate voltages, time-dependent barriers, ultrafast laser fields, and photonic coupling [1505.03567, 2201.02646, 1907.05252].

Specific applications include:
- **Ultrafast electron scattering by laser pulses**: TD-NEGF captures dynamical polarization, induced density modulations, and the outgoing wave packets generated by pulse-driven scattering [1505.03567].
- **Photon-assisted and AC transport**: Closed-form analytic formulas capture transient and long-time currents, with non-adiabatic features such as ringing, frequency mixing, and sideband generation emerging naturally in the formalism [1412.4118, 1907.05252].
- **Correlated excitation transport**: Wave-packet dynamics of electron-hole pairs, polaronic evolution, and their decoherence can be tracked in real-time [2211.15635].
- **Spin and charge pumping**: TD-NEGF+LLG predicts AC and DC current generation by moving or pulse-driven magnetic domain walls, and quantifies the breakdown of perturbative spin-motive force theory [1802.05682].
- **Molecular and nanoscale devices**: TD-NEGF implemented with muffin-tin Poisson solvers enables full-scale time-resolved simulation including dynamical charge transfer and screening in metal–semiconductor junctions [1810.02577].

## 7. Efficiency, Scalability, and Limitations

Numerical advances in pole expansion, hierarchical ODE solvers, and memory kernel reduction techniques have enabled TD-NEGF to reach system sizes and time scales inaccessible to brute-force two-time integration [2211.15635, 1307.6419, 1810.02577]. For typical nanodevice calculations, runtime is reduced by up to one to two orders of magnitude, and the method remains stable for tens of picoseconds and beyond.

Limitations persist for systems with strong nonequilibrium correlations beyond the single-particle sector, for extremely narrow-band reservoirs (breakdown of WBL), for very large density distortions, or for strong dynamical exchange-correlation effects not captured at the Hartree level [1810.02577, 2211.15635].

Potential extensions under active investigation include the incorporation of fully dynamical exchange-correlation kernels, energy-dependent embedding, and time-dependent quantum capacitance, as well as further integration with stochastic or bath models for strongly open quantum systems [2201.02646, 1810.02577].

---

**Principal References**: [1505.03567], [1707.03524], [2201.02646], [2211.15635], [1412.4118], [1307.6419], [1810.02577], [1907.05252], [1802.05682], [1304.0192].

Source: https://www.emergentmind.com/topics/time-dependent-negf-td-negf