Experimental and theoretical study of the low-temperature kinetics of the reaction of CN with CH2O and implications for interstellar environments.

Rate coefficients for the reaction of CN with CH2O were measured for the first time below room temperature in the range 32-103 K using a pulsed Laval nozzle apparatus together with the Pulsed Laser Photolysis-Laser-Induced Fluorescence technique. The rate coefficients exhibited a strong negative temperature dependence, reaching (4.62 ± 0.84) × 10-11 cm3 molecule-1 s-1 at 32 K, and no pressure dependence was observed at 70 K. The potential energy surface (PES) of the CN + CH2O reaction was calculated at the CCSD(T)/aug-cc-pVTZ//M06-2X/aug-cc-pVTZ level of theory, with the lowest energy channel to reaction characterized by the formation of a weakly-bound van der Waals complex, bound by 13.3 kJ mol-1, prior to two transition states with energies of -0.62 and 3.97 kJ mol-1, leading to the products HCN + HCO or HNC + HCO, respectively. For the formation of formyl cyanide, HCOCN, a large activation barrier of 32.9 kJ mol-1 was calculated. Reaction rate theory calculations were performed with the MESMER (Master Equation Solver for Multi Energy well Reactions) package on this PES to calculate rate coefficients. While this ab initio description provided good agreement with the low-temperature rate coefficients, it was not capable of describing the high-temperature experimental rate coefficients from the literature. However, increasing the energies and imaginary frequencies of both transition states allowed MESMER simulations of the rate coefficients to be in good agreement with data spanning 32-769 K. The mechanism for the reaction is the formation of a weakly-bound complex followed by quantum mechanical tunnelling through the small barrier to form HCN + HCO products. MESMER calculations showed that channel generating HNC is not important. MESMER simulated the rate coefficients from 4-1000 K which were used to recommend best-fit modified Arrhenius expressions for use in astrochemical modelling. The UMIST Rate12 (UDfa) model yielded no significant changes in the abundances of HCN, HNC, and HCO for a variety of environments upon inclusion of rate coefficients reported here. The main implication from this study is that the title reaction is not a primary formation route to the interstellar molecule formyl cyanide, HCOCN, as currently implemented in the KIDA astrochemical model.

Comparison of intermolecular energy transfer from vibrationally excited benzene in mixed nitrogen-benzene baths at 140 K and 300 K.

Gas phase intermolecular energy transfer (IET) is a fundamental component of accurately explaining the behavior of gas phase systems in which the internal energy of particular modes of molecules is greatly out of equilibrium. In this work, chemical dynamics simulations of mixed benzene/N2 baths with one highly vibrationally excited benzene molecule (Bz*) are compared to experimental results at 140 K. Two mixed bath models are considered. In one, the bath consists of 190 N2 and 10 Bz, whereas in the other bath, 396 N2 and 4 Bz are utilized. The results are compared to results from 300 K simulations and experiments, revealing that Bz*-Bz vibration-vibration IET efficiency increased at low temperatures consistent with longer lived "chattering" collisions at lower temperatures. In the simulations, at the Bz* excitation energy of 150 kcal/mol, the averaged energy transferred per collision, ⟨ΔEc⟩, for Bz*-Bz collisions is found to be ∼2.4 times larger in 140 K than in 300 K bath, whereas this value is ∼1.3 times lower for Bz*-N2 collisions. The overall ⟨ΔEc⟩, for all collisions, is found to be almost two times larger at 140 K compared to the one obtained from the 300 K bath. Such an enhancement of IET efficiency at 140 K is qualitatively consistent with the experimental observation. However, the possible reasons for not attaining a quantitative agreement are discussed. These results imply that the bath temperature and molecular composition as well as the magnitude of vibrational energy of a highly vibrationally excited molecule can shift the overall timescale of rethermalization.

Non-statistical intermolecular energy transfer from vibrationally excited benzene in a mixed nitrogen-benzene bath.

A chemical dynamics simulation was performed to model experiments [N. A. West et al., J. Chem. Phys. 145, 014308 (2016)] in which benzene molecules are vibrationally excited to 148.1 kcal/mol within a N2-benzene bath. A significant fraction of the benzene molecules are excited, resulting in heating of the bath, which is accurately represented by the simulation. The interesting finding from the simulations is the non-statistical collisional energy transfer from the vibrationally excited benzene C6H6 * molecules to the bath. The simulations find that at ∼10-7 s and 1 atm pressure there are four different final temperatures for C6H6 * and the bath. N2 vibration is not excited and remains at the original bath temperature of 300 K. Rotation and translation degrees of freedom of both N2 and C6H6 in the bath are excited to a final temperature of ∼340 K. Energy transfer from the excited C6H6 * molecules is more efficient to vibration of the C6H6 bath than its rotation and translation degrees of freedom, and the final vibrational temperature of the C6H6 bath is ∼453 K, if the average energy of each C6H6 vibration mode is assumed to be RT. There is no vibrational equilibration between C6H6 * and the C6H6 bath molecules. When the simulations are terminated, the vibrational temperatures of the C6H6 * and C6H6 bath molecules are ∼537 K and ∼453 K, respectively. An important question is the time scale for complete energy equilibration of the C6H6 * and N2 and C6H6 bath system. At 1 atm and 300 K, the experimental V-T (vibration-translation) relaxation time for N2 is ∼10-4 s. The simulation time was too short for equilibrium to be attained, and the time for complete equilibration of C6H6 * vibration with translation, rotation, and vibration of the bath was not determined.

Resolving the energy and temperature dependence of C6H6 (∗) collisional relaxation via time-dependent bath temperature measurements.

The relaxation of highly vibrationally excited benzene, generated by 193 nm laser excitation, was studied using the transient rotational-translational temperature rise of the N2 bath, which was measured by proxy using two-line laser induced fluorescence of seeded NO. The resulting experimentally measured time-dependent N2 temperature rises were modeled with MultiWell based simulations of Collisional Energy Transfer (CET) from benzene vibration to N2 rotation-translation. We find that the average energy transferred in benzene deactivating collisions depends linearly on the internal energy of the excited benzene molecules and depends approximately linearly on the N2 bath temperature between 300 K and 600 K. The results are consistent with experimental studies and classical trajectory calculations of CET in similar systems.

A method for the determination of speed-dependent semi-classical vector correlations from sliced image anisotropies.

We present analytical expressions relating the bipolar moment ß(Q)(K)(k(1)k(2)) parameters of Dixon to the measured anisotropy parameters of different pump/probe geometry sliced ion images. In the semi-classical limit, when there is no significant coherent contribution from multiple excited states to fragment angular momentum polarization, the anisotropy of the images alone is sufficient to extract the ß(Q)(K)(k(1)k(2)) parameters with no need to reference relative image intensities. The analysis of sliced images is advantageous since the anisotropy can be directly obtained from the image at any radius without the need for 3D-deconvolution, which is not applicable for most pump/probe geometries. This method is therefore ideally suited for systems which result in a broad distribution of fragment velocities. The bipolar moment parameters are obtained for NO(2) dissociation at 355 nm using these equations, and are compared to the bipolar moment parameters obtained from a proven iterative fitting technique for crushed ion images. Additionally, the utility of these equations in extracting speed-dependent bipolar moments is demonstrated on the recently investigated NO(3) system.