To simulate disorder effects in linear spectra (Figs. ), we fix the collective coupling strength 2 √ N ℏg c at 0.2 eV for N = 500 molecules. For the molecules, we used the average excitation energy of ε = 2.0 eV. Throughout this work, we consider the resonance condition where