Density functional theory (in the form of the PW91, BP86, OLYP, and B3LYP exchangecorrelation functionals) has been used to map out the low-energy states of a series of eight-coordinate square-antiprismatic (D(2)d) first-row transition metal complexes, involving Mn(II), Fe(II), Co(II), Ni(II), and Cu(II), along with a pair of tetradentate N4 ligands. Of the five complexes, the Mn(II) and Fe(II) complexes have been synthesized and characterized structurally and spectroscopically, whereas the other three are as yet unknown. Each N4 ligand consists of a pair of terminal imidazole units linked by an o-phenylenediimine unit. The imidazole units are the strongest ligands in these complexes and dictate the spatial disposition of the metal three-dimensional orbitals. Thus, the dx(-)(2)y(2) orbital, whose lobes point directly at the coordinating imidazole nitrogens, has the highest orbital energy among the five d orbitals, whereas the d(xy) orbital has the lowest orbital energy. In general, the following orbital ordering (in order of increasing orbital energy) was found to be operative: d(xy) < d(xz) = d(yz) = d(z)(2) < d(x)(-y)(2)(2). The square-antiprism geometry does not lead to large energy gaps between the d orbitals, which leads to an S = 2 ground state for the Fe(II) complex. Nevertheless, the d(xy) orbital has significantly lower energy relative to that of the dxz and dyz orbitals. Accordingly, the ground state of the Fe(II) complex corresponds unambiguously to a d(xy)(2)d(xz)(1)d(yz)(1)d(z)(2)1d(x)(-)(2)y(2)1 electronic configuration. Unsurprisingly, the Mn(II) complex has an S = 5/2 ground state and no low-energy d-d excited states within 1.0 eV of the ground state. The Co(II) complex, on the other hand, has both a low-lying S = 1/2 state and multiple low-energy S = 3/2 states. Very long metal-nitrogen bonds are predicted for the Ni(II) and Cu(II) complexes; these bonds may be too fragile to survive in solution or in the solid state, and the complexes may therefore not be isolable. Overall, the different exchange-correlation functionals provided a qualitatively consistent and plausible picture of the low-energy d-d excited states of the complexes.