The key geometrical parameters, the singlet-triplet energy splitting (ΔEST), the natural charge densities (QAl and QSb), the binding energies (BE), the HOMO-LUMO energy gaps, and the Wiberg bond index (WBI) for RAl☰SbR using the M06-2X/Def2-TZVP, B3PW91/Def2-TZVP (in round brackets) and B3LYP/LANL2DZ + dp (in square brackets) levels of theory.
The effect of substitution on the potential energy surfaces of RAl☰SbR (R = F, OH, H, CH3, SiH3, SiMe(SitBu3)2, SiiPrDis2, Tbt, and Ar*) is investigated using density functional theories (M06-2X/Def2-TZVP, B3PW91/Def2-TZVP, and B3LYP/LANL2DZ + dp). The theoretical results demonstrated that all the triply bonded RAl☰SbR compounds with small substituents are unstable and can spontaneously rearrange to other doubly bonded isomers. That is, the smaller groups, such as R = F, OH, H, CH3 and SiH3, neither kinetically nor thermodynamically stabilize the triply bonded RAl☰SbR compounds. However, the triply bonded R’Al☰SbR´ molecules that feature bulkier substituents (R´ = SiMe(SitBu3)2, SiiPrDis2, Tbt, and Ar*) are found to possess the global minimum on the singlet potential energy surface and are both kinetically and thermodynamically stable. In particular, the bonding characters of the R’Al☰SbR´ species agree well with the valence-electron bonding model (model) as well as several theoretical analyses (the natural bond orbital, the natural resonance theory, and the charge decomposition analysis). That is to say, R’Al☰SbR´ molecules that feature groups are regarded as R′─Al Sb─R′. Their theoretical evidence shows that both the electronic and the steric effects of bulkier substituent groups play a decisive role in making triply bonded R′Al☰SbR′ species synthetically accessible and isolable in a stable form.
- group 13 elements
- group 13 elements
- triple bond
The chemical synthesis and structural characterization of molecules that feature triple bonds  between heavier group 14 elements (E14 = Si, Ge, Sn and Pb) are of interest because of their interesting structural chemistry and their potential applications in organic and inorganic synthesis [1, 2, 3, 4, 5, 6, 7, 8, 9, 10]. Although understanding of these RE14☰E14R molecules that feature heavier group 14 atoms has increased during the last two decades, the understanding of the RE13☰E15R compounds, which are isoelectronic to acetylene from a valence electron viewpoint, is still limited. The reason for this limited knowledge of acetylene analogues, RE13☰E15R, could be due to the fact that there has been limited preparation and the isolation of these species in a stable form [11, 12]. Theoretical methods allow a theoretical design of the RE13☰E15R molecules to be made that increases understanding of their potential properties.
The III-V semiconductors that contain antimony have several important applications in optoelectronic devices that operate in the infrared region and in high-speed devices, which has prompted widespread studies of promising precursor systems for these materials . In particular, the chemical synthesis and structural characterization of AlSb single-source precursors of the type R3Al-SbR´3 has attracted much attention, owing to their importance in CVD procedures , which is a developing industry for the production of thin films of the corresponding semiconducting materials . As far as the authors are aware, only a handful of group 13 antimonides that contain Al─Sb σ-bonds have been discovered , No triply bonded RAl☰SbR species, which is isoelectronic to HC☰CH, has been reported both experimentally and theoretically.
Density functional theory (DFT) is sued to determine the structures, the kinetic stability and bonding properties of various RAl☰SbR triply bonded forms on the singlet ground state, in order to obtain a better understanding of aluminum☰antimony triple bonds. This work reports the possible existence of triply bonded RAl☰SbR molecules, from the viewpoint of the effect of substituents, using DFT . That is, M06-2X/Def2-TZVP, B3PW91/Def2-TZVP and B3LYP/LANL2DZ + dp are used for small substituents (R = H, F, OH, CH3, and SiH3) and M06-2X/Def2-TZVP  for large substituents (R = SiMe(SitBu3)2, SiiPrDis2, Tbt, and Ar*; see Scheme 1) .
2. General considerations
The valence-bond bonding model is a well-known satisfactory method, which is an approximate theory to explain the electron pair or chemical bond by quantum mechanics, for predicting molecular geometries . Two valence-bond bonding models (Figure 1) are thus used to interpret the bonding properties of triply bonded RAl☰SbR species. In model , the RAl☰SbR molecule is partitioned into two units: a singlet R─Al and a singlet R─Sb. In model , the RAl☰SbR compound is divided into two moieties: a triplet R─Al and a triplet R─Sb. As a result, the choice of the bonding model that is used to explain the bonding characters of RAl☰SbR depends on the promotion energies (ΔEST = Etriplet − Esinglet) of the R─Al and R─Sb fragments. According to current theoretical calculations (see below), it is known that R─Al occupies the singlet ground state, but R─Sb occupies the triplet ground state. In consequence, if the value of ΔEST for R─Al is much larger than that for R─Sb, the latter easily jumps to the singlet excited state. Hence, model  can be used to explain the bonding nature of the RAl☰SbR molecule. In contrast, if the value of ΔEST for R─Al is smaller than that for R─Sb, the former is readily promoted to the excited triplet state. Therefore, model  is used to interpret the bond constitutions of the RAl☰SbR compound.
Two points are worthy of note. The first is that since aluminum and antimony respectively belong to group 13 and group 15 and both elements have different atomic radii (covalent radii: 118 pm and 140 for Al and Sb, respectively) , the overlapping populations between Al and Sb should not be strong. The second is that the lone pairs of both aluminum and antimony feature the valence s character. This, in turn, makes the overlap integrals between the lone pair orbital and the pure p orbital small. These two factors mean that the triple bond between aluminum and antimony is weak, unlike the traditional triple bond in acetylene.
Bearing the above bonding analyses in mind, theoretical evidences are given in the following sections.
3. Results and discussion
3.1. Small ligands on substituted RAl☰SbR
Five small substituents (R = F, OH, H, CH3 and SiH3) are chosen, which include electronegative and electropositive groups, to determine their stability and bonding properties on the triply bonded RAl☰SbR molecules using the three types of DFT calculations (i.e., M06-2X/Def2-TZVP, B3PW91/Def2-TZVP and B3LYP/LANL2DZ + dp). Figure 2 shows the potential energy surfaces of the intra-molecular 1,2-migration reactions for five triply bonded RAl☰SbR compounds that feature small substituents. That is to say, the triply bonded RAl☰SbR species can undergo a 1,2-shift to give either R2Al〓Sb: or: Al〓SbR2 doubly bonded isomers.
As seen in Figure 2, the three DFT computational results demonstrate that the triply bonded RAl☰SbR species that feature small substituents are all both kinetically and thermodynamically unstable on the intra-molecular 1,2-migration reaction potential energy surfaces. In other words, once the triply bonded RAl☰SbR with small substituents is formed, it can easily proceed along the 1,2-migration to give the thermodynamically stable doubly bonded isomer, either R2Al〓Sb: or: Al〓SbR2. The theoretical findings give strong evidence that the triply bonded RAl☰SbR molecules that feature the small ligands are highly unlikely to be detected experimentally.
Although current theoretical observations show that the formation of RAl☰SbR involving small ligands is not likely, some of their physical properties, which are shown in Table 1, must be theoretically determined in order to design much more stable aluminum☰antimony acetylene analogues.
|ΔEST for Al-R (kcal/mol)3||79.78|
|ΔEST for Sb-R (kcal/mol)4||−32.40|
As seen in Table 1, the three DFT computational results predict that the Al☰Sb triple bond distance (Å) is in the ranges 2.388–2.539 (M06-2X/Def2-TZVP), 2.397–2.536 (B3PW91/Def2-TZVP) and 2.436–2.565 (B3LYP/LANL2DZ + dp). Table 1 also shows that all of the geometrical structures of RAl☰SbR adopt the bent form, as demonstrated in Scheme 2. That is, ∠R─Al─Sb ≈ 180.0° and ∠Al─Sb─R ≈ 90.0°. The reason for this vertical angle at the Sb center can be ascribed to the relativistic effect, as discussed previously . The three DFT calculations shown in Table 1 all indicate that the electronic ground states for R─Al and the R─Sb fragments are singlet and triplet, respectively. In particular, all of the DFT results shown in Table 1 show that most of the singlet-triplet energy splitting (ΔEST) of R-Al is larger than that of the corresponding R─Sb. This strongly implies that the bonding characters of the triply bonded RAl☰SbR species that feature small substituents are better described by model , as shown in Figure 1. In other words, the triple bond consists of one donor-acceptor σ bond and two donor-acceptor π bonds, which are schematically represented as R─Al
In brief, the three DFT calculations shown in this work show that irrespective of their electronegativity, the triply bonded RAl☰SbR molecules that feature small ligands are highly unlikely to exist, even in the low-temperature matrices. In particular, the bond orders of these Al☰Sb triple bonds are theoretically predicted to be a weak double bond, rather than a triple bond.
3.2. Large ligands on substituted R’Al☰SbR´
Three bulky groups were then used to search for kinetically stable triple-bonded R’Al☰SbR´ molecules: R´(〓SiMe(SitBu3)2, SiiPrDis2, Tbt, and Ar*) . These are shown in Scheme 1. It is known that London dispersion (nonvalent interactions) plays a prominent role in both chemical and physical properties of inorganic molecules . As a result, the dispersion-corrected M06-2X/Def2-TZVP method is used in the present study to investigate the behaviors of the triply bonded R’Al☰SbR´ compounds bearing bulky substituents. Similarly to the cases for small ligands on substituted RAl☰SbR, the dispersion-corrected M06-2X/Def2-TZVP level of theory is used to determine the potential energy surfaces for the intra-molecular 1,2-migration reactions of R’Al☰SbR´, as shown in Scheme 3. The computed relative energies are listed in Table 2. The reaction enthalpies for both the 1,2-shift reactions (R’Al☰SbR´ → R2’Al〓Sb and R’Al☰SbR´ → R2’Sb〓Al) are apparently too high. They are estimated to be at least 80 kcal/mol. The reason that both doubly bonded R2’Al〓Sb and R2’Sb〓Al isomers occupy such high energy points is simply because two bulky groups can cause steric overcrowding. As a consequence, the theoretical findings strongly suggest that the triply bonded R’Al☰SbR´, which is attached by two bulkier substituents, is kinetically stabilized.
|ΔEST for Al─R′ (kcal/mol)3||28.89||27.30||42.50||40.21|
|ΔEST for Sb─R′ (kcal/mol)4||−16.89||−24.80||−30.51||−15.92|
Table 2 shows that the Al☰Sb triple bond distance is predicted to be 2.422–2.477 Å. Since no experimental results for the Al☰Sb triple bond length have been reported, these values are estimates. These theoretical calculations also show that the geometrical structures of R′Al☰SbR′ molecules that feature bulky groups adopt a bent structure; i.e., ∠R′─Al─Sb ≈ 160.0° and ∠Al─Sb─R′ ≈ 120.0°. As stated previously, the triply bonded R′Al☰SbR′ species feature this bent geometry because of the relativistic effect .
In addition, the bonding energy (BE) that is shown in Table 2 shows that the central aluminum and antimony atoms in the substituted R′Al☰SbR′ compounds are strongly bonded, since the BE values are in the range 71–97 kcal/mol for R′ = SiMe(SitBu3)2, SiiPrDis2, Tbt, and Ar*. Table 2 also shows that the modulus ΔEST (kcal/mol) for Al─R′ and Sb─R′ fragments are predicted to be 43–27 and 31–16. These theoretical values allow two interpretations. Firstly, even when attached by bulkier groups, it is theoretically verified that both the Al─R′ and the Sb─R′ units occupy the ground singlet state and the ground triplet state, respectively. Since the ΔEST values for Al─R′ are so small (compared with those for Al─R, as shown in Table 1), model  in Figure 1 is most suitable to interpret the triple bonding characters in the R′Al☰SbR′ species that feature bulky substituents. As schematically shown in Figure 1, the nature of the Al☰Sb triple bond can be considered as one conventional σ bond, one conventional π bond and one donor-acceptor π bond. That is, R′─Al
Besides these, Dapprich and Frenking developed a useful method , which is called the introduced charge decomposition analysis (CDA), from which one may analyze donor-acceptor interactions of a A-B molecule. From CDA, one may obtain three parts. The first part is the number of electrons donated from the R′─Al unit to the R′─Sb monomer, which can be considered as (R′─Al) → (R′─Sb). The second part is the number of electrons back donated from the R′─Sb component to the R′─Al moiety, which can be represented as (R′─Al) ← (R′─Sb). The third part is the repulsive interactions between (R′─Al) and (R′─Sb), which can be described as (R′─Al) ↔ (R′─Sb). The CDA results about the (SiMe(SitBu3)2)Al☰Sb(SiMe(SitBu3)2) molecule based on the dispersion-corrected M06-2X/Def2-TZVP method are given in Table 3. As seen in Table 3, for the (R′─Sb) fragment, its largest contribution is No. 267 (HOMO) orbital, displaying that a R′─Sb component donates electrons to a R′─Ga unit mainly through the HOMO orbital. In consequence, the net amount of electron transfer is estimated to be −0.207, implying that the R′─Sb part donates more electrons to the R′─Al moiety. This theoretical finding agrees well with the valence-electron bonding model shown in Figure 1 (i.e., model ). Namely, the bonding character of R′Al☰SbR′ can be recognized as R′Al
The bonding characters of the Al☰Sb triple bond in R′Al☰SbR′ molecules were examined using the natural bond orbital (NBO)  and the natural resonance theory (NRT)  analysis, whose results are given in Table 4, are used to determine the bonding properties. For instance, Table 4 shows that for (SiMe(SitBu3)2)Al☰Sb(SiMe(SitBu3)2), the NBO model shows that the Al-Sb σ bonding orbital contains about 23% natural Al orbitals and 77% natural Sb orbitals. Also, the Al☰Sb π bonding orbital contains averagely about 25% natural Al orbitals and 75% natural Sb orbitals (Figure 3). These values give strong evidence that the Al☰Sb π bond is polarized. Table 4 also shows that the Al☰Sb π bonding interaction: π⊥(Al☰Sb) = 0.529(3s3p1.98)Al + 0.849(5s5p12.43)Sb and π‖(Al☰Sb) = 0.475(3s3p99.99)Al + 0.880(5s5p99.99)Sb, which again implies that the predominant bonding interaction between the Al─R and the Sb─R moieties originates from 3p(Al) ← 5p(Sb) donation. In other words, the electron deficiency on Al and the π bond polarity are partially balanced by the donation of the Sb lone pair to the empty Al p orbital (Figure 3). Table 4 also shows that, on the basis of the NRT analyses of the electron density for (SiMe(SitBu3)2)Al☰Sb(SiMe(SitBu3)2), its Al☰Sb triple bond has a greater covalent character, as shown by the greater covalent part of the NRT bond order (1.25), compared to its ionic part (0.81). The reason for this may be due to the fact that the difference between the electronegativity values for the Al and Sb elements is small (Al: 1.5 and Sb: 1.8) .
|R’Al☰SbR’||WBI||NBO analysis||NRT analysis|
|R’ = SiMe(SitBu3)2||2.17||σ: 1.91||σ: 0.4799 Al (sp3.23) + 0.8773 Sb (sp0.60)||23.03% (Al)|
|π⊥: 1.81||π⊥: 0.5288 Al (sp1.98) + 0.8487 Sb (sp12.43)||27.96% (Al)|
|π‖: 1.89||π‖: 0.4753 Al (sp99.99) + 0.8798 Sb (sp99.99)||22.59% (Al)|
|R’ = SiiPrDis2||2.18||σ: 1.91||σ: 0.5525 Al (sp1.71) + 0.8335 Sb (sp1.15)||30.53% (Al)|
|π⊥: 1.86||π⊥: 0.4723 Al (sp3.67) + 0.8815 Sb (sp3.68)||22.30% (Al)|
|π‖: 1.89||π‖: 0.4476 Al (sp99.99) + 0.8943 Sb (sp99.99)||20.03% (Al)|
|R’ = Tbt||2.07||σ: 1.95||σ: 0.6923 Al (sp0.18) + 0.7216 Sb (sp12.38)||47.93% (Al)|
|π⊥: 1.88||π⊥: 0.4488 Al (sp47.14) + 0.8936 Sb (sp99.99)||20.14% (Al)|
|π‖: 1.91||π‖: 0.4772 Al (sp99.99) + 0.8788 Sb (sp99.99)||22.78% (Al)|
|R’ = Ar*||2.02||σ: 1.96||σ: 0.6946 Al (sp0.16) + 0.7194 Sb (sp18.14)||48.25% (Al)|
|π⊥: 1.83||π⊥: 0.4543 Al (sp99.99) + 0.8908 Sb (sp40.30)||20.64% (Al)|
|π‖: 1.92||π‖: 0.4266 Al (sp99.99) + 0.9044 Sb (sp99.99)||18.20% (Al)|
This study uses DFT computations to theoretically design substituted RAl☰SbR molecules that feature the Al☰Sb triple bond, that are stable from the kinetic viewpoint. The theoretical observations show that only bulky substituents (R′) can significantly stabilize the triply bonded R′Al☰SbR′ compounds, and not small substituents. The theoretical findings also show that the bonding characters of the R′Al☰SbR′ species that feature bulky groups can be represented as R′─Al
The authors are grateful to the National Center for High-Performance Computing of Taiwan for generous amounts of computing time, and the Ministry of Science and Technology of Taiwan for the financial support.