Non-Born-Oppenheimer molecular dynamics of NaFH photodissociation

The accuracy of non-Born-Oppenheimer (electronically nonadiabatic) semiclassical trajectory methods for simulations of "deep quantum" systems is reevaluated in light of recent quantum mechanical calculations of the photodissociation of the NaFH van der Waals complex. In contrast to the conclusion arrived at in an earlier study, semiclassical trajectory methods are shown to be qualitatively accurate for this system, thus further validating their use for systems with large electronic energy gaps. Product branching in semiclassical surface hopping and decay-of-mixing calculations is affected by a region of coupling where the excited state is energetically forbidden. Frustrated hops in this region may be attributed to a failure of the treatment of decoherence, and a stochastic model for decoherence is introduced into the surface hopping method and is shown to improve the agreement with the quantum mechanical results. A modification of the decay-of-mixing method resulting in faster decoherence in this region is shown to give similarly improved results.