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 of the laboratory or astronomical detection of c-(CH)C 3 H 2 + .  more » « less
Award ID(s):
1757220
NSF-PAR ID:
10251876
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Physical Chemistry Chemical Physics
Volume:
22
Issue:
23
ISSN:
1463-9076
Page Range / eLocation ID:
12951 to 12958
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Rovibrational spectral data for several tetra-atomic silicon carbide clusters (TASCCs) are computed in this work using a CCSD(T)-F12b/cc-pCVTZ-F12 quartic force field. Accurate theoretical spectroscopic data may facilitate the observation of TASCCs in the interstellar medium which may lead to a more complete understanding of how the smallest silicon carbide (SiC) solids are formed. Such processes are essential for understanding SiC dust grain formation. Due to SiC dust prevalence in the interstellar medium, this may also shed light on subsequent planetary formation. Rhomboidal Si2C2is shown here to have a notably intense (247 km mol−1) anharmonic vibrational frequency at 988.1 cm−1(10.1 μm) forν2, falling into one of the spectral emission features typically associated with unknown infrared bands of various astronomical regions. Notable intensities are also present for several of the computed anharmonic vibrational frequencies including the cyclic forms of C4, SiC3, Si3C, and Si4. These features in the 6–10 μm range are natural targets for infrared observation with theJames Webb Space Telescope(JWST)’s MIRI instrument. Additionally,t-Si2C2,d-Si3C, andr-SiC3each possess dipole moments of greater than 2.0 D making them interesting targets for radioastronomical searches especially sinced-SiC3is already known in astrophysical media.

     
    more » « less
  2. 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 → 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. 
    more » « less
  3. High-level rovibrational characterization of methanediol, the simplest geminal diol, using state-of-the-art, purely ab initio techniques unequivocally confirms previously reported gas phase preparation of this simplest geminal diol in its C 2 conformation. The F12-TZ-cCR and F12-DZ-cCR quartic force fields (QFFs) utilized in this work are among the largest coupled cluster-based anharmonic frequencies computed to date, and they match the experimental band origins of the spectral features in the 980–1100 cm −1 range to within 3 cm −1 , representing a significant improvement over previous studies. The simulated spectrum also matches the experimental spectrum in the strong Q branch feature and qualitative shape of the 980–1100 cm −1 region. Additionally, the full set of rotational constants, anharmonic vibrational frequencies, and quartic and sextic distortion constants are provided for both the lowest energy C 2 conformer as well as the slightly higher C s conformer. Several vibrational modes have intensities of 60 km mol −1 or higher, facilitating potential astronomical or atmospheric detection of methanediol or further identification in laboratory work especially now that gas phase synthesis of this molecule has been established. 
    more » « less
  4. 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 single 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. 
    more » « less
  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 2 + 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. 
    more » « less