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.

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.

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,≷

Abstract. The electrostatic potential of a moving dust grain in a complex plasma with magnetized ions is computed using linear response theory, thereby extending our previous work for unmagnetized plasmas [P. Ludwig et al., New J. Phys. 14, 053016 (2012)]. In addition to the magnetic field, our approach accounts for a finite ion temperature as well as ion-neutral collisions. Our recently introduced code Kielstream is used for an efficient calculation of the dust potential. Increasing the magnetization of the ions, we find that the shape of the potential crucially depends on the Mach number M . In the regime of subsonic ion flow (M < 1), a strong magnetization gives rise to a potential distribution that is qualitatively different from the unmagnetized limit, while for M > 1 the magnetic field effectively suppresses the plasma wakefield.