To get the spectra, we have to calculate the time evolution of transition dipole moment on the excited state of aggregate μ(t) = Û (t) † μ Û (t) and Û (t) = exp(-i Ĥagg t). Here, we adopt the NMSSE proposed in our group to propagate. Within the NMSSE framework, the correlation function for absorption spectra is given by