Crossed molecular beam experiments and theoretical simulations on the multichannel reaction of toluene with atomic oxygen.

Despite extensive experimental and theoretical studies on the kinetics of the O(3P) + C7H8 (toluene) reaction and a pioneering crossed molecular beam (CMB) investigation, the branching fractions (BFs) of the CH3C6H4O(methylphenoxy) + H, C6H5O(phenoxy) + CH3, and spin-forbidden C5H5CH3 (methylcyclopentadiene) + CO product channels remain an open question, which has hampered the proper inclusion of this important reaction in the chemical modelling of various chemical environments. We report a CMB study with universal soft electron-ionization mass-spectrometric detection of the reactions O(3P,1D) + toluene at the collision energy of 34.7 kJ mol-1. From CMB data we have inferred the reaction dynamics and quantified the BFs of the primary products and the role of intersystem crossing (ISC). The CH3-elimination channel dominates (BF = 0.69 ± 0.22) in the O(3P) reaction, while the H-displacement and CO-formation channels are minor (BF = 0.22 ± 0.07 and 0.09 ± 0.05, respectively), with ISC accounting for more than 50% of the reactive flux. Synergistic transition-state theory (TST)-based master equation simulations including nonadiabatic TST on ab initio coupled triplet/singlet potential energy surfaces were employed to compute the product BFs and assist in the interpretation of the CMB results. In the light of the good agreement between the theoretical predictions for the O(3P) + toluene reaction and the CMB results as well as the absolute rate constant as a function of temperature (T) (from literature), the so-validated computational methodology was used to predict channel-specific rate constants as a function of T at 1 atm.

Unveiling the Reaction Mechanism of the N(2D) + Pyridine Reaction: Ring-Contraction versus 7-Membered-Ring Formation Channels.

Despite the relevance of the reactions of the prototypical nitrogen-containing six-membered aromatic molecule (N-heterocyclic) of pyridine (C6H5N) in environmental science, astrochemistry, planetary science, prebiotic chemistry, and materials science, few experimental/theoretical studies exist on the bimolecular reactions involving pyridine and neutral atomic/molecular radicals. We report a combined experimental and theoretical study on the elementary reaction of pyridine with excited nitrogen atoms, N(2D), aimed at providing information about the primary reaction products and their branching fractions (BFs). From previous crossed molecular beam (CMB) experiments with mass-spectrometric detection and present synergistic calculations of the reactive potential energy surface (PES) and product BFs we have unveiled the reaction mechanism. It is found that the reaction proceeds via N(2D) barrierless addition to pyridine that, via bridged intermediates followed by N atom "sliding" into the ring, leads to 7-membered-ring structures. They further evolve, mainly via ring-contraction mechanisms toward 5-membered-ring radical products and, to a smaller extent, via H-displacement mechanisms toward 7-membered-ring isomeric products and their isomers. Using the theoretical statistical estimates, an improved fit of the experimental data previously reported has been obtained, leading to the following results for the dominant product channels: C4H4N (pyrrolyl) + HCN (BF = 0.61 ± 0.20), C3H3N2 (1H-imidazolyl/1H-pyrazolyl) + C2H2 (BF = 0.11 ± 0.06), and C5H4N2 (7-membered-ring molecules or pyrrole carbonitriles) + H (BF = 0.28 ± 0.10). The ring-contraction product channels C4H4N (pyrrolyl) + HCN, C3H3N2 (1H-imidazolyl) + C2H2, C3H3N2 (1H-pyrazolyl) + C2H2, and C5H5 (cyclopentadienyl) + N2 have statistical BFs of 0.54, 0.09, 0.11, and 0.07, respectively. Among the H-displacement channels, the cyclic-CHCHCHCHNCN + H channel and cyclic-CHCHCHCHCN2 + H are theoretically predicted to have a comparable BF (0.07 and 0.06, respectively), while the other isomeric 7-membered-ring molecule + H channel has a BF of 0.03. Pyrrole-carbonitriles and 1H-ethynyl-1H-imidazole (+ H) isomeric channels have an overall BF of 0.03. Implications for the chemistry of Saturn's moon Titan and prebiotic chemistry, as well as for understanding the N-doping of graphene or carbon nanotubes, are noted.

An experimental and theoretical investigation of the N(2D) + C6H6 (benzene) reaction with implications for the photochemical models of Titan.

We report on a combined experimental and theoretical investigation of the N(2D) + C6H6 (benzene) reaction, which is of relevance in the aromatic chemistry of the atmosphere of Titan. Experimentally, the reaction was studied (i) under single-collision conditions by the crossed molecular beams (CMB) scattering method with mass spectrometric detection and time-of-flight analysis at the collision energy (Ec) of 31.8 kJ mol-1 to determine the primary products, their branching fractions (BFs), and the reaction micromechanism, and (ii) in a continuous supersonic flow reactor to determine the rate constant as a function of temperature from 50 K to 296 K. Theoretically, electronic structure calculations of the doublet C6H6N potential energy surface (PES) were performed to assist the interpretation of the experimental results and characterize the overall reaction mechanism. The reaction is found to proceed via barrierless addition of N(2D) to the aromatic ring of C6H6, followed by formation of several cyclic (five-, six-, and seven-membered ring) and linear isomeric C6H6N intermediates that can undergo unimolecular decomposition to bimolecular products. Statistical estimates of product BFs on the theoretical PES were carried out under the conditions of the CMB experiments and at the temperatures relevant for Titan's atmosphere. In all conditions the ring-contraction channel leading to C5H5 (cyclopentadienyl) + HCN is dominant, while minor contributions come from the channels leading to o-C6H5N (o-N-cycloheptatriene radical) + H, C4H4N (pyrrolyl) + C2H2 (acetylene), C5H5CN (cyano-cyclopentadiene) + H, and p-C6H5N + H. Rate constants (which are close to the gas kinetic limit at all temperatures, with the recommended value of 2.19 ± 0.30 × 10-10 cm3 s-1 over the 50-296 K range) and BFs have been used in a photochemical model of Titan's atmosphere to simulate the effect of the title reaction on the species abundances as a function of the altitude.

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

Acrylonitrile (CH2CHCN) 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(3P) + CH2CHCN 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 3P and excited 1D atomic oxygen. From product angular and time-of-flight distributions in crossed molecular beam experiments with mass spectrometric detection at a collision energy, Ec, 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(3P) reaction leads to two main product channels: (i) CH2CNH (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 CH2CHCN and a possible formation route of CH2CNH in the interstellar medium. The O(1D) + CH2CHCN reaction mainly leads to the formation of CH2CNH + CO adiabatically on the singlet PES. This result can improve models related to the chemistry of interstellar ice and cometary comas, where O(1D) reactions can play a role. Overall, our results are expected to be useful for improving the models of combustion and extraterrestrial environments.

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 (HCnN, 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 3P and excited 1D 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(3P, 1D) with HCCCN(X1Σ+) 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, Ec, 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(3P) 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 1HCCN (cyanomethylene) + CO. The O(1D) 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(3P) + HC3N at 300 K is in the 10-12 cm3 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.

OH(2Π) + C2H4 Reaction: A Combined Crossed Molecular Beam and Theoretical Study.

The reaction between the ground-state hydroxyl radical, OH(2Π), and ethylene, C2H4, has been investigated under single-collision conditions by the crossed molecular beam scattering technique with mass-spectrometric detection and time-of-flight analysis at the collision energy of 50.4 kJ/mol. Electronic structure calculations of the underlying potential energy surface (PES) and statistical Rice-Ramsperger-Kassel-Marcus (RRKM) calculations of product branching fractions on the derived PES for the addition pathway have been performed. The theoretical results indicate a temperature-dependent competition between the anti-/syn-CH2CHOH (vinyl alcohol) + H, CH3CHO (acetaldehyde) + H, and H2CO (formaldehyde) + CH3 product channels. The yield of the H-abstraction channel could not be quantified with the employed methods. The RRKM results predict that under our experimental conditions, the anti- and syn-CH2CHOH + H product channels account for 38% (in similar amounts) of the addition mechanism yield, the H2CO + CH3 channel for ∼58%, while the CH3CHO + H channel is formed in negligible amount (<4%). The implications for combustion and astrochemical environments are discussed.

The N(2D) + CH2CHCN (Vinyl Cyanide) Reaction: A Combined Crossed Molecular Beam and Theoretical Study and Implications for the Atmosphere of Titan.

The reaction of electronically excited nitrogen atoms, N(2D), with vinyl cyanide, CH2CHCN, has been investigated under single-collision conditions by the crossed molecular beam (CMB) scattering method with mass spectrometric detection and time-of-flight (TOF) analysis at the collision energy, Ec, of 31.4 kJ/mol. Synergistic electronic structure calculations of the doublet potential energy surface (PES) have been performed to assist in the interpretation of the experimental results and characterize the overall reaction micromechanism. Statistical (Rice-Ramsperger-Kassel-Marcus, RRKM) calculations of product branching fractions (BFs) on the theoretical PES have been carried out at different values of temperature, including the one corresponding to the temperature (175 K) of Titan's stratosphere and at a total energy corresponding to the Ec of the CMB experiment. According to our theoretical calculations, the reaction is found to proceed via barrierless addition of N(2D) to the carbon-carbon double bond of CH2═CH-CN, followed by the formation of cyclic and linear intermediates that can undergo H, CN, and HCN elimination. In competition, the N(2D) addition to the CN group is also possible via a submerged barrier, leading ultimately to N2 + C3H3 formation, the most exothermic of all possible channels. Product angular and TOF distributions have been recorded for the H-displacement channels leading to the formation of a variety of possible C3H2N2 isomeric products. Experimentally, no evidence of CN, HCN, and N2 forming channels was observed. These findings were corroborated by the theory, which predicts a variety of competing product channels, following N(2D) addition to the double bond, with the main ones, at Ec = 31.4 kJ/mol, being six isomeric H forming channels: c-CH(N)CHCN + H (BF = 35.0%), c-CHNCHCN + H (BF = 28.1%), CH2NCCN + H (BF = 26.3%), c-CH2(N)CCN(cyano-azirine) + H (BF = 7.4%), trans-HNCCHCN + H (BF = 1.6%), and cis-HNCCHCN + H (BF = 1.3%), while C-C bond breaking channels leading to c-CH2(N)CH(2H-azirine) + CN and c-CH2(N)C + HCN are predicted to be negligible (0.02% and 0.2%, respectively). The highly exothermic N2 + CH2CCH (propargyl) channel is also predicted to be negligible because of the very high isomerization barrier from the initial addition intermediate to the precursor intermediate able to lead to products. The predicted product BFs are found to have, in general, a very weak energy dependence. The above cyclic and linear products containing an additional C-N bond could be potential precursors of more complex, N-rich organic molecules that contribute to the formation of the aerosols on Titan's upper atmosphere. Overall, the results are expected to have a significant impact on the gas-phase chemistry of Titan's atmosphere and should be properly included in the photochemical models.

Unsaturated Dinitriles Formation Routes in Extraterrestrial Environments: A Combined Experimental and Theoretical Investigation of the Reaction between Cyano Radicals and Cyanoethene (C2H3CN).

The reaction between cyano radicals (CN, X2Σ+) and cyanoethene (C2H3CN) has been investigated by a combined approach coupling crossed molecular beam (CMB) experiments with mass spectrometric detection and time-of-flight analysis at a collision energy of 44.6 kJ mol-1 and electronic structure calculations to determine the relevant potential energy surface. The experimental results can be interpreted by assuming the occurrence of a dominant reaction pathway leading to the two but-2-enedinitrile (1,2-dicyanothene) isomers (E- and Z-NC-CH═CH-CN) in a H-displacement channel and, to a much minor extent, to 1,1-dicyanoethene, CH2C(CN)2. In order to derive the product branching ratios under the conditions of the CMB experiments and at colder temperatures, including those relevant to Titan and to cold interstellar clouds, we have carried out RRKM statistical calculations using the relevant potential energy surface of the investigated reaction. We have also estimated the rate coefficient at very low temperatures by employing a semiempirical method for the treatment of long-range interactions. The reaction has been found to be barrierless and fast also under the low temperature conditions of cold interstellar clouds and the atmosphere of Titan. Astrophysical implications and comparison with literature data are also presented. On the basis of the present work, 1,2-dicyanothene and 1,1-dicyanothene are excellent candidates for the search of dinitriles in the interstellar medium.

The Reaction N(2D) + CH3CCH (Methylacetylene): A Combined Crossed Molecular Beams and Theoretical Investigation and Implications for the Atmosphere of Titan.

The reaction of excited nitrogen atoms N(2D) with CH3CCH (methylacetylene) was investigated under single-collision conditions by the crossed molecular beams (CMB) scattering method with mass spectrometric detection and time-of-flight analysis at the collision energy (Ec) of 31.0 kJ/mol. Synergistic electronic structure calculations of the doublet potential energy surface (PES) were performed to assist the interpretation of the experimental results and characterize the overall reaction micromechanism. Theoretically, the reaction is found to proceed via a barrierless addition of N(2D) to the carbon-carbon triple bond of CH3CCH and an insertion of N(2D) into the CH bond of the methyl group, followed by the formation of cyclic and linear intermediates that can undergo H, CH3, and C2H elimination or isomerize to other intermediates before unimolecularly decaying to a variety of products. Kinetic calculations for addition and insertion mechanisms and statistical (Rice-Ramsperger-Kassel-Marcus) computations of product branching fractions (BFs) on the theoretical PES were performed at different values of total energy, including the one corresponding to the temperature (175 K) of Titan's stratosphere and that of the CMB experiment. Up to 14 competing product channels were statistically predicted, with the main ones, at Ec = 31.0 kJ/mol, being the formation of CH2NH (methanimine) + C2H (ethylidyne) (BF = 0.41), c-C(N)CH + CH3 (BF = 0.32), CH2CHCN (acrylonitrile) + H (BF = 0.12), and c-CH2C(N)CH + H (BF = 0.04). Of the 14 possible channels, seven correspond to H displacement channels of different exothermicity, for a total H channel BF of ∼0.25 at Ec = 31.0 kJ/mol. Experimentally, dynamical information could only be obtained about the overall H channels. In particular, the experiment corroborates the formation of acrylonitrile + H, which is the most exothermic of all 14 reaction channels and is theoretically calculated to be the dominant H-forming channel (BF = 0.12). The products containing a novel C-N bond could be potential precursors to form other nitriles (C2N2, C3N) or more complex organic species containing N atoms in planetary atmospheres, such as those of Titan and Pluto. Overall, the results are expected to have a potentially significant impact on the understanding of the gas-phase chemistry of Titan's atmosphere and the modeling of that atmosphere.

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(3P). During the past decade, we have performed in-depth investigations of the reactions of O(3P) 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 (C6H6) molecule. Despite extensive experimental and theoretical studies on the kinetics and dynamics of the O(3P) + C6H6 reaction, the relative importance of the C6H5O (phenoxy) + H open-shell products and of the spin-forbidden C5H6 (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(3P) + benzene at a collision energy (Ec) of 8.2 kcal/mol, focusing on the occurrence of the phenoxy + H and spin-forbidden C5H6 + 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(1D) + benzene at the same Ec. 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(3P) + benzene reaction, the so-validated computational methodology has been used to predict (i) the BF between the C6H5O + H and C5H6 + 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 (C5H6 + 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.

Combined Experimental and Theoretical Studies of the O(3P) + 1-Butene Reaction Dynamics: Primary Products, Branching Fractions, and Role of Intersystem Crossing.

Information on the detailed mechanism and dynamics (primary products, branching fractions (BFs), and channel specific rate constants as a function of temperature) for many important combustion reactions of O(3P) with unsaturated hydrocarbons is still lacking. We report synergistic experimental/theoretical studies on the mechanism and dynamics of the O(3P) + 1-C4H8 (1-butene) reaction by combining crossed molecular beam (CMB) experiments with soft electron ionization mass-spectrometric detection and time-of-flight analysis at 10.5 kcal/mol collision energy (Ec) to high-level ab initio electronic structure calculations of the underlying triplet and singlet potential energy surfaces (PESs) and statistical Rice-Ramsperger-Kassel-Marcus/Master Equation (RRKM/ME) computations of BFs including intersystem crossing (ISC). The reactive interaction of O(3P) with 1-butene is found to mainly break apart the 4-carbon atom chain, leading to the radical product channels ethyl + vinoxy (BF = 0.34 ± 0.11), methyl + C3H5O (BF = 0.28 ± 0.09), formyl + propyl (BF = 0.17 ± 0.05), ethyl + acetyl (BF = 0.014 ± 0.007), and butanal radical (ethylvinoxy) + H (BF = 0.013 ± 0.004), and molecular product channels formaldehyde + propenylidene/propene (BF = 0.15 ± 0.05) and butenone (ethyl ketene) + H2 (BF = 0.037 ± 0.015). As some of these products can only be formed via ISC from triplet to singlet PESs, from BFs an extent of ISC of 50% is inferred. This value is significantly larger than that recently observed for O(3P) + propene (22%) at similar Ec, underlying the question of how important it is to consider nonadiabatic effects for these and similar combustion reactions. Comparison of the derived BFs with those of statistical (RRKM/ME) simulations on the ab initio coupled triplet/singlet PESs shows good agreement, warranting the use of the RRKM/ME approach to provide information on the variation of BFs with temperature and to derive channel specific rate constants as a function of temperature (T) and pressure (p). Notably, ISC is predicted to decrease strongly with increasing temperature (from about 70% at 300 K to 46% at Ec = 10.5 kcal/mol, and about 1% at 2000 K). The present results lead to a detailed understanding of the complex reaction mechanism of O(3P) + 1-butene and, by providing channel specific rate constants as a function of T and p, should facilitate the improvement of current fossil-fuel (1-butene) as well as biofuel (1-butanol) combustion models.

Crossed beam polyatomic reaction dynamics: recent advances and new insights.

Over the past ten years or so, great advances in our understanding of the dynamics of elementary (bimolecular) polyatomic reactions in the gas-phase have occurred. This has been made possible by critical improvements (a) in crossed molecular beam (CMB) instruments with rotating mass spectrometric detection and time-of-flight analysis, especially following the implementation of soft ionization (by tunable low energy electrons or vacuum-ultraviolet synchrotron radiation) for product detection with increased sensitivity and universal detection power, and (b) in REMPI-slice velocity map ion imaging (VMI) detection techniques in pulsed CMB experiments for obtaining product pair-correlated information through high-resolution measurements directly in the center of mass system. The improved universal CMB method is permitting us to identify all primary reaction products, characterize their formation dynamics, and determine the branching ratios (BRs) for multichannel non-adiabatic reactions, such as those of ground state oxygen atoms, O(3P), with unsaturated hydrocarbons (alkynes, alkenes, dienes). The improved slice VMI CMB technique is permitting us to explore at an unprecedented level of detail, through pair-correlated measurements, the reaction dynamics of a prototype polyatomic molecule such as CH4 (and isotopologues) in its ground state with a variety of important X radicals such as F, Cl, O, and OH. In this review, we highlight this recent progress in the field of CMB reaction dynamics, with an emphasis on the experimental side, but with the related theoretical work, at the level of state-of-the-art calculations of both the underlying potential energy surfaces and the reaction dynamics, noted throughout. In particular, the focus is (a) on the effect of molecular complexity and structure on product distributions, branching ratios and role of intersystem crossing for the multichannel, addition-elimination reactions of unsaturated hydrocarbons with O atoms, and (b) on the very detailed dynamics of the abstraction reactions of ground-state methane (and isotopologues) with atoms (F, Cl, O) and diatoms (OH), with inclusion of also rotational mode specificity in the vibrationally excited methane reactions.

Reaction Dynamics of O((3)P) + Propyne: II. Primary Products, Branching Ratios, and Role of Intersystem Crossing from Ab Initio Coupled Triplet/Singlet Potential Energy Surfaces and Statistical Calculations.

The mechanism of the O((3)P) + CH3CCH reaction was investigated using a combined experimental/theoretical approach. Experimentally the reaction dynamics was studied using crossed molecular beams (CMB) with mass-spectrometric detection and time-of-flight analysis at 9.2 kcal/mol collision energy. Theoretically master equation (ME) simulations were performed on a potential energy surface (PES) determined using high-level ab initio electronic structure calculations. In this paper (II) the theoretical results are described and compared with experiments, while in paper (I) are reported and discussed the results of the experimental study. The PES was investigated by determining structures and vibrational frequencies of wells and transition states at the CASPT2/aug-cc-pVTZ level using a minimal active space. Energies were then determined at the CASPT2 level increasing systematically the active space and at the CCSD(T) level extrapolating to the complete basis set limit. Two separate portions of the triplet PES were investigated, as O((3)P) can add either on the terminal or the central carbon of the unsaturated propyne bond. Minimum energy crossing points (MECPs) between the triplet and singlet PESs were searched at the CASPT2 level. The calculated spin-orbit coupling constants between the T1 and S0 electronic surfaces were ∼25 cm(-1) for both PESs. The portions of the singlet PES that can be accessed from the MECPs were investigated at the same level of theory. The system reactivity was predicted integrating stochastically the one-dimensional ME using Rice-Ramsperger-Kassel-Marcus theory to determine rate constants on the full T1/S0 PESs, accounting explicitly for intersystem crossing (ISC) using the Landau-Zener model. The computational results are compared both with the branching ratios (BRs) determined experimentally in the companion paper (I) and with those estimated in a recent kinetic study at 298 K. The ME results allow to interpret the main system reactivity: CH3CCO + H and CH3 + HCCO are the major channels active on the triplet PES and are formed from the wells accessed after O addition to the terminal and central C, respectively; (1)CH3CH + CO and C2H3 + HCO are the major channels active on the singlet PES and are formed from the methylketene and acrolein wells after ISC. However, also a large number of minor channels (∼15) are active, so that the system reactivity is quite complicated. The comparison between computational and experimental BRs is quite good for the kinetic study, while some discrepancy with the CMB estimations suggests that dynamic non-ergodic effects may influence the system reactivity. Channel specific rate constants are calculated in the 300-2250 K and 1-30 bar temperature and pressure ranges. It is found that as the temperature increases the H abstraction reaction, whose contribution is negligible in the experimental conditions, increases in relevance, and the extent of ISC decreases from ∼80% at 300 K to less than 2% at 2250 K.

Reaction Dynamics of O((3)P) + Propyne: I. Primary Products, Branching Ratios, and Role of Intersystem Crossing from Crossed Molecular Beam Experiments.

We performed synergic experimental/theoretical studies on the mechanism of the O((3)P) + propyne reaction by combining crossed molecular beams experiments with mass-spectrometric detection and time-of-flight analysis at 9.2 kcal/mol collision energy (Ec) with ab initio electronic structure calculations at a high level of theory of the relevant triplet and singlet potential energy surfaces (PESs) and statistical calculations of branching ratios (BRs) taking into account intersystem crossing (ISC). In this paper (I) we report the results of the experimental investigation, while the accompanying paper (II) shows results of the theoretical investigation with comparison to experimental results. By exploiting soft electron ionization detection to suppress/mitigate the effects of the dissociative ionization of reactants, products, and background gases, product angular and velocity distributions at different charge-to-mass ratios were measured. From the laboratory data angular and translational energy distributions in the center-of-mass system were obtained for the five competing most important product channels, and product BRs were derived. The reactive interaction of O((3)P) with propyne under single collision conditions is mainly leading to the rupture of the three-carbon atom chain, with production of the radical products methylketenyl + atomic hydrogen (BR = 0.04), methyl + ketenyl (BR = 0.10), and vinyl + formyl (BR = 0.11) and the molecular products ethylidene/ethylene + carbon monoxide (BR = 0.74) and propandienal + molecular hydrogen (BR = 0.01). Because some of the products can only be formed via ISC from the entrance triplet to the low-lying singlet PES, we infer from their BRs an amount of ISC larger than 80%. This value is dramatically large when compared to the negligible ISC reported for the O((3)P) reaction with the simplest alkyne, acetylene. At the same time, it is much larger than that (∼20%) recently observed in the related reaction of the three-carbon atom alkene, O((3)P) + propene at a comparable Ec. This poses the question of how important it is to consider nonadiabatic effects and their dependence on molecular structure for this kind of combustion reactions. The prevalence of the CH3 over the H displacement channels is not explained by invoking a preference for the addition on the methyl-substituted acetylenic carbon atom, but rather it is believed to be due to the different tendencies of the two addition triplet intermediates CH3CCHO (preferentially leading to H elimination) and CH3COCH (preferentially leading to CH3 elimination) to undergo ISC to the underlying singlet PES. It is concluded that the main coproduct of the CO forming channel is singlet ethylidene ((1)CH3CH) rather than ground-state ethylene. By comparing the derived BRs with those very recently derived from kinetics studies at room temperature and 4 Torr we obtained information on how BRs vary with collision energy. The extent of ISC is estimated to remain essentially constant (∼85%) with increasing Ec from ∼1 to ∼10 kcal/mol. The present experimental results shed light on the mechanism of the title reaction at energies comparable to those involved in combustion and, when compared with the results from the statistical calculations on the ab initio coupled PESs (see accompanying paper II), lead to an in-depth understanding of the rather complex reaction mechanism of O + propyne. The overall results are expected to contribute to the development of more refined models of hydrocarbon combustion.

Crossed Molecular Beams and Quasiclassical Trajectory Surface Hopping Studies of the Multichannel Nonadiabatic O((3)P) + Ethylene Reaction at High Collision Energy.

The combustion relevant O((3)P) + C2H4 reaction stands out as a prototypical multichannel nonadiabatic reaction involving both triplet and singlet potential energy surfaces (PESs), which are strongly coupled. Crossed molecular beam (CMB) scattering experiments with universal soft electron ionization mass spectrometric detection have been used to characterize the dynamics of this reaction at the relatively high collision energy Ec of 13.7 kcal/mol, attained by crossing the reactant beams at an angle of 135°. This work is a full report of the data at the highest Ec investigated for this reaction. From laboratory product angular and velocity distribution measurements, angular and translational energy distributions in the center-of-mass system have been obtained for the five observed exothermic competing reaction channels leading to H + CH2CHO, H + CH3CO, CH3 + HCO, CH2 + H2CO, and H2 + CH2CO. The product branching ratios (BRs) have been derived. The elucidation of the reaction dynamics is assisted by synergic full-dimensional quasiclassical trajectory surface-hopping calculations of the reactive differential cross sections on coupled ab initio triplet/singlet PESs. This joint experimental/theoretical study extends and complements our previous combined CMB and theoretical work at the lower collision energy of 8.4 kcal/mol. The theoretically derived BRs and extent of intersystem crossing (ISC) are compared with experimental results. In particular, the predictions of the QCT results for the three main channels (those leading to vinoxy + H, methyl + HCO and methylene + H2CO formation) are compared directly with the experimental data in the laboratory frame. Good overall agreement is noted between theory and experiment, although some small, yet significant shortcomings of the theoretical differential cross section are noted. Both experiment and theory find almost an equal contribution from the triplet and singlet surfaces to the reaction, with a clear tendency of the degree of ISC to decrease with increasing Ec and with theory slightly overestimating the extent of ISC.

Intersystem crossing and dynamics in O(3P) + C2H4 multichannel reaction: experiment validates theory.

The O((3)P) + C(2)H(4) reaction, of importance in combustion and atmospheric chemistry, stands out as a paradigm reaction involving triplet- and singlet-state potential energy surfaces (PESs) interconnected by intersystem crossing (ISC). This reaction poses challenges for theory and experiments owing to the ruggedness and high dimensionality of these potentials, as well as the long lifetimes of the collision complexes. Primary products from five competing channels (H + CH(2)CHO, H + CH(3)CO, H(2) + CH(2)CO, CH(3) + HCO, CH(2) + CH(2)O) and branching ratios (BRs) are determined in crossed molecular beam experiments with soft electron-ionization mass-spectrometric detection at a collision energy of 8.4 kcal/mol. As some of the observed products can only be formed via ISC from triplet to singlet PESs, from the product BRs the extent of ISC is inferred. A new full-dimensional PES for the triplet state as well as spin-orbit coupling to the singlet PES are reported, and roughly half a million surface hopping trajectories are run on the coupled singlet-triplet PESs to compare with the experimental BRs and differential cross-sections. Both theory and experiment find almost equal contributions from the two PESs to the reaction, posing the question of how important is it to consider the ISC as one of the nonadiabatic effects for this and similar systems involved in combustion chemistry. Detailed comparisons at the level of angular and translational energy distributions between theory and experiment are presented for the two primary channel products, CH(3) + HCO and H + CH(2)CHO. The agreement between experimental and theoretical functions is excellent, implying that theory has reached the capability of describing complex multichannel nonadiabatic reactions.

Dynamics of the O((3)P) + C2H2 reaction from crossed molecular beam experiments with soft electron ionization detection.

The reaction between ground state oxygen atoms, O((3)P), and the acetylene molecule, C2H2, has been investigated in crossed molecular beam experiments with mass-spectrometric detection and time-of-flight analysis at three different collision energies, Ec = 34.4, 41.1 and 54.6 kJ mol(-1). From product angular and velocity distribution measurements of the HCCO and CH2 products in the laboratory frame, product angular and translational energy distributions in the center-of-mass frame were determined. Measurements on the CH2 product were made possible by employing for product detection the recently implemented soft electron-ionization (EI) technique with low-energy, tunable electrons, which has permitted suppressing interference coming from the dissociative ionization of reactants, products and background gases. It has been found that the title reaction leads only to two competing channels: H + HCCO (ketenyl) and CO + CH2 (triplet methylene). The branching ratio of cross sections between the two competing channels has been determined to be σ(HCCO)/[σ(HCCO) + σ(CH2)] = 0.79 ± 0.05, independent of collision energy within the experimental uncertainty. This value is in line with that obtained in the most recent and accurate kinetics determination at room temperature as well as with that predicted from recent theoretical calculations based on statistical rate theory and weak-collision master equation analysis and on dynamics surface-hopping quasiclassical trajectory calculations on-the-fly on coupled triplet/singlet ab initio potential energy surfaces. The firm assessment of the branching ratio as a function of translational energy for this important reaction, besides its fundamental significance, is of considerable relevance for the implementation of theoretical models of hydrocarbon combustion.

Combined crossed beam and theoretical studies of the C(1D) + CH4 reaction.

The reaction involving atomic carbon in its first electronically excited state (1)D and methane has been investigated in crossed molecular beam experiments at a collision energy of 25.3 kJ mol(-1). Electronic structure calculations of the underlying potential energy surface (PES) and Rice-Ramsperger-Kassel-Marcus (RRKM) estimates of rates and branching ratios have been performed to assist the interpretation of the experimental results. The reaction proceeds via insertion of C((1)D) into one of the C-H bonds of methane leading to the formation of the intermediate HCCH(3) (methylcarbene or ethylidene), which either decomposes directly into the products C(2)H(3) + H or C(2)H(2) + H(2) or isomerizes to the more stable ethylene, which in turn dissociates into C(2)H(3) + H or H(2)CC + H(2). The experimental results indicate that the H-displacement and H(2)-elimination channels are of equal importance and that for both channels the reaction mechanism is controlled by the presence of a bound intermediate, the lifetime of which is comparable to its rotational period. On the contrary, RRKM estimates predict a very short lifetime for the insertion intermediate and the dominance of the H-displacement channel. It is concluded that the reaction C((1)D) + CH(4) cannot be described statistically and a dynamical treatment is necessary to understand its mechanism. Possibly, nonadiabatic effects are responsible for the discrepancies, as triplet and singlet PES of methylcarbene cross each other and intersystem crossing is possible. Similarities with the photodissociation of ethylene and with the related reactions N((2)D) + CH(4), O((1)D) + CH(4) and S((1)D) + CH(4) are also commented on.