arXiv:0906.2136 · Physical Review B

Kadanoff-Baym approach to quantum transport through interacting nanoscale systems: From the transient to the steady-state regime

Petri Myöhänen, Adrian Stan, Gianluca Stefanucci, Robert van Leeuwen

Published Sep 8, 2009• 222 citations

Abstract

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.

DOI: 10.1103/physrevb.80.115107

Related papers

Many-body approach to time-dependent quantum transport

Anna-Maija Uimonen

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.

Jan 1, 20090 citations

A many-body approach to quantum transport dynamics: Initial correlations and memory effects

Petri Myöhänen, Adrian Stan, Gianluca Stefanucci +1 more

We study time-dependent quantum transport through a correlated double quantum dot (DQD) model system by means of time propagation of the nonequilibrium many-body Green's function. The theory is an extension of the Kadanoff-Baym approach for finite inhomogeneous systems ( Phys. Rev. Lett. , 98 (2007) 153004) to open inhomogeneous systems and generalizes the Meir-Wingreen formula to include initial correlations and memory effects. Important features of the theory are 1) the possibility to study the ultrafast dynamics of transients and other time-dependent regimes and 2) the inclusion of exchange and correlation effects in a conserving approximation scheme. We calculate time-dependent local currents and densities for different many-body approximations and highlight the role of initial correlations and memory effects on the transient dynamics. Furthermore we show that coherent charge oscillations on the DQD are strongly affected by the confined Coulomb interaction and can be directly related to the local equilibrium spectral density.

arXiv:0808.3483Europhysics Letters (EPL)Nov 13, 2008146 citations

Kadanoff-Baym approach to time-dependent quantum transport in AC and DC fields

Petri Myöhänen, Adrian Stan, Gianluca Stefanucci +1 more

We have developed a method based on the embedded Kadanoff-Baym equations to study the time evolution of open and inhomogeneous systems. The equation of motion for the Green's function on the Keldysh contour is solved using different conserving many-body approximations for the self-energy. Our formulation incorporates basic conservation laws, such as particle conservation, and includes both initial correlations and initial embedding effects, without restrictions on the time-dependence of the external driving field. We present results for the time-dependent density, current and dipole moment for a correlated tight binding chain connected to one-dimensional non-interacting leads exposed to DC and AC biases of various forms. We find that the self-consistent 2B and GW approximations are in extremely good agreement with each other at all times, for the long-range interactions that we consider. In the DC case we show that the oscillations in the transients can be understood from interchain and lead-chain transitions in the system and find that the dominant frequency corresponds to the HOMO-LUMO transition of the central wire. For AC biases with odd inversion symmetry odd harmonics to high harmonic order in the driving frequency are observed in the dipole moment, whereas for asymmetric applied bias also even harmonics have considerable intensity. In both cases we find that the HOMO-LUMO transition strongly mixes with the harmonics leading to harmonic peaks with enhanced intensity at the HOMO-LUMO transition energy.

arXiv:1006.2912Journal of Physics Conference SeriesApr 1, 201027 citations

Developments in many-body theory of quantum transport and spectroscopy with non-equilibrium Green's functions and time-dependent density functional theory

Anna-Maija Uimonen

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 ﬁrst 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 deﬁciencies 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 identiﬁed 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.

Jan 1, 20152 citations

Kadanoff-Baym description of Hubbard clusters out of equilibrium: performance of many-body schemes, correlation-induced damping and multiple steady states

M. Puig von Friesen, Claudio Verdozzi, C.‐O. Almbladh

We present in detail a method we recently introduced (PRL. 103, 176404 (2009)) to describe finite systems in and out of equilibrium, where the evolution in time is performed via the Kadanoff-Baym Equations (KBE) within Many-Body Perturbation Theory (MBPT). The main property we analyze is the time-dependent density. We also study is the exchange-correlation potential of TDDFT, obtained via reverse engineering from the time-dependent density. Our systems consist of small, strongly correlated clusters, described by a Hubbard Hamiltonian within the Hartree-Fock, second Born, GW and T-matrix approximations. We compare the results from the KBE dynamics to those from exact numerical solutions. The outcome of our comparisons is that, among the many-body schemes considered, the T-matrix approximation is overall superior at all electron densities. Such comparisons permit a general assessment of the whole idea of applying MBPT, in the KBE sense, to finite systems. A striking outcome of our analysis is that when the system evolves under a strong external field, the KBE develop a steady-state solution as a consequence of a correlation-induced damping. This damping is present both in isolated (finite) systems, where it is purely artificial, as well as in clusters contacted to (infinite) macroscopic leads. To illustrate this point we present selected results for a system coupled to contacts within the T-matrix and second Born approximation. The extensive characterization we performed indicates that this behavior is present whenever approximate self energies, based upon infinite partial summations, are used. A second important result is that, for isolated clusters, the steady state reached is not unique but depends on how one switches on the external field. This may also true for clusters connected to leads.

arXiv:1004.3662Apr 21, 20100 citations

Real-time switching between multiple steady-states in quantum transport

A.-M. Uimonen, E Khosravi, Gianluca Stefanucci +3 more

We study transport through an interacting model system consisting of a central correlated site coupled to finite bandwidth tight-binding leads, which are considered as effectively noninteracting. Its nonequilibrium properties are determined by real-time propagation of the Kadanoff-Baym equations after applying a bias voltage to the system. The electronic interactions on the central site are incorporated by means of self-energy approximations at Hartree-Fock, second Born and GW level. We investigate the conditions under which multiple steady-state solutions occur within different self-energy approximations, and analyze in detail the nature of these states from an analysis of their spectral functions. At the Hartree-Fock level at least two stable steady-state solutions with different densities and currents can be found. By applying a gate voltage-pulse at a given time we are able to switch between these solutions. With the same parameters we find only one steady-state solution when the self-consistent second Born and GW approximations are considered. We therefore conclude that treatment of many-body interactions beyond mean-field can destroy bistability and lead to qualitatively different results as compared those at mean-field level.

arXiv:1106.5631Journal of Physics Conference SeriesApr 1, 201011 citations