We note that because ĤQ = 0 in Models 1 & 2, the L-MFE method in Eq. 72 becomes c(t + dt) = T (dt) • c(t), where the dynamics are completely dictated by the Lindblad decay process governed by L L. The dynamics are insensitive to the choice of dt, which means that one can choose an arbitrarily large dt and obtain identical results. This is not the case for Ehrenfest+R approach, where the "+R" dynamics will be sensitive to the choice of dt even when ĤQ = 0.