“…In such systems, the direct coupling between gain molecules is maximal, which leads to random Coulomb shifts of molecule excitation energies and, hence, to dephasing [59,60]. For small systems, the spasing eigenstates were found via exact numerical diagonalization [58]; however, for larger systems with realistic random dipole orientations, direct numerical determination of many-body eigenstates of interacting QEs is not feasible. At the same time, for large ensembles of randomlyoriented QEs, the ensemble-averaged dipole coupling between individual QEs vanishes, while weak fluctuations of QEs' excitation energies do not significantly affect the collective system eigenstates strongly coupled to radiation [61].…”