With the technical progress of radio-frequency setups, high frequency quantum transport experiments have moved from theory to the lab. So far the standard theoretical approach used to treat such problems numerically - known as Keldysh or NEGF (Non Equilibrium Green's Functions) formalism - has not been very successful mainly because of a prohibitive computational cost. We propose a reformulation of the non-equilibrium Green's function technique in terms of the electronic wave functions of the system in an energy-time representation. The numerical algorithm we obtain scales now linearly with the simulated time and the volume of the system, and makes simulation of systems with 105-106atoms/sites feasible. We illustrate our method with the propagation and spreading of a charge pulse in the quantum Hall regime. We identify a classical and a quantum regime for the spreading, depending on the number of particles contained in the pulse. This numerical experiment is the condensed matter analogue to the spreading of a Gaussian wavepacket discussed in quantum mechanics textbooks.
With the technical progress of radio-frequency setups, high frequency quantum transport experiments have moved from theory to the lab. So far the standard theoretical approach used to treat such problems numerically - known as Keldysh or NEGF (Non Equilibrium Green's Functions) formalism - has not been very successful mainly because of a prohibitive computational cost. We propose a reformulation of the non-equilibrium Green's function technique in terms of the electronic wave functions of the system in an energy-time representation. The numerical algorithm we obtain scales now linearly with the simulated time and the volume of the system, and makes simulation of systems with 105-106atoms/sites feasible. We illustrate our method with the propagation and spreading of a charge pulse in the quantum Hall regime. We identify a classical and a quantum regime for the spreading, depending on the number of particles contained in the pulse. This numerical experiment is the condensed matter analogue to the spreading of a Gaussian wavepacket discussed in quantum mechanics textbooks.
As experimental probes have matured to observe ultrafast transient and high frequency responses of materials and devices, so to have the theoretical methods to numerically and analytically simulate time- and frequency-resolved transport. In this review article, we discuss recent progress in the development of the time-dependent and frequency-dependent non-equilibrium Green function (NEGF) technique. We begin with an overview of the theoretical underpinnings of the underlying Kadanoff-Baym equations and derive the fundamental NEGF equations in the time and frequency domains. We discuss how these methods have been applied to a variety of condensed matter systems such as molecular electronics, nanoscale transistors, and superconductors. In addition, we survey the application of NEGF in fields beyond condensed matter, where it has been used to study thermalization in ultra-cold atoms and to understand leptogenesis in the early universe. Throughout, we pay special attention to the challenges of incorporating contacts and interactions, as the NEGF method is uniquely capable of accounting for such features.
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).
Real-time nonequilibrium Green functions (NEGFs) have been very successfully used to simulate the dynamics of correlated many-particle systems far from equilibrium. However, NEGF simulations are computationally expensive since the effort scales cubically with the simulation duration. Recently, we introduced the G1--G2 scheme that allows for a dramatic reduction to time-linear scaling [N. Schl\"unzen et al., Phys. Rev. Lett. 124, 076601 (2020); J.-P. Joost et al., Phys. Rev. B 101, 245101 (2020)]. Here we tackle another problem: the rapid growth of the computational effort with the system size. In many situations where the system of interest is coupled to a bath, to electric contacts, or to similar macroscopic systems for which a microscopic resolution of the electronic properties is not necessary, efficient simplifications are possible. This is achieved by the introduction of an embedding self-energy---a concept that has been successful in standard NEGF simulations. Here, we demonstrate how the embedding concept can be introduced into the G1--G2 scheme, allowing us to drastically accelerate NEGF embedding simulations. The approach is compatible with all advanced self-energies that can be represented by the G1--G2 scheme [as described in J.-P. Joost et al., Phys. Rev. B 105, 165155 (2022)] and retains the memoryless structure of the equations and their time-linear scaling. As a numerical illustration, we investigate the charge transfer between a Hubbard nanocluster and an additional site which is of relevance for the neutralization of ions in matter.
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.
The Non-Equilibrium Green's Function (NEGF) is the standard formalism for nano-scale electron transport. By recasting the NEGF scattering problem as a linear system of equations whose solution encodes the transmission and reflection amplitudes, we present the first quantum computerized implementation of NEGF. We apply both the Harrow--Hassidim--Lloyd and Variational Quantum Linear Solver algorithms to compute the transmission coefficient T(E) of quantum dragon nanodevices within the single-band tight-binding model. Quantum dragon devices exhibit perfect transmission across the full conducting band regardless of internal disorder. The problem maps onto compact circuits of 3 and 4 total physical qubits for the 2-site and 6-site dragon devices, respectively. A similarity transformation block-diagonalizes the NEGF linear system reducing the Pauli decomposition of the block-encoded matrix. We demonstrate the feasibility of quantum computation by performing ideal and noise-aware simulations and computations on physical IBM quantum processor.