Fig. presents the population dynamics of Model 1 (Eq. 76a), obtained from the L-MFE approach (dots) as well as exact Lindblad dynamics (solid lines). In this case, the only dynamics present are the Lindblad exponential decay from state |1⟩ to state |0⟩. In Fig. , the diagonal populations of both states |1⟩ and |0⟩ exhibit the expected exponential decay/growth at a rate of Γ = 0.05, while the coherence between the two states stays at 0. This is expected because, in the absence of any Hamiltonian dynamics, the coherences between states |1⟩ and |0⟩ should only monotonically decrease to 0 from the initial time coherence ρ 01 (0). Fig. presents the population dynamics of Model 2 (Eq. 76b), which has an initial condition of a superposition of state |1⟩ and |0⟩. Similar to Model 1, only Lindblad exponential decay/growth is present. The diagonal populations match the Lindblad dynamics, and the coherence shows the expected Γ/2 decay rate from Lindblad dynamics.