Abstract
The time evolution in quantum many-body systems after external excitations is attracting high interest in many fields. 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^2$-scaling 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 [Schl\"unzen \textit{et al.}, Phys. Rev. Lett. \textbf{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$.
Highlights
Nonequibrium Green functions (NEGF) [1,2,3] have proven highly successful in simulations of the dynamics of correlated many-body systems
If the HF-generalized Kadanoff-Baym ansatz (GKBA) is applied to improved self energies, such as the T -matrix self energy [19,27], which is required for strongly correlated systems [18], or the GKBA with SOA (GW) self energy [28] which is required to capture dynamical-screening effects, the CPU-time scaling is again increased to Nt3
In a recent Letter we reported a breakthrough for NEGF simulations within the HF-GKBA scheme: We demonstrated that time-linear scaling, i.e., a CPU time that is of order Nt1, can be achieved if the equations of motion are properly reformulated, without any approximations
Summary
Nonequibrium Green functions (NEGF) [1,2,3] have proven highly successful in simulations of the dynamics of correlated many-body systems This is due to a number of attractive properties that include conservation laws and the existence of systematic approximation schemes that are based on Feynman diagrams. If the HF-GKBA is applied to improved self energies, such as the T -matrix self energy [19,27], which is required for strongly correlated systems [18], or the GW self energy [28] which is required to capture dynamical-screening effects, the CPU-time scaling is again increased to Nt3. In a recent Letter we reported a breakthrough for NEGF simulations within the HF-GKBA scheme: We demonstrated that time-linear scaling, i.e., a CPU time that is of order Nt1, can be achieved if the equations of motion are properly reformulated, without any approximations.
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have