where ρ(t) is the full density matrix of the system, and the trace Tr x and Tr X are performed on the phonon bath {x j } and the effective photon bath { Xj }, respectively (see Eq. 18). In order to evaluate the forward rate constant, we follow the previous work by defining the timedependent reactant (R) and product state (P) populations as