Previous research on the present complexes has shown that DFT has problems in predicting the accurate spin-state, equilibrium geometry, and spectrum unless specialized functionals are employed. In Ref. 51 it was shown that the OPBE functional yields acceptable results for these systems. For this reason, equilibrium geometries were taken from that reference.