For clarity, we denote the reduced density matrix in the adiabatic-Fock basis as ρ af ij (t), and ρ pl IJ (t) is the reduced density matrix in the polariton basis (expressed in Eq. 36). To obtain the density matrix in the polaritonic basis during the dynamics, we use the following unitary transformation