A time-stepping scheme with adaptivity in both the step size and the integration order is presented in the context of non-equilibrium dynamics described via Kadanoff-Baym equations. The accuracy and effectiveness of the algorithm are analysed by obtaining numerical solutions of exactly solvable models. We find a significant reduction in the number of time-steps compared to fixed-step methods. Due to the at least quadratic scaling of Kadanoff-Baym equations, reducing the amount of steps can dramatically increase the accessible integration time, opening the door for the study of long-time dynamics in interacting systems. A selection of illustrative examples is provided, among them interacting and open quantum systems as well as classical stochastic processes. An open-source implementation of our algorithm in the scientific-computing language Julia is made available.
A time-stepping scheme with adaptivity in both the step size and the integration order is presented in the context of non-equilibrium dynamics described via Kadanoff-Baym equations. The accuracy and effectiveness of the algorithm are analysed by obtaining numerical solutions of exactly solvable models. We find a significant reduction in the number of time-steps compared to fixed-step methods. Due to the at least quadratic scaling of Kadanoff-Baym equations, reducing the amount of steps can dramatically increase the accessible integration time, opening the door for the study of long-time dynamics in interacting systems. A selection of illustrative examples is provided, among them interacting and open quantum systems as well as classical stochastic processes. An open-source implementation of our algorithm in the scientific-computing language Julia is made available.
The nonequilibrium Green's function gives access to one-body observables for quantum systems. Of particular interest are quantities such as density, currents, and absorption spectra which are important for interpreting experimental results in quantum transport and spectroscopy. We present an integration scheme for the Green's function's equations of motion, the Kadanoff-Baym equations (KBE), which is both adaptive in the time integrator step size and method order as well as the history integration order. We analyze the importance of solving the KBE self-consistently and show that adapting the order of history integral evaluation is important for obtaining accurate results. To examine the efficiency of our method, we compare runtimes to a state-of-the-art fixed time step integrator for several test systems and show an order of magnitude speedup at similar levels of accuracy. Published by the American Physical Society 2024
In this thesis we present methodologies that have been developed to efficiently solve for the non-equilibrium Green’s Function (NEGF) of interacting electron systems, and demonstrate their applications to several different condensed matter systems. These solvers are based on the numerical integration of the equations of motion for the NEGF, which are known as the Kadanoff-Baym Equations (KBE). The KBE are a general set of equations which allow for the study of a wide range of systems that are of interest to condensed matter physicists, including superconductivity, nuclear matter, molecules and atoms, and qubits, among many others. They are also useful for solving impurity models and systems in the dynamical mean-field theory (DMFT) approximation through their ability to incorporate dynamical hybridization functions. The generality of these equations comes from their connection to time-dependent Many-Body Perturbation Theory (MBPT), which introduces non-Markovian dynamical correlations through the self-energy functional. MBPT ensures important conservation laws, such as those for density and energy, are obeyed no matter the type of time-dependent perturbations that are applied to the Hamiltonian. The non-Markovian correlations manifest in a kernel that couples the dynamics to the entire history of the system. This coupling leads to computational costs that become prohibitively expensive as longer-time dynamics are probed, this is especially true for systems which decohere slowly, leading to undamped dynamics and kernels. Due to the large amount of interest in non-equilibrium quantum dynamics, the availability of efficient solvers for the KBE are of great importance. This thesis will focus on the implementation and analysis of integration methods for the KBE, including compression algorithms which lower computational complexity and adaptive integration techniques which allow for resolution of dynamics which vary over wide ranges of timescales. Chapter 1 gives an introduction to quantum systems and the problems we will be studying. Chapter 2 provides a derivation of the NEGF and its representation on the Keldysh contour. Chapter 3 derives the KBE from MBPT and introduces the self- energy and its conservation laws. Chapter 4 introduces the numerical techniques used to integrate the KBE, as well as an analysis of its causality structure. In Chapter 5, we use compression techniques to lower the complexity of the KBE, and provide details on its implementation. In Chapter 6 we apply the compression algorithm to study the Higgs mode in a superconducting system. In Chapter 7, we couple our solver to an adaptive integrator and analyze its performance against fixed-timestep integrators. Finally, Chapter 8 provides a summary of the results and an outlook for future work on efficient solvers.
Although the Kadanoff-Baym equations are typically solved using time-stepping methods, iterative global-in-time solvers offer potential algorithmic advantages, particularly when combined with compressed representations of two-time objects. We examine the computational complexity and stability of several global-in-time iterative methods, including multiple variants of fixed point iteration, Jacobian-free methods, and a Newton-Krylov method using automatic differentiation. We consider the ramped and periodically-driven Falicov-Kimball and Hubbard models within time-dependent dynamical mean-field theory. Although we observe that several iterative methods yield stable convergence at large propagation times, a standard forward fixed point iteration does not. We find that the number of iterations required to converge to a given accuracy with a fixed time step size scales roughly linearly with the number of time steps. This scaling is associated with the formation of a propagating front in the residual error, whose velocity is method-dependent. We identify key challenges which must be addressed in order to make global solvers competitive with time-stepping methods.
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.
(Classical) Boltzmann equations suffer from several shortcomings as compared to their quantum mechanical generalizations, the so-called Kadanoff-Baym equations. Nevertheless, in practical calculations approximations to Boltzmann equations are widely used. A prominent example is the computation of the baryon asymmetry of the universe. In this work, we present a detailed comparison between the Kadanoff-Baym and Boltzmann equations in the framework of a scalar phi 4 quantum field theory. We solve the Kadanoff-Baym and the Boltzmann equations numerically in 3 + 1 dimensions and compare their predictions on the evolution of systems out of thermal equilibrium for various initial conditions. The obtained numerical solutions reveal significant discrepancies between both types of equations. Most notably, the universality respected by Kadanoff-Baym equations is severely restricted in the case of Boltzmann equations.