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.
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.
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. ACKNO
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.
The method of nonequilibrium Green functions (NEGFs) constitutes a solid framework to describe closed and open quantum systems. It captures situations far from equilibrium including arbitrary perturbations and correlations - effects higher than first order in the binary interactions. On the contrary, in the absence of external fields, it reduces to the imaginary-time Matsubara formalism. The abilities (and successes) of the method originate from its self-consistency and the conserving character of many-body approximations that are introduced by diagram technique. On this account, NEGFs have, now for a long time, brought insight into the dynamics of nuclear matter, non-ideal quantum gases and (semiconductor) plasmas. However, applications to localized, finite and strongly interacting systems have emerged only within the recent decade. The dynamical properties follow from quantum kinetic equations - the two-time Kadanoff-Baym equations for the one-particle NEGF - the numerical solution of which is the central topic of this thesis. The non-Markovian structure of these equations thereby inhibits an unlimited, gradual propagation of the nonequilibrium Green function: Any self-consistent, future dynamics of the system is tightly linked to the (complete) past by non-local self-energies. Further, the success in obtaining solutions is being diminished if the system under investigation is not isotropic or inhomogeneous and, in addition, the spatial topology (e.g., the charge carrier or spin density) needs to be resolved in time. On top of this, small isolated systems are generally not candidates showing fast decoherence which offers little hope for smooth and well-damped integration kernels. Despite these difficulties, the thesis at hand is devoted to the computation of the two-time NEGF for finite, closed systems of which ensembles of electrons bound by nuclei forming atoms or molecules are the most natural representatives. Here, in the spirit of the well-studied one-dimensional helium atom, the electron dynamics in model atoms and molecules is addressed in ab initio fashion. To render calculations possible in the first place, (i) a novel hybrid representation of the NEGF is introduced based on the conception of finite elements (FE) and the discrete variable representation (DVR), and (ii) a parallel algorithm is developed that allows for an efficient, well scalable distributed-memory computation of the full two-time Green function. The first point attributes plenty of flexibility and makes the method readily adjustable to different Hamiltonians avoiding the need to, by hand or in intricate numerical manner, calculate the whole bunch of matrix elements of one- and two-particle operators. Also, in consequence, it drastically simplifies the evaluation of self-energies. The second point makes the NEGF approach ready for large-scale, high-performance computing facilities, the resources (of computer power: mainly memory!) of which are indispensable for (some of) the calculations presented. The applicability of the developed methods is demonstrated inter alia on the basis of helium, beryllium, hydrogen and lithium hydride modeled in one spatial dimension whereby electron-electron correlations are treated in the second(-order) Born approximation. After a thorough discussion on the self-consistent ground states, the thesis focuses on the electron dynamics in the linear response regime as well as in the presence of strong external fields. Overall, to assess the the performance of the second Born approximation, the results are compared to the Hartree-Fock approximation, on the one hand, and to the exact solution of the time-(in)dependent Schrödinger equation, on the other hand. In addition, aiming at the possibility of extended time-dependent calculations, the behavior of the generalized Kadanoff-Baym ansatz is discussed.