Frustrated hops in Ring Polymer Surface Hopping: Real-time dynamics and detailed balance

Abstract

Ring Polymer Surface-Hopping (RPSH) has been recently introduced as a well-tailored method for incorporating nuclear quantum effects (NQEs), such as zero-point energy and tunneling, into non-adiabatic molecular dynamics simulations. The practical widespread usage of RPSH demands a comprehensive benchmarking of different reaction regimes and conditions with equal emphasis on demonstrating both the cons and pros of the method. Here, we investigate the fundamental questions related to the conservation of energy and detailed balance in the context of RPSH. Using Tully's avoided crossing model as well as a 2-level system coupled to a classical bath undergoing Langevin dynamics, we probe the critical problem of the proper treatment of the classically forbidden transitions stemming from the surface hopping algorithm. We show that proper treatment of these frustrated hops is key to the accurate description of real-time dynamics as well as reproducing the exact quantum Boltzmann population.

arXiv comment8 pages, 6 Figures, 1 Table

References in corpus (7)

• Efficient stochastic thermostatting of path integral molecular dynamics
• How to remove the spurious resonances from ring polymer molecular dynamics
• Mapping Variable Ring Polymer Molecular Dynamics: A Path-Integral Based Method for Nonadiabatic Processes
• Direct simulation of electron transfer using ring polymer molecular dynamics: Comparison with semiclassical instanton theory and exact quantum methods
• Path-Integral Isomorphic Hamiltonian for Including Nuclear Quantum Effects in Non-adiabatic Dynamics
• Coherent State Mapping Ring-Polymer Molecular Dynamics for Non-Adiabatic quantum propagations
• State Dependent Ring Polymer Molecular Dynamics for Investigating Excited Nonadiabatic Dynamics