Ehrenfest+R approach that uses the modified γ ′ R . This modified Ehrenfest+R approach produces results identical to those obtained from the L-MFE method for this model. Consequentially, in Fig. , the population dynamics of the modified Ehrenfest+R approach are identical to the Lindblad population dynamics. The corresponding coherence dynamics are provided in Fig. , which again match perfectly with the Lindblad coherence dynamics. This demonstrates that it is important to include the effect that modifying the magnitudes of the coefficients has on the coherences between the states (see analysis in Appendix B). We note that using the modified γ ′ R in the Ehrenfest+R method yields the same expectation values for the populations and coherences as the L-MFE method as Γdt → 0, but begins to lose accuracy for larger Γdt due to the use of a Poisson process for decoherence (through Eq. A28), which is sensitive to the choice of Γdt (Eq. B7).