The quest for the carbene bent-pentadiynylidene isomer of C5H2

Electron ionization induced fragmentation pathways of trichloroanisole

Trichloroanisole (TCA) is one of the most significant contaminants in cork stoppers. The presence of TCA leads to an unpleasant odor known as "cork taint", resulting in high economic losses for the cork and wine industries. Hence, the detection, quantification, and characterization of TCA are essential to address this concern. The present study investigates the electron-driven fragmentation pathways of TCA through electron ionization mass spectrometry as a function of electron energy (0-100 eV), and the results are supported by theoretical characterization of ionization potentials, dissociation thresholds, and electron ionization cross sections. The appearance energies of ten cations were measured, including the first experimental evaluation of the molecule's ionization energy at 8.8 ± 0.3 eV, in excellent agreement with the calculations (8.83 eV). For lower energies, around 20 eV, the parent cation accounted for more than 60% of the total ion signal, followed by its demethylated fragment. Taken together, these ion signals could be used as fingerprints of TCA in industrial quality control by low-energy electron ionization mass spectrometry. Fifty other fragments have been identified at higher electron energies, revealing the very rich fragmentation pattern of TCA.

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.

Why Are MgC3H Isomers Missing in the Interstellar Medium?

Considering the recent findings of linear doublet (²Σ⁺) MgC_nH isomers (n = 2, 4, and 6) in the evolved carbon star IRC+10216, various structural isomers of MgC₃H and MgC₃H⁺ are theoretically investigated here. For MgC₃H, 11 doublet and 8 quartet stationary points ranging from 0.0 to 71.8 and 0.0 to 110.1 kcal mol^-1, respectively, have been identified initially at the UωB97XD/6-311++G(2d,2p) level. To get accurate relative energies, further energy evaluations are carried out for all isomers with coupled cluster methods and thermochemical modules such as G3//B3LYP, G4MP2, and CBS-QB3 methods. Unlike the even series, where the global minima are linear molecules with a Mg atom at one end, in the case of MgC₃H, the global minimum geometry turns out to be a cyclic isomer, 2-magnesabicyclo[1.1.0]but-1,3,4-triyl (1, C_2v, ²A₁). In addition, five low-lying isomers, magnesium-substituted cyclopropenylidene (2, C_s, ²A'), 1-magnesabut-2,3-dien-1-yl-4-ylidene (3, C_s, ²A″), 1-magnesabut-2-yn-1-yl-4-ylidene (4, C_s, ²A″), 2λ³-magnesabicyclo[1.1.0]but-1,3-diyl-4-ylidene (5, C_2v;, ²A₁), and 1-magnesabut-2,3-dien-2-yl-4-ylidene (6, C_∞v, ²Σ⁺), were also identified. The doublet linear isomer of MgC₃H, 1-magnesabutatrienyl (10, C_∞v, ²Σ⁺) turns out to be a minimum but lies 54.1 kcal mol^-1 above 1 at the ROCCSD(T)/cc-pVTZ level. The quartet (⁴Σ⁺) electronic state of 10 was also found to be a minimum, but it lies 8.0 kcal mol^-1 above 1 at the same level. Among quartets, isomer 10 is the most stable molecule. The next quartet electronic state (of isomer 11) is 34.4 kcal mol^-1 above 10, and all other quartet electronic states of other isomers are not energetically close to low-lying doublet isomers 2 to 6. Overall, the chemical space of MgC₃H contains more cyclic isomers (1, 2, and 3) on the low-energy side unlike their even-numbered MgC_nH counterparts (n = 2, 4, and 6). Though the quartet electronic state of 10 is linear, it is not the global minimum geometry on the MgC₃H potential energy surface. Isomerization pathways among the low-lying isomers (doublets of 1-4 and a quartet of 10) reveal that these molecules are kinetically stable. For the cation, MgC₃H⁺, the cyclic isomers (1⁺, 2⁺, and 3⁺) are on the low-energy side. The singlet linear isomer, 10⁺, is a fourth-order saddle point. The low-lying cations are quite polar, with dipole moment values of >7.00 D. The current theoretical data would be helpful to both laboratory astrophysicists and radioastronomers for further studies on the MgC₃H^0/+ isomers.

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.

Six Low-Lying Isomers of C11H8 Are Unidentified in the Laboratory-A Theoretical Study

Isomers of C₁₁H₈ have been theoretically examined using density functional theory and coupled-cluster methods. The current investigation reveals that 2aH-cyclopenta[cd]indene (2), 7-ethynyl-1H-indene (6), 4-ethynyl-1H-indene (7), 6-ethynyl-1H-indene (8), 5-ethynyl-1H-indene (9), and 7bH-cyclopenta[cd]indene (10) remain elusive till date in the laboratory. The puckered low-lying isomer 2 lies at 9 kJ mol^-1 below the experimentally known molecule, cyclobuta[de]naphthalene (3), at the fc-CCSD(T)/cc-pVTZ//fc-CCSD(T)/cc-pVDZ level of theory. 2 lies at 36 kJ mol^-1 above the thermodynamically most stable and experimentally known isomer, 1H-cyclopenta[cd]indene (1), at the same level. It is identified that 1,2-H transfer from 1 yields 2H-cyclopenta[cd]indene (14) and subsequent 1,2-H shift from 14 yields 2. Appropriate transition states have been identified, and intrinsic reaction coordinate calculations have been carried out at the B3LYP/6-311+G(d,p) level of theory. Recently, 1-ethynyl-1H-indene (11) has been detected using synchrotron-based vacuum ultraviolet ionization mass spectrometry. 2-Ethynyl-1H-indene (4) and 3-ethynyl-1H-indene (5) have been synthetically characterized in the past. While the derivatives of 7bH-cyclopenta[cd]indene (10) have been isolated elsewhere, the parent compound remains unidentified till date in the laboratory. Although C₁₁H₈ is a key elemental composition of astronomical interest for the formation of polycyclic aromatic hydrocarbons in the interstellar medium, none of its low-lying isomers have been characterized by rotational spectroscopy though they are having a permanent dipole moment (μ ≠ 0). Therefore, energetic and spectroscopic properties have been computed, and the present investigation necessitates new synthetic studies on C₁₁H₈, in particular 2, 6-10, and also rotational spectroscopic studies on all low-lying isomers.

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.

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.

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.

Pentadiynylidene and Its Methyl-Substituted Derivates: Threshold Photoelectron Spectroscopy of R1-C5-R2 Triplet Carbon Chains

Mass-selective threshold photoelectron spectroscopy in the gas phase was employed to characterize the dialkynyl triplet carbenes pentadiynylidene (HC₅H), methylpentadiynylidene (MeC₅H), and dimethylpentadiynylidene (MeC₅Me). Diazo compounds were employed as precursors to generate the carbenes by flash pyrolysis. The R₁-C₅-R₂ carbon chains were photoionized by vacuum ultraviolet (VUV) synchrotron radiation in photoelectron photoion coincidence (PEPICO) experiments. High-level ab initio computations were carried out to support the interpretation of the experiments. For the unsubstituted pentadiynylidene (R₁ = R₂ = H) the recorded spectrum yields an adiabatic ionization energy (IE_ad) of 8.36 ± 0.03 eV. In addition, a second carbene isomer, 3-(didehydrovinylidene)cyclopropene, with a singlet electronic ground state, was identified in the spectrum based on the IE_ad of 8.60 ± 0.03 eV and Franck-Condon simulations. We found that multireference computations are required to reliably calculate the IE_ad for this molecule. CASPT2 computations predicted an IE_ad = 8.55 eV, while coupled-cluster computations significantly overestimate the IE. The cyclic isomer is most likely formed from another isomer of the precursor present in the sample. Stepwise methyl-substitution of the carbene leads to a reduction of the IE to 7.77 ± 0.04 eV for methylpentadiynylidene and 7.27 ± 0.06 eV for dimethylpentadiynylidene. The photoionization and dissociative photoionization of the precursors is investigated as well.