Abstract:Accurate prediction of spin-state energy difference is crucial for understanding the spin crossover (SCO) phenomena and is very challenging for the density functional approximations, especially for the local and semi-local approximations, due to delocalization errors. Here, we investigate the effect of self-interaction error removal from the local spin density approximation (LSDA) and PBE generalized gradient approximation (GGA) on the spin-state gaps of Fe(II) complexes with various ligands using recently developed locally scaled self-interaction correction (LSIC) by Zope et al. [J. Chem. Phys. 151, 214108 (2019)]. The LSIC method is exact for one-electron density, which recovers uniform electron gas limit of underlying functional and approaches the well-known Perdew-Zunger self-interaction correction [Phys. Rev. B, 23, 5048 (1981)] (PZSIC) as a special case when the scaling factor is constant. Our results, when compared with reference diffusion Monte Carlo (DMC) results, show that the PZSIC method significantly overestimates spin-state gaps favoring low spin states for all ligands and does not improve upon DFAs. The perturbative LSIC-LSDA using PZSIC densities significantly improves the gaps with a mean absolute error of 0.51 eV but slightly overcorrects for the stronger CO ligands. The quasi-self-consistent LSIC-LSDA, like CCSD(T), gives a correct sign of spin-state gaps for all ligands with MAE of 0.56 eV, comparable to that of CCSD(T) (0.49 eV).

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: **Accurately predict the spin - state energy gap of iron(II) octahedral complexes** in order to better understand the Spin Crossover (SCO) phenomenon. Specifically, the paper focuses on how to improve the prediction accuracy of the spin - state energy gaps of these complexes by improving Density Functional Approximations (DFA), especially by removing the Self - Interaction Error (SIE). ### Problem Background 1. **Importance of Spin - State Energy Gap**: The spin - state energy gap is crucial for understanding the spin - crossover phenomenon. Under certain temperature changes, the spin - state of iron(II) complexes will change, which has an important impact on their physical and chemical properties. 2. **Limitations of Existing Methods**: Traditional Density Functional Theory (DFT) methods, especially local and semi - local approximations (such as LSDA and PBE), are not accurate enough in predicting the spin - state energy gap due to the existence of the Self - Interaction Error (SIE). SIE will lead to excessive delocalization of electrons, thus affecting the accuracy of spin - state ordering. ### Research Methods The authors used two self - interaction correction methods: 1. **Perdew - Zunger Self - Interaction Correction (PZSIC)**: This is a classical self - interaction correction method. 2. **Local Scaling Self - Interaction Correction (LSIC)**: This is a relatively new method that can handle the self - interaction error in multi - electron systems more precisely. ### Main Research Contents - **Comparing the Effects of Different Methods**: The authors compared the performance of LSDA, PBE, PZSIC and LSIC methods by calculating the spin - state energy gaps of a series of iron(II) octahedral complexes with different ligands. - **Evaluating the Accuracy of the Results**: The calculated results were compared with the reference Diffusion Monte Carlo (DMC) results to evaluate the Mean Absolute Error (MAE) of different methods. ### Main Findings 1. **Deficiencies of the PZSIC Method**: Although the PZSIC method can theoretically remove the self - interaction error, in practical applications, it tends to over - stabilize the low - spin state, resulting in inaccurate prediction of the spin - state energy gap, even worse than the uncorrected DFA method. 2. **Advantages of the LSIC Method**: The LSIC method significantly improves the prediction accuracy of the spin - state energy gap, especially for strong - field ligands (such as CO), and its Mean Absolute Error is significantly lower than other methods. In addition, the LSIC method shows better balance when dealing with complex systems. ### Conclusions - **Effectiveness of the LSIC Method**: The LSIC method can more effectively remove the self - interaction error, thereby providing more accurate prediction of the spin - state energy gap. - **Future Research Directions**: It is necessary to further study the performance of the LSIC method in other spin - crossover complexes to verify its wide applicability. Through these studies, the authors hope to provide a more reliable method for the accurate prediction of the spin - state energy gap, thereby promoting the research progress in related fields.

Spin-state Gaps and Self-Interaction-Corrected Density Functional Approximations: Octahedral Fe(II) Complexes as Case Study

Improving the Accuracy in the Prediction of Transition-Metal Spin-State Energetics Using a Robust Variation-Based Approach: Density Functional Theory, CASPT2 and MC-PDFT Applied to the Case Study of Tris-Diimine Fe(II) Complexes

Spin gaps and spin-flip energies in density-functional theory

Thermal spin crossover in Fe(ii) and Fe(iii). Accurate spin state energetics at the solid state

Conceptual Insights into DFT Spin-State Energetics of Octahedral Transition-Metal Complexes through a Density Difference Analysis.

Predicting Spin States of Iron Porphyrins with DFT Methods Including Crystal Packing Effects and Thermodynamic Corrections

Benchmarks and reliable DFT results for spin-crossover complexes

Appraising spin-state energetics in transition metal complexes using double-hybrid models: accountability of SOS0-PBESCAN0-2(a) as a promising paradigm

Symmetry breaking and self-interaction correction in the chromium atom and dimer

How Evenly Can Approximate Density Functionals Treat the Different Multiplicities and Ionization States of 4d Transition Metal Atoms?

Self-interaction correction schemes for non-collinear spin-density-functional theory

Importance of self-interaction-error removal in density functional calculations on water cluster anions

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

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

The rise and fall of stretched bond errors: Extending the analysis of Perdew–Zunger self-interaction corrections of reaction barrier heights beyond the LSDA

Performance of CASPT2 and DFT for Relative Spin-State Energetics of Heme Models

Spin-flip, Tensor Equation-of-motion Configuration Interaction with a Density-Functional Correction: A Spin-Complete Method for Exploring Excited-State Potential Energy Surfaces.

Spin–spin interactions in defects in solids from mixed all-electron and pseudopotential first-principles calculations

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