where ω α are the frequencies for the α th phonon mode, b † α,n and bα,n are the bath phonon creation and annihilation operators that satisfy the bosonic commutation relations. Ĥsb describes the systembath interaction, where c α denotes the coupling strength between the molecules and the α-th bath phonon mode. The system-bath interactions are determined by the spectral density