The nonequilibrium interplay of electrons and phonons plays an important role in the thermalization of solids, yet the microscopic picture of transient states and relaxation pathways remains incomplete. Previous nonequilibrium Green's function (NEGF) studies with full two-time dependence were restricted to local phonons and local self-energy approximations, leaving momentum-dependent phonon dynamics largely unexplored. Here, we demonstrate the recently developed quantics-tensor-train (QTT) NEGF framework via large-scale lattice simulations of models with arbitrary phonon dispersions. QTTs provide a memory-efficient representation of full two-time Green's functions, enabling momentum-resolved simulations with full electron-phonon feedback on lattices up to 256x256 sites. Comparing optical and acoustic phonon models, we reveal a hierarchy of relaxation bottlenecks that extends the well-known phonon-window bottleneck effect. For optical phonons, we confirm the main phonon-energy window and uncover an additional reduced window separating momentum-space regions with excess and deficit electronic population. We also identify a separate phonon-thermalization bottleneck rooted in momentum-dependent coupling to the particle-hole continuum. High momentum and frequency resolution makes the phonon--charge correspondence directly visible. For acoustic phonons, the phonon-energy window acquires pronounced momentum dependence dictated by energy and momentum conservation. The reduced window becomes asymmetric, and directional scattering between Brillouin-zone regions produces a persistent bottleneck for low-momentum phonon modes. Our results establish diagrammatic QTT-NEGF simulations as a scalable framework for quantitative nonequilibrium electron-phonon dynamics, overcoming previous lattice-size and propagation-time limitations and providing accurate reference data for time-resolved spectroscopies.
One of the challenges in diagrammatic simulations of nonequilibrium phenomena in lattice models is the large memory demand for storing momentum-dependent two-time correlation functions. This problem can be overcome with the recently introduced quantics tensor train (QTT) representation of multivariable functions. Here, we demonstrate nonequilibrium Green's function simulations within the GW and Migdal approximations with high momentum resolution, up to times which exceed the capabilities of standard implementations and are long enough to study, e.g., transient Floquet physics during multicycle electric field pulses and thermalization dynamics. The self-consistent calculation on the three-leg Kadanoff-Baym contour is fully self-contained, employing only QTT-compressed functions and input functions which are either generated directly in QTT form or obtained via quantics tensor cross interpolation.
Correlated real-time dynamics in large, spatially inhomogeneous quantum systems remain difficult to access with nonequilibrium many-body methods. Two-time nonequilibrium Green functions (NEGF) retain dynamical correlations but their computational runtime grows cubically with the number of time steps
Nt. This scaling bottleneck could recently be overcome by introducing the G1--G2 scheme that is linear in
Nt, but requires propagation of a two-particle correlation function and may suffer from numerical instabilities. This has restricted simulations to small systems with
Nb∼102 basis states. Here we introduce a quantum-fluctuation formulation of nonequilibrium Green functions, denoted
δNEGF, that represents dynamical two-particle correlations through fluctuations of field-operator products,
δG^. This guarantees stable dynamics by preserving the positivity of the reduced density matrices, avoids the explicit storage of the two-particle Green function, and reduces the propagation to a finite ensemble of Hartree-Fock-like trajectories. Combined with a stochastic low-rank decomposition of the correlation functions, the method retains time-linear scaling while extending dynamical
GW and particle-particle and particle-hole
T-matrix simulations to basis sizes of order
Nb∼104. We benchmark
δNEGF against exact and HF-GKBA results for lattice systems, finding stable correlated dynamics also at strong coupling. We further demonstrate large-scale simulations of diffusion in two-dimensional Hubbard lattices and ultrafast relaxation in graphene nanoribbon heterostructures with long-range Coulomb interactions. These results establish
δNEGF as a scalable route to dynamical self-energy simulations of large, spatially inhomogeneous correlated quantum systems beyond the reach of existing NEGF implementations.
We investigate the nonequilibrium dynamics of the weak-coupling Hubbard-Holstein model after a sudden switch-on of the electron-phonon interaction within nonequilibrium dynamical mean-field theory (DMFT). Using the self-consistent Migdal approximation for the electron-phonon coupling together with second-order perturbation theory for the electron-electron interaction, we show that the relaxation dynamics exhibits a crossover between electron-dominated and phonon-dominated regimes, extending to finite Hubbard interaction the scenario previously identified in the Holstein model. To investigate the microscopic buildup of the thermal state, we analyze the dynamics within the Step-by-Step DMFT framework. In the plane of real time and DMFT iteration number, thermalization is marked by a sharp propagating front. This front appears in electronic observables already for weak quenches within the simulated time window, whereas the phononic sector exhibits a visible front only at sufficiently strong coupling. Thus, at weak coupling the local dispersionless phonons show a delayed onset of front formation, while near and beyond the crossover the front develops on comparable timescales in both the electronic and phononic sectors. Whenever both fronts are resolved, they propagate with the same velocity, showing that thermalization spreads coherently through the coupled electron-phonon system.
The nonequilibrium Green's function formalism provides a versatile and powerful framework for numerical studies of nonequilibrium phenomena in correlated many-body systems. For calculations starting from an equilibrium initial state, a standard approach consists of discretizing the Kadanoff-Baym contour and implementing a causal time-stepping scheme in which the self-energy of the system plays the role of a memory kernel. This approach becomes computationally expensive at long times, because of the convolution integrals and the large amount of computer memory needed to store the Green's functions. A recent idea for the compression of nonequilibrium Green's functions is the quantics tensor train representation. Here, we explore this approach by implementing equilibrium and nonequilibrium simulations of the two-dimensional Hubbard model with a second-order weak-coupling approximation to the self-energy. We show that calculations with compressed two-time functions are possible without any loss of accuracy, and that the quantics tensor train implementation shows a much improved scaling of the computational effort and memory demand with the length of the time contour.
This article presents an overview on recent progress in the theory of nonequilibrium Green functions (NEGF). We discuss applications of NEGF simulations to describe the femtosecond dynamics of various finite fermionic systems following an excitation out of equilibrium. This includes the expansion dynamics of ultracold atoms in optical lattices following a confinement quench and the excitation of strongly correlated electrons in a solid by the impact of a charged particle. NEGF, presently, are the only ab initio quantum approach that is able to study the dynamics of correlations for long times in two and three dimensions. However, until recently, NEGF simulations have mostly been performed with rather simple selfenergy approximations such as the second-order Born approximation (SOA). While they correctly capture the qualitative trends of the relaxation towards equilibrium, the reliability and accuracy of these NEGF simulations has remained open, for a long time. Here we report on recent tests of NEGF simulations for finite lattice systems against exact-diagonalization and density-matrix-renormalization-group benchmark data. The results confirm the high accuracy and predictive capability of NEGF simulations-provided selfenergies are used that go beyond the SOA and adequately include strong correlation and dynamical-screening effects. With an extended arsenal of selfenergies that can be used effectively, the NEGF approach has the potential of becoming a powerful simulation tool with broad areas of new applications including strongly correlated solids and ultracold atoms. The present review aims at making such applications possible. To this end we present a selfcontained introduction to the theory of NEGF and give an overview on recent numerical applications to compute the ultrafast relaxation dynamics of correlated fermions. In the second part we give a detailed introduction to selfenergies beyond the SOA. Important examples are the third-order approximation, the [Formula: see text] approximation, the T-matrix approximation and the fluctuating-exchange approximation. We give a comprehensive summary of the explicit selfenergy expressions for a variety of systems of practical relevance, starting from the most general expressions (general basis) and the Feynman diagrams, and including also the important cases of diagonal basis sets, the Hubbard model and the differences occuring for bosons and fermions. With these details, and information on the computational effort and scaling with the basis size and propagation duration, readers will be able to choose the proper basis set and straightforwardly implement and apply advanced selfenergy approximations to a broad class of systems.
We study thermalization dynamics in a fermion-phonon variant of the Sachdev-Ye-Kitaev model coupled to an external cold thermal bath of harmonic oscillators. We find that quantum critical fermions thermalize more efficiently than phonons, in sharp contrast to the behavior in the Fermi liquid regime. In addition, after a short prethermal stage, the system acquires a quasithermal distribution given by a time-dependent effective temperature, reminiscent of ``hydrodynamic'' relaxation. All physical observables relax at the same rate which scales with the final temperature through an exponent that depends universally on the low-energy spectrum of the system and the bath. Such relaxation rate is derived using a hydrodynamic approximation in full agreement with the numerical solution of a set quantum kinetic equations derived from the Keldysh formalism for nonequilibrium Green's functions. Our results hint toward further research on the applicability of the hydrodynamic picture in the description of the late time dynamics of open quantum systems despite the absence of conserved quantities in regimes dominated by conserving collisions.