We have developed a time-propagation scheme for the Kadanoff-Baym equations for general inhomogeneous systems. These equations describe the time evolution of the nonequilibrium Green function for interacting many-body systems in the presence of time-dependent external fields. The external fields are treated nonperturbatively whereas the many-body interactions are incorporated perturbatively using Phi-derivable self-energy approximations that guarantee the satisfaction of the macroscopic conservation laws of the system. These approximations are discussed in detail for the time-dependent Hartree-Fock, the second Born, and the GW approximation.
We have developed a time-propagation scheme for the Kadanoff-Baym equations for general inhomogeneous systems. These equations describe the time evolution of the nonequilibrium Green function for interacting many-body systems in the presence of time-dependent external fields. The external fields are treated nonperturbatively whereas the many-body interactions are incorporated perturbatively using Phi-derivable self-energy approximations that guarantee the satisfaction of the macroscopic conservation laws of the system. These approximations are discussed in detail for the time-dependent Hartree-Fock, the second Born, and the GW approximation.
We implement time propagation of the nonequilibrium Green function for atoms and molecules by solving the Kadanoff-Baym equations within a conserving self-energy approximation. We here demonstrate the usefulness of time propagation for calculating spectral functions and for describing the correlated electron dynamics in a nonperturbative electric field. We also demonstrate the use of time propagation as a method for calculating charge-neutral excitation energies, equivalent to highly advanced solutions of the Bethe-Salpeter equation.
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.
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 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.
We have implemented time-propagation of the non-equilibrium Green function for atoms and molecules, by solving the Kadanoff-Baym equations within a conserving self-energy approximation. We here demonstrate the usefulnes of time-propagation for calculating spectral functions and for describing the correlated electron dynamics in a non-perturbative electric field. We also demonstrate the use of time-propagation as a method for calculating charge-neutral excitation energies, equivalent to highly advanced solutions of the Bethe-Salpeter equation.
We present in detail a method we recently introduced [Phys. Rev. Lett. 103, 176404 (2009)] to describe finite systems in and out of equilibrium, where the evolution in time is performed via the Kadanoff-Baym equations within many-body perturbation theory. Our systems consist of small, strongly correlated clusters, described by a Hubbard Hamiltonian within the Hartree-Fock, second Born, GW, and T-matrix approximations. We compare the results from the Kadanoff-Baym dynamics to those from exact numerical solutions. The outcome of our comparisons is that, among the many-body schemes considered, the T-matrix approximation is superior at low electron densities while none of the tested approximations stands out at half filling. Such comparisons permit a general assessment of the whole idea of applying many-body perturbation theory, in the Kadanoff-Baym sense, to finite systems. A striking outcome of our analysis is that when the system evolves under a strong external field, the Kadanoff-Baym equations develop a steady-state solution as a consequence of a correlation-induced damping. This damping is present both in isolated (finite) systems, where it is purely artificial, as well as in clusters contacted to (infinite) macroscopic leads. The extensive numerical characterization we performed indicates that this behavior is present whenever approximate self-energies, which include correlation effects, are used. Another important result is that, for isolated clusters, the steady state reached is not unique but depends on how one switches on the external field. When the clusters are coupled to macroscopic leads, one may reach multiple quasisteady states with arbitrarily long lifetimes.