Dynamical Coulomb blockade and the derivative discontinuity of time-dependent density functional theory

A projection-based reduced-order method for electron transport problems with long-range interactions

Long-range interactions play a central role in electron transport. At the same time, they present a challenge for direct computer simulations since sufficiently large portions of the bath have to be included in the computation to accurately compute the Coulomb potential. This article presents a reduced-order approach by deriving an open quantum model for the reduced density matrix. To treat the transient dynamics, the problem is placed in a reduced-order framework. The dynamics described by the Liouville-von Neumann equation is projected to subspaces using a Petrov-Galerkin projection. In order to recover the global electron density profile as a vehicle to compute the Coulomb potential, we propose a domain decomposition approach, where the computational domain also includes segments of the bath that are selected using logarithmic grids. This approach leads to a multi-component self-energy that enters the effective Hamiltonian. We demonstrate the accuracy of the reduced model using a molecular junction built from lithium chains.

Electronic transport in molecular junctions: The generalized Kadanoff-Baym ansatz with initial contact and correlations

The generalized Kadanoff-Baym ansatz (GKBA) offers a computationally inexpensive approach to simulate out-of-equilibrium quantum systems within the framework of nonequilibrium Green's functions. For finite systems, the limitation of neglecting initial correlations in the conventional GKBA approach has recently been overcome [Karlsson et al., Phys. Rev. B 98, 115148 (2018)]. However, in the context of quantum transport, the contacted nature of the initial state, i.e., a junction connected to bulk leads, requires a further extension of the GKBA approach. In this work, we lay down a GKBA scheme that includes initial correlations in a partition-free setting. In practice, this means that the equilibration of the initially correlated and contacted molecular junction can be separated from the real-time evolution. The information about the contacted initial state is included in the out-of-equilibrium calculation via explicit evaluation of the memory integral for the embedding self-energy, which can be performed without affecting the computational scaling with the simulation time and system size. We demonstrate the developed method in carbon-based molecular junctions, where we study the role of electron correlations in transient current signatures.

Time-resolved impurity-invisibility in graphene nanoribbons

We investigate time-resolved charge transport through graphene nanoribbons supplemented with adsorbed impurity atoms. Depending on the location of the impurities with respect to the hexagonal carbon lattice, the transport properties of the system may become invisible to the impurity due to the symmetry properties of the binding mechanism. This motivates a chemical sensing device since dopants affecting the underlying sublattice symmetry of the pristine graphene nanoribbon introduce scattering. Using the time-dependent Landauer-Büttiker formalism, we extend the stationary current-voltage picture to the transient regime, where we observe how the impurity invisibility takes place at sub-picosecond time scales further motivating ultrafast sensor technology. We further characterize time-dependent local charge and current profiles within the nanoribbons, and we identify rearrangements of the current pathways through the nanoribbons due to the impurities. We finally study the behavior of the transients with ac driving which provides another way of identifying the lattice-symmetry breaking caused by the impurities.

Communication: Gibbs phenomenon and the emergence of the steady-state in quantum transport

Simulations are increasingly employing explicit reservoirs-internal, finite regions-to drive electronic or particle transport. This naturally occurs in simulations of transport via ultracold atomic gases. Whether the simulation is numerical or physical, these approaches rely on the rapid development of the steady state. We demonstrate that steady state formation is a manifestation of the Gibbs phenomenon well-known in signal processing and in truncated discrete Fourier expansions. Each particle separately develops into an individual steady state due to the spreading of its wave packet in energy. The rise to the steady state for an individual particle depends on the particle energy-and thus can be slow-and ringing oscillations appear due to filtering of the response through the electronic bandwidth. However, the rise to the total steady state-the one from all particles-is rapid, with time scale π/W, where W is the bandwidth. Ringing oscillations are now also filtered through the bias window, and they decay with a higher power. The Gibbs constant-the overshoot of the first ring-can appear in the simulation error. These results shed light on the formation of the steady state and support the practical use of explicit reservoirs to simulate transport at the nanoscale or using ultracold atomic lattices.

Nonadiabatic Dynamics in Single-Electron Tunneling Devices with Time-Dependent Density-Functional Theory

We simulate the dynamics of a single-electron source, modeled as a quantum dot with on-site Coulomb interaction and tunnel coupling to an adjacent lead in time-dependent density-functional theory. Based on this system, we develop a time-nonlocal exchange-correlation potential by exploiting analogies with quantum-transport theory. The time nonlocality manifests itself in a dynamical potential step. We explicitly link the time evolution of the dynamical step to physical relaxation timescales of the electron dynamics. Finally, we discuss prospects for simulations of larger mesoscopic systems.

Extreme electron transport suppression in siloxane ring-based molecular devices

Siloxane ring-based molecular devices possess excessive transport suppression and size-dependent transport decay, based on an analysis of electronic coupling.

Transport through correlated systems with density functional theory

We present recent advances in density functional theory (DFT) for applications in the field of quantum transport, with particular emphasis on transport through strongly correlated systems. We review the foundations of the popular Landauer-Büttiker(LB)  +  DFT approach. This formalism, when using approximations to the exchange-correlation (xc) potential with steps at integer occupation, correctly captures the Kondo plateau in the zero bias conductance at zero temperature but completely fails to capture the transition to the Coulomb blockade (CB) regime as the temperature increases. To overcome the limitations of LB  +  DFT, the quantum transport problem is treated from a time-dependent (TD) perspective using TDDFT, an exact framework to deal with nonequilibrium situations. The steady-state limit of TDDFT shows that in addition to an xc potential in the junction, there also exists an xc correction to the applied bias. Open shell molecules in the CB regime provide the most striking examples of the importance of the xc bias correction. Using the Anderson model as guidance we estimate these corrections in the limit of zero bias. For the general case we put forward a steady-state DFT which is based on one-to-one correspondence between the pair of basic variables, steady density on and steady current across the junction and the pair local potential on and bias across the junction. Like TDDFT, this framework also leads to both an xc potential in the junction and an xc correction to the bias. Unlike TDDFT, these potentials are independent of history. We highlight the universal features of both xc potential and xc bias corrections for junctions in the CB regime and provide an accurate parametrization for the Anderson model at arbitrary temperatures and interaction strengths, thus providing a unified DFT description for both Kondo and CB regimes and the transition between them.

Achieving Predictive Description of Molecular Conductance by Using a Range-Separated Hybrid Functional

The conductance of molecular bridges tends to be overestimated by computational studies in comparison to measured values. While this well-established trend may be related to difficulties for achieving robust bridges, the employed computational scheme can also contribute to this tendency. In particular, caveats of the traditional functionals employed in first-principles-based calculations can lead to discrepancies reflected in exaggerated conductance. Here, we show that by employing a range-separated hybrid functional the calculated values are within the same order as the measured conductance for all four considered cases. On the other hand, with B3LYP, which is a widely used functional, the calculated values greatly overestimate the conductance (by about 1-2 orders of magnitude). The improved description of the conductance with a RSH functional builds on achieving a physically meaningful treatment of the quasi particles associated with the frontier orbitals.

Steady-State Density Functional Theory for Finite Bias Conductances

In the framework of density functional theory, a formalism to describe electronic transport in the steady state is proposed which uses the density on the junction and the steady current as basic variables. We prove that, in a finite window around zero bias, there is a one-to-one map between the basic variables and both local potential on as well as bias across the junction. The resulting Kohn-Sham system features two exchange-correlation (xc) potentials, a local xc potential, and an xc contribution to the bias. For weakly coupled junctions the xc potentials exhibit steps in the density-current plane which are shown to be crucial to describe the Coulomb blockade diamonds. At small currents these steps emerge as the equilibrium xc discontinuity bifurcates. The formalism is applied to a model benzene junction, finding perfect agreement with the orthodox theory of Coulomb blockade.

The Hubbard dimer: a density functional case study of a many-body problem

This review explains the relationship between density functional theory and strongly correlated models using the simplest possible example, the two-site Hubbard model. The relationship to traditional quantum chemistry is included. Even in this elementary example, where the exact ground-state energy and site occupations can be found analytically, there is much to be explained in terms of the underlying logic and aims of density functional theory. Although the usual solution is analytic, the density functional is given only implicitly. We overcome this difficulty using the Levy-Lieb construction to create a parametrization of the exact function with negligible errors. The symmetric case is most commonly studied, but we find a rich variation in behavior by including asymmetry, as strong correlation physics vies with charge-transfer effects. We explore the behavior of the gap and the many-body Green's function, demonstrating the 'failure' of the Kohn-Sham (KS) method to reproduce the fundamental gap. We perform benchmark calculations of the occupation and components of the KS potentials, the correlation kinetic energies, and the adiabatic connection. We test several approximate functionals (restricted and unrestricted Hartree-Fock and Bethe ansatz local density approximation) to show their successes and limitations. We also discuss and illustrate the concept of the derivative discontinuity. Useful appendices include analytic expressions for density functional energy components, several limits of the exact functional (weak- and strong-coupling, symmetric and asymmetric), various adiabatic connection results, proofs of exact conditions for this model, and the origin of the Hubbard model from a minimal basis model for stretched H2.

Nonequilibrium density matrix in quantum open systems: generalization for simultaneous heat and charge steady-state transport

We suggest a generalization of the expression of the nonequilibrium (NE) density matrix obtained by Hershfield's method for the cases where both heat and charge steady-state currents are present in a quantum open system. The finite-size quantum system, connected to two temperature and particle reservoirs, is driven out of equilibrium by the presence of both a temperature gradient and a chemical potential gradient between the two reservoirs. We show that the NE density matrix is given by a generalized Gibbs-like ensemble and is in full agreement with the general results of the McLennan-Zubarev nonequilibrium ensembles. The extra nonequilibrium terms are related to the entropy production in the system and characterize the fluxes of heat and particle. An explicit example, for the lowest-order expansion, is provide for a model system of noninteracting fermions.

First-principles investigation of quantum transport through an endohedral N@C60 in the Coulomb blockade regime

We report an investigation of Coulomb blockade transport through an endohedral N@C60 weakly coupled with aluminum leads, employing the first-principles method combined with the Keldysh non-equilibrium Green's function derived from the equation of motion beyond the Hartree-Fock approximation. The differential conductance characteristics of the molecular device are calculated within the Coulomb blockade regime, which shows the Coulomb diamond as observed experimentally. When the gate voltage is less than that of the degeneracy point, there are two peaks in the differential conductance with an excited state induced by the change of the exchange interaction between the spin of C60 and the encapsulated nitrogen atom due to the transition from N@C(1-)(60) to N@C(2-)(60), while for a gate voltage larger than that of the degeneracy point, no excited state is available due to the quenching of exchange energy. As a result, there is only one Coulomb blockade peak in the differential conductance from the electron tunneling through the highest energy level below the Fermi level. Our first-principles results are in good agreement with experimental data obtained by an endohedral N@C60 molecular device.

Derivative discontinuity in the strong-interaction limit of density-functional theory

We generalize the exact strong-interaction limit of the exchange-correlation energy of Kohn-Sham density functional theory to open systems with fluctuating particle numbers. When used in the self-consistent Kohn-Sham procedure on strongly interacting systems, this functional yields exact features crucial for important applications such as quantum transport. In particular, the steplike structure of the highest-occupied Kohn-Sham eigenvalue is very well captured, with accurate quantitative agreement with exact many-body chemical potentials. While it can be shown that a sharp derivative discontinuity is present only in the infinitely strongly correlated limit, at finite correlation regimes we observe a slightly smoothened discontinuity, with qualitative and quantitative features that improve with increasing correlation. From the fundamental point of view, our results obtain the derivative discontinuity without making the assumptions used in its standard derivation, offering independent support for its existence.

Kondo memory in driven strongly correlated quantum dots

We investigate the real-time current response of strongly correlated quantum dot systems under sinusoidal driving voltages. By means of an accurate hierarchical equations of motion approach, we demonstrate the presence of prominent memory effects induced by the Kondo resonance on the real-time current response. These memory effects appear as distinctive hysteresis line shapes and self-crossing features in the dynamic current-voltage characteristics, with concomitant excitation of odd-number overtones. They emerge as a cooperative effect of quantum coherence-due to inductive behavior-and electron correlations-due to the Kondo resonance. We also show the suppression of memory effects and the transition to classical behavior as a function of temperature. All these phenomena can be observed in experiments and may lead to novel quantum memory applications.

Dynamical correction to linear Kohn-Sham conductances from static density functional theory

For molecules weakly coupled to leads the exact linear Kohn-Sham (KS) conductance can be orders of magnitude larger than the true linear conductance due to the lack of dynamical exchange-correlation (xc) corrections. In this work we show how to incorporate dynamical effects in KS transport calculations. The only quantity needed is the static xc potential in the molecular junction. Our scheme provides a comprehensive description of Coulomb blockade without breaking the spin symmetry. This is explicitly demonstrated in single-wall nanotubes where the corrected conductance is in good agreement with experimental data whereas the KS conductance fails dramatically.

Exact exchange-correlation potential of an ionic Hubbard model with a free surface

In Kohn-Sham density functional theory (DFT) the interacting electron problem is mapped into a noninteracting problem in an effective potential v_KS. It is known that the charge gap of the interacting system is different from the gap of the effective problem due to a jump Δ_xc in v_KS when an electron is added but its magnitude and its role in the ubiquitous discrepancy between the experimental gaps and approximate DFT computations is poorly understood. Here we compute the exact v_KS of a strongly interacting one-dimensional lattice model which can be driven from an ionic to a Mott insulating state. Presence of a “vacuum” region allows to determine the absolute value of v_KS. We show that in the ionic regime Δ_xc is determined by nearest-neighbor interaction, while in the Mott regime Δ_xc is determined by on-site Hubbard interaction.

Time-dependent electron transport through a strongly correlated quantum dot: multiple-probe open-boundary conditions approach

We present a time-dependent study of electron transport through a strongly correlated quantum dot, which combines adiabatic lattice density functional theory in the Bethe ansatz local-density approximation (BALDA) to the Hubbard model, with the multiple-probe battery method for open-boundary simulations in the time domain. In agreement with the recently proposed dynamical picture of Coulomb blockade, a characteristic driven regime, defined by regular current oscillations, is demonstrated for a certain range of bias voltages. We further investigate the effects of systematically improving the approximation for the electron-electron interaction at the dot site (going from non-interacting, through Hartree-only to adiabatic BALDA) on the transmission spectrum and the I-V characteristics. In particular, a negative differential conductance is obtained at large bias voltages and large Coulomb interaction strengths. This is attributed to the combined effect of the electron-electron interaction at the dot and the finite bandwidth of the electrodes.

Hierarchical Liouville-space approach for accurate and universal characterization of quantum impurity systems

A hierarchical equations of motion based numerical approach is developed for accurate and efficient evaluation of dynamical observables of strongly correlated quantum impurity systems. This approach is capable of describing quantitatively Kondo resonance and Fermi-liquid characteristics, achieving the accuracy of the latest high-level numerical renormalization group approach, as demonstrated on single-impurity Anderson model systems. Its application to a two-impurity Anderson model results in differential conductance versus external bias, which correctly reproduces the continuous transition from Kondo states of individual impurity to singlet spin states formed between two impurities. The outstanding performance on characterizing both equilibrium and nonequilibrium properties of quantum impurity systems makes the hierarchical equations of motion approach potentially useful for addressing strongly correlated lattice systems in the framework of dynamical mean-field theory.

Universal dynamical steps in the exact time-dependent exchange-correlation potential

We show that the exact exchange-correlation potential of time-dependent density-functional theory displays dynamical step structures that have a spatially nonlocal and time nonlocal dependence on the density. Using one-dimensional two-electron model systems, we illustrate these steps for a range of nonequilibrium dynamical situations relevant for modeling of photochemical or physical processes: field-free evolution of a nonstationary state, resonant local excitation, resonant complete charge transfer, and evolution under an arbitrary field. A lack of these steps in the usual approximations yields inaccurate dynamics, for example, predicting faster dynamics and incomplete charge transfer.

Strong correlation in Kohn-Sham density functional theory

We use the exact strong-interaction limit of the Hohenberg-Kohn energy density functional to approximate the exchange-correlation energy of the restricted Kohn-Sham scheme. Our approximation corresponds to a highly nonlocal density functional whose functional derivative can be easily constructed, thus transforming exactly, in a physically transparent way, an important part of the electron-electron interaction into an effective local one-body potential. We test our approach on quasi-one-dimensional systems, showing that it captures essential features of strong correlation that restricted Kohn-Sham calculations using the currently available approximations cannot describe.

Persistent current and Drude weight for the one-dimensional Hubbard model from current lattice density functional theory

The Bethe ansatz local density approximation (LDA) to lattice density functional theory (LDFT) for the one-dimensional repulsive Hubbard model is extended to current-LDFT (CLDFT). The transport properties of mesoscopic Hubbard rings threaded by a magnetic flux are then systematically investigated by this scheme. In particular we present calculations of ground state energies, persistent currents and Drude weights for both a repulsive homogeneous and a single impurity Hubbard model. Our results for the ground state energies in the metallic phase compare favorably well with those obtained with numerically accurate many-body techniques. Also the dependence of the persistent currents on the Coulomb and the impurity interaction strength, and on the ring size are all well captured by LDA-CLDFT. Our study demonstrates the value of CLDFT in describing the transport properties of one-dimensional correlated electron systems. As its computational overheads are rather modest, we propose this method as a tool for studying problems where both disorder and interaction are present.

Towards a description of the Kondo effect using time-dependent density-functional theory

We demonstrate that the zero-temperature conductance of the Anderson model can be calculated within the Landauer formalism combined with static density-functional theory. The proposed approximate functional is based on finite-temperature density-functional theory and yields the exact Kohn-Sham potential at the particle-hole symmetric point. Furthermore, in the limit of zero temperature it correctly exhibits a derivative discontinuity which is shown to be essential to reproduce the conductance plateau. On the other hand, at the Kondo temperature the exact Kohn-Sham conductance overestimates the real one by an order of magnitude. To understand the failure of density-functional theory, we resort to its time-dependent version and conclude that the suppression of the Kondo resonance must be attributed to dynamical exchange-correlation corrections.

Stability analysis of multiple nonequilibrium fixed points in self-consistent electron transport calculations

We present a method to perform stability analysis of nonequilibrium fixed points appearing in self-consistent electron transport calculations. The nonequilibrium fixed points are given by the self-consistent solution of stationary, nonlinear kinetic equation for single-particle density matrix. We obtain the stability matrix by linearizing the kinetic equation around the fixed points and analyze the real part of its spectrum to assess the asymptotic time behavior of the fixed points. We derive expressions for the stability matrices within Hartree-Fock and linear response adiabatic time-dependent density functional theory. The stability analysis of multiple fixed points is performed within the nonequilibrium Hartree-Fock approximation for the electron transport through a molecule with a spin-degenerate single level with local Coulomb interaction.

Active Brownian motion in threshold distribution of a Coulomb blockade model

Randomly distributed offset charges affect the nonlinear current-voltage property via the fluctuation of the threshold voltage above which the current flows in an array of a Coulomb blockade (CB). We analytically derive the distribution of the threshold voltage for a model of one-dimensional locally coupled CB arrays and propose a general relationship between conductance and distribution. In addition, we show that the distribution for a long array is equivalent to the distribution of the number of upward steps for aligned objects of different heights. The distribution satisfies a novel Fokker-Planck equation corresponding to active Brownian motion. The feature of the distribution is clarified by comparing it with the Wigner and Ornstein-Uhlenbeck processes. It is not restricted to the CB model but is instructive in statistical physics generally.