Crossed beam investigations of the reaction dynamics of O(3P) with allyl radical, C3H5

ChemDyME: Kinetically Steered, Automated Mechanism Generation through Combined Molecular Dynamics and Master Equation Calculations

In many scientific fields, there is an interest in understanding the way in which chemical networks evolve. The chemical networks which researchers focus upon have become increasingly complex, and this has motivated the development of automated methods for exploring chemical reactivity or conformational change in a "black-box" manner, harnessing modern computing resources to automate mechanism discovery. In this work, we present a new approach to automated mechanism generation which couples molecular dynamics and statistical rate theory to automatically find kinetically important reactions and then solve the time evolution of the species in the evolving network. The key to this chemical network mapping through combined dynamics and ME simulation approach is the concept of "kinetic convergence", whereby the search for new reactions is constrained to those species which are kinetically favorable at the conditions of interest. We demonstrate the capability of the new approach for two systems, a well-studied combustion system and a multiple oxygen addition system relevant to atmospheric aerosol formation.

A combined crossed-beam and theoretical study of the reaction dynamics of O(3P) + C2H3 → C2H2 + OH: analysis of the nascent OH products with the preferential population of the Π(A') component

The gas-phase reaction dynamics of ground-state atomic oxygen [O((3)P) from the photo-dissociation of NO(2)] with vinyl radicals [C(2)H(3) from the supersonic flash pyrolysis of vinyl iodide, C(2)H(3)I] has been investigated using a combination of high-resolution laser-induced fluorescence spectroscopy in a crossed-beam configuration and ab initio calculations. Unlike the previous gas-phase bulk kinetic experiments by Baulch et al. [J. Phys. Chem. Ref. Data 34, 757 (2005)], a new exothermic channel of O((3)P) + C(2)H(3) → C(2)H(2) + OH (X (2)Π: υ" = 0) has been identified for the first time, and the population analysis shows bimodal nascent rotational distributions of OH products with low- and high-N" components with a ratio of 2.4:1. No spin-orbit propensities were observed, and the averaged ratios of Π(A('))∕Π(A") were determined to be 1.66 ± 0.27. On the basis of computations at the CBS-QB3 theory level and comparison with prior theory, the microscopic mechanisms responsible for the nascent populations can be understood in terms of two competing dynamical pathways: a direct abstraction process in the low-N" regime as the major pathway and an addition-complex forming process in the high-N" regime as the minor pathway. Particularly, during the bond cleavage process of the weakly bound van der Waals complex C(2)H(2)-OH, the characteristic pathway from the low dihedral-angle geometry was consistent with the observed preferential population of the Π(A') component in the nascent OH products. A molecular-level discussion of the reactivity, mechanism, and dynamical features of the title reaction are presented together with a comparison to gas-phase oxidation reactions of a series of prototypical hydrocarbon radicals.

Probing the dynamics of polyatomic multichannel elementary reactions by crossed molecular beam experiments with soft electron-ionization mass spectrometric detection

In this Perspective we highlight developments in the field of chemical reaction dynamics. Focus is on the advances recently made in the investigation of the dynamics of elementary multichannel radical-molecule and radical-radical reactions, as they have become possible using an improved crossed molecular beam scattering apparatus with universal electron-ionization mass spectrometric detection and time-of-flight analysis. These improvements consist in the implementation of (a) soft ionization detection by tunable low-energy electrons which has permitted us to reduce interfering signals originating from dissociative ionization processes, usually representing a major complication, (b) different beam crossing-angle set-ups which have permitted us to extend the range of collision energies over which a reaction can be studied, from very low (a few kJ mol(-1), as of interest in astrochemistry or planetary atmospheric chemistry) to quite high energies (several tens of kJ mol(-1), as of interest in high temperature combustion systems), and (c) continuous supersonic sources for producing a wide variety of atomic and molecular radical reactant beams. Exploiting these new features it has become possible to tackle the dynamics of a variety of polyatomic multichannel reactions, such as those occurring in many environments ranging from combustion and plasmas to terrestrial/planetary atmospheres and interstellar clouds. By measuring product angular and velocity distributions, after having suppressed or mitigated, when needed, the problem of dissociative ionization of interfering species (reactants, products, background gases) by soft ionization detection, essentially all primary reaction products can be identified, the dynamics of each reaction channel characterized, and the branching ratios determined as a function of collision energy. In general this information, besides being of fundamental relevance, is required for a predictive description of the chemistry of these environments via computer models. Examples are taken from recent on-going work (partly published) on the reactions of atomic oxygen with acetylene, ethylene and allyl radical, of great importance in combustion. A reaction of relevance in interstellar chemistry, as that of atomic carbon with acetylene, is also discussed briefly. Comparison with theoretical results is made wherever possible, both at the level of electronic structure calculations of the potential energy surfaces and dynamical computations. Recent complementary CMB work as well as kinetic work exploiting soft photo-ionization with synchrotron radiation are noted. The examples illustrated in this article demonstrate that the type of dynamical results now obtainable on polyatomic multichannel radical-molecule and radical-radical reactions might well complement reaction kinetics experiments and hence contribute to bridging the gap between microscopic reaction dynamics and thermal reaction kinetics, enhancing significantly our basic knowledge of chemical reactivity and understanding of the elementary reactions which occur in real-world environments.

Crossed molecular beam studies of the reactions of allyl radicals, C3H5(X2A2), with methylacetylene (CH3CCH(X1A1)), allene (H2CCCH2(X1A1)), and their isotopomers

The chemical dynamics of the reaction of allyl radicals, C(3)H(5)(X(2)A(2)), with two C(3)H(4) isomers, methylacetylene (CH(3)CCH(X(1)A(1))) and allene (H(2)CCCH(2)(X(1)A(1))) together with their (partially) deuterated counterparts, were unraveled under single-collision conditions at collision energies of about 125 kJ mol(-1) utilizing a crossed molecular beam setup. The experiments indicate that the reactions are indirect via complex formation and proceed via an addition of the allyl radical with its terminal carbon atom to the terminal carbon atom of the allene and of methylacetylene (alpha-carbon atom) to form the intermediates H(2)CCHCH(2)CH(2)CCH(2) and H(2)CCHCH(2)CHCCH(3), respectively. The lifetimes of these intermediates are similar to their rotational periods but too short for a complete energy randomization to occur. Experiments with D4-allene and D4-methylacetylene verify explicitly that the allyl group stays intact: no hydrogen emission was observed but only the release of deuterium atoms from the perdeuterated reactants. Further isotopic substitution experiments with D3-methylacetylene combined with the nonstatistical nature of the reaction suggest that the intermediates decompose via hydrogen atom elimination to 1,3,5-hexatriene, H(2)CCHCH(2)CHCCH(2), and 1-hexen-4-yne, H(2)CCHCH(2)CCCH(3), respectively, via tight exit transition states located about 10-15 kJ mol(-1) above the separated products. The overall reactions were found to be endoergic by 98 +/- 4 kJ mol(-1) and have characteristic threshold energies to reaction between 105 and 110 kJ mol(-1). Implications of these findings to combustion and interstellar chemistry are discussed.

Crossed beam studies of radical–radical reactions: O(3P) + C3H5(allyl)

The dynamics of the radical-radical reaction O((3)P) + C(3)H(5) has been investigated by means of the crossed molecular beam technique with mass spectrometric detection at a collision energy of 73.0 kJ mol(-1); the reaction mechanism of the H-displacement channel has been elucidated, while experimental evidence of the occurrence of one or more C-C bond-breaking channels at this collision energy has been obtained.

A Theoretical Investigation of the Gas-Phase Oxidation Reaction of the Saturated tert-Butyl Radical

The radical-radical reaction mechanisms and dynamics of ground-state atomic oxygen [O(3P)] with the saturated tert-butyl radical (t-C4H9) are investigated using the density functional method and the complete basis set model. Two distinctive reaction pathways are predicted to be in competition: addition and abstraction. The barrierless addition of O(3P) to t-C4H9 leads to the formation of an energy-rich intermediate (OC4H9) on the lowest doublet potential energy surface, which undergoes subsequent direct elimination or isomerization-elimination leading to various products: C3H6O + CH3, iso-C4H8O + H, C3H7O + CH2, and iso-C4H8 + OH. The respective microscopic reaction processes examined with the aid of statistical calculations, predict that the major addition pathway is the formation of acetone (C3H6O) + CH3 through a low-barrier, single-step cleavage. For the direct, barrierless H-atom abstraction mechanism producing iso-C4H8 (isobutene) + OH, which was recently reported in gas-phase crossed-beam investigations, the reaction is described in terms of both an abstraction process (major) and a short-lived addition dynamic complex (minor).

Crossed-beam radical-radical reaction dynamics of O(3P)+C3H3-->H(2S)+C3H2O

The radical-radical oxidation reaction, O(3P)+C3H3 (propargyl)-->H(2S)+C3H2O (propynal), was investigated using vacuum-ultraviolet laser-induced fluorescence spectroscopy in a crossed-beam configuration, together with ab initio and statistical calculations. The barrierless addition of O(3P) to C3H3 is calculated to form energy-rich addition complexes on the lowest doublet potential energy surface, which subsequently undergo direct decomposition steps leading to the major reaction products, H+C3H(2)O (propynal). According to the nascent H-atom Doppler-profile analysis, the average translational energy of the products and the fraction of the average transitional energy to the total available energy were determined to be 5.09+/-0.36 kcal/mol and 0.077, respectively. On the basis of a comparison with statistical prior calculations, the reaction mechanism and the significant internal excitation of the polyatomic propynal product can be rationalized in terms of the formation of highly activated, short-lived addition-complex intermediates and the adiabaticity of the excess available energy along the reaction coordinate.

A study of the radical-radical reaction dynamics of O(P3)+t-C4H9→OH+iso-C4H8

The radical-radical reaction dynamics of ground-state atomic oxygen [O(3P)] with t-butyl radicals (t-C4H9) in the gas phase were investigated using high-resolution laser spectroscopy in a crossed-beam configuration, together with ab initio theoretical calculations. The radical reactants, O(3P) and t-C4H9, were produced by the photodissociation of NO2 and the supersonic flash pyrolysis of the precursor, azo-t-butane, respectively. A new exothermic channel, O(3P)+t-C4H9 --> OH+iso-C4H8, was identified and the nascent rovibrational distributions of the OH (X 2Pi: upsilon" = 0,1,2) products were examined. The population analyses for the two spin-orbit states of F1(2Pi32) and F2(2Pi12) showed that the upsilon" = 0 level is described by a bimodal feature composed of low- and high-N" rotational components, whereas the upsilon" = 1 and 2 levels exhibit unimodal distributions. No noticeable spin-orbit or Lambda-doublet propensities were observed in any vibrational state. The partitioning ratio of the vibrational populations (Pupsilon") with respect to the low-N" components of the upsilon" = 0 level was estimated to be P0:P1:P2 = 1:1.17+/-0.24:1.40+/-0.11, indicating that the nascent internal distributions are highly excited. On the basis of the comparison of the experimental results with the statistical theory, the reaction mechanism at the molecular level can be described in terms of two competing dynamic pathways: the major, direct abstraction process leading to the inversion of the vibrational populations, and the minor, short-lived addition-complex process responsible for the hot rotational distributions. After considering the reaction exothermicity, the barrier height, and the number of intermediates along the addition reaction pathways on the lowest doublet potential energy surface, the formation of CH3COCH3(acetone)+CH3 was predicted to be dominant in the addition mechanism.

Ab initio investigations of the radical-radical reaction of O(3P) + C3H3

We present ab initio calculations of the reaction of ground-state atomic oxygen [O((3)P)] with a propargyl (C(3)H(3)) radical based on the application of the density-functional method and the complete basis-set model. It has been predicted that the barrierless addition of O((3)P) to C(3)H(3) on the lowest doublet potential-energy surface produces several energy-rich intermediates, which undergo subsequent isomerization and decomposition steps to generate various exothermic reaction products: C(2)H(3)+CO, C(3)H(2)O+H, C(3)H(2)+OH, C(2)H(2)+CHO, C(2)H(2)O+CH, C(2)HO+CH(2), and CH(2)O+C(2)H. The respective reaction pathways are examined extensively with the aid of statistical Rice-Ramsperger-Kassel-Marcus calculations, suggesting that the primary reaction channel is the formation of propynal (CHCCHO)+H. For the minor C(3)H(2)+OH channel, which has been reported in recent gas-phase crossed-beam experiments [H. Lee et al., J. Chem. Phys. 119, 9337 (2003); 120, 2215 (2004)], a comparison on the basis of prior statistical calculations is made with the nascent rotational state distributions of the OH products to elucidate the mechanistic and dynamic characteristics at the molecular level.

A combined crossed-beam and theoretical investigation of radical-radical reaction dynamics of O(P3)+t-C4H9→OH+iso-C4H8

The oxidation reaction dynamics of a saturated hydrocarbon radical t-butyl leading to the isobutene +OH (X 2Pi:v"=0, 1, 2) products in the gas phase were first investigated by applying a combination of high-resolution spectroscopy in a crossed-beam configuration and ab initio calculations. By comparing the nascent OH populations with the statistical theory, the reaction mechanism at the molecular level can be described in terms of two competing dynamic pathways: the major direct abstraction process leading to the inversion of vibrational populations, and the minor short-lived addition-complex process for hot rotational distributions.

Exploring the dynamics of hydrogen atom release from the radical–radical reaction of O(3P) with C3H5

The gas-phase radical-radical reaction dynamics of O(3P) + C3H5 --> H(2S) + C3H4O was studied at an average collision energy of 6.4 kcal/mol in a crossed beam configuration. The ground-state atomic oxygen [O(3P)] and allyl radicals (C3H5) were generated by the photolysis of NO2 and the supersonic flash pyrolysis of allyl iodide, respectively. Nascent hydrogen atom products were probed by the vacuum-ultraviolet-laser induced fluorescence spectroscopy in the Lyman-alpha region centered at 121.6 nm. With the aid of the CBS-QB3 level of ab initio theory, it has been found that the barrierless addition of O(3P) to C3H5 forms the energy-rich addition complexes on the lowest doublet potential energy surface, which are predicted to undergo a subsequent direct decomposition step leading to the reaction products H + C3H4O. The major counterpart C3H4O of the probed hydrogen atom is calculated to be acrolein after taking into account the factors of barrier height, reaction enthalpy, and the number of intermediates involved along the reaction pathway. The nascent H-atom Doppler profile analysis shows that the average center-of-mass translational energy of the H + C3H4O products and the fraction of the total available energy released as the translational energy were determined to be 3.83 kcal/mol and 0.054, respectively. On the basis of comparison with statistical calculations, the reaction proceeds through the formation of short-lived addition complexes rather than statistical, long-lived intermediates, and the polyatomic acrolein product is significantly internally excited at the moment of the decomposition.

A combined crossed beam and theoretical investigation of O(3P)+C3H3→C3H2+OH

The radical-radical reaction dynamics of ground-state atomic oxygen [O(3P)] with propargyl radicals (C3H3) has first been investigated in a crossed beam configuration. The radical reactants O(3P) and C3H3 were produced by the photodissociation of NO2 and the supersonic flash pyrolysis of precursor propargyl bromide, respectively. A new exothermic channel of O(3P) + C3H3 --> C3H2 + OH was identified and the nascent distributions of the product OH in the ground vibrational state (X 2Pi:nu" = 0) showed bimodal rotational excitations composed of the low- and high-N" components without spin-orbit propensities. The averaged ratios of Pi(A')/Pi(A") were determined to be 0.60 +/- 0.28. With the aid of ab initio theory it is predicted that on the lowest doublet potential energy surface, the reaction proceeds via the addition complexes formed through the barrierless addition of O(3P) to C3H3. The common direct abstraction pathway through a collinear geometry does not occur due to the high entrance barrier in our low collision energy regime. In addition, the major reaction channel is calculated to be the formation of propynal (CHCCHO) + H, and the counterpart C3H2 of the probed OH product in the title reaction is cyclopropenylidene (1c-C3H2) after considering the factors of barrier height, reaction enthalpy and structural features of the intermediates formed along the reaction coordinate. On the basis of the statistical prior and rotational surprisal analyses, the ratio of population partitioning for the low- and high-N" is found to be about 1:2, and the reaction is described in terms of two competing addition-complex mechanisms: a major short-lived dynamic complex and a minor long-lived statistical complex. The observed unusual reaction mechanism stands in sharp contrast with the reaction of O(3P) with allyl radical (C3H5), a second significant conjugated hydrocarbon radical, which shows totally dynamic processes [J. Chem. Phys. 117, 2017 (2002)], and should be understood based upon the characteristic electronic structures and reactivity of the intermediates on the potential energy surface.