Evaluation of tight-binding DFT performance for the description of organic photochromes properties

Photochromic molecules are widely studied and developed for their many potential applications. To optimize the required properties through theoretical models, a considerable chemical space is to be explored, and their environment in devices is to be accounted for.. To this end, cheap and reliable computational methods can be powerful tools to steer synthetic developments. As ab initio methods remain costly for extensive studies (in terms of the size of the system and/or number of molecules), semiempirical methods such as density functional tight-binding (TB) could offer a good compromise between accuracy computational cost. However, these approaches necessitate benchmarking on the families of compounds of interest. Thus, the aim of the present study is to evaluate the accuracy of several key features calculated with TB methods (DFTB2, DFTB3, GFN2-xTB, and LC-DFTB2) for three sets of photochromic organic molecules: azobenzene (AZO), norbornadiene/quadricyclane (NBD/QC), and dithienylethene (DTE) derivatives. The features considered here are the optimized geometries, the difference in energy between the two isomers (ΔE), and of the energies of the first relevant excited states. All the TB results are compared to those obtained with DFT methods and state-of-the-art electronic structure calculation methods: DLPNO-CCSD(T) for ground states and DLPNO-STEOM-CCSD for excited states. Our results show that, overall, DFTB3 is the TB method leading to the best results for the geometries and the ΔE values and can be used alone for these purposes for NBD/QC and DTE derivatives. Single point calculations at the r2SCAN-3c level using TB geometries allow circumventing the deficiencies of the TB methods in the AZO series. For electronic transition calculations, the range separated LC-DFTB2 method is the most accurate TB method tested for AZO and NBD/QC derivatives, in close agreement with the reference.