Abstract Correlated classical and quantum many‐particle systems out of equilibrium are of high interest in many fields, including dense plasmas, correlated solids, and ultracold atoms. Accurate theoretical description of these systems is challenging both, conceptionally and with respect to computational resources. While for classical systems, in principle, exact simulations are possible via molecular dynamics, this is not the case for quantum systems. Alternatively, one can use many‐particle approaches such as hydrodynamics, kinetic theory, or nonequilibrium Green functions (NEGF). However, NEGF exhibit a very unfavorable cubic scaling of the CPU time with the number of time steps. An alternative is the G1–G2 scheme [N. Schlünzen et al., Phys. Rev. Lett. 124, 076601 (2020)] which allows for NEGF simulations with time linear scaling, however, at the cost of large memory consumption. The reason is the need to store the two‐particle correlation function. This problem can be overcome for a number of approximations by reformulating the kinetic equations in terms of fluctuations – an approach that was developed, for classical systems, by Yu.L. Klimontovich [JETP 33, 982 (1957)]. Here, we present an overview of his ideas and extend them to quantum systems. In particular, we demonstrate that this quantum fluctuations approach can reproduce the nonequilibrium GW approximation [E. Schroedter et al., Cond. Matt. Phys. 25, 23401 (2022)] promising high accuracy at low computational cost which arises from an effective semiclassical stochastic sampling procedure. We also demonstrate how to extend the approach to the two‐time exchange‐correlation functions and the density response properties. [E. Schroedter et al., Phys. Rev. B 108, 205109 (2023)]. The results are equivalent to the Bethe–Salpeter equation for the two‐time exchange‐correlation function when the generalized Kadanoff‐Baym ansatz with Hartree‐Fock propagators is applied [E. Schroedter and M. Bonitz, phys. stat. sol. (b) 2024, 2300564].
Abstract Correlated classical and quantum many‐particle systems out of equilibrium are of high interest in many fields, including dense plasmas, correlated solids, and ultracold atoms. Accurate theoretical description of these systems is challenging both, conceptionally and with respect to computational resources. While for classical systems, in principle, exact simulations are possible via molecular dynamics, this is not the case for quantum systems. Alternatively, one can use many‐particle approaches such as hydrodynamics, kinetic theory, or nonequilibrium Green functions (NEGF). However, NEGF exhibit a very unfavorable cubic scaling of the CPU time with the number of time steps. An alternative is the G1–G2 scheme [N. Schlünzen et al., Phys. Rev. Lett. 124, 076601 (2020)] which allows for NEGF simulations with time linear scaling, however, at the cost of large memory consumption. The reason is the need to store the two‐particle correlation function. This problem can be overcome for a number of approximations by reformulating the kinetic equations in terms of fluctuations – an approach that was developed, for classical systems, by Yu.L. Klimontovich [JETP 33, 982 (1957)]. Here, we present an overview of his ideas and extend them to quantum systems. In particular, we demonstrate that this quantum fluctuations approach can reproduce the nonequilibrium GW approximation [E. Schroedter et al., Cond. Matt. Phys. 25, 23401 (2022)] promising high accuracy at low computational cost which arises from an effective semiclassical stochastic sampling procedure. We also demonstrate how to extend the approach to the two‐time exchange‐correlation functions and the density response properties. [E. Schroedter et al., Phys. Rev. B 108, 205109 (2023)]. The results are equivalent to the Bethe–Salpeter equation for the two‐time exchange‐correlation function when the generalized Kadanoff‐Baym ansatz with Hartree‐Fock propagators is applied [E. Schroedter and M. Bonitz, phys. stat. sol. (b) 2024, 2300564].
Correlated quantum many‐particle systems out of equilibrium are of high interest in many fields, including correlated solids, ultracold atoms, or dense plasmas. Accurate theoretical description of these systems is challenging both, conceptionally and with respect to computational resources. A quantum fluctuations approach is recently presented, which is equivalent to the nonequilibrium GW approximation that promises high accuracy at low computational cost. The method exhibits process time scaling that is linear in the number of time steps, like the G1–G2 scheme, however, at a much reduced computer memory cost. In a second publication, this approach is extended to the two‐time exchange–correlation functions and the dynamic density response properties. Herein, the properties of this approach are analyzed in more detail. The physical meaning of the central approximation, the quantum polarization approximation, is established. It is demonstrated that the method is equivalent to the Bethe–Salpeter equation for the two‐time exchange–correlation function when the generalized Kadanoff–Baym ansatz with Hartree–Fock propagators is applied.
The quantum dynamics of fermionic or bosonic many-body systems following external excitation can be successfully studied using two-time nonequilibrium Green's functions (NEGF) or single-time reduced density matrix methods. Approximations are introduced via a proper choice of the many-particle self-energy or decoupling of the BBGKY hierarchy. These approximations are based on Feynman's diagram approaches or on cluster expansions into single-particle and correlation operators. Here, we develop a different approach where, instead of equations of motion for the many-particle NEGF (or density operators), single-time equations for the correlation functions of fluctuations are analyzed. We present a derivation of the first two equations of the alternative hierarchy of fluctuations and discuss possible decoupling approximations. In particular, we derive the polarization approximation (PA) which is shown to be equivalent to the single-time version [following by applying the generalized Kadanoff-Baym ansatz (GKBA)] of the nonequilibrium GW approximation with exchange effects of NEGF theory, for weak coupling. The main advantage of the quantum fluctuations approach is that the standard ensemble average can be replaced by a semiclassical average over different initial realizations, as was demonstrated before by Lacroix and co-workers [see e.g. D. Lacroix et al., Phys. Rev. B, 2014, 90, 125112]. Here, we introduce the stochastic GW (SGW) approximation and the stochastic polarization approximation (SPA) which are demonstrated to be equivalent to the single-time GW approximation without and with exchange, respectively, in the weak coupling limit. In addition to the standard stochastic approach to sample initial configurations, we also present a deterministic approach. Our numerical tests confirm that our approach has the same favorable linear scaling with the computation time as the recently developed G1-G2 scheme [Schluenzen et al., Phys. Rev. Lett., 2020, 124, 076601]. At the same time, the SPA and SGW approaches scale more favorably with the system size than the G1-G2 scheme, allowing to extend nonequilibrium GW calculations to bigger systems.
The theory of nonequilibrium Green functions (NEGF) has seen a rapid development over the recent three decades. Applications include diverse correlated many‐body systems in and out of equilibrium. Very good agreement with experiments and available exact theoretical results could be demonstrated if the proper selfenergy approximations were used. However, full two‐time NEGF simulations are computationally costly, as they suffer from a cubic scaling of the computation time with the simulation duration. Recently the G1–G2 scheme that exactly reformulates the generalized Kadanoff–Baym ansatz with Hartree–Fock propagators (HF‐GKBA) into time‐local equations is introduced, which achieves time‐linear scaling and allows for a dramatic speedup and extension of the simulations (Schluenzen et al. Phys. Rev. Lett. 2020, 124, 076601). Remarkably, this scaling is achieved quickly, and also for high‐level selfenergies, including the nonequilibrium GW and T‐matrix approximations (Joost et al. Phys. Rev. B 2020, 101, 245101). Even the dynamically screened ladder approximation is now feasible (Joost et al. Phys. Rev. B 2022, 105, 165155), and also applications to electron‐boson systems are demonstrated. Herein, an overview on recent results that are achieved with the G1–G2 scheme is presented. Problems and open questions are discussed and further ideas of how to overcome the current limitations of the scheme and present are presented. The G1–G2 scheme is illustrated by presenting applying it to the excitation dynamics of Hubbard clusters, to optical excitation of graphene, and to charge transfer during stopping of ions by correlated materials.
The quantum dynamics of correlated fermionic or bosonic many-body systems following external excitation can be successfully studied using nonequilibrium Green functions (NEGFs) or reduced density matrix methods. Approximations are introduced via a proper choice of the many-particle self-energy or decoupling of the BBGKY-hierarchy, respectively. These approximations are based on Feynman's diagram approaches or on cluster expansions into single-particle and correlation operators. In a recent paper [E. Schroedter, J.-P. Joost, and M. Bonitz, Condens. Matter Phys. 25, 23401 (2022)] we presented a different approach in which, instead of equations of motion for the many-particle NEGF (or density operators), equations for the correlation functions of fluctuations are analyzed. In particular, we derived the stochastic GW and polarization approximations that are closely related to the nonequilibrium GW approximation. Here, we extend this approach to the computation of two-time observables depending on the specific ordering of the underlying operators. In particular, we apply this extension to the calculation of the density correlation function and dynamic structure factor of correlated Hubbard clusters in and out of equilibrium.
In a previous publication [H. S. K\"ohler, Phys. Rev. C 51, 3232 (1995)] equilibration rates were studied for an extended system of nonequilibrium nuclear matter. The two-time Green's functions in the Kadanoff-Baym (KB) formalism were solved numerically. These quantum-mechanical calculations were compared with the Markovian (semi)classical Uehling-Uhlenbeck collision term often used in the analysis of collisions between heavy ions. The equilibration was found to be substantially slowed down by quantum effects. The theoretical as well as experimental study of nonequilibrium systems is of course of great interest in many areas of physics, e.g., within the realm of quantum transport in solid state devices. On the theoretical side various methods and approximations are the subject of intensive study. By the KB ansatz, for example, the equations can be simplified to a one-time kinetic equation while preserving the non-Markovian character. The generalized KB (GKB) ansatz of Lipavsky, Spicka, and Velicky [Phys. Rev. B 34, 6933 (1986)] is often preferred. In this paper approximations based upon the GKB ansatz are discussed in relation to the exact KB formalism. Numerical comparisons are made for both the strongly interacting nuclear medium as well as an electron plasma. A non-Markovian collision integral, including memory effects in a quasiparticle approximation, is often discussed in the literature. It is found that this approximation appreciably overestimates the collision rates even more than the classical (Boltzmann) approximation. The inclusion of a width in the spectral function, i.e., going beyond the quasiparticle approximation, gives considerable improvement. The competition between the memory and correlation effects is found to lead to an anomalous saturation of the relaxtion rate. \textcopyright{} 1996 The American Physical Society.
The method of nonequilibrium Green functions (NEGFs) constitutes a solid framework to describe closed and open quantum systems. It captures situations far from equilibrium including arbitrary perturbations and correlations - effects higher than first order in the binary interactions. On the contrary, in the absence of external fields, it reduces to the imaginary-time Matsubara formalism. The abilities (and successes) of the method originate from its self-consistency and the conserving character of many-body approximations that are introduced by diagram technique. On this account, NEGFs have, now for a long time, brought insight into the dynamics of nuclear matter, non-ideal quantum gases and (semiconductor) plasmas. However, applications to localized, finite and strongly interacting systems have emerged only within the recent decade. The dynamical properties follow from quantum kinetic equations - the two-time Kadanoff-Baym equations for the one-particle NEGF - the numerical solution of which is the central topic of this thesis. The non-Markovian structure of these equations thereby inhibits an unlimited, gradual propagation of the nonequilibrium Green function: Any self-consistent, future dynamics of the system is tightly linked to the (complete) past by non-local self-energies. Further, the success in obtaining solutions is being diminished if the system under investigation is not isotropic or inhomogeneous and, in addition, the spatial topology (e.g., the charge carrier or spin density) needs to be resolved in time. On top of this, small isolated systems are generally not candidates showing fast decoherence which offers little hope for smooth and well-damped integration kernels. Despite these difficulties, the thesis at hand is devoted to the computation of the two-time NEGF for finite, closed systems of which ensembles of electrons bound by nuclei forming atoms or molecules are the most natural representatives. Here, in the spirit of the well-studied one-dimensional helium atom, the electron dynamics in model atoms and molecules is addressed in ab initio fashion. To render calculations possible in the first place, (i) a novel hybrid representation of the NEGF is introduced based on the conception of finite elements (FE) and the discrete variable representation (DVR), and (ii) a parallel algorithm is developed that allows for an efficient, well scalable distributed-memory computation of the full two-time Green function. The first point attributes plenty of flexibility and makes the method readily adjustable to different Hamiltonians avoiding the need to, by hand or in intricate numerical manner, calculate the whole bunch of matrix elements of one- and two-particle operators. Also, in consequence, it drastically simplifies the evaluation of self-energies. The second point makes the NEGF approach ready for large-scale, high-performance computing facilities, the resources (of computer power: mainly memory!) of which are indispensable for (some of) the calculations presented. The applicability of the developed methods is demonstrated inter alia on the basis of helium, beryllium, hydrogen and lithium hydride modeled in one spatial dimension whereby electron-electron correlations are treated in the second(-order) Born approximation. After a thorough discussion on the self-consistent ground states, the thesis focuses on the electron dynamics in the linear response regime as well as in the presence of strong external fields. Overall, to assess the the performance of the second Born approximation, the results are compared to the Hartree-Fock approximation, on the one hand, and to the exact solution of the time-(in)dependent Schrödinger equation, on the other hand. In addition, aiming at the possibility of extended time-dependent calculations, the behavior of the generalized Kadanoff-Baym ansatz is discussed.