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 non-perturbative 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 selfconsistent 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 compares 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

No abstract

This article presents an overview on recent progress in the theory of nonequilibrium Green functions (NEGF). We discuss applications of NEGF simulations to describe the femtosecond dynamics of various finite fermionic systems following an excitation out of equilibrium. This includes the expansion dynamics of ultracold atoms in optical lattices following a confinement quench and the excitation of strongly correlated electrons in a solid by the impact of a charged particle. NEGF, presently, are the only ab initio quantum approach that is able to study the dynamics of correlations for long times in two and three dimensions. However, until recently, NEGF simulations have mostly been performed with rather simple selfenergy approximations such as the second-order Born approximation (SOA). While they correctly capture the qualitative trends of the relaxation towards equilibrium, the reliability and accuracy of these NEGF simulations has remained open, for a long time. Here we report on recent tests of NEGF simulations for finite lattice systems against exact-diagonalization and density-matrix-renormalization-group benchmark data. The results confirm the high accuracy and predictive capability of NEGF simulations—provided selfenergies are used that go beyond the SOA and adequately include strong correlation and dynamical-screening effects. With an extended arsenal of selfenergies that can be used effectively, the NEGF approach has the potential of becoming a powerful simulation tool with broad areas of new applications including strongly correlated solids and ultracold atoms. The present review aims at making such applications possible. To this end we present a selfcontained introduction to the theory of NEGF and give an overview on recent numerical applications to compute the ultrafast relaxation dynamics of correlated fermions. In the second part we give a detailed introduction to selfenergies beyond the SOA. Important examples are the third-order approximation, the approximation, the T-matrix approximation and the fluctuating-exchange approximation. We give a comprehensive summary of the explicit selfenergy expressions for a variety of systems of practical relevance, starting from the most general expressions (general basis) and the Feynman diagrams, and including also the important cases of diagonal basis sets, the Hubbard model and the differences occuring for bosons and fermions. With these details, and information on the computational effort and scaling with the basis size and propagation duration, readers will be able to choose the proper basis set and straightforwardly implement and apply advanced selfenergy approximations to a broad class of systems.

In non-equilibrium Green's function calculations the use of the Generalized Kadanoff-Baym Ansatz (GKBA) allows for a simple approximate reconstruction of the two-time Green's function from its time-diagonal value. With this a drastic reduction of the computational needs is achieved in time-dependent calculations, making longer time propagation possible and more complex systems accessible. This paper gives credit to the GKBA that was introduced 25 years ago. After a detailed derivation of the GKBA, we recall its application to homogeneous systems and show how to extend it to strongly correlated, inhomogeneous systems. As a proof of concept, we present results for a 2electron quantum well, where the correct treatment of the correlated electron dynamics is crucial for the correct description of the equilibrium and dynamic properties.

Finite lattice models are a prototype for strongly correlated quantum systems and capture essential properties of condensed matter systems. With the dramatic progress in ultracold atoms in optical lattices, finite fermionic Hubbard systems have become directly accessible in experiments, including their ultrafast dynamics far from equilibrium. Here, we present a theoretical approach that is able to treat these dynamics in any dimension and fully includes inhomogeneity effects. The method consists in stochastic sampling of mean-field trajectories and is found to be more accurate and efficient than current nonequilibrium Green functions approaches. This is demonstrated for Hubbard clusters with up to 512 particles in one, two and three dimensions.