where bα,n , b † α,n are the α th bath phonon annihilation and creation operators for the n th molecule, who satisfy the bosonic commutation relations, with ω α the phonon frequency. The last term ĤSB in Eq. 1 characterizes the system-bath interaction which describe the excitonphonon interaction, which we assume takes the linear form as follows