In particular, the oscillation in Im[ρ 20 ] is due to the presence of electronic coupling between state |0⟩ and |2⟩, without further decoherence from L (see Eq. 24), suggesting that the L-MFE method can correctly describe dynamics involving {|j⟩} states even in the presence of random phases in the coefficient of state |0⟩.