Simulations of spectroscopies for large-scale systems are challenging in several ways that energy calculations are not, at least if they are to be made practical as an analytical tool for exploring ideas and interpreting observations. In time-dependent DFT, the canonical computational task is the construction of Fock matrices and therefore, the calculation of response properties scales identically as the SCF optimization. The key difference, however, is that in the former case there is typically a need to determine about two orders of magnitude more Fock matrices as compared to the latter case. To put an emphasis on the needed responsiveness of the spectrum simulation toolbox, we have introduced a KPI in terms of the number of Fock matrices produced per hour, and shown that by means of GPU acceleration, we can reach an indicator value above 26 h -1 for a system with close to 12,000 basis functions and with employment of a range-separated functional that has a higher computational cost for ERIs as compared to the corresponding hybrid GGA-functional.