Comparisons of classical chemical dynamics simulations of the unimolecular decomposition of classical and quantum microcanonical ensembles

Previous studies have shown that classical trajectory simulations often give accurate results for short-time intramolecular and unimolecular dynamics, particularly for initial non-random energy distributions. To obtain such agreement between experiment and simulation, the appropriate distributions must be sampled to choose initial coordinates and momenta for the ensemble of trajectories. If a molecule's classical phase space is sampled randomly, its initial decomposition will give the classical anharmonic microcanonical (RRKM) unimolecular rate constant for its decomposition. For the work presented here, classical trajectory simulations of the unimolecular decomposition of quantum and classical microcanonical ensembles, at the same fixed total energy, are compared. In contrast to the classical microcanonical ensemble, the quantum microcanonical ensemble does not sample the phase space randomly. The simulations were performed for CH(4), C(2)H(5), and Cl(-)---CH(3)Br using both analytic potential energy surfaces and direct dynamics methods. Previous studies identified intrinsic RRKM dynamics for CH(4) and C(2)H(5), but intrinsic non-RRKM dynamics for Cl(-)---CH(3)Br. Rate constants calculated from trajectories obtained by the time propagation of the classical and quantum microcanonical ensembles are compared with the corresponding harmonic RRKM estimates to obtain anharmonic corrections to the RRKM rate constants. The relevance and accuracy of the classical trajectory simulation of the quantum microcanonical ensemble, for obtaining the quantum anharmonic RRKM rate constant, is discussed.     

Non-RRKM dynamics in the CH3O2 + NO reaction system

Quasi-classical trajectory (QCT) calculations on a model potential energy surface (PES) show strong deviations from statistical Rice-Ramsperger-Kassel-Marcus (RRKM) rate theory for the decomposition reaction (1) CH3OONO* --> CH3O + NO2, where the highly excited CH3OONO* was formed by (2) CH3O2 + NO --> CH3OONO*. The model PES accurately describes the vibrational frequencies, structures, and thermochemistry of the cis- and trans-CH3OONO isomers; it includes cis-trans isomerization in addition to reactions 1 and 2 but does not include nitrate formation, which is too slow to affect the decay rate of CH3OONO*. The QCT results give a strongly time-dependent rate constant for decomposition and damped oscillations in the decomposition rate, not predicted by statistical rate theory. Anharmonicity is shown to play an important role in reducing the rate constant by a factor of 10 smaller than predicted using classical harmonic RRKM theory (microcanonical variational transition state theory). Master equation simulations of organic nitrate yields published previously by two groups assumed that RRKM theory is accurate for reactions 1 and 2 but required surprising parametrizations to fit experimental nitrate yield data. In the present work, it is hypothesized that the non-RRKM rate of reaction (1) and vibrational anharmonicity are at least partly responsible for the surprising parameters.     

Trajectory studies of S(N)2 nucleophilic substitution. 8. Central barrier dynamics for gas phase Cl(-) + CH(3)Cl

Quasiclassical direct dynamics trajectories, calculated at the MP2/6-31G level of theory, are used to study the central barrier dynamics for the C1(-) + CH(3)Cl S(N)2 reaction. Extensive recrossings of the central barrier are observed in the trajectories. The dynamics of the Cl(-)-CH(3)Cl complex is non-RRKM and transition state theory (TST) is predicted to be an inaccurate model for calculating the Cl(-) + CH(3)Cl S(N)2 rate constant. Direct dynamics trajectories also show that Cl(-) + CH(3)Cl trajectories, which collide backside along the S(N)2 reaction path, do not form the Cl(-)-CH(3)Cl complex. This arises from weak coupling between the Cl(-)-CH(3)Cl intermolecular and CH(3)Cl intramolecular modes. The trajectory results are very similar to those of a previous trajectory study, based on a HF/6-31G* analytic potential energy function, which gives a less accurate representation of the central barrier region of the Cl(-) + CH(3)Cl reaction than does the MP2/6-31G* level of theory used here. Experiments are suggested for investigating the non-RRKM and non-TST dynamics predicted by the trajectories.     

Insight into selected reactions in low-temperature dimethyl ether combustion from Born-Oppenheimer molecular dynamics

Dimethyl ether is under consideration as an alternative diesel fuel. Its combustion chemistry is as yet ill-characterized. Here we use Born-Oppenheimer molecular dynamics (BOMD) based on DFT-B3LYP forces to investigate the short-time dynamics of selected features of the low-temperature dimethyl ether (DME) oxidation potential energy surface. Along the chain propagation pathway, we run BOMD simulations from the transition state involving the decomposition of (*)CH(2)OCH(2)OOH to two CH(2)=O and an (*)OH radical. We predict that formaldehyde C-O stretch overtones are excited, consistent with laser photolysis experiments. We also predict that O-H overtones are excited for the (*)OH formed from (*)CH(2)OCH(2)OOH dissociation. We also investigate short-time dynamics involved in chain branching. First, we examine the isomerization transition state of (*)OOCH(2)OCH(2)OOH --> HOOCH(2)OCHOOH. The latter species is predicted to be a short-lived metastable radical that decomposes within 500 fs to hydroperoxymethyl formate (HPMF; HOOCH(2)OC(=O)H) and the first (*)OH of chain branching. The dissociation of HOOCH(2)OCHOOH exhibits non-RRKM behavior in its lifetime profile, which may be due to conformational constraints or slow intramolecular vibrational energy transfer (IVR) from the nascent H-O bond to the opposite end of the radical, where O-O scission occurs to form HPMF and (*)OH. In a few trajectories, we see HOOCH(2)OCHOOH recross back to (*)OOCH(2)OCH(2)OOH because the isomerization is endothermic, with only an 8 kcal/mol barrier to recrossing. Therefore, some inhibition of chain-branching may be due to recrossing. Second, trajectories run from the transition state leading to the direct decomposition of HPMF (an important source of the second (*)OH radical in chain branching) to HCO, (*)OH, and HC(=O)OH show that these products can recombine to form many other possible products. These products include CH(2)OO + HC(=O)OH, H(2)O + CO + HC(=O)OH, HC(=O)OH + HC(=O)OH, and HC(=O)C(=O)H + H(2)O, which (save CH(2)OO + HC(=O)OH) are all more thermodynamically stable than the original HCO + (*)OH + HC(=O)OH products. Moreover, the multitude of extra products suggest that standard statistical rate theories cannot completely describe the reaction kinetics of significantly oxygenated compounds such as HPMF. These secondary products consume the second (*)OH required for explosive combustion, suggesting an inhibition of DME fuel combustion is likely.     

Ab initio energies and product branching ratios for the O+C3H6 reaction

Intermediate and transition-state energies have been calculated for the O+C3H6 (propene) reaction using the compound ab initio CBS-QB3 and G3 methods in combination with density functional theory. The lowest-lying triplet and singlet potential energy surfaces of the O-C3H6 system were investigated. RRKM statistical theory was used to predict product branching fractions over the 300-3000 K temperature and 0.001-760 Torr pressure ranges. The oxygen atom adds to the C3H6 terminal olefinic carbon in the primary step to generate a nascent triplet biradical, CH3CHCH2O. On the triplet surface, unimolecular dissociation of CH3CHCH2O to yield H+CH3CHCHO is favored over the entire temperature range, although the competing H2CO+CH3CH product channel becomes significant at high temperature. Rearrangement of triplet CH3CHCH2O to CH3CH2CHO (propanal) via a 1,2 H-atom shift has a barrier of 122.3 kJ mol(-1), largely blocking this reaction channel and any subsequent dissociation products. Intersystem crossing of triplet CH3CHCH2O to the singlet surface, however, leads to facile rearrangement to singlet CH3CH2CHO, which dissociates via numerous product channels. Pressure was found to have little influence over the branching ratios under most conditions, suggesting that the vibrational self-relaxation rates for p相似文献   

Dissociation of the ground state vinoxy radical and its photolytic precursor chloroacetaldehyde: electronic nonadiabaticity and the suppression of the H+ketene channel

This work is a study of the competition between the two unimolecular reaction channels available to the vinoxy radical (CH(2)CHO), C-H fission to form H+ketene, and isomerization to the acetyl radical (CH(3)CO) followed by C-C fission to form CH(3) + CO. Chloroacetaldehyde (CH(2)ClCHO) was used as a photolytic precursor to the vinoxy radical in its ground state; photodissociation of chloroacetaldehyde at 193 nm produces vinoxy radicals with internal energies spanning the G3//B3LYP calculated barriers to the two available unimolecular reaction channels. The onset of the CH(3) + CO channel, via isomerization to the acetyl radical, was found to occur at an internal energy of 41 +/- 2 kcal/mol, agreeing well with our calculated isomerization barrier of 40.8 kcal/mol. Branching to the H+ketene channel was too small to be detected; we conclude that the branching to the H+ketene channel must be at least a factor of 200 lower than what is predicted by a RRKM analysis based on our electronic structure calculations. This dramatic result may be explained in part by the presence of a conical intersection at planar geometries along the reaction coordinate leading to H+ketene, which results in electronically nonadiabatic recrossing of the transition state.     

Ab initio molecular dynamics study on the electron capture processes of protonated methane (CH5+)

Electron capture dynamics of protonated methane (CH5(+)) have been investigated by means of a direct ab initio molecular dynamics (MD) method. First, the ground and two low-lying state structures of CH5 (+) with eclipsed Cs , staggered Cs and C2v symmetries were examined as initial geometries in the dynamics calculation. Next, the initial structures of CH5 (+) in the Franck-Condon (FC) region were generated by inclusion of zero point energy and then trajectories were run from the selected points on the assumption of vertical electron capture. Two competing reaction channels were observed: CH5 (+) + e (-)--> CH4 + H (I) and CH5 (+) + e (-) --> CH3 + H2 (II). Channel II occurred only from structures very close to the s- Cs geometry for which two protons with longer C-H distances are electronically equivalent in CH5 (+). These protons have the highest spin density as hydrogen atoms following vertical electron capture of CH5 (+) and are lost as H2. On the other hand, channel I was formed from a wide structural region of CH5 (+). The mechanism of the electron capture dynamics of CH5 is discussed on the basis of the theoretical results.     

A direct dynamics trajectory study of F- + CH(3)OOH reactive collisions reveals a major non-IRC reaction path

A direct dynamics simulation at the B3LYP/6-311+G(d,p) level of theory was used to study the F- + CH3OOH reaction dynamics. The simulations are in excellent agreement with a previous experimental study (J. Am. Chem. Soc. 2002, 124, 3196). Two product channels, HF + CH2O + OH- and HF + CH3OO-, are observed. The former dominates and occurs via an ECO2 mechanism in which F- attacks the CH3- group, abstracting a proton. Concertedly, a carbon-oxygen double bond is formed and OH- is eliminated. Somewhat surprisingly this is not the reaction path, predicted by the intrinsic reaction coordinate (IRC), which leads to a deep potential energy minimum for the CH2(OH)2...F- complex followed by dissociation to HF + CH2(OH)O-. None of the direct dynamics trajectories followed this path, which has an energy release of -63 kcal/mol and is considerably more exothermic than the ECO2 path whose energy release is -27 kcal/mol. Other product channels not observed, and which have a lower energy than that for the ECO2 path, are F- + CO + H2 + H2O (-43 kcal/mol), F- + CH2O + H2O (-51 kcal/mol), and F- + CH2(OH)2 (-60 kcal/mol). Formation of the CH3OOH...F- complex, with randomization of its internal energy, is important, and this complex dissociates via the ECO2 mechanism. Trajectories which form HF + CH3OO- are nonstatistical events and, for the 4 ps direct dynamics simulation, are not mediated by the CH3OOH...F- complex. Dissociation of this complex to form HF + CH3OO- may occur on longer time scales.     

烯丙基自由基(·C3H5)与氧气(O2)反应机理的理论计算

用量子化学计算方法对CH3CH=·CH与O2气的反应机理进行了理论研究, 在B3LYP/6-311G(d,p) 水平下优化稳定分子结构和寻找过渡态, 并在此构型的基础上, 采用CCSD(T)/6-311G(d,p)方法得到各驻点的高级单点能量. 找到主要路径R(CH3CH=·CH+O2)→m1(trans-CH3CH=CHOO)→m2(形成COO三元环)→m3(C—C键断裂,同时生成CH3CH—O—CHO)→P2(C—O键断裂生成CH3CHO+CHO); 并与C2H3等共轭体系进行了对比.     

Quantifying the non-RRKM effect in the H + O2 ⇄ OH + O reaction

In the present investigation the non-RRKM behavior in the title reaction is quantified in two different ways: (1) Quasiclassical trajectory calculations of the thermal rate coefficient are compared with results from a microcanonical variational transition-state theory/RRKM model. Results on both the Varandas DMBE IV and Melius-Blint potentials indicate that the non-RRKM behavior acts to reduce the thermal rate coefficient by about a factor of two, independent of temperature from 250 K to 5500 K. The QCT thermal rate coefficients on the two potentials are in remarkably good agreement with experiment and with each other over the entire temperature range. (2) The non-RRKM behavior as a classical phenomenon is demonstrated and quantified on both potentials by a direct test of the fundamental assumption. Complex-forming classical trajectories, started as either O + OH or H + O₂, are shown preferentially to return to the region of configuration space from which they were started. This test is discussed in detail in the text. The transition of the non-RRKM behavior from classical to quantum mechanics is also discussed. © 1997 John Wiley & Sons, Inc. Int J Chem Kinet 29: 275–287, 1997.     

E−Z Isomerization in Guanidine: Second-order Saddle Dynamics,Non-statisticality,and Time-frequency Analysis

Our recent work on the E−Z isomerization reaction of guanidine using ab initio chemical dynamics simulations [Rashmi et al., Regul. Chaotic Dyn. 2021 , 26, 119] emphasized the role of second-order saddle ( SOS ) in the isomerization reaction; however, we could not unequivocally establish the non-statistical nature of the dynamics followed in the reaction. In the present study, we performed thousands of on-the-fly trajectories using forces computed at the MNDO level to investigate the influence of second-order saddle in the E−Z isomerization reaction of guanidine and the role of intramolecular vibrational energy redistribution (IVR) on the reaction dynamics. The simulations reveal that while majority of the trajectories follow the traditional transition state pathways, 15 % of the trajectories follow the SOS path. The dynamics was found to be highly non-statistical with the survival probabilities of the reactants showing large deviations from those obtained within the RRKM assumptions. In addition, a detailed analysis of the dynamics using time-dependent frequencies and the frequency ratio spaces reveal the existence of multiple resonance junctions that indicate the existence of regular dynamics and long-lived quasi-periodic trajectories in the phase space associated with non-RRKM behavior.     

Implementation of a microcanonical variational transition state theory for direct dynamics calculations of rate constants

Calculation of microcanonical rate constants has been an important field in chemical dy-namic studies for many years because it can be used not only to give good prediction of rate con-stants in microcanonical assembly, but also to calculate rate constants with certain conserved quantum numbers such as the total angular momentum, and in turn, can be easily converted into thermal rate constants[1—3]. The widely used method for calculating microcanonical rate constants of unimolecular reac-tions…     

Experimental and theoretical studies of the reactions Y (a2D) + H2CO and Y (a2D) + CH3CHO

The reactions of ground state Y (a(2)D) with H(2)CO and CH(3)CHO were studied at a range of collision energies in crossed molecular beams. For reaction with H(2)CO, three product channels were observed: formation of YH(2) + CO, YCO + H(2), and YHCO + H. Reaction with CH(3)CHO led to three analogous product channels involving formation of HYCH(3) + CO, YCH(2)CO + H(2), and YCH(3)CO + H. The calculated CCSD(T) energetics and DFT geometries for key intermediates in both reactions, together with RRKM theory, are used to calculate a priori the branching ratios between various product channels. These calculated values are compared to those obtained experimentally.     

The products of the thermal decomposition of CH3CHO

We have used a heated 2 cm × 1 mm SiC microtubular (μtubular) reactor to decompose acetaldehyde: CH(3)CHO + Δ → products. Thermal decomposition is followed at pressures of 75-150 Torr and at temperatures up to 1675 K, conditions that correspond to residence times of roughly 50-100 μs in the μtubular reactor. The acetaldehyde decomposition products are identified by two independent techniques: vacuum ultraviolet photoionization mass spectroscopy (PIMS) and infrared (IR) absorption spectroscopy after isolation in a cryogenic matrix. Besides CH(3)CHO, we have studied three isotopologues, CH(3)CDO, CD(3)CHO, and CD(3)CDO. We have identified the thermal decomposition products CH(3) (PIMS), CO (IR, PIMS), H (PIMS), H(2) (PIMS), CH(2)CO (IR, PIMS), CH(2)=CHOH (IR, PIMS), H(2)O (IR, PIMS), and HC≡CH (IR, PIMS). Plausible evidence has been found to support the idea that there are at least three different thermal decomposition pathways for CH(3)CHO; namely, radical decomposition: CH(3)CHO + Δ → CH(3) + [HCO] → CH(3) + H + CO; elimination: CH(3)CHO + Δ → H(2) + CH(2)=C=O; isomerization∕elimination: CH(3)CHO + Δ → [CH(2)=CH-OH] → HC≡CH + H(2)O. An interesting result is that both PIMS and IR spectroscopy show compelling evidence for the participation of vinylidene, CH(2)=C:, as an intermediate in the decomposition of vinyl alcohol: CH(2)=CH-OH + Δ → [CH(2)=C:] + H(2)O → HC≡CH + H(2)O.     

Energy-flow dynamics in the molecular channel of propanal photodissociation, C2H5CHO --> C2H6 + CO: direct ab initio molecular dynamics study

Direct ab initio molecular dynamics calculations have been carried out for the molecular channel of the photodissociation of propanal, C2H5CHO --> C2H6 + CO, at the RMP2(full)/cc-pVDZ level of ab initio molecular orbital theory. The initial conditions were generated using the microcanonical sampling to put the excess energy randomly into all vibrational modes of the TS. Starting from the TS, a total of approximately 700 trajectories were numerically integrated for 100 fs. The obtained final energy distributions for the C2H6 and CO fragments and their relative translational motion were found to be quite similar to those obtained for the acetaldehyde reaction, CH3CHO --> CH4 + CO, in our previous study (Chem. Phys. Lett. 2006, 421, 549) despite the fact that the number of degree of freedom for C2H6 is larger than that for CH4. The coupling between the intrinsic reaction coordinate and one of the generalized normal modes orthogonal to it was predicted substantially strong around s = 1.4 amu(1/2) bohr, and it is expected that the energy flow out of C2H6 proceeds through this coupling. However, the obtained energy distributions strongly suggest that the coupling among the modes in C2H6 is quite small and the intramolecular energy redistribution does not occur efficiently in this molecule.     

Ab initio and direct dynamics studies of the reaction of singlet methylene with acetylene and the lifetime of the cyclopropene complex

The energetics of the (1)CH(2) + C(2)H(2) --> H + C(3)H(3) reaction are accurately calculated using an extrapolated coupled-cluster/complete basis set (CBS) method based on the cc-pVDZ, cc-pVTZ, and cc-pVQZ basis sets. The reaction enthalpy (0 K) is predicted to be -20.33 kcal/mol. This reaction has no classical barrier in either the entrance or exit channel. However, there are several stable intermediates-cyclopropene (c-C(3)H(4)), allene (CH(2)CCH(2)), and propyne (CH(3)CCH)-along the minimum energy path. These intermediates with zero-point energy corrections lie below the reactants by 87.11 (c-C(3)H(4)), 109.69 (CH(2)CCH(2)), and 110.78 kcal/mol (CH(3)CCH). The vibrationally adiabatic ground-state (VAG) barrier height for c-C(3)H(4) isomerization to allene is obtained as 45.2 kcal/mol, and to propyne as 37.2 kcal/mol. In addition, the (1)CH(2) + C(2)H(2) reaction is investigated utilizing the dual-level "scaling all correlation" (SAC) ab initio method of Truhlar et al., i.e., the UCCSD(SAC)/cc-pVDZ theory. Results show that the reaction occurs via long-lived complexes. The lifetime of the cyclopropene intermediate is obtained as 3.2 +/- 0.4 ps. It is found that the intermediate propyne can be formed directly from reactants through the insertion of (1)CH(2) into a C-H bond of C(2)H(2). However, compared to the major mechanism in which the propyne is produced through a ring-opening of the cyclopropene complex, this reaction pathway is much less favorable. Finally, the theoretical thermal rate constant exhibits a negative temperature dependence, which is in excellent agreement with the previous results. The temperature dependence is consistent with the earlier RRKM results but weaker than the experimental observations at high temperatures.     

Exploring the dynamics of reaction N((2)D)+C2H4 with crossed molecular-beam experiments and quantum-chemical calculations

We conducted the title reaction using a crossed molecular-beam apparatus, quantum-chemical calculations, and RRKM calculations. Synchrotron radiation from an undulator served to ionize selectively reaction products by advantage of negligibly small dissociative ionization. We observed two products with gross formula C(2)H(3)N and C(2)H(2)N associated with loss of one and two hydrogen atoms, respectively. Measurements of kinetic-energy distributions, angular distributions, low-resolution photoionization spectra, and branching ratios of the two products were carried out. Furthermore, we evaluated total branching ratios of various exit channels using RRKM calculations based on the potential-energy surface of reaction N((2)D)+C(2)H(4) established with the method CCSD(T)/6-311+G(3df,2p)//B3LYP/6-311G(d,p)+ZPE[B3LYP/6-311G(d,p)]. The combination of experimental and computational results allows us to reveal the reaction dynamics. The N((2)D) atom adds to the C=C π-bond of ethene (C(2)H(4)) to form a cyclic complex c-CH(2)(N)CH(2) that directly ejects a hydrogen atom or rearranges to other intermediates followed by elimination of a hydrogen atom to produce C(2)H(3)N; c-CH(2)(N)CH+H is the dominant product channel. Subsequently, most C(2)H(3)N radicals, notably c-CH(2)(N)CH, further decompose to CH(2)CN+H. This work provides results and explanations different from the previous work of Balucani et al. [J. Phys. Chem. A, 2000, 104, 5655], indicating that selective photoionization with synchrotron radiation as an ionization source is a good choice in chemical dynamics research.     

High-dimensional ab initio potential energy surfaces for reaction dynamics calculations

There has been great progress in the development of potential energy surfaces (PESs) for reaction dynamics that are fits to ab initio energies. The fitting techniques described here explicitly represent the invariance of the PES with respect to all permutations of like atoms. A review of a subset of dynamics calculations using such PESs (currently 16 such PESs exist) is then given. Bimolecular reactions of current interest to the community, namely, H + CH(4) and F + CH(4), are focused on. Unimolecular reactions are then reviewed, with a focus on the photodissociation dynamics of H(2)CO and CH(3)CHO, where so-called "roaming" pathways have been discovered. The challenges for electronically non-adiabatic reactions, and associated PESs, are presented with a focus on the OH* + H(2) reaction. Finally, some thoughts on future directions and challenges are given.     

Trajectory surface hopping study of the O(3P) + ethylene reaction dynamics

Spin-orbit coupling (SOC) induced intersystem crossing (ISC) has long been believed to play a crucial role in determining the product distributions in the O(3P) + C2H4 reaction. In this paper, we present the first nonadiabatic dynamics study of the title reaction at two center-of-mass collision energies: 0.56 eV, which is barely above the H-atom abstraction barrier on the triplet surface, and 3.0 eV, which is in the hyperthermal regime. The calculations were performed using a quasiclassical trajectory surface hopping (TSH) method with the potential energy surface generated on the fly at the unrestricted B3LYP/6-31G(d,p) level of theory. To simplify our calculations, nonadiabatic transitions were only considered when the singlet surface intersects the triplet surface. At the crossing points, Landau-Zener transition probabilities were computed assuming a fixed spin-orbit coupling parameter, which was taken to be 70 cm-1 in most calculations. Comparison with a recent crossed molecular beam experiment at 0.56 eV collision energy shows qualitative agreement as to the primary product branching ratios, with the CH3 + CHO and H + CH2CHO channels accounting for over 70% of total product formation. However, our direct dynamics TSH calculations overestimate ISC so that the total triplet/singlet ratio is 25:75, compared to the observed 43:57. Smaller values of SOC reduce ISC, resulting in better agreement with the experimental product relative yields; we demonstrate that these smaller SOC values are close to being consistent with estimates based on CASSCF calculations. As the collision energy increases, ISC becomes much less important and at 3.0 eV, the triplet to singlet branching ratio is 71:29. As a result, the triplet products CH2 + CH2O, H + CH2CHO and OH + C2H3 dominate over the singlet products CH3 + CHO, H2 + CH2CO, etc.     

H atom attack on propene

The reaction of propene (CH(3)CH═CH(2)) with hydrogen atoms has been investigated in a heated single-pulsed shock tube at temperatures between 902 and 1200 K and pressures of 1.5-3.4 bar. Stable products from H atom addition and H abstraction have been identified and quantified by gas chromatography/flame ionization/mass spectrometry. The reaction for the H addition channel involving methyl displacement from propene has been determined relative to methyl displacement from 1,3,5-trimethylbenzene (135TMB), leading to a reaction rate, k(H + propene) → H(2)C═CH(2) + CH(3)) = 4.8 × 10(13) exp(-2081/T) cm(3)/(mol s). The rate constant for the abstraction of the allylic hydrogen atom is determined to be k(H + propene → CH(2)CH═CH(2) + H(2)) = 6.4 × 10(13) exp(-4168/T) cm(3)/(mol s). The reaction of H + propene has also been directly studied relative to the reaction of H + propyne, and the relationship is found to be log[k(H + propyne → acetylene + CH(3))/k(H + propene → ethylene + CH(3))] = (-0.461 ± 0.041)(1000/T) + (0.44 ± 0.04). The results showed that the rate constant for the methyl displacement reaction with propene is a factor of 1.05 ± 0.1 larger than that for propyne near 1000 K. The present results are compared with relevant earlier data on related compounds.