Nuclear Dynamics at Molecule-Metal Interfaces: A Pseudoparticle Perspective

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.

Nonadiabatic Dynamics at Metal Surfaces: Fewest Switches Surface Hopping with Electronic Relaxation

A new scheme is proposed for modeling molecular nonadiabatic dynamics near metal surfaces. The charge-transfer character of such dynamics is exploited to construct an efficient reduced representation for the electronic structure. In this representation, the fewest switches surface hopping (FSSH) approach can be naturally modified to include electronic relaxation (ER). The resulting FSSH-ER method is valid across a wide range of coupling strengths as supported by tests applied to the Anderson-Holstein model for electron transfer. Future work will combine this scheme with ab initio electronic structure calculations.

Non-adiabatic effects of nuclear motion in quantum transport of electrons: A self-consistent Keldysh-Langevin study

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.

Chemical Reaction Rates for Systems with Spin-Orbit Coupling and an Odd Number of Electrons: Does Berry's Phase Lead to Meaningful Spin-Dependent Nuclear Dynamics for a Two State Crossing?

Within the context of a simple avoided crossing, we investigate the effect of a complex-valued diabatic coupling in determining spin-dependent rate constants and scattering states. We find that, if the molecular geometry is not linear and the Berry force is not zero, one can find significant spin polarization of the products. This study emphasizes that, when analyzing nonadiabatic reactions with spin orbit coupling (and a complex-valued Hamiltonian), one must consider how Berry force affects nuclear motion-at least in the context of gas phase reactions. Work is currently ongoing as far as extrapolating these conclusions to the condensed phase, where interesting spin selection has been observed in recent years.

A practical ansatz for evaluating the electronic friction tensor accurately, efficiently, and in a nearly black-box format

It is well-known that under conditions of fast electronic equilibration and weak nonadiabaticity, nonadiabatic effects induced by electron-hole pair excitations can be partly incorporated through a frictional force. However, ab initio computation of the electronic friction tensor suffers from numerical instability and usually demands a convergence check. In this study, we present an efficient and accurate interpolation method for computing the electronic friction tensor in a nearly black-box manner as appropriate for molecular dynamics. In almost all cases, our method agrees quite well with the exact friction tensor which is available for several quadratic Hamiltonians. As such, we outperform more conventional approaches that are based on the introduction of a broadening parameter. Future work will implement this interpolation approach within ab initio software packages.

Perspective: How to understand electronic friction

Electronic friction is a correction to the Born-Oppenheimer approximation, whereby nuclei in motion experience a drag in the presence of a manifold of electronic states. The notion of electronic friction has a long history and has been (re-)discovered in the context of a wide variety of different chemical and physical systems including, but not limited to, surface scattering events, surface reactions or chemisorption, electrochemistry, and conduction through molecular-(or nano-) junctions. Over the years, quite a few different forms of electronic friction have been offered in the literature. In this perspective, we briefly review these developments of electronic friction, highlighting the fact that we can now isolate a single, unifying form for (Markovian) electronic friction. We also focus on the role of electron-electron interactions for understanding frictional effects and offer our thoughts on the strengths and weaknesses of using electronic friction to model dynamics in general.

When is electronic friction reliable for dynamics at a molecule–metal interface?

Conditions under which electronic friction dynamics are applicable in the nonadiabatic limit are determined by examination of three model systems.

Coupled wave-packets for non-adiabatic molecular dynamics: a generalization of Gaussian wave-packet dynamics to multiple potential energy surfaces

Coupled wave-packets for non-adiabatic dynamics is a new method for simulation of molecular dynamics on coupled potential energy surfaces, which efficiency and correctly accounts for decoherence and interferences effects.

Accurate simulation of the non-adiabatic dynamics of molecules in excited electronic states is key to understanding molecular photo-physical processes. Here we present a novel method, based on a semi-classical approximation, that is as efficient as the commonly used mean field Ehrenfest or ad hoc surface hopping methods and properly accounts for interference and decoherence effects. This novel method is an extension of Heller's thawed Gaussian wave-packet dynamics that includes coupling between potential energy surfaces. By studying several standard test problems we demonstrate that the accuracy of the method can be systematically improved while maintaining high efficiency. The method is suitable for investigating the role of quantum coherence in the non-adiabatic dynamics of many-atom molecules.