“… , However, this localized unpaired electron can easily leak across the two dioxolene ligands during geometry optimization due to the common self-interaction error (SIE), − an effect where the electrons display unphysical interaction with themselves. Selection of a density functional that prevents such artificial delocalization of the spin density is vital, something not properly considered in the literature despite knowledge of this problem, with methods applied that did not localize the spin. ,,,, In general, generalized gradient approximations (GGA) functionals, such as BP86, , give reliable geometries and are widely used in many computational studies including for valence tautomeric cobalt complexes. ,, However, for [Co(dbdiox)(dbsq)(phen)], the BP86 functional results in delocalized spin distribution in the doublet LS-Co III -(cat)(SQ) tautomer resulting from high amounts of SIE (Figure ). The 2010 study by Yoshikawa et al employed B3LYP, B3LYP*, BLYP, , BP86 and OPBE , to obtain the relative energy of the VT process and Co–O bond lengths .…”