Si3C2H2 isomers with a planar tetracoordinate carbon or silicon atom(s)

Ligand-enforced geometric constraints and associated reactivity in p-block compounds

The geometry at an element centre can generally be predicted based on the number of electron pairs around it using valence shell electron pair repulsion (VSEPR) theory. Strategies to distort p-block compounds away from these predicted geometries have gained considerable interest due to the unique structural outcomes, spectroscopic properties or reactivity patterns engendered by such distortion. This review presents an up-to-date group-wise summary of this exciting and rapidly growing field with a focus on understanding how the ligand employed unlocks structural features, which in turn influences the associated reactivity. Relevant geometrically constrained compounds from groups 13-16 are discussed, along with selected stoichiometric and catalytic reactions. Several areas for advancement in this field are also discussed. Collectively, this review advances the notion of geometric tuning as an important lever, alongside electronic and steric tuning, in controlling bonding and reactivity at p-block centres.

Theoretical Studies on the Isomerization Kinetics of Low-Lying Isomers of the SiC4H2 System

The low-lying isomers of SiC4H2 are investigated to understand the kinetics of isomerization pathways using density functional theory. In our earlier work, we studied the various possible isomers (J. Phys. Chem. A, 2020, 124, 987-1002) and the chemical bonding of low-lying isomers of SiC4H2 (J. Phys. Chem. A, 2022, 126, 9366-9374). Among them, four isomers, 1-ethynyl-3-silacycloprop-1-en-3-ylidene (1), 3-silapent-1,4-diyn-3-ylidene (2), 1-silapent-1,2,3,4-tetraen-1-ylidene (4), and 1-silapent-2,4-diyn-1-ylidene (5) have already been identified in the laboratory. The previously known theoretical isomer 2-methylene-1-silabicyclo[1.1.0]but-1(3)-en-4-ylidene (3) and the newly identified unknown isomer through the present kinetic studies 5-silabicyclo[2.1.0]pent-1(4),2-dien-5-ylidene (N6) remain elusive in the laboratory to date. The isomerization pathways of the low-lying isomers of SiC4H2 are predicted through the transition state structures. Intrinsic reaction coordinate analysis identifies the minimum energy reaction pathways connecting the transition state from one isomer to another of the investigated system. The present kinetic data reveal the isomerization of global minimum energy isomer 1 to thermodynamically stable low-lying isomers, 2 and 5. Interestingly, isomer 3 interconverts to the experimentally known low-energy isomer 4, the second most thermodynamically stable isomer among them. The thermodynamic and kinetic parameters of the low-lying isomers of SiC4H2 are also documented in this work. The rate coefficient and equilibrium constant for isomerization reactions are calculated using the Rice-Ramsperger-Kassel-Marcus theory. The equilibrium constant delineates the difficulties in forming N6 and 3 through the isomerization pathways. Furthermore, ab initio molecular dynamics studies dictate the stability of low-lying isomers of SiC4H2 within the time scale of the simulation.

CCSD(T) Rotational Constants for Highly Challenging C5H2 Isomers-A Comparison between Theory and Experiment

We evaluate the accuracy of CCSD(T) and density functional theory (DFT) methods for the calculation of equilibrium rotational constants (Ae, Be, and Ce) for four experimentally detected low-lying C5H2 isomers (ethynylcyclopropenylidene (2), pentatetraenylidene (3), ethynylpropadienylidene (5), and 2-cyclopropen-1-ylidenethenylidene (8)). The calculated rotational constants are compared to semi-experimental rotational constants obtained by converting the vibrationally averaged experimental rotational constants (A0, B0, and C0) to equilibrium values by subtracting the vibrational contributions (calculated at the B3LYP/jun-cc-pVTZ level of the theory). The considered isomers are closed-shell carbenes, with cumulene, acetylene, or strained cyclopropene moieties, and are therefore highly challenging from an electronic structure point of view. We consider both frozen-core and all-electron CCSD(T) calculations, as well as a range of DFT methods. We find that calculating the equilibrium rotational constants of these C5H2 isomers is a difficult task, even at the CCSD(T) level. For example, at the all-electron CCSD(T)/cc-pwCVTZ level of the theory, we obtain percentage errors ≤0.4% (Ce of isomer 3, Be and Ce of isomer 5, and Be of isomer 8) and 0.9-1.5% (Be and Ce of isomer 2, Ae of isomer 5, and Ce of isomer 8), whereas for the Ae rotational constant of isomers 2 and 8 and Be rotational constant of isomer 3, high percentage errors above 3% are obtained. These results highlight the challenges associated with calculating accurate rotational constants for isomers with highly challenging electronic structures, which is further complicated by the need to convert vibrationally averaged experimental rotational constants to equilibrium values. We use our best CCSD(T) rotational constants (namely, ae-CCSD(T)/cc-pwCVTZ for isomers 2 and 5, and ae-CCSD(T)/cc-pCVQZ for isomers 3 and 8) to evaluate the performance of DFT methods across the rungs of Jacob's Ladder. We find that the considered pure functionals (BLYP-D3BJ, PBE-D3BJ, and TPSS-D3BJ) perform significantly better than the global and range-separated hybrid functionals. The double-hybrid DSD-PBEP86-D3BJ method shows the best overall performance, with percentage errors below 0.5% in nearly all cases.

Chemical Bonding Perspective on Low-Lying SiC4H2 Isomers: Conceptual Quantum Chemical Views

The nature of the chemical bonding in seven low-lying isomers of SiC₄H₂ is analyzed through quantum chemical concepts. Out of the seven, four isomers, 1-ethynyl-3-silacycloprop-1(2)-en-3-ylidene (1), diethynylsilylidene (2), 1-sila-1,2,3,4-pentatetraenylidene (4), and 1,3-butadiynylsilylidene (5), have already been identified in the laboratory. The other three isomers, 2-methylenesilabicyclo[1.1.0]but-1(3)-en-4-ylidene (3), 4-sila-2-methylenebicyclo[1.1.0]but-1(3)-en-4-ylidene (6), and 3-ethynyl-1-silapropadienylidene (7) remain elusive in the laboratory to date (J. Phys. Chem. A, 2020, 124, 987-1002). Deep insight into the characteristics of chemical bonding is explored with different bonding analysis tools. Quantum theory of atoms in molecules (QTAIM), interaction quantum atoms analysis, natural bond orbital analysis, adaptive natural density partitioning, electron localization function (ELF), Laplacian of electron density, energy decomposition analysis, atomic charge analysis, bond order analysis, and frontier molecular orbital analysis are employed in the present work to gain a better understanding of the chemical bonding perspective in SiC₄H₂ isomers. Different quantum chemical topology approaches (QTAIM, ELF, and Laplacian of electron density) are employed to complement each other. The obtained results dictate that the lone pair of the silicon atom participate in delocalization and influences the structural stability of isomers.

Pentacoordinate Carbon Atoms in a Ferrocene Dication Derivative—[Fe(Si2-η5-C5H2)2]2+

Pentacoordinate carbon atoms are theoretically predicted here in a ferrocene dication derivative in the eclipsed-(1; C2v), gauche-(2; C2) and staggered-[Fe(Si2-η5-C5H2)2]2+(3; C2h) forms for the first time. Energetically, the relative energy gaps for 2 and 3 range from −3.06 to 16.74 and −2.78 to 40.34 kJ mol−1, respectively, when compared to the singlet electronic state of 1 at different levels. The planar tetracoordinate carbon (ptC) atom in the ligand Si2C5H2 becomes a pentacoordinate carbon upon complexation. The ligand with a ptC atom was predicted to be both a thermodynamically and kinetically stable molecule by some of us in our earlier theoretical works. Natural bond orbital and adaptive natural density partitioning analyses confirm the pentacoordinate nature of carbon in these three complexes (1–3). Although they are hypothetical at the moment, they support the idea of “hypercoordinate metallocenes” within organometallic chemistry. Moreover, ab initio molecular dynamics simulations carried out at 298 K temperature for 2000 fs suggest that these molecules are kinetically stable.

From Molecules with a Planar Tetracoordinate Carbon to an Astronomically Known C5H2 Carbene

Ethynylcyclopropenylidene (2), an isomer of C₅H₂, is a known molecule in the laboratory and has recently been identified in Taurus Molecular Cloud-1 (TMC-1). Using high-level coupled-cluster methods up to the CCSDT(Q)/CBS level of theory, it is shown that two isomers of C₅H₂ with a planar tetracoordinate carbon (ptC) atom, (SP-4)-spiro[2.2]pent-1,4-dien-1,4-diyl (11) and (SP-4)-spiro[2.2]pent-1,4-dien-1,5-diyl (13), serve as the reactive intermediates for the formation of 2. Here, a theoretical connection has been established between molecules containing ptC atoms (11 and 13) and a molecule (2) that is present nearly 430 light years away, thus providing evidence for the existence of ptC species in the interstellar medium. The reaction pathways connecting the transition states and the reactants and products have been confirmed by intrinsic reaction coordinate calculations at the CCSDT(Q)/CBS//B3LYP-D3BJ/cc-pVTZ level. While isomer 11 is non-polar (μ = 0), isomers 2 and 13 are polar, with dipole moment values of 3.52 and 5.17 Debye at the CCSD(T)/cc-pVTZ level. Therefore, 13 is also a suitable candidate for both laboratory and radioastronomical studies.

The inversion of tetrahedral p-block element compounds: general trends and the relation to the second-order Jahn-Teller effect

The tetrahedron is the primary structural motif among the p-block elements and determines the architecture of our bio- and geosphere. However, a broad understanding of the configurational inversion of tetrahedral compounds is missing. Here, we report over 250 energies (DLPNO-CCSD(T)) for square planar inversion of third- and fourth-period element species of groups 13, 14, and 15. Surprisingly low inversion barriers are identified for compounds of industrial relevance (e.g., ≈100 kJ mol-1 for Al(OH)4 -). More fundamentally, the second-order Jahn-Teller theorem is disclosed as suitable to rationalize substituent and central element effects. Bond analysis tools give further insights into the preference of eight valence electron systems with four substituents to be tetrahedral. Hence, this study develops a model to understand, memorize, and predict the angular flexibility of tetrahedral species. Perceiving the tetrahedron not as forcingly rigid but as a dynamic structural entity might leverage new approaches and visions for adaptive matter.

BAl4Mg−/0/+: Global Minima with a Planar Tetracoordinate or Hypercoordinate Boron Atom

We have explored the chemical space of BAl4Mg−/0/+ for the first time and theoretically characterized several isomers with interesting bonding patterns. We have used chemical intuition and a cluster building method based on the tabu-search algorithm implemented in the Python program for aggregation and reaction (PyAR) to obtain the maximum number of possible stationary points. The global minimum geometries for the anion (1a) and cation (1c) contain a planar tetracoordinate boron (ptB) atom, whereas the global minimum geometry for the neutral (1n) exhibits a planar pentacoordinate boron (ppB) atom. The low-lying isomers of the anion (2a) and cation (3c) also contain a ppB atom. The low-lying isomer of the neutral (2n) exhibits a ptB atom. Ab initio molecular dynamics simulations carried out at 298 K for 2000 fs suggest that all isomers are kinetically stable, except the cation 3c. Simulations carried out at low temperatures (100 and 200 K) for 2000 fs predict that even 3c is kinetically stable, which contains a ppB atom. Various bonding analyses (NBO, AdNDP, AIM, etc.) are carried out for these six different geometries of BAl4Mg−/0/+ to understand the bonding patterns. Based on these results, we conclude that ptB/ppB scenarios are prevalent in these systems. Compared to the carbon counter-part, CAl4Mg−, here the anion (BAl4Mg−) obeys the 18 valence electron rule, as B has one electron fewer than C. However, the neutral and cation species break the rule with 17 and 16 valence electrons, respectively. The electron affinity (EA) of BAl4Mg is slightly higher (2.15 eV) than the electron affinity of CAl4Mg (2.05 eV). Based on the EA value, it is believed that these molecules can be identified in the gas phase. All the ptB/ppB isomers exhibit π/σ double aromaticity. Energy decomposition analysis predicts that the interaction between BAl4−/0/+ and Mg is ionic in all these six systems.

Stereomutation in Tetracoordinate Centers via Stabilization of Planar Tetracoordinated Systems

The quest for stabilizing planar forms of tetracoordinate carbon started five decades ago and intends to achieve interconversion between [R]- and [S]-stereoisomers without breaking covalent bonds. Several strategies are successful in making the planar tetracoordinate form a minimum on its potential energy surface. However, the first examples of systems where stereomutation is possible were reported only recently. In this study, the possibility of neutral and dications of simple hydrocarbons (cyclopentane, cyclopentene, spiropentane, and spiropentadiene) and their counterparts with the central carbon atom replaced by elements from groups 13, 14, and 15 are explored using ab initio MP2 calculations. The energy difference between the tetrahedral and planar forms decreases from row II to row III or IV substituents. Additionally, aromaticity involving the delocalization of the lone pair on the central atom appears to help in further stabilizing the planar form compared to the tetrahedral form, especially for the row II substituents. We identified 11 systems where the tetrahedral state is a minimum on the potential energy surface, and the planar form is a transition state corresponding to stereomutation. Interestingly, the planar structures of three systems were found to be minimum, and the corresponding tetrahedral states were transition states. The energy profiles corresponding to such transitions involving both planar and tetrahedral states without the breaking of covalent bonds were examined. The systems showcased in this study and research in this direction are expected to realize molecules that experimentally exhibit stereomutation.

CAl4Mg0/−: Global Minima with a Planar Tetracoordinate Carbon Atom

Isomers of CAl4Mg and CAl4Mg− have been theoretically characterized for the first time. The most stable isomer for both the neutral and anion contain a planar tetracoordinate carbon (ptC) atom. Unlike the isovalent CAl4Be case, which contains a planar pentacoordinate carbon atom as the global minimum geometry, replacing beryllium with magnesium makes the ptC isomer the global minimum due to increased ionic radii of magnesium. However, it is relatively easier to conduct experimental studies for CAl4Mg0/− as beryllium is toxic. While the neutral molecule containing the ptC atom follows the 18 valence electron rule, the anion breaks the rule with 19 valence electrons. The electron affinity of CAl4Mg is in the range of 1.96–2.05 eV. Both the global minima exhibit π/σ double aromaticity. Ab initio molecular dynamics simulations were carried out for both the global minima at 298 K for 10 ps to confirm their kinetic stability.

Planar Tetracoordinated Silicon (ptSi): Room-Temperature Stable Compounds Containing Anti-van't Hoff/Le Bel Silicon

While a variety of compounds containing planar tetracoordinated carbon (ptC), the so-called anti-van't Hoff/Le Bel carbon, are known experimentally, stable systems containing planar tetracoordinated silicon (ptSi) are barely known. As part of our studies on the application of stereoelectronically well-defined transition-metal fragments to stabilize silicon in unprecedented bonding modes, we report herein the synthesis and full characterization of a series of thermally stable complexes of the general formula [Tp'(CO)₂MSiC(R¹)C(R²)M(CO)₂Tp'] (M = Mo, W; R¹ = R² = Me or R¹ = H, R² = SiMe₃, Ph; Tp' = κ³-N,N',N″-hydridotris(3,5-dimethylpyrazolyl)borate), which incorporate a ptSi atom in addition to two ptC atoms. The complexes were obtained by reacting the metallasilylidyne complexes [Tp'(CO)₂M≡Si-M(CO)₂(PMe₃)Tp'] with alkynes R¹C≡CR² and were comprehensively analyzed by experimental studies and quantum chemical calculations. The analyses revealed that the ptSi atom is embedded in a tricyclic trapezoidal core featuring one internal SiC₂ and two outer M-Si-C three-membered rings, which are fused via two Si-C bonds. The structural peculiarities evoked by the presence of an anti-van't Hoff/Le Bel ptSi center, such as the short M-Si bonds, a nearly linear M-Si-M spine, long M-C bonds, and the presence of two planar tetracoordinated carbon atoms were elucidated by a detailed analysis of the electronic structure, suggesting that one factor for the stabilization of the ptSi geometry is the aromaticity of the central SiC₂ ring having two delocalized π electrons. Remarkably, the results further indicate the existence of both anti-van't Hoff/Le Bel carbon and silicon centers next to each other in the isolated complexes.

Kinetic Stability of Si2C5H2 Isomer with a Planar Tetracoordinate Carbon Atom

Dissociation pathways of the global minimum geometry of Si2C5H2 with a planar tetracoordinate carbon (ptC) atom, 2,7-disilatricyclo[4.1.0.01,3]hept-2,4,6-trien-2,7-diyl (1), have been theoretically investigated using density functional theory and coupled-cluster (CC) methods. Dissociation of Si-C bond connected to the ptC atom leads to the formation of 4,7-disilabicyclo[4.1.0]hept-1(6),4(5)-dien-2-yn-7-ylidene (4) through a single transition state. Dissociation of C-C bond connected to the ptC atom leads to an intermediate with two identical transition states and leads back to 1 itself. Simultaneous breaking of both Si-C and C-C bonds leads to an acyclic transition state, which forms an acyclic product, cis-1,7-disilahept-1,2,3,5,6-pentaen-1,7-diylidene (19). Overall, two different products, four transition states, and an intermediate have been identified at the B3LYP/6-311++G(2d,2p) level of theory. Intrinsic reaction coordinate calculations have also been done at the latter level to confirm the isomerization pathways. CC calculations have been done at the CCSD(T)/cc-pVTZ level of theory for all minima. Importantly, all reaction profiles for 1 are found be endothermic in Si2C5H2. These results are in stark contrast compared to the structurally similar and isovalent lowest-energy isomer of C7H2 with a ptC atom as the overall reaction profiles there have been found to be exothermic. The activation energies for Si-C, C-C, and Si-C/C-C breaking are found to be 30.51, 64.05, and 61.85 kcal mol−1, respectively. Thus, it is emphasized here that 1 is a kinetically stable molecule. However, it remains elusive in the laboratory to date. Therefore, energetic and spectroscopic parameters have been documented here, which may be of relevance to molecular spectroscopists in identifying this key anti-van’t-Hoff-Le Bel molecule.

How to Accomplish a Square C(N)4 Substructure of the Planar Tetracoordinate Carbon

Nitrogen-based groups are usually not used as ligands to coordinate to the ptC atom. However, here we reported only nitrogen-based ligands to accomplish a theoretically successful square planar C(N)4 substructure. The first difficulty in accomplishing a square ptC(N)4 substructure is to conquer the tremendous strain from the planar to tetrahedral arrangements, and the second is to restrict it in a suitable system with the right symmetry. We designed several neutral molecules with the square ptC(N)4 substructures, and the molecules were studied using the density functional theory method at the B3LYP/6-311++G(3df,3pd) and TPSSh/6-311++G(3df,3pd) level of theory. The results of this work show that the molecules are all real minima on the potential energy surface and successfully achieved the square ptC(N)4 substructure in the theoretical method. The group orbitals among the square ptC(N)4 arrangement in the D 2d symmetry have been discussed and used to investigate the bonding interactions among all atoms in the square ptC(N)4 substructure. Usually, the ptC systems have 18 valence electrons, but the present ptC systems mentioned in this work have 24 valence electrons, which is unusual for ptC.

MgC6H2 Isomers: Potential Candidates for Laboratory and Radioastronomical Studies

Eighty three stationary points of MgC₆H₂ isomers spanning from 0 to 215 kcal mol^-1 have been theoretically identified using density functional theory at the B3LYP/6-311++G(2d,2p) level of theory. Among them, four low-lying isomers lying within 23.06 kcal mol^-1 (1 eV) have been further characterized in detail using high-level coupled-cluster (CC) methods. The thermodynamically most stable isomer turns out to be 1-magnesacyclohepta-4-en-2,6-diyne (1). The other three isomers, 3-magnesahepta-1,4,6-triyne (2), 1-magnesacyclohepta-2,3,4-trien-6-yne (3), and 1-magnesahepta-2,4,6-triyne (4) lie 8.24, 19.76, and 21.36 kcal mol^-1, respectively, above 1 at the ae-CCSD(T)/cc-pCVTZ level of theory. All the four isomers are polar with a permanent electric dipole moment (μ ≠ 0). Hence, they are potential candidates for rotational spectroscopic studies. Considering the recent identification of magnesium-bearing hydrocarbons such as, MgC₂H and MgC₄H in IRC+10216, it is believed that the current theoretical data may be of relevance to laboratory molecular spectroscopic and radioastronomical studies on MgC₆H₂ isomers. The energetic and spectroscopic information gathered in this study would aid the detection of low-lying MgC₆H₂ isomers in the laboratory, which are indispensable for radioastronomical studies. It is also noted here that neither the National Institute of Standards and Technology Chemistry WebBook nor the Kinetic Database for Astrochemistry lists any isomer of MgC₆H₂ at the moment. Therefore, these isomers are studied here theoretically for the very first time.

Si2C5H2 isomers - search algorithms versus chemical intuition

The pros and cons of using search algorithms alone in identifying new geometries have been discussed by using the Si₂C₅H₂ elemental composition as an example. Within 30 kcal mol^-1 at the CCSD(T)/def2-TZVP//PBE0/def2-TZVP level of theory, the coalescence kick and cuckoo methods postulate merely four isomers (1, 3, 6, and 7) for Si₂C₅H₂ (O. Yañez et. al., Chem. Commun., 2017, 53, 12112). On the contrary, chemical intuition yields fourteen (2, 4, 5, and 8-18) new isomers within the same energy range at the B3LYP/6-311++G(2d,2p) level of theory. Based on the relative energies of the first eleven isomers of Si₂C₅H₂ (1, C_2v, 0.00; 2, C_s, 21.39; 3, C_s, 21.95; 4, C_s, 22.76; 5, C_s, 24.74; 6, C_s, 25.34; 7, C_s, 25.64; 8, C_s, 25.79; 9, C_s, 27.20; 10, C_2v, 28.59; and 11, C_2v, 29.16 kcal mol^-1) calculated at the CCSD(T)/cc-pVTZ level of theory, it is evident that the search algorithms had missed at least seven isomers in the same energy range. The relative energy gaps of isomers 12-18 fall in the range of 30-40 kcal mol^-1 at the latter level of theory. Consequentially, this scenario triggers a speculation going forward with search algorithms alone in the search of all new isomers. While one cannot underestimate the power of these algorithms, the role of chemical intuition may not be completely neglected. Retrospectively, the fourteen new isomers found by chemical intuition may help in writing better search algorithms. All eighteen isomers - including the most stable isomer with a planar tetracoordinate carbon atom 1- remain elusive in the laboratory to date. Thus, structural and spectroscopic parameters have been presented here, which may possibly aid the future experimental studies.

Theoretical Studies of SiC4H2 Isomers Delineate Three Low-Lying Silylidenes Are Missing in the Laboratory

Eleven isomers of SiC₄H₂ lying within 50 kcal mol^-1 have been theoretically investigated using density functional theory and high-level coupled-cluster methods. Among them, four isomers, 1-ethynyl-3-silacycloprop-1(2)-en-3-ylidene (1), diethynylsilylidene (2), 1-sila-1,2,3,4-pentatetraenylidene (4), and 1,3-butadiynylsilylidene (5), have already been identified in the laboratory. The current investigation reports three low-lying (<1 eV) silylidenes [2-methylenesilabicyclo[1.1.0]but-1(3)-en-4-ylidene (3), 4-sila-2-methylenebicyclo[1.1.0]but-1(3)-en-4-ylidene (6), and 3-ethynyl-1-silapropadienylidene (7)] and three high-lying (>1 eV) silylidenes [2-sila-(didehydrovinylidene)cyclopropene (8), an isomer with a planar tetracoordinate carbon (ptC) atom (10), and 1-ethynyl-1-silapropadienylidene (11)], which remain elusive in the laboratory to date. Isomer 9 also contains a ptC atom, which turned out to be a transition state at all levels. Though all isomers are polar (μ ≠ 0), rotational spectrum is available only for 4. Using matrix isolation, three isomers (1, 2, and 5) have been trapped in the laboratory at 10 K. Considering the astrochemical relevance of silicon-carbide clusters in the interstellar medium, the current theoretical data demand new molecular spectroscopic studies on SiC₄H₂. Surprisingly, unlike the isovalent C₅H₂ isomers, where the bent carbenes are yet to be identified in the laboratory, the bent silylidenes (2 and 5) have been trapped in the case of SiC₄H₂. In both the cases, molecules with transannular C-C and/or Si-C bonds remain elusive, though they lie in the low-lying region. Using suitable precursors, whether these peculiar geometries (especially 3 and 6) would be identified or not in the laboratory needs to be addressed by molecular spectroscopists. The present investigation documents structural and spectroscopic information of SiC₄H₂ isomers, which may compliment future molecular spectroscopic observations including radioastronomical searches.