Spin-orbit coupling has been introduced into our newly developed ligand field density functional theory (LFDFT), using the zero-order regular approximation as , is used to manifest further the effect of bonding changes on the sign and magnitude of the spin-orbit constant. Ligand field and spin-orbit coupling matrices are found to be correlated, with the higher erxtent of antibonding being accompanied by lower values of the spin-orbit coupling constant. In cases of little or no symmetry, this leads to situations in which ligand field and spin-orbit coupling cannot be neatly separated in the mathematical description. Using these results, the electronic energy levels of this series of compounds are predicted to be in good agreement with available spectral and magnetic data from literature.