Si Jia Isabel
Phang‡
a,
Zheng-Feng
Zhang‡
b,
Chi-Shiun
Wu
b,
Zhen Xuan
Wong
a,
Ming-Der
Su
*bc and
Cheuk-Wai
So
*a
aSchool of Chemistry, Chemical Engineering and Biotechnology, Nanyang Technological University, Singapore 637371, Singapore. E-mail: CWSo@ntu.edu.sg
bDepartment of Applied Chemistry, National Chiayi University, Chiayi 60004, Taiwan
cDepartment of Medicinal and Applied Chemistry, Kaohsiung Medical University, Kaohsiung 80708, Taiwan
First published on 27th January 2025
The replacement of all carbon atoms in aromatic rings with main-group elements to afford inorganic ring systems is highly desirable due to their distinct aromatic character. However, fused polycyclic main-group element rings are rare and the feasibility of aromaticity in such compounds has yet to be explored. To explore aromaticity in fused polycyclic main-group element rings, a stable di-silicon analogue of fused bicyclic borirene, namely bicyclo[1.1.0]-2,4-diborylenyldisil-1(3)-ene 5 was synthesized from an N-phosphinoamidinato chlorosilylene 3. Compound 5 consists of a bridgehead SiSi double bond bonded with two bridging borons resulting in an unsaturated fused bicyclic skeleton. The bridgehead Si
Si σ- and π-electrons and bridging Si–B σ-electrons are stabilized by both σ- and π-aromatic delocalization on the Si2B2 fused bicyclic ring.
It is noteworthy that compounds I–IX are all mono-cyclic aromatic rings, and are representative main group element analogues of aromatic carbon molecules, such as cyclobutane-1,3-diyl and the cyclopropenyl cation. Based on fused polycyclic aromatic hydrocarbons such as naphthalene, it is plausible to extend the concept of aromaticity to fused polycyclic inorganic rings, which have yet to be discovered.
Among fused polycyclic aromatic rings, bicyclic borirene X is the smallest neutral fused aromatic ring.17 It features bicyclic 2π-aromatic delocalization, resulting from the overlapping of the bridgehead CC π orbital with two empty p orbitals on the bridging boron center. Compared to typical alkenes such as ethylene and propylene, high strain energy owing to the exceptionally strained bridgehead C
C double bond in the bicyclic scaffold, arising from an inversion of the conformation about the C
C double bond, is expected. The isolation of such a ring remains unknown experimentally. Thus, in an effort to isolate a more stable analogue of bicyclic borirene X, we replaced the carbon atoms in the ring with silicon to offer more chemical adaptability and reduced ring strain due to the larger atomic size of silicon. We further used a bidentate ligand bonded with the boron atoms to provide enhanced kinetic stabilization. Herein, we report a phosphinoamidinate-stabilized bicyclo[1.1.0]-2,4-diborylenyldisil-1(3)-ene, where the bridgehead Si
Si σ- and π-electrons and bridging Si–B σ-electrons are aromatically delocalized on the Si2B2 fused bicyclic ring, accounting for the stability of the fused bicyclic borirene ring structure.
Compound 3 underwent oxidative addition with BI3 in toluene at 65 °C for 16 h to afford the N-phosphinoamidinate-bridged borylsilane 4 (Scheme 2), which was isolated as a colorless crystalline solid (yield: 55%). The 29Si{1H} NMR signal of 4 (−8.8 ppm) is broad due to quadrupolar coupling with the B nucleus. The 31P{1H} and 11B{1H} NMR signals are found at 26.6 ppm and −45.6 ppm, respectively. X-ray crystallography shows a six-membered ring, where the ligand is an N,P-chelate, bridged between the Si–B bond (Fig. 1b). The Si1–B1 bond (2.021(5) Å) is typical of a single bond. A similar B–Cl bond oxidative addition was observed in the reactivity of a disilicon(I) compound.20
Compound 4 was reacted with excess KC8 in THF at room temperature for 3 h to afford N-phosphinoamidinato bicyclo[1.1.0]-2,4-diborylenyldisil-1(3)-ene 5, which was isolated as a reddish-orange crystalline solid (yield: 10%). Compound 5 is considered a bridging boron analogue of the bicyclo[1.1.0]tetrasil-1(3)-ene reported by Iwamoto et al.,21–25 but the presence of the bridging N-phosphinoamidinato boron moieties in compound 5 induces some degree of difference in electronic properties compared with the bicyclo[1.1.0]tetrasil-1(3)-ene, as indicated by NMR spectroscopy and X-ray crystallography (Scheme 3). The 31P{1H} NMR signal (40.9 ppm) is broad and the 11B{1H} NMR signal (30.5 ppm) is in the low-field region, even though the boron centers in compound 5 are four-coordinate. It is upfield shifted in comparison with that of the endocyclic boron center in VIII (50.1 ppm).12 The 29Si{1H} NMR resonance (233.0 ppm) is intermediate between cyclotrisilenylium cations (284.6–288.1 ppm)6 and bicyclo[1.1.0]tetrasil-1(3)-ene (217 ppm).21 The broad 29Si{1H} NMR signal of compound 5 is also downfield shifted in comparison with those of disilenes (50–155 ppm)26 and VIII (112.2 ppm). The molecular structure of 5 obtained by X-ray crystallography shows the planar fused bicyclo-Si2B2 ring being orthogonal to two N-phosphinoamidinate ligands (Fig. 2). It is expected that the exceptionally strained bridgehead SiSi double bond in the bicyclic scaffold distorts the structural parameters. The Si–B (2.034(4), 2.048(3) Å) and Si1–Si1A (2.3583(19) Å) bond lengths are longer than those in compound VIII (Si–B: 1.911(7)–1.952(3); Si–Si: 2.133(2)– 2.1469(11) Å).12 The Si1–Si1A bond is intermediate between those in disilenes (2.14–2.29 Å)26 and bicyclo[1.1.0]tetrasil-1(3)-ene (2.4873(10) Å).21 The Si–Si bond length in compound 5 is comparable with the Si–Si single bond length (ca. 2.3 Å). The UV-vis spectrum of compound 5 in THF (dark orange solution) shows an absorption band at λmax = 412 nm corresponding to the
(HOMO−2 → LUMO) and
(HOMO → LUMO+2) with equal contribution, and an intense absorption band at 474 nm corresponding to the
(HOMO−1 → LUMO).
![]() | ||
Scheme 3 Comparison of NMR signals (in red) and Si![]() |
The downfield 11B{1H} and 29Si{1H} NMR signals together with long SiSi double bond length indicate that the bridgehead Si
Si σ and π electrons in compound 5 could delocalize in the fused bicyclic scaffold. However, the multinuclear NMR signals and the distorted structural parameters of the Si2B2 ring cannot definitively confirm the presence of aromaticity. 1H NMR spectroscopy is a widely used experimental tool to assess whether a molecule is aromatic, although determining the degree of aromaticity remains a challenge.27 In the case of compound 5, there are no 1H NMR signals for the Si2B2 ring. Therefore, the degree of electron delocalization as σ- and π-aromaticity can only be evaluated through DFT calculations, which is discussed later in the manuscript.
The reaction of 4 with KC8 was performed in a less polar solvent, namely toluene, at room temperature (Scheme 2). The reaction was traced by 31P{1H} and 11B{1H} NMR spectroscopy after stirring for 16 h, where a mixture of 1,3-dibora-2,4-disilacyclobutanes 6a and 6b was observed (31P{1H} NMR: 102.4 and 108.4 ppm; 11B{1H} NMR: −11.3 and −34.8 ppm). They were isolated as a co-crystalline colorless solid (Fig. 3, yield: 10%) from the reaction mixture. The formation of 6a and 6b indicates that halogen scrambling occurred during the reduction. X-ray crystallography of 6a shows that the N-phosphinoamidinate ligands are in a P,N-chelate fashion, bridging across the B–Si bonds (Fig. 3). The B–Si bond lengths (2.033(10) and 2.025(11) Å) are comparable to those of 5. The Si1–N1 (1.804(7) Å) and P1–B1 bonds (1.979(10) Å) are typical single bonds. The C–N bond lengths (C1–N1: 1.416(10) and C1–N2: 1.281(10) Å) are unequal, showing that the exocyclic N2 atom is an imine moiety. When the reaction of 4 with KC8 in toluene was performed for 3 days, compound 5 was observed, indicating that compounds 6a and 6b are intermediates for the formation of compound 5. To support this, a mixture of compounds 6a and 6b were reacted with excess KC8 in toluene for 16 h, leading to the formation of compound 5.
From compounds 6a and 6b, it is plausible that the reduction of compound 4 with KC8 proceeds through the formation of compounds 6a and 6b first, followed by further reduction to form a diboradisilacyclobutadiene intermediate, which is anti-aromatic (Scheme S1†). Subsequently, it undergoes rearrangement by coordinating the N-phosphinoamidinate ligands with the boron centers to form the bridgehead SiSi double bond and bridging boron in compound 5.
To understand the nature of the Si–Si σ-orbital, the electron localization function (ELF) was calculated (Fig. 5a). The ELF plot of the fused bicyclic Si2B2 ring reveals four Si–B bond critical points (blue dots) and one Si2B2 ring critical point (orange dot). Notably, a weak σ-interaction is observed between the two bridgehead Si atom as reflected by the small ELF values (in green). It is because some of the Si–Si σ-electrons are delocalized with the B–Si σ-electrons (annular red region). Upon the removal of σ electrons from HOMO−2 (Fig. 5b), the distribution of energy density is altered, revealing a nearly non-existent σ bond between the two bridgehead Si atoms (in blue). This confirms the presence of a weak Si–Si σ bond which cannot be completely disregarded.
A quantum theory of atoms in molecules (QTAIM) analysis was conducted to quantify the Si–Si σ-interaction. QTAIM revealed the formation of 4 Si–B σ bonds, as seen from the bond critical point (BCP: blue dot 1–4, Fig. 6a) along each Si–B bond path. The mean Laplacian distribution ∇2ρ(rc) and energy density H(r) of the Si–B bonds at the BCP are 0.009 e a−5 and –0.060 hartree a−3, respectively. It should be noted that the π-electron density at the ring critical point (RCP) becomes zero, which is attributed to the RCP being located precisely on the nodal plane of the π orbital. As such, the values of ∇2ρ(rc) (0.008 e a−5) and H(r) (−0.032 hartree a−3) observed at the Si2B2 RCP (orange dot 5) indicate the presence of a distinct and weak σ bond between the two Si atoms. Upon the removal of σ electrons from HOMO−2 (Fig. 6b), the H(r) of the RCP reduced from −0.032 hartree a−3 to −0.017 hartree a−3, suggesting that HOMO−2 contributes to the weak σ bond between the two bridgehead Si atoms. These results are in line with ELF calculations. The feasibility of weak σ interaction in the strained and saturated fused bicyclic ring system has also been illustrated by Foroutan-Nejad's theoretical calculations.28
To understand the overall σ and π electronic delocalization in the entire Si2B2 ring, adaptive natural density partitioning (AdNDP) analysis was performed to elucidate the electronic structure of 5 in terms of the classical Lewis elements (lone-pairs and two-center-two-electron bonds) and delocalized n-center two-electron (nc-2e) bonds (Fig. 7) in order to identify possible delocalized electron-pair bonding.31 A simplified truncated model 5-H was used since all the bonding features of interest are located at the Si2B2 core in 5, which is preserved and clearly presented in the model. First, a 2-center-2-electron π-bond with occupation number (ON) = 1.41 e was found between two Si atoms (Fig. 7, top left). Such low occupancy indicates that the Si–Si π-bond is delocalized over the Si2B2 bicyclic ring, supporting the presence of 2π-aromaticity. Second, a 4-center-2-electron σ bond with ON = 1.99 e is found between two silicon atoms (Fig. 7, top right). In addition, four 4-center-2-electron σ bonds with ON = 1.88–1.99 e are observed at the Si2B2 ring (Fig. 7, bottom). These five 4-center-2-electron σ bonds indicate the presence of σ-aromatization,32 which is also illustrated by the contour plot (annular red region) in ELF calculations (Fig. 5a).
![]() | ||
Fig. 7 Adaptive natural population density (AdNDP) analysis of a simplified truncated model 5-H, where substituents (Dipp, Ph, and tBu) in compound 5 are substituted with hydrogen atoms for clarity. |
Anisotropy of the induced current density (ACID) was calculated to elucidate the σ- and π-aromaticity. The current density vectors plotted onto the ACID isosurface indicate a strong diatropic ring current in the π system above and underneath the ring, showing π-aromaticity in the Si2B2 ring (Fig. 8a). It is supported by the strong diatropic ring current in the HOMO (Fig. 8c), where the Si–Si π bond exhibits 2π-aromatic delocalization. Besides π delocalization, the current density vectors plotted onto the ACID isosurface indicate a diatropic ring current in the molecular plane mainly localized inside the four-membered ring (Fig. 8a). This is supported by the diatropic ring current in HOMO-1 (Fig. 8b), which illustrates the presence of σ-tangential aromatic delocalization33 among the Si–B bonds. It is noteworthy that the diatropic ring current is restricted to the periphery. A similar observation can be found in naphthalene.34 In addition, the positive natural population analysis (NPA) charge on the bridgehead Si atoms (0.36e) and the negative NPA charge on the B atoms (−0.55e) suggest that the electron density in the ring flows from the silicon to boron atoms. The highly negative Nucleus Independent Chemical Shift (NICS) values at the center of each three-membered ring [NICS(0): −38.7, NICS(1): −13.8, NICS(1)zz: −22.4 ppm] and at the center of the whole Si2B2 four-membered ring [NICS(0): −38.5, NICS(1): −22.1, NICS(1)zz: −34.0 ppm] indicate considerable inorganic aromaticity in the fused bicyclic 2π-Si2B2 ring.
To further quantify the aromaticity of the B2Si2 four-membered ring in compound 5, the electron density of delocalized bonds (EDDB) was employed. The EDDB value for compound 5 is 1.901 e. The EDDB value comprises the EDDBσ component (1.193 e) and the EDDBπ component (0.708 e) within the B2Si2 four-membered ring of compound 5. These results illustrate that compound 5 possesses both σ- and π-aromaticity in the Si2B2 ring of compound 5. The EDDBπ component (0.708 e) is comparable with that of VI and related 2π-aromatic Si4 ring analogues (0.586–0.979 e, CAM-B3LYP/def2-TZVP).10 Other ring systems containing both σ and π aromatic delocalization have been reported by Saito and Berndt.35
5-NMe2 with two -NMe2 substituted tricoordinate boron centers was used as a model molecule (Fig. 9), and its NICS and EDDB calculations were performed for comparison with compound 5. Similar to compound 5, the 5-NMe2 model possesses both σ- and π-aromatic delocalization. The HOMO of the 5-NMe2 model is a π orbital delocalized over the entire fused bicyclic Si2B2 ring. It can be considered as the overlapping of pπ orbitals on the Si atoms with the B–N π* orbitals. The HOMO-1 of the 5-NMe2 model represents the σ-interaction between two bridgehead silicon atoms, which significantly extends over the fused bicyclic Si2B2 ring. The NICS value at the center of the Si2B2 ring [NICS(0): −37.5, NICS(1): −23.2, NICS(1)zz: −35.8 ppm], the EDDBσ value (2.071 e) and the EDDBπ value (1.396 e) illustrate the σ- and π-aromaticity in the Si2B2 ring of the 5-NMe2 model. The higher EDDBπ value of 5-NMe2 is attributed to the empty p orbital on the boron center which allows for more effective overlapping with the SiSi π orbital. The lower EDDBπ value in compound 5 is due to the less effective overlapping of pπ orbitals on the Si atoms with the B–N and B–P σ* orbitals, having pseudo π symmetry at the four-coordinate boron atoms.
Footnotes |
† Electronic supplementary information (ESI) available: Experimental procedures, complete characterization data, computational details, and NMR spectra. CCDC 2266104–2266108 and 2380825 for 2–7. For ESI and crystallographic data in CIF or other electronic format see DOI: https://doi.org/10.1039/d4sc05867d |
‡ These authors contributed equally. |
This journal is © The Royal Society of Chemistry 2025 |