Nguyen Minh Tam*ab,
Long Van Duongc,
Ngo Tuan Cuongd and
Minh Tho Nguyene
aComputational Chemistry Research Group, Ton Duc Thang University, Ho Chi Minh City, Vietnam. E-mail: nguyenminhtam@tdtu.edu.vn
bFaculty of Applied Sciences, Ton Duc Thang University, Ho Chi Minh City, Vietnam
cInstitute for Computational Science and Technology (ICST), Quang Trung Software City, Ho Chi Minh City, Vietnam
dFaculty of Chemistry and Center for Computational Science, Hanoi National University of Education, Hanoi, Vietnam
eDepartment of Chemistry, KU Leuven, Celestijnenlaan 200F, B-3001 Leuven, Belgium
First published on 30th August 2019
Structures of the binary AlnSim clusters in both neutral and cationic states were investigated using DFT and TD-DFT (B3LYP/6-311+G(d)) and (U)CCSD(T)/cc-pvTZ calculations. Silicon-doped aluminum clusters are characterized by low spin ground states. For small sizes, the Si dopant prefers to be located at vertices having many edges. For larger sizes, the Si atom prefers to be endohedrally doped inside an Aln cage. Relative stability, adiabatic ionization energy and dissociation energies of each cluster size were evaluated. A characteristic of most Si doped Al clusters is the energetic degeneracy of two lowest-lying isomers. Calculated results confirm the high stability of the sizes Al4Si2, Al12Si and Al11Si2+ as “magic” clusters, that exhibit 20 or 40 shell electrons and are thermodynamically more stable as compared to their neighbors. Electronic absorption spectra of isoelectronic magic clusters Al13−, Al12Si, and Al11Si2+ that have two pronounced bands corresponding to blue and violet lights, have been rationalized by using the electron shell model. The magnetically included ring current density (MICD) analyses suggest that they are also aromatic structures as a result of the “magic” 40 shell electrons.
Along with pure Al clusters, several doped Al clusters have also been investigated. For alkali and alkali-earth dopants, while Rao et al.18 presented a theoretical study on their geometries and energetics including adsorption energies, ionization potentials, and electron affinities of AlnLi and AlnK, Sun and co-workers19 examined AlnBe using the high accuracy method CCSD(T). Wang and coworkers20,21 investigated the Al Zintl anion moieties within the mixed sodium as well as magnesium aluminum clusters NamAln−/0 and MgmAln− using a combination of anion photoelectron spectroscopy and density functional theory (DFT) calculations. Studies on transition metal (TM) doped Al clusters provide us with a considerable amount of results. For copper derivatives AlnCu with n = 11–14, unlike the alkali atom-doped counterparts in which the alkali dopant prefers an exohedral position, the copper atom tends to reside inside an Al cage.22 For platinum doped aluminum clusters, Zhang et al. reported a joint photoelectron spectroscopic and theoretical study of the PtAl− and PtAl2− anions.23 An experimental study on the argon physisorption ability of TM indicated that AlnTM+ clusters attach one argon atom up to a critical cluster size ncrit, with ncrit = 16 and ncrit = 19–21 for TM = V, Cr, and TM = Ti, respectively, and undergo a geometrical transition in going from exohedrally to endohedrally doped clusters in which the TM atom is located inside an Al cage.24 Experimental results on AlnTi were confirmed by subsequent theoretical DFT computations.25
More recently, the influence of spin on the properties of TM doped aluminum clusters in a series of low spin clusters was found in which a prominent odd–even oscillation in all calculated properties supports the presence of jellium shell structure. The electron jellium model considers a cluster of atoms as a superatom where the cores of constituent atoms form a constant positive background, the “jellium density”, and only valence electrons are treated explicitly. Under the model, the shell molecular orbitals of the clusters are denoted as 1S, 1P, 1D, 2S, 2P, 1F, … (usually the notation 1s, 1p, 1d, 2s, 2p, 1f, … have also been used), rather than as the atomic orbitals 1s, 2s, 2p, 3s, 3p, 3d, … of hydrogen-like atoms. High spin ground state isomers show more smooth trends similar to the bulk materials.26 Geometrical and electronic structures as well as physical and chemical properties of Al clusters doped by various common nonmetal dopants have also abundantly been investigated.27–45
Interaction of hydrogen molecules with AlnRh2+ with n = 10–13 clusters was recently studied by mass spectrometry and infrared multiple photon dissociation (IRMPD) spectroscopy.46 A comparison of the IRMPD spectra with predictions obtained using DFT calculations showed that for n = 10 and 11, a single H2 molecule dissociates upon binding to the cluster, whereas for n = 12 and 13, it adsorbs molecularly. Upon adsorption of a second H2 molecule to the cation Al12Rh2+, both hydrogen molecules dissociate. These results point out a peculiar catalytic effect of doped Al clusters.
Silicon, being contiguous to Al in the periodic table and having one valence electron more than Al, has been and still is massively used in the semiconductors and optoelectronic industries. A mixture of aluminum–silicon atoms gives rise to some quite interesting compounds as previous studies indicated some remarkable properties of Al–Si nanomaterials. For instance, interaction between a Si nanowire and one Al-atom was found to increase the electrical conductivity as compared to pristine Si-nanowire.47 Al atoms form an ordered array of magic clusters on the surfaces of Si(111)48 and a persistent formation of Al–Si nanowires was reported.49 However, understanding of geometrical and electronic structures and physico-chemical properties of binary Al–Si clusters remains limited. Earlier studies mainly examined the Al doped Si clusters.50–52 whereas only a few studies on the Si doped Al clusters were performed and concentrated on high symmetry structures53,54 that are not necessarily the most thermodynamically stable structural motifs.
In view of the scarcity of reliable results on mixed Al–Si clusters, we set out to carry out a systematic investigation, using both density functional theory (DFT) and wavefunction methods, on a series of small singly and doubly silicon doped aluminum cluster AlnSim, with n = 3–16 and m = 1–2, in both neutral and cationic states. Geometries of the most stable equilibrium structures and their energetic parameters are determined by using both DFT and coupled-cluster theory CCSD(T) methods. It is of importance to emphasize the similarities and differences between the behavior of Al–Si clusters that have a similar size and number of electrons. As for a typical example, as the isoelectronic species Al13−, Al12Si and Al11Si2+ have closed electronic shells in the jellium model (40 valence electrons), we thus model their optical absorption spectra by using the time-dependent DFT (TD-DFT).
Geometries and harmonic vibrational frequencies of the lower-lying isomers having relative energies of <10 eV of SinAlm at the small sizes (n ≤ 4) and their corresponding cations are then reoptimized with the higher spin states (triplet and quintet for closed shell clusters, and quartet and sextet for open shell clusters). However, no high spin states having a remarkable stability are found, due to the fact that both Al and Si are basically non-magnetic atoms. It can thus be predicted that the larger sizes of AlnSim also exist as non-magnetic clusters, and therefore we can neglect the calculations of their high spin states due to the limitations of our computational resource.
The local minima with relative energies of <5 eV with respect to the corresponding lowest-lying isomer of each size are considered again and reoptimized using the same B3LYP functional but with the larger 6-311+G(d) basis set. Harmonic vibrational frequencies are again calculated at this level in order to identify the obtained structures as local minima on the potential energy surface and to evaluate zero-point energy corrections. To obtain more reliable energetic parameters, the lower-lying isomers of each size having relative energies within 1.0 eV are selected for single point electronic energies using the coupled-cluster theory CCSD(T) with the correlation consistent cc-pVTZ basis set, based on B3LYP optimized geometries. The unrestricted formalism (UHF, UB3LYP, UCCSD) is employed for open-shell structures.
In an attempt to probe further the effects induced by the Si dopants on the electronic structure, we perform time-dependent density functional theory (TD-DFT) calculations using the same functional/basis set and the optimized geometries of the closed-shell isoelectronic Al13−, Al12Si and Al11Si2+ species. These calculations also allow us to construct their absorption spectra. For each species, the number of electronically excited states considered is large enough for the spectra spread out up to a wavelength of 250 nm. The shape of molecular orbitals involved in the electron shells and the density of states of the above clusters are also plotted assisting the analysis of their spectra.
In addition, we use the magnetically included ring current density (MICD) approach which is computed by integrating the electronic current passing through interatomic surfaces located by the atom-in-molecule theory (QTAIM) between neighboring atoms by using the AIMALL suite of programs.56 The purpose of the latter computations is to assess the aromatic character of the Al13−, Al12Si, and Al11Si2+ species. The magnetic current density derived from the MICD method is performed with the magnetic field vector perpendicular the molecular plane and out of that plane. As for a convention, clockwise flow of the magnetic current stands for a diamagnetic current, and anticlockwise flow for a paramagnetic current.
As for a convention, label X.n.m.Y is used to denote the isomers in which X = N and C stands for a neutral and cationic state, respectively, n is the number of Al atoms, m is the number of Si atoms, and Y = A, B, C… referring to the different isomers with increasing relative energy. Accordingly, the notation X.n.m.A invariably refers to the most stable isomer A of the X.n.m series.
The corresponding cations of N.6.1.A, N.7.1.A, and N.8.1.A remain also the lowest-lying isomers of the cationic counterparts. Again, calculated results show that the two most stable isomers exist as energetically degenerate for both cations Al6Si+ and Al7Si+ with energy gap of only 0.03 and 0.01 eV, respectively. For the cation Al9Si+, C.9.1.B which is the corresponding cation of N.9.1.A, now lies 0.17 eV higher in energy than the C2v symmetric structure C.9.1.A although N.9.1.C is the corresponding neutral of C.9.1.A and is 0.40 eV less stable than N.9.1.A.
For the neutral Al14Si, calculations continue to emphasize an energetic degeneracy with an energy gap of 0.03 eV between the high symmetry structure N.14.1.A (D2h) and N.14.1.B in which a Si and an Al atom are exohedrally attached on the icosahedron Al13. In other words, it can be considered that N.14.1.B is formed by substituting a Si atom at an Al-top of the high symmetrical (D3d) anion Al15− and this may lead to the stability of N.14.1.B.12 The AlnSi can be regarded to be generated upon substitution of the central Al position of the icosahedron Al13 of the Aln+1 species, with n = 12–16, by a Si atom.12
Generally, since the Si element has one valence electron more than the Al, the neutral Si is isoelectronic with the anion Al−. Therefore, it is observed in Fig. 1 and 2 that in small neutral AlnSi clusters (n ≤ 9), the Si atom prefers a substitution at an Al position, which exhibits a higher coordination number, of the pure aluminum framework of the isoelectronic anion Aln+1−. The neutral Al10Si represents however an exception because its most stable structure significantly differs from the Al11− framework. The ground state structure of a cationic AlnSi+ cluster turns out to be also quite similar to those of the pure neutral Aln+1. Moreover, for the cluster AlnSi with n ≥ 10, the Si dopant prefers to be located at the inner cage formed by Al atoms. According to Weigend et al.,57 this can be rationalized by the fact that the Si atom has larger nuclear charge as compared to Al atom, and its valence orbitals (3s, 3p) are lower-lying in energy, more compact to its nucleus as compared to Al atom, causing a higher electron density at the center as compared to that in the surface of the AlnSi cluster.
Similar to the growth pattern of singly silicon doped AlnSi0/+, it can be considered that the structures of doubly silicon doped AlnSi20/+ are obtained following substitution at two Al positions of the pure corresponding aluminum framework Aln+2 by two Si atoms. For n = 12–16, AlnSi20/+ are built up on the basis of an icosahedral growth path in which one Si atom is situated at the centre of the icosahedral Al12Si, whereas the additional Si atom is now externally attached.
Eb(AlnSi) = [nE(Al) + E(Si) − E(AlnSi)]/(n + 1) | (1) |
Eb(AlnSi+) = [(n − 1)E(Al) + E(Al+) + E(Si) − E(AlnSi+)]/(n + 1) | (2) |
Eb(AlnSi2) = [nE(Al) + 2xE(Si) − E(AlnSi2)]/(n + 2) | (3) |
Eb(AlnSi2+) = [(n − 1)E(Al) + E(Al+) + 2xE(Si) − E(AlnSi2+)]/(n + 2) | (4) |
Eb(Aln) = [nE(Al) − E(Aln)]/n | (5) |
Eb(Aln+) = [(n − 1)E(Al) + E(Al+) − E(Aln+)]/n | (6) |
The Eb value tends to rise with increasing cluster sizes. For singly doped AlnSi clusters, the neutral Al12Si reveals the largest Eb value as compared to the remaining singly doped neutral species. Of the neutral AlnSi2 clusters, Al4Si2 presents an enhanced Eb value as compared to those of the remaining doubly doped species. In the series of cations considered, Al11Si2+ gets a maximum peak in the Eb plot which demonstrates its higher thermodynamic stability.
The enhanced relative stability of these clusters can simply be rationalized in terms of “magic clusters” in which the number of valence electrons amounts to 2, 8, 18, 20, 34, 40 and 58… etc. Accordingly, the closed electron shells [1S2 1P6 1D10 2S2 1F14 2P6…] of the neutral Al4Si2 contains the same number of 20 valence electrons but an atom fewer in comparison to the pure cation Al7+, while both Al12Si and Al11Si2+ have each 40 valence electrons and the same icosahedral structure with a total number of 13 atoms such as the anion Al13−. As a result, they behave as “magic clusters” with an enhanced thermodynamic stability.
Fig. 4 also displays a comparison of the binding energies of the pure Aln, AlnSi, and AlnSi2 using (U)CCSD(T)/cc-pVTZ + ZPE calculated results in which the Eb value of AlnSi is larger than the Eb value of Aln and smaller than AlnSi2 at the same number of atoms. Accordingly, it appears that the Si dopant tends to increase the cluster stability with respect to various fragmentation processes.
Number of atoms, k | Alk−1Si | Alk−2Si2 | Alk |
---|---|---|---|
N.k − 1.1.A | N.k − 2.2.A | ||
4 | 6.75 | — | 6.43 |
5 | 6.67 | 6.39 | 6.36 |
6 | 6.80 | 7.10 | 6.69 |
7 | 6.30 | 6.17 | 5.68 |
8 | 5.94 | 6.40 | 6.42 |
9 | 6.55 | 6.16 | 6.09 |
10 | 5.83 | 6.31 | 6.49 |
11 | 6.59 | 5.71 | 6.13 |
12 | 5.88 | 6.43 | 6.23 |
13 | 6.85 | 5.46 | 5.63 |
14 | 5.60 | 5.88 | 5.84 |
15 | 5.89 | 5.84 | 5.66 |
16 | 5.71 | 6.03 | 5.97 |
17 | 5.79 | 5.36 | 5.63 |
18 | — | 5.87 | 6.05 |
n | De(1) remove a Si | De(2) remove a Al | De(3) remove a Si+ | De(4) remove a Al+ | ||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|
AlnSi | AlnSi+ | AlnSi2 | AlnSi2+ | AlnSi | AlnSi+ | AlnSi2 | AlnSi2+ | AlnSi+ | AlnSi2+ | AlnSi+ | AlnSi2+ | |
3 | 3.56 | 3.22 | 3.82 | 4.18 | n/a | n/a | n/a | n/a | 4.91 | 5.53 | n/a | n/a |
4 | 4.36 | 4.12 | 4.74 | 4.31 | 2.71 | 2.79 | 3.63 | 2.92 | 5.79 | 5.74 | 1.98 | 2.47 |
5 | 4.75 | 4.31 | 3.81 | 4.44 | 2.95 | 2.82 | 2.02 | 2.96 | 6.06 | 5.75 | 2.09 | 1.79 |
6 | 4.49 | 4.88 | 3.64 | 3.54 | 2.74 | 3.24 | 2.57 | 2.34 | 6.29 | 5.35 | 2.38 | 2.11 |
7 | 3.59 | 3.33 | 4.17 | 3.95 | 2.25 | 2.61 | 2.77 | 3.02 | 5.76 | 6.12 | 2.25 | 2.56 |
8 | 4.22 | 4.09 | 3.81 | 4.06 | 2.97 | 2.36 | 2.62 | 2.46 | 5.77 | 5.61 | 2.36 | 2.25 |
9 | 3.97 | 4.23 | 4.16 | 4.27 | 2.26 | 2.99 | 2.61 | 3.2 | 6.25 | 6.55 | 2.37 | 2.83 |
10 | 4.68 | 4.58 | 4.16 | 4.32 | 3.39 | 2.62 | 3.39 | 2.67 | 6.2 | 5.83 | 2.73 | 2.9 |
11 | 4.49 | 4.75 | 5.32 | 5.73 | 2.53 | 3.24 | 3.68 | 4.65 | 6.72 | 7.96 | 2.59 | 4.16 |
12 | 6.28 | 5.66 | 2.7 | 3.67 | 4.9 | 3.93 | 2.28 | 1.86 | 7.53 | 4.92 | 3.99 | 2.34 |
13 | 4.69 | 4.72 | 3.87 | 3.63 | 1.66 | 2.91 | 2.83 | 2.88 | 7.2 | 6.14 | 2.00 | 2.93 |
14 | 4 | 3.94 | 4.15 | 4.01 | 2.58 | 2.29 | 2.86 | 2.67 | 6.22 | 6.23 | 2.63 | 2.77 |
15 | 4.07 | 4.01 | 3.66 | 4.02 | 2.55 | 2.73 | 2.06 | 2.74 | 6.46 | 6.41 | 2.78 | 2.64 |
16 | 4.19 | 4.37 | 3.66 | 3.59 | 2.93 | 2.85 | 2.93 | 2.42 | 6.5 | 5.9 | 3.08 | 3.01 |
The De values of the neutrals AlnSi for the Si removal channel (1) AlnSi → Aln + Si turn out to be larger than those for the Al-loss channel (2) AlnSi → Aln−1Si + Al proving the stronger Si–Al bond. From the obtained knowledge in solid state, the loss of Al is a lower energy channel than the loss of Si, and this can be understood by the fact that the cohesive energy of Si (4.63 eV per atom) is larger than that of Al (3.39 eV per atom). Similar observations can be made for the positively charged species in that AlnSi+ tend to be fragmented generating one Al element or its cation Al+ plus a smaller Aln−1Si+/0 along the fragmentation channels (2) and (4).
Again, for doubly doped neutral AlnSi2 clusters, dissociation energies for Si-loss channels (1) AlnSi2 → AlnSi + Si are constantly larger than those for Al-elimination pathways (2) AlnSi2 → Aln−1Si2 + Al. Similarly, the cations AlnSi2+ follow a preferential fragmentation to form one Al0/+ plus a smaller Aln−1Si+/0 as described in the channels (2) and (4) (cf. Table 2).
Moreover, from the values given in Table 2, we can infer two exchange reactions:
Aln + Si → Aln−1Si + Al; | (R1) |
AlnSi + Si → Aln−1Si2 + Al | (R2) |
The energies the exchange reaction (R1) and (R2) can easily be calculated via the differences between De(1) and De(2) of AlnSi and AlnSi2. Generally, all energetic values of (R1) and (R2) are approximately in a range of 0.5–2 eV. Remarkably, the maximum value of 3.03 eV, obtained at n = 12 for (R1), corresponds to that of Al12Si. It proves the very particular case of neutral Al12Si for which the energy of (R1) is much higher than for any other case.
Of the whole series of the Aln−, Aln−1Si, and Aln−2Si2+ clusters, let us consider a typical series of high symmetry isoelectronic structures Al13−, Al12Si and Al11Si2+ to emphasize on how an assemblage of degenerate or nearly degenerate levels is manifested in their optical absorption spectra.
In order to investigate the electronic structure of Al13−, Al12Si, Al11Si2+, their MO diagrams are first plotted and displayed in Fig. 5. The electron configuration of the outer electron shell of each cluster could then be written as […1S2 1P6 1D10 2S2 2P6 1F14…], in which the high fold degeneracy levels such as the 1F, 1D, even the 1P and 2P could be split and mixed with other levels. To confirm that the number of 40 electrons in each cluster shell are arising from the valence electrons of Al atoms, we calculate the contributions of atomic orbitals (AO) of Al atoms to the shell MOs by using Natural Bonding Orbital (NBO) analysis.
The contributions of s, p and d AOs of Al atoms into the shell MOs of the anion Al13− are found as follows: 1S = 94% s + 5% p + 0.1% d; 1P = 81% s + 18% p + 0.22% d; 1D = 89% s + 10% p + 0.49% d; 2S = 97% s + 3% p + 0.25% d; 1F (three lower 1F shell orbitals) = 94% s + 5% p + 0.22% d; 2P = 10% s + 89% p + 0.74% d and 1F (four higher ones) = 0% s + 98% p + 1.6% d. This result shows that the lower shell orbitals 1S, 1P, 1D, 2S, 1F (three lower ones) are mainly constructed from the s(Al) AOs and the higher shell orbitals 2P and 1F (four higher orbitals) are mainly constructed from p(Al) AOs. In other words, the 26 s-electrons of 13 Al atoms form the lower shells [1S2 1P6 1D10 2S2 1F6…] and 14 p-electrons including the adding one (to form the negative charge) form the higher shell [2P6 1F8]. It is worth to note that while for the single Al atom, the 3s–3p energy gap is 3.60 eV,50 in the Al13− cluster, according to our calculations, the gap between the 1F lower shell, which is formed from 3s(Al) AOs and the 2P shell which is formed from 3p(Al) AOs, is now reduced to 0.36 eV.
Similarly, the contributions of s, p and d AOs of Al and Si atoms into the shell MOs of Al12Si are found as follows: 1S = 90% s(Si) + 9% s(Al); 1P = 22% p(Si) + 12% p(Al) + 66% s(Al); 2S = 96% s(Si) + 4% s(Al); 1D = 8% p(Al) + 91%s(Al); 2P = 40% p(Si) + 8% p(Al) + 52% s(Al); 1F = 6% p (Al) + 94% s(Al) (three lower orbitals); and 1F = 2% d(Al) + 98% p(Al) (four higher orbitals). Introduction of a Si atom into the center of Al12 also splits the 1F level into two sub-levels in which one sub-level has four-fold degeneracy (higher 1F) and the other is three-fold (lower 1F), as it is shown in Fig. 5a and b. The energy gap between them is ∼0.47 eV for Al13− and ∼0.36 eV for Al12Si. Unlike the 1F shell orbitals, the shells 1P, 1D and 2P do not split upon introduction of a Si atom into the center of the Al12 icosahedron. The shapes of higher 1F shell-orbitals (four-fold degeneracy) are similar to those of the 4f(y(3x2 − y2)), 4f(yz2), 4f(xz2) and 4f(x(3y2 − x2)), and the shapes of lower 1F shell orbitals are similar to those of 4f(xyz), 4f(z(x2 − y2)) and 4f(z3) AOs. The splitting of 1F shell orbitals described above can be understood on the basis that the central ions (Al3+ and Si4+) cause the orbitals with small angular momentum (they have larger densities in the inner part) more energetically favorable.
Herein the three lower 1F shell orbitals are constructed mainly from s(Al) AOs having small angular momentum, and the four higher 1F shell orbitals are composed of p(Al) AOs having larger angular momentum. It is also worth noting that in the formation of 2P shell orbitals, the 3p AOs of the Si atom contribute up to 40% and the contribution of 3s AOs of Al atoms drops to 52%, while in the formation of the three lower 1F shell orbitals, the 3s(Al) AOs are the dominant contributor (94%). This constitutes a reason that in the Al12Si cluster, the 2P orbitals are located below the three lower 1F orbitals as it is shown in Fig. 5b. Similarly, the MO diagram of the Al11Si2+, also having 40 valence electrons, is shown in Fig. 5c and described in detail in ESI.†
As for a comparison, orbital energies of the Al13−, Al12Si, and Al11Si2+ clusters are also illustrated in Fig. 5. Let us briefly state again that, for the Al13− anion and the Al12Si neutral, owing to an icosahedral structure with Ih symmetry, orbital degeneracy is observed for each of the 1P, 1D shells, and the 1F shell which is split into two subshells. As compared to the case of neutral Al12Si cluster, orbital energy levels of the anion Al13− are significantly higher due to the repulsion induced by the extra electron.
For Al12Si, the level ordering of 2S and 1D is reversed (relative to those in Al13−), and the 2P levels obviously get lower energy. For Al11Si2+, the 2P, 1D and 1F states markedly split due to the Si substitution and the lower symmetry, and a mixture of both 2P and 1F orbitals, as well as a mixture of the 2S level with the 1D levels are observed. However, the widths of the 1D/2S and 1F/2P sub-bands are still narrower than the spacings between other bands. In fact, this is the main feature for keeping the jellium model valid. The ellipsoidal shell structures in metal clusters were previously discussed by Brack and Clemenger61,62 and crossovers of the subshells in the deformed structures are common.
Calculated results also show that the clusters considered probably become semiconductors as suggested by the HOMO–LUMO energy gap ∼2.6 eV, though aluminum in its bulk form is a good conductor. The HOMOs of these clusters are the 1F and 2P shells, while the LUMOs are the 2D shells allow us to expect that electronic transitions from HOMO to LUMO is allowed Δ = ±1.
Moreover, our TD-DFT calculations indicate that the absorption spectra of the clusters considered, as illustrated in Fig. 6 and listed in Table S1 in ESI,† have two main characteristic absorption bands in the visible region. The absorption spectrum of the Al13− cluster is characterized by two strong bands which are centered at ∼345 (peak A) and ∼435 nm (peak B). The first peak emerges at ∼345 nm due to electronic transitions from 2P to 3S shells. The second band corresponds to electronic transitions collectively from the 1F shell orbital (HOMOs) as well as 2P shell orbitals to the 2D shell orbitals (LUMOs).
It is worth noting that the absorption spectra of small (n ≤ 5) Aln clusters are dominated by discrete orbital transition, while for the larger size clusters, more collective electron transitions are responsible for their absorption spectra.63
The absorption spectra of both clusters Al12Si and Al11Si2+ are very similar to each other and characterized by two absorption bands centered at 316–403 and 403–394 nm, respectively. Analogous to the case of Al13− anion, the longer wavelength absorption band at 403–394 nm is also due to electronic transitions collectively occurred from the 2P and 1F-the higher shell orbitals to the 2D shell orbitals. The 2P → 3S electronic transition is responsible for the shorter wavelength peak. As compared to the Al13− anion, the corresponding absorption peaks of these two Si-doped clusters shift to the shorter wavelength (Fig. 6). This is reasonable because the metal Al atom(s) are substituted by the semi-metal Si atom(s) and the net charge of the cluster is changed from negative to neutral and positive values.
For the purpose of comparison, the absorption spectra of the second isomers of the Al12Si, isomer N.12.1.B, (D5h) and Al11Si2+ C.11.2.B, (C1) clusters are also computed and plotted in Fig. 6. Similar to other isomers, the absorption spectra of N.12.1.B and C.11.2.B are also characterized by two absorption bands (A and B) in the visible region. For these isomers, the electronic transition from 1F → 2D (HOMO → LUMO) gives rise to a low intensity absorption band centered at around 450 nm. While the absorption band A remains almost unchanged (cf. Fig. 6), the absorption band B shifts from higher to lower intensity and from shorter to longer wave length, as compared to the isomers N.12.1.A and C.11.2.A. In view of the fact that N.12.1.B is 0.66 eV less stable than N.12.1.A and C.11.2.B is 0.43 eV higher in energy than C.11.2.A, their absorption bands (Fig. 6) are probably not observed, in case that the isomers are thermally populated. However, for other sizes where both lowest-lying isomers are energetically degenerate, each resulting UV-Vis absorption spectrum is likely produced by a superposition of both spectra of both corresponding A and B isomers.
The magnetic maps of Al13−, Al12Si and Al11Si2+, obtained by using the MICD approach, are shown in Fig. 7. The clockwise flows of the magnetic current point out that Al13−, Al12Si, and Al11Si2+ are aromatic structures thanks to the presence of 40 valence electrons. Furthermore, the MICD maps for each shell of Al13−, Al12Si and Al11Si2+, that are shown in Fig. 8 and 9, indicate that they are all aromatic. The blue arrows correspond to the strongest current density in each map. Fig. 8 displays the MICD maps for the 1S, 1P, 1D, and 1F shells of Al13−, Al12Si, and Al11Si2+. The positions of blue arrows in these structures move from the center to the outer of each structure when the electron shell change from the 1S to 1P to 1D and to 1F. The MICD maps confirm the validity of the shell model for these structures. As for a comparison, the MICD maps for both 2S and 2P shells are shown in Fig. 9. Obviously, the profiles of MICD maps for all the shells considered in Fig. 8 and 9 denote the electrons of 2S and 2P shells effect at both inside and outside the cage, whereas the electrons of 1S, 1P, 1D and 1F shells only induce effect within the cage of clusters.
All the doped clusters considered are characterized by low spin ground states, singlet for closed-shell species, and doublet for open-shell counterparts. For small size clusters, the Si atom prefers to be located at vertices having many edges. For larger sizes, the Si atom prefers to be endohedrally doped inside the corresponding Aln cages.
A typical characteristic in most Si doped Al clusters is a consistent degeneracy of energies of the two lowest-lying isomers, either in the neutral or the cationic state. Such a characteristic suggests the coexistence of two distinct isomeric equilibrium structures, which is expected to give rise to complicated spectra, such as vibrational, electronic or photoelectron spectra.
Relative stabilities of the clusters considered, as compared to their neighbours, were also determined. Adiabatic ionization energy and dissociation energy of each cluster size were evaluated. Calculated results confirm the high stability of the Al4Si2, Al12Si and Al11Si2+ species as “magic clusters”, that are clusters having 20 or 40 valence electrons distributed over their outer electron shells.
The UV-Vis absorption spectra of the isoelectronic magic clusters Al13−, Al12Si, and Al11Si2+ have also been modelled, showing that they have two pronounced bands centered in the regions of ∼440 and ∼350 nm, corresponding to blue light and violet light, respectively. A rationalization for the absorption spectra was assisted by an analysis of the jellium shell model of the clusters involved. Moreover, MICD maps reveal that the isoelectronic series Al13−, Al12Si and Al11Si2+, having an aromatic shell, are aromatic structures possessing 40 valence electrons.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c9ra04004h |
This journal is © The Royal Society of Chemistry 2019 |