It is well known 151 that FSSH suffers from producing overly coherent results (or, equivalently, a lack of proper decoherence between quantum states) within the expansion electronic coefficients and will subsequently be problematic for the polaritonic coefficients. Many ad hoc corrections exist to modify the expansion coefficients in FSSH to account for decoherence, such as the instantaneous decoherence correction (IDC), 155;127 the energy-based decoherence correction (EDC), 156 etc., as well as other forms of the surface hopping scheme, such as the augmented surface hopping (A-FSSH), 157 the decoherence-induced surface hopping (DISH), and the global flux surface hopping 159 schemes. More recently, a spin-mapping 145;144;140;137;138;139;160;141;142;143 analogue of the surface hopping approach was constructed and mitigates much of the shortcomings of the standard FSSH procedure without ad hoc adjustments. A major simplicity afforded by the FSSH method is that the derivative coupling vectors d j k (R(t ) ) are not explicitly required as the nuclear forces (unlike in the mean-field MQC methods) do not require this quantity for time-evolution of the electronic or nuclear DOFs except at the hops for the nuclear velocity rescaling. To be clear, the electronic propagation only requires the scalar non-adiabatic coupling terms