where Λ (defined in Eq. 3) characterizes the couplings between N RPV modes (e.g., solvent DOF) {Q j } with R 0 , and ω Q is the solvent frequency which we assume to be identical for all {Q j }. In addition, Γ Q characterizes the excitation decay rate in the {Q j } modes