Quantum Chemistry Study on the Structures and Electronic Properties of Bimetallic Ca2-Doped Magnesium Ca2Mgn (n = 1-15) Clusters.

Here, by utilizing crystal structure analysis through the particle swarm optimization (CALYPSO) structural searching method with density functional theory (DFT), we investigate the systemic structures and electronic properties of Ca2Mgn (n = 1-15) clusters. Structural searches found that two Ca atoms prefer to occupy the external position of magnesium-doped systems at n = 2-14. Afterward, one Ca atom begins to move from the surface into the internal of the caged skeleton at n = 15. Calculations of the average binding energy, second-order difference of energies, and HOMO-LUMO gaps indicated that the pagoda construction Ca2Mg8 (as the magic cluster) has higher stability. In addition, the simulated IR and Raman spectra can provide theoretical guidance for future experimental and theoretical investigation. Last, further electronic properties were determined, including the charge transfer, density of states (DOS) and bonding characteristics. We hope that our work will provide theoretical and experimental guidance for developing magnesium-based nanomaterials in the future.

1. Introduction

Nanomaterials with small particle sizes, specific surface areas and high surface energies possess wide applications in chemistry, physics, biology, medicine, materials and nanodevices. Magnesium atoms contain s2 closed-shell electron configuration similarly to helium, which plays an essential role in aerospace, mobile electronics, automobile and biomedical applications [1,2,3,4]. Currently, metal-doped magnesium clusters with unique geometries and fascinating electronic properties have received considerable attention in magnesium-based multi-function materials.

In the past decades, many experimental techniques and theoretical studies have been reported for the related structures and properties of pure magnesium clusters [5,6,7,8,9,10,11,12,13,14,15,16,17,18,19]. For instance, on the experimental side, the transition points of Mg clusters were determined at n = 20 by Diederich’s group [5]. By measuring the photoelectron spectra and observing the sp band gaps of Mg (n = 3–35) clusters, Thomas et al. found that the anion magnesium clusters exhibited a metallic character from n = 18 [6]. On the theoretical side, the ground state van der Waals potential of a magnesium dimer was described by five essential parameters of the Tang–Toennies potential model [7].

Janecek et al. studied the structures of neutral and cationic Mg (n to 30) clusters using the local spin density functional [8]. Based on the MP4 (SDTQ) and CCSD(T) levels, the electron affinities were calculated for magnesium dimers and trimers [9]. By using the spin-unrestricted density functional theory with a local density approximation, Gong and co-workers studied the electronic structures of Mg (n < 57) clusters [10]. Their results found that more 3p electrons will be involved in the sp hybridization with the cluster size increasing. More recently, Akola and co-workers focused on the structural and electronic properties of Mg (n < 13) clusters [11]. Their investigations showed that the onset of metallization of Mg clusters is difficult to assign due to the energy gap and sp hybridization.

Subsequently, electron binding energies, structural and electronic properties and the nonmetal-to-metal transition were studied systemically by Jellinek and Acioli [12,13]. In addition, by using ab initio theoretical methods (B3PW91, B3LYP and MP4), Lyalin et al. investigated the structure and electronic properties of neutral and singly charged Mg (n = 2–21) clusters [14]. The above study by Lyalin et al. suggested that the hexagonal ring structure determines the cluster growth from Mg15. Moreover, the electronic shell effects and the Jellium-like behavior manifest themselves in the formation of geometrical properties; however, the shell effects do not determine the geometry of the Mg clusters completely.

In addition, the energetic structural properties of neutral magnesium clusters Mg (n = 2–22) were investigated utilizing density functional theory [15]. More recently, Heidari and Xia’s groups presented further details about the structural transition and electronic properties of Mg (n = 10–56) and Mg0/− (n = 3–20) clusters, respectively [16,17]. Recently, Duanmu’s group investigated the geometries, binding energies, adiabatic ionization potentials and adiabatic electron affinities of Mg+/0/− (n = 1–7) clusters [18] and cohesive energies by using the CCSD(T) scheme with MP2/CBS correction for Mg (n = 2–19 and 28), respectively [19].

Most recently, research found that heteroatom doping is an effective strategy to stabilize geometrical structures or to tune electronic properties. Up to now, based on the different quantum chemistry calculations, studies on metal-doped (Be, Al, Ge, Sn, 3d and 4d TM atoms) and nonmetal-doped (B, C, N, O, F and Si) magnesium clusters have harvested many great achievements [20,21,22,23,24,25,26,27,28,29,30,31,32]. More importantly, based on the CALYPSO structural searching method and DFT, the structures and electronic properties of Be-, Be2-, Sr2- and Ba2-atom-doped differently sized magnesium clusters have been systemically discussed by Zeng’s, Zhao’s and our groups [33,34,35,36,37]. For example, Zeng and co-workers found that BeMg90, BeMg9+1 and BeMg8−1 clusters possess relative highly stability out of the studied systems [33].

Moreover, the stability mainly originating from the σ-type covalent bond, is formed by the interaction between Be-s and Mg-p orbitals. Zhao’s group performed the structural and electronic properties of BeMg (n = 10–20) clusters and their anions [34]. The research concluded that the position of Be atom changes from completely encapsulated sites to surface sites after reverting to the caged magnesium motif. Subsequently, the structures and electronic properties have been investigated for two Be-, Sr- and Ba-atom-doped small-sized magnesium clusters by our group [35,36,37]. For Be2Mg (n = 1–20) clusters, from n = 10, the structures transfer from 3D to filled cage-like frameworks [35].

Furthermore, in the small size, one Be atom prefers the surface sites, and the other Be atom tends to embed inside magnesium motif. However, for the large size clusters (n > 18), two Be atoms were completely encapsulated into magnesium cages. In addition, the Be2Mg8 cluster possesses robust stability, and Be-2p and Mg-3p orbitals revealed increasing metallic behavior. More interestingly, based on the same calculated method, studies on the structural evolution and electronic properties were performed for SrMg0/− (n = 1–12) clusters [36]. As a result, the tower-like framework of the Sr2Mg8 cluster possesses higher stability out of the studied systems.

Moreover, the stronger sp hybridization leads to stronger Sr-Mg bonds, which is supported by the analysis of the multi-center bonds. Subsequently, our groups systemically reported the structures and electronic properties of two-barium-atom-doped magnesium in both neutral and anionic species [37]. A pagoda-like Ba2Mg8 was determined by analyzing the relative stability.

Analysis of the molecular orbitals indicated that the high stability comes mainly from the interaction between Ba-6s and Mg-3p orbitals. In conclusion, the stronger sp hybridization leads to stronger M-Mg bonds and metallic behavior. Last, the geometric structures and electronic properties have been investigated systemically for lithium-doped magnesium clusters [38]. The results indicated that lithium atoms prefer to occupy the convex sites of LiMg structures. The LiMg9 cluster possesses relatively higher stability. In addition, the charges transfer from the Li to Mg atoms, and there exists strong hybridization among sp orbitals.

Thus, metal-atom-doped magnesium clusters provide an effective approach to creating novel structures and electronic properties. As the same group of alkaline earth metals, Ca and Be, Mg, Sr and Ba have the identical valence electronic configuration of ns2. Nevertheless, due to the different electronegativity and atomic radius of Ca atoms compared with Be, Mg, Sr and Ba atoms, are there similar frameworks for two-calcium-atom-doped magnesium clusters? If yes, do these clusters possess novel electronic and bonding properties? How does hybridization change?

To date, minor investigations on the structures and electronic properties of two-calcium-atom-doped magnesium clusters have been reported in theoretical calculations and experimental works. Thus, in the present work, motivated by Be-, Mg-, Sr- and Ba-doped magnesium, we performed a systematic investigation for two-calcium-atom-doped magnesium Ca2Mg (n = 1–15) clusters.

First, we conducted wide structural searching and precise structural optimization to explore the structural evolution rule. Second, determining the stable configuration of Ca2Mg (n = 1–15) clusters was conducted by analyzing the stability properties. Finally, some electronic properties, such as the charge transfer, IR and Raman spectra, DOS and bonding characteristics, are discussed for Mg-doped alkaline-earth clusters. We hope that our investigations will provide a theoretical and experimental basis for studying the microscopic mechanism of magnesium doped with alkaline-earth nanomaterials.

2. Computational Detail

In the present section, a stochastic global search algorithm based on the CALYPSO structural searching method was used to obtain the lowest and low-lying isomers of Ca2Mg (n = 1–15) clusters and Mg+2 (n = 1–15) clusters [39,40,41]. Based on this method, the corresponding stable or metastable structures can be successful with the chemical composition and given external conditions [42,43,44,45,46,47,48,49,50], and the detailed searching process can be found in our reported papers [42,43,44,45,46,47]. The structural optimization and energy calculations employed the B3PW91 functional and 6−311+G(d) basis set for Ca and Mg atoms, respectively [51,52].

In particular, the B3PW91 functional has been widely tested for magnesium and magnesium-based clusters [35,36]. In addition, the spin multiplicity (1, 3, 5 and 7) is included, no imaginary frequencies are validated. All calculations were performed using the Gaussian09 program package [53]. In the following works, the relative stabilities of the ground state Ca2Mg (n = 1–15) and Mg+2 (n = 1–15) clusters were studied by computing the average binding energy (E) and second-order difference energy (). Subsequently, the IR and Raman spectra, DOS, molecular orbitals and AdNDP were systemically investigated using the Multiwfn software for the studied clusters [54]. To ensure the reliability of our computational method, the bond length (re), vibration frequencies (ωe) and dissociation energies (De) were calculated for Mg2, Ca2 and CaMg dimers, respectively.

For the Mg2 dimer, our calculated values were re = 3.651 Å, De = 0.0790 eV, which are in good agreement with the experimental results (3.891 Å, 0.0866 eV) [55]. For the Ca2 and CaMg dimers, there are no experimental values available. Our calculated results for the bond lengths, vibration frequencies and dissociation energies are 4.2667 Å, 72.32 cm−1 and 0.1478 V for the Ca2 dimer and 3.909 Å, 82.63 cm−1 and 0.1114 V for the CaMg dimer, respectively. Moreover, the bond length and frequency of Ca2 dimer are also in excellent agree with Soltani’s theoretical values (4.285 Å and 65.2 cm−1), respectively [56]. This indicated the reliability of the proposed method in this work.

3. Results and Discussions

3.1. Geometric Structures

In Figure 1, the lowest and lower-lying energy structures are present for Ca2Mg (n = 1–15) and Mg+2 (n = 1–15) clusters. These isomers with energies from low to high are designated by na, nb, nc and nd. Moreover, the electronic states and point symmetry are also collected in Table 1. Simultaneously, the Cartesian coordinates of the lowest energy structures of Ca2Mg (n = 1–15) clusters are given in Table S1 in the Supplementary Materials.

Figure 1

Optimized geometrical structures of lowest and low-lying isomers of Ca2Mg (n = 1–15) clusters and Mg+2 (n = 1–15) clusters along with the point group symmetry, electronic states and relative energy (eV). The pinkish and orange balls are magnesium and calcium atoms, respectively.

Table 1

Electronic states, symmetries, average binding energies E (eV), HOMO–LUMO energy gaps E (eV), the second-order difference energy (∆2E) and charges on the Ca atoms of the most stable Ca2Mg (n = 1–15) clusters.

Clusters	State	Sym.	Eb (eV)	Δ2E (eV)	Eg (eV)	Charge (e)
						Ca1	Ca2
Ca2Mg	1A1	C2V	0.20	-	1.91	0.11	0.11
Ca2Mg2	1A1	C2V	0.40	0.57	2.01	0.31	0.31
Ca2Mg3	1A1	C2V	0.40	−0.21	1.85	0.37	0.37
Ca2Mg4	1A1	C2V	0.44	−0.25	1.71	0.41	0.41
Ca2Mg5	1A1’	D5H	0.50	0.33	1.39	0.53	0.53
Ca2Mg6	1A	D5H	0.51	−0.37	1.51	0.37	0.37
Ca2Mg7	1A1	C2V	0.56	−0.19	1.31	0.68	0.68
Ca2Mg8	1A’	CS	0.61	0.66	1.77	0.84	0.94
Ca2Mg9	1A	C2	0.60	−0.20	1.77	0.71	0.71
Ca2Mg10	1A	C1	0.60	−0.19	1.16	0.86	0.71
Ca2Mg11	1A	CS	0.62	0.14	1.13	0.87	0.71
Ca2Mg12	1A’	CS	0.63	−0.42	1.12	0.85	0.85
Ca2Mg13	1A’	CS	0.66	0.06	1.19	0.86	0.86
Ca2Mg14	1A’	CS	0.69	−0.06	0.93	1.06	1.06
Ca2Mg15	1A1	C3V	0.71	-	1.41	−0.18	0.85

First, the lowest energy structures of the Mg clusters agree with previous research by Zhang and Li et al. [35,37], which indicates that the present functional and basis sets are reliable. Moreover, the Ca2Mg (n = 1–8) clusters possess similar geometric structures compared with the Mg+2 clusters. Second, Ca2Mg possesses a triangular plane structure, and the lowest energy structures of Ca2Mg (n = 2–15) clusters maintain a three-dimensional (3D) configuration.

Specifically, Ca2Mg (n = 2–6) clusters can be obtained by adding one Mg atom to the Ca2Mg−1 clusters, and Ca2Mg (n = 7–9) clusters can be generated from the substitution of Mg+2 clusters by two Ca atoms, respectively. Interestingly, the lowest energy structure of the Ca2Mg8 cluster possesses the same geometrical form as those of X2Mg8 (X = Be, Sr and Ba) clusters. Finally, from 2 to 14, the doped two Ca atoms prefer to locate outside the host Mg+2 cluster. At n = 15, one Ca atom starts to move from the surface into the internal of the caged skeleton.

3.2. Relative Stability

To determine the relative stabilities of Ca2Mg and Mg+2 (n = 1–15) clusters, the average binding energy (E) and second-order difference of energy are calculated as follows:

E denotes the total energy of the corresponding clusters or atoms. The calculated results are plotted in Table 1 and Figure 2A,B. The following information can be concluded: (1) The values of E(Mg+2) clusters are lower than those of E(Ca2Mg) clusters, indicating that the doped Ca atom can improve the stability of pure magnesium clusters. (2) With an increase of Mg atoms, the values of E (Mg+2 and Ca2Mg) monotonically increase, which indicates an enhanced effect on the stabilities for Ca2Mg and Mg+2 (n = 1–15) clusters. (3) values show a certain degree of oscillation, and there exist four visible peaks in the curves at n = 2, 5, 8, 11 and 13, implying that the Ca2Mg2,5,8,11,13 and Mg4,7,10,13,15 clusters possess higher relative stability. Most interestingly, the Ca2Mg8 cluster with 20 valence electrons possesses an notably high value (0.66 eV), which indicates superior relative stability.

Figure 2

(A) E, (B) Δ2E, (C) E and (D) the total charges on Ca atoms in the ground state of Ca2Mg (n = 1–15) clusters.

The HOMO–LUMO energy gap (Eg) is the other powerful tool to study the relative stabilities. Generally speaking, large values indicate stronger stability. In the present section, the HOMO–LUMO energy gaps are presented in Table 2 and Figure 2C. First, values of Eg(Mg+2) are larger than those of Eg(Ca2Mg) clusters, indicating that Ca2Mg clusters are more stable, which is in reasonable agreement with previous research on averaged binding energies.

Table 2

The HOMO and LUMO energy of Ca2Mg and Mg+2 (n = 1–15).

Clusters	HOMO (eV)	LUMO (eV)	Clusters	HOMO (eV)	LUMO (eV)
Ca2Mg	−3.834	−1.920	Mg3	−4.799	−1.937
Ca2Mg2	−3.981	−1.973	Mg4	−4.961	−2.076
Ca2Mg3	−3.831	−1.979	Mg5	−4.258	−2.099
Ca2Mg4	−3.860	−2.148	Mg6	−4.294	−2.311
Ca2Mg5	−3.531	−2.140	Mg7	−4.436	−2.347
Ca2Mg6	−3.897	−2.384	Mg8	−4.284	−2.507
Ca2Mg7	−3.917	−2.611	Mg9	−4.393	−2.868
Ca2Mg8	−3.954	−2.188	Mg10	−4.306	−2.350
Ca2Mg9	−4.048	−2.275	Mg11	−4.385	−2.338
Ca2Mg10	−3.750	−2.588	Mg12	−4.012	−2.583
Ca2Mg11	−3.671	−2.543	Mg13	−4.105	−2.652
Ca2Mg12	−3.746	−2.623	Mg14	−4.092	−2.881
Ca2Mg13	−3.730	−2.539	Mg15	−3.973	−2.750
Ca2Mg14	−3.518	−2.592	Mg16	−3.985	−2.521
Ca2Mg15	−3.965	−2.558	Mg17	−3.917	−2.823

Second, Ca2Mg2,6,8,9,15 clusters with local maxima of Egap suggests that those clusters are more stable than their neighbors. In summary, combining the conclusions of E, Δ2E and E values, the Ca2Mg8 cluster corresponds to the magic numbers and exhibits robust stability.

3.3. Charge Transfer

In this section, the charge-transfer information is analyzed by natural population analysis (NPA) in Table 1 and Figure 2D. Clearly, the doped Ca atoms possess positive charges in doped systems, meaning the charges transfer from calcium to magnesium atoms.

Thus, the role of magnesium atoms is as charge acceptors, and Ca atoms are the charge donors. This is expected as Ca (1.00) has a smaller electronegativity compared with that of Mg (1.31) [57]. Second, the transferred charges increase with increasing cluster size. However, the value falls sharply at n = 15, and this situation may be the result of the position of Ca atom in the caged skeleton. In addition, the transfer charges in the range of n = 5–14 are all greater than 1.0 eV except for Ca2Mg6.

3.4. Infrared (IR) and Raman Spectra

In order to facilitate the characterization of spectra, we computationally simulated the infrared (IR) and Raman spectra of the lowest energy structures of Ca2Mg8 cluster. The simulated spectra with atomic labels are shown in Figure 3. The results found the highest intense IR frequency was located at 197.74 cm−1 with 4Mg-7Mg-10Mg bond tensile vibration; however, two Ca atoms were almost silent. The second and third strongest peaks can be found at 129.62 and 218.78 cm−1. For Raman spectra, the strongest peak at 181.51 cm−1 corresponds to the breathing vibration of all atoms. The second- and third-strongest Raman frequencies at 138.28 and 129.62 cm−1 correspond to the swing vibration of all atoms. In addition, IR and Raman spectra revealed that the strongest spectral frequencies are displayed in the range of 100–200 cm−1.

Figure 3

The Infrared and Raman spectra of the most stable cluster of the Ca2Mg8 cluster.

3.5. The Density of States

To understand the nature of the chemical bonding, the total density of states (TDOS) and partial density of states (PDOS) of the lowest energy structures of Ca2Mg (n = 1–15) are displayed in Figure 4. The TDOS is represented by the khaki shade; PDOS of Ca-s, Ca-p, Mg-s and Mg-p AOs (atomic orbitals) are represented by the red and blue solid lines as well as magenta and green dotted lines, respectively. We found that the contribution to TDOS mainly comes from the PDOS of Ca-s, Mg-s and Mg-p AOs in the region of occupied orbitals.

Figure 4

The calculated total densities of states (TDOS) and partial densities of states (PDOS) of Ca2Mg (n = 1–15) clusters.

This indicates that sp hybridization has occurred in Ca-Mg atoms and Mg-Mg atoms. In fact, the sp hybridization of Mg-Mg promotes the formation of Mg frames in Ca2Mg clusters. The sp hybridization of Ca-Mg can promote the interaction between the two doped calcium atoms and magnesium frames of Ca2Mg clusters, which is also the main reason why the stability of Ca2Mg is higher than that of their corresponding pure magnesium clusters.

3.6. Bonding Characters

Based on the above analyses, a pagoda-like Ca2Mg8 structure possesses superior stability. To illustrate the source of higher stability, the bonding nature, such as the MOs (molecular orbitals) and multi-center bonds, are discussed for the lowest energy structure of Ca2Mg8 cluster. The molecular energy levels and corresponding orbitals are presented in Figure 5. First, calcium and magnesium atoms are composed of the same valence configuration of s2, and Ca2Mg8 with 20 valence electrons meets the requirement of the Jellium model in terms of the valence electron number.

Figure 5

Molecular orbitals and the corresponding energy levels of the Ca2Mg8 cluster. The HOMO–LUMO gap is indicated (in azure).

Moreover, the shell structures consist of one 1S orbital, three 1P orbital, five 1D orbitals and one 2S orbital, all of which are occupied by the paired electrons. The energy of 1S, 1P, 1D and 2S states are arranged in order from low to high without energy levels overlapping, and the energy levels are also relatively concentrated. In addition, all the splitting energy levels of 1D orbital are lower in energy than the 2S orbital. As a result, Ca2Mg8 cluster is a closed shell 1S21P61D102S2 filled with 20 valence electrons.

In addition, the contributions of molecular orbital for Ca2Mg8 cluster were probed utilizing Multiwfn 3.8 program. The HOMO corresponding to 2S state involves Ca-s (30.98%), Mg-s (19.53%) and Mg-p (44.86%). The HOMO-m (m = 1–5) exhibit 1D state, in which HOMO-1 and HOMO-4 are mainly composed of Ca-s (10.62%, 26.60%), Mg-s (41.87%, 16.99% for) and Mg-p (43.70%, 53.33%), respectively, and the remaining 1D state orbital is mainly contributed to by Mg atoms. The compositions of HOMO-6, 7, 8 corresponding to three 1P orbitals comes mainly from the Mg-s and p AOs as well as to a small extent from the Ca-s and p AOs. In the case of HOMO-9 (1S), Mg atoms provide more than 90% of the orbital contribution. Hence, the sp hybridization between the Ca and Mg atoms could promote the interaction between the doped-Ca and host-Mg atoms and form stronger Ca-Mg bonds.

To further elucidate the bonding patterns of the Ca2Mg8 cluster, adaptive natural density partitioning (AdNDP) is included, representing the bonding of a molecule in terms of n-center two-electron (nc-2e) bonds. Notably, 1c-2e and 2c-2e bonds mean a localized character; the nc-2e bond (n ≥ 3) belongs to the delocalized character. For the Ca2Mg8 cluster, ten delocalized bonds with different occupation number (ON) are shown in Figure 6, which corresponds to three delocalized 3c-2e σ-bonds, two delocalized 4c-2e σ-bonds, two delocalized 5c-2e π-bonds, two delocalized 7c-2e π-bonds and one delocalized 9c-2e π-bond, respectively.

Figure 6

AdNDP chemical bonds and the corresponding structural units for the Ca2Mg8 cluster. (ON denotes the occupation number).

Three 3c-2e σ-bonds with an ON of 1.745–1.766 |e| are formed by two opposite trigonal 5Mg-8Mg-10Mg and 3Mg-6Mg-7Mg units and one 2Ca-5Mg-3Mg unit. Two 4c-2e σ-bonds are accountable for the bonding between the quadrangle 4Mg-6Mg-8Mg-9Mg units (ON = 1.89 |e|) and pyramid 1Ca-7Mg-9Mg-10Mg units (ON = 1.774 |e|), respectively. There are two 5c-2e π-type bonds with ON = 1.80 and 1.792 |e|, derived from two Ca-Mg units of 1Ca-2Ca-3Mg-5Mg-9Mg and 1Ca-3Mg-4Mg-5Mg-9Mg units, respectively. Two 7c-2e ON = 1.945, 1.903 |e| (two symmetric units) and one 9c-2e ON = 1.994 |e| exhibit π-type bonding character. In addition, to deeply understand the nature of bonding, the Wiberg bond orders are calculated and listed in Table 3. The results indicate that the strong Ca-Mg bonds are greater than Ca-Ca and most Mg-Mg bonds.

Table 3

The Wiberg bond orders of the Ca2Mg8 cluster.

Atoms	Ca-1	Ca-2	Mg-3	Mg-4	Mg-5	Mg-6	Mg-7	Mg-8	Mg-9
Ca-2	0.130
Mg-3	0.383	0.462
Mg-4	0.151	0.143	0.233
Mg-5	0.383	0.462	0.524	0.233
Mg-6	0.139	0.131	0.517	0.529	0.228
Mg-7	0.398	0.149	0.528	0.439	0.237	0.523
Mg-8	0.139	0.131	0.228	0.529	0.517	0.182	0.199
Mg-9	0.179	0.410	0.465	0.458	0.465	0.487	0.231	0.487
Mg-10	0.398	0.149	0.237	0.439	0.528	0.199	0.523	0.523	0.231

4. Conclusions

In summary, a detailed investigation of the structures and electronic properties of Ca2Mg (n = 1–15) was performed using the CALYPSO searching method and DFT calculations. Structural searching found that two Ca atoms prefer to occupy the external position of magnesium-doped systems at n = 2–14 and that one Ca atom tends to move from the surface into the internal of the caged skeleton at n = 15. The size-dependent binding energies, second-order difference of energies, and HOMO–LUMO gaps found a pagoda-like Ca2Mg8 as a magic cluster that possessed higher stability.

Upon charge transfer analysis, charges transferred from calcium to magnesium atoms. The simulated IR and Raman spectra of the magic cluster revealed that the strongest spectral frequencies were displayed in the range of 100~200 cm−1. In addition, the high stability of Ca2Mg8 with a 20 valence electron cluster possessed a closed-shell electron configuration of 1S21P61D102S2 in terms of the Jellium model. Last, the sp hybridization of Ca-Mg and Mg-Mg bonds was confirmed by the molecular orbitals and AdNDP, which contribute to the high stability of the Ca2Mg8 cluster.