The Generalized Kadanoff--Baym Ansatz with mean-field propagators is increasingly used to simulate non-equilibrium quantum fermionic and bosonic systems. Compared with the full Kadanoff--Baym equations, its mean-field propagators substantially reduce computational cost, enabling simulations with computational effort that scales linearly with propagation time. However, the time-diagonal structure of the Generalized Kadanoff--Baym Ansatz obscures spectral properties within the collision integral of the transport equation. Here, we recover and investigate these hidden spectral properties using the extended Generalized Kadanoff--Baym Ansatz. For a Hubbard cluster, we compare ground-state spectral functions obtained with the extended Generalized Kadanoff--Baym Ansatz against those from standard Kadanoff--Baym equations. At moderate interaction strengths, the spectral functions show good agreement. At large interaction strengths, however, significant deviations emerge, and the spectral function obtained with the extended Generalized Kadanoff--Baym Ansatz becomes negative. These results demonstrate that, despite its computational advantages, the Generalized Kadanoff--Baym Ansatz has important limitations in strongly interacting regimes and may produce unphysical spectral properties when interactions become sufficiently large.
A systematic and unambiguous method of deriving generalized transport equations, i.e., equations for distribution functions having a single-time structure, on the basis of the nonequilibrium Green function is obtained, if the common Kadanoff-Baym (KB) ansatz is replaced by a modification which we call the generalized KB ansatz. This new ansatz is fully consistent with the dynamical structure of the theory and is independent of any specific representation. The resulting equations appear to be the zeroth-order approximation of a systematic expansion in terms of the collision duration. In the case of the electron transport in a strong homogeneous electric field, the generalized ansatz is shown to be in agreement with the superoperator methods, whereas the KB ansatz is known to fail in this case.
For a minimal Hubbard-type system at different interaction strengths U , we investigate the density-response for an excitation beyond the linear regime using the generalized Kadanoff-Baym ansatz (GKBA) and the second Born (2B) approximation. We find strong correlation features in the response spectra and establish the connection to an involved double excitation process. By comparing approximate and exact Green's function results, we also observe an anomalous U -dependence of the energy of this double excitation in 2B+GKBA. This is in accordance with earlier findings [K. Balzer et al., EPL 98 , 67002 (2012)] on double excitations in quantum wells.
The time evolution in quantum many-body systems after external excitations is attracting high interest in many fields, including dense plasmas, correlated solids, laser-excited materials, or fermionic and bosonic atoms in optical lattices. The theoretical modeling of these processes is challenging, and the only rigorous quantum-dynamics approach that can treat correlated fermions in two and three dimensions is nonequilibrium Green functions (NEGF). However, NEGF simulations are computationally expensive due to their
T3 scaling with the simulation duration
T. Recently,
T2 scaling was achieved with the generalized Kadanoff-Baym ansatz (GKBA), for the second-order Born (SOA) self energy, which has substantially extended the scope of NEGF simulations. In a recent Letter [Schl\"unzen et al., Phys. Rev. Lett. 124, 076601 (2020).] we demonstrated that GKBA-NEGF simulations can be efficiently mapped onto coupled time-local equations for the single-particle and two-particle Green functions on the time diagonal, hence the method has been called the G1-G2 scheme. This allows one to perform the same simulations with order
T1 scaling, both for SOA and
GW self energies giving rise to a dramatic speedup. Here we present more details on the G1-G2 scheme, including derivations of the basic equations including results for a general basis, for Hubbard systems, and for jellium. Also, we demonstrate how to incorporate initial correlations into the G1-G2 scheme. Further, the derivations are extended to a broader class of self energies, including the
T matrix in the particle-particle and particle-hole channels and the dynamically-screened-ladder approximation. Finally, we demonstrate that, for all self energies, the CPU-time scaling of the G1-G2 scheme with the basis dimension
Nb can be improved compared to our first report: The overhead compared to the original GKBA is not more than an additional factor
Nb, even for Hubbard systems.
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.
In the non-equilibrium Green function calculations, the use of the generalized Kadanoff–Baym ansatz (GKBA) allows for a simple approximate reconstruction of the two-time Green function from its time-diagonal value. With this, a drastic reduction of the computational needs is achieved in time-dependent calculations, making longer time propagation possible and more complex systems accessible. This paper gives credit to the GKBA that was introduced 25 years ago. After a detailed derivation of the GKBA, we recall its application to homogeneous systems and show how to extend it to strongly correlated, inhomogeneous systems. As a proof of concept, we present the results for a two-electron quantum well, where the correct treatment of the correlated electron dynamics is crucial for a correct description of the equilibrium and dynamic properties.
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.