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.

Prediction of the inner-shell contribution to the correlation energy through DLPNO-CEPA/1 and Scaled same-spin second order Møller-Plesset perturbation theory

Accurate Prediction of Noncovalent Interaction Energies with the Effective Fragment Potential Method: Comparison of Energy Components to Symmetry-Adapted Perturbation Theory for the S22 Test Set

Efficient and accurate local approximations to coupled-electron pair approaches: An attempt to revive the pair natural orbital method

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

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

Toward Pair Atomic Density Fitting for Correlation Energies with Benchmark Accuracy

Assessing the domain-based local pair natural orbital (DLPNO) approximation for non-covalent interactions in sizable supramolecular complexes

Local electronic correlation in multicomponent Møller-Plesset perturbation theory

Staggered mesh method for correlation energy calculations of solids: Second order Møller-Plesset perturbation theory

Spinless formulation of linearized adiabatic connection approximation and its comparison with second order N-electron valence state perturbation theory

Factorized Quadruples and a Predictor of Higher-Level Correlation in Thermochemistry

Simultaneous Attenuation of Both Self-Interaction Error and Nondynamic Correlation Error in Density Functional Theory: A Spin-Pair Distinctive Adiabatic-Connection Approximation.

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

MP2-based composite extrapolation schemes can predict core-ionization energies for first-row elements with coupled-cluster level accuracy

Global Method for Electron Correlation

A regularized second-order correlation method from Green's function theory

Global Study of Beyond-Mean-Field Correlation Energies in Covariant Energy Density Functional Theory Using A Collective Hamiltonian Method

Density Functional Study Of Multiplicity-Changing Valence And Rydberg Excitations Of P-Block Elements: Delta Self-Consistent Field, Collinear Spin-Flip Time-Dependent Density Functional Theory (Dft), And Conventional Time-Dependent Dft

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

Combining density based dynamical correlation with a reduced density matrix strong correlation description