In a recent paper, we presented a road map for how Tully's fewest switches surface hopping (FSSH) algorithm can be derived, under certain circumstances, from the mixed quantum-classical Liouville equation. In this communication, we now demonstrate how this new interpretation of surface hopping can yield significantly enhanced results for electronic properties in nonadiabatic calculations. Specifically, we calculate diabatic populations for the spin-boson problem using FSSH trajectories. We show that, for some Hamiltonians, without changing the FSSH algorithm at all but rather simply reinterpreting the ensemble of surface hopping trajectories, we recover excellent results and remove any and all ambiguity about the initial condition problem.
We compare the dynamics of Fewest Switches Surface Hopping (FSSH) in different parameter regimes of the spin-boson model. We show that for exceptional regions of the spin-boson parameter space, FSSH dynamics are in fact time-reversible. In these rare instances, FSSH does recover the correct Marcus rate scaling (as a function of diabatic coupling) without the addition of decoherence. In regions where dynamics are irreversible, however, FSSH does not recover the correct Marcus rate scaling. Finally, by comparing the friction dependence of rates predicted by various decoherence schemes to an analytic result by Zusman, we provide yet more evidence that the method of introducing decoherence has a qualitative effect on the accuracy of results and this effect must be treated carefully.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.