skip to main content

Title: Overcoming the out-of-plane bending issue in an aromatic hydrocarbon: the anharmonic vibrational frequencies of c-(CH)C 3 H 2 +
The challenges associated with the out-of-plane bending problem in multiply-bonded hydrocarbon molecules can be mitigated in quartic force field analyses by varying the step size in the out-of-plane coordinates. Carbon is a highly prevalent element in astronomical and terrestrial environments, but this major piece of its spectra has eluded theoretical examinations for decades. Earlier explanations for this problem focused on method and basis set issues, while this work seeks to corroborate the recent diagnosis as a numerical instability problem related to the generation of the potential energy surface. Explicit anharmonic frequencies for c-(CH)C 3 H 2 + are computed using a quartic force field and the CCSD(T)-F12b method with cc-pVDZ-F12, cc-pVTZ-F12, and aug-cc-pVTZ basis sets. The first of these is shown to offer accuracy comparable to that of the latter two with a substantial reduction in computational time. Additionally, c-(CH)C 3 H 2 + is shown to have two fundamental frequencies at the onset of the interstellar unidentified infrared bands, at 5.134 and 6.088 μm or 1947.9 and 1642.6 cm −1 , respectively. This suggests that the results in the present study should assist in the attribution of parts of these aromatic bands, as well as provide data in support more » of the laboratory or astronomical detection of c-(CH)C 3 H 2 + . « less
; ; ;
Award ID(s):
Publication Date:
Journal Name:
Physical Chemistry Chemical Physics
Page Range or eLocation-ID:
12951 to 12958
Sponsoring Org:
National Science Foundation
More Like this
  1. Context. Formaldehyde is a potential biogenic precursor involved in prebiotic chemical evolution. The cold conditions of the interstellar medium (ISM) allow H 2 CO to be reactive, playing a significant role as a chemical intermediate in formation pathways leading to interstellar complex organic molecules. However, gas-phase molecular formation mechanisms in cold regions of the ISM are poorly understood. Aims. We computationally determine the most favored gas-phase molecular formation mechanisms at local thermodynamic equilibrium conditions that can produce the detected amounts of H 2 CO in diffuse molecular clouds (DMCs), in dark, cold, and dense molecular clouds (DCDMCs), and in three regions of circumstellar envelopes of low-mass protostars (CELMPs). Methods. The potential energy surfaces, thermodynamic functions, and single-point energies for transition states were calculated at the CCSD(T)-F12/cc-pVTZ-F12 and MP2/aug-cc-pVDZ levels of theory and basis sets. Molecular thermodynamics and related partition functions were obtained by applying the Maxwell-Boltzmann quantum statistics theory from energies computed at CCSD(T)-F12/cc-pVTZ-F12 with corrections for zero-point energy. A literature review on detected abundances of reactants helped us to propose the most favorable formation routes. Results. The most probable reactions that produce H 2 CO in cold astrophysical regions are: 1 CH 2 + ⋅ 3 O 2 →more »1 H 2 CO + O⋅( 3 P) in DMCs, ⋅ 3 CH 2 + ⋅ 3 O 2 → 1 H 2 CO + ⋅O( 3 P) in DCDMCs, and ⋅CH 3 + ⋅O( 3 P) → 1 H 2 CO + ⋅H in region III, ⋅CH 3 +⋅O( 1 D) → 1 H 2 CO + ⋅H in region II, and 1 CH 2 + ⋅ 3 O 2 → 1 H 2 CO + ⋅O( 3 P) in region I belonging to CELMPs. Conclusions. Quantum chemical calculations suggest that the principal carbonaceous precursors of H 2 CO in cold regions for the gas-phase are CH 2 (a 1 A 1 ), and ⋅CH 2 (X 3 B 1 ) combined with ⋅O 2 ( 3 Σ g ) and ⋅CH 3 ( 2 A ” ) + ⋅O( 3 P) / O( 1 D). Reactions based on more complex reagents yield less effective thermodynamics in the gas-phase H 2 CO molecular formation.« less
  2. The gas-phase reaction of the methylidyne (CH; X 2 Π) radical with dimethylacetylene (CH 3 CCCH 3 ; X 1 A 1g ) was studied at a collision energy of 20.6 kJ mol −1 under single collision conditions with experimental results merged with ab initio calculations of the potential energy surface (PES) and ab initio molecule dynamics (AIMD) simulations. The crossed molecular beam experiment reveals that the reaction proceeds barrierless via indirect scattering dynamics through long-lived C 5 H 7 reaction intermediate(s) ultimately dissociating to C 5 H 6 isomers along with atomic hydrogen with atomic hydrogen predominantly released from the methyl groups as verified by replacing the methylidyne with the D1-methylidyne reactant. AIMD simulations reveal that the reaction dynamics are statistical leading predominantly to p28 (1-methyl-3-methylenecyclopropene, 13%) and p8 (1-penten-3-yne, 81%) plus atomic hydrogen with a significant amount of available energy being channeled into the internal excitation of the polyatomic reaction products. The dynamics are controlled by addition to the carbon–carbon triple bond with the reaction intermediates eventually eliminating a hydrogen atom from the methyl groups of the dimethylacetylene reactant forming 1-methyl-3-methylenecyclopropene (p28). The dominating pathways reveal an unexpected insertion of methylidyne into one of the six carbon–hydrogen singlemore »bonds of the methyl groups of dimethylacetylene leading to the acyclic intermediate, which then decomposes to 1-penten-3-yne (p8). Therefore, the methyl groups of dimethylacetylene effectively ‘screen’ the carbon–carbon triple bond from being attacked by addition thus directing the dynamics to an insertion process as seen exclusively in the reaction of methylidyne with ethane (C 2 H 6 ) forming propylene (CH 3 C 2 H 3 ). Therefore, driven by the screening of the triple bond, one propynyl moiety (CH 3 CC) acts in four out of five trajectories as a spectator thus driving an unexpected, but dominating chemistry in analogy to the methylidyne – ethane system.« less
  3. Abstract

    Five substituted cyclopropenylidene derivatives (c-C3HX, X=CN, OH, F, NH2), all currently undetected in the interstellar medium (ISM), are found herein to have mechanistically viable, gas-phase formation pathways through neutral–neutral additions of ·X ontoc-C3H2. The detection and predicted formation mechanism ofc-C3HC2H introduces a need for the chemistry ofc-C3H2and any possible derivatives to be more fully explored. Chemically accurate CCSD(T)-F12/cc-pVTZ-F12 calculations provide exothermicities of additions of various radical species toc-C3H2, alongside energies of submerged intermediates that are crossed to result in product formation. Of the novel reaction mechanisms proposed, the addition of the cyano radical is the most exothermic at -16.10 kcal mol−1. All five products are found to or are expected to have at least one means of associating barrierlessly to form a submerged intermediate, a requirement for the cold chemistry of the ISM. The energetically allowed additions arise as a result of the strong electrophilicity of the radical species as well as the product stability gained through substituent-ring conjugation.

  4. We report and interpret recently recorded high-resolution infrared spectra for the fundamentals of the CH 2 scissors and CH stretches of gas phase cyclopentane at −26.1 and −50 °C, respectively. We extend previous theoretical studies of this molecule, which is known to undergo barrierless pseudorotation due to ring puckering, by constructing local mode Hamiltonians of the stretching and scissor vibrations for which the frequencies, couplings, and linear dipoles are calculated as functions of the pseudorotation angle using B3LYP/6-311++(d,p) and MP2/cc-pVTZ levels of theory. Symmetrization ( D 5 h ) of the vibrational basis sets leads to simple vibration/pseudorotation Hamiltonians whose solutions lead to good agreement with the experiment at medium resolution, but which miss interesting line fractionation when compared to the high-resolution spectra. In contrast to the scissor motion, pseudorotation leads to significant state mixing of the CH stretches, which themselves are Fermi coupled to the scissor overtones.
  5. The kinetic energy dependent reactions of Re + with SO 2 were studied with guided ion beam tandem mass spectrometry. ReO + , ReO 2 + , and OReS + species were observed as products, all in endothermic reactions. Modeling of the kinetic energy dependent cross sections yields 0 K bond dissociation energies (BDEs, in eV) of 4.78 ± 0.06 (Re + –O), 5.75 ± 0.02 (Re + –O 2 ), and 4.35 ± 0.14 (Re + –SO). The latter two values can be combined with other information to derive the additional values 6.05 ± 0.05 (ORe + –O) and 4.89 ± 0.19 (ORe + –S). BDEs of ReO + and ReO 2 + agree with literature values whereas the values for OReS + are the first measurements. The former result is obtained even though formation of ground state ReO + is spin-forbidden. Quantum mechanical calculations at the B3LYP level of theory with a def2-TZVPPD basis set yield results that agree reasonably well with experimental values. Additional calculations at the BP86 and CCSD(T) levels of theory using def2-QZVPPD and aug-cc-pVxZ (x = T, Q, and 5) basis sets were performed to compare thermochemistry with experiment to determine that ReO 2more »+ rather than the isobaric ReS + is formed. Product ground states are 3 Δ 3 (ReO + ), 3 B 1 (OReO + ), 5 Π −1 (ReS + ), and 3 A′′(OReS + ) after including empirical spin–orbit corrections, which means that formation of ground state products is spin-forbidden for all three product channels. The potential energy surfaces for the ReSO 2 + system were also explored at the B3LYP/def2-TZVPPD level and exhibited no barriers in excess of the endothermicities for all products. BDEs for rhenium oxide and sulfide diatomics and triatomics are compared and discussed. The present result for formation of ReO + is compared to that for formation of ReO + in the reactions of Re + + O 2 and CO, where the former system exhibited interesting dual cross section features. Results are consistent with the hypothesis that the distinction of in-plane and out-of-plane C S symmetry in the triatomic systems might be the explanation for the two endothermic features observed in the Re + + O 2 reaction.« less