Molecular electronic states near metal surfaces at equilibrium using potential of mean force and numerical renormalization group methods: Hysteresis revisited

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.

Assessing Mixed Quantum-Classical Molecular Dynamics Methods for Nonadiabatic Dynamics of Molecules on Metal Surfaces

Mixed quantum-classical (MQC) methods for simulating the dynamics of molecules at metal surfaces have the potential to accurately and efficiently provide mechanistic insight into reactive processes. Here, we introduce simple two-dimensional models for the scattering of diatomic molecules at metal surfaces based on recently published electronic structure data. We apply several MQC methods to investigate their ability to capture how nonadiabatic effects influence molecule-metal energy transfer during the scattering process. Specifically, we compare molecular dynamics with electronic friction, Ehrenfest dynamics, independent electron surface hopping, and the broadened classical master equation approach. In the case of independent electron surface hopping, we implement a simple decoherence correction approach and assess its impact on vibrationally inelastic scattering. Our results show that simple, low-dimensional models can be used to qualitatively capture experimentally observed vibrational energy transfer and provide insight into the relative performance of different MQC schemes. We observe that all approaches predict similar kinetic energy dependence but return different vibrational energy distributions. Finally, by varying the molecule-metal coupling, we can assess the coupling regime in which some MQC methods become unsuitable.

Nonequilibrium reaction rate theory: Formulation and implementation within the hierarchical equations of motion approach

The study of chemical reactions in environments under nonequilibrium conditions has been of interest recently in a variety of contexts, including current-induced reactions in molecular junctions and scanning tunneling microscopy experiments. In this work, we outline a fully quantum mechanical, numerically exact approach to describe chemical reaction rates in such nonequilibrium situations. The approach is based on an extension of the flux correlation function formalism to nonequilibrium conditions and uses a mixed real and imaginary time hierarchical equations of motion approach for the calculation of rate constants. As a specific example, we investigate current-induced intramolecular proton transfer reactions in a molecular junction for different applied bias voltages and molecule-lead coupling strengths.

On the meaning of Berry force for unrestricted systems treated with mean-field electronic structure

We show that the Berry force as computed by an approximate, mean-field electronic structure can be meaningful if properly interpreted. In particular, for a model Hamiltonian representing a molecular system with an even number of electrons interacting via a two-body (Hubbard) interaction and a spin-orbit coupling, we show that a meaningful nonzero Berry force emerges whenever there is spin unrestriction-even though the Hamiltonian is real-valued and formally the on-diagonal single-surface Berry force must be zero. Moreover, if properly applied, this mean-field Berry force yields roughly the correct asymptotic motion for scattering through an avoided crossing. That being said, within the context of a ground-state calculation, several nuances do arise as far interpreting the Berry force correctly, and as a practical matter, the Berry force diverges near the Coulson-Fischer point (which can lead to numerical instabilities). We do not address magnetic fields here.

Nonadiabatic Molecular Dynamics at Metal Surfaces

Dynamics at molecule-metal interfaces are a subject of intense current interest and come in many different flavors of experiments: gas-phase scattering, chemisorption, electrochemistry, nanojunction transport, and heterogeneous catalysis, to name a few. These dynamics involve nuclear degrees of freedom entangled with many electronic degrees of freedom (in the metal), and as such there is always the possibility for nonadiabatic phenomena to appear: the nuclei do not necessarily need to move slower than the electrons to break the Born-Oppenheimer (BO) approximation. In this Feature Article, we review a set of dynamical methods developed recently to deal with such nonadiabatic phenomena at a metal surface, methods that serve as alternatives to Tully's independent electron surface hopping (IESH) model. In the weak molecule-metal coupling regime, a classical master equation (CME) can be derived and a simple surface hopping approach is proposed to propagate nuclear and electronic dynamics stochastically. In the strong molecule-metal interaction regime, a Fokker-Planck equation can be derived for the nuclear dynamics, with electronic DoFs incorporated into the overall friction and random force. Lastly, a broadened classical master equation (BCME) can interpolate between the weak and strong molecule-metal interactions. Here, we briefly review these methods and the relevant benchmarking data, showing in particular how the methods can be used to calculate nonequilibrium transport properties. We highlight several open questions and pose several avenues for future study.

A comparison of surface hopping approaches for capturing metal-molecule electron transfer: A broadened classical master equation versus independent electron surface hopping

Within a generalized Anderson-Holstein model, we investigate electron transfer rates using two different surface hopping algorithms: a broadened classical master equation (BCME) and independent electron surface hopping (IESH). We find that for large enough bandwidth and density of one electron states, and in the presence of external friction, the IESH results converge to the BCME results for impurity-bath model systems, recovering both relaxation rates and equilibrium populations. Without external friction, however, the BCME and IESH results can strongly disagree, and preliminary evidence suggests that IESH does not always recover the correct equilibrium state. Finally, we also demonstrate that adding an electronic thermostat to IESH does help drive the metallic substrate to the correct equilibrium state, but this improvement can sometimes come at the cost of worse short time dynamics. Overall, our results should be of use for all computational chemists looking to model either gas phase scattering or electrochemical dynamics at a metal interface.

Covalently bonded single-molecule junctions with stable and reversible photoswitched conductivity

Through molecular engineering, single diarylethenes were covalently sandwiched between graphene electrodes to form stable molecular conduction junctions. Our experimental and theoretical studies of these junctions consistently show and interpret reversible conductance photoswitching at room temperature and stochastic switching between different conductive states at low temperature at a single-molecule level. We demonstrate a fully reversible, two-mode, single-molecule electrical switch with unprecedented levels of accuracy (on/off ratio of ~100), stability (over a year), and reproducibility (46 devices with more than 100 cycles for photoswitching and ~10(5) to 10(6) cycles for stochastic switching).