Theoretical study on conformational energies of transition metal complexes

• Published in: Physical chemistry chemical physics : PCCP 2021-01, Vol.23 (1), p.287-299
• Main Authors: Bursch, Markus, Hansen, Andreas, Pracht, Philipp, Kohn, Julia T, Grimme, Stefan
• Format: Article
• Language: English
• Subjects: Benchmarks; Cartesian coordinates; Coordination compounds; Correlation coefficients; Crystal structure; Density functional theory; Low cost; Performance assessment; Production methods; Protocol (computers); Transition metal compounds
• ISSN: 1463-9076; 1463-9084
• DOI: 10.1039/d0cp04696e

Conformational energies are an important chemical property for which a performance assessment of theoretical methods is mandatory. Existing benchmark sets are often limited to biochemical or main group element containing molecules, while organometallic systems are generally less studied. A key problem herein is to routinely generate conformers for these molecules due to their complexity and manifold of possible coordination patterns. In this study we used our recently published CREST protocol [Pracht et al. , Phys. Chem. Chem. Phys. , 2020, 22 , 7169-7192] to generate conformer ensembles for a variety of 40 challenging transition metal containing molecules, which were then used to form a comprehensive conformational energy benchmark set termed TMCONF40. Several low-cost semiempirical, density functional theory (DFT) and force-field methods were compared to high level DLPNO-CCSD(T1) and double-hybrid DFT reference values. Close attention was paid to the energetic ordering of the conformers in the statistical evaluation. With respect to the double-hybrid references, both tested low-cost composite DFT methods produce high Pearson correlation coefficients of r p,mean,B97-3c//B97-3c = 0.922 and r p,mean,PBEh-3c//B97-3c = 0.890, with mean absolute deviations close to or below 1 kcal mol −1 . This good performance also holds for a comparison to DLPNO-CCSD(T1) reference energies for a smaller subset termed TMCONF5. Based on DFT geometries, the GFN n -xTB methods yield reasonable Pearson correlation coefficients of r p,mean,GFN1-xTB//B97-3c = 0.617 (MAD mean = 2.15 kcal mol −1 ) and r p,mean,GFN2-xTB//B97-3c = 0.567 (MAD mean = 2.68 kcal mol −1 ), outperforming the widely used PM x methods on the TMCONF40 test set. Employing the low-cost composite DFT method B97-3c on GFN2-xTB geometries yields an slightly improved correlation of r p,mean,B97-3c//GFN2-xTB = 0.632. Furthermore, for 68% of the investigated complexes at least one low-energy conformer was found that is more stable than the respective crystal structure conformation, which signals the importance of conformational studies. General recommendations for the application of the CREST protocol and DFT methods for transition metal conformational energies are given. Conformational energies are an important chemical property for which a performance assessment of theoretical methods is mandatory. Efficient low-cost methods are valuable for the generation and energetic ranking of conformers.