Semiempirical local spin: Theory and implementation of the ZILSH method for predicting Heisenberg exchange constants of polynuclear transition metal complexes