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.
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.
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 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.
In this work we discuss the application of nonequilibrium Green functions theory to atomic and molecular systems with the aim to study charge and energy transport in these systems. We apply the Kadanoff-Baym equations to atoms and diatomic molecules initially in the ground state. The results obtained for the correlated initial states are used to analyze variational energy functionals of the Green function which are shown to perform very well. We further show an application of the Kadanoff-Baym equations to a molecule exposed to an external laser field. Finally we discuss the connection between nonequilibrium Green function theory and time-dependent density-functional theory with the aim to develop better density functionals in order to treat larger systems than those attainable with the nonequilibrium Green function method.
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.
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.