Fast vs Accurate: Evaluating TD-DFTB for Large-Scale Screening of Organic Dyes
Abstract
The effectiveness of the semi-empirical TD-DFTB (Time-Dependent Density Functional Tight Binding) method in reproducing the spectral properties of organic dyes was investigated using the example of a library of isomers of the thiophene-containing donor-acceptor Effenberger dye, known for its pronounced solvatochromism. The aim of the work was to find out how suitable the accelerated TD-DFTB approach is for modern molecular design tasks, with a necessity to quickly and reliably identify compounds with intense electronic transitions in the long-wavelength region of the UV-Vis spectrum. The library contained 60 structures in which the positions of the donor (N,N-dimethylamine) and acceptor (NO₂) substituents, as well as the degree of planarity of the π-framework, were systematically varied. For each isomer, the geometry was first optimized at the DFTB level, after which the excitation energies were calculated using the TD-DFTB method. The obtained values were compared with TD-DFT calculations (B3LYP and CAM-B3LYP functionals) performed with geometries, obtained both at DFT and DFTB methods. Such a hybrid scheme significantly reduces the computational costs, allowing screening of large libraries without losing the accuracy. The correlation between excitation energies calculated by TD-DFTB and TD-DFT is given. As obtained, TD-DFTB tends to systematically underestimate the excitation energies, but largely reflects compounds with minimal excitation energies and large oscillator strengths, which makes it a reliable tool for initial screening. Several isomers with long-wavelength absorption and sufficient transition intensity were identified, which are promising for further modification. Thus, TD-DFTB in combination with TD-DFT on optimized DFTB geometries demonstrates an optimal balance between accuracy and speed for prescreening donor-acceptor dyes with given spectral parameters, which significantly enhances the capabilities of rational design of functional organic materials.
Downloads
References
Martinez-Mayorga K., Madariaga-Mazon A., Medina-Franco J. L., Maggiora G. The impact of chemoinformatics on drug discovery in the pharmaceutical industry. Expert Opin. Drug Discov. 2020, 15(3), 293-306.
Humbeck L., Koch O. What can we learn from bioactivity data? Chemoinformatics tools and applications in chemical biology research. ACS Chem. Biol. 2017, 12(1), 23-35.
Runge E., Gross E. K. Density-functional theory for time-dependent systems. Phys. Rev. Lett. 1984, 52(12), 997.
Rüger R., Van Lenthe E., Heine T., Visscher L. Tight-binding approximations to time-dependent density functional theory—A fast approach for the calculation of electronically excited states. J. Chem. Phys. 2016, 144(18).
Effenberger F., Wuerthner F., Steybe F. Synthesis and solvatochromic properties of donor-acceptor-substituted oligothiophenes. J. Org. Chem. 1995, 60(7), 2082-2091.
https://dftbplus-recipes.readthedocs.io/en/latest/
Zakharov A. B., Tsarenko D. K., Ivanov V. V. Topological characteristics of iterated line graphs in the QSAR problem: A multigraph in the description of properties of unsaturated hydrocarbons. Struct. Chem. 2021, 32(4), 1629-1639.
Zakharov A. B., Ivanov V. V. A new approach in topological descriptors usage. Iterated line graphs in the theoretical prediction of physico-chemical properties of saturated hydrocarbons. Kharkiv Univ. Bull. Chem. Ser. 2019, (32), 38-45.
Berdnyk M. I., Zakharov A. B., Ivanov, V. V. Application of L1-regularization approach in QSAR problem. Linear regres sion and artificial neural networks. Methods Objects Chem. Anal. 2019, 14(2), 79-90.
Zakharov A. B., Dyachenko A. V., Ivanov V. V. Topological characteristics of iterated line graphs in QSAR problem: Octane numbers of saturated hydrocarbons. J. Chemom. 2019, 33(9), e3169.
Zakharov A. B., Kyrpa M., Kyrychenko A. V., Kovalenko S. M., Kalugin O. N., Ivanov V. V., Adamowicz L. Towards the computational design of organic molecules with specified properties. Struct. Chem. 2025, 36(2), 723-738.
Papaconstantopoulos D. A., Mehl M. J. The Slater–Koster tight-binding method: a computationally efficient and accurate approach. J. Phys. Condens. Matter. 2003, 15(10), R413.
https://dftb.org/parameters/download.html
Gaussian 09, Revision A.02, M. J. Frisch, G. W. Trucks, H. B. Schlegel, G. E. Scuseria, M. A. Robb, J. R. Cheeseman, G. Scalmani, V. Barone, G. A. Petersson, H. Nakatsuji, X. Li, M. Caricato, A. Marenich, J. Bloino, B. G. Janesko, R. Gomperts, B. Mennucci, H. P. Hratchian, J. V. Ortiz, A. F. Izmaylov, J. L. Sonnenberg, D. Williams-Young, F. Ding, F. Lipparini, F. Egidi, J. Goings, B. Peng, A. Petrone, T. Henderson, D. Ranasinghe, V. G. Zakrzewski, J. Gao, N. Rega, G. Zheng, W. Liang, M. Hada, M. Ehara, K. Toyota, R. Fukuda, J. Hasegawa, M. Ishida, T. Nakajima, Y. Honda, O. Kitao, H. Nakai, T. Vreven, K. Throssell, J. A. Montgomery, Jr., J. E. Peralta, F. Ogliaro, M. Bearpark, J. J. Heyd, E. Brothers, K. N. Kudin, V. N. Staroverov, T. Keith, R. Kobayashi, J. Normand, K. Raghavachari, A. Rendell, J. C. Burant, S. S. Iyengar, J. Tomasi, M. Cossi, J. M. Millam, M. Klene, C. Adamo, R. Cammi, J. W. Ochterski, R. L. Martin, K. Morokuma, O. Farkas, J. B. Foresman, and D. J. Fox, Gaussian, Inc., Wallingford CT, 2016.