Surface hopping, transition state theory, and decoherence. II. Thermal rate constants and detailed balance

We investigate a simple approach to compute a non-adiabatic thermal rate constant using the fewest switches surface hopping (FSSH) dynamics. We study the effects of both decoherence (using our augmented-FSSH (A-FSSH) algorithm) and forbidden hops over a large range of parameters, including high and low friction regimes, and weak and strong electronic coupling regimes. Furthermore, when possible, we benchmark our results against exact hierarchy equations of motion results, where we usually find a maximum error of roughly a factor of two (at reasonably large temperatures). In agreement with Hammes-Schiffer and Tully, we find that a merger of transition state theory and surface hopping can be both accurate and efficient when performed correctly. We further show that detailed balance is followed approximately by A-FSSH dynamics.