Communication: The correct interpretation of surface hopping trajectories: how to calculate electronic properties

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.