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 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 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 framework of nonequilibrium Green functions (NEGFs) is developed to accurately describe interacting quantum systems evolving in the presence of strongly exciting ultra-fast fields.At its core stands the selfconsistent determination of the selfenergy, a functional of the single-particle Green function, which incorporates all relevant many-body physics and, via diagrammatic techniques, guarantees the conservation of the constants of motion.Over the last decades, NEGFs have been successfully applied to a broad range of domains ranging from homogeneous systems such as non-ideal quantum gases and plasmas over condensed matter to highly inhomogeneous systems such as atoms and molecules photo-ionized by laser-fields or colliding nuclear matter.Only recently, the first accurate simulation of the expansion dynamics of clouds of ultra-cold fermionic atoms in a two-dimensional optical lattice could be achieved using NEGFs.The single-particle Green function is the solution of a coupled set of kinetic equations.Its intrinsic dependence on two independent time-arguments, ranging over the complete past of the system, not only enables the determination of all singleparticle observables but also provides access to important two-particle observables and spectral information.With this coupling to the history, though, comes the practical necessity to apply physically motivated approximation techniques to the structure of the selfenergy and the Green functions, to achieve numerical tractability.This thesis is devoted to provide a thorough overview of the possible approximation techniques including perturbative and non-perturbative approaches.Starting from the exact coupled set of equations of motion for the Green functions, the expressions for all selfenergy approximations are derived in the most general basis representation, i.e., for interaction tensor elements w ijkl .From this, the corresponding expressions in a diagonal basis, where the interaction can be described as a matrix w ij , and for the special case of the Hubbard model, where the interaction is just given as a scalar U , are deduced.
In this thesis we present methodologies that have been developed to efficiently solve for the non-equilibrium Green’s Function (NEGF) of interacting electron systems, and demonstrate their applications to several different condensed matter systems. These solvers are based on the numerical integration of the equations of motion for the NEGF, which are known as the Kadanoff-Baym Equations (KBE). The KBE are a general set of equations which allow for the study of a wide range of systems that are of interest to condensed matter physicists, including superconductivity, nuclear matter, molecules and atoms, and qubits, among many others. They are also useful for solving impurity models and systems in the dynamical mean-field theory (DMFT) approximation through their ability to incorporate dynamical hybridization functions. The generality of these equations comes from their connection to time-dependent Many-Body Perturbation Theory (MBPT), which introduces non-Markovian dynamical correlations through the self-energy functional. MBPT ensures important conservation laws, such as those for density and energy, are obeyed no matter the type of time-dependent perturbations that are applied to the Hamiltonian. The non-Markovian correlations manifest in a kernel that couples the dynamics to the entire history of the system. This coupling leads to computational costs that become prohibitively expensive as longer-time dynamics are probed, this is especially true for systems which decohere slowly, leading to undamped dynamics and kernels. Due to the large amount of interest in non-equilibrium quantum dynamics, the availability of efficient solvers for the KBE are of great importance. This thesis will focus on the implementation and analysis of integration methods for the KBE, including compression algorithms which lower computational complexity and adaptive integration techniques which allow for resolution of dynamics which vary over wide ranges of timescales. Chapter 1 gives an introduction to quantum systems and the problems we will be studying. Chapter 2 provides a derivation of the NEGF and its representation on the Keldysh contour. Chapter 3 derives the KBE from MBPT and introduces the self- energy and its conservation laws. Chapter 4 introduces the numerical techniques used to integrate the KBE, as well as an analysis of its causality structure. In Chapter 5, we use compression techniques to lower the complexity of the KBE, and provide details on its implementation. In Chapter 6 we apply the compression algorithm to study the Higgs mode in a superconducting system. In Chapter 7, we couple our solver to an adaptive integrator and analyze its performance against fixed-timestep integrators. Finally, Chapter 8 provides a summary of the results and an outlook for future work on efficient solvers.
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 nonequilibrium description of quantum systems requires, for more than two or three particles, the use of a reduced description to be numerically tractable. Two possible approaches are based on either reduced density matrices or nonequilibrium Green functions (NEGF). Both concepts are formulated in terms of hierarchies of coupled equations—the Bogoliubov-Born-Green-Kirkwood-Yvon (BBGKY) hierarchy for the reduced density operators and the Martin-Schwinger-hierarchy (MS) for the Green functions, respectively. In both cases, similar approximations are introduced to decouple the hierarchy, yet still many questions regarding the correspondence of both approaches remain open. Here we analyze this correspondence by studying the generalized Kadanoff–Baym ansatz (GKBA) that reduces the NEGF to a single-time theory. Starting from the BBGKY-hierarchy we present the approximations that are necessary to recover the GKBA result both, with Hartree-Fock propagators (HF-GKBA) and propagators in second Born approximation. To test the quality of the HF-GKBA, we study the dynamics of a 4-electron Hubbard nanocluster starting from a strong nonequilibrium initial state and compare to exact results and the Wang-Cassing approximation to the BBGKY hierarchy presented recently by Akbari et al. [1].
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.