A DFT computational study of spin crossover in iron(III) and iron(II) tripodal imidazole complexes. A comparison of experiment with calculations

B3LYP* functionals were used to model the sixteen iron(II) (1A, LS and 5T, HS) and iron(III) (2T, LS and 6A, HS) complexes of the 1:3 Schiff base condensate of tris(2-aminoethyl)amine and imidazole-4-carboxaldehyde, H3L1, and its deprotonated forms, H2L1]1-, HL1]2-, and L1]3-. This ligand system is unusual in that FeH3L1]3+, FeH3L1]2+ and FeL1]- all exhibit a spin crossover between 100-300 K. This makes these complexes ideal for a hybrid DFT computational approach and provides an opportunity to refine the value of the exact exchange admixture parameter, c3, and to predict properties of partially protonated complexes that are not experimentally available. The accepted value of 0.20 is larger than the value of approximately 0.13 that was found to best reproduce experimental data in terms of spin state predictions. With iron(III) B3LYP calculations showed that all of the complexes were low spin at 298 K with the exception of FeH3L1]3+ which is spin crossover in agreement with experimental results. It was also shown for iron(III) that the ligand field increased as the number of protons decreased. In contrast all of the iron(II) complexes were close to the spin crossover region regardless of protonation state. Experimental structures are fairly well modeled by this system in regard to the key structural indicators of spin state, which are the bite and trans angles. The calculated iron to nitrogen atom distances are always larger in the high spin form than the low spin form but all iron to nitrogen bond distances are larger than the experimental values. In general non-bonded interactions are not well modeled by this methodology.