Research Article  Open Access
Theoretical Prediction for Synthetic Realization: Pyramidal Systems ClE[E′_{4}R_{4}] (E = B–Ga, E′ = C–Ge, R = SiMe_{3}, SiMe^{t}Bu_{2}): A DFT Study
Abstract
Structure and bonding of hybrid group 13/14 pyramidal molecules ClE[R_{4}] (E = B–Ga, E′ = C–Ge, R = SiMe_{3}, SiMe^{t}Bu_{2}) were studied by DFT calculations. Six pyramidal structures are suggested for their potential synthetic realization.
1. Introduction
Among the great variety of polyhedral cluster compounds, pyramidanes are of particular importance due to their unusual structures and bonding situations, especially nonclassical interaction mode between their apexes and bases [1–4]. However, despite such evident interest, pyramidanes eluded their synthetic realization until very recently and became synthetically accessible only several years ago [5–13]. The closest approximation to the most challenging but still hypothetical allcarbon pyramidanes C[C_{4}R_{4}] represents derivatives in which apical carbon is replaced with the main group element. To date, neutral pyramidanes of the type E[C_{4}(SiMe_{3})_{4}] (E = Ge, Sn, Pb) [5, 14] were synthesized and fully structurally characterized. Moreover, pyramidanes in which the basal carbon atoms are replaced with the heavier group 14 elements (Si, Ge, and Sn) were also approached both experimentally and theoretically [14, 15]. Since the first investigations on the pyramidanes, it quickly became apparent that the structural and chemical peculiarities of the latter are totally dictated by the particular nature of their apical and basal atoms, whereas stability of pyramidanes is governed by the substituents at the basal atoms. Thus, for example, for the model E[H_{4}] (E = C–Pb, E′ = C, Si, Ge) pyramidanes, squarepyramidal C_{4v} energy minima structures were found only for pyramidanes with the carbon base (E′ = C; E = C–Pb), whereas with the heavier atoms at the base (E′ = Si, Ge; E = C–Pb) the optimized molecule represents a distorted structure [14]. However, substituents at the basal atoms play a decisive role in the overall stabilization of the squarepyramidal structures, which actually takes place upon the successive increase in the size of substituents on going from H to H_{3}Si to Me_{3}Si to ^{t}Bu_{2}MeSi. Thus, for example, most pyramidal structures E[Si_{4}(SiMe_{2})_{4}] (E = Si–Sn) correspond to the energy minima on the potential energy surface (PES), whereas for the structures E[Ge_{4}(SiMe_{2})_{4}] stationary points can be found as either energy minima or transition states (TS) with low inversion barriers [14]. Further extending the range of pyramidal structures family, we also synthesized cationic pyramidanes of the type E[C_{4}(SiMe_{3})_{4}]}^{+}∙[A]^{–} (E = P, Sb; [A]^{–} = [B(C_{6}F_{5})_{4}]^{–}, [Sb_{2}F_{7}]^{–})[16, 17] Moreover, very recently the first (chloro)borapyramidane ClB[C_{4}(SiMe_{3})_{4}] was also reported [18], being the first and still remaining the single representative of the neutral pyramidanes with the group 13 elements at the top of the square pyramid. Given very promising structural and chemical properties of the borapyramidane, it would be very challenging to expand the scope of the reaction and prepare a whole series of the hybrid group 13/14 elements pyramidanes (group 13 elements at the apex and group 14 elements at the base) and study the general trends within this series. However, it seems reasonable before pursuing experimental attempts to study the accessibility of such hybrid pyramidanes from a computational viewpoint. But such highly desirable theoretical treatment is still lacking; therefore we decided to undertake a systematic study of the general trends in stability, structural, and bonding features of hybrid pyramidanes depending on the nature of the apical and basal atoms and substituents. In this work, we report on the systematic DFT study of structure, stability, and bonding features of the series of hybrid neutral pyramidanes ClE[R_{4}] (E = B–Ga; E′ = C–Ge; R = H, SiMe_{3}, SiMe) possessing C_{4}symmetry.
2. Computational Details
All computations were performed with the Gaussian 16 [19] suite of programs at the B3LYP/Def2TZVP level [20–24]. Our previous works [14–18] showed that this level of theory is wellsuited to reproduce the Xray data of synthesized pyramidanes. The stationary points on the PES were located by full geometry optimization with calculations of force constants. XYZcoordinates for reported structures are given in the Supporting Information (SI). Natural resonance theory (NRT) [25–27] analysis was carried out by the NBO 6.0 program package [28]. NRT calculations were performed for the corresponding model systems with R = H. NICS(1)_{zz} indices [29] were calculated below the bases, that is, on the opposite side to the apex. The views of the optimized geometries were generated using Chemcraft 1.8 program [30].
3. Results and Discussions
At the beginning, the simplest hydrogensubstituted model structures were studied. Only one model system ClB[C_{4}H_{4}] 1 has symmetrical pyramidal structure as local minimum on the PES (Table 1). Other structures represent second or thirdorder saddle points on the PES (Table 1) and their minima structures have distorted geometries (for the minima geometries of systems 4–9 as well as pyramidal nonminima structures 2–9; see the SI). It is interesting that carbonbased systems 2 and 3 have housenetype structures as local minima (Figure 1). Small relative energy [0.8 kcal mol^{−1}, here and elsewhere – with zero point energy] of pyramidal (С_{4v}) 2 towards the C_{s}–minimum suggests flattened PES of ClAl[C_{4}H_{4}], and there is a lowlying (0.9 kcal mol^{−1}) TS between two degenerated C_{s}–minima (Figure 2). Thus in case of ClAl[C_{4}H_{4}] 2 there is a fast migration of ClAl fragment around C_{4}H_{4} base providing pseudopyramidal structure. For housene system ClGa[C_{4}H_{4}] 3 we also found corresponding TS with relative energy 9.8 kcal mol^{−1}.

Our previous work [14] showed the critical role of substituents for the stabilization of pyramidal structures, especially for the heavy ones with silicon and germanium bases. Thus we performed the DFT calculations of compounds 1′–9′ with real substituents that are typically used in our experimental works (Table 2, R = SiMe_{3}, SiMe^{t}Bu_{2}). As reported earlier [18], structure 1′ ClB[C_{4}(SiMe_{3})_{4}] has slightly distorted geometry with planar C_{4}base. Calculated geometrical parameters (Figure 3) are in good agreement with the experimental Xray data (C–C = 1.502, 1.493; B–C = 1.653, 1.643, 1.634; B–Cl = 1.766 Å) [18]. Geometrical trends of 1′ are in good agreement with the NRT calculations. Thus, NRT bond order for C–C bonds is 1.03 manifesting some little degree of the double bond character that corresponds to the calculated bond length value of 1.492 Å being intermediate between the standard single (1.54 Å) and double (1.34 Å) bond. Calculated NICS(1)_{zz} index [29] of –7.9 for C_{4} base indicates electron delocalization in basal cycle what keeps planarity of base and provides the pyramidal geometry for the system. For the comparison, NICS(1)_{zz} for the benzene is –29.9 at the same level of theory. Pyramidal boroncarbon bonds in 1′ of 1.637 and 1.642 Å are longer than the sum of the single bond covalent radii of B and C atoms (1.60 Å) [31]. NRT bond order for B–C bonds is only 0.63. Nevertheless, B–C bonds are mostly covalent: the NRT covalent character is 58.8 %. The short length of B–Cl bond of 1.775 Å (the sum of the single bond covalent radii of boron and chlorine atoms is 1.84 Å) [31] well corresponds to the NRT bond order of 1.06 indicating some participation of double bonding.

Calculations of 2′ (E = Al) also predict stable pyramidal structure with symmetrical C_{4}geometry (Figure 3). C–C bonds in 2′ of 1.506 Å are longer that in 1′, but still intermediate between the standard single and double bonds, and have the greater NRT bond orders of 1.16. Aromaticity of the C_{4}base is proved by the NICS(1)_{zz} value of –10.9. Pyramidal Al–C bonds of 2.013 Å are only slightly longer than the sum of the single bond covalent radii of Al and C atoms (2.01 Å) [31]. NRT bond order for these bonds is 0.49, and the ionic character is 61.8 %. Apical ClAlgroup has the Al–Cl bond value of 2.086 Å (the sum of the single bond covalent radii of Al and Cl atoms is 2.25 Å) [31], NRT Al–Cl bond order is 1.05, and the bond is by mostly ionic (NRT ionic character is 62.4 %).
In contrast with 1′ and 2′, system 3′ (E = Ga) is represented by the housenetype structure, as shown in Figure 3. Symmetric pyramidal structure represents secondorder saddle point (Table 2, for geometry see Figure S4 in the SI). As we mentioned above, there is a quite lowlying TS for the model housene structure 3 ClGa[C_{4}H_{4}]. Bulkier substituents might further decrease the energy of TS for such migrations [16]. The calculated barrier to migration is only 5.3 kcal mol^{−1} for H_{3}Sisubstituted model system 3 ClGa[C_{4}(SiH_{3})_{4}], as shown in Figure 2. Such value is reasonable to consider a room temperature process for the migration of ClGagroup in 3′ that finally provides pseudopyramidal structure for it. Decrease in the migration barrier is accompanied by the decreasing of deviation angle of gallium towards C_{4}base (see Figure 1): 102.7° for R = H, 99.3° for R = SiH_{3}, and 97.9° for R = SiMe_{3}.
Siliconbased structures 4′ (E = B, R = SiMe^{t}Bu_{2}), 5′ (E = Al, R = SiMe^{t}Bu_{2}), and 6′ (E = Ga, R = SiMe^{t}Bu_{2}) all have perfect pyramidal C_{4}symmetry geometry at their energy minima points on the PES (Figure 4, Table 2) and are therefore good candidates for the synthetic realization. Si–Si bonds in 4′ of 2.279 Å are intermediate between the standard single (2.32 Å) [31] and double bond (2.14 Å) [32], and NICS(1)_{zz} index of –9.6 approves electron delocalization in basal cycle of 4′. Nevertheless NRT bond order of 0.97 indicate single Si–Si bonding in 4′. Structures 5′ and 6′ have Si–Si bonds of 2.32–2.33 Å corresponding to the single bond values. But NRT bond orders of 1.05 and 1.05 describe Si–Si bonds as having some degree of double bond character, and NICS(1)_{zz} values of –10.2 and –10.9 indicate electron delocalization in Si_{4}bases of 5′ and 6′, respectively. As expected, pyramidal E–Si bonds in 4′, 5′, 6′ are longer than the sum of the single bond covalent radii of E and Si atoms (2.01, 2.42, 2.40 [31] Å for E = B, Al, Ga, respectively) and have low NRT bond orders of 0.69, 0.60, 0.59. Covalent characters calculated within the NRT are 75.6, 54.5, 61.9 % for the E–Si bonds (E = B, Al, Ga) in 4′, 5′, 6′, respectively. E–Cl bonds in 4′, 5′, 6′ are shorter than the corresponding standard single bonds [31] by 0.02, 0.12, and 0.04 Å, respectively, and their NRT results point to a little extent of double bonding: the NRT bond orders are 1.02, 1.04, 1.00 for B–Cl, Al–Cl, and Ga–Cl, respectively.
Among the systems with the heaviest Ge_{4}base there are no true pyramidanes: 7′ (E = B, R = SiMe^{t}Bu_{2}), 8′ (E = Al, R = SiMe^{t}Bu_{2}), and 9′ (E = Ga, R = SiMe^{t}Bu_{2}) minima structures all have distorted bases (like structure 7′ in Figure 4, for the minima structures 8′ and 9′ see Figure S3 in the SI). But among them the system ClB[Ge_{4}(SiMe^{t}Bu_{2})_{4}] 7′ could be interpreted as effective pyramidal because of its lowlying (+2.7 kcal mol^{−1}) TS representing true squarepyramidal geometry (Figure S4 in the SI). Such phenomenon was previously observed in pentagermapyramidane Ge[Ge_{4}(SiMe^{t}Bu_{2})_{4}] [15] and can be explained by the fast conformational interconversion between the folded minima and squareplanar pyramidal TS (Scheme 1).
Four Ge–Ge bonds of 2.458 Å in minimum structure 7′ are longer than the standard single bonds (2.42 Å) [31] but have NRT bond order of 1.01 manifesting little contribution of the double bond character. The dihedral angle involving four Ge atoms is 26°. The TS structure characterized by planar Ge_{4}base, shorter Ge–Ge bonds of 2.410 Å and NICS(1)_{zz} index of –11.8. Two B–Ge bonds of 2.084 Å in minimum structure 7′ are close to the standard single bonds (2.06 Å) [31]. The other two B–Ge distances of 2.539 Å are much longer than the normal single bonds. B–Cl bond of 1.806 Å is shorter than the standard single bond (1.84 Å) [31], its NRT bond order is 1.00, and covalent character is 67 %.
Described structures could also be realized with other halogen atoms, and we calculated a structure BrB[C_{4}(SiMe_{3})_{4}] 10′, as such example. System 10′ corresponds to the energy minimum on the PES and has the same pyramidal geometry, as ClB[C_{4}(SiMe_{3})_{4}] 1′ (Figure 3). C–C bonds in 10′ of 1.494 Å have NRT bond orders of 1.03. The NICS(1)_{zz} index of –4.4 for C_{4} base is nearly twice smaller than that for C_{4} base in 1′ (–7.9). Pyramidal B–C bonds in 10′ of 1.637 and 1.640 Å are elongated comparing with single B–C bonds and have NRT bond orders of only 0.64. B–Br bond of 1.930 Å is shorter than the sum of the single bond covalent radii of boron and bromine atoms (1.99 Å [31]) and has the corresponding NRT bond order of 1.05 manifesting some extent of the double bonding.
4. Conclusions
According to our DFT calculations, synthetic realization of six pyramidal structures with the group 13 elements apical atom seems to be feasible: 1′ ClB[C_{4}(SiMe_{3})_{4}] (experimentally realized), 2′ ClAl[C_{4}(SiMe_{3})_{4}], 3′ ClGa[C_{4}(SiMe_{3})_{4}] (effective pyramidal), 4′ ClB[Si_{4}(SiMe^{t}Bu_{2})_{4}], 5′ ClAl[Si_{4}(SiMe^{t}Bu_{2})_{4}], 6′ ClGa[Si_{4}(SiMe^{t}Bu_{2})_{4}], and 7′ ClB[Ge_{4}(SiMe^{t}Bu_{2})_{4}] (effective pyramidal). Structures with other halogen atoms could also be realized.
Data Availability
The geometry optimization data used to support the findings of this study are included within the supplementary information files.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was financially supported by the Russian Ministry of Education and Science for Research (Project Part, Project 4.844.2017/4.6), and by the JSPS KAKENHI Grants program (No. JP18K05137) from the Ministry of Education, Science, Sports, and Culture of Japan.
Supplementary Materials
Additional Supporting Information may be found online in the supporting information tab for this article. (Supplementary Materials)
References
 V. I. Minkin, R. M. Minyaev, and R. Hoffmann, “Nonclassical structures of organic compounds: unusual stereochemistry and hypercoordination,” Russian Chemical Reviews, vol. 71, pp. 869–892, 2002. View at: Google Scholar
 V. I. Minkin, “Current trends in the development of A. M. Butlerov’s theory of chemical structure,” Russian Chemical Bulletin, vol. 61, pp. 1265–1290, 2012. View at: Google Scholar
 E. Lewars, Computational Chemistry (Introduction to the Theory and Applications of Molecular and Quantum Mechanics), Springer, 2nd edition, 2011.
 E. Lewars, Modeling Marvels (Computational Anticipation of Novel Molecules), Springer, 2008.
 V. Ya. Lee, Y. Ito, A. Sekiguchi et al., “Pyramidanes,” Journal of the American Chemical Society, vol. 135, no. 24, pp. 8794–8797, 2013. View at: Publisher Site  Google Scholar
 M. D. Francis and P. B. Hitchcock, “Synthesis, crystal and molecular structure of [Sn(η^{4})]: the first non transition metal 1,3diphosphacyclobutadienyl compound,” Chemical Communications, pp. 8687, 2002. View at: Google Scholar
 C. Condick, M. A. Fox, R. Greatrex, C. Jones, and D. L. Ormsby, “A new nido5vertex cluster, phosphacarbanidopentaborane, 21,2PCB_{3}H_{5},” Chemical Communications, vol. 2, no. 14, pp. 14481449, 2002. View at: Publisher Site  Google Scholar
 J. M. Lynam, M. C. Copsey, M. Green et al., “Selective preparation of the [3,5‐tBu_{2}‐1,2,4‐C_{2}P_{3}]^{−} ion and synthesis and structure of the cationic species nido‐[3,5‐tBu_{2}‐1,2,4‐C_{2}P_{3}]^{+}, isoelectronic with [C_{5}R_{5}]^{+},” Angewandte Chemie International Edition, vol. 42, no. 24, pp. 2778–2782, 2003. View at: Publisher Site  Google Scholar
 M. D. Francis and P. B. Hitchcock, “Studies of the 24Electron nido5Vertex Clusters [M(η^{4})] (M = Ge, Pb, Sn): Synthesis and Characterization of [M(η^{4})] (M = Ge, Pb), [Sn(η^{4})W(CO)_{5}_{n}] (n = 1, 2), and [I_{2}Ge] and the Mössbauer Spectrum of [Sn(η^{4})],” Organometallics, vol. 14, pp. 2891–2896, 2003. View at: Google Scholar
 G. Anderson, J. C. Green, and M. D. Francis, “Electronic structure of M(η^{4}), where M = Ge, Sn, and Pb: photoelectron spectroscopy and density functional studies,” Organometallics, vol. 14, pp. 2897–2901, 2003. View at: Google Scholar
 D. A. Pantazis, J. E. McGrady, J. M. Lynam, C. A. Russell, and M. Green, “Structure and bonding in the isoelectronic series : is phosphorus a carbon copy?” Dalton Transactions, no. 14, pp. 2080–2086, 2004. View at: Google Scholar
 C. Fish, M. Green, J. Jeffery et al., “Cationic phosphorus–carbon–pnictogen cages isolobal to [C_{5}R_{5}]^{+},” Chemical Communications, no. 13, pp. 1375–1377, 2006. View at: Google Scholar
 C. Fish, M. Green, J. C. Jeffery et al., “A maingroup analogue of housene: the subtle influence of the inertpair effect in group 15 clusters,” Angewandte Chemie International Edition, vol. 45, no. 40, pp. 6685–6689, 2006. View at: Publisher Site  Google Scholar
 V. Ya. Lee, O. A. Gapurenko, Y. Ito et al., “Pyramidanes: the covalent form of the ionic compounds,” Organometallics, vol. 35, no. 3, pp. 346–356, 2016. View at: Publisher Site  Google Scholar
 V. Ya. Lee, Y. Ito, O. A. Gapurenko et al., “Pentagermapyramidane: crystallizing the "transitionstate" structure,” Angewandte Chemie International Edition, vol. 54, no. 19, pp. 5654–5657, 2015. View at: Google Scholar
 V. Ya. Lee, H. Sugasawa, O. A. Gapurenko et al., “A cationic phosphapyramidane,” Chemistry  A European Journal, vol. 22, no. 49, pp. 17585–17589, 2016. View at: Publisher Site  Google Scholar
 V. Ya. Lee, K. Ota, Y. Ito et al., “Bis(stibahousene),” Journal of the American Chemical Society, vol. 139, no. 39, pp. 13897–13902, 2017. View at: Publisher Site  Google Scholar
 V. Ya. Lee, H. Sugasawa, O. A. Gapurenko et al., “From borapyramidane to borole dianion,” Journal of the American Chemical Society, vol. 140, no. 19, pp. 6053–6056, 2018. View at: Publisher Site  Google Scholar
 M. J. Frisch, G. W. Trucks, H. B. Schlegel et al., Gaussian 16, Revision A.03, Gaussian, Inc., Wallingford, Conn, USA, 2016.
 A. D. Becke, “Densityfunctional exchangeenergy approximation with correct asymptotic behavior,” Physical Review A: Atomic, Molecular and Optical Physics, vol. 38, no. 6, pp. 3098–3100, 1988. View at: Publisher Site  Google Scholar
 C. Lee, W. Yang, and R. G. Parr, “Development of the ColleSalvetti correlationenergy formula into a functional of the electron density,” Physical Review B: Condensed Matter and Materials Physics, vol. 37, no. 2, pp. 785–789, 1988. View at: Publisher Site  Google Scholar
 B. Miehlich, A. Savin, H. Stoll, and H. Preuss, “Results obtained with the correlation energy density functionals of becke and Lee, Yang and Parr,” Chemical Physics Letters, vol. 157, no. 3, pp. 200–206, 1989. View at: Publisher Site  Google Scholar
 F. Weigend and R. Ahlrichs, “Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: design and assessment of accuracy,” Physical Chemistry Chemical Physics, vol. 7, no. 18, pp. 3297–3305, 2005. View at: Publisher Site  Google Scholar
 F. Weigend, “Accurate Coulombfitting basis sets for H to Rn,” Physical Chemistry Chemical Physics, vol. 8, no. 9, pp. 1057–1065, 2006. View at: Publisher Site  Google Scholar
 E. D. Glendening and F. J. Weinhold, “Natural resonance theory: I. General formalism,” Journal of Computational Chemistry, vol. 19, no. 6, pp. 593–609, 1998. View at: Publisher Site  Google Scholar
 E. D. Glendening and F. J. Weinhold, “Natural resonance theory: II. Natural bond order and valency,” Journal of Computational Chemistry, vol. 19, no. 6, pp. 610–627, 1998. View at: Publisher Site  Google Scholar
 E. D. Glendening, J. K. Badenhoop, and F. J. Weinhold, “Natural resonance theory: III. Chemical applications,” Journal of Computational Chemistry, vol. 19, no. 6, pp. 628–646, 1998. View at: Google Scholar
 E. D. Glendening, J. K. Badenhoop, A. E. Reed et al., NBO 6.0, Theoretical Chemistry Institute, University of Wisconsin, Madison, Wis, USA, 2013.
 H. FallahBagherShaidaei, C. S. Wannere, C. Corminboeuf, R. Puchta, and P. V. R. Schleyer, “Which NICS aromaticity index for planar π rings is best?” Organic Letters, vol. 8, no. 5, pp. 863–866, 2006. View at: Publisher Site  Google Scholar
 G. A. Zhurko, “ChemCraft software,” http://www.chemcraftprog.com. View at: Google Scholar
 P. Pyykkö and M. Atsumi, “Molecular single‐bond covalent radii for elements 1–118,” Chemistry – A European Journal, vol. 15, no. 1, pp. 186–197, 2009. View at: Google Scholar
 P. Pyykkö and M. Atsumi, “Molecular double‐bond covalent radii for elements Li–E112,” Chemistry – A European Journal, vol. 15, no. 46, pp. 12770–12779, 2009. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2019 Olga A. Gapurenko 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.