The Reaction of Tricarbon with Acetylene: An Ab Initio/RRKM Study of the Potential Energy Surface and Product Branching Ratios

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.

New Carbenes and Cyclic Allenes Energetically Comparable to Experimentally Known 1-Azulenylcarbene

1-Azulenylcarbene (18; 0 kJ mol^-1) is experimentally known as the key reactive intermediate for the rearrangement reactions of aryl carbenes in the laboratory. Here, using coupled-cluster methods up to the fc-CCSD(T)/cc-pVTZ//fc-CCSD(T)/cc-pVDZ level, thirteen new carbenes and one new cyclic allene are theoretically identified within the C₁₁H₈ elemental composition that either energetically lie below or very close to 18. While the cyclic allene, bicyclo[5.4.0]undeca-2,3,5,7,9,11-hexene (1; -166 kJ mol^-1), is the experimentally known lowest energy isomer, three other cyclic allenes, bicyclo[5.4.0]undeca-1,2,4,6,8,10-hexene (2; -100 kJ mol^-1), bicyclo[5.4.0]undeca-1,3,4,6,8,10-hexene (3; -97 kJ mol^-1), and bicyclo[6.3.0]undeca-1,2,4,6,8,10-hexene (13; -42 kJ mol^-1), demand new experimental studies. In total, thirty-one isomers are studied in this work (within -166 to +15 kJ mol^-1 from 18) and all are found to be polar (μ ≠ 0). Among these, 1H-benzo[7]annulen-1-ylidene (17; -4 kJ mol^-1; μ = 5.24 D), bicyclo[5.4.0]undeca-2,4,6,8,11-pentaene-10-ylidene (24; 13 kJ mol^-1; μ = 7.59 D), 5-methylene-naphthalen-1-ylidene (26; 15 kJ mol^-1; μ = 5.32 D), 6-methylene-naphthalen-2-ylidene (27; -43 kJ mol^-1; μ = 6.60 D), and 8-methylene-naphthalen-2-ylidene (28; -39 kJ mol^-1; μ = 5.55 D) are competitively polar compared to 18 (μ = 5.39 D). Therefore, these carbene molecules are potential targets for rotational spectroscopists and radioastronomers. Considering the importance of naphthyl and azulenylcarbenes in reactive intermediate chemistry, mechanisms of different rearrangement reactions and plausible formation pathways of some of these new carbenes are studied in this work using density functional theory.

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.

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.

Formation mechanism and spectroscopy of C6H radicals in extreme environments: a theoretical study

This study examined the reaction mechanisms of singlet (rhombic) and triplet (linear) C4 with acetylene by using accurate ab initio CCSD(T)/CBS//B3LYP/6-311G(d,p) calculations followed by a kinetic analysis of various reaction pathways and computations of relative product yields in combustion and planetary atmospheres. These calculations were combined with the Rice-Ramsperger-Kassel-Marcus (RRKM) calculations of reaction rate constants for predicting product-branching ratios, which depend on the collision energy under single-collision conditions. The results demonstrate that the initial reaction begins with the formation of an intermediate 3i2 with an entrance barrier of 3.0 kcal mol-1 and an intermediate 1i1 without entrance barriers. The product-branching ratios obtained by solving kinetic equations with individual rate constants calculated using the RRKM and variational transition-state theories for determining the collision energies between 5 kcal mol-1 and 25 kcal mol-1 demonstrate that l-C6H + H is the dominant reaction product, whereas HC3C3 + H, l-C6 + H2, c-C6H + H, and c-C6 + H2 are minor products. The electronic absorption spectra of solid neon matrices in the range of 17 140-22 200 cm-1 were obtained by Maier et al., and the optimized ground and excited state structures of C6H were used to simulate the absorption spectra by one-photon excitation equations. The displaced harmonic oscillator approximation and the Franck-Condon approximation were used to simulate the absorption spectrum of the B2Π ← X2Π transition of C6H. This indicates that the vibronic structures were dominated by one of the six active completely symmetric modes, with v3 being the most crucial.

Theoretical Studies of Two Key Low-Lying Carbenes of C5H2 Missing in the Laboratory

The equilibrium geometries and spectroscopic properties of two key singlet carbenes, buta-1,3-diynylcarbene (6) and 2-methylenebicyclo[1.1.0]but-1(3)-en-4-ylidene (9), which have not been experimentally observed to date, are investigated using high-level coupled-cluster (CC) methods. The current theoretical study necessitates new experimental data on C₅H₂ isomers considering the relevance of these molecules to interstellar chemistry. Bent-pentadiynylidene (4) has been missing in the laboratory and the prime focus of our earlier theoretical work. The present theoretical study indicates that isomers 6 and 9 are also viable experimental targets. Apart from ethynylcyclopropenylidene (2), pentatetraenylidene (3), ethynylpropadienylidene (5), and 3-(didehydrovinylidene)cyclopropene (8), which are identified by Fourier transform microwave spectroscopy, the dipole moments of elusive 4, 6, and 9 are also nonzero (μ ≠ 0). The relative energies of these isomers, calculated at the CCSDT(Q)/CBS level of theory, with respect to linear triplet pentadiynylidene (1) reveal that they all lie within 25.1 kcal mol^-1. Therefore, geometric, energetic, aromatic, and spectroscopic parameters are reported here, which may assist the efforts of molecular spectroscopists in the future. Anharmonic vibrational calculations on isomers 6 and 9 indicate that the former is loosely bound and would be challenging to be detected experimentally. Among the undetected carbenes, 9 may be considered as a potential target molecule considering its higher polarity and aromatic nature.

Dynamics of the reaction of C₃(a³Πu) radicals with C₂H₂: a new source for the formation of C₅H

The reaction C3(a(3)Πu) + C2H2 → C5H + H was investigated at collision energy 10.9 kcal mol(-1) that is less than the enthalpy of ground-state reaction C3(X(1)Σg (+)) + C2H2 → C5H + H. C3(a(3)Πu) radicals were synthesized from 1% C4F6/He by pulsed high-voltage discharge. The title reaction was conducted in a crossed molecular-beam apparatus equipped with a quadrupole-mass filter. Product C5H was interrogated with time-of-flight spectroscopy and synchrotron vacuum-ultraviolet ionization. Reactant C3(a(3)Πu) and product C5H were identified using photoionization spectroscopy. The ionization thresholds of C3(X(1)Σg(+)) and C3(a(3)Πu) are determined as 11.6 ± 0.2 eV and 10.0 ± 0.2 eV, respectively. The C5H product is identified as linear pentynylidyne that has an ionization energy 8.4 ± 0.2 eV. The title reaction releases translational energy 10.6 kcal mol(-1) in average and has an isotropic product angular distribution. The quantum-chemical calculation indicates that the C3(a(3)Πu) radical attacks one of the carbon atoms of C2H2 and subsequently a hydrogen atom is ejected to form C5H + H, in good agreement with the experimental observation. As far as we are aware, the C3(a(3)Πu) + C2H2 reaction is investigated for the first time. This work gives an implication for the formation of C5H from the C3(a(3)Πu) + C2H2 reaction occurring in a combustion or discharge process of C2H2.

A crossed molecular beam study on the synthesis of the interstellar 2,4-pentadiynylidyne radical (HCCCCC)

Crossed molecular beam experiments are performed to elucidate the synthesis of the 2,4-penta-diynylidyne [HCCCCC(X (2)Pi)] radical under single collision conditions--a crucial reaction intermediate to form polycyclic aromatic hydrocarbons and carbonaceous nanostructures in the interstellar medium and in combustion flames. The experiments demonstrate that the chemical dynamics of ground state carbon reacting with diacetylene [HCCCCH(X (1)Sigma(g)(+))] are indirect and proceed via addition of the electrophilic carbon atom to the pi electron density of the diacetylene molecule yielding ultimately the carbenelike HCCCCCH(X (3)Sigma(g)(-)) molecule. This intermediate fragments via hydrogen atom emission to yield the 2,4-pentadiynylidyne [HCCCCC(X (2)Pi)] radical. The chemical dynamics elucidated also allows us to predict that reaction of carbon atoms with polyynes of the generic formula H(C[triple bond]C)(n)H leads to the formation of hydrogen-terminated carbon clusters of the generic form HC(2n+1) in extreme environments. The acetylene-related reactivity and electronic structure of the diacetylene molecule also allow us to project that reactions of the diacetylene molecule with cyano and ethynyl radicals result in a stepwise extension of the carbon skeleton forming cyanodiacetylene (HCCCCCN) and triacetylene (HCCCCCCH) plus atomic hydrogen. These predictions open the door to extensive laboratory studies involving hitherto poorly understood reactions of the diacetylene molecule under single collision conditions.