For the last case of Lorentzian distribution for N = 1 in Fig. , the UP and LP collapse is nearly linear in the energy disorder magnitude while both the Gaussian-and rectangulardisorder cases were non-linear. Specifically, we are tracking the maximal peak location in the spectra as a function of the disorder. Further, the collapse is much more rapid compared to the other two cases, with LP and UP merging around 2 ∼ 4 σ / Ã0 compared to 16 + σ / Ã0 for rectangular and 10 σ / Ã0 for Gaussian. This could be due to the intrinsic long-tailed nature of the Lorentzian distribution compared to the short-tail Gaussian and no-tail rectangular distributions. This implies that for few-molecule spectroscopy, for example, in plasmonic cavity designs , the emergence of UP and LP spectral features may depend heavily on the shape and magnitude of the molecular disorder. These features of the molecular ensemble are dictated by the internal processes of the molecules as well as by their local environment.