Structure and bonding in mixed silicon–carbon clusters and their anions

Theoretical Rovibrational Spectroscopy of Magnesium Tricarbide-Multireference Character Thwarts a Full Analysis of All Isomers

Magnesium tricarbide isomers are studied herein with coupled cluster theory and multireference configuration interaction to support their possible detection in astrochemical environments such as the circumstellar envelope surrounding the star IRC +10216 or in terrestrial laboratories. Magnesium-bearing species may abound in the interstellar medium (ISM), but only eight (MgNC, MgCN, HMgNC, MgC₂H, MgC₃N, MgC₄H, MgC₅N, and MgC₆H) have been directly identified thus far. Several possible isomers for the related MgC₃ system are explored in their singlet and triplet spin multiplicities. Overall, this work offers quantum chemical insight of rovibrational spectroscopic data for MgC₃ using quartic force fields (QFFs) based on the CCSD(T) and CCSD(T)-F12 levels of theory at the complete basis set (CBS) limit. Additional corrections with small basis set CCSDT(Q) and scalar relativistic effects are also included in the analysis. Salient multireference character is found in the singlet diamond electronic state, which makes a definitive assignment of the ground state challenging. Nevertheless, coupled cluster-based composite energies and multireference configuration interaction both predict that the ¹A₁ diamond isomer is 1.6-2.2 kcal mol^-1 lower in energy than the ³A₁ diamond isomer. Furthermore, highly accurate binding energies of various isomers MgC₃ are provided for comparison to photodetachment experiments. Dipole moments along with harmonic infrared intensities will guide efforts for astronomical and spectroscopic characterization.

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.

A template for a planar tetracoordinate heavier group 14 atom: a global study of C2Si2X(q) (X = C, Si, Ge, Sn, Pb; q = +1, 0, -1)

Through a global isomeric study, we computationally identified the first structural template C2Si2X that could encompass a planar tetracoordinate X for all the heavier group 14 elements X in the 0, +1 or -1 charge state. We thus significantly expanded the traditional 16/17/18ve rules to 19/20/21ve for ptX.

Discovery of a Missing Link: Detection and Structure of the Elusive Disilicon Carbide Cluster

The rotational spectrum of the elusive but fundamentally important silicon carbide SiCSi has been detected using sensitive microwave techniques aided by high-level ab initio methods. Its equilibrium structure has been determined to very high precision using isotopic substitution and vibrational corrections calculated quantum-chemically: it is an isosceles triangle with a Si-C bond length of 1.693(1) Å, and an apex angle of 114.87(5)°. Now that all four Si(m)C(n) clusters with m + n = 3 have been observed experimentally, their structure and chemical bonding can be rigorously compared. Because Si2C is so closely linked to other Si-bearing molecules that have been detected in the evolved carbon star IRC+10216, it is an extremely promising candidate for detection with radio telescopes.

Vibrational spectra and structures of SinC clusters (n = 3–8)

The geometries of C-doped silicon clusters determined from infrared spectroscopy and computational chemistry reveal the stable Si₃C unit as a common structural motif.

Smallest deltahedra silicon dicarbide: C2Si32−

In this paper, we reported the smallest main-group dicarbide with all deltahedras, which is also the first main-group dicarbide with (n + 1) polyhedral skeletal electron pairs (PSEPs).

An ab initio calculation study of silicon and carbon binary clusters C7Si(n) (n = 1-7)

Binary C(7)Si(n) (n = 1-7) clusters are studied using density functional calculations at the level of B3LYP/6-311G(d). Lowest-energy structures have been determined theoretically and their properties such as binding energies, second differences in energy and highest-occupied and lowest-unoccupied molecular orbital gaps have been analyzed. It is found that the lowest-energy structures of the C(7)Si(n) (n = 1-7) clusters change from linear to planar when n ≥ 3, and in the planar structures C atoms prefer to form five- and six-membered rings surrounded by extra Si atoms in the form of the C(2)Si units.

Unusual isomers of disilacyclopropenylidene (Si2CH2)

Nine electronic singlet state structures of Si(2)CH(2) have been systematically investigated by high level theoretical methods. This research employed coupled cluster (CC) methods with single and double excitations (CCSD) and CCSD with perturbative triple excitations [CCSD(T)] using the correlation-consistent polarized valence cc-pVXZ/cc-pV(X+d)Z (X = D, T, and Q) basis sets. Full valence complete active space self-consistent-field (CASSCF) wave functions were used for the interpretation of geometries and physical properties. Among the nine singlet stationary points, six structures (1S-6S) are found to be minima, two structures (7S and 8S) are transition states, and one structure (9S) is a second-order saddle point. The existence of the two peculiar hydrogen bridged isomers, 1S (Si...H...Si) and 4S (agostic CH...Si) is established. Extensive focal point analyses are used to obtain complete basis set (CBS) limit energies. For the six lowest-lying singlet minima, after focal point analyses, the energy ordering and energy differences (in kcal mol(-1), with the zero-point vibrational energy corrected values in parentheses) are predicted to be 1S [0.0 (0.0)] < 3S [14.7 (14.5)] < 4S [25.1 (25.3)] < 5S [28.2 (26.0)] < 6S [45.0 (45.4)] < 2S [73.8 (72.0)]. Their relative energies are strikingly different from those for the isovalent parent C(3)H(2) molecule. Geometries, dipole moments, harmonic vibrational frequencies, and associated infrared (IR) intensities are reported for all equilibrium structures.

A theoretical study of the cyclization processes of energized CCCSi and CCCP

Calculations at the CCSD(T)/aug-cc-pVDZ//B3LYP/6-31+G(d) level of theory have shown that cyclization of both the ground state triplet and the corresponding singlet state of CCCSi may rearrange to give cyclic isomers which upon ring opening may reform linear C(3)Si isomers in which the carbon atoms are scrambled. The cyclization processes are energetically favorable with barriers to the transition states from 13 to 16 kcal mol(-1). This should be contrasted with the analogous process of triplet CCCC to triplet rhombic C(4), which requires an excess energy of 25.8 kcal mol(-1). A similar cyclization of doublet CCCP requires 50.4 kcal mol(-1) of excess energy; this should be contrasted with the same process for CCCN, which requires 54.7 kcal mol(-1) to effect cyclization.

A new class of silicon-carbon clusters: a full study of the hydrogenated SinC2H2, n=3,4,5, clusters in comparison with their isoelectronic carboranes C2BnHn+2

The structural and electronic characteristics of the Si(n)C(2)H(2), n=3,4,5, clusters are studied by ab initio calculations based on coupled cluster and density functional theory using the hybrid B3LYP functional. It is demonstrated that all three clusters are structurally and electronically homologous to the corresponding isoelectronic organometallic carboranes C(2)B(n)H(n+2). This homology, which is in full agreement with the analogy of Si(6) (2-) and B(6)H(6) (2-) demonstrated recently by the author [J. Chem. Phys. 127, 014314 (2007)], includes not only the ground states but also the lower-lying isomers as well. These lowest lying isomers can be obtained by ortho, para, and meta substitutions from the corresponding Si(n) (2-), n=3,4,5, dianions. The energetic ordering of the low-lying isomers is in full agreement with the known valence and topological charge stability rules developed for carboranes. The hydrogenated clusters are much more stable than their nonhydrogenated counterparts. It is suggested that Si(3)C(2)H(2), Si(4)C(2)H(2), and Si(5)C(2)H(2), which can be probably found in interstellar space, are special examples of a general class of silicon-carbon clusters of the form Si(n)C(2)H(2), with analogous properties and similarities to the corresponding carboranes C(2)B(n)H(n+2). It is furthermore illustrated that the lowest energy structures of the Si(n)C(2) clusters can be obtained through a systematic and straightforward procedure from the Si(n)C(2)H(2) clusters. The present results could hopefully make possible the exploitation of the rich borane and carborane chemistry for the design and development of novel silicon and silicon-carbon composite nanomaterials.

A Density Functional Study on Beryllium-Doped Carbon Dianion Clusters CnBe2- (n = 4−14)

Making use of the software of molecular graphics, we designed numerous models of C(n)()Be(2-) (n = 4-14). We carried out geometry optimization and calculation on vibration frequency by means of the B3LYP density functional method. After comparison of structure stability, we found that the ground-state isomers of C(n)()Be(2-) (n = 4-14) are linear with the beryllium atom located inside the C(n)() chain. When a side carbon chain is with an even number of carbon atoms, it is polyacetylene-like, whereas when a side chain is with an odd number of carbon atoms, it is cumulene-like. The C(n)Be(2-) (n = 4-14) clusters with an even number of carbon atoms are more stable than that with an odd number of carbon atoms, matching the peak pattern observed in accelerator mass spectrometry (AMS) and Coulomb Explosion Imaging (CEI) investigations of C(n)()Be(2-) (n = 4-14). The trend of such odd/even alternation is explained based on concepts of bonding characteristics, electronic configuration, electron detachment, and incremental binding energy.

A systematic multireference perturbation-theory study of the low-lying states of SiC3

The three known lowest-energy isomers of SiC(3), two cyclic singlets (2s and 3s) and a linear triplet (1t), have been reinvestigated using multireference second-order perturbation theory (MRPT2). The dependence of the relative energies of the isomers upon the quality of the basis sets and the sizes of the reference active spaces is explored. When using a complete-active-space self-consistent-field reference wave function with 12 electrons in 11 orbitals [CASSCF (12, 11)] together with basis sets that increase in size up to the correlation-consistent polarized core-valence quadruple zeta basis set (cc-pCVQZ), the MRPT2 method consistently predicts the linear triplet to be the most stable isomer. A new parallel direct determinant MRPT2 code has been used to systematically explore reference spaces that vary in size from CASSCF (8,8) to full optimized reaction space [FORS or CASSCF (16,16)] with the cc-pCVQZ basis. It is found that the relative energies of the isomers change substantially as the active space is increased. At the best level of theory, MRPT2 with a full valence FORS reference, the 2s isomer is predicted to be more stable than 3s and 1t by 4.7 and 2.2 kcal/mol, respectively.

Investigating bonding in small silicon–carbon clusters: Exploration of the potential energy surfaces of Si3C4, Si4C3, and Si4C4 using ab initio molecular dynamics

A theoretical investigation of the properties of the Si3C4, Si4C3, and Si4C4 clusters is reported. Systematic explorations of the potential energy surfaces of the three clusters are performed using a combination of ab initio molecular dynamics and local energy minimizations using density functional theory. A large number of isomers with a large variety of geometries has been found. The geometries, energies, and vibrational frequencies yielded are discussed. Furthermore, a quantitative analysis of the interatomic distances, angles, and coordination numbers observed, as well as the conclusions on the bonding properties, are presented. The cluster properties are then compared to those of solid SiC and of the smaller Si-C clusters (with size up to 6) obtained in a previous study. Analysis of our results and comparison with bulk properties show that even clusters as small as Si3C4, Si4C3, and Si4C4 exhibit properties similar to those of the amorphous bulk, in particular as for the structures and bonds formed by C atoms.

Linear Triplet SiC3 and SiC5: Results of Coupled Cluster Calculations

Partially restricted coupled cluster calculations involving connected triple substitutions have been carried out for linear SiC