where N is the total number of trajectories. Instead of calculating the polaritonic state density matrix [ρ pl (R l (t))] using Eq. 36 directly, in the estimator used in this work, we calculate the diagonal elements of ρ pl using the active state index, and calculate the offdiagonal elements using the polaritonic state expansion coefficients {c I (t)}