Atom-radical reaction dynamics of O(3P)+C3H5→C3H4+OH: Nascent rovibrational state distributions of product OH

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.

Vibrational predissociation of the phenol-water dimer: a view from the water

The vibrational predissociation (VP) dynamics of the phenol-water (PhOH-H₂O) dimer were studied by detecting H₂O fragments and using velocity map imaging (VMI) to infer the internal energy distributions of PhOH cofragments, pair-correlated with selected rotational levels of the H₂O fragments. Following infrared (IR) laser excitation of the hydrogen-bonded OH stretch fundamental of PhOH (Pathway 1) or the asymmetric OH stretch localized on H₂O (Pathway 2), dissociation to H₂O + PhOH was observed. H₂O fragments were monitored state-selectively by using 2+1 Resonance-Enhanced Multiphoton Ionization (REMPI) combined with time-of-flight mass spectrometry (TOF-MS). VMI of H₂O in selected rotational levels was used to derive center-of-mass (c.m.) translational energy (E_T) distributions. The pair-correlated internal energy distributions of the PhOH cofragments derived via Pathway 1 were well described by a statistical prior distribution. On the other hand, the corresponding distributions obtained via Pathway 2 show a propensity to populate higher-energy rovibrational levels of PhOH than expected from a statistical distribution and agree better with an energy-gap model. The REMPI spectra of the H₂O fragments from both pathways could be fit by Boltzmann plots truncated at the maximum allowed energy, with a higher temperature for Pathway 2 than that for Pathway 1. We conclude that the VP dynamics depends on the OH stretch level initially excited.

Exploring the chemical kinetics of partially oxidized intermediates by combining experiments, theory, and kinetic modeling

Partially oxidized intermediates play a central role in combustion and atmospheric chemistry. In this perspective, we focus on the chemical kinetics of alkoxy radicals, peroxy radicals, and Criegee intermediates, which are key species in both combustion and atmospheric environments. These reactive intermediates feature a broad spectrum of chemical diversity. Their reactivity is central to our understanding of how volatile organic compounds are degraded in the atmosphere and converted into secondary organic aerosol. Moreover, they sensitively determine ignition timing in internal combustion engines. The intention of this perspective article is to provide the reader with information about the general mechanisms of reactions initiated by addition of atomic and molecular oxygen to alkyl radicals and ozone to alkenes. We will focus on critical branching points in the subsequent reaction mechanisms and discuss them from a consistent point of view. As a first example of our integrated approach, we will show how experiment, theory, and kinetic modeling have been successfully combined in the first infrared detection of Criegee intermediates during the gas phase ozonolysis. As a second example, we will examine the ignition timing of n-heptane/air mixtures at low and intermediate temperatures. Here, we present a reduced, fuel size independent kinetic model of the complex chemistry initiated by peroxy radicals that has been successfully applied to simulate standard n-heptane combustion experiments.

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.

Gas-phase kinetics of the hydroxyl radical reaction with allene: absolute rate measurements at low temperature, product determinations, and calculations

The gas phase reaction of the hydroxyl radical with allene has been studied theoretically and experimentally in a continuous supersonic flow reactor over the range 50 ≤ T/K ≤ 224. This reaction has been found to exhibit a negative temperature dependence over the entire temperature range investigated, varying between (0.75 and 5.0) × 10(-11) cm(3) molecule(-1) s(-1). Product formation from the reaction of OH and OD radicals with allene (C(3)H(4)) has been investigated in a fast flow reactor through time-of-flight mass spectrometry, at pressures between 0.8 and 2.4 Torr. The branching ratios for adduct formation (C(3)H(4)OH) in this pressure range are found to be equal to 34 ± 16% and 48 ± 16% for the OH and OD + allene reactions, respectively, the only other channel being the formation of CH(3) or CH(2)D + H(2)CCO (ketene). Moreover, the rate constant for the OD + C(3)H(4) reaction is also found to be 1.4 times faster than the rate constant for the OH + C(3)H(4) reaction at 1.5 Torr and at 298 K. The experimental results and implications for atmospheric chemistry have been rationalized by quantum chemical and RRKM calculations.

Crossed beam study of the atom-radical reaction of ground state carbon atoms (C(3P)) with the vinyl radical (C2H3(X2A'))

The atom-radical reaction of ground state carbon atoms (C((3)P)) with the vinyl radical (C(2)H(3)(X(2)A')) was conducted under single collision conditions at a collision energy of 32.3 ± 2.9 kJ mol(-1). The reaction dynamics were found to involve a complex forming reaction mechanism, which is initiated by the barrier-less addition of atomic carbon to the carbon-carbon-double bond of the vinyl radical forming a cyclic C(3)H(3) radical intermediate. The latter has a lifetime of at least 1.5 times its rotational period and decomposes via a tight exit transition state located about 45 kJ mol(-1) above the separated products through atomic hydrogen loss to the cyclopropenylidene isomer (c-C(3)H(2)) as detected toward cold molecular clouds and in star forming regions.

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.

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.

Radical reaction C3H+NO: a mechanistic study

Although a number of hydrocarbon radicals including the heavier C(3)-radicals C(3)H(3) and C(3)H(5) have been experimentally shown to deplete NO effectively, no theoretical or experimental attempts have been made on the reactivity of the simplest C(3)-radical towards NO. In this article, we report our detailed mechanistic study on the C(3)H+NO reaction at the Gussian-3//B3LYP/6-31G(d) level by constructing the singlet and triplet electronic state [H,C(3),N,O] potential energy surfaces (PESs). The l-C(3)H+NO reaction is shown to barrierlessly form the entrance isomer HCCCNO followed by the direct O-elimination leading to HCCCN+(3)O on triplet PES, or by successive O-transfer, N-insertion, and CN bond-rupture to generate the product (1)HCCN+CO on singlet PES. The possible singlet-triplet intersystem crossings are also discussed. Thus, the novel reaction l-C(3)H+NO can proceed effectively even at low temperatures and is expected to play an important role in both combustion and interstellar processes. For the c-C(3)H+NO reaction, the initially formed H-cCCC-NO can most favorably isomerize to HCCCNO, and further evolution follows that of the l-C(3)H+NO reaction. Quantitatively, the c-C(3)H+NO reaction can take place barrierlessly on singlet PES, yet it faces a small barrier 2.7 kcal/mol on triplet PES. The results will enrich our understanding of the chemistry of the simplest C(3)-radical in both combustion and interstellar processes, which to date have received little attention despite their importance and available abundant studies on its structural and spectroscopic properties.

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.

Probing the barrier for CH2CHCO→CH2CH+CO by the velocity map imaging method

This work determines the dissociation barrier height for CH2CHCO --> CH2CH + CO using two-dimensional product velocity map imaging. The CH2CHCO radical is prepared under collision-free conditions from C-Cl bond fission in the photodissociation of acryloyl chloride at 235 nm. The nascent CH2CHCO radicals that do not dissociate to CH2CH + CO, about 73% of all the radicals produced, are detected using 157-nm photoionization. The Cl(2P(3/2)) and Cl(2P(1/2)) atomic fragments, momentum matched to both the stable and unstable radicals, are detected state selectively by resonance-enhanced multiphoton ionization at 235 nm. By comparing the total translational energy release distribution P(E(T)) derived from the measured recoil velocities of the Cl atoms with that derived from the momentum-matched radical cophotofragments which do not dissociate, the energy threshold at which the CH2CHCO radicals begin to dissociate is determined. Based on this energy threshold and conservation of energy, and using calculated C-Cl bond energies for the precursor to produce CH2CHC*O or C*H2CHCO, respectively, we have determined the forward dissociation barriers for the radical to dissociate to vinyl + CO. The experimentally determined barrier for CH2CHC*O --> CH2CH + CO is 21+/-2 kcal mol(-1), and the computed energy difference between the CH2CHC*O and the C*H2CHCO forms of the radical gives the corresponding barrier for C*H2CHCO --> CH2CH + CO to be 23+/-2 kcal mol(-1). This experimental determination is compared with predictions from electronic structure methods, including coupled-cluster, density-functional, and composite Gaussian-3-based methods. The comparison shows that density-functional theory predicts too low an energy for the C*H2CHCO radical, and thus too high a barrier energy, whereas both the Gaussian-3 and the coupled-cluster methods yield predictions in good agreement with experiment. The experiment also shows that acryloyl chloride can be used as a photolytic precursor at 235 nm of thermodynamically stable CH2CHC*O radicals, most with an internal energy distribution ranging from approximately 3 to approximately 21 kcal mol(-1). We discuss the results with respect to the prior work on the O(3P) + propargyl reaction and the analogous O(3P) + allyl system.

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.