In this work, we examine the importance of nuclear quantum effects on capturing the line broadening and vibronic structure of optical spectra. We determine the absorption spectra of three aromatic molecules indole, pyridine, and benzene using time dependent density functional theory with several molecular dynamics sampling protocols: force-field based empirical potentials, ab initio simulations, and finally path-integrals for the inclusion of nuclear quantum effects. We show that the absorption spectrum for all these chromophores are similarly broadened in the presence of nuclear quantum effects regardless of the presence of hydrogen bond donor or acceptor groups. We also show that simulations incorporating nuclear quantum effects are able to reproduce the heterogeneous broadening of the absorption spectra even with empirical force fields. The spectral broadening associated with nuclear quantum effects can be accounted for by the broadened distribution of chromophore size as revealed by a particle in the box model. We also highlight the role that nuclear quantum effects have on the underlying electronic structure of aromatic molecules as probed by various electrostatic properties.