Search

In part I, we presented the ring-polymer instanton with explicit friction (RPI-EF) method and showed how it can be connected to the \textit{ab initio} electronic friction formalism. This framework allows 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 second part, 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 cross-over 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 part I, these results advance the calculations of dissipative tunneling rates from first principles.

Reactions involving adsorbates on metallic surfaces and impurities in bulk metals are ubiquitous in a wide range of technological applications. The theoretical modelling 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. In more detail, 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 \textit{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.