The role of vibrational anharmonicity in the computational study of thermal spin crossover

Abstract

Spin crossover in transition metal complexes can be studied in great detail with computational chemistry. Over the years, the understanding has grown that the relative stability of high-spin (HS) versus low-spin (LS) states is a subtle balance of many factors that all need to be taken into account for a reliable description. Among the different contributions, the zero-point energy (ZPE) and the entropy play key roles. These quantities are usually calculated assuming a harmonic oscillator model for the molecular vibrations. We investigated the impact of including anharmonic corrections on the ZPE and the entropy and indirectly on the critical temperature of spin crossover. As test systems, we used a set of ten Fe(II) complexes and one Fe(III) complex, covering different coordination modes (mono-, bi-, and tri-dentate ligands), decreasing coordination number upon spin crossover, coordination by second- and third-row atoms, and changes in the oxidation state. The results show that the anharmonicity has a measurable effect, but it is in general rather small, and tendencies are not easily recognized. As a conclusion, we put forward that for high precision results, one should be aware of the anharmonic effects, but as long as computational chemistry is still struggling with other larger factors like the influence of the environment and the accurate determination of the electronic energy difference between HS and LS, the anharmonicity of the vibrational modes is a minor concern.