This dissertation investigates the possibilities and limitations of time-dependent many-body perturbation theory by studying small Hubbard clusters for which the exact solution is available. The first part of the thesis is comprised of a short introduction to the concepts and methodologies used. The second part consists of an review of the main findings of the thesis and a short summary of four original papers. In paper I we study the dynamics of short Hubbard chains within many-body perturbation theory and compare to the exact solution. The main outcomes are that the T-matrix approximation works well in the low filling regime and that all approximation which involve correlation effects develop an unphysical steady state. In paper II we present the methodology used in paper I. We find that the correlation-induced damping appears even in the presence of leads and that there are multiple steady and quasi steady states. In paper III we make a review of the status of time dependent density functional theory for lattice models. In particular we find that a non-perturbative adiabatic local density approximation describes strong correlations well while many-body perturbation theory accounts better for non-adiabatic effects. In paper IV we propose a method to obtain the double occupancy from the Kadanoff-Baym equations. We show that the positiveness condition may be violated in the GW or the second Born approximation but fulfiled in the T-matrix approximation. We apply this method to obtain the local entanglement entropy.
This dissertation investigates the possibilities and limitations of time-dependent many-body perturbation theory by studying small Hubbard clusters for which the exact solution is available. The first part of the thesis is comprised of a short introduction to the concepts and methodologies used. The second part consists of an review of the main findings of the thesis and a short summary of four original papers. In paper I we study the dynamics of short Hubbard chains within many-body perturbation theory and compare to the exact solution. The main outcomes are that the T-matrix approximation works well in the low filling regime and that all approximation which involve correlation effects develop an unphysical steady state. In paper II we present the methodology used in paper I. We find that the correlation-induced damping appears even in the presence of leads and that there are multiple steady and quasi steady states. In paper III we make a review of the status of time dependent density functional theory for lattice models. In particular we find that a non-perturbative adiabatic local density approximation describes strong correlations well while many-body perturbation theory accounts better for non-adiabatic effects. In paper IV we propose a method to obtain the double occupancy from the Kadanoff-Baym equations. We show that the positiveness condition may be violated in the GW or the second Born approximation but fulfiled in the T-matrix approximation. We apply this method to obtain the local entanglement entropy.
We present in detail a method we recently introduced (PRL. 103, 176404 (2009)) to describe finite systems in and out of equilibrium, where the evolution in time is performed via the Kadanoff-Baym Equations (KBE) within Many-Body Perturbation Theory (MBPT). The main property we analyze is the time-dependent density. We also study is the exchange-correlation potential of TDDFT, obtained via reverse engineering from the time-dependent density. 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 KBE 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 overall superior at all electron densities. Such comparisons permit a general assessment of the whole idea of applying MBPT, in the KBE sense, to finite systems. A striking outcome of our analysis is that when the system evolves under a strong external field, the KBE 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. To illustrate this point we present selected results for a system coupled to contacts within the T-matrix and second Born approximation. The extensive characterization we performed indicates that this behavior is present whenever approximate self energies, based upon infinite partial summations, are used. A second important result is that, for isolated clusters, the steady state reached is not unique but depends on how one switches on the external field. This may also true for clusters connected to leads.
We perform a comparative study of exact and approximate time-evolved densities in small Hubbard chains. The approximate densities are obtained via many-body perturbation theory (Hartree-Fock, 2nd Born, GW and T-matrix approximations) within the framework of the time-dependent Kadanoff-Baym equations. Benchmarking approximate results against exact ones allows us to address two rather fundamental issues in the non equilibrium dynamics of strongly correlated systems. I) A characterisation of the performance of several standard MBAs in the non-equilibrium regime. Having a definite notion of how good a specific MBA can be is highly relevant to its application to cases (typically, infinite systems) where exact solutions are not available. Our results show that the T-matrix approximation is overall superior to the other MBAs, at all electron densities. II) A scrutiny of the whole idea of Many Body Perturbation Theory in the Kadanoff-Baym sense, when applied to finite systems. The surprising outcome of our study is that during the time evolution, the KBE develop an unphysical steady state solution. This is a genuinely novel feature of the time-dependent KBE, i.e. is not inherited from possible limitations/approximations in the calculation of the initial state. Our extensive numerical characterisation gives robust evidence that the problem occurs in general, whenever MBPT is applied to finite systems, and approximate self energies based upon infinite partial summations are used. We also offer some more conceptual and general consideration on the dependence of this behaviour on the number of particles and system size. This is followed by our conclusions and glimpses of future work.
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.
In this thesis, an approximate method introduced to solve time-dependent many-body problems known as time-dependent many-body perturbation theory is studied. Many-body perturbation theory for interacting electrons and phonons is reviewed. In particular, the electron propagator G and an unconventional two-component phonon propagator, which satisfy coupled integral Dyson equations, are introduced. In practice, the associated integral kernels known as the electron Σ and phonon self-energies need to be approximated. The conserving approximations known as the Hartree (-Fock) and the first and second Born approximations, which respect the continuity equation between the electron density and current, are considered in this work. Time-dependent systems of interest are studied in this thesis by using the integro-differential forms of the Dyson equations referred to as the Kadanoff-Baym Equations (KBE). The Kadanoff-Baym equations are introduced for the electron and phonon propagators unconventionally as coupled firstorder integro-differential equations. It is reviewed how these equations are solved numerically by describing an integration rule, a class of practical methods and a parallel implementation of the numerical method. In addition, documentation of how the Kadanoff-Baym equations allow to solve the Bethe-Salpeter Equation (BSE) with the kernel δΣ/δG for the density response function, is provided. In two of the enclosed publications, we benchmarked observables obtained by using the Hartree and partially and fully self-consistent Born approximations against numerically exact results for the two- site, two-electron Holstein model. In this model, the two electrons which are constrained to move between two lattice sites interact indirectly via the electron-phonon coupling. It was observed that only the fully self-consistent Born approximation could cope qualitatively correctly with the competition between the delocalizing and localizing effects of the kinetic and interaction energies. For the other two approximations, spurious bifurcative symmetry breaking with an associated unbounded density response was observed. Nevertheless, also the self-consistent Born approximation was concluded to fail in describing the bipolaronic behavior of the true system. In the third publication, we benchmarked the Hartree-Fock and second Born approximations against an exact method for the few-site Hubbard and Pariser-Parr-Pople models in which the underlying lattice is inert and the electrons interact amongst themselves directly. It was found that the second Born approximation is capable of describing the so-called correlation induced doubly-excited states. This is not possible for time-local approximations such as Hartree-Fock. In addition to the qualitative results, which highlight successes of the applied simple self-energy approximations, the approximate and exact results were also compared on a more quantitative level. It is the quantitative and qualitative results combined which are used in this thesis to assess the quality of the many-body approximations used, with the aim to better understand when these methods are predictive.
We study the nonequilibrium dynamics of small, strongly correlated clusters, described by a Hubbard Hamiltonian, by propagating in time the Kadanoff-Baym equations within the Hartree-Fock, second Born, GW, and T-matrix approximations. We compare the results to exact numerical solutions. We find that the time-dependent T matrix is overall superior to the other approximations, and is in good agreement with the exact results in the low-density regime. In the long time limit, the many-body approximations attain an unphysical steady state which we attribute to the implicit inclusion of infinite-order diagrams in a few-body system.
In this thesis the electron transport through interacting quantum wires is studied.The time-dependent properties of the system are investigated by using the non-equilibrium Green function technique.This technique provides a way to include ultra-fast timedependent processes, arbitrary strong external fields and also brings forward an effective treatment of electron-electron and electron-phonon interactions.The electronelectron interactions are included via various self-energy approximations at Hartree-Fock, second-order Born and GW levels.These approximations are called conserving since they fulfill the physical conservation laws like the charge conservation law.The interacting central system is coupled to one-dimensional tight-binding leads, which are considered effectively as noninteracting.The system is perturbed by a bias voltage after which the time-evolution of the system is followed by Kadanoff-Baym equations.The non-equilibrium steady-state value of the current depends on the used self-energy approximation, which indicates the importance of the history-dependence in the calculation of current.This is an important observation in the aspect of time-dependent density functional theory, since it is one indication for the need of functionals with memory.The steady-state current is independent of the initial state, which is in accordance with the memory-loss theorem.The transient current shows small oscillations, which are consequence of the density oscillations especially in the first and last site.The density oscillations in the sites are effect of the electron transition from the lead to the central system and vice versa.