The extraordinary quantum properties of nonequilibrium systems governed by dissipative dynamics have become a focal point in contemporary scientific inquiry. The nonequilibrium Green's functions (NEGF) theory provides a versatile method for addressing driven nondissipative systems, utilizing the powerful diagrammatic technique to incorporate correlation effects. We here present a second-quantization approach to the dissipative NEGF theory, reformulating Keldysh ideas to accommodate Lindbladian dynamics and extending the Kadanoff-Baym equations accordingly. Generalizing diagrammatic perturbation theory for many-body Lindblad operators, the formalism enables correlated and dissipative real-time simulations for the exploration of transient and steady-state changes in the electronic, transport, and optical properties of materials.
The extraordinary quantum properties of nonequilibrium systems governed by dissipative dynamics have become a focal point in contemporary scientific inquiry. The nonequilibrium Green's functions (NEGF) theory provides a versatile method for addressing driven nondissipative systems, utilizing the powerful diagrammatic technique to incorporate correlation effects. We here present a second-quantization approach to the dissipative NEGF theory, reformulating Keldysh ideas to accommodate Lindbladian dynamics and extending the Kadanoff-Baym equations accordingly. Generalizing diagrammatic perturbation theory for many-body Lindblad operators, the formalism enables correlated and dissipative real-time simulations for the exploration of transient and steady-state changes in the electronic, transport, and optical properties of materials.
We present a unified many-body perturbation theory for open quantum systems that treats dissipation, correlations, and external driving on equal footing. Using a Keldysh-Lindblad formalism, we introduce a diagrammatic treatment of dissipative interaction lines representing quasiparticle flows and fluctuations. Two new Feynman rules render the evaluation of dissipative diagrams compact and systematically improvable, while preserving the Keldysh and anti-Hermitian symmetries of the closed-system theory. Consequently, the structure of the Kadanoff-Baym equations remains unchanged, enabling existing numerical methods to be directly applied. To illustrate this, we derive dissipative versions of the second Born and GW approximations and identify the physical content of the self-energy components. Moreover, we demonstrate that time-linear approximations to the full Kadanoff-Baym equations retain their closed-form structure and can be efficiently used to simulate relaxation and decoherence dynamics. The impact of dissipation-induced correlations is illustrated in the driven Haldane model, where quasiparticles exhibit nontrivial stabilization and acquire lifetimes that far exceed those of the bare system. This framework establishes a general route toward first-principles modeling of correlated, driven, and dissipative quantum materials.
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.
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.
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.
Abstract We review one of the most versatile theoretical approaches to the study of time-dependent correlated quantum transport in nano-systems: the non-equilibrium Green’s function (NEGF) formalism. Within this formalism, one can treat, on the same footing, inter-particle interactions, external drives and/or perturbations, and coupling to baths with a (piece-wise) continuum set of degrees of freedom. After a historical overview on the theory of transport in quantum systems, we present a modern introduction of the NEGF approach to quantum transport. We discuss the inclusion of inter-particle interactions using diagrammatic techniques, and the use of the so-called embedding and inbedding techniques which take the bath couplings into account non-perturbatively. In various limits, such as the non-interacting limit and the steady-state limit, we then show how the NEGF formalism elegantly reduces to well-known formulae in quantum transport as special cases. We then discuss non-equilibrium transport in general, for both particle and energy currents. Under the presence of a time-dependent drive—encompassing pump–probe scenarios as well as driven quantum systems—we discuss the transient as well as asymptotic behavior, and also how to use NEGF to infer information on the out-of-equilibrium system. As illustrative examples, we consider model systems general enough to pave the way to realistic systems. These examples encompass one- and two-dimensional electronic systems, systems with electron–phonon couplings, topological superconductors, and optically responsive molecular junctions where electron–photon couplings are relevant.
Within the Landauer approach to conduction discussed in the previous chapter, electron interactions have been included only at the mean-field level. This is quite a strong approximation, especially in nanojunctions, where large current densities are common. I therefore want to go beyond this level of description. I can follow two different routes. (1) I can develop a functional theory of quantum correlations via time-dependent effective single-particle equations, as I will do in Chapter 7 where I will use dynamical density-functional theories within the micro-canonical approach to conduction. (2) I can try to solve the time-dependent Schrödinger equation 1.16 – or its mixed-state version, the Liouville-von Neumann equation 1.60 – directly. Written this way, this last proposition seems hopeless. In reality, we can employ a many-body technique, known as the non-equilibrium Green's function formalism (NEGF), also referred to as the Keldysh formalism (Keldysh, 1964; Kadanoff and Baym, 1962), which allows us, at least in principle, to do just that: solve the time-dependent Schrödinger equation for an interacting many-body system exactly, from which one can, in principle, calculate the time-dependent current. This is done by solving equations of motion for specific time-dependent single-particle Green's functions, from which the physical properties of interest, such as the charge and current densities, can be obtained. I have stressed the term “time-dependent” several times, because, as I will show in a moment, the NEGF is “exact” only when one solves the time-dependent Schrödinger equation for a closed quantum system, subject to deterministic perturbations: the system is closed but not necessarily isolated. These perturbations may drive the system far away from its initial state of thermodynamic equilibrium (Keldysh, 1964).