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.