Quantifying the Performances of DFT for Predicting Vibrationally Resolved Optical Spectra: Asymmetric Fluoroborate Dyes as Working Examples

This article aims at a quantitative assessment of the performances of a panel of exchange-correlation functionals, including semilocal (BLYP and PBE), global hybrids (B3LYP, PBE0, M06, BHandHLYP, M06-2X, and M06-HF), and range-separated hybrids (CAM-B3LYP, LC-omega PBE, LC-BLYP, omega B97X, and omega B97X-D), in predicting the vibrationally resolved absorption spectra of BF2-carrying compounds. To this end, for 19 difluoroborates as examples, we use, as a metric, the vibrational reorganization energy (lambda(vib)) that can be determined based on the computationally efficient linear coupling model (a.k.a. vertical gradient method). The reference values of were determined by employing the CC2 method combined with the cc-pVTZ basis set for a representative subset of molecules. To validate the performances of CC2, comparisons with experimental data have been carried out as well. This study shows that the vibrational reorganization energy, involving Huang-Rhys factors and normal-mode frequencies, can indeed be used to quantify the reliability of functionals in the calculations of the vibrational fine structure of absorption bands, i.e., an accurate prediction of the vibrational reorganization energy leads to absorption band shapes better fitting the selected reference. The CAM-B3LYP, M06-2X, wB97X-D, wB97X, and BHandHLYP functionals all deliver vibrational reorganization energies with absolute relative errors smaller than 20% compared to CC2, whereas 10% accuracy can be achieved with the first three functionals. Indeed, the set of examined exchange-correlation functionals can be divided into three groups: (i) BLYP, B3LYP, PBE, PBE0, and M06 yield inaccurate band shapes (lambda(vib,TDDFT) lambda(vib,CC2)). This study also demonstrates that lambda(vib), can be reliably estimated using the CC2 model and the relatively small cc-pVDZ basis set. Therefore, the linear coupling model combined with the CC2/cc-pVDZ level of theory can be used as a very efficient approach to determine lambda(vib) values that can be used to select the most adequate functional for more accurate vibronic calculations, e.g., including more refined models and environmental effects