The dynamics governed by L Eh (Eq. A17) do not match the Lindblad dynamics governed by L (Eq. A5) since the decay rates k Eh and γ Eh from Eqs. A18 and A19 are not constant and depend on the density matrix elements, among other reasons. This is not surprising because it is impossible to recover open system Lindblad dynamics through the deterministic dynamics of a Hamiltonian (Eq. A10). Thus, an additional relaxation propagation (the "+R" stage of the dynamics) is introduced, which is designed to correct both the diagonals and offdiagonals of the reduced density matrix elements in the {|0⟩, |1⟩} subspace. This additional decay process (the +R dynamics) is governed by