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

The nonequilibrium dynamics of strongly-correlated fermions in lattice systems have attracted considerable interest in the condensed matter and ultracold atomic-gas communities. While experiments have made remarkable progress in recent years, there remains a need for the further development of theoretical tools that can account for both the nonequilibrium conditions and strong correlations. For instance, time-dependent theoretical quantum approaches based on the density matrix renormalization group (DMRG) methods have been primarily applied to one-dimensional setups. Recently, two-dimensional quantum simulations of the expansion of fermions based on nonequilibrium Green functions (NEGF) have been presented [Schlünzen et al., Phys. Rev. B 93, 035107 (2016)] that showed excellent agreement with the experiments. Here we present an extensive comparison of the NEGF approach to numerically accurate DMRG results. The results indicate that NEGF are a reliable theoretical tool for weak to intermediate coupling strengths in arbitrary dimensions and make long simulations possible. This is complementary to DMRG simulations which are particularly efficient at strong coupling.

Quantum dynamics in strongly correlated systems are of high current interest in many fields including dense plasmas, nuclear matter and condensed matter and ultracold atoms. An important model case are fermions in lattice systems that is well suited to analyze, in detail, a variety of electronic and magnetic properties of strongly correlated solids. Such systems have recently been reproduced with fermionic atoms in optical lattices which allow for a very accurate experimental analysis of the dynamics and of transport processes such as diffusion. The theoretical analysis of such systems far from equilibrium is very challenging since quantum and spin effects as well as correlations have to be treated non‐perturbatively. The only accurate method that has been successful so far are density matrix renormalization group (DMRG) simulations. However, these simulations are presently limited to one‐dimensional (1D) systems and short times. Extension of quantum dynamics simulations to two and three dimensions is commonly viewed as one of the major challenges in this field. Recently we have reported a breakthrough in this area [N. Schlünzen et al., Phys. Rev. B (2016)] where we were able to simulate the expansion dynamics of strongly correlated fermions in a Hubbard lattice following a quench of the confinement potential in 1D, 2D and 3D. The results not only exhibited excellent agreement with the experimental data but, in addition, revealed new features of the short‐time dynamics where correlations and entanglement are being build up. The method used in this work are nonequilibrium Green functions (NEGF) which are found to be very powerful in the treatment of fermionic lattice systems filling the gap presently left open by DMRG in 2D and 3D. In this paper we present a detailed introduction in the NEGF approach and its application to inhomogeneous Hubbard clusters. In detail we discuss the proper strong coupling approximation which is given by T ‐matrix selfenergies that sum up two‐particle scattering processes to infinite order. The efficient numerical implemen‐tation of the method is discussed in detail as it has allowed us to achieve dramatic performance gains. This has been the basis for the treatment of more than 100 particles over large time intervals. The numerical results presented in this paper concentrate on the diffusion in 1D to 3D lattices. We find that the expansion dynamics consist of three different phases that are linked with the build‐up of correlations. In the long time limit, a universal scaling with the particle number is revealed. By extrapolating the expansion velocities to the macroscopic limit, the obtained results show excellent agreement with recent experiments on ultracold fermions in optical lattices. Moreover we present results for the site‐resolved behavior of correlations and entanglement that can be directly compared with experiments using the recently developed atomic microscope technique. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)

The dynamics of strongly correlated fermions following an external excitation reveals extremely rich collective quantum effects. Examples are fermionic atoms in optical lattices, electrons in correlated materials, and dense quantum plasmas. Presently, the only quantum-dynamics approach that rigorously describes these processes in two and three dimensions is nonequilibrium Green functions (NEGF). However, NEGF simulations are computationally expensive due to their T 3 scaling with the simulation duration T . Recently, T 2 scaling was achieved with the generalized Kadanoff-Baym ansatz (GKBA), for second order Born (SOA) selfenergies, which has substantially extended the scope of NEGF simulations. Here we demonstrate that GKBA-NEGF simulations can be performed with order T 1 scaling, both for SOA and GW selfenergies, and point out the remarkable capabilities of this approach.

No abstract

A combined nonequilibrium Green functions-Ehrenfest dynamics approach is developed that allows for a time-dependent study of the energy loss of a charged particle penetrating a strongly correlated system at zero and finite temperature. Numerical results are presented for finite inhomogeneous two-dimensional Fermi-Hubbard models, where the many-electron dynamics in the target are treated fully quantum mechanically and the motion of the projectile is treated classically. The simulations are based on the solution of the two-time Dyson (Keldysh-Kadanoff-Baym) equations using the second-order Born, third-order and T-matrix approximations of the self-energy. As application, we consider protons and helium nuclei with a kinetic energy between 1 and 500 keV/u passing through planar fragments of the two-dimensional honeycomb lattice and, in particular, examine the influence of electron-electron correlations on the energy exchange between projectile and electron system. We investigate the time dependence of the projectile's kinetic energy (stopping power), the electron density, the double occupancy and the photoemission spectrum. Finally, we show that, for a suitable choice of the Hubbard model parameters, the results for the stopping power are in fair agreement with ab-initio simulations for particle irradiation of single-layer graphene.Comment: 13 pages, 12 figure

1] V.E. Fortov, R.I. Ilkaev, V.A. Arinin, V.V. Burtzev, V.A. Golubev, I.L. Iosilevskiy, V.V. Khrustalev, A.L. Mikhailov, M.A. Mochalov, V.Ya. Ternovoi, and M.V. Zhernokletov, Phase Transition in a Strongly Non-ideal Deuterium Plasma Generated by Quasi-Isentropical Compression at Megabar Pressures, Phys. Rev. Lett. 99, 185001 (2007). [2] T.

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 T 3 -scaling with the simulation duration T . Recently, T 2scaling was achieved with the generalized Kadanoff-Baym ansatz (GKBA), for the second-order Born (SOA) selfenergy, which has substantially extended the scope of NEGF simulations. In a recent Letter [Schlnzen 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 G1-G2 scheme. This allows one to perform the same simulations with order T 1 -scaling, both for SOA and GW selfenergies 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 selfenergies, including the T matrix in the particle-particle and particle-hole channels, and the dynamically screened-ladder approximation. Finally, we demonstrate that, for all selfenergies, the CPU time scaling of the G1-G2 scheme with the basis dimension, N b , can be improved compared to our first report: the overhead compared to the original GKBA, is not more than an additional factor N b , even for Hubbard systems.real-time components [28], that we define as follows(2),F ijkl (z, z , z, z ) = G il (z, z )G jk (z , z) , G corr ijkl (z, z ) := G(2),corr ijkl (z, z, z , z + ) , G H,≷