Communication: Global flux surface hopping in Liouville space

Recent years have witnessed substantial progress in the surface hopping (SH) formulation of non-adiabatic molecular dynamics. A generalization of the traditional fewest switches SH (FSSH), global flux SH (GFSH) utilizes the gross population flow between states to derive SH probabilities. The Liouville space formulation of FSSH puts state populations and coherences on equal footing, by shifting the hopping dynamics from Hilbert to Liouville space. Both ideas have shown superior results relative to the standard FSSH in Hilbert space, which has been the most popular approach over the past two and a half decades. By merging the two ideas, we develop GFSH in Liouville space. The new method is nearly as straightforward as the standard FSSH, and carries comparable computational expense. Tested with a representative super-exchange model, it gives the best performance among all existing techniques in the FSSH series. The obtained numerical results match almost perfectly the exact quantum mechanical solutions. Moreover, the results are nearly invariant under the choice of a basis state representation for SH, in contrast to the earlier techniques which exhibit notable basis set dependence. Unique to the developed approach, this property is particularly encouraging, because exact quantum dynamics is representation independent. GFSH in Liouville space significantly improves accuracy and applicability of SH for a broad range of chemical and physical processes.