We develop nonequilibribrium Green's function based transport theory, which includes effects of nonadiabatic nuclear motion in the calculation of the electric current in molecular junctions. Our approach is based on the separation of slow and fast timescales in the equations of motion for the Green's functions by means of the Wigner representation. Time derivatives with respect to central time serves as a small parameter in the perturbative expansion enabling the computation of nonadiabatic corrections to molecular Green's functions. Consequently, we produce series of analytic expressions for non-adiabatic electronic Green's functions (up to the second order in the central time derivatives); which depend not solely on instantaneous molecular geometry but likewise on nuclear velocities and accelerations. Extended formula for electric current is derived which accounts for the non-adiabatic corrections. This theory is concisely illustrated by the calculations on a model molecular junction.

Molecule-electrode interfaces in molecular electronic junctions are prone to chemical reactions, structural changes, and localized heating effects caused by electric current. These can be exploited for device functionality or may be degrading processes that limit performance and device lifetime. We develop a nonequilibrium Green's function based transport theory in which the central region atoms and, more importantly, atoms on molecule-electrode interfaces are allowed to move. The separation of time-scales of slow nuclear motion and fast electronic dynamics enables the algebraic solution of the Kadanoff-Baym equations in the Wigner space. As a result, analytical expressions for dynamical corrections to the adiabatically computed Green's functions are produced. These dynamical corrections depend not only on the instantaneous molecular geometry but also on the nuclear velocities. To make the theoretical approach fully self-consistent, the same time-separation approach is used to develop expressions for the adiabatic, dissipative, and stochastic components of current-induced forces in terms of adiabatic Green's functions. Using these current induced forces, the equation of motion for the nuclear degrees of freedom is cast in the form of a Langevin equation. The theory is applied to model molecular electronic junctions. We observe that the interplay between the value of the spring constant for the molecule-electrode chemical bond and electronic coupling strength to the corresponding electrode is critical for the appearance of structural instabilities and, consequently, telegraphic switching in the electric current. The range of model parameters is identified to observe structurally stable molecular junctions as well as various different kinds of current-induced telegraphic switching. The interfacial structural instabilities are also quantified based on current noise calculations.

We present quantum electron transport theory that incorporates dynamical effects of motion of atoms on electrode-molecule interfaces in the calculations of the electric current. The theory is based on non-equilibrium Green's functions. We separate time scales in the Green's functions as fast relative time and slow central time. The derivative with respect to the central time serves as a small parameter in the theory. We solve the real-time Kadanoff-Baym equations for molecular Green's functions using Wigner representation and keep terms up to the second order with respect to the central time derivatives. Molecular Green's functions and consequently the electric current are expressed as functions of molecular junction coordinates as well as velocities and accelerations of molecule-electrode interface nuclei. We apply the theory to model a molecular system and study the effects of non-adiabatic nuclear motion on molecular junction conductivity.

Non-equilibrium Green's function theory for non-adiabatic effects in quantum transport Kosov, J.Chem. Phys. 2017, 147, 224109 and J. Chem. Phys. 2018, 149, 044121] is extended to the case of interacting electrons. We consider a general problem of quantum transport of interacting electrons through a central region with dynamically changing geometry. The approach is based on the separation of time scales in the non-equilibrium Green's functions and the use of the Wigner transformation to solve the Kadanoff-Baym equations. The Green's functions and correlation self-energy are non-adiabatically expanded up to the second order central time derivatives. We produce expressions for Green's functions with non-adiabatic corrections and a modified formula for electric current; both depend not only on instantaneous molecular junction geometry but also on nuclear velocities and accelerations. The theory is illustrated by the study of electron transport through a model single-resonant level molecular junction with local electron-electron repulsion and a dynamically changing geometry.

The molecular junction geometry is modeled in terms of nuclear degrees of freedom that are embedded in a stochastic quantum environment of non-equilibrium electrons. The time-evolution of the molecular geometry is governed via a mean force, a frictional force, and a stochastic force, forces arising from many electrons tunneling across the junction for a given nuclear vibration. Conversely, the current-driven nuclear dynamics feed back to the electronic current, which can be captured according to the extended expressions for the current that have explicit dependences on classical nuclear velocities and accelerations. Current-induced nuclear forces and the non-adiabatic electric current are computed using non-equilibrium Green’s functions via a timescale separation solution of Keldysh–Kadanoff–Baym equations in the Wigner space. Applying the theory to molecular junctions demonstrated that non-adiabatic corrections play an important role when nuclear motion is considered non-equilibrium and, in particular, showed that non-equilibrium and equilibrium descriptions of nuclear motion produce significantly different current characteristics. It is observed that non-equilibrium descriptions generally produce heightened conductance profiles relative to the equilibrium descriptions and provide evidence that the effective temperature is an effective measure of the steady-state characteristics. Finally, we observe that the non-equilibrium descriptions of nuclear motion can give rise to the Landauer blowtorch effect via the emergence of multi-minima potential energy surfaces in conjunction with non-uniform temperature profiles. The Landauer blowtorch effect and its impact on the current characteristics, waiting times, and the Fano factor are explored for an effective adiabatic potential that morphs between a single, double, and triple potential as a function of voltage.