To account for the influence of the solvent and the continuum of far-field electromagnetic modes outside the cavity, we consider an environment composed of two bosonic baths, H E = H ph + H pt , namely, a phonon bath for dissipative solvent DoFs and a photon bath for the cavity lossy environment. Each bath can be modeled as an ensemble of infinite harmonic oscillators,