Experimental and theoretical studies of the O(3P) + C2H4reaction dynamics: Collision energy dependence of branching ratios and extent of intersystem crossing

High-Resolution Spectroscopic Investigation of the CH2CHO Radical in the Sub-Millimeter Region

In this work, the pure rotational spectrum of the vinoxy radical (CH2CHO) has been studied at millimeter and sub-millimeter wavelengths (110-860 GHz). CH2CHO was produced by H-abstraction from acetaldehyde (CH3CHO) using atomic fluorine in a double-pass absorption cell at room temperature. A Zeeman-modulation spectrometer, in which an external magnetic field generated inside the absorption cell is amplitude-modulated, was used to record the pure rotational transitions of the radical. The recorded spectra are devoid of signals from closed-shell species, allowing for relatively fast acquisitions over wide spectral windows. Transitions involving values of the rotational quantum numbers N″ and Ka″ up to 41 and 18, respectively, were measured and combined with all available high-resolution literature data (both pure rotation and ground-state combination differences from ro-vibration) to greatly improve the modeling of the CH2CHO spectrum. The combined experimental line list is fit using a semirigid rotor Hamiltonian, and the results are compared to quantum chemical calculations. This laboratory study provides the spectroscopic information needed to search for CH2CHO in various interstellar environments, from cold (e.g., typically 10 K for dense molecular clouds) to warm (e.g., ∼200 K for hot corinos) objects.

A combined crossed molecular beam and theorerical study of the O(3P,1D) + acrylonitrile (CH2CHCN) reactions and implications for combustion and extraterrestrial environments

Acrylonitrile (CH₂CHCN) is ubiquitous in space (molecular clouds, solar-type star forming regions, and circumstellar envelopes) and is also abundant in the upper atmosphere of Titan. The reaction O(³P) + CH₂CHCN can be of relevance in the chemistry of the interstellar medium because of the abundance of atomic oxygen. The oxidation of acrylonitrile is also important in combustion as the thermal decomposition of pyrrolic and pyridinic structures present in fuel-bound nitrogen generates many nitrogen-bearing compounds, including acrylonitrile. Despite its relevance, limited information exists on this reaction. We report a combined experimental and theoretical investigation of the reactions of acrylonitrile with both ground ³P and excited ¹D atomic oxygen. From product angular and time-of-flight distributions in crossed molecular beam experiments with mass spectrometric detection at a collision energy, E_c, of 31.4 kJ mol^-1, we have identified the primary reaction products and determined their branching fractions (BFs). Theoretical calculations of the relevant triplet and singlet potential energy surfaces (PESs) were performed to interpret the experimental results and elucidate the reaction mechanism. Adiabatic statistical calculations of product BFs for the decomposition of the main triplet and singlet intermediates have been carried out. Combining the experimental and theoretical results, we conclude that the O(³P) reaction leads to two main product channels: (i) CH₂CNH (ketenimine) + CO (dominant with a BF of 0.87 ± 0.05), formed via efficient intersystem crossing from the entrance triplet PES to the underlying singlet PES, and (ii) HCOCHCN + H (minor, with a BF of 0.13 ± 0.05), occurring adiabatically on the triplet PES. Our study suggests the inclusion of this reaction as a possible destruction pathway of CH₂CHCN and a possible formation route of CH₂CNH in the interstellar medium. The O(¹D) + CH₂CHCN reaction mainly leads to the formation of CH₂CNH + CO adiabatically on the singlet PES. This result can improve models related to the chemistry of interstellar ice and cometary comas, where O(¹D) reactions can play a role. Overall, our results are expected to be useful for improving the models of combustion and extraterrestrial environments.

A highly accurate full-dimensional ab initio potential surface for the rearrangement of methylhydroxycarbene (H3C-C-OH)

We report here a full-dimensional machine learning global potential surface (PES) for the rearrangement of methylhydroxycarbene (H₃C-C-OH, 1t). The PES is trained with the fundamental invariant neural network (FI-NN) method on 91 564 ab initio energies calculated at the UCCSD(T)-F12a/cc-pVTZ level of theory, covering three possible product channels. FI-NN PES has the correct symmetry properties with respect to permutation of four identical hydrogen atoms and is suitable for dynamics studies of the 1t rearrangement. The averaged root mean square error (RMSE) is 11.4 meV. Six important reaction pathways, as well as the energies and vibrational frequencies at the stationary geometries on these pathways are accurately preproduced by our FI-NN PES. To demonstrate the capacity of the PES, we calculated the rate coefficient of hydrogen migration in -CH₃ (path A) and hydrogen migration of -OH (path B) with instanton theory on this PES. Our calculations predicted the half-life of 1t to be 95 min, which is excellent in agreement with experimental observations.

Reactions O(3P, 1D) + HCCCN(X1Σ+) (Cyanoacetylene): Crossed-Beam and Theoretical Studies and Implications for the Chemistry of Extraterrestrial Environments

Cyanoacetylene (HCCCN), the first member of the cyanopolyyne family (HC_nN, where n = 3, 5, 7, ...), is of particular interest in astrochemistry being ubiquitous in space (molecular clouds, solar-type protostars, protoplanetary disks, circumstellar envelopes, and external galaxies) and also relatively abundant. It is also abundant in the upper atmosphere of Titan and comets. Since oxygen is the third most abundant element in space, after hydrogen and helium, the reaction O + HCCCN can be of relevance in the chemistry of extraterrestrial environments. Despite that, scarce information exists not only on the reactions of oxygen atoms with cyanoacetylene but with nitriles in general. Here, we report on a combined experimental and theoretical investigation of the reactions of cyanoacetylene with both ground ³P and excited ¹D atomic oxygen and provide detailed information on the primary reaction products, their branching fractions (BFs), and the overall reaction mechanisms. More specifically, the reactions of O(³P, ¹D) with HCCCN(X¹Σ⁺) have been investigated under single-collision conditions by the crossed molecular beams scattering method with mass spectrometric detection and time-of-flight analysis at the collision energy, E_c, of 31.1 kJ/mol. From product angular and time-of-flight distributions, we have identified the primary reaction products and determined their branching fractions (BFs). Theoretical calculations of the relevant triplet and singlet potential energy surfaces (PESs) were performed to assist the interpretation of the experimental results and clarify the reaction mechanism. Adiabatic statistical calculations of product BFs for the decomposition of the main triplet and singlet intermediates have also been carried out. Merging together the experimental and theoretical results, we conclude that the O(³P) reaction is characterized by a minor adiabatic channel leading to OCCCN (cyanoketyl) + H (experimental BF = 0.10 ± 0.05), while the dominant channel (BF = 0.90 ± 0.05) occurs via intersystem crossing to the underlying singlet PES and leads to formation of ¹HCCN (cyanomethylene) + CO. The O(¹D) reaction is characterized by the same two channels, with the relative CO/H yield being slightly larger. Considering the recorded reactive signal and the calculated entrance barrier, we estimate that the rate coefficient for reaction O(³P) + HC₃N at 300 K is in the 10^-12 cm³ molec^-1 s^-1 range. Our results are expected to be useful to improve astrochemical and photochemical models. In addition, they are also relevant in combustion chemistry, because the thermal decomposition of pyrrolic and pyridinic structures present in fuel-bound nitrogen generates many nitrogen-bearing compounds, including cyanoacetylene.

Crossed-beam and theoretical studies of multichannel nonadiabatic reactions: branching fractions and role of intersystem crossing for O(3P)+1,3-butadiene

Atomic oxygen reactions can contribute significantly to the oxidation of unsaturated aliphatic and aromatic hydrocarbons. The reaction mechanism is started by electrophilic O atom addition to the unsaturated bond(s) to...

Crossed-Beam and Theoretical Studies of the O(3P, 1D) + Benzene Reactions: Primary Products, Branching Fractions, and Role of Intersystem Crossing

Reliable modeling of hydrocarbon oxidation relies critically on knowledge of the branching fractions (BFs) as a function of temperature (T) and pressure (p) for the products of the reaction of the hydrocarbon with atomic oxygen in its ground state, O(³P). During the past decade, we have performed in-depth investigations of the reactions of O(³P) with a variety of small unsaturated hydrocarbons using the crossed molecular beam (CMB) technique with universal mass spectrometric (MS) detection and time-of-flight (TOF) analysis, combined with synergistic theoretical calculations of the relevant potential energy surfaces (PESs) and statistical computations of product BFs, including intersystem crossing (ISC). This has allowed us to determine the primary products, their BFs, and extent of ISC to ultimately provide theoretical channel-specific rate constants as a function of T and p. In this work, we have extended this approach to the oxidation of one of the most important species involved in the combustion of aromatics: the benzene (C₆H₆) molecule. Despite extensive experimental and theoretical studies on the kinetics and dynamics of the O(³P) + C₆H₆ reaction, the relative importance of the C₆H₅O (phenoxy) + H open-shell products and of the spin-forbidden C₅H₆ (cyclopentadiene) + CO and phenol adduct closed-shell products are still open issues, which have hampered the development of reliable benzene combustion models. With the CMB technique, we have investigated the reaction dynamics of O(³P) + benzene at a collision energy (E_c) of 8.2 kcal/mol, focusing on the occurrence of the phenoxy + H and spin-forbidden C₅H₆ + CO and phenol channels in order to shed further light on the dynamics of this complex and important reaction, including the role of ISC. Concurrently, we have also investigated the reaction dynamics of O(¹D) + benzene at the same E_c. Synergistic high-level electronic structure calculations of the underlying triplet/singlet PESs, including nonadiabatic couplings, have been performed to complement and assist the interpretation of the experimental results. Statistical (RRKM)/master equation (ME) computations of the product distribution and BFs on these PESs, with inclusion of ISC, have been performed and compared to experiment. In light of the reasonable agreement between the CMB experiment, literature kinetic experimental results, and theoretical predictions for the O(³P) + benzene reaction, the so-validated computational methodology has been used to predict (i) the BF between the C₆H₅O + H and C₅H₆ + CO channels as a function of collision energy and temperature (at 0.1 and 1 bar), showing that their increase progressively favors radical (phenoxy + H)-forming over molecule (C₅H₆ + CO and phenol stabilization)-forming channels, and (ii) channel-specific rate constants as a function of T and p, which are expected to be useful for improved combustion models.

Electronic Structure and Excited States of the Collision Reaction O(3P) + C2H4: A Multiconfigurational Perspective

We present a study of the O(³P) + C₂H₄ scattering reaction, a process that takes place in the interstellar medium and is of relevance in atmospheric chemistry as well. A comprehensive investigation of the electronic properties of the system has been carried out based on multiconfigurational ab initio CASSCF/CASPT2 calculations, using a robust and consistent active space that can deliver accurate potential energy surfaces in the key regions visited by the system. The paper discloses detailed description of the primary reaction pathways and the relevant singlet and triplet excited states at the CASSCF and CASPT2 level, including an accurate description of the critical configurations, such as minima and transition states. The chosen active space and the CASSCF/CASPT2 computational protocol are assessed against coupled-cluster calculations to further check the stability and reliability of the entire multiconfigurational procedure.

The reaction of O(3P) with alkynes: a dynamic and computational study focusing on formyl radical production

Production of formyl radical, HCO, from reactions of O(³P) with alkynes (acetylene, propyne, 1-butyne, and 1-pentyne) has been investigated using cavity ringdown laser absorption spectroscopy (CRDLAS) and computational methods. No HCO was detected from reaction with acetylene, while the amount of HCO increased for propyne and 1-butyne, dropping off somewhat for 1-pentyne. These results differ from trends previously observed for reactions of O(³P) with alkenes, which exhibit the largest HCO production for the smallest alkene and drop off as the alkene size increases. Computational studies employing density functional and coupled cluster methods have been employed to investigate the triplet and singlet state pathways for HCO production. Because intersystem crossing (ISC) has been shown to be important in these processes, the minimum energy crossing point (MECP) between the triplet and singlet surfaces has been studied. We find the MECP for propyne to possess C₁ symmetry and to lie lower in energy than previous studies have found. Natural Bond Orbital and Natural Resonance Theory analyses have been performed to investigate the changes in spin density and bond order along the reaction pathways for formation of HCO. Explanations are suggested for the trend in HCO formation observed for the alkynes. The trend in alkyne HCO yield also is compared and contrasted with the trend previously observed for the alkenes.

Roaming Reactions and Dynamics in the van der Waals Region

Roaming reactions were first clearly identified in photodissociation of formaldehyde 15 years ago, and roaming dynamics are now recognized as a universal aspect of chemical reactivity. These reactions typically involve frustrated near-dissociation of a quasibound system to radical fragments, followed by reorientation at long range and intramolecular abstraction. The consequences can be unexpected formation of molecular products, depletion of the radical pool in chemical systems, and formation of products with unusual internal state distributions. In this review, I examine some current aspects of roaming reactions with an emphasis on experimental results, focusing on possible quantum effects in roaming and roaming dynamics in bimolecular systems. These considerations lead to a more inclusive definition of roaming reactions as those for which key dynamics take place at long range.

A global ab initio potential energy surface and dynamics of the proton-transfer reaction: OH− + D2 → HOD + D−

The reaction mechanisms of OH⁻ + D₂ → HOD + D⁻ were first revealed by theory, based on an accurate full-dimensional PES.

An accurate full-dimensional potential energy surface and quasiclassical trajectory dynamics of the H + H2O2 two-channel reaction

Quasiclassical trajectory calculations reveal interesting dynamics features based on an accurate FI-NN PES for the H + H₂O₂ two-channel reaction.

Crossed beam polyatomic reaction dynamics: recent advances and new insights

This review summarizes the developments in polyatomic reaction dynamics, focusing on reactions of unsaturated hydrocarbons with O-atoms and methane with atoms/radicals.

Theoretical Investigation of the Radical-Radical Reaction of O((3)P) + C2H3 and Comparison with Gas-Phase Crossed-Beam Experiments

Herein, we present an ab initio study of the prototypal radical-radical reactions of ground-state atomic oxygen [O((3)P)] with the vinyl (C2H3) radical using density functional theory and a complete basis set model. Two distinctive pathways on the lowest doublet potential energy surfaces (PESs) were predicted to be in competition: addition and abstraction. The barrierless addition of O((3)P) to the hydrocarbon radicals leads to energy-rich intermediate formation followed by subsequent isomerization and decomposition to yield various products: CH2CO (ketene) + H, CO + CH3, C2HOH (acetylenol) + H, (3,1)CCHOH + H, H2O + C2H, (3,1)CH2 + HCO, H2CO (formaldehyde) + CH, C2H2 (acetylene) + OH, and (3,1)CCH2 + OH. The competing but minor H-atom abstraction mechanisms produce C2H2 + OH and (1,3)CCH2 + OH. The optimized structures of the reactants, products, intermediates, and transition states and the reaction mechanisms were obtained on the lowest doublet PESs. The major pathway was predicted to be the formation of CH2CO + H through the low-barrier, single-step cleavages of the addition intermediates. The Levine-Bernstein prior method, statistical surprisal approach, and microcanonical Rice-Ramsperger-Kassel-Marcus theory were applied to deduce the energy distributions of H atoms and OH products and quantitative rate constants. On the basis of the statistical theory and the population analysis, the predicted energy distributions were compared to the kinetic energy release of H and the preferential population of the Π(A') component of OH products reported in recent gas-phase crossed-beam investigations (Park, M. J.; Jang, S. C.; Choi, J. H. J. Chem. Phys. 2012, 137, 204311), and their kinetic and dynamic characteristics were discussed.

Trajectory surface hopping study of the O((3)P) + C2H2 reaction dynamics: effect of collision energy on the extent of intersystem crossing

Intersystem crossing (ISC) dynamics plays an important role in determining the product branching in the O((3)P) + C2H2 reaction despite the necessarily small spin-orbit coupling constant values. In this study we investigate the effect of collision energy on the extent of the contribution of a spin non-conserving route through ISC dynamics to the product distributions at the initial collision energies 8.2, 9.5, and 13.1 kcal/mol. A direct dynamics trajectory surface hopping method is employed with potential energy surfaces generated at the unrestricted B3LYP/6-31G(d,p) level of theory to perform nonadiabatic dynamics. To make our calculation simpler, nonadibatic transitions were only considered at the triplet-singlet intersections. At the crossing points, Landau-Zener transition probabilities were calculated using spin-orbit coupling constant values computed at the same geometry. The Landau-Zener model for the title reaction is validated against a more rigorous Tully's fewest switches method and found to be working reasonably well as expected because of weak spin-orbit coupling. We have compared our results with the recent crossed molecular beam experiments and observed a very good agreement with respect to the primary product branching ratios. Our calculation revealed that there is no noticeable effect of the initial collision energy on the overall product distributions that corroborates the recent experimental findings. Our calculation indicates, however, that the extent of intersystem crossing contributions varies significantly with collision energy, needed to be verified, experimentally.

A global full-dimensional potential energy surface and quasiclassical trajectory study of the O(1D) + CH4 multichannel reaction

The QCT calculations based on an accurate global full-dimensional PES are capable of reproducing the experimental dynamic features for O(¹D) + CH₄.