Linear and angular momentum conservation in surface hopping methods

We demonstrate that, for systems with spin-orbit coupling and an odd number of electrons, the standard fewest switches surface hopping algorithm does not conserve the total linear or angular momentum. This lack of conservation arises not so much from the hopping direction (which is easily adjusted) but more generally from propagating adiabatic dynamics along surfaces that are not time reversible. We show that one solution to this problem is to run along eigenvalues of phase-space electronic Hamiltonians H(R, P) (i.e., electronic Hamiltonians that depend on both nuclear position and momentum) with an electronic-nuclear coupling Gamma center dot P [see Eq. (25)], and we delineate the conditions that must be satisfied by the operator Gamma. The present results should be extremely useful as far as developing new semiclassical approaches that can treat systems where the nuclear, electronic orbital, and electronic spin degrees of freedom altogether are all coupled together, hopefully including systems displaying the chiral-induced spin selectivity effect.