Developments in many-body theory of quantum transport and spectroscopy with non-equilibrium Green's functions and time-dependent density functional theory
The problem of quantum dynamics in open systems has gained attention in recent decades and not the least due to the advances made in quantum transport in molecular systems. The main motivation behind quantum transport and molecular electronics is the futuristic goal to be able at some point to replace, or to complement, the silicon-based technology and to make the electronic devices faster. On a fundamental level, one has to deal with time-dependent processes where electron-electron or electron- phonon interactions are of great importance, and they can cause profound quantitative and qualitative changes on the physical and dynamical properties of electronic systems compared to the non-interacting case. Most of the studies of quantum transport have been focused on the steady-state description while neglecting the short-time dynamics. However, the dynamical effects are of great importance since fast- switching processes play a pivotal role in the operation of future devices. We studied the problem of time-dependent electron transport through the Anderson impurity model by using many-body perturbation theory (MBPT) together with Keldysh Green’s functions as well as with time-dependent density functional theory (TDDFT). These methods were compared with numerically exact time-dependent density-matrix renormalization group (tDMRG) method. We found that the many-body perturbation theory results beyond Hartree-Fock approximation were in close agreement with tDMRG results. In addition we studied the possibility of multistablity in the density and current of an interacting nanoscale junction as well as how to reversibly switch between the multiple solutions in time domain. An accurate theoretical treatment of electron correlation even in as simple model as an interacting electron gas at metallic densities still continues to be a challenge; especially description of features in the photoemission spectra due to electron correlations provides a theoretical challenge. The many-body perturbation theory yields a systematic way to study electron-electron (electron-phonon) correlations in various systems. One of the widely used approximations in MBPT is the GW approximation in which the bare interaction line is replaced with screened interaction line in the first order exchange diagram. The GW approximation gives good estimates for the band gap values close to experimental ones but especially the self-consistent GW approximation has a number of deficiencies like washing out of plasmon features and overestimation of bandwidths compared to experiment. One way to improve GW calculations is to include vertex corrections. Unfortunately, the straightforward inclusion of vertex corrections yields negative spectra in some frequency regions. We developed a diagrammatic approach to construct self-energy approximations with positive spectral properties. Our approach consists of expressing a self-energy of response diagram as a product of half-diagrams after which a minimal set of additional diagrams is identified to construct a perfect square. We applied this method to study vertex corrections in a homogeneous electron gas. In addition we analyzed the diagrammatic content of photocurrent with density functional theory. The expression for the photocurrent was obtained as an integral over the Kohn-Sham spectral function renormalized by effective potentials that depend on the exchange correlation kernel of current density functional theory. The expression for the photocurrent gives us the angular dependence of the photocurrent but it does not provide a direct access to the kinetic energy distribution of the photoelectrons.
The problem of quantum dynamics in open systems has gained attention in recent decades and not the least due to the advances made in quantum transport in molecular systems. The main motivation behind quantum transport and molecular electronics is the futuristic goal to be able at some point to replace, or to complement, the silicon-based technology and to make the electronic devices faster. On a fundamental level, one has to deal with time-dependent processes where electron-electron or electron- phonon interactions are of great importance, and they can cause profound quantitative and qualitative changes on the physical and dynamical properties of electronic systems compared to the non-interacting case. Most of the studies of quantum transport have been focused on the steady-state description while neglecting the short-time dynamics. However, the dynamical effects are of great importance since fast- switching processes play a pivotal role in the operation of future devices. We studied the problem of time-dependent electron transport through the Anderson impurity model by using many-body perturbation theory (MBPT) together with Keldysh Green’s functions as well as with time-dependent density functional theory (TDDFT). These methods were compared with numerically exact time-dependent density-matrix renormalization group (tDMRG) method. We found that the many-body perturbation theory results beyond Hartree-Fock approximation were in close agreement with tDMRG results. In addition we studied the possibility of multistablity in the density and current of an interacting nanoscale junction as well as how to reversibly switch between the multiple solutions in time domain. An accurate theoretical treatment of electron correlation even in as simple model as an interacting electron gas at metallic densities still continues to be a challenge; especially description of features in the photoemission spectra due to electron correlations provides a theoretical challenge. The many-body perturbation theory yields a systematic way to study electron-electron (electron-phonon) correlations in various systems. One of the widely used approximations in MBPT is the GW approximation in which the bare interaction line is replaced with screened interaction line in the first order exchange diagram. The GW approximation gives good estimates for the band gap values close to experimental ones but especially the self-consistent GW approximation has a number of deficiencies like washing out of plasmon features and overestimation of bandwidths compared to experiment. One way to improve GW calculations is to include vertex corrections. Unfortunately, the straightforward inclusion of vertex corrections yields negative spectra in some frequency regions. We developed a diagrammatic approach to construct self-energy approximations with positive spectral properties. Our approach consists of expressing a self-energy of response diagram as a product of half-diagrams after which a minimal set of additional diagrams is identified to construct a perfect square. We applied this method to study vertex corrections in a homogeneous electron gas. In addition we analyzed the diagrammatic content of photocurrent with density functional theory. The expression for the photocurrent was obtained as an integral over the Kohn-Sham spectral function renormalized by effective potentials that depend on the exchange correlation kernel of current density functional theory. The expression for the photocurrent gives us the angular dependence of the photocurrent but it does not provide a direct access to the kinetic energy distribution of the photoelectrons.
We study time-dependent electron transport through an Anderson model. The electronic interactions on the impurity site are included via the self-energy approximations at Hartree-Fock (HF), second Born (2B), GW, and T-matrix levels as well as within a time-dependent density functional (TDDFT) scheme based on the adiabatic Bethe-ansatz local density approximation (ABALDA) for the exchange-correlation potential. The Anderson model is driven out of equilibrium by applying a bias to the leads, and its nonequilibrium dynamics is determined by real-time propagation. The time-dependent currents and densities are compared to benchmark results obtained with the time-dependent density matrix renormalization group (tDMRG) method. Many-body perturbation theory beyond HF gives results in close agreement with tDMRG, especially within the 2B approximation. We find that the TDDFT approach with the ABALDA approximation produces accurate results for the densities on the impurity site, but overestimates the currents. This problem is found to have its origin in an overestimation of the lead densities, which indicates that the exchange-correlation potential must attain nonzero values in the leads.
In this thesis the electron transport through interacting quantum wires is studied.The time-dependent properties of the system are investigated by using the non-equilibrium Green function technique.This technique provides a way to include ultra-fast timedependent processes, arbitrary strong external fields and also brings forward an effective treatment of electron-electron and electron-phonon interactions.The electronelectron interactions are included via various self-energy approximations at Hartree-Fock, second-order Born and GW levels.These approximations are called conserving since they fulfill the physical conservation laws like the charge conservation law.The interacting central system is coupled to one-dimensional tight-binding leads, which are considered effectively as noninteracting.The system is perturbed by a bias voltage after which the time-evolution of the system is followed by Kadanoff-Baym equations.The non-equilibrium steady-state value of the current depends on the used self-energy approximation, which indicates the importance of the history-dependence in the calculation of current.This is an important observation in the aspect of time-dependent density functional theory, since it is one indication for the need of functionals with memory.The steady-state current is independent of the initial state, which is in accordance with the memory-loss theorem.The transient current shows small oscillations, which are consequence of the density oscillations especially in the first and last site.The density oscillations in the sites are effect of the electron transition from the lead to the central system and vice versa.
In this thesis we focus on improvements of the description of the electron-electron correlation effects in nonequilibrium nanosystems. We mainly focus on developments of two nonequilibrium methods, namely the formalism of Nonequilibrium Green’s Function and Time Dependent Density Functional Theory and we explore the possibility to improve existing approximations in these theories. A smaller part of the thesis is devoted to the Exact Diagonalization method which provides a numerically exact description of small systems.Paper I: We review the current methods for description of correlated materials in nonequilibrium and their connection to pump-probe spectroscopy.Paper II: We propose a hybrid method for the real time dynamics of strongly correlated materials which includes memory effects beyond the adiabatic local density approximation.Paper III: We study the dynamics of desorption of a molecule from a surface with different levels of approximation for both the nuclear and the electronic part. We compare a full quantum mechanical treatment to the Ehrenfestapproximation for the molecule and perturbative approximations for the electrons.Paper IV: We develop a theory of current-induced forces within Adiabatic Ehrenfest Dynamics which includes effects of electron-electron interactions. We study a dependence of the electronic friction on interaction strength.We also benchmark it against nonadiabatic Ehrenfest dynamics.Paper V: We study the competition of interaction and disorder in systems with steady state currents - in transport and ring geometries. We exactly define the exchange-correlation screening of the disorder by the interactioneffects via Kohn–Sham construction of DFT.Paper VI: We study a competition between Kondo and RKKY interaction in small clusters of Periodic Anderson Model (ring geometries), we construct a nonequilibrium Doniach-phase like diagram. We then determine anoptimal pulse to induce transitions with the highest fidelity.
We propose a time-dependent many-body approach to study the short-time dynamics of correlated electrons in quantum transport through nanoscale systems contacted to metallic leads. This approach is based on the time propagation of the Kadanoff-Baym equations for the nonequilibrium many-body Green's function of open and interacting systems out of equilibrium. An important feature of the method is that it takes full account of electronic correlations and embedding effects in the presence of time-dependent external fields, while at the same time satisfying the charge conservation law. The method further extends the Meir-Wingreen formula to the time domain for initially correlated states. We study the electron dynamics of a correlated quantum wire attached to two-dimensional leads exposed to a sudden switch on of a bias voltage using conserving many-body approximations at Hartree-Fock, second Born and GW level. We obtain detailed results for the transient currents, dipole moments, spectral functions, charging times, and the many-body screening of the quantum wire as well as for the time-dependent density pattern in the leads, and we show how the time dependence of these observables provides a wealth of information on the energy level structure of the quantum wire out of equilibrium. For moderate interaction strengths the second Born and GW results are in excellent agreement at all times. We find that many-body effects beyond the Hartree-Fock approximation have a large effect on the qualitative behavior of the system and lead to a bias-dependent gap closing and quasiparticle broadening, shortening of the transient times and washing out of the step features in the current-voltage curves.
The most of the experimental measurable quantities in many-body electron system, such as the optical absorption, luminescent spectra and exciton effect, are related to the excitation states description. The density functional theory with local density approximation (DFT-LDA) as a ground state theory, i.e. the solutions based on the Kohn-Sham equation is a powerful tool for studying the ground state properties in many-particle systems. However, the first principles computation of excited states is more complexity than ground-state calculations. A key problem is that the exchange-correlation interaction in excited states is differ from the ground-states. In recent years, nevertheless, several electronic-excitation theories have been developed. The most important theoretical and computational method include the many-body perturbation theory which is based on the quasi-particle concept and the Green function equation, the time-dependent DFT and the Bethe-Salpeter equation for describing the electron-hole interaction. Among them the central ingredient is an electron's self-energy operator Σ that describes the exchange correlation effect beyond Hartree approximation. The implementations of the above theories are unavoidable to introduce some of approximations. A good approximation for Σ is the GW approach by Hedin. It is shown from computation simulations for many real condensed matter systems that the GW approach is a successful method for electronic-excitation problems. Another effective method for excited-states is so-called screen-exchange local-density approximation (sX-LDA). It is a combination of LDA and Hartree-Fock (HF)theory, in which the intrinsic local screen-exchange interaction has been replaced by a non-local one and a generalized Kohn-Sham equation (GKS) is given. Based on a PWscf package, the sX-LDA-PWscf method has been tested and indicates that an available approach for treating band-gap problem in semiconductors. In this paper, the many-body theories for the excitation process, its development and signification are reviewed. The relation and differences consist in different theoretical methods are discussed. On the based of above discussions, the applications of many-body theory to the band-band transition (band-gap underestimation problem) and exciton effect in semiconductors are presented.
For the calculation of neutral excitations, time-dependent density functional theory (TDDFT) is an exact reformulation of the many-body time-dependent Schrödinger equation, based on knowledge of the density instead of the many-body wavefunction. The density can be determined in an efficient scheme by solving one-particle non-interacting Schrödinger equations—the Kohn–Sham equations. The complication of the problem is hidden in the—unknown—time-dependent exchange and correlation potential that appears in the Kohn–Sham equations and for which it is essential to find good approximations. Many approximations have been suggested and tested for finite systems, where even the very simple adiabatic local-density approximation (ALDA) has often proved to be successful. In the case of solids, ALDA fails to reproduce optical absorption spectra, which are instead well described by solving the Bethe–Salpeter equation of many-body perturbation theory (MBPT). On the other hand, ALDA can lead to excellent results for loss functions (at vanishing and finite momentum transfer). In view of this and thanks to recent successful developments of improved linear-response kernels derived from MBPT, TDDFT is today considered a promising alternative to MBPT for the calculation of electronic spectra, even for solids. After reviewing the fundamentals of TDDFT within linear response, we discuss different approaches and a variety of applications to extended systems.