where The initial conditions are similar to that of the EH approach; however, if there exists a distribution of polaritonic coefficients at initial time {c α (0)}, then initial active state should also be sampled independently for each trajectory (similarly to sampling of nuclear DOFs) from the probability distribution defined by P S (0) = {|c a (0)| 2 }. It is well known, 262 FSSH suffers from producing overly coherent (or lack of proper electronic decoherence) 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), the energy-based decoherence correction (EDC), 265 etc., as well as other forms of the surface hopping scheme, such as the augmented surface hopping (A-FSSH), 266 the decoherence-induced surface hopping (DISH), 267 and the global flux surface hopping 268 schemes.