Abstract:Reliable prediction of spin-state energetics for transition metal (TM) complexes is recognized as a challenging and compelling problem in quantum chemistry, with implications for modeling catalytic reaction mechanisms and computational discovery of materials. The computed spin–state energetics are highly variable with the choice of method and credible reference data are scarce, making it difficult to conduct conclusive computational studies of open-shell TM systems. Here, we present a novel benchmark set of first-row TM spin–state energetics, which is derived from curated experimental data of 17 representative complexes containing Fe(II), Fe(III), Co(II), Co(III), Mn(II), and Ni(II) with chemically diverse ligands. The reference values of adiabatic or vertical energy differences are derived from spin-crossover enthalpies (9 complexes) or energies of spin-forbidden absorption bands in reflectance spectra (8 complexes). These are carefully back-corrected for relevant vibrational and environmental effects (due to solvation or crystal lattice) in order to provide the reference values directly comparable with computed electronic energy differences. The new benchmark set makes it possible to assess the accuracy of spin–state energetics from approximate density functional theory (DFT) and wave function methods with a level of statistical reliability not attained in earlier studies. The lowest mean absolute error (MAE) of 1.5 kcal/mol and maximum error of −3.5 kcal/mol are found for the coupled-cluster CCSD(T) method, which outperforms all tested multireference methods: CASPT2, MRCI+Q, CASPT2/CC and CASPT2+δMRCI. Contrary to earlier claims in the literature, the use of Kohn–Sham instead of Hartree–Fock orbitals in the reference determinant is not found to consistently improve the accuracy of the CCSD(T) spin–state energetics. The best performing DFT methods are double-hybrids (PWPB95-D3(BJ), B2PLYP-D3(BJ)) with the MAEs below 3 kcal/mol and maximum errors within 6 kcal/mol, whereas DFT methods traditionally recommended for spin states (e.g., B3LYP*-D3(BJ) and TPSSh-D3(BJ)) are found to perform much worse with the MAEs of 5–7 kcal/mol and maximum errors beyond 10 kcal/mol. The results of this work are relevant for the proper choice of methods to characterize TM systems in computational catalysis and (bio)inorganic chemistry, and may also stimulate new developments in quantum-chemical or machine learning approaches.

Systematic Error Estimation for Chemical Reaction Energies

Quantifying the error propagation in microkinetic modeling of catalytic reactions with model-predicted binding energies

Kinetics and Thermodynamics of Reactions Involving Criegee Intermediates: An Assessment of Density Functional Theory and Ab Initio Methods Through Comparison with CCSDT(Q)/CBS Data

A framework for quantifying uncertainty in DFT energy corrections

The performance of exchange—correlation functionals in describing electron density parameters of saddle point structures along chemical reactions

Error Assessment of Computational Models in Chemistry

Insight into Organic Reactions from the Direct Random Phase Approximation and Its Corrections.

Developing Orbital-Dependent Corrections for the Non-Additive Kinetic Energy in Subsystem Density Functional Theory

Understanding the Source of Error in First-Principles-Based Micro-Kinetic Modeling: Density Functional Theory Calculations Versus the Mean-Field Approximation

Accuracy of reaction coordinate based rate theories for modelling chemical reactions: insights from the thermal isomerization in retinal

Understanding density driven errors for reaction barrier heights

Putting error bars on density functional theory

Accurate transition state generation with an object-aware equivariant elementary reaction diffusion model

Computed vs Experimental Energy Barriers in Solution: Influence of the type of the density functional approximation.

A Systematic DFT Approach for Studying Mechanisms of Redox Active Enzymes

Computed versus experimental energy barriers in solution: Influence of the type of the density functional approximation

Performance of Quantum Chemistry Methods for Benchmark Set of Spin–State Energetics Derived from Experimental Data of 17 Transition Metal Complexes (SSE17)

Calculation of reaction constants using Transition Path Sampling with a local Lyapunov bias

Correction of Residual Errors in Configuration Interaction Electronic Structure Calculations

Refit to numerically problematic UMIST reaction rate coefficients

A Potential Energy Surface Construction Scheme For Accurate Reaction Rate Calculations: General Approach And A Test For The H+Ch4 -> H-2+Ch3 Reaction