being identical for all excitons (for n ∈ [1, N ]). We use a Drude-Lorentz form for the spectral density in our investigations, where ω f is the bath characteristic frequency, and λ is the reorganization energy, which can also be expressed in terms of the coupling strength and the phonon frequencies as