Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userThe G1--G2 Scheme: Dramatic Acceleration of Nonequilibrium Green Functions Simulations Within the Hartree--Fock-GKBA
66
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic 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 [Schlunzen 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$.
rate research
Read More
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) which has substantially extended the scope of NEGF simulations. Here we present a novel approach to GKBA-NEGF simulations that is of order $T$, and demonstrate its remarkable capabilities.
We study the anisotropic quantum Heisenberg antiferromagnet for spin-1/2 that interpolates smoothly between the one-dimensional (1D) and the two-dimensional (2D) limits. Using the spin Hartree-Fock approach we construct a quantitative theory of heat capacity in the quasi-1D regime with a finite coupling between spin chains. This theory reproduces closely the exact result of Bethe Ansatz in the 1D limit and does not produces any spurious phase transitions for any anisotropy in the quasi-1D regime at finite temperatures in agreement with the Mermin-Wagner theorem. We study the static spin-spin correlation function in order to analyse the interplay of lattice geometry and anisotropy in these systems. We compare the square and triangular lattice. For the latter we find that there is a quantum transition point at an intermediate anisotropy of $sim0.6$. This quantum phase transition establishes that the quasi-1D regime extends upto a particular point in this geometry. For the square lattice the change from the 1D to 2D occurs smoothly as a function of anisotropy, i.e. it is of the crossover type. Comparing the newly developed theory to the available experimental data on the heat capacity of $rm{Cs}_2rm{CuBr}_4$ and $rm{Cs}_2rm{CuCl}_4$ we extract the microscopic constants of the exchange interaction that previously could only be measured using inelastic neutron scattering in high magnetic fields.
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 [Schluenzen 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.
This article presents an overview on recent progress in the theory of nonequilibrium Green functions (NEGF). NEGF, presently, are the only textit{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. 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 GWAx, the TMA 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 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, an easy use of these approximations in numerical applications is made possible.
The energy gap of correlated Hubbard clusters is well studied for one-dimensional systems using analytical methods and density-matrix-renormalization-group (DMRG) simulations. Beyond 1D, however, exact results are available only for small systems by quantum Monte Carlo. For this reason and, due to the problems of DMRG in simulating 2D and 3D systems, alternative methods such as Green functions combined with many-body approximations (GFMBA), that do not have this restriction, are highly important. However, it has remained open whether the approximate character of GFMBA simulations prevents the computation of the Hubbard gap. Here we present new GFMBA results that demonstrate that GFMBA simulations are capable of producing reliable data for the gap which agrees well with the DMRG benchmarks in 1D. An interesting observation is that the accuracy of the gap can be significantly increased when the simulations give up certain symmetry restriction of the exact system, such as spin symmetry and spatial homogeneity. This is seen as manifestation and generalization of the symmetry dilemma introduced by Lowdin for Hartree--Fock wave function calculations.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
