Comparison of Molecular Mechanics with GPU-Accelerated, DFT-Accuracy Quantum Mechanical Optimization via Differentiable GFN1-xTB and GFN2-xTB Methods
Abstract
Geometry optimization is the most critical initial step in computational chemistry. This study presents a GPU-accelerated geometry optimization workflow built upon dxtb, a fully differentiable PyTorch implementation of the semi-empirical GFN1-xTB method. The structural accuracy of the developed workflow was benchmarked against high-resolution X-ray crystal structures obtained from the Cambridge Structural Database (CSD), and compared with the MMFF94 force field and the GFN2-xTB method. For this purpose, a test set of 25 organic molecules representing distinct physical and chemical foundations was utilized: push-pull systems (p-nitroaniline), halogen-rich structures (hexachlorobenzene), rigid cage systems (cubane, adamantane, dodecahedrane, hexamethylenetetramine), flexible macrocyclic rings (18-crown-6), and polar disaccharides (sucrose). In the overall comparison based on heavy-atom RMSD values, GFN1-xTB exhibited the best mean performance (0.2686 Å), followed by MMFF94 (0.3119 Å) and GFN2-xTB (0.3818 Å). The superiority of GFN1-xTB is particularly pronounced in conjugated systems such as biphenyl, yielding a twofold increase in accuracy compared to MMFF94, and in highly flexible structures like cholesterol. Despite its more sophisticated electronic structure definition, GFN2-xTB lagged in mean RMSD, though it showed specific strengths in certain polycyclic and heteroaromatic systems; notably, its computational cost increased disproportionately for large molecules compared to GFN1-xTB. The results demonstrate that GFN1-xTB provides structural accuracy approaching the Density Functional Theory (DFT) level while significantly reducing computational costs compared to standard quantum mechanical methods. Although it may appear to lag behind MMFF94 in terms of absolute speed, GFN1-xTB offers highly practical utility for novel molecules, as it eliminates the need for predefined parameterization.
Keywords
GFN1-xTB, GFN2-xTB, DFT, MMFF94
Supporting Institution
Ethical Statement
Thanks
References
- Bannwarth, C., Caldeweyher, E., Ehlert, S., Hansen, A., Pracht, P., Seibert, J., Spicher, S., & Grimme, S. (2021). Extended tight-binding quantum chemistry methods. WIREs Computational Molecular Science, 11(2), Article e1493. https://doi.org/10.1002/wcms.1493
- Bannwarth, C., Ehlert, S., & Grimme, S. (2019). GFN2-xTB—An accurate and broadly parametrized self-consistent tight-binding quantum chemical method with multipole electrostatics and density-dependent dispersion contributions. Journal of Chemical Theory and Computation, 15(3), 1652–1671. https://doi.org/10.1021/acs.jctc.8b01176
- Bursch, M., Neugebauer, H., & Grimme, S. (2019). Structure optimisation of large transition-metal complexes with extended tight-binding methods. Angewandte Chemie International Edition, 58(32), 11078–11087. https://doi.org/10.1002/anie.201904021
- Bursch, M., Mewes, J.-M., Hansen, A., & Grimme, S. (2022). Best-practice DFT protocols for basic molecular computational chemistry. Angewandte Chemie International Edition, 61(42), Article e202205735. https://doi.org/10.1002/anie.202205735
- Ehlert, S., Stahn, M., Spicher, S., & Grimme, S. (2021). Robust and efficient implicit solvation model for fast semiempirical methods. Journal of Chemical Theory and Computation, 17(7), 4250–4261. https://doi.org/10.1021/acs.jctc.1c00471
- Friede, M., Hölzer, C., Ehlert, S., & Grimme, S. (2024). dxtb—An efficient and fully differentiable framework for extended tight-binding. The Journal of Chemical Physics, 161(6), Article 062501. https://doi.org/10.1063/5.0216715
- Grimme, S., Bannwarth, C., & Shushkov, P. (2017). A robust and accurate tight-binding quantum chemical method for structures, vibrational frequencies, and noncovalent interactions of large molecular systems parametrized for all spd-block elements (Z = 1–86). Journal of Chemical Theory and Computation, 13(5), 1989–2009. https://doi.org/10.1021/acs.jctc.7b00118
- Groom, C. R., Bruno, I. J., Lightfoot, M. P., & Ward, S. C. (2016). The cambridge structural database. Acta Crystallographica Section B: Structural Science, Crystal Engineering and Materials, 72(2), 171–179. https://doi.org/10.1107/S2052520616003954
- Halgren, T. A. (1996). Merck molecular force field. I. Basis, form, scope, parameterization, and performance of MMFF94. Journal of Computational Chemistry, 17(5–6), 490–519. https://doi.org/10.1002/(SICI)1096-987X(199604)17:5/6%3C490::AID-JCC1%3E3.0.CO;2-P
- Kabsch, W. (1976). A solution for the best rotation to relate two sets of vectors. Acta Crystallographica Section A: Foundations and Advances, 32(5), 922–923. https://doi.org/10.1107/S0567739476001873