We address the performance of the vertical and adiabatic Franck-Condon (VFC/AFC) approaches combined with time-independent or time-dependent (TI/TD) formalisms in simulating the one-photon absorption spectra of three flavin compounds with distinct structural features. Calculations were done in the gas phase and in two solvents (water, benzene) for which experimental reference measurements are available. We utilized the independent mode displaced harmonic oscillator model without or with frequency alteration (IMDHO/IMDHO-FA) and also accounted for Duschinsky mixing effects. In the initial validation on the first excited singlet state of riboflavin, the range-separated functionals, CAM-B3LYP and ωB97xD, showed the best performance, but B3LYP also gave a good compromise between peak positions and spectral topology. Large basis sets were not mandatory to obtain high-quality spectra for the selected systems. The presence of a symmetry plane facilitated the computation of vibrationally broadened spectra, since different FC variants yield similar results and the harmonic approximation holds rather well. Compared with the AFC approach, the VFC approach performed equally well or even better for all three flavins while offering several advantages, such as avoiding error-prone geometry optimization procedures on excited-state surfaces. We also explored the advantages of curvilinear displacements and of a Duschinsky treatment for the AFC spectra in cases when a rotatable group is present on the chromophore. Taken together, our findings indicate that the combination of the VFC approach with the TD formalism and the IMDHO-FA model offers the best overall performance.
Read full abstract