Search papers, labs, and topics across Lattice.
This paper addresses the inconsistency between stochastic surface hops and wavefunction propagation in Fewest-Switches Surface Hopping (FSSH) by developing a fully-consistent hopping-probability formula specifically for trivial crossings. This novel approach ensures that attempted hops accurately reflect wavefunction time-evolution across arbitrary time steps, mitigating errors and enabling larger integration steps. The method is validated on challenging systems with trivial crossings, demonstrating high accuracy even with significantly increased computational efficiency.
Trivial crossings, a notorious FSSH bottleneck, are finally tamed by a hopping scheme that guarantees wavefunction consistency across any time step.
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.