where ρ is the density matrix describing the state of the closed system. To model the dimer system, we use Ĥdimer as the Hamiltonian in Eq. 16. The expression above leads to a system of ODEs for the matrix elements of the density matrix, which can be solved to obtain the population dynamics of the exciton (i.e., the diagonal elements of the density matrix ρ).