Efficient simulation of resonance Raman spectra with tight-binding approximations to density functional theory

Resonance Raman spectroscopy has long been established as one of the most sensitive techniques for detection, structure characterization, and probing the excited-state dynamics of biochemical systems. However, the analysis of resonance Raman spectra is much facilitated when measurements are accompanied by Density Functional Theory (DFT) calculations that are expensive for large biomolecules. In this work, resonance Raman spectra are therefore computed with the Density Functional Tight-Binding (DFTB) method in the time-dependent excited-state gradient approximation. To test the accuracy of the tight-binding approximations, this method is first applied to typical resonance Raman benchmark molecules, such as β-carotene, and compared to results obtained with pure and range-separated exchange-correlation functionals. We then demonstrate the efficiency of the approach by considering a computationally challenging heme variation. Overall, we find that the vibrational frequencies and excited-state properties (energies and gradients) that are needed to simulate the spectra are reasonably accurate and suitable for interpretation of experiments. We can therefore recommend DFTB as a fast computational method to interpret resonance Raman spectra.