Perspective: How to understand electronic friction

Nonadiabatic dynamics of molecules interacting with metal surfaces: Extending the hierarchical equations of motion and Langevin dynamics approach to position-dependent metal-molecule couplings

Electronic friction and Langevin dynamics is a popular mixed quantum-classical method for simulating the nonadiabatic dynamics of molecules interacting with metal surfaces, as it can be computationally more efficient than fully quantum approaches. In this work, we extend the theory of electronic friction within the hierarchical equations of motion formalism to models with a position-dependent metal-molecule coupling. We show that the addition of a position-dependent metal-molecule coupling adds new contributions to the electronic friction and other forces, which are highly relevant for many physical processes. Our expressions for the electronic forces within the Langevin equation are valid both in and out of equilibrium and for molecular models containing strong interactions. We demonstrate the approach by applying it to different models of interest.

i-PI 3.0: A flexible and efficient framework for advanced atomistic simulations

Atomic-scale simulations have progressed tremendously over the past decade, largely thanks to the availability of machine-learning interatomic potentials. These potentials combine the accuracy of electronic structure calculations with the ability to reach extensive length and time scales. The i-PI package facilitates integrating the latest developments in this field with advanced modeling techniques thanks to a modular software architecture based on inter-process communication through a socket interface. The choice of Python for implementation facilitates rapid prototyping but can add computational overhead. In this new release, we carefully benchmarked and optimized i-PI for several common simulation scenarios, making such overhead negligible when i-PI is used to model systems up to tens of thousands of atoms using widely adopted machine learning interatomic potentials, such as Behler-Parinello, DeePMD, and MACE neural networks. We also present the implementation of several new features, including an efficient algorithm to model bosonic and fermionic exchange, a framework for uncertainty quantification to be used in conjunction with machine-learning potentials, a communication infrastructure that allows for deeper integration with electronic-driven simulations, and an approach to simulate coupled photon-nuclear dynamics in optical or plasmonic cavities.

First-Principles Nonadiabatic Dynamics of Molecules at Metal Surfaces with Vibrationally Coupled Electron Transfer

Accurate description of nonadiabatic dynamics of molecules at metal surfaces involving electron transfer has been a long-standing challenge for theory. Here, we tackle this problem by first constructing high-dimensional neural network diabatic potentials including state crossings determined by constrained density functional theory, then applying mixed quantum-classical surface hopping simulations to evolve coupled electron-nuclear motion. Our approach accurately describes the nonadiabatic effects in CO scattering from Au(111) without empirical parameters and yields results agreeing well with experiments under various conditions for this benchmark system. We find that both adiabatic and nonadiabatic energy loss channels have important contributions to the vibrational relaxation of highly vibrationally excited CO(v_{i}=17), whereas relaxation of low vibrationally excited states of CO(v_{i}=2) is weak and dominated by nonadiabatic energy loss. The presented approach paves the way for accurate first-principles simulations of electron transfer mediated nonadiabatic dynamics at metal surfaces.

Quantum rates in dissipative systems with spatially varying friction

We investigate whether making the friction spatially dependent on the reaction coordinate introduces quantum effects into the thermal reaction rates for dissipative reactions. Quantum rates are calculated using the numerically exact multi-configuration time-dependent Hartree method, as well as the approximate ring-polymer molecular dynamics (RPMD), ring-polymer instanton methods, and classical molecular dynamics. By conducting simulations across a wide range of temperatures and friction strengths, we can identify the various regimes that govern the reactive dynamics. At high temperatures, in addition to the spatial-diffusion and energy-diffusion regimes predicted by Kramer's rate theory, a (coherent) tunneling-dominated regime is identified at low friction. At low temperatures, incoherent tunneling dominates most of Kramer's curve, except at very low friction, when coherent tunneling becomes dominant. Unlike in classical mechanics, the bath's influence changes the equilibrium time-independent properties of the system, leading to a complex interplay between spatially dependent friction and nuclear quantum effects even at high temperatures. More specifically, a realistic friction profile can lead to an increase (or decrease) of the quantum (classical) rates with friction within the spatial-diffusion regime, showing that classical and quantum rates display qualitatively different behaviors. Except at very low frictions, we find that RPMD captures most of the quantum effects in the thermal reaction rates.

Viscous Current-Induced Forces

We study the motion (translational, vibrational, and rotational) of a diatomic impurity immersed in an electron liquid and exposed to electronic current. An approach based on the linear response time-dependent density functional theory combined with the Ehrenfest dynamics leads to a system of linear algebraic equations, which account for the competing and counteracting effects of the current-induced force (electron wind) and the electronic friction. We find and emphasize the coupling between the center of mass motion and that of the nuclei relative to each other, the feature due to the mediation of the two-body interaction by the environment. The current-induced forces, by means of the dynamic exchange-correlation (xc) kernel f_{xc}(r,r^{'},ω), include the electronic viscosity contribution. Starting from the ground state at the equilibrium internuclear distance and applying a current pulse, we observe three phases of the motion: (i) acceleration due to the prevalence of the current-induced force, (ii) stabilization upon balancing of the two forces, and (iii) deceleration due to the friction after the end of the pulse. At lower, but still metallic, electron densities, the dynamic xc contribution to the force significantly affects the acceleration (deceleration) at the first (third) phase of the process. For the Cs density (r_{s}≈6  a.u.), this correction amounts up to 40% in the rotation regime.

Vibrational energy transfer in collisions of molecules with metal surfaces

The Born-Oppenheimer approximation (BOA), which serves as the basis for our understanding of chemical bonding, reactivity and dynamics, is routinely violated for vibrationally inelastic scattering of molecules at metal surfaces. The title-field therefore represents a fascinating challenge to our conventional wisdom calling for new concepts that involve explicit electron dynamics occurring in concert with nuclear motion. Here, we review progress made in this field over the last decade, which has witnessed dramatic advances in experimental methods, thereby providing a much more extensive set of diverse observations than has ever before been available. We first review the experimental methods used in this field and then provide a systematic tour of the vast array of observations that are currently available. We show how these observations - taken together and without reference to computational simulations - lead us to a simple and intuitive picture of BOA failure in molecular dynamics at metal surfaces, one where electron transfer between the molecule and the metal plays a preeminent role. We also review recent progress made in the theory of electron transfer mediated BOA failure in molecule-surface interactions, describing the most important methods and their ability to reproduce experimental observation. Finally, we outline future directions for research and important unanswered questions.

Time-dependent electron transfer and energy dissipation in condensed media

We study a moving adsorbate interacting with a metal electrode immersed in a solvent using the time-dependent Newns-Anderson-Schmickler model Hamiltonian. We have adopted a semiclassical trajectory treatment of the adsorbate to discuss the electron and energy transfers that occur between the adsorbate and the electrode. Using Keldysh Green's function scheme, we found a non-adiabatically suppressed electron transfer caused by the motion of the adsorbate and coupling with bath phonons that model the solvent. The energy is thus dissipated into electron-hole pair excitations, which are hindered by interacting with the solvent modes and facilitated by the applied electrode potential. The average energy transfer rate is discussed in terms of the electron friction coefficient and given an analytical expression in the slow-motion limit.

Electronic friction near metal surface: Incorporating nuclear quantum effect with ring polymer molecular dynamics

The molecular dynamics with electronic friction (MDEF) approach can accurately describe nonadiabatic effects at metal surfaces in the weakly nonadiabatic limit. That being said, the MDEF approach treats nuclear motion classically such that the nuclear quantum effects are completely missing in the approach. To address this limitation, we combine Electronic Friction with Ring Polymer Molecular Dynamics (EF-RPMD). In particular, we apply the averaged electronic friction from the metal surface to the centroid mode of the ring polymer. We benchmark our approach against quantum dynamics to show that EF-RPMD can accurately capture zero-point energy as well as transition dynamics. In addition, we show that EF-RPMD can correctly predict the electronic transfer rate near metal surfaces in the tunneling limit as well as the barrier crossing limit. We expect that our approach will be very useful to study nonadiabatic dynamics near metal surfaces when nuclear quantum effects become essential.

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.

Theoretical study of nonadiabatic hydrogen atom scattering dynamics on metal surfaces using the hierarchical equations of motion method

Hydrogen atom scattering on metal surfaces is investigated based on a simplified Newns-Anderson model. Both the nuclear and electronic degrees of freedom are treated quantum mechanically. By partitioning all the surface electronic states as the bath, the hierarchical equations of motion method for the fermionic bath is employed to simulate the scattering dynamics. It is found that, with a reasonable set of parameters, the main features of the recent experimental studies of hydrogen atom scattering on metal surfaces can be reproduced. Vibrational states on the chemisorption state whose energies are close to the incident energy are found to play an important role, and the scattering process is dominated by a single-pass electronic transition forth and back between the diabatic physisorption and chemisorption states. Further study on the effects of the atom-surface coupling strength reveals that, upon increasing the atom-surface coupling strength, the scattering mechanism changes from typical nonadiabatic transitions to dynamics in the electronic friction regime.

Nonadiabatic Potential Energy Surfaces for a Molecule on a Surface as Found by Constrained Complete Active Space Theory

In order to study electron-transfer mediated chemical processes on a metal surface, one requires not one but two potential energy surfaces (one ground state and one excited state) as in Marcus theory. In this letter, we report that a novel, dynamically weighted, state-averaged constrained CASSCF(2,2) (DW-SA-cCASSCF(2,2)) can produce such surfaces for the Anderson impurity model. Both ground and excited state potentials are smooth, they incorporate states with a charge transfer character, and the accuracy of the ground state surface can be verified for some model problems by renormalization group theory. Future development of gradients and nonadiabatic derivative couplings should allow for the study of nonadiabatic dynamics for molecules near metal surfaces.

Classical Quantum Friction at Water-Carbon Interfaces

Friction at water-carbon interfaces remains a major puzzle with theories and simulations unable to explain experimental trends in nanoscale waterflow. A recent theoretical framework─quantum friction (QF)─proposes to resolve these experimental observations by considering nonadiabatic coupling between dielectric fluctuations in water and graphitic surfaces. Here, using a classical model that enables fine-tuning of the solid's dielectric spectrum, we provide evidence from simulations in general support of QF. In particular, as features in the solid's dielectric spectrum begin to overlap with water's librational and Debye modes, we find an increase in friction in line with that proposed by QF. At the microscopic level, we find that this contribution to friction manifests more distinctly in the dynamics of the solid's charge density than that of water. Our findings suggest that experimental signatures of QF may be more pronounced in the solid's response rather than liquid water's.

Quantitative molecular simulations

All-atom simulations can provide molecular-level insights into the dynamics of gas-phase, condensed-phase and surface processes. One important requirement is a sufficiently realistic and detailed description of the underlying intermolecular interactions. The present perspective provides an overview of the present status of quantitative atomistic simulations from colleagues' and our own efforts for gas- and solution-phase processes and for the dynamics on surfaces. Particular attention is paid to direct comparison with experiment. An outlook discusses present challenges and future extensions to bring such dynamics simulations even closer to reality.

Dissipative tunneling rates through the incorporation of first-principles electronic friction in instanton rate theory. II. Benchmarks and applications

In Paper I [Litman et al., J. Chem. Phys. (in press) (2022)], we presented the ring-polymer instanton with explicit friction (RPI-EF) method and showed how it can be connected to the ab initio electronic friction formalism. This framework allows for the calculation of tunneling reaction rates that incorporate the quantum nature of the nuclei and certain types of non-adiabatic effects (NAEs) present in metals. In this paper, we analyze the performance of RPI-EF on model potentials and apply it to realistic systems. For a 1D double-well model, we benchmark the method against numerically exact results obtained from multi-layer multi-configuration time-dependent Hartree calculations. We demonstrate that RPI-EF is accurate for medium and high friction strengths and less accurate for extremely low friction values. We also show quantitatively how the inclusion of NAEs lowers the crossover temperature into the deep tunneling regime, reduces the tunneling rates, and, in certain regimes, steers the quantum dynamics by modifying the tunneling pathways. As a showcase of the efficiency of this method, we present a study of hydrogen and deuterium hopping between neighboring interstitial sites in selected bulk metals. The results show that multidimensional vibrational coupling and nuclear quantum effects have a larger impact than NAEs on the tunneling rates of diffusion in metals. Together with Paper I [Litman et al., J. Chem. Phys. (in press) (2022)], these results advance the calculations of dissipative tunneling rates from first principles.

Dissipative tunneling rates through the incorporation of first-principles electronic friction in instanton rate theory. I. Theory

Reactions involving adsorbates on metallic surfaces and impurities in bulk metals are ubiquitous in a wide range of technological applications. The theoretical modeling of such reactions presents a formidable challenge for theory because nuclear quantum effects (NQEs) can play a prominent role and the coupling of the atomic motion with the electrons in the metal gives rise to important non-adiabatic effects (NAEs) that alter atomic dynamics. In this work, we derive a theoretical framework that captures both NQEs and NAEs and, due to its high efficiency, can be applied to first-principles calculations of reaction rates in high-dimensional realistic systems. More specifically, we develop a method that we coin ring polymer instanton with explicit friction (RPI-EF), starting from the ring polymer instanton formalism applied to a system-bath model. We derive general equations that incorporate the spatial and frequency dependence of the friction tensor and then combine this method with the ab initio electronic friction formalism for the calculation of thermal reaction rates. We show that the connection between RPI-EF and the form of the electronic friction tensor presented in this work does not require any further approximations, and it is expected to be valid as long as the approximations of both underlying theories remain valid.

Quantum Dynamics with Electronic Friction

A theory of electronic friction is developed using the exact factorization of the electronic-nuclear wave function. No assumption is made regarding the electronic bath, which can be made of independent or interacting electrons, and the nuclei are treated quantally. The ensuing equation of motion for the nuclear wave function is a nonlinear Schrödinger equation including a friction term. The resulting friction kernel agrees with a previously derived mixed quantum-classical result by Dou et al., [Phys. Rev. Lett. 119, 046001 (2017)]PRLTAO0031-900710.1103/PhysRevLett.119.046001, except for a pseudomagnetic contribution in the latter that is here removed. More specifically, it is shown that the electron dynamics generally washes out the gauge fields appearing in the adiabatic dynamics. However, these are fully re-established in the typical situation where the electrons respond rapidly on the slow time scale of the nuclear dynamics (Markov limit). Hence, we predict Berry's phase effects to be observable also in the presence of electronic friction. Application to a model vibrational relaxation problem proves that the proposed approach represents a viable way to account for electronic friction in a fully quantum setting for the nuclear dynamics.

Independent Trajectory Mixed Quantum-Classical Approaches Based on the Exact Factorization

Mixed quantum-classical dynamics based on the exact factorization exploits the "derived" electron-nuclear correlation (ENC) term aiming for the description of quantum coherences. The ENC contains interactions between the phase of electronic states and nuclear quantum momenta which depend on the spatial shape of the nuclear density.The original surface hopping based on the exact factorization (SHXF) [\textit{J. Phys. Chem. Lett.} \textbf{2018}, \textit{9}, 1097] exploits frozen Gaussian functions to construct the nuclear density in the ENC term while the phase of electronic states is approximated as a fictitious nuclear momentum change.However, in reality, the width of nuclear wave packets varies in time depending on the shape of potential energy surfaces.In this work, we present a modified SHXF approach and a newly-developed Ehrenfest dynamics based on the exact factorization (EhXF) with time-dependent Gaussian functions and phases by enforcing total energy conservation.We perform numerical tests for various one-dimensional two-state model Hamiltonians.Overall, the time-dependent width of Gaussian functions and the energy conserving phase show a reliable decoherence compared to the original frozen Gaussian-based SHXF and the exact quantum mechanical calculation.Especially, the energy conserving phase is crucial for EhXF to reproduce the correct quantum dynamics.

Fluctuation-induced quantum friction in nanoscale water flows

The flow of water in carbon nanochannels has defied understanding thus far¹, with accumulating experimental evidence for ultra-low friction, exceptionally high water flow rates and curvature-dependent hydrodynamic slippage^2-5. In particular, the mechanism of water-carbon friction remains unknown⁶, with neither current theories⁷ nor classical^8,9 or ab initio molecular dynamics simulations¹⁰ providing satisfactory rationalization for its singular behaviour. Here we develop a quantum theory of the solid-liquid interface, which reveals a new contribution to friction, due to the coupling of charge fluctuations in the liquid to electronic excitations in the solid. We expect that this quantum friction, which is absent in Born-Oppenheimer molecular dynamics, is the dominant friction mechanism for water on carbon-based materials. As a key result, we demonstrate a marked difference in quantum friction between the water-graphene and water-graphite interface, due to the coupling of water Debye collective modes with a thermally excited plasmon specific to graphite. This suggests an explanation for the radius-dependent slippage of water in carbon nanotubes⁴, in terms of the electronic excitations of the nanotubes. Our findings open the way for quantum engineering of hydrodynamic flows through the electronic properties of the confining wall.

Random Force in Molecular Dynamics with Electronic Friction

Originally conceived to describe thermal diffusion, the Langevin equation includes both a frictional drag and a random force, the latter representing thermal fluctuations first seen as Brownian motion. The random force is crucial for the diffusion problem as it explains why friction does not simply bring the system to a standstill. When using the Langevin equation to describe ballistic motion, the importance of the random force is less obvious and it is often omitted, for example, in theoretical treatments of hot ions and atoms interacting with metals. Here, friction results from electronic nonadiabaticity (electronic friction), and the random force arises from thermal electron-hole pairs. We show the consequences of omitting the random force in the dynamics of H-atom scattering from metals. We compare molecular dynamics simulations based on the Langevin equation to experimentally derived energy loss distributions. Despite the fact that the incidence energy is much larger than the thermal energy and the scattering time is only about 25 fs, the energy loss distribution fails to reproduce the experiment if the random force is neglected. Neglecting the random force is an even more severe approximation than freezing the positions of the metal atoms or modelling the lattice vibrations as a generalized Langevin oscillator. This behavior can be understood by considering analytic solutions to the Ornstein-Uhlenbeck process, where a ballistic particle experiencing friction decelerates under the influence of thermal fluctuations.

Determining the Effect of Hot Electron Dissipation on Molecular Scattering Experiments at Metal Surfaces

Nonadiabatic effects that arise from the concerted motion of electrons and atoms at comparable energy and time scales are omnipresent in thermal and light-driven chemistry at metal surfaces. Excited (hot) electrons can measurably affect molecule-metal reactions by contributing to state-dependent reaction probabilities. Vibrational state-to-state scattering of NO on Au(111) has been one of the most studied examples in this regard, providing a testing ground for developing various nonadiabatic theories. This system is often cited as the prime example for the failure of electronic friction theory, a very efficient model accounting for dissipative forces on metal-adsorbed molecules due to the creation of hot electrons in the metal. However, the exact failings compared to experiment and their origin from theory are not established for any system because dynamic properties are affected by many compounding simulation errors of which the quality of nonadiabatic treatment is just one. We use a high-dimensional machine learning representation of electronic structure theory to minimize errors that arise from quantum chemistry. This allows us to perform a comprehensive quantitative analysis of the performance of nonadiabatic molecular dynamics in describing vibrational state-to-state scattering of NO on Au(111) and compare directly to adiabatic results. We find that electronic friction theory accurately predicts elastic and single-quantum energy loss but underestimates multiquantum energy loss and overestimates molecular trapping at high vibrational excitation. Our analysis reveals that multiquantum energy loss can potentially be remedied within friction theory whereas the overestimation of trapping constitutes a genuine breakdown of electronic friction theory. Addressing this overestimation for dynamic processes in catalysis and surface chemistry will likely require more sophisticated theories.

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.

Electronic spin separation induced by nuclear motion near conical intersections

Though the concept of Berry force was proposed thirty years ago, little is known about the practical consequences of this force as far as chemical dynamics are concerned. Here, we report that when molecular dynamics pass near a conical intersection, a massive Berry force can appear as a result of even a small amount of spin-orbit coupling (<10⁻³ eV), and this Berry force can in turn dramatically change pathway selection. In particular, for a simple radical reaction with two outgoing reaction channels, an exact quantum scattering solution in two dimensions shows that the presence of a significant Berry force can sometimes lead to spin selectivity as large as 100%. Thus, this article opens the door for organic chemists to start designing spintronic devices that use nuclear motion and conical intersections (combined with standard spin-orbit coupling) in order to achieve spin selection. Vice versa, for physical chemists, this article also emphasizes that future semiclassical simulations of intersystem crossing (which have heretofore ignored Berry force) should be corrected to account for the spin polarization that inevitably arises when dynamics pass near conical intersections.

Spin polarization is at the basis of quantum information and underlies some natural processes, but many aspects still need to be explored. Here, the authors, by quantum mechanical computations, show that even a weak spin-orbit coupling near a conical intersection can induce large spin selection, with consequences for spin manipulation in photochemical or electrochemical reactions.

A simple approximation to the electron-phonon interaction in population dynamics

The modeling of coupled electron-ion dynamics including a quantum description of the nuclear degrees of freedom has remained a costly and technically difficult practice. The kinetic model for electron-phonon interaction provides an efficient approach to this problem, for systems evolving with low amplitude fluctuations, in a quasi-stationary state. In this work, we propose an extension of the kinetic model to include the effect of coherences, which are absent in the original approach. The new scheme, referred to as Liouville-von Neumann + Kinetic Equation (or LvN + KE), is implemented here in the context of a tight-binding Hamiltonian and employed to model the broadening, caused by the nuclear vibrations, of the electronic absorption bands of an atomic wire. The results, which show close agreement with the predictions given by Fermi's golden rule (FGR), serve as a validation of the methodology. Thereafter, the method is applied to the electron-phonon interaction in transport simulations, adopting to this end the driven Liouville-von Neumann equation to model open quantum boundaries. In this case, the LvN + KE model qualitatively captures the Joule heating effect and Ohm's law. It, however, exhibits numerical discrepancies with respect to the results based on FGR, attributable to the fact that the quasi-stationary state is defined taking into consideration the eigenstates of the closed system rather than those of the open boundary system. The simplicity and numerical efficiency of this approach and its ability to capture the essential physics of the electron-phonon coupling make it an attractive route to first-principles electron-ion dynamics.

Spintronics Meets Nonadiabatic Molecular Dynamics: Geometric Spin Torque and Damping on Dynamical Classical Magnetic Texture due to an Electronic Open Quantum System

We analyze a quantum-classical hybrid system of steadily precessing around the fixed axis slow classical localized magnetic moments (LMMs), forming a head-to-head domain wall, surrounded by fast electrons driven out of equilibrium by LMMs and residing within a metallic wire whose connection to macroscopic reservoirs makes electronic quantum system an open one. The model captures the essence of dynamical noncollinear magnetic textures encountered in spintronics, while making it possible to obtain the exact time-dependent nonequilibrium density matrix of electronic systems and split it into four contributions. The Fermi surface contribution generates dissipative (or dampinglike in spintronics terminology) spin torque on LMMs, as the counterpart of electronic friction in nonadiabatic molecular dynamics (MD). Among two Fermi sea contributions, one generates geometric torque dominating in the adiabatic regime, which remains as the only nonzero contribution in a closed system with disconnected reservoirs. Locally geometric torque can have nondissipative (or fieldlike in spintronics terminology) component, acting as the counterpart of geometric magnetism force in nonadiabatic MD, as well as a much smaller dampinglike component acting as "geometric friction." Such current-independent geometric torque is absent from widely used micromagnetics or atomistic spin dynamics modeling of magnetization dynamics based on the Landau-Lifshitz-Gilbert equation, while previous analyses of how to include our Fermi-surface dampinglike torque have severely underestimated its total magnitude.

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.

Nature of Non-Adiabatic Electron-Ion Forces in Liquid Metals

An accurate description of electron-ion interactions in materials is crucial for our understanding of their equilibrium and nonequilibrium properties. Here we assess the properties of frictional forces experienced by ions in noncrystalline metallic systems, including liquid metals and warm dense plasmas, that arise from electronic excitations driven by the nuclear motion due to the presence of a continuum of low-lying electronic states. To this end, we perform detailed ab initio calculations of the full friction tensor that characterizes the set of friction forces. The non-adiabatic electron-ion interactions introduce hydrodynamic couplings between the ionic degrees of freedom, which are sizable between nearest neighbors. The friction tensor is generally inhomogeneous, anisotropic, and nondiagonal, especially at lower densities.

Cooling molecular electronic junctions by AC current

Electronic current flowing in a molecular electronic junction dissipates significant amounts of energy to vibrational degrees of freedom, straining and rupturing chemical bonds and often quickly destroying the integrity of the molecular device. The infamous mechanical instability of molecular electronic junctions critically limits performance and lifespan and raises questions as to the technological viability of single-molecule electronics. Here, we propose a practical scheme for cooling the molecular vibrational temperature via application of an AC voltage over a large, static operational DC voltage bias. Using nonequilibrium Green's functions, we computed the viscosity and diffusion coefficient experienced by nuclei surrounded by a nonequilibrium "sea" of periodically driven, current-carrying electrons. The effective molecular junction temperature is deduced by balancing the viscosity and diffusion coefficients. Our calculations show the opportunity of achieving in excess of 40% cooling of the molecular junction temperature while maintaining the same average current.

High-Fidelity Potential Energy Surfaces for Gas-Phase and Gas-Surface Scattering Processes from Machine Learning

In this Perspective, we review recent advances in constructing high-fidelity potential energy surfaces (PESs) from discrete ab initio points, using machine learning tools. Such PESs, albeit with substantial initial investments, provide significantly higher efficiency than direct dynamics methods and/or high accuracy at a level that is not affordable by on-the-fly approaches. These PESs not only are a necessity for quantum dynamical studies because of delocalization of wave packets but also enable the study of low-probability and long-time events in (quasi-)classical treatments. Our focus here is on inelastic and reactive scattering processes, which are more challenging than bound systems because of the involvement of continua. Relevant applications and developments for dynamical processes in both the gas phase and at gas-surface interfaces are discussed.

Ion energy-loss characteristics and friction in a free-electron gas at warm dense matter and nonideal dense plasma conditions

We investigate the energy-loss characteristics of an ion in warm dense matter (WDM) and dense plasmas concentrating on the influence of electronic correlations. The basis for our analysis is a recently developed ab initio quantum Monte Carlo- (QMC) based machine learning representation of the static local field correction (LFC) [Dornheim et al., J. Chem. Phys. 151, 194104 (2019)JCPSA60021-960610.1063/1.5123013], which provides an accurate description of the dynamical density response function of the electron gas at the considered parameters. We focus on the polarization-induced stopping power due to free electrons, the friction function, and the straggling rate. In addition, we compute the friction coefficient which constitutes a key quantity for the adequate Langevin dynamics simulation of ions. Considering typical experimental WDM parameters with partially degenerate electrons, we find that the friction coefficient is of the order of γ/ω_{pi}=0.01, where ω_{pi} is the ionic plasma frequency. This analysis is performed by comparing QMC-based data to results from the random-phase approximation (RPA), the Mermin dielectric function, and the Singwi-Tosi-Land-Sjölander (STLS) approximation. It is revealed that the widely used relaxation time approximation (Mermin dielectric function) has severe limitations regarding the description of the energy loss of ions in a correlated partially degenerate electrons gas. Moreover, by comparing QMC-based data with the results obtained using STLS, we find that the ion energy-loss properties are not sensitive to the inaccuracy of the static local field correction (LFC) at large wave numbers, k/k_{F}>2 (with k_{F} being the Fermi wave number), but that a correct description of the static LFC at k/k_{F}≲1.5 is important.

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.

Fundamental mechanisms for molecular energy conversion and chemical reactions at surfaces

The dream of theoretical surface chemistry is to predict the outcome of reactions in order to find the ideal catalyst for a certain application. Having a working ab initio theory in hand would not only enable these predictions but also provide insights into the mechanisms of surface reactions. The development of theoretical models can be assisted by experimental studies providing benchmark data. Though for some reactions a quantitative agreement between experimental observations and theoretical calculations has been achieved, theoretical surface chemistry is in general still far away from gaining predictive power. Here we review recent experimental developments towards the understanding of surface reactions. It is demonstrated how quantum-state resolved scattering experiments on reactive and nonreactive systems can be used to test front-running theoretical approaches. Two challenges for describing dynamics at surfaces are addressed: nonadiabaticity in diatomic molecule surface scattering and the increasing system size when observing and describing the dynamics of polyatomic molecules at surfaces. Finally recent experimental studies on reactive systems are presented. It is shown how elementary steps in a complex surface reaction can be revealed experimentally.

Quantum Trajectory Mean-Field Method for Nonadiabatic Dynamics in Photochemistry

The mixed quantum-classical dynamical approaches have been widely used to study nonadiabatic phenomena in photochemistry and photobiology, in which the time evolutions of the electronic and nuclear subsystems are treated based on quantum and classical mechanics, respectively. The key issue is how to deal with coherence and decoherence during the propagation of the two subsystems, which has been the subject of numerous investigations for a few decades. A brief description on Ehrenfest mean-field and surface-hopping (SH) methods is first provided, and then different algorithms for treatment of quantum decoherence are reviewed in the present paper. More attentions were paid to quantum trajectory mean-field (QTMF) method under the picture of quantum measurements, which is able to overcome the overcoherence problem. Furthermore, the combined QTMF and SH algorithm is proposed in the present work, which takes advantages of the QTMF and SH methods. The potential to extend the applicability of the QTMF method was briefly discussed, such as the generalization to other type of nonadiabatic transitions, the combination with multiscale computational models, and possible improvements on its accuracy and efficiency by using machine-learning techniques.

Nonadiabatic Quantum Molecular Dynamics in Dense Manifolds of Electronic States

Most nonadiabatic molecular dynamics methods require the determination of a basis of adiabatic or diabatic electronic states at every time step, but in dense manifolds of electronic states, such approaches become intractable. A notable exception is Ehrenfest molecular dynamics, which can be implemented without explicit determination of such a basis but suffers from unphysical behavior when propagation on a mean-field potential energy surface (PES) does not accurately reflect the true dynamics on multiple electronic states. Here we introduce the multiple cloning for dense manifolds of states (MCDMS) method, a systematically improvable approximation to the multiple cloning method. MCDMS avoids both the mean-field PES problem and the need to compute the full electronic spectrum. This is achieved by reformulating multiple cloning to use a subspace of approximate eigenstates constructed from the time-dependent Ehrenfest electronic wave function. By application to model systems, we show that this approach allows a substantial reduction in the size of the required electronic basis without significant loss in accuracy.

Dynamics in reactions on metal surfaces: A theoretical perspective

Recent advances in theoretical characterization of reaction dynamics on metal surfaces are reviewed. It is shown that the widely available density functional theory of metals and their interactions with molecules have enabled first principles theoretical models for treating surface reaction dynamics. The new theoretical tools include methods to construct high-dimensional adiabatic potential energy surfaces, to characterize nonadiabatic processes within the electronic friction models, and to describe dynamics both quantum mechanically and classically. Three prototypical surface reactions, namely, dissociative chemisorption, Eley-Rideal reactions, and recombinative desorption, are surveyed with a focus on some representative examples. While principles governing gas phase reaction dynamics may still be applicable, the presence of the surface introduces a higher level of complexity due to strong interaction between the molecular species and metal substrate. Furthermore, most of these reactive processes are impacted by energy exchange with surface phonons and/or electron-hole pair excitations. These theoretical studies help to interpret and rationalize experimental observations and, in some cases, guide experimental explorations. Knowledge acquired in these fundamental studies is expected to impact many practical problems in a wide range of interfacial processes.

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.

Hot-electron effects during reactive scattering of H2 from Ag(111): the interplay between mode-specific electronic friction and the potential energy landscape

The breakdown of the Born-Oppenheimer approximation gives rise to nonadiabatic effects in gas-surface reactions at metal surfaces. However, for a given reaction, it remains unclear which factors quantitatively determine whether these effects measurably contribute to surface reactivity in catalysis and photo/electrochemistry. Here, we systematically investigate hot electron effects during H2 scattering from Ag(111) using electronic friction theory. We combine first-principles calculations of tensorial friction by time-dependent perturbation theory based on density functional theory and an analytical neural network representation, to overcome the limitations of existing approximations and explicitly simulate mode-specific nonadiabatic energy loss during molecular dynamics. Despite sizable hot-electron-induced energy loss, no measurable nonadiabatic effects can be found for H2 scattering on Ag(111). This is in stark contrast to previous reports for vibrationally excited H2 scattering on Cu(111). By detailed analysis of the two systems, we attribute this discrepancy to a subtle interplay between the magnitude of electronic friction along intramolecular vibration and the shape of the potential energy landscape that controls the molecular velocity at impact. On the basis of this characterization, we offer guidance for the search of highly nonadiabatic surface reactions.

Hot electron effects during reactive scattering of H2 from Ag(111): assessing the sensitivity to initial conditions, coupling magnitude, and electronic temperature

We use an analytical representation of electronic friction for H₂ on Ag(111) to assess the validity and robustness of the MDEF method based on TDPT.