where ρ(0) is the initial state for the exciton-photonic subsystem (the quantum subsystem), and ρR is the thermal initial condition for the phonon, Q = Tr R [e -β ĤR ] is the phonon bath partition function, and β = 1/k B T with T as the temperature and k B as the Boltzmann factor.