In the nonequilibrium Green's function approach, the approximation of the correlation self-energy at the second-Born level is of particular interest, since it allows for a maximal speed-up in computational scaling when used together with the generalized Kadanoff-Baym ansatz for the Green's function. The present day numerical time-propagation algorithms for the Green's function are able to tackle first principles simulations of atoms and molecules, but they are limited to relatively small systems due to unfavorable scaling of self-energy diagrams with respect to the basis size. We propose an efficient computation of the self-energy diagrams by using tensor-contraction operations to transform the internal summations into functions of external low-level linear algebra libraries. We discuss the achieved computational speed-up in transient electron dynamics in selected molecular systems.
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 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.
While the use of Green's function techniques has a long tradition in quantum chemistry, the possibility of propagating the Kadanoff-Baym equations has remained largely unexplored. We have implemented the time-propagation for atoms and diatomic molecules, starting from a system in the groundstate. The initial stage of the calculation requires solving the Dyson equation self-consistently for the equilibrium Green's function. This Green's function contains a huge amount of information, and we have found it particularly interesting to compare the self-consistent total energies to the results of variational energy functionals of the Green's function. We also use time-propagation for calculating linear response functions, as a means for obtaining the excitation energies of the system. We have presently implemented the propagation for the second Born approximation, while the GW approximation has now been implemented for the ground state calculations.
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.
Simulations of interacting electrons and bosons out of equilibrium, starting from first principles and aiming at realistic multiscale scenarios, is a grand theoretical challenge. Here, using the formalism of nonequilibrium Green's functions and relying in a crucial way on the recently discovered time-linear formulation of the Kadanoff-Baym equations, we present a versatile toolbox for the simulation of correlated electron-boson dynamics. A large class of methods are available, from the Ehrenfest to the dressed
GD for the treatment of electron-boson interactions in combination with perturbative, i.e., Hartree-Fock and second-Born, or nonperturbative, i.e.,
GW and
T matrices either without or with exchange effects, for the treatment of the Coulomb interaction. In all cases the numerical scaling is linear in time and the equations of motion satisfy all fundamental conservation laws.