Assessment of Electronic Structure Methods for the Determination of Ground Spin States of Fe(II), Fe(III), and Fe(IV) Complexes

Verma, P.; Varga, Z.; Klein, J. E. M. N.; Cramer, C. J.; Que, L., Jr.; Truhlar, D. G.
Phys. Chem. Chem. Phys. 2017, 19, 13049 (doi:10.1039/c7cp01263b).

Our ability to understand and simulate the reactions catalyzed by iron depends strongly on our ability to predict the relative energetics of spin states. In this work, we studied the electronic structures of Fe²⁺ ion, gaseous FeO and 14 iron complexes using Kohn-Sham density functional theory with particular focus on determining the ground spin state of these species as well as the magnitudes of relevant spin-state energy splittings. The 14 iron complexes investigated in this work have hexacoordinate geometries of which seven are Fe(II), five are Fe(III) and two are Fe(IV) complexes. These are calculated using 20 exchange-correlation functionals. In particular, we use a local spin density approximation (LSDA) — GVWN5, four generalized gradient approximations (GGAs) — BLYP, PBE, OPBE and OLYP, two non-separable gradient approximations (NGAs) — GAM and N12, two meta-GGAs — M06-L and M11-L, a meta-NGA — MN15-L, five hybrid GGAs — B3LYP, B3LYP*, PBE0, B97-3 and SOGGA11-X, four hybrid meta-GGAs — M06, PW6B95, MPW1B95 and M08-SO and a hybrid meta-NGA — MN15. The density functional results are compared to reference data, which include experimental results as well as the results of diffusion Monte Carlo (DMC) calculations and ligand field theory estimates from the literature. For the Fe²⁺ ion, all functionals except M11-L correctly predict the ground spin state to be quintet. However, quantitatively, most of the functionals are not close to the experimentally determined spin-state splitting energies. For FeO all functionals predict quintet to be the ground spin state. For the 14 iron complexes, the hybrid functionals B3LYP, MPW1B95 and MN15 correctly predict the ground spin state of 13 out of 14 complexes and PW6B95 gets all the 14 complexes right. The local functionals, OPBE, OLYP and M06-L, predict the correct ground spin state for 12 out of 14 complexes. Two of the tested functionals are not recommended to be used for this type of study, in particular M08-SO and M11-L, because M08-SO systematically overstabilizes the high spin state, and M11-L systematically overstabilizes the low spin state.