An Examination of the Electron Densities in a Series of Tripodal Cobalt Complexes Bridged by Magnesium , Calcium , Strontium and Barium

X-ray crystallographic and theoretical charge density data for a series of compounds [(Co(Ts3tren))M(Co(Ts3tren))] (M = Mg, Ca, Sr and Ba) is examined. The crystal structures are isostructural and the alkaline earth metal ions have the same arrangement of donor oxygen atoms despite the large variation in ionic radii. The isomorphism of these molecules is surprising and a theoretical examination of the electronic structures, with the different metal ions along the series, provides detailed insight into their stabilities. The theoretical and experimental data are consistent and agree well. The local properties of the Co(II) ion and its donor atoms are relatively independent of the alkali earth metal.


Introduction
Spectroscopic and electrochemical data, and preliminary X-ray crystallographic data for a series of compounds-[(Co(Ts 3 tren))M(Co(Ts 3 tren))], (M = Mg, Ca, Sr and Ba)-were reported [1]].The results showed that the cobalt(II) anions were in an unusual four-coordinate environment, which was trigonal pyramidal.The fourth nitrogen atom above the triangle of three nitrogen atoms formed the apex of the trigonal pyramid, but the fifth site in the potential trigonal bipyramid was vacant; instead, a redox-inert alkaline earth metal was nearby.The alkaline earth metal bridged two trigonal-pyramidal cobalt(II) anions.The bridging symmetry was enforced by a crystallographic inversion center.Six oxygen atoms, three from three Ts groups of one [(Co(Ts 3 tren))] 1− complex and another three from the second [(Co(Ts 3 tren))] 1− complex, were bound to the alkaline earth metal.The triangles of the three oxygen atoms above and below the alkaline earth metal were twisted through an angle of 45 • relative to one another about the vertical axis (the Co-M-Co axis), forming a trigonal antiprism, and an overall octahedral arrangement of oxygen atoms.The crystal structures were isostructural, and the M 2+ ions had the same octahedral arrangement of six sulfonamide oxygen atoms.The M 2+ atoms were placed at the Wyckoff position −3, where the 3-fold axis and the inversion center intersect, for all the alkali-earth metals.The O atoms occupy the same crystallographic sites in the four structures.The isomorphism of these molecules was surprising because the magnesium and barium atoms are so different.Barium atoms favor much higher coordination numbers with weaker interactions.Because a six-coordinate Mg 2+ ion is much more stable than a six-coordinate Ba 2+ ion, it is reasonable to assume there was distortion and stress in the Ba-crystal.Thus, appreciable differences in the energies were expected across this series of compounds.Crystals of the Mg-compound were of the best quality, but most crystals of the Sr-and Ba-compounds were twinned, and their structures were disordered.The disordered crystals were made up of intimate components that were rotated about each other along the c-axis.We hypothesize that stress in compounds with larger ions may distort crystal lattices, causing crystallographic defects.Interestingly, although the neighboring domains possessed different orientations, splitting of the domains within the crystal coincided exactly with the c-axis, giving a hexagonal-packing motif.The disorder was borne out by twinning by inversion.Differences in the stabilities of the isomorphic compounds were also borne out by reactivity studies which showed the Ba-compound to be the least stable, and the most vulnerable to hydrolysis [1][2][3].The poor solubility of the Mg-complex precluded its inclusion in the hydrolysis experiments, but the observed trend was Ba < Sr < Ca, so it is reasonable to expect that the Mg-complex is also the least vulnerable to hydrolysis.We were interested to see whether the electronic nature of the cobalt ions, and whether the donor characteristics of the sulfonamide oxygen donors changed when the small, acidic Mg 2+ ion was replaced by bigger ions.This was especially of interest because not only does the identity of the alkaline earth metal play a prominent role in catalysis [2][3][4][5][6][7][8][9], but spectroscopic profiling of interactions between these redox-inert metals and cobalt-supporting ligands through electron-rich sulfonamide oxygen donors is challenging.Elucidating the origin of the redox-inert effect of metal ions is critical to the potential development of these compounds as catalysts for industrially important transformations.Spectroscopic and electrochemical techniques were employed to examine the cobalt chromophore, but these studies did not reveal any differences in the electronic structure of the cobalt ion, as a function of the identity of the alkaline earth metal.Therefore, charge-density studies of these compounds were performed using diffraction data from single crystals, and quantum-mechanical calculations.Electron densities were analyzed using the quantum theory of atoms in molecules (QTAIM) [10].Specifically, local properties at the bond critical points, and atomic charges for the cobalt ions, sulfonamide oxygen donors, and alkaline earth metals were compared and trends examined.Delocalization indices were calculated to examine the basic pattern of interactions between the alkaline earth metals and the oxygen donor atoms, and of interactions between Co atoms separated by the alkaline earth metal.The amount of electron sharing between atom basins (i.e., the degree of delocalization) was derived from the exchange density.The QTAIM atomic charges and delocalization indices between two atomic basins gave reasonable indications of the formal charges and bond orders between atoms in the complexes, except for the ionic M-O bonds, as evidenced by their similarity to the values derived from the simplified Lewis model of the chemical bonding.These results showed the effects of the M (Mg, Ca, Sr, and Ba) substitutions on the electron density distribution, and hence, their potential catalytic activity.We found that the alkaline earth metals had a very small observable effect on the transition metals, despite them having a dramatic effect on the reactivity and stability of the compounds.The energy differences of the isolated molecules were found to be small, with similar electron distributions about the cobalt atoms (except for the Ba-compound which had a stronger axial Co-N bond).

Materials and Methods
The theoretical electron densities of [C 36 H 78 MCo 2 N 8 O 12 S 6 ], where M was Mg, Ca, Sr, or Ba, were obtained using the M06-2X functional with a double-zeta plus polarization Douglas-Kroll-Hess (DZP-DKH) basis set.The M06-2X functional performed reasonably well when benchmarked against both coupled-cluster-single-double-perturbative-triple (CCSD(T)) and experimental results for structures, dissociation energies, and vibrational frequencies of alkaline-earth compounds [11].The DZP-DKH is a contracted Gaussian basis set of double-zeta-valence qualities plus polarization functions, with contraction coefficients optimized using the relativistic DKH Hamiltonian for elements H to Kr [12], and Cs to Rn [13].This theoretical level allowed us to analyze and compare the complexes featuring the four alkaline earth metals.The Kohn-Sham molecular orbitals of each complex were obtained with Gaussian 09 software [14], and the local and integrated properties of electron density were calculated with the AIMAll suite of programs [15].Quality X-ray charge-density data were accessible for the Mg-and Ca-complexes, but not for the other compounds due to issues related to crystal quality and twinning.Multipolar refinements [16] were applied with a multipole expansion of the valence shell, up to hexadecapoles, for Co, S, Mg, Ca, O, and N atoms.Hartree-Fock functions were used for the core-and spherical-valence densities.Conventional refinements, where the charge-density parameters-P v , P lm , and κ-were optimized on the experimentally measured X-ray structure factors, were carried out [17].The charge-density parameters-P v , P lm , κ-were also obtained for Co, S, O, and N atoms from non-relativistic wave-function calculations through a published procedure [18].In this procedure, density-functional calculations are performed on the experimental charge-density-optimized atomic coordinates.The program TONTO [19] was used to calculate structure factors from the theoretical electron density.The differences in orbital populations derived from the experimental data and from the theoretical densities were within 1%.

The Crystal Structures of [(Co(Ts 3 tren))M(Co(Ts 3 tren))]
Although the dimeric cobalt(II) complexes bridged by alkaline earth metals were isostructural, the smaller metals were more ordered, with the Mg-and Ca-compounds being the most ordered and well-crystallized of the series.The molecular structure of the Mg-compound, shown in Figure 1, will be discussed in the most detail.The molecule was located on a 3-fold rotation axis through atoms Mg(1), Co(1), and N (1).Each magnesium ion served as a bridge between two Co atoms, and formed a very interesting column along the c-axis.As these atoms were situated on the rotation axis, the entire N(1)-Co(1)-Mg-Co(1) -N(1) system was linear, within experimental error, by crystallographic symmetry.Both Co-Mg distances were equal to 3.408(1) Å.Each Co atom was coordinated to four N atoms in a trigonal-pyramidal arrangement with C 3v symmetry.The implication was that the symmetry of orbitals on the two cobalt atoms (and the apical N atoms) were constrained by crystallographic symmetry, and also had local C 3v symmetry.However, we did not impose this 3-fold symmetry in the quantum-mechanical calculations.On the other hand, we did retain the overall inversion symmetry of the molecule during the calculations.In the multipolar model, the spherical-harmonics functions with m = 1, 2, 4 were forbidden, but m = 0, 3 were allowed for this symmetry.An interesting feature, from a bonding and charge-density perspective, was that the cobalt atom was missing an apical ligand (i.e., there was no ligand forming a bipyramid for dative bonds with cobalt).The electron-rich sulfonamide oxygen atoms were a long distance from the Co atom (2.966(3) Å in the Ca-compound, and 2.863(3) Å in the Mg-compound).Instead, there was a M 2+ ion at the apical position.Interactions between Co and these sulfonamide oxygen atoms could not be ruled out, and one of the aims of this charge-density study was to determine the nature and the extent of these interactions.The M-O bonds were expected to be completely ionic, and it was unlikely for there to be any spin-spin coupling between Co(II) cations because of their large separation and with an absence of an spin-coupling pathway between them.The spin state was investigated using optimizations of unrestricted wave functions, and distributions of the spin densities in the molecule were visualized.The spin state was found to be nearly exactly S = 3/2 at each cobalt atom, corresponding to a high-spin d 7 configuration.The unrestricted value of S was 3.0023 for the entire complex, and nuclear magnetic resonance (NMR) spectroscopy confirmed the assignment of S = 3/2 [1].

Theoretical Charge-Density Analysis of [C 36 H 78 MCo 2 N 8 O 12 S 6 ]
Table 1 shows the molecular graphs [10,20] of [C 36 H 78 MCo 2 N 8 O 12 S 6 ] for each alkaline earth metal, where it was possible to observe, beside the covalent structure, an octahedral coordination sphere with six gradient paths between the alkaline earth metals and the sulfonamide oxygen atoms.
The small coordination sphere of Mg allowed the presence of a path between adjacent oxygen atoms, which were also involved in non-conventional hydrogen bonds [21].
Table 2 presents the local properties at the bond critical points of selected interactions included in the molecular graph.In general, the alkaline-earth-metal-oxygen interactions were classified as closed-shell interactions with low values of electron densities and of electronic energy densities, along with positive Laplacian operators of electron density [22].For closed shell interactions the curvature of the electron density along the interaction line is positive because charge is depleted at the interatomic surface.These local properties agreed with the ionic characteristics of the alkaline-earth-metal-oxygen interactions, and the configuration of the coordination sphere around the metal followed Gillespie's rules [23].There were some trends in the local properties.The charge depletion (∇ 2 ρ bcp ) decreased with atomic number (Z) of the alkaline earth metals, reducing the ionic characteristics of the interaction (∇ 2 ρ bcp = −0.0015Z + 0.2462, R 2 = 0.93).The same conclusion could be obtained from the energy density, where potential energy density increased with atomic number.Electron delocalization between alkaline earth metals and the oxygen atom increased with covalent contribution.The O-O contacts observed in the MgL 2 complex presented features of very weak, closed-shell interactions.On the other hand, the coordination sphere of the cobalt center was composed of four bond paths, with nitrogen atoms of the ligand in a trigonal-pyramidal configuration.On the basis of the properties of the bond critical points, Co-N interactions presented features of directional metal-ligand interactions [24].In general, equatorial ligands were strongly bonded to the cobalt atom (larger values of electron density, and its Laplacian operator), except in complexes with barium, which presented stronger axial metal-ligand interactions.
The small coordination sphere of Mg allowed the presence of a path between adjacent oxygen atoms, which were also involved in non-conventional hydrogen bonds [21].
Table 2 presents the local properties at the bond critical points of selected interactions included in the molecular graph.In general, the alkaline-earth-metal-oxygen interactions were classified as closed-shell interactions with low values of electron densities and of electronic energy densities, along with positive Laplacian operators of electron density [22].For closed shell interactions the curvature of the electron density along the interaction line is positive because charge is depleted at the interatomic surface.These local properties agreed with the ionic characteristics of the alkaline-earth-metal-oxygen interactions, and the configuration of the coordination sphere around the metal followed Gillespie's rules [23].There were some trends in the local properties.The charge depletion (∇ 2 ρbcp) decreased with atomic number (Z) of the alkaline earth metals, reducing the ionic characteristics of the interaction (∇ 2 ρbcp= −0.0015 Z + 0.2462, R 2 = 0.93).The same conclusion could be obtained from the energy density, where potential energy density increased with atomic number.Electron delocalization between alkaline earth metals and the oxygen atom increased with covalent contribution.The O-O contacts observed in the MgL2 complex presented features of very weak, closed-shell interactions.On the other hand, the coordination sphere of the cobalt center was composed of four bond paths, with nitrogen atoms of the ligand in a trigonal-pyramidal configuration.On the basis of the properties of the bond critical points, Co-N interactions presented features of directional metal-ligand interactions [24].In general, equatorial ligands were strongly bonded to the cobalt atom (larger values of electron density, and its Laplacian operator), except in complexes with barium, which presented stronger axial metal-ligand interactions.

Ba
The contour map of the Laplacian operators of electron density of MgL2, displayed in Figure 2, showed the spherical distribution of charge concentration within the valence shell of the magnesium atom.It was also possible to observe the polarization of the valence shell of the cobalt metal center, with the charge depletions directed toward the charge concentrations of the nitrogen atoms.An interesting feature of the molecular graphs was the interatomic region between the magnesium and cobalt atoms, where an axial charge depletion of the cobalt valence shell faced a (3, +3) critical point, which was a local minimum of the electron density.This feature was also observed in the other complexes (Table 3).

Ba
The contour map of the Laplacian operators of electron density of MgL2, displayed in Figure 2, showed the spherical distribution of charge concentration within the valence shell of the magnesium atom.It was also possible to observe the polarization of the valence shell of the cobalt metal center, with the charge depletions directed toward the charge concentrations of the nitrogen atoms.An interesting feature of the molecular graphs was the interatomic region between the magnesium and cobalt atoms, where an axial charge depletion of the cobalt valence shell faced a (3, +3) critical point, which was a local minimum of the electron density.This feature was also observed in the other complexes (Table 3).The contour map of the Laplacian operators of electron density of MgL 2 , displayed in Figure 2, showed the spherical distribution of charge concentration within the valence shell of the magnesium atom.It was also possible to observe the polarization of the valence shell of the cobalt metal center, with the charge depletions directed toward the charge concentrations of the nitrogen atoms.An interesting feature of the molecular graphs was the interatomic region between the magnesium and cobalt atoms, where an axial charge depletion of the cobalt valence shell faced a (3, +3) critical point, which was a local minimum of the electron density.This feature was also observed in the other complexes (Table 3).

Sr
Ba Table 4 presents selected atomic properties useful in showing the influence of the alkaline earth metals on the cobalt metal center.There was a nonlinear correlation between the atomic number of the alkaline earth metals and their electron delocalization indexes (DI) within the molecule as evidenced by the slope (y = −0.0003DI 2 + 0.0256 DI + 0.0044, R 2 = 0.99).The delocalization indexes of linked atoms through bond paths (DIB) presented the same correlation.Also, it was possible to observe correlations between non-bonded delocalization indexes (DINB) and the atomic number (DINB = 1 × 10 −05 Z 2 − 0.0015 Z + 0.2082, R 2 = 0.99).There was a correlation between the DIB of the alkaline earth metals (DIB(M)) and the DINB of the cobalt atom (DINB(Co)), where an increase in the DIB(M) was associated with a decrease in the DINB(Co).

Visualisation and Qualitative Analysis of the Deformation Density of the [(Co(Ts3tren))M(Co(Ts3tren))] Complexes
Visually, the qualitative characteristics of the deformation electron density (ED) appeared to be very similar for the Mg-, Ca-, Sr-, and Ba-compounds by visual inspection of their difference density maps.A deformation electron density map in the plane of the Ca 2+ atom and two of the oxygen Sr Ba Table 4 presents selected atomic properties useful in showing the influence of the alkaline earth metals on the cobalt metal center.There was a nonlinear correlation between the atomic number of the alkaline earth metals and their electron delocalization indexes (DI) within the molecule as evidenced by the slope (y = −0.0003DI 2 + 0.0256 DI + 0.0044, R 2 = 0.99).The delocalization indexes of linked atoms through bond paths (DIB) presented the same correlation.Also, it was possible to observe correlations between non-bonded delocalization indexes (DINB) and the atomic number (DINB = 1 × 10 −05 Z 2 − 0.0015 Z + 0.2082, R 2 = 0.99).There was a correlation between the DIB of the alkaline earth metals (DIB(M)) and the DINB of the cobalt atom (DINB(Co)), where an increase in the DIB(M) was associated with a decrease in the DINB(Co).

Visualisation and Qualitative Analysis of the Deformation Density of the [(Co(Ts3tren))M(Co(Ts3tren))] Complexes
Visually, the qualitative characteristics of the deformation electron density (ED) appeared to be very similar for the Mg-, Ca-, Sr-, and Ba-compounds by visual inspection of their difference density maps.A deformation electron density map in the plane of the Ca 2+ atom and two of the oxygen Sr Ba Table 4 presents selected atomic properties useful in showing the influence of the alkaline earth metals on the cobalt metal center.There was a nonlinear correlation between the atomic number of the alkaline earth metals and their electron delocalization indexes (DI) within the molecule as evidenced by the slope (y = −0.0003DI 2 + 0.0256 DI + 0.0044, R 2 = 0.99).The delocalization indexes of linked atoms through bond paths (DIB) presented the same correlation.Also, it was possible to observe correlations between non-bonded delocalization indexes (DINB) and the atomic number (DINB = 1 × 10 −05 Z 2 − 0.0015 Z + 0.2082, R 2 = 0.99).There was a correlation between the DIB of the alkaline earth metals (DIB(M)) and the DINB of the cobalt atom (DINB(Co)), where an increase in the DIB(M) was associated with a decrease in the DINB(Co).

Visualisation and Qualitative Analysis of the Deformation Density of the [(Co(Ts3tren))M(Co(Ts3tren))] Complexes
Visually, the qualitative characteristics of the deformation electron density (ED) appeared to be very similar for the Mg-, Ca-, Sr-, and Ba-compounds by visual inspection of their difference density maps.A deformation electron density map in the plane of the Ca 2+ atom and two of the oxygen Sr Ba Table 4 presents selected atomic properties useful in showing the influence of the alkaline earth metals on the cobalt metal center.There was a nonlinear correlation between the atomic number of the alkaline earth metals and their electron delocalization indexes (DI) within the molecule as evidenced by the slope (y = −0.0003DI 2 + 0.0256 DI + 0.0044, R 2 = 0.99).The delocalization indexes of linked atoms through bond paths (DIB) presented the same correlation.Also, it was possible to observe correlations between non-bonded delocalization indexes (DINB) and the atomic number (DINB = 1 × 10 −05 Z 2 − 0.0015 Z + 0.2082, R 2 = 0.99).There was a correlation between the DIB of the alkaline earth metals (DIB(M)) and the DINB of the cobalt atom (DINB(Co)), where an increase in the DIB(M) was associated with a decrease in the DINB(Co).

Visualisation and Qualitative Analysis of the Deformation Density of the [(Co(Ts3tren))M(Co(Ts3tren))] Complexes
Visually, the qualitative characteristics of the deformation electron density (ED) appeared to be very similar for the Mg-, Ca-, Sr-, and Ba-compounds by visual inspection of their difference density maps.A deformation electron density map in the plane of the Ca 2+ atom and two of the oxygen Table 4 presents selected atomic properties useful in showing the influence of the alkaline earth metals on the cobalt metal center.There was a nonlinear correlation between the atomic number of the alkaline earth metals and their electron delocalization indexes (DI) within the molecule as evidenced by the slope (y = −0.0003DI 2 + 0.0256 DI + 0.0044, R 2 = 0.99).The delocalization indexes of linked atoms through bond paths (DI B ) presented the same correlation.Also, it was possible to observe correlations between non-bonded delocalization indexes (DI NB ) and the atomic number (DI NB = 1 × 10 −5 Z 2 − 0.0015 Z + 0.2082, R 2 = 0.99).There was a correlation between the DI B of the alkaline earth metals (DI B (M)) and the DI NB of the cobalt atom (DI NB (Co)), where an increase in the DI B (M) was associated with a decrease in the DI NB (Co).Visually, the qualitative characteristics of the deformation electron density (ED) appeared to be very similar for the Mg-, Ca-, Sr-, and Ba-compounds by visual inspection of their difference density maps.A deformation electron density map in the plane of the Ca 2+ atom and two of the oxygen donor atoms, and in the plane of the Co 2+ atom and two of the nitrogen donor atoms is shown in Figure 3. Four charge concentrations corresponding to the d electrons about the cobalt atom were visible.These were asymmetric, and appeared to be more dominant on the side facing the Mg 2+ ion.A representation of the three-dimensional deformation density, showing the characteristics of the chemical bonding for the Ca-compound, is shown in Figure 4.There were concentrations of charge corresponding to the d and lone-pair electrons.The deformation density about M 2+ was very nearly spherical, showing the purely ionic characteristics of this atom.The contour map in a cross section through the Ca atom and the two O atoms showed the minimal perturbation of the ED by the oxygen donor atoms.However, there was a complex distribution of charge about the Co 2+ atom.There was an annular region of charge depletion in the xy-plane, containing the three equatorial N atoms.According to simple crystal field theory, d electrons in this region were repulsed by the electrostatic field of the N atoms, and the energy levels of the d electrons were expected to be higher, and their populations were expected to be lower in this region of the crystal field.In contrast, there were two nearly equivalent annular regions, with local charge concentrations at angles subtended by the equatorial xy-plane and the axial z-axis (containing the axial N atom).In each of these two regions, there were three charge maxima.Two nearly equivalent regions of charge depletion were present in the axial direction, directed both toward the axial N atom and toward the Ca atom.The charge depletions appeared to be smaller than the equatorial depletions.The two axial charge depletions were almost the same, although one was directed toward an axial ligand.This distribution of charge was what was expected on the basis of a simplified classical crystal-field picture of bonding.In trigonal-monopyramidal (C 3v ) geometry, we would expect the dx 2 − y 2 and dxy orbitals to be the most destabilized, with the dxz and dyz orbitals being degenerate and the lowest in energy, and the dz 2 orbital being in the middle.This arrangement does depend on the relative donor strength of the sulfonamide and amine donors.There was an annular region of charge depletion in the xy-plane, containing the three equatorial N atoms.According to simple crystal field theory, d electrons in this region were repulsed by the electrostatic field of the N atoms, and the energy levels of the d electrons were expected to be higher, and their populations were expected to be lower in this region of the crystal field.In contrast, there were two nearly equivalent annular regions, with local charge concentrations at angles subtended by the equatorial xy-plane and the axial z-axis (containing the axial N atom).In each of these two regions, there were three charge maxima.Two nearly equivalent regions of charge depletion were present in the axial direction, directed both toward the axial N atom and toward the Ca atom.The charge depletions appeared to be smaller than the equatorial depletions.The two axial charge depletions were almost the same, although one was directed toward an axial ligand.This distribution of charge was what was expected on the basis of a simplified classical crystal-field picture of bonding.In trigonal-monopyramidal (C3v) geometry, we would expect the dx 2 − y 2 and dxy orbitals to be the most destabilized, with the dxz and dyz orbitals being degenerate and the lowest in energy, and the dz 2 orbital being in the middle.This arrangement does depend on the relative donor strength of the sulfonamide and amine donors.

Comparison of Some Experimentally and Theoretically Derived Parameters
Theoretical and experimental multipolar parameters for cobalt are compared in Table 5. D orbital populations were obtained from multipole populations, and are given in Table 6.

Comparison of Some Experimentally and Theoretically Derived Parameters
Theoretical and experimental multipolar parameters for cobalt are compared in Table 5. D orbital populations were obtained from multipole populations, and are given in Table 6.

Conclusions
Accurate mapping of the electron density provided an experimental probe into the interactions around alkaline earth metals.Delocalization indexes between two atomic basins gave reasonable indications of formal charges and bond orders between atoms in the complexes.The differences in common atomic properties across the series were small, despite large differences in the stabilities of these compounds.The net atomic charges (Table 4) of the M 2+ ions were ranged from +1.78 to +1.93, and there was no clear trend in the variation of charges.Atomic charges for the cobalt(II) atoms were nearly identical across the series (+1.31).It is worth noting that there was an uneven increase in bond delocalization indexes across the series (Mg (0.2378), Ca (0.4058), Sr (0.5027), and Ba (0.4892)), with a sudden increase between magnesium and calcium.These indexes gave a measure of the amount of electron-sharing between the magnesium and oxygen donor atoms.The atomic charges of the oxygen donor atoms with respect to the M 2+ ions were very similar as evidenced by the nearly invariant value of −1.46(1) across the series.We found that the alkaline earth metals had a very small observable effect on the transition metals, despite having a dramatic effect on the reactivity and stability of the compounds.The energy differences of the isolated molecules were found to be small, with similar electron distributions existing about the cobalt atoms (except for the Ba-compound which had a stronger axial Co-N bond).The differences in reactivity need to be understood by examining the amount of energy released upon forming the reaction products (appreciably less energy may be released upon forming the hexa-aqua Mg 2+ ion than upon forming the fully hydrated Ba 2+ ion from the [(Co(Ts 3 tren))M(Co(Ts 3 tren))] compounds).However, there were subtle differences in the electron-density distributions in the Co-O-Mg-O-Co planes, in the bond critical points around the Co atom, and in the molecular graph of the Mg-compound.These slight variations indicated that a complete study of the reactivities of these compounds should include an analysis of the reaction products, and that QTAIM does not necessarily predict the relative stabilities of molecules a priori.

ρ
bcp = density at the bond critical point.∇ 2 ρ bcp = Laplacian of ρ bcp .ε = bond ellipticity.H bcp = the local total energy density at the the bond critical point (H bcp = G bcp + V bcp ).DI(M|O) = Electron delocalization index between atoms M and O.

Figure 2 .
Figure 2. A contour map of the of the Laplacian operators of electron density in the plane of the Mg, Co, and oxygen donor atoms in the dimeric Co 2+ complex, [(Co(Ts3tren))Mg(Co(Ts3tren))].The two small blue spheres are (3, +3) critical points, which are local minima of the electron density between the magnesium and cobalt atoms.

Figure 2 .
Figure 2. A contour map of the of the Laplacian operators of electron density in the plane of the Mg, Co, and oxygen donor atoms in the dimeric Co 2+ complex, [(Co(Ts 3 tren))Mg(Co(Ts 3 tren))].The two small blue spheres are (3, +3) critical points, which are local minima of the electron density between the magnesium and cobalt atoms.

Figure 2 .
Figure 2. A contour map of the of the Laplacian operators of electron density in the plane of the Mg, Co, and oxygen donor atoms in the dimeric Co 2+ complex, [(Co(Ts3tren))Mg(Co(Ts3tren))].The two small blue spheres are (3, +3) critical points, which are local minima of the electron density between the magnesium and cobalt atoms.

Figure 2 .
Figure 2. A contour map of the of the Laplacian operators of electron density in the plane of the Mg, Co, and oxygen donor atoms in the dimeric Co 2+ complex, [(Co(Ts3tren))Mg(Co(Ts3tren))].The two small blue spheres are (3, +3) critical points, which are local minima of the electron density between the magnesium and cobalt atoms.

Figure 2 .
Figure 2. A contour map of the of the Laplacian operators of electron density in the plane of the Mg, Co, and oxygen donor atoms in the dimeric Co 2+ complex, [(Co(Ts3tren))Mg(Co(Ts3tren))].The two small blue spheres are (3, +3) critical points, which are local minima of the electron density between the magnesium and cobalt atoms.

Figure 2 .
Figure 2. A contour map of the of the Laplacian operators of electron density in the plane of the Mg, Co, and oxygen donor atoms in the dimeric Co 2+ complex, [(Co(Ts3tren))Mg(Co(Ts3tren))].The two small blue spheres are (3, +3) critical points, which are local minima of the electron density between the magnesium and cobalt atoms.
q(A) = charge on atom A. N(A) = average electron population of atom A. LI(A) = average number of electrons localized on atom A. DI = average number of electrons in atom A that are shared between atom A and other atoms in the molecule.DIB = bonding contribution to the delocalization index DI.DINB = non-bonding contribution to the delocalization index DI.E(A) = electronic energy of atom A.

Crystals 2018, 8 ,
x FOR PEER REVIEW 8 of 12 q(A) = charge on atom A. N(A) = average electron population of atom A. LI(A) = average number of electrons localized on atom A. DI = average number of electrons in atom A that are shared between atom A and other atoms in the molecule.DIB = bonding contribution to the delocalization index DI.DINB = non-bonding contribution to the delocalization index DI.E(A) = electronic energy of atom A.

Ba
Crystals 2018, 8, x FOR PEER REVIEW 8 of 12 q(A) = charge on atom A. N(A) = average electron population of atom A. LI(A) = average number of electrons localized on atom A. DI = average number of electrons in atom A that are shared between atom A and other atoms in the molecule.DIB = bonding contribution to the delocalization index DI.DINB = non-bonding contribution to the delocalization index DI.E(A) = electronic energy of atom A.

Crystals 2018, 8 ,
x FOR PEER REVIEW 8 of 12 q(A) = charge on atom A. N(A) = average electron population of atom A. LI(A) = average number of electrons localized on atom A. DI = average number of electrons in atom A that are shared between atom A and other atoms in the molecule.DIB = bonding contribution to the delocalization index DI.DINB = non-bonding contribution to the delocalization index DI.E(A) = electronic energy of atom A.
q(A) = charge on atom A. N(A) = average electron population of atom A. LI(A) = average number of electrons localized on atom A. DI = average number of electrons in atom A that are shared between atom A and other atoms in the molecule.DI B = bonding contribution to the delocalization index DI.DI NB = non-bonding contribution to the delocalization index DI.E(A) = electronic energy of atom A.

Crystals 2018, 8 ,
x FOR PEER REVIEW 9 of 12 donor atoms, and in the plane of the Co 2+ atom and two of the nitrogen donor atoms is shown in Figure3.Four charge concentrations corresponding to the d electrons about the cobalt atom were visible.These were asymmetric, and appeared to be more dominant on the side facing the Mg 2+ ion.A representation of the three-dimensional deformation density, showing the characteristics of the chemical bonding for the Ca-compound, is shown in Figure4.There were concentrations of charge corresponding to the d and lone-pair electrons.The deformation density about M 2+ was very nearly spherical, showing the purely ionic characteristics of this atom.The contour map in a cross section through the Ca atom and the two O atoms showed the minimal perturbation of the ED by the oxygen donor atoms.However, there was a complex distribution of charge about the Co 2+ atom.

Figure 3 .
Figure 3. (a) A deformation electron density map in the plane of the Ca 2+ atom and two of the oxygen donor atoms.The density is almost perfectly spherical about the calcium, and two lobes of density corresponding to the two pairs of lone-pair electrons about each oxygen atom are visible.Red lines represent positive contours, and blue lines represent negative contours.Contours are drawn at 0.1 eÅ −3 intervals.(b) A deformation electron density map in the plane of the Co 2+ atom and two of the nitrogen donor atoms.Four charge concentrations corresponding to the d electrons about the cobalt atoms are visible, and two of these are larger towards the Mg 2+ ion.Red lines represent positive contours, and blue lines represent negative contours.Contours are drawn at 0.1 eÅ −3 intervals.

Figure 3 .
Figure 3. (a) A deformation electron density map in the plane of the Ca 2+ atom and two of the oxygen donor atoms.The density is almost perfectly spherical about the calcium, and two lobes of density corresponding to the two pairs of lone-pair electrons about each oxygen atom are visible.Red lines represent positive contours, and blue lines represent negative contours.Contours are drawn at 0.1 eÅ −3 intervals.(b) A deformation electron density map in the plane of the Co 2+ atom and two of the nitrogen donor atoms.Four charge concentrations corresponding to the d electrons about the cobalt atoms are visible, and two of these are larger towards the Mg 2+ ion.Red lines represent positive contours, and blue lines represent negative contours.Contours are drawn at 0.1 eÅ −3 intervals.

Figure 4 .
Figure 4. Positive and negative deformation-density isosurfaces (0.02 eÅ −3 ) of the dimeric Co 2+ complex, [C36H78CaCo2N8O12S6], from refinements with the predicted multipolar parameters.Charge concentrations are shown in red, and charge depletions are shown in blue.

Figure 4 .
Figure 4. Positive and negative deformation-density isosurfaces (0.02 eÅ −3 ) of the dimeric Co 2+ complex, [C 36 H 78 CaCo 2 N 8 O 12 S 6 ], from refinements with the predicted multipolar parameters.Charge concentrations are shown in red, and charge depletions are shown in blue.

Table 1 .
Molecular graph of the complexes.

Full Molecular Graph Partial Molecular Graph (Zoom) MgTable 1 .
Molecular graph of the complexes.
l Full Molecular Graph Partial Molecular Graph (Zoom) Mg Ca Sr Crystals 2018, 8, x FOR PEER REVIEW 5 of

Table 1 .
Molecular graph of the complexes.

Table 1 .
Molecular graph of the complexes.

Table 1 .
Molecular graph of the complexes.

Table 1 .
Molecular graph of the complexes.

Table 1 .
Molecular graph of the complexes.
Ba Crystals 2018, 8, x FOR PEER REVIEW 6 of 12

Table 2 .
Properties of the bond critical points, and delocalization of M-O and Co-N interactions.

Table 2 .
Properties of the bond critical points, and delocalization of M-O and Co-N interactions.

Table 2 .
Properties of the bond critical points, and delocalization of M-O and Co-N interactions.

Table 3 .
Molecular graph of the complexes.

Table 3 .
Molecular graph of the complexes.
MetalFull Molecular Graph Partial Molecular Graph (Zoom) Mg

Table 3 .
Molecular graph of the complexes.

Table 3 .
Molecular graph of the complexes.

Table 3 .
Molecular graph of the complexes.

Table 3 .
Molecular graph of the complexes.

Table 4 .
Integrated properties of the M 2+ and Co atoms.

Table 4 .
Integrated properties of the M 2+ and Co atoms.

Table 4 .
Integrated properties of the M 2+ and Co atoms.

Table 4 .
Integrated properties of the M 2+ and Co atoms.

Table 4 .
Integrated properties of the M 2+ and Co atoms.