Classical nuclear motion in quantum transport

Electron Dynamics in Open Quantum Systems: The Driven Liouville-von Neumann Methodology within Time-Dependent Density Functional Theory

A first-principles approach to describe electron dynamics in open quantum systems driven far from equilibrium via external time-dependent stimuli is introduced. Within this approach, the driven Liouville-von Neumann methodology is used to impose open boundary conditions on finite model systems whose dynamics is described using time-dependent density functional theory. As a proof of concept, the developed methodology is applied to simple spin-compensated model systems, including a hydrogen chain and a graphitic molecular junction. Good agreement between steady-state total currents obtained via direct propagation and those obtained from the self-consistent solution of the corresponding Sylvester equation indicates the validity of the implementation. The capability of the new computational approach to analyze, from first principles, non-equilibrium dynamics of open quantum systems in terms of temporally and spatially resolved current densities is demonstrated. Future extensions of the approach toward the description of dynamical magnetization and decoherence effects are briefly discussed.

How an electrical current can stabilize a molecular nanojunction

The stability of molecular junctions under transport is of the utmost importance for the field of molecular electronics. This question is often addressed within the paradigm of current-induced heating of nuclear degrees of freedom or current-induced forces acting upon the nuclei. At the same time, an essential characteristic of the failure of a molecular electronic device is its changing conductance - typically from a finite value for the intact device to zero for a device that lost its functionality. In this publication, we focus on the current-induced changes in the molecular conductance, which are inherent to molecular junctions at the limit of mechanical stability. We employ a numerically exact framework based on the hierarchical equations of motion approach, which treats both electronic and nuclear degrees of freedom on an equal footing and does not impose additional assumptions. Studying generic model systems for molecular junctions with dissociative potentials for a wide range of parameters spanning the adiabatic and the nonadiabatic regime, we find that molecular junctions that exhibit a decrease in conductance upon dissociation are more stable than junctions that are more conducting in their dissociated state. This represents a new mechanism that stabilizes molecular junctions under current. Moreover, we identify characteristic signatures in the current of breaking junctions related to the interplay between changes in the conductance and the nuclear configuration and show how these are related to properties of the leads rather than characteristics of the molecule itself.

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.

Nonequilibrium Bond Forces in Single-Molecule Junctions

Passing a current across two touching C₆₀ molecules imposes a nonequilibrium population of bonding and antibonding molecular orbitals, which changes the equilibrium bond character and strength. A current-induced bond force therefore contributes to the total force at chemical-bond distances. The combination of first-principles calculations with scanning probe experiments exploring currents and forces in a wide C₆₀-C₆₀ distance range consistently evidences the presence of current-induced attraction that occurs when the two molecules are on the verge of forming a chemical bond. The unique opportunity to arrange matter at the atomic scale with the atomic force and scanning tunneling microscope tip has enabled closely matching molecular junctions in theory and experiment. The findings consequently represent the first report of current-induced bond forces at the single-molecule level and further elucidate the intimate relation between charge transport and force. The results are relevant to molecular electronics and chemical reactions in the presence of a current.

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.

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.

Nonequilibrium thermodynamics of interacting tunneling transport: variational grand potential, density functional formulation and nature of steady-state forces

The standard formulation of tunneling transport rests on an open-boundary modeling. There, conserving approximations to nonequilibrium Green function or quantum statistical mechanics provide consistent but computational costly approaches; alternatively, the use of density-dependent ballistic-transport calculations (e.g., Lang 1995 Phys. Rev. B 52 5335), here denoted 'DBT', provides computationally efficient (approximate) atomistic characterizations of the electron behavior but has until now lacked a formal justification. This paper presents an exact, variational nonequilibrium thermodynamic theory for fully interacting tunneling and provides a rigorous foundation for frozen-nuclei DBT calculations as a lowest-order approximation to an exact nonequilibrium thermodynamic density functional evaluation. The theory starts from the complete electron nonequilibrium quantum statistical mechanics and I identify the operator for the nonequilibrium Gibbs free energy which, generally, must be treated as an implicit solution of the fully interacting many-body dynamics. I demonstrate a minimal property of a functional for the nonequilibrium thermodynamic grand potential which thus uniquely identifies the solution as the exact nonequilibrium density matrix. I also show that the uniqueness-of-density proof from a closely related Lippmann-Schwinger collision density functional theory (Hyldgaard 2008 Phys. Rev. B 78 165109) makes it possible to express the variational nonequilibrium thermodynamic description as a single-particle formulation based on universal electron-density functionals; the full nonequilibrium single-particle formulation improves the DBT method, for example, by a more refined account of Gibbs free energy effects. I illustrate a formal evaluation of the zero-temperature thermodynamic grand potential value which I find is closely related to the variation in the scattering phase shifts and hence to Friedel density oscillations. This paper also discusses the difference between the here-presented exact thermodynamic forces and the often-used electrostatic forces. Finally the paper documents an inherent adiabatic nature of the thermodynamic forces and observes that these are suited for a nonequilibrium implementation of the Born-Oppenheimer approximation.

Non-adiabatic simulations of current-related structural transformations in metallic nanodevices

One of the less explored aspects of molecular electronics is the effect of current on the mechanical stability of the conducting molecule: charge flow can alter both the geometry and electronic properties of the device, modifying the conductance and giving rise to nonlinear conduction characteristics or conductance switching. We performed a fundamental study on the correlation between the geometry and evolving electronic structure of small Au clusters that were embedded in finite Au wires and subjected to periodic transient currents. Both the current-carrying electronic states and the local electronic structure of the model system were described away from the ground state within a time-dependent Ehrenfest formalism. Non-adiabatic molecular dynamics simulations revealed that clusters undergo structural transformations between several representative geometries that coincide with patterns in cluster charging. The shape changes were enabled by the fluctuations in cluster band filling associated with the current and assisted by current-related forces. Metastable configurations of stable clusters were linked to events of charge trapping on the cluster.

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

The role of the discontinuity of the exchange-correlation potential of density functional theory is studied in the context of electron transport and shown to be intimately related to Coulomb blockade. By following the time evolution of an interacting nanojunction attached to biased leads, we find that, instead of evolving to a steady state, the system reaches a dynamical state characterized by correlation-induced current oscillations. Our results establish a dynamical picture of Coulomb blockade manifesting itself as a periodic sequence of charging and discharging of the nanostructure.

Blowing the fuse: Berry's phase and runaway vibrations in molecular conductors

We examine a molecular bridge connecting two metallic electrodes. We find that an electronic current passing across the bridge can cause a vibrational instability of the molecule, which ultimately can lead to a breakdown of the bridge. This instability is generated by a hitherto never considered mechanism, which surprisingly involves the quantum mechanical phase of the electronic waves, the "Berry phase". This mechanism works for highly conducting bridges, and contrary to breakdown by traditional Joule heating, this instability is deterministic and occurs at certain critical voltages. We demonstrate the new mechanism using state-of-the-art ab initio calculations on realistic molecular bridges.

Current-driven atomic waterwheels

A current induces forces on atoms inside the conductor that carries it. It is now possible to compute these forces from scratch, and to perform dynamical simulations of the atomic motion under current. One reason for this interest is that current can be a destructive force--it can cause atoms to migrate, resulting in damage and in the eventual failure of the conductor. But one can also ask, can current be made to do useful work on atoms? In particular, can an atomic-scale motor be driven by electrical current, as it can be by other mechanisms? For this to be possible, the current-induced forces on a suitable rotor must be non-conservative, so that net work can be done per revolution. Here we show that current-induced forces in atomic wires are not conservative and that they can be used, in principle, to drive an atomic-scale waterwheel.

Coupled electron-phonon transport from molecular dynamics with quantum baths

Based on generalized quantum Langevin equations for the tight-binding wavefunction amplitudes and lattice displacements, electron and phonon quantum transport are obtained exactly using molecular dynamics (MD) in the ballistic regime. The electron-phonon interactions can be handled with a quasi-classical approximation. Both charge and energy transport and their interplay can be studied. We compare the MD results with those of a fully quantum mechanical nonequilibrium Green's function (NEGF) approach for the electron currents. We find a ballistic to diffusive transition of the electron conduction in one-dimensional chains as the chain length increases.

The non-linear response of molecular junctions: the polaron model revisited

A polaron model proposed as a possible mechanism for non-linear conductance (Galperin et al 2005 Nano Lett. 5 125-30) is revisited with the focus on the differences between the weak and strong molecule-lead coupling cases. Within the one-molecule-level model we present an approximate expression for the electronic Green function corresponding to the inelastic transport case, which in the appropriate limits reduces to expressions presented previously for the isolated molecule and for a molecular junction coupled to a slow vibration (static limit). The relevance of considerations based on the isolated molecule limit to understanding properties of molecular junctions is discussed.