On the other hand, our previous theoretical work suggests that if the relaxation are dictated by a "golden rule" type of law, then even the dynamics are largely dictated by collective quantity √ NA 0 . In the future, we will extend this exploration in three directions: (i) non-linear response via a modification of the observable in Eq. 15, (ii) the inclusion of multi-mode cavity structures and (iii) dynamical effects through Chebyshev propagation , all with varied molecular distributions and number of collectively coupled molecules. Further, the inclusion of the dipole self-energy (DSE), while small, may yield non-trivial results in the large-N limit due to the ∝ N scaling of the DSE term compared to the ∝ √ N scaling of the bilinear coupling term.