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.
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.
The quantum dynamics of fermionic or bosonic many-body systems following external excitation can be successfully studied using two-time nonequilibrium Green's functions (NEGF) or single-time reduced density matrix methods. Approximations are introduced via a proper choice of the many-particle self-energy or decoupling of the BBGKY hierarchy. These approximations are based on Feynman's diagram approaches or on cluster expansions into single-particle and correlation operators. Here, we develop a different approach where, instead of equations of motion for the many-particle NEGF (or density operators), single-time equations for the correlation functions of fluctuations are analyzed. We present a derivation of the first two equations of the alternative hierarchy of fluctuations and discuss possible decoupling approximations. In particular, we derive the polarization approximation (PA) which is shown to be equivalent to the single-time version [following by applying the generalized Kadanoff-Baym ansatz (GKBA)] of the nonequilibrium GW approximation with exchange effects of NEGF theory, for weak coupling. The main advantage of the quantum fluctuations approach is that the standard ensemble average can be replaced by a semiclassical average over different initial realizations, as was demonstrated before by Lacroix and co-workers [see e.g. D. Lacroix et al., Phys. Rev. B, 2014, 90, 125112]. Here, we introduce the stochastic GW (SGW) approximation and the stochastic polarization approximation (SPA) which are demonstrated to be equivalent to the single-time GW approximation without and with exchange, respectively, in the weak coupling limit. In addition to the standard stochastic approach to sample initial configurations, we also present a deterministic approach. Our numerical tests confirm that our approach has the same favorable linear scaling with the computation time as the recently developed G1-G2 scheme [Schluenzen et al., Phys. Rev. Lett., 2020, 124, 076601]. At the same time, the SPA and SGW approaches scale more favorably with the system size than the G1-G2 scheme, allowing to extend nonequilibrium GW calculations to bigger systems.
The theory of nonequilibrium Green functions (NEGF) has seen a rapid development over the recent three decades. Applications include diverse correlated many‐body systems in and out of equilibrium. Very good agreement with experiments and available exact theoretical results could be demonstrated if the proper selfenergy approximations were used. However, full two‐time NEGF simulations are computationally costly, as they suffer from a cubic scaling of the computation time with the simulation duration. Recently the G1–G2 scheme that exactly reformulates the generalized Kadanoff–Baym ansatz with Hartree–Fock propagators (HF‐GKBA) into time‐local equations is introduced, which achieves time‐linear scaling and allows for a dramatic speedup and extension of the simulations (Schluenzen et al. Phys. Rev. Lett. 2020, 124, 076601). Remarkably, this scaling is achieved quickly, and also for high‐level selfenergies, including the nonequilibrium GW and T‐matrix approximations (Joost et al. Phys. Rev. B 2020, 101, 245101). Even the dynamically screened ladder approximation is now feasible (Joost et al. Phys. Rev. B 2022, 105, 165155), and also applications to electron‐boson systems are demonstrated. Herein, an overview on recent results that are achieved with the G1–G2 scheme is presented. Problems and open questions are discussed and further ideas of how to overcome the current limitations of the scheme and present are presented. The G1–G2 scheme is illustrated by presenting applying it to the excitation dynamics of Hubbard clusters, to optical excitation of graphene, and to charge transfer during stopping of ions by correlated materials.
The Hartree-Fock generalized Kadanoff-Baym ansatz (HF-GKBA) offers an approximate numerical procedure for propagating the two-time nonequilibrium Green's function (NEGF). Here, using the
GW self-energy, we compare the HF-GKBA to exact results for a variety of systems with long- and short-range interactions, different two-body interaction strengths, and various nonequilibrium preparations. We find excellent agreement between the HF-GKBA and exact time evolution in models when more realistic long-range exponentially decaying interactions are considered. This agreement persists for long times and for intermediate to strong interaction strengths. In large systems, HF-GKBA becomes prohibitively expensive for long-time evolutions. For this reason, we look at the use of dynamical mode decomposition (DMD) to reconstruct long-time NEGF trajectories from a sample of the initial trajectory. Using no more than 16% of the total time evolution, we reconstruct the total trajectory with high fidelity. Our results show the potential for DMD to be used in conjunction with HF-GKBA to calculate long-time trajectories in large-scale systems.
The dynamics of strongly correlated fermions following an external excitation reveals extremely rich collective quantum effects. Examples are fermionic atoms in optical lattices, electrons in correlated materials, and dense quantum plasmas. Presently, the only quantum-dynamics approach that rigorously describes these processes in two and three dimensions is the nonequilibrium Green functions (NEGF) method. However, NEGF simulations are computationally expensive due to their T^{3} scaling with the simulation duration T. Recently, T^{2} scaling was achieved with the generalized Kadanoff-Baym ansatz (GKBA), for second-order Born (SOA) selfenergies, which has substantially extended the scope of NEGF simulations. Here we demonstrate that GKBA-NEGF simulations can be performed with order T^{1} scaling, both for SOA and GW selfenergies, and point out the remarkable capabilities of this approach.