where ω k are the phonon frequencies, and ν † k,n and νk,n are the k th phonon mode's creation and annihilation operators for the n th molecule that satisfy the bosonic commutation relations. The last term ĤSB in Eq. 1. characterizes the system-bath interaction, which we assume takes the linear form,