  Mean-Field Ring Polymer Molecular Dynamics (MF-RPMD) is a powerful,
efficient, and accurate method for approximate quantum dynamic simulations of
multi-level system dynamics. Initial efforts to compute nonadiabatic reaction
rates using MF-RPMD were not successful; recent work showed that this can be
remedied by including a simple, if adhoc, correction term that accounts for the
formation of `kinked' or mixed electronic state ring polymer configurations.
Here, we build on this idea, introducing a electronic state population based
reaction coordinate and novel skew dividing surface that constrains nuclear
positions to configurations where the reactant and product state potentials are
near-degenerate and that samples kinked electronic state configurations. We
then demonstrate the numerical accuracy of this method in computing rates for a
series of nonadiabatic model systems.
