The proton transfer reaction H3 + + CO → HCO+/HOC+ + H2 has gained considerable attention in the literature due to its importance in interstellar chemistry. The reaction productsformyl cation (HCO+) and isoformyl cation (HOC+)are known to initiate multiple chemical reaction networks, resulting in complex molecules found in space. Several experimental and theoretical studies probing the structure and energetics of the [H3CO]+ system, HCO+/HOC+ product branching ratios, reaction mechanisms, etc., have been reported in the literature. In the present work, we investigated the H3 + + CO bimolecular reaction in the gas phase using direct dynamics methodology. The simulation conditions were chosen to mimic recently reported velocity map imaging experiments on the same reaction. The calculations were performed using the density functional PBE0/aug-cc-pVDZ level of electronic structure theory. Internal energy and scattering angle distributions of reaction products found from the simulations are in qualitative agreement with the experiment. However, the product branching ratios at low collision energies were in contrast with the experimental predictions. Interesting dynamical features were observed in the simulations, and detailed atomic level mechanisms are presented.

γ-Ketohydroperoxide [3-(hydroperoxy)propanal] is an important reagent in synthetic chemistry and, in particular, oxidation reactions. It is considered to be a precursor for secondary organic aerosol formation in the troposphere. Due to enhanced reactivity and limitations associated with analytical techniques, theoretical methods have been employed to study the unimolecular reactivity of hydroperoxides. A number of automated reaction discovery techniques have been used to study the reactivity of γ-ketohydroperoxide, and a large number of reactions have been reported in such studies. In the present work, we have investigated the unimolecular reaction dynamics of this molecule using electronic structure theory calculations and direct chemical dynamics simulations to assess the relevance of different reaction pathways. Classical trajectories were launched from the reactant well with fixed amounts of total energies and integrated on-the-fly using density functional B3LYP/6-31+G* model chemistry. Three dissociation channels among the previously reported reactions were identified as important. Korcek decomposition, which was proposed earlier as a source of carbonyl compounds from thermal decomposition of γ-ketohydroperoxide, was not observed in the present high-temperature simulations. However, trajectories showed the formation of carbonyl compounds such as aldehydes via other pathways. Results are compared with previous studies, and detailed atomic-level reaction mechanisms are presented.

Thiophene is an organo-sulfur aromatic molecule present in fossil fuels and alternate fuels such as shale oils and contributes to air pollution via fuel burning. Hence, it is essential to remove thiophene and its derivatives during the refining process. In this regard, experimental and electronic structure theory studies investigating the thermal decomposition of thiophene have been reported in the literature. In the present work, high temperature thermal decomposition of thiophene was investigated using Born–Oppenheimer direct dynamics simulations. The trajectory integrations were performed on-the-fly at the density functional B3LYP/6-31+G* level of electronic structure theory to investigate the atomic level decomposition mechanisms. Simulation results show that C–S cleavage accompanied by an intramolecular proton transfer to C is the dominant initial dissociation step. Acetylene was observed as primary decomposition product and the results are in agreement with previous experimental studies.