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.

Can Spin-Component Scaled MP2 Achieve kJ/mol Accuracy for Cohesive Energies of Molecular Crystals?

Can spin-component scaled MP2 achieve kJ/mol accuracy for cohesive energies of molecular crystals?

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

Spin-free formulation of the multireference driven similarity renormalization group: A benchmark study of first-row diatomic molecules and spin-crossover energetics

Predicting Organic Crystal Lattice Energies with Chemical Accuracy

How Accurate are Simulations and Experiments for the Lattice Energies of Molecular Crystals?

How do spin-scaled double hybrids designed for excitation energies perform for noncovalent excited-state interactions? An investigation on aromatic excimer models

Crystal Polymorphism in Oxalyl Dihydrazide: Is Empirical DFT-D Accurate Enough?

Efficient Calculation of the Lattice Thermal Conductivity by Atomistic Simulations with Ab Initio Accuracy

Regularized and Opposite spin-scaled functionals from Møller-Plesset adiabatic connection -- higher accuracy at lower cost

Accurate interaction energies by spin component scaled Möller-Plesset second order perturbation theory calculations with optimized basis sets (SCS-MP2mod): Development and application to aromatic heterocycles

How do spin‐scaled double hybrids designed for excitation energies perform for noncovalent excited‐state interactions? An investigation on aromatic excimer models

Reaching high accuracy for energetic properties at second-order perturbation cost by merging self-consistency and spin-opposite scaling

A Multimer Embedding Approach for Molecular Crystals up to Harmonic Vibrational Properties

Effective spin-orbit models using correlated first-principles wave functions

Attaining high accuracy for charge-transfer excitations in non-covalent complexes at second-order perturbation cost: the importance of state-specific self-consistency

Towards accurate spin-orbit splittings from relativistic multireference electronic structure theory

Toward Accurate Spin-Orbit Splittings from Relativistic Multireference Electronic Structure Theory

Quantification of the thermodynamic effects of the low-spin – high-spin interaction in molecular crystals of a mononuclear iron(II) spin crossover complex