Ensuring Complete Consistency between Hops and Wavefunction Propagation in Fewest-Switches Surface Hopping

Fewest-switches surface hopping (FSSH) is a popular and versatile approach for performing dynamical simulations beyond the Born-Oppenheimer approximation, in which nonadiabatic effects are incorporated through stochastic hops of the active potential energy surface. However, the FSSH algorithm does not ensure that hops are completely consistent with the underlying electronic wavefunction dynamics. This can lead to significant errors as well as requiring small propagation time steps, which can significantly compromise the computational efficiency of the method. These problems are known to be particularly severe at so-called trivial crossings, for which the nonadiabatic coupling is entirely localized at the crossing point. Here, we solve the trivial-crossing problem by developing a fully-consistent hopping-probability formula, which generates attempted hops that match the time-evolution of the wavefunction across any size time step by construction. A state-tracking algorithm is also utilized to ensure that attempted hops are always accepted at trivial crossings. Through application to systems that contain extremely challenging trivial crossings, we demonstrate that this hopping scheme is capable of obtaining highly accurate results, even when large propagation time steps are used.

Source: arXiv:2610.18743v1 - http://arxiv.org/abs/2610.18743v1 PDF: https://arxiv.org/pdf/2610.18743v1 Original Link: http://arxiv.org/abs/2610.18743v1