Hidden optical nonlinearities in linear spectra of quantum emitter arrays

Classical optical frameworks such as the discrete dipole approximation (DDA) assume that the linear spectrum of coupled quantum emitters can be computed solely from the linear susceptibilities of individual constituents. However, recent polariton studies show that cavity linear response can encode nonlinear optical susceptibilities. Here, we demonstrate that this phenomenon is more general: emitter-emitter interactions allow nonlinearities of individual emitters to emerge in the linear response of arrays, without cavities or permutational symmetry. To illustrate this phenomenon, we show linear spectra for coupled heterodimers and linear chains, and demonstrate that Raman features of individual monomers show up as vibrational sidebands of collective resonances. Moreover, tuning Raman-type anharmonicities enables systematic control of spectral features, establishing a genuine quantum optical effect in molecular aggregates and quantum emitter arrays, which goes beyond mean-field descriptions in light-matter interactions.