Abstract We apply a computationally efficient approach to study the time- and energy-resolved spectral properties of a two-site Hubbard model using the nonequilibrium Green’s function formalism. By employing the iterative generalized Kadanoff-Baym ansatz ( i GKBA) within a time-linear framework, we avoid the computational cost of solving the full two-time Kadanoff-Baym equations. Spectral information is extracted by coupling the system to multiple narrow-band leads, establishing a direct analogy to photoemission experiments. Our results reveal correlation-induced shifts and broadening of spectral features, along with a suppression of transient current oscillations. This approach provides a promising avenue for analyzing correlated electron dynamics in open quantum systems.
Abstract We apply a computationally efficient approach to study the time- and energy-resolved spectral properties of a two-site Hubbard model using the nonequilibrium Green’s function formalism. By employing the iterative generalized Kadanoff-Baym ansatz ( i GKBA) within a time-linear framework, we avoid the computational cost of solving the full two-time Kadanoff-Baym equations. Spectral information is extracted by coupling the system to multiple narrow-band leads, establishing a direct analogy to photoemission experiments. Our results reveal correlation-induced shifts and broadening of spectral features, along with a suppression of transient current oscillations. This approach provides a promising avenue for analyzing correlated electron dynamics in open quantum systems.
Electron dynamics in a two-sites Hubbard model is studied using the nonequilibrium Green's function approach. The study is motivated by the empirical observation that a full solution of the integro-differential Kadanoff-Baym equation (KBE) is more stable and often accompanied by artificial damping [Marc Puig von Friesen, C. Verdozzi, and C.-O. Almbladh (2009)] than its time-linear reformulations relying on the generalized Kadanoff-Baym ansatz (GKBA). Additionally, for conserving theories, numerical simulations suggest that KBE produces natural occupations bounded by one and zero in agreement with the Pauli exclusion principle, whereas, in some regimes, GKBA-based theories violate this principle. As the first step for understanding these issues, the electron dynamics arising in the adiabatic switching scenario is studied. Many-body approximations are classified according to the channel of the Bethe-Salpeter equation in which electronic correlations are explicitly treated. They give rise to the so-called second Born, T-matrix and GW approximations. In each of these cases, the model is reduced to a system of ordinary differential equations, which resemble equations of motion for a driven harmonic oscillator with time-dependent frequencies. A more complete treatment of electronic correlations is achieved by combining different correlation channels, with parquet theory serving as a starting point.
Electron dynamics in a two‐sites Hubbard model is studied using the nonequilibrium Green's function approach. The study is motivated by the empirical observation that a full solution of the integro‐differential Kadanoff–Baym equation (KBE) is more stable and often accompanied by artificial damping [Marc Puig von Friesen, C. Verdozzi, and C.‐O. Almbladh (2009)] than its time‐linear reformulations relying on the generalized Kadanoff–Baym ansatz (GKBA). Additionally, for conserving theories, numerical simulations suggest that KBE produces natural occupations bounded by one and zero in agreement with the Pauli exclusion principle, whereas, in some regimes, GKBA‐based theories violate this principle. As the first step for understanding these issues, the electron dynamics arising in the adiabatic switching scenario is studied. Many‐body approximations are classified according to the channel of the Bethe–Salpeter equation in which electronic correlations are explicitly treated. They give rise to the so‐called second Born, T‐matrix, and GW approximations. In each of these cases, the model is reduced to a system of ordinary differential equations, which resemble equations of motion for a driven harmonic oscillator with time‐dependent frequencies. A more complete treatment of electronic correlations is achieved by combining different correlation channels, with parquet theory serving as a starting point.
The Hubbard model is a prototype for strongly correlated many-particle systems, including electrons in condensed matter and molecules, as well as for fermions or bosons in optical lattices. While the equilibrium properties of these systems have been studied in detail, the nonequilibrium dynamics following a strong nonperturbative excitation only recently came into the focus of experiments and theory. It is of particular interest how the dynamics depend on the coupling strength and on the particle number and whether there exist universal features in the time evolution. Here, we present results for the dynamics of finite Hubbard clusters based on a self-consistent nonequilibrium Green functions (NEGF) approach invoking the generalized Kadanoff-Baym ansatz (GKBA). We discuss the conserving properties of the GKBA with Hartree-Fock propagators in detail and present a generalized form of the energy conservation criterion of Baym and Kadanoff for NEGF. Furthermore, we demonstrate that the HF-GKBA cures some artifacts of prior two-time NEGF simulations. Besides, this approach substantially speeds up the numerical calculations and thus presents the capability to study comparatively large systems and to extend the analysis to long times allowing for an accurate computation of the excitation spectrum via time propagation. Our data obtained within the second Born approximation compare favorably with exact diagonalization results (available for up to 13 particles) and are expected to have predictive capability for substantially larger systems in the weak-coupling limit.
We present a systematic comparison of the real-time Dyson expansion (RTDE) with established non-equilibrium Green's function approaches for simulating driven, interacting quantum systems. Focusing on density matrix dynamics, time-off-diagonal Green's functions, and time-resolved photoemission spectra, we benchmark RTDE against fully self-consistent Kadanoff-Baym equation (KBE) calculations, the generalized Kadanoff-Baym ansatz (GKBA), and exact diagonalization for small systems using second order many-body perturbation theory. Using a driven two-band Hubbard model, we show that mean-field single particle density matrix trajectories provide a reliable baseline for RTDE across a broad range of interaction strengths and excited-carrier populations. Further, RTDE accurately captures correlation effects in the Green's functions, including long-lived oscillations and revivals that are strongly suppressed by the overdamping inherent to self-consistent KBE schemes. As a consequence, RTDE resolves rich non-equilibrium spectral structure in time-resolved photoemission, such as interaction- and population-dependent quasiparticle splittings and bandgap renormalization, which are largely washed out in self-consistent approaches, yet are present in the exact solutions. Our results demonstrate that RTDE bridges the gap between mean-field propagation and full two-time KBE simulations, retaining favorable linear scaling while capturing essential dynamical correlations relevant for ultrafast spectroscopy.
We investigate the spectral properties of an open interacting system by solving the Generalized Kadanoff-Baym Ansatz (GKBA) master equation for the single-particle density matrix, namely the time-diagonal lesser Green's function. To benchmark its validity, we compare the solution obtained within the GKBA with the solution of the Dyson equation (equivalently the full Kadanoff-Baym equations). In both approaches, we treat the interaction within the self-consistent second-order Born approximation, whereas the GKBA still retains the retarded propagator calculated at the Hartree-Fock level. We consider the case of two leads connected through a central correlated region where particles can interact and exploit the stationary particle current at the boundary of the junction as a probe of the spectral features of the system. In this work, as an example, we take the central region to be a one-dimensional quantum wire and a two-dimensional carbon nanotube and show that the solution of the GKBA master equation well captures their spectral features. Our result demonstrates that, even when the propagator used is at the Hartree-Fock level, the GBKA solution retains the main spectral features of the self-energy used.
Nonequilibrium heat transport in quantum systems coupled to wide-band embeddings provides a striking example of the limitations of the generalized Kadanoff-Baym ansatz (GKBA), while solving the full two-time Kadanoff-Baym equations remains computationally prohibitive. To address this challenge, we propose an iterated solution to the reconstruction problem, resulting in a time-linear evolution scheme involving 14 correlators for systems with finite bandwidth embeddings. This approach eliminates artifacts related to the GKBA and resolves convergence issues that arise in the wide-band limit. Furthermore, it enables the calculation of energy- and time-resolved currents, facilitating the modeling of heat flows in quantum systems and energy- and time-resolved photoemission experiments, all at significantly reduced computational cost.