Research Article  Open Access
Prinka Batra, Ritu Gaba, Upasana Issar, Rita Kakkar, "Structures and Stabilities of Alkaline Earth Metal Oxide Nanoclusters: A DFT Study", Journal of Theoretical Chemistry, vol. 2013, Article ID 720794, 14 pages, 2013. https://doi.org/10.1155/2013/720794
Structures and Stabilities of Alkaline Earth Metal Oxide Nanoclusters: A DFT Study
Abstract
The stability orders of a number of alkaline earth oxide cluster isomers , M = Mg, Ca, Sr, Ba and have been determined by means of density functional theory studies using the LDAPWC functional. Among the candidate structures, the hexagonalringbased isomers and the slab shapes are found to display similar stabilities. Stacks of hexagonal (MO)_{3} rings are found to be the slightly preferred growth strategy among the (MgO)_{6}, isomers. In contrast, the slab structures are slightly preferred for the other alkaline metal oxide (MO)_{6} clusters. An explanation based on packing and aromaticity arguments has been proposed. This study may have important implications for modeling and understanding the initial growth patterns of small nanostructures of alkaline earth metals.
1. Introduction
In the last few years, considerable effort has been directed to the understanding of metallic and semiconductor clusters. Clusters are aggregates of atoms or molecules intermediate in size between individual atoms and bulk matter, and their studies provide an interesting way to develop materials with varying properties by changing size and shape. Hence, studies of cluster properties as a function of size have received prominence in recent years. While much progress has been made on clusters of metals and semiconductors, metal oxide particles are often considered to be bulk fragments. However, their structure and properties could be entirely different in small clusters [1–4].
In this work, we have performed a comparative study of the structures, stabilities, and properties of some alkaline earth metal oxides (, , , and ). Magnesium oxide crystallizes in the rocksalt structure and has some typical semiconducting properties, such as wide valence band (~6 eV), large dielectric constant (9.8), and small exciton binding energy (<0.1 eV). For bulk MgO, the experimental value of the band gap is 7.8 eV [5]. It is close to an ideal insulating ionic solid with a valence band structure dominated by the strong potential of the ionic cores. Studies of the electronic properties of MgO are motivated by its technological applications, such as in catalysis, microelectronics, and electrochemistry. Bulk MgO is relatively inert, but its reactivity is greatly enhanced in the nanoscale. The high surface area and the intrinsically high surface reactivity of MgO nanocrystals make these materials especially effective as adsorbents [6]. In fact, they have been called “destructive adsorbents” because of their tendency to adsorb and simultaneously destroy by bond breaking processes a series of toxic chemicals [6–9].
It is interesting to study a similar system like calcium oxide in order to assess whether those trends are a general feature of alkaline earth oxide clusters or not. From the theoretical point of view, Ca^{2+} is larger than Mg^{2+}, so we can expect ionic size effects to play an important role in determining structural differences. Besides, Ca^{2+} is approximately six times as polarizable as Mg^{2+}, and the polarizabilities of the oxide anions are also larger in CaO because the bonding is weaker than in MgO. Calcium oxide also crystallizes in the closepacked “rocksalt” structure and is primarily an ionic material, with some degree of covalency in its bonding. It is considered as a prototype oxide from the theoretical point of view, with a wide band gap (7.1 eV) [10] and a high dielectric constant (11.8). Furthermore, local density approximation band structure calculations predicted a halfmetallic ferromagnetic ground state for CaO [11]. Nanocrystalline CaO is used as an absorbent to remove COD from paper mill effluent [12].
Barium oxide is an oxide with interesting electronic and structural properties. It is also a precursor to the ferroelectric perovskite oxide BaTiO_{3} and a component of the earth’s mantle. Barium strontium oxide coated carbon nanotubes serve as field emitters [13].
Theoretical work on ionic materials has been centered mostly in the family of alkali metal halides, and studies of metal oxide clusters have been comparatively scarce, despite their importance in many branches of surface physics, such as heterogeneous catalysis or corrosion. The mass spectra and collision induced fragmentation data for stoichiometric and cluster ions have been reported [14, 15]. The mass spectra of cluster ions [16] and experimental measurements of several singly and doubly ionized cluster ions of MgO and CaO by laser ionization timeofflight mass spectrometry [17–20] have also been published. Simple ionic models based on phenomenological pair potentials have been used to explain the global trends found in these experiments [17–20]. Several ab initio calculations on stoichiometric MgO clusters have been presented [21–32], but the growth of these clusters is still not well understood.
We aim to study the electronic properties of the clusters of these alkaline earth metal oxides using the density functional approach.
2. Computational Details
In the calculations reported in the paper, firstprinciples density functional (DF) calculations were performed using the DMol [5] code [33–36], available from Accelrys Inc. in the Materials Studio 3.2 package. The DFT calculations were carried out employing both the generalized gradient approximation (GGA) with the PW91 functional [37], as well as the local density approximation (LDA) with the PWC (PerdewWang local correlation) functional [38]. Hybrid functionals such as B3LYP, though more accurate for metal oxide dissociation energies, cannot be efficiently used with plane waves and are hence not useful for calculations on solid materials [39]. Our calculations employed numerical basis sets of doubleζ quality plus polarization functions (DNP) to describe the valence orbitals. This basis set is the numerical equivalent of the Gaussian basis, 631G**. The cores of Ba^{2+} and Sr^{2+} were treated with the allelectron approach.
Complete geometry optimizations for all structures were carried out. The atomic positions were relaxed to achieve minimum energy, until the system energy converged to 2 × 10^{−5} Ha and the gradient to 0.004 Ha Å^{−1}. The SCF tolerance was set at 1 × 10^{−5} and the maximum displacement at 0.005 Å. The binding energies (BEs), the highest occupied molecular orbital (HOMO)lowest unoccupied molecular orbital (LUMO) gaps, Fermi energies, and density of states were also computed. The reported binding energy values were corrected for zeropoint vibrational energies.
3. Results
Various structures, including the slab, hexagonal, octagonal, ladder, and other types, were studied for various numbers of formula units of the four alkaline earth metal oxides. Various theoretical studies at different levels of calculations have been reported in the literature [1, 2, 29, 32, 40–42], but there is no clear consensus regarding the suitability of LDA, GGA, or hybrid functionals for calculations on metal oxide nanoclusters. We therefore first compared results for the MgO molecular form obtained by different methods with the experimental quantities. The calculated LDAPWC, GGAPW91, B3LYP [40, 41], and MP4 [40, 41] values for the binding energy are 3.69, 3.22, 2.03, and 3.22 eV, respectively. The latter three fall short of the experimental value [43] of 3.57 eV. The LDA result is far superior to the other calculations, although its tendency to overbind is clear from the result. Similarly, the computed Mg–O bond lengths are 1.743, 1.767, and 1.756 Å, respectively, for LDAPWC, GGAPW91, and B3LYP [40, 41] calculations in comparison with the experimental [30, 31, 43] value of 1.749 Å. Here, again, our LDA result shows the best correspondence with experiment. The calculated vibrational frequencies are 751 cm^{−1} (LDA) and 719 cm^{−1} (GGA), in comparison with the observed [43] value of 785 cm^{−1}.
For CaO, too, the LDAPWC calculated binding energy (5.08 eV) is in better agreement with the experimental [43] value of 4.76 eV than the B3LYP/6311G(2d) [40, 41] value (4.28 eV). Though the GGAPW91 value (4.55 eV) is in slightly better agreement with experiment, the LDA Ca–O bond distance (1.818 Å) is in excellent agreement with the experimental [43] value (1.822 Å), while the GGA value is considerably larger (1.843 Å), reflecting the tendency of GGA to underbind. In the rest of the paper, therefore, we report the LDA results, but we also offer comparison with our calculated GGA results.
3.1. Stabilities of Structures
(MO)_{2} Where M = Mg, Ca, Sr, Ba. The optimized structures are rhombus shaped and planar (Figure 1). The angle about the magnesiums is obtuse (95.7°) in (MgO)_{2}, indicating the overlap repulsion between the large oxygen ions in close proximity since the cation size is small. However, for CaO, the bond angles are acute (∠OCaO = 86.4°). The bond angle about the metal ion decreases with increasing atomic number of the metal ion (∠OSrO = 82.7°, ∠OBaO = 79.3°), in parallel with the increasing ionic radii of the metal ions. The metal radii for Mg, Ca, Sr, and Ba are 1.60, 1.97, 2.15, and 2.17 Å, respectively, whereas the radii for the corresponding M^{2+} cations are 0.65, 0.97, 1.15, and 1.35 Å [44]. The atomic and ionic radii of O and lattice O^{2−} are 0.66 and 1.40 Å, respectively. The optimized M–O bond distances are 1.858, 2.005, 2.140, and 2.277 Å, respectively, for M = Mg, Ca, Sr, and Ba, which are all smaller than the sum of the ionic radii of M^{2+} and O^{2−}, indicating that these bonds are not purely ionic. In contrast, the observed M–O bond distances in the ionic crystals of the metal oxides are 2.106, 2.405, 2.565, and 2.762 Å, respectively. The corresponding gas phase values [43] for the diatomic species are 1.749, 1.822, 1.920, and 1.940 Å, respectively. The optimized M–O bond distances are slightly closer to the gas phase values than to the ionic values. The calculated Mulliken charges on the metal ions are 0.930, 1.253, 0.959, and 0.943, respectively. In view of the fact that the metal and oxygen charges are close to +1 and −1, respectively, the actual ionic radii are expected to lie between those for the neutral state and the divalent ions. The values of the partial charges indicate the slightly higher ionic character of (CaO)_{2} compared to the other (MO)_{2} systems.
Figure 1 gives the optimized structures for all the (MO)_{2} moieties. In this figure, as well as all other figures, the metal ions are represented by green balls and the oxide ions by red balls.
(MO)_{3} Where M = Mg, Ca, Sr, Ba. Table 1 summarizes the results for the various (MO)_{3} systems studied here. Here, MO stands for the four alkali metal oxides, MgO, CaO, BaO, and SrO. For , the two possible structures, namely, ladder and hexagonal, were studied for all the metal oxides.

(MgO)_{3}. Both ladder and hexagonal ring starting structures optimized to a distorted hexagonal form with a binding energy of −23.33 eV (see Table 1 and Figure 2(a)). It is well known that (MgO)_{3} ring structures are competitive building blocks in the growth of very small MgO clusters [29]. This result is also in agreement with our earlier calculations [22, 23] for (MgO)_{12} clusters, which optimized to stacked rings from initial cubic rock salt structures.
(a)
(b)
(c)
In this ring structure (Figure 2(a)), the optimized (LDAPWC) bond angle about each oxygen is 103.9°, and that about each Mg is 136.1°. The Mg–O bond distance is 1.819 Å. The increase of the ∠OMgO signifies repulsion between the oxygens. In fact, the O–O distances are 3.374 Å, compared to 2.866 Å for the Mg–Mg distance. The Mg–O bond orders are 0.872, implying significant covalent character. The Mulliken charge on each oxygen is −0.915. This signifies increased size of the oxygen ion in the system relative to the oxygen atom. Therefore, the energy is lowered by keeping the oxide ions away from each other. The GGAPW91 values are similar, although the computed Mg–O bond distance is much larger (1.843 Å), reflecting the tendency of GGA to underbind atoms. This structure has a high band gap (3.05 eV), which is much lower than that reported for bulk MgO [45–47].
(CaO)_{3}. For (CaO)_{3}, however, it was found that both initial structures optimized to different geometries. LDA calculations indicate that the ladder structure is slightly preferred, while GGA calculations predict a slight tilt in favor of the hexagonal structure (~0.1 eV). However, the energy differences are too small for one to make a definite statement regarding the relative stabilities. While the band gap for the hexagonal structure is 2.3 eV, that of the ladder structure is much smaller (1.8 eV).
In the ladder structure, there are two types of atoms—the central ones having a coordination number of 3, while the outer atoms having a coordination of 2 only and are more unsaturated. As a result, the Mulliken charge densities on the outer Ca atoms are 1.249, compared to 1.292 for the central Ca atom. Likewise, the terminal twocoordinate oxygens have a smaller negative charge (−1.240), while the central one has a partial charge of −1.311. The central bond length is also longer (2.294 Å) compared to the outer ones (1.952 Å). The increase in the central bond length can be understood in terms of the increased coordination of the central ions. The external field produced by the larger number of surrounding ions increases the ionic character of the central Ca–O bond, which resembles the Ca–O lattice limit (2.405 Å), while the terminal atoms are closer to the molecular limit (1.822 Å). The bond orders are 0.837 and 0.610, respectively. The Ca–O bond lengths and bond orders are 2.059 Å and 0.566, respectively, in the Ca–O–Ca face, and 2.108 Å and 0.662 in the O–Ca–O face. Again, this difference is due to the substantial ionic radius of O^{2−}, and the O–Ca–O face has two of these ions, compared to only one in the opposite face.
Figure 2(b) depicts the optimized structures. It can be seen that the bond angles about the oxygens in the ladder structure are obtuse, while those about the calciums are acute.
(SrO)_{3}. For (SrO)_{3}, the ladder structure is preferred by 0.16 eV. In this structure (Figure 2(c)), the outer Sr–O–Sr angle is obtuse (101.0°), but the O–Sr–O bond angle is acute (86.9°). The shortest Sr–O bond distances are the outer ones, that is, between two 2coordinate sites (2.107 Å), and the corresponding bond order is 1.128, while the longest bond is the central bond (2.450 Å) with a bond order of only 0.397. The Mulliken charges on the atoms show a behavior similar to that observed for (CaO)_{3}; that is, the charges on , , and are 0.915, 1.055, −0.952 and −0.981, respectively.
(BaO)_{3}. For (BaO)_{3}, again, two different structures are obtained. The ladder structure is slightly preferred (by 0.12 eV; see Table 1). The Mulliken charges on the atoms show a behavior similar to that observed for (CaO)_{3} and (SrO)_{3}; that is, the charges on , , , and are 0.880, 1.035, −0.915, and −0.965, respectively. It is observed that the charge on the central (3coordinate) metal ion decreases with increasing atomic number of the metal ion (−1.311, −0.981, and −0.965, for ladder (CaO)_{3}, (SrO)_{3}, and (BaO)_{3}, resp.), indicating the increasing involvement of the metal ion d orbitals and hence increasing covalency. As for SrO, the inner Ba–O bonds are longer (2.512 Å), and have a bond order of only 0.436, but the outer Ba–O bonds are shorter (2.179 Å), and the corresponding bond order is 1.160.
(MO)_{4} Where M = Mg, Ca, Sr, Ba. In this case, three structures, namely, slab, octagonal, and ladder, were studied. Table 2 gives the calculated energies and the HOMOLUMO gaps for the various structures. It is found that in all cases, the slab form is preferred over the other two. The energy differences are higher for this case, and the slab form is more emphatically preferred. Comparison of the calculated results is also made with B3LYP calculations [40, 41].
 
From reference [40, 41]. 
(MgO)_{4}. In the slab structure, all atoms are equivalent. We find that, from (MgO)_{4} onward, threedimensional structures are favored. (MgO)_{4} has a cubic structure with rhombohedral distortion (Figure 3(a)), each atom being tricoordinated. The Mulliken charge on Mg is 0.930, and all the Mg–O bond orders are 0.570. The bonding is therefore primarily ionic.
(a)
(b)
(c)
(CaO)_{4}. For CaO, the slab structure is again preferred. Although both LDA and GGA indicate this preference, it is interesting to see that the optimized structures for the octagonal initial structure are different for the two cases (Figure 3(b)) and different from the initial ring structure, unlike the case of (MgO)_{4} (Figure 3(a)). In the LDA case, the optimized structure consists of three fused rhombi, while the GGA optimized structure consists of fused six and fourmembered rings.
(SrO)_{4} and (BaO)_{4}. (SrO)_{4} and (BaO)_{4} show behavior similar to (CaO)_{4}, except that both LDA and GGA give similar optimized structures for the octagonal form. For this reason, only the optimized structures for (SrO)_{4} are shown in Figure 3(c).
An interesting result is that, although the slab structure is preferred in all cases, the next important structure is the ring for (MgO)_{4}, but for the other metal oxides, it is the ladder structure. Unlike (MgO)_{4}, the initial octagonal structure undergoes considerable distortion in all other cases.
(MO)_{5} Where M = Mg, Ca, Sr, Ba. The largest number of structures is possible in this case, namely, ladder, hexagonal, decagonal, chair, and many others (Table 3). In the case of MgO, both LDA and GGA predict that the most stable structure is the chair form. Similar is the case for CaO, that is, the chair form is the most stable structure, but in the case of SrO and BaO, the ladder form is found to be slightly preferred over the other forms. However, it may be noted that the energy differences are not very large for this stoichiometry.

(MgO)_{5}. Although MgO seems to prefer hexagonal structures, interestingly, in this case, the hexagonal fused ring optimizes to the ladder structure (Figure 4(a)). Similarly, the MgOI and MgOII starting structures optimize to the same geometry, which is a hybrid of one hexagonal ring fused with two fourmembered rings. The lowest energy structure for (MgO)_{5} is obtained from (MgO)_{4} by capping an edge such that the capping atoms are bicoordinated, while the rest of the atoms are tricoordinated. The corresponding edge on (MgO)_{4} opens up because of the increase in coordination number to three and consequent increase in ionic character. Similar is the case for the (CaO)_{5} structures, except that LDA predicts two distinct structures, CaOI and CaOII (Figure 4(b)). For (SrO)_{5} (Figure 4(c)) and (BaO)_{5} (Figure 4(d)), the ladder structures are preferred, and the chair structures could not be optimized.
(a)
(b)
(c)
(d)
(MO)_{6} Where M = Mg, Ca, Sr, Ba. In this case, three structures, namely, slab, hexagonal, and ladder, were studied (Table 4).
 
From reference [40, 41]. 
Table 4 indicates that in the case of MgO, the stacked hexagonal form is the most stable, and the slab structure is higher in energy by 0.21 eV. This energy difference increases to 0.57 eV for the (MgO)_{12} cluster [22, 23]. For the other metal oxides, the slab structure is preferred. The (MO)_{6} system is the first system for which both the slab structure and the hexagonal structure are possible, and we can make a comparison of the two. We observe that the Mg–O bond distances in the terminal rings of (MgO)_{6} are reduced when going from the rhombic (slab) to the hexagonal structure from 1.919 Å in the former to 1.891 Å in the latter. As for the ladder structures, the bond distance in the inner ring of the threering slab cluster (2.123 Å) is markedly elongated as compared to the terminal bond distance. As expected, the outer Mg–O distance is closer to the molecular 1.822 Å, while the inner ring distance is closer to the 2.106 Å lattice limit. Two distinctly different interplanar distances (1.936 Å and 1.898 Å) are also observed, depending on which atom, Mg or O, sits on the terminal ring of the threering stack. It is larger when O is on the terminal ring because of the larger O^{2−} radius. The interplanar distance in the hexagonal stacked structure (1.980 Å) is much larger than the Mg–O bond distance. The increased charge separation in the interior has only a minor effect on the terminal rings beyond that already seen for the doublering (MgO)_{4} system (1.943 Å).
For CaO, we do not observe any Ca–O bond compression with increasing number of atoms in the terminal rings; that is, the bond lengths do not vary too much when going from the slab (2.108 Å) to the hexagonal ring (2.106 Å) structure. There is, however, an interior ring expansion (2.290 Å) in the slab structure, similar to the MgO system, due to increased polarization of the Ca–O bond under the influence of the terminal rings. The terminal Ca–O bond distance in the threering stack is observed to be the mean of the 2.405 Å lattice value and the molecular 1.822 Å distance [43]. The innerring Ca–O bond distance (2.288 Å) is considerably closer to the bulk limit.
This behavior continues down the series. For (SrO)_{6}, the optimized interior Sr–O bond distance (2.422 Å) is similar to the lattice value (2.565 Å), while the terminal ring distance (2.245 Å) is smaller and closer to the gas phase value (1.920 Å). For (BaO)_{6}, the inner and outer Ba–O bond distances are 2.564 Å and 2.379 Å, respectively, compared to the gas phase and bulk values of 1.940 Å and 2.762 Å, respectively. An interesting trend is also observed. As the atomic number of the metal ion increases, the interplanar distance in the stacked hexagonal ring approaches the M–O distance in the ring, suggesting that spherical clusters become more important for the heavier alkaline earth metal oxides. For example, for (MgO)_{6}, the ring and interplanar distances are 1.891 and 1.980 Å, respectively. The corresponding distances are 2.118 Å and 2.159 Å (CaO), 2.272 Å and 2.290 Å (SrO), and 2.422 Å and 2.431 Å (BaO).
The inner ring bond distance is also only slightly larger than the bulk value for the (MgO)_{6} slab, and the deviation from the bulk value increases with increasing atomic number of the metal, suggesting a faster overall convergence to bulk properties for MgO clusters than for other alkaline metal oxide clusters.
Figure 5 depicts the optimized structures for the (MO)_{6} clusters for different starting geometries.
(a)
(b)
(c)
(d)
3.2. Vibrational Frequencies
For small MgO clusters, the experimental [48] vibrational frequency at 640 cm^{−1} matches with a strong resonance observed at 651 cm^{−1} in the highresolution electron energy loss spectroscopy (HREELS) surface phonon spectrum of the solid surface [49, 50]. The present calculations indicate that the vibrational frequency for this optically allowed stretching motion of oxygen atoms perpendicular to the surface of the (MgO)_{6} hexagonal cluster is at 691 cm^{−1} (intensity = 512 km mol^{−1}), while for the slab structure, this band is at 740 cm^{−1} (intensity = 445 km mol^{−1}). For the (MgO)_{12} nanocluster, the corresponding values [22, 23] are 657 cm^{−1} (intensity = 1132 km mol^{−1}) (hexagonal) and 677 cm^{−1} (intensity = 705 km mol^{−1}) (slab).
3.3. Electronic Structures
As noted previously, the preferred geometry for the (MgO)_{6} cluster is tubular, which becomes more stable than the rectangular bulklike cluster due to stabilization of the occupied levels. The isosurfaces of the HOMO and LUMO reveal that these comprise mainly the 2p orbitals of oxygen and 3s and 3p orbitals of Mg, respectively, but the LUMO also has an oxygen 3s component (see Figure 6). Thus, the anioncentered nature of the HOMO indicates that its energy depends strongly on the O–O distances. In the hexagonal structure, the O–O distance is 3.45 Å, compared to 2.86 Å for the rectangular cluster. The smaller anionanion repulsion in the nanotube stabilizes the HOMO, increasing the HOMOLUMO gap.
The electronic density of states (DOS) near the Fermi level for the two structures is shown in Figure 7. The separation of bands for the hexagonal structure is about 3.3 eV, which is higher than that obtained for the slab cluster (2.9 eV), but considerably smaller than the value of 4.8 eV computed for bulk MgO.
(a)
(b)
The DOS plots for the nanotube and cubelike structure are qualitatively similar, but one important difference is noticeable. In the hexagonal structure, there is greater involvement of the Mg^{2+} ion 3d orbitals near the Fermi level, although the population analysis reveals that the Mg electron configurations in the two are similar (3s^{0.34}2p^{6}3p^{0.46}3d^{0.26} for the hexagonal structure and 3s^{0.36}2p^{6}3p^{0.45}3d^{0.26} for the terminal slab ions). The peak height for the slab structure is also smaller. The computed electron configuration of the magnesium ion differs from the expected 2p^{6}. The 3s, 3p, and 3d orbitals are occupied, leading to a formal charge on Mg closer to one than the expected two.
The densities of states were also calculated for the slab structures of the other (MO)_{6} systems. These are shown in Figure 8.
(a)
(b)
(c)
It can be seen from Figure 8 that the involvement of d orbitals increases with increasing atomic number of the metal ion. Indeed, the participation of d orbitals increases from 0.263 in the terminal atoms of the MgO slab to 0.406 in CaO, 0.715 in SrO, and 0.835 for BaO. This happens at the expense of the p orbital population, so that the overall charge on the metal ion remains close to +1. For the inner metal ions, the participation of d orbitals is a little smaller. The additional electron comes from a 2p orbital of oxygen, and the oxygen atoms have electron configurations close to 1s^{2}2s^{2}2p^{5} in all cases. Another noticeable feature is that, while the bands are sharp for the other (MO)_{6} systems, they are broad for (CaO)_{6}. In fact, CaO appears to exhibit some anomalous behavior, as the Mulliken charge on Ca is slightly higher than that calculated for the metal ions in the other metal oxides.
4. Discussion
The bond lengths and total energy per molecule for the lowest energy structures are given in Table 5. The binding energy for a single MgO molecule is quite low in comparison to other clusters. The bond length (1.743 Å) is, however, small. As the cluster size increases, the bond lengths and binding energies increase in an oscillatory manner. The bond length is elongated to 1.882 Å in the most favorable structure of (MgO)_{2}, which is a rhombus. The binding energy per molecule also increases significantly, though this value is still much smaller than that for larger clusters. The lowest energy structure for (MgO)_{3} is planar and ring type, with each atom being bicoordinated. The Mg–O bond length in this case is shorter compared with the (MgO)_{2} cluster, and there is a significant increase in the binding energy. The bond angles about oxygen ions are smaller than those about magnesium ions, but the MgO bond lengths remain the same within both planar structures. In some of the structures, notably the chair structures, a large variation in M–O distances within the cluster is discernible from Table 5.

The increased stability of slab structures (both and ) is obvious from the binding energies per molecule for all systems (except , for which the hexagonal ringbased structures turn out to be more stable), in agreement with previous studies [29, 32, 40, 41]. For , our results that ring structures lie higher in energy are in contrast with reported results [17–20], in which and 5 were reported to have ring structures. For all (MO)_{4}, the slab structures are preferred over the octagonal ring and ladder structures. However, while for all the other (MO)_{4} system, the ring structure is the least preferred, for (MgO)_{4} it is preferred over the ladder structure. It may be noted that the ring structure undergoes considerable distortion in all cases except for MgO.
An important observation from the optimized threedimensional isomers is the chairtype structures of clusters. This is very interesting, as it indicates the existence of some covalent bonding in MgO clusters similar to covalently bonded silicon, although MgO is considered to be ionic in the bulk. Covalent bonding in alkalineearth oxides increases [51] as one goes from MgO to CaO, SrO, and BaO due to increasing involvement of d orbitals in bonding (Figure 8). While in bulk this is negligible for MgO, a reduction in the bond lengths in clusters could be responsible for the increased covalent character.
In the case of , there is a consistent preference for slabshaped structures, but the stability difference between the most stable and the second most stable structure is always small. Thus, hexagonal rings are slightly more stable than the slabshaped structures in the case of , whereas the opposite is true for . Consequently, we consider whether this possible trend continues, that is, toward increased relative stability of the slab structures relative to the hexagonal. Whereas the overall patterns of the Sr and Ba compounds are similar to Ca, we note in the case of (BaO)_{6} the lesser difference in stability between the slab and hexagonal ringbased isomers. Based on these observations, it can be concluded that there is no definite trend towards increasing preference for the slab shape with increasing atomic number on the metal for the small alkaline earth oxide clusters.
The similar values of the calculated Mulliken charges, approximately +1 on the metal ions in the various metal oxide clusters, make it difficult to assign an electronic cause for the slight preference for hexagonal structures in the case of MgO and slab structures for the other metal oxides. This leads us to believe that it is a packing effect rather than an electronic one. As stated in the sections above, due to the small cation size in MgO, the Mg–O bond is short, and, consequently, the fourmembered ring in the slab structure is too strained. In order to accommodate the small cation and the large anion in the fourmembered ring, the Mg–O bond length increases, leading to a weakening of the bonding and consequent instability. On the other hand, the octagonal ring structure is too open and is hence not favored for any of the metal oxides. This leaves the sixmembered ring as a compromise for MgO systems.
An alternate explanation for the preference for the sixmembered ring structure in the case of MgO could be the existence of aromaticity. In order to quantify aromaticity, we used the nucleusindependent chemical shift (NICS) method proposed by Schleyer et al. [52]. The NICS values were calculated at the center of the sixmembered ring (NICS(0)) and at a plane 1 Å above it (NICS(1)) and compared with the corresponding values calculated for benzene. In this method, negative NICS values indicate aromaticity and positive values antiaromaticity. We had earlier [53, 54] concluded that the NICS(1) value is the best measure of aromaticity for benzene. The calculated NICS(1) value for benzene is −10.84, whereas the corresponding values for (MgO)_{3} and (CaO)_{3} ring structures are −2.05 and 2.75, respectively, clearly indicating that the MgO ring is about 20% aromatic, while the CaO ring is antiaromatic, accounting for the increased stability of the MgO ring.
5. Conclusions
An important finding of the present study is that hexagonal tubelike structures are preferred for clusters, while slablike structures are preferred for the other alkaline earth metal oxide clusters. Explanations based on ionic size effects and aromaticity have been proposed in this work. It is gratifying to note that experimental observations of mass spectra [14, 15, 48, 55] indicate the existence and stabilities of such stacked hexagonal rings, at least for small gasphase clusters of MgO. Other experimental and theoretical [22, 23, 56] studies also provide evidence for the existence and stabilities of MgO nanotubes.
An outstanding result of the present study is the similar stabilities of the hexagonalringbased structures and the rocksaltlike slabshaped isomers. While this observation is important as such and contradicts the exclusive nature of the latter structural shapes proposed previously [17–20], it is noteworthy how the stability ordering changes as the metal atomic number increases among the alkaline earth elements. In the case of , the hexagonalringbased structure is the more stable one, although the energy difference between the two structures is small. Going to , the situation is reversed in that the slab structure prevails. For , the trend towards increasing stabilization of the slab structures continues, and the slab structure is evidently the more stable one. For , the slab structure is still the preferred one, but to a lesser extent. An explanation based on simple packing arguments has been proposed to explain the variation in relative stabilities. Aromaticity in the (MgO)_{3} ring also accounts for its stability.
It is difficult to find experimental verification for our results, as neutral clusters are difficult to study experimentally. Their structures are usually inferred indirectly from the mass spectra of ionized clusters, the more abundant species being interpreted as the more stable. However, the results from such studies on alkaline earth metal oxides are contradictory and depend on the process of formation of the clusters. Two conclusions, however, result from these studies. Firstly, for small clusters, hexagonal stacked rings are preferred for , but these give way to rocksalt cubic structures for large values of . Secondly, fragmentation clusters are found for both even values of and when is a multiple of 3. These results suggest that the basic clusterbuilding blocks are different for the two materials, as observed from the present calculations. The results of the experimental studies and our calculations can be reconciled if we assume that the neutral stoichiometric and clusters show structural differences: the basic building block is an (MgO)_{3} hexagonal fragment in the case of MgO and a (CaO)_{3} rectangular 2 × 3 (or doublechain) fragment for CaO, as the one found in the present studies. This difference is just a packing effect due to the larger overlap repulsion between anions in MgO since the cation size is very small and is also due to the aromaticity of the (MgO)_{3} ring.
These experiments also suggest that the hexagonal ring and rectangular slab structures are topologically equivalent. Deformation along one of the directions orthogonal to the rings stack transforms it into the slab structure. As noted earlier, this intense vibration mode for (MgO)_{6} occurs at a low wavenumber (691 cm^{−1}). Thus, experimental knowledge of abundance of masses alone cannot distinguish between the two structures, and sophisticated calculations such as the present ones can only decide the relative stabilities. Our earlier studies [22, 23] on the (MgO)_{12} cluster had indicated that the (MgO)_{12} nanotube, consisting of four stacked hexagonal (MgO)_{3} rings, is more stable than the bulklike cubic structure by 0.48 eV. Moreover, the calculated energy barrier for the rearrangement of the cubic structure to the tube was found to be only 0.13 eV, and hence the two structures are easily interconvertible.
Conflict of Interests
The authors declare no conflict of interests with any financial organization regarding the material presented in the paper.
Acknowledgments
The authors thank the Council of Scientific and Industrial Research (CSIR) for financial support (Grant no. 01(2554)/12/EMR1). Ritu Gaba and Upasana Issar thank the CSIR and the University Grants Commission (UGC), respectively, for Senior Research Fellowships.
References
 A. Jain, V. Kumar, M. Sluiter, and Y. Kawazoe, “First principles studies of magnesium oxide clusters by parallelized Tohoku University MixedBasis program TOMBO,” Computational Materials Science, vol. 36, no. 12, pp. 171–175, 2006. View at: Google Scholar
 N. Sharma and R. Kakkar, “Recent advancements in warfare agents/metal oxides surface chemistry and their simulant studies,” Advanced Materials Letters, 2013. View at: Google Scholar
 P. N. Kapoor, A. K. Bhagi, R. S. Mulukutla, and K. J. Klabunde, Dekker Encyclopedia of Nanoscience & Technology, Marcel Dekker, New York, NY, USA, 2004.
 A. Khaleel, P. N. Kapoor, and K. J. Klabunde, “Nanocrystalline metal oxides as new adsorbents for air purification,” Nanostructured Materials, vol. 11, no. 4, pp. 459–468, 1999. View at: Google Scholar
 J. Heyd, J. E. Peralta, and G. E. Scuseria, “Energy band gaps and lattice parameters evaluated with the HeydScuseriaErnzerhof screened hybrid functional,” Journal of Chemical Physics, vol. 123, no. 17, Article ID 174101, 8 pages, 2005. View at: Google Scholar
 S. Utamapanya, K. J. Klabunde, and J. R. Schlup, “Nanoscale metal oxide particles/clusters as chemical reagents,” Chemistry of Materials, vol. 3, pp. 175–181, 1991. View at: Google Scholar
 O. Koper, Y. X. Li, and K. J. Klabunde, “Destructive adsorption of chlorinated hydrocarbons on ultrafine (nanoscale) particles of calcium oxide,” Chemistry of Materials, vol. 5, no. 4, pp. 500–505, 1993. View at: Google Scholar
 Y. X. Li and K. J. Klabunde, “Nanoscale metal oxide particles as chemical reagents. Destructive adsorption of a chemical agent simulant, dimethyl methylphosphonate, on heattreated magnesium oxide,” Langmuir, vol. 7, no. 7, pp. 1388–1393, 1991. View at: Google Scholar
 Y.X. Li, J. R. Schlup, and K. J. Klabunde, “Fourier transform infrared photoacoustic spectroscopy study of the adsorption of organophosphorus compounds on heattreated magnesium oxide,” Langmuir, vol. 7, no. 7, pp. 1394–1399, 1991. View at: Google Scholar
 R. C. Whited, C. J. Flaten, and W. C. Walker, “Exciton thermoreflectance of MgO and CaO,” Solid State Communications, vol. 13, no. 11, pp. 1903–1905, 1973. View at: Google Scholar
 I. S. Elfimov, S. Yunoki, and G. A. Sawatzky, “Possible path to a new class of ferromagnetic and halfmetallic ferromagnetic materials,” Physical Review Letters, vol. 89, no. 21, Article ID 216403, 4 pages, 2002. View at: Google Scholar
 B. Nagappa and G. T. Chandrappa, “Nanocrystalline CaO as adsorbent to remove COD from paper mill effluent,” Journal of Nanoscience and Nanotechnology, vol. 7, no. 3, pp. 1039–1042, 2007. View at: Publisher Site  Google Scholar
 F. Jin, Y. Liu, and M. D. Christopher, “Barium strontium oxide coated carbon nanotubes as field emitters,” Applied Physics Letters, vol. 90, no. 14, Article ID 143114, 3 pages, 2007. View at: Google Scholar
 W. A. Saunders, “Structural dissimilarities between small IIVI compound clusters: MgO and CaO,” Physical Review B, vol. 37, no. 11, pp. 6583–6586, 1988. View at: Google Scholar
 W. A. Saunders, “Molecules and clusters,” Zeitschrift für Physik D, vol. 12, pp. 601–603, 1989. View at: Google Scholar
 T. P. Martin and T. Bergmann, “Mass spectra of CaO and BaO clusters,” The Journal of Chemical Physics, vol. 90, no. 11, pp. 6664–6667, 1989. View at: Google Scholar
 P. J. Ziemann and A. W. Castleman, “Massspectrometric study of the formation, evaporation, and structural properties of doubly charged MgO clusters,” Physical Review B, vol. 44, no. 12, pp. 6488–6499, 1991. View at: Publisher Site  Google Scholar
 P. J. Ziemann and A. W. Castleman, “Structures and bonding properties of calcium oxide clusters inferred from mass spectral abundance patterns,” The Journal of Physical Chemistry, vol. 96, no. 11, pp. 4271–4276, 1992. View at: Google Scholar
 P. J. Ziemann and A. W. Castleman Jr., “Mass spectrometric study of MgO clusters produced by the gas aggregation technique,” Zeitschrift für Physik D, vol. 20, pp. 97–99, 1991. View at: Google Scholar
 P. J. Ziemann and A. W. Castleman Jr., “Stabilities and structures of gas phase MgO clusters,” Journal of Chemical Physics, vol. 94, no. 1, pp. 718–728, 1991. View at: Google Scholar
 M. Gutowski, P. Skurski, X. Li, and L.S. Wang, “${(\text{MgO})}_{\text{n}}{\text{\hspace{0.17em}}}^{}$ (n = 1–5) clusters: multipolebound anions and photodetachment spectroscopy,” Physical Review Letters, vol. 85, no. 15, pp. 3145–3148, 2000. View at: Publisher Site  Google Scholar
 R. Kakkar, P. N. Kapoor, and K. J. Klabunde, “Theoretical study of the adsorption of formaldehyde on magnesium oxide nanosurfaces: size effects and the role of lowcoordinated and defect sites,” The Journal of Physical Chemistry B, vol. 108, no. 47, pp. 18140–18148, 2004. View at: Google Scholar
 R. Kakkar, P. N. Kapoor, and K. J. Klabunde, “First principles density functional study of the adsorption and dissociation of carbonyl compounds on magnesium oxide nanosurfaces,” Journal of Physical Chemistry B, vol. 110, no. 51, pp. 25941–25949, 2006. View at: Publisher Site  Google Scholar
 R. Dong, X. Chen, X. Wang, and W. Lu, “Structural transition of hexagonal tube to rocksalt for (MgO)_{3n}, $2\le n\le 10$,” Journal of Chemical Physics, vol. 129, no. 4, Article ID 044705, 5 pages, 2008. View at: Google Scholar
 M. Srnec and R. Zahradník, “Small group IIaVIa clusters and related systems: a theoretical study of physical properties, reactivity, and electronic spectra,” European Journal of Inorganic Chemistry, vol. 2007, no. 11, pp. 1529–1543, 2007. View at: Google Scholar
 M. Wilson, “Stability of small MgO nanotube clusters: predictions of a transferable ionic potential model,” The Journal of Physical Chemistry B, vol. 101, no. 25, pp. 4917–4924, 1997. View at: Google Scholar
 M. Wilson, P. A. Madden, N. C. Pyper, and J. H. Harding, “Molecular dynamics simulations of compressible ions,” Journal of Chemical Physics, vol. 104, no. 20, pp. 8068–8081, 1996. View at: Google Scholar
 E. De La Puente, A. Aguado, A. Ayuela, and J. M. López, “Structural and electronic properties of small neutral (MgO)_{n} clusters,” Physical Review B, vol. 56, no. 12, pp. 7607–7614, 1997. View at: Google Scholar
 M.J. Malliavin and C. Coudray, “Ab initio calculations on (MgO)_{n}, (CaO)_{n}, and (NaCl)_{n} clusters (n = 1–6),” Journal of Chemical Physics, vol. 106, no. 6, pp. 2323–2330, 1997. View at: Google Scholar
 J. M. Recio, R. Pandey, A. Ayuela, and A. B. Kunz, “Molecular orbital calculations on (MgO)_{n} and ${(\text{MgO})}_{\text{n}}{\text{\hspace{0.17em}}}^{+}$ clusters (n = 1–13),” The Journal of Chemical Physics, vol. 98, no. 6, pp. 4783–4792, 1993. View at: Google Scholar
 J. M. Recio and R. Pandey, “Ab initio study of neutral and ionized microclusters of MgO,” Physical Review A, vol. 47, no. 3, pp. 2075–2082, 1993. View at: Google Scholar
 A. Aguado and J. M. López, “Structures and stabilities of CaO and MgO clusters and cluster ions: an alternative interpretation of the experimental mass spectra,” The Journal of Physical Chemistry B, vol. 104, no. 35, pp. 8398–8405, 2000. View at: Google Scholar
 B. Delley, “An all electron numerical method for solving the local density functional for polyatomic molecules,” Journal of Chemical Physics, vol. 92, no. 1, pp. 508–517, 1990. View at: Google Scholar
 B. Delley, “Analytic energy derivatives in the numerical local density functional approach,” Journal of Chemical Physics, vol. 94, no. 11, pp. 7245–7250, 1991. View at: Google Scholar
 B. Delley, “Fast calculation of electrostatics in crystals and large molecules,” The Journal of Physical Chemistry, vol. 100, no. 15, pp. 6107–6110, 1996. View at: Google Scholar
 B. Delley, “From molecules to solids with the DMol^{3} approach,” Journal of Chemical Physics, vol. 113, no. 18, pp. 7756–7764, 2000. View at: Google Scholar
 J. P. Perdew, J. A. Chevary, S. H. Vosko et al., “Atoms, molecules, solids, and surfaces: applications of the generalized gradient approximation for exchange and correlation,” Physical Review B, vol. 46, no. 11, pp. 6671–6687, 1992. View at: Publisher Site  Google Scholar
 J. P. Perdew and Y. Wang, “Accurate and simple analytic representation of the electrongas correlation energy,” Physical Review B, vol. 45, no. 23, pp. 13244–13249, 1992. View at: Publisher Site  Google Scholar
 R. Kakkar, R. Grover, and P. Gahlot, “Density functional study of the properties of isomeric aminophenylhydroxamic acids and their copper (II) complexes,” Polyhedron, vol. 25, no. 3, pp. 759–766, 2006. View at: Publisher Site  Google Scholar
 F. Bawa and I. Panas, “Competing pathways for MgO, CaO, SrO, and BaO nanocluster growth,” Physical Chemistry Chemical Physics, vol. 4, pp. 103–108, 2002. View at: Google Scholar
 Á. Vibók and G. J. Halász, “Parametrization of complex absorbing potentials for timedependent quantum dynamics using multistep potentials,” Physical Chemistry Chemical Physics, vol. 3, pp. 3042–3047, 2001. View at: Google Scholar
 S. Veliah, R. Pandey, Y. S. Li, J. M. Newsam, and B. Vessal, “Density functional study of structural and electronic properties of cubelike MgO clusters,” Chemical Physics Letters, vol. 235, no. 12, pp. 53–57, 1995. View at: Google Scholar
 K. P. Huber and G. Herzberg, Molecular Spectra and Molecular Structure, IV. Constants of Diatomic Molecules, van Nostrand Reinhold, New York, NY, USA, 1979.
 M. C. Day and J. Selbin, Theoretical Inorganic Chemistry, chapter 4, Reinhold, New York, NY, USA; Chapman & Hall, London, UK, 1962.
 D. J. Driscoll, W. Martir, J. X. Wang, and J. H. Lunsford, “Formation of gasphase methyl radicals over magnesium oxide,” Journal of the American Chemical Society, vol. 107, no. 1, pp. 58–63, 1985. View at: Google Scholar
 T. Ito and J. H. Lunsford, “Synthesis of ethylene and ethane by partial oxidation of methane over lithiumdoped magnesium oxide,” Nature, vol. 314, pp. 721–722, 1985. View at: Google Scholar
 T. Ito, J. X. Wang, C. H. Lin, and J. H. Lunsford, “Oxidative dimerization of methane over a lithiumpromoted magnesium oxide catalyst,” Journal of the American Chemical Society, vol. 107, no. 18, pp. 5062–5068, 1985. View at: Google Scholar
 D. Van Heijnsbergen, G. Von Helden, G. Meijer, and M. A. Duncan, “Infrared resonanceenhanced multiphoton ionization spectroscopy of magnesium oxide clusters,” Journal of Chemical Physics, vol. 116, no. 6, pp. 2400–2406, 2002. View at: Publisher Site  Google Scholar
 P. A. Cox and A. A. Williams, “HREELS studies of simple ionic solids,” Journal of Electron Spectroscopy and Related Phenomena, vol. 39, pp. 45–48, 1986. View at: Google Scholar
 V. E. Henrich and P. A. Cox, The Surface Science of Metal Oxides, Cambridge University Press, Cambridge, Mass, USA, 1994.
 G. Pacchioni, C. Sousa, F. Illas, P. S. Bagus, and F. Parmigiani, “Measures of ionicity of alkalineearth oxides from the analysis of ab initio cluster wave functions,” Physical Review B, vol. 48, no. 16, pp. 11573–11582, 1993. View at: Publisher Site  Google Scholar
 P. V. R. Schleyer, C. Maerker, A. Dransfeld, H. Jiao, and N. J. R. Van Eikema Hommes, “Nucleusindependent chemical shifts: a simple and efficient aromaticity probe,” Journal of the American Chemical Society, vol. 118, no. 26, pp. 6317–6318, 1996. View at: Publisher Site  Google Scholar
 R. Kakkar and C. Singh, “Theoretical study of the kojic acid structure in gas phase and aqueous solution,” International Journal of Quantum Chemistry, vol. 111, no. 15, pp. 4318–4329, 2011. View at: Google Scholar
 R. Kakkar, M. Bhandari, and R. Gaba, “Tautomeric transformations and reactivity of alloxan,” Computational and Theoretical Chemistry, vol. 986, pp. 14–24, 2012. View at: Publisher Site  Google Scholar
 G. W. Wang and H. Hattori, “Reaction of adsorbed carbon monoxide with hydrogen on magnesium oxide,” Journal of the Chemical Society, Faraday Transactions 1, vol. 80, pp. 1039–1047, 1984. View at: Publisher Site  Google Scholar
 G. Bilalbegović, “Structural and electronic properties of MgO nanotube clusters,” Physical Review B, vol. 70, no. 4, Article ID 045407, 6 pages, 2004. View at: Google Scholar
Copyright
Copyright © 2013 Prinka Batra et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.