Thermochemistry and accurate quantum reaction rate calculations for H2/HD/D2 + CH3

Accurate quantum-mechanical results for thermodynamic data, cumulative reaction probabilities (for J = 0), thermal rate constants, and kinetic isotope effects for the three isotopic reactions H2 + CH3 --> CH4 + H, HD + CH3 --> CH4 + D, and D2 + CH3 --> CH(3)D + D are presented. The calculations are performed using flux correlation functions and the multiconfigurational time-dependent Hartree (MCTDH) method to propagate wave packets employing a Shephard interpolated potential energy surface based on high-level ab initio calculations. The calculated exothermicity for the H2 + CH3 --> CH4 + H reaction agrees to within 0.2 kcal/mol with experimentally deduced values. For the H2 + CH3 --> CH4 + H and D2 + CH3 --> CH(3)D + D reactions, experimental rate constants from several groups are available. In comparing to these, we typically find agreement to within a factor of 2 or better. The kinetic isotope effect for the rate of the H2 + CH3 --> CH4 + H reaction compared to those for the HD + CH3 --> CH4 + D and D2 + CH3 --> CH(3)D + D reactions agree with experimental results to within 25% for all data points. Transition state theory is found to predict the kinetic isotope effect accurately when the mass of the transferred atom is unchanged. On the other hand, if the mass of the transferred atom differs between the isotopic reactions, transition state theory fails in the low-temperature regime (T < 400 K), due to the neglect of the tunneling effect.     

Kinetics of the thiazolium ion-catalyzed benzoin condensation.

The formation of benzoin (Ph-CHOH-CO-Ph) from two molecules of benzaldehyde, catalyzed by 3-benzyl-5-(2-hydroxyethyl)-4-methylthiazolium bromide in methanol buffered with Et(3)N/Et(3)NH(+)Cl(-) has been studied. Initial-rate studies at various concentrations of PhCHO (0.1-1.7 M) showed that the reaction is close to being first order in PhCHO. Following the reaction in deuteriomethanol, (1)H NMR spectroscopy allowed rate constants for all three kinetically significant steps to be determined. These show that all three steps are partially rate-determining. A normal deuterium kinetic isotope effect for the overall reaction (k(H)/k(D) approximately 3.4) is observed using PhCDO, and a large inverse solvent isotope effect (k(D)/k(H) approximately 5.9) is observed using deuteriomethanol, consistent with the kinetic scheme presented here.     

Theoretical dynamic studies on the reactions of CH3C(O)CH3-nCl(n) (n = 0-3) with the chlorine atom

The theoretical investigations were performed on the reaction mechanisms for the title reactions CH(3)C(O)CH(3) + Cl --> products (R1), CH(3)C(O)CH(2)Cl + Cl --> products (R2), CH(3)C(O)CHCl(2) + Cl --> products (R3), and CH(3)C(O)CCl(3) + Cl --> products (R4) by ab initio direct dynamics approach. Two different reaction channels have been found: abstract of the H atom from methyl (--CH(3)) group or chloromethyl (--CH(3-n)Cl(n)) group of chloroacetone and addition of a Cl atom to the carbon atom of the carbonyl group of chloroacetone followed by methyl or chloromethyl eliminations. Because of the higher potential energy barrier, the contribution of addition-elimination reaction pathway to the total rate constants is very small and thus this pathway is insignificant in atmospheric conditions. The rate constants for the H-abstraction reaction channels are evaluated by using canonical variational transition state theory incorporating with the small-curvature tunneling correction. Theoretical overall rate constants are in good agreement with the available experimental values and decrease in the order of k(1) > k(2) > k(3) > k(4). The results indicate that for halogenated acetones the substitution of halogen atom (F or Cl) leads to the decrease in the C--H bond reactivity and more decrease of reactivity is caused by F-substitution.     

A systematic computational study of the reactions of HO2 with RO2: the HO2 + C2H5O2 reaction

The reaction of HO2 with C2H5O2 has been studied using the density functional theory (B3LYP) and the coupled-cluster theory [CCSD(T)]. The reaction proceeds on the triplet potential energy surface via hydrogen abstraction to form ethyl hydroperoxide and oxygen. On the singlet potential energy surface, the addition-elimination mechanism is revealed. Variational transition state theory is used to calculate the temperature-dependent rate constants in the range 200-1000 K. At low temperatures (e.g., below 300 K), the reaction takes place predominantly on the triplet surface. The calculated low-temperature rate constants are in good agreement with the experimental data. As the temperature increases, the singlet reaction mechanism plays more and more important role, with the formation of OH radical predominantly. The isotope effect of the reaction (DO2 + C2D5O2 vs HO2 + C2H5O2) is negligible. In addition, the triplet abstraction energetic routes for the reactions of HO2 with 11 alkylperoxy radicals (CnHmO2) are studied. It is shown that the room-temperature rate constants have good linear correlation with the activation energies for the hydrogen abstraction.     

Theoretical study on the rate constants for the C2H5 + HBr --> C2H6 + Br reaction

The reaction C(2)H(5) + HBr --> C(2)H(6) + Br has been theoretically studied over the temperature range from 200 to 1400 K. The electronic structure information is calculated at the BHLYP/6-311+G(d,p) and QCISD/6-31+G(d) levels. With the aid of intrinsic reaction coordinate theory, the minimum energy paths (MEPs) are obtained at the both levels, and the energies along the MEP are further refined by performing the single-point calculations at the PMP4(SDTQ)/6-311+G(3df,2p)//BHLYP and QCISD(T)/6-311++G(2df,2pd)//QCISD levels. The calculated ICVT/SCT rate constants are in good agreement with available experimental values, and the calculate results further indicate that the C(2)H(5) + HBr reaction has negative temperature dependence at T < 850 K, but clearly shows positive temperature dependence at T > 850 K. The current work predicts that the kinetic isotope effect for the title reaction is inverse in the temperature range from 200 to 482 K, i.e., k(HBr)/k(DBr) < 1.     

Unimolecular rate constants for HX or DX elimination (X = F, Cl) from chemically activated CF3CH2CH2Cl, C2H5CH2Cl, and C2D5CH2Cl: threshold energies for HF and HCl elimination

Vibrationally activated CF(3)CH(2)CH(2)Cl molecules were prepared with 94 kcal mol(-1) of vibrational energy by the combination of CF(3)CH(2) and CH(2)Cl radicals and with 101 kcal mol(-1) of energy by the combination of CF(3) and CH(2)CH(2)Cl radicals at room temperature. The unimolecular rate constants for elimination of HCl from CF(3)CH(2)CH(2)Cl were 1.2 x 10(7) and 0.24 x 10(7) s(-1) with 101 and 94 kcal mol(-1), respectively. The product branching ratio, k(HCl)/k(HF), was 80 +/- 25. Activated CH(3)CH(2)CH(2)Cl and CD(3)CD(2)CH(2)Cl molecules with 90 kcal mol(-1) of energy were prepared by recombination of C(2)H(5) (or C(2)D(5)) radicals with CH(2)Cl radicals. The unimolecular rate constant for HCl elimination was 8.7 x 10(7) s(-1), and the kinetic isotope effect was 4.0. Unified transition-state models obtained from density-functional theory calculations, with treatment of torsions as hindered internal rotors for the molecules and the transition states, were employed in the calculation of the RRKM rate constants for CF(3)CH(2)CH(2)Cl and CH(3)CH(2)CH(2)Cl. Fitting the calculated rate constants from RRKM theory to the experimental values provided threshold energies, E(0), of 58 and 71 kcal mol(-1) for the elimination of HCl or HF, respectively, from CF(3)CH(2)CH(2)Cl and 54 kcal mol(-1) for HCl elimination from CH(3)CH(2)CH(2)Cl. Using the hindered-rotor model, threshold energies for HF elimination also were reassigned from previously published chemical activation data for CF(3)CH(2)CH(3,) CF(3)CH(2)CF(3), CH(3)CH(2)CH(2)F, CH(3)CHFCH(3), and CH(3)CF(2)CH(3). In an appendix, the method used to assign threshold energies was tested and verified using the combined thermal and chemical activation data for C(2)H(5)Cl, C(2)H(5)F, and CH(3)CF(3).     

Understanding the complex dissociation dynamics of energy selected dichloroethylene ions: neutral isomerization energies and heats of formation by imaging photoelectron-photoion coincidence

The dissociative photoionization of 1,1-C(2)H(2)Cl(2), (E)-1,2-C(2)H(2)Cl(2), and (Z)-1,2-C(2)H(2)Cl(2) has been investigated at high energy and mass resolution using the imaging photoelectron photoion coincidence instrument at the Swiss Light Source. The asymmetric Cl-atom loss ion time-of-flight distributions were fitted to obtain the dissociation rates in the 10(3) s(-1) < k < 10(7) s(-1) range as a function of the ion internal energy. The results, supported by ab initio calculations, show that all three ions dissociate to the same C(2v) symmetry ClC═CH(2)(+) product ion. The 0 K onset energies thus establish the relative heats of formation of the neutral isomers, that is, the isomerization energies. The experimental rate constants, k(E), as well as ab initio calculations indicate an early isomerization transition state and no overall reverse barrier to dissociation. The major high energy channels are the parallel HCl loss and the sequential ClC═CH(2)(+) → HCCH(+) + Cl process, the latter in competition with a ClC═CH(2)(+) → ClCCH(+) + H reaction. A parallel C(2)H(2)Cl(2)(+) → C(2)HCl(2)(+) + H channel also weakly asserts itself. The 0 K onset energy for the sequential Cl loss reaction suggests no barrier to the production of the most stable acetylene ion product; thus the sequential Cl-atom loss is preceded by a ClC═CH(2)(+) → HC(Cl)CH(+) reorganization step with a barrier lower than that of the second Cl-atom loss. The breakdown diagram corresponding to this sequential dissociation reveals the internal energy distribution of the first C(2)H(2)Cl(+) daughter ion, which is determined by the kinetic energy release in the first, Cl loss reaction at high excess energies. At low kinetic energy release, this distribution corresponds to the predicted two translational degrees of freedom, whereas at higher energies, the excess energy partitioning is characteristic of only one translational degree of freedom. New Δ(f)H(o)(298K) of 3.7, 2.5, and 0.2 ± 1.75 kJ mol(-1) are proposed for 1,1-C(2)H(2)Cl(2), (E)-1,2-C(2)H(2)Cl(2), and (Z)-1,2-C(2)H(2)Cl(2), respectively, and the proton affinity of ClCCH is found to be 708.6 ± 2.5 kJ mol(-1).     

The addition reaction between silylene and ethyne: further isotope studies, pressure dependence studies, and quantum chemical calculations

Time-resolved kinetic studies of the reaction of dideutero-silylene, SiD 2, generated by laser flash photolysis of phenylsilane-d 3, have been carried out to obtain rate constants for its bimolecular reaction with C 2H 2. The reaction was studied in the gas phase over the pressure range 1-100 Torr in SF 6 bath gas, at five temperatures in the range 297-600 K. The second-order rate constants obtained by extrapolation to the high-pressure limits at each temperature fitted the Arrhenius equation log( k (infinity)/cm (3) molecule (-1) s (-1)) = (-10.05 +/- 0.05) + (3.43 +/- 0.36 kJ mol (-1))/ RT ln 10. The rate constants were used to obtain a comprehensive set of isotope effects by comparison with earlier obtained rate constants for the reactions of SiH 2 with C 2H 2 and C 2D 2. Additionally, pressure-dependent rate constants for the reaction of SiH 2 with C 2H 2 in the presence of He (1-100 Torr) were obtained at 300, 399, and 613 K. Quantum chemical (ab initio) calculations of the SiC 2H 4 reaction system at the G3 level support the initial formation of silirene, which rapidly isomerizes to ethynylsilane as the major pathway. Reversible formation of vinylsilylene is also an important process. The calculations also indicate the involvement of several other intermediates, not previously suggested in the mechanism. RRKM calculations are in semiquantitative agreement with the pressure dependences and isotope effects suggested by the ab initio calculations, but residual discrepancies suggest the possible involvement of the minor reaction channel, SiH 2 + C 2H 2 --> Si( (3)P 1) + C 2H 4. The results are compared and contrasted with previous studies of this reaction system.     

Study of the carbon-13 and deuterium kinetic isotope effects in the Cl and OH reactions of CH4 and CH3Cl

Relative rate experiments have been carried out for three isotopologues of chloromethane and their reactions with Cl atoms and OH radicals. The OH and Cl reaction rates of CH2DCl and CHD2Cl were measured by long-path FTIR spectroscopy relative to CH3Cl at 298+/-2 K and 1013+/-10 hPa in purified air. The FTIR spectra were fitted using a nonlinear least squares spectral fitting method including measured high-resolution infrared spectra as references. The relative reaction rates defined by alpha=klight/kheavy were determined to be kOH+CH3Cl/kOH+CH2DCl=1.41+/-0.05, kOH+CH3Cl/kOH+CHD2Cl=2.03+/-0.05, kCl+CH3Cl/kCl+CH2DCl=1.42+/-0.04, and kCl+CH3Cl/kCl+CHD2Cl=2.27+/-0.04. The carbon-13 and deuterium kinetic isotope effects in the OH and Cl reactions of CH3Cl were investigated further using variational transition state theory, and the results were compared to similar calculations performed for the CH4+OH/Cl reaction systems. The calculations show that the order of magnitude difference for the carbon-13 kinetic isotope effect in the OH reaction of CH3Cl compared to CH4 reported by Gola et al. (Atmos. Chem. Phys. 2005, 5, 2395) can be explained by the lower barrier to internal rotation of the OH radical in the transition state of the CH4+OH reaction than in the CH3Cl+OH reaction. The deuterium kinetic isotope effects can be explained in terms of combined variational effects and tunneling.     

Quantum-instanton evaluation of the kinetic isotope effects

A general quantum-mechanical method for computing kinetic isotope effects is presented. The method is based on the quantum-instanton approximation for the rate constant and on the path-integral Metropolis-Monte Carlo evaluation of the Boltzmann operator matrix elements. It computes the kinetic isotope effect directly, using a thermodynamic integration with respect to the mass of the isotope, thus avoiding the more computationally expensive process of computing the individual rate constants. The method should be more accurate than variational transition-state theories or the semiclassical instanton method since it does not assume a single tunneling path and does not use a semiclassical approximation of the Boltzmann operator. While the general Monte Carlo implementation makes the method accessible to systems with a large number of atoms, we present numerical results for the Eckart barrier and for the collinear and full three-dimensional isotope variants of the hydrogen exchange reaction H + H2 --> H2 + H. In all seven test cases, for temperatures between 250 and 600 K, the error of the quantum instanton approximation for the kinetic isotope effects is less than approximately 10%.     

Isotope effect-mapping of the ionization of glucose demonstrates unusual charge sharing

Isotopic substitution is known to affect kinetic rate constants and equilibrium constants in chemistry. In this study, we have used tritium substitution and high pH to probe the glucose left harpoon over right harpoon glucose(-) + H(+) equilibrium. Passing partially ionized mixtures of [(3)H]- and [(14)C]glucose over anionic exchange resin has permitted the detection of subtle differences in pK(a). We have found that, at pH 11.7 in an anionic exchange system, [(3)H]glucose always elutes ahead of the [(14)C]glucose, and we report isotope effects of 1.051 +/- 0.0007, 1.012 +/- 1.0005, 1.014 +/- 0.0004, 1.024 +/- 0.0003, 1.014 +/- 0.0004, and 1.015 +/- 0.0014 for [1-(3)H]-, [2-(3)H]-, [3-(3)H]-, [4-(3)H]-, [5-(3)H]-, and [6,6-(3)H(2)]glucose, respectively, as compared to either [2-(14)C]-or [6-(14)C]glucose. The elevated isotope effects at H1 and H4 imply unusual charge sharing in anionic aqueous glucose. Base titration of (13)C-chemical shift changes demonstrates the dominance of pyranose forms at elevated pH, and ab initio isotope effect computations on gas-phase glucose anions are presented.     

Laser photolysis studies of the photochemical reactions of chromium(III) octaethylporphyrin and tetramesitylporphyrin complexes in toluene solution

A nanosecond laser photolysis study was carried out for the Cr(III) porphyrin complexes of 2,3,7,8,12,13,17,18-octaethylporphyrin, [Cr(OEP)(Cl)(L)], and of 5,10,15,20-tetramesitylporphyrin, [Cr(TMP)(Cl)(L)], in toluene containing water and an excess amount of L (L = axial ligand). The laser photolysis generates the triplet excited state of the parent complex as well as a five-coordinate complex, [Cr(porphyrin)(Cl)], produced by the photodissociation of the axial ligand L. The yields for the formation of the triplet state and the photodissociation of L are found to markedly depend on the nature of both L and porphyrin ligand. The five-coordinate [Cr(porphyrin)(Cl)] readily reacts with both H(2)O and L in the bulk solution to give [Cr(porphyrin)(Cl)(H(2)O)] and [Cr(porphyrin)(Cl)(L)], respectively. The axial H(2)O ligand in [Cr(porphyrin)(Cl)(H(2)O)] is then substituted by the ligand L to regenerate the original complex [Cr(porphyrin)(Cl)(L)]. In principle, the substitution reaction takes place by the dissociative mechanism: the first step is the dissociation of H(2)O from [Cr(porphyrin)(Cl)(H(2)O)], followed by the reaction of the five-coordinate [Cr(porphyrin)(Cl)] with the ligand L to regenerate [Cr(porphyrin)(Cl)(L)]. The rate constants for this ligand substitution reaction are found to exhibit bell-shaped ligand concentration dependence. The detailed kinetic analysis revealed that both ligands L and H(2)O in toluene make a hydrogen bond with the axial H(2)O ligand in [Cr(porphyrin)(Cl)(H(2)O)] to yield dead-end complexes for the substitution reaction. The reaction mechanisms are discussed on the basis of the substituent effects of the porphyrin peripheral groups and the kinetic parameters determined from the temperature dependence of the rate constants.     

Kinetics and mechanism of hydrogen-atom abstraction from rhodium hydrides by alkyl radicals in aqueous solutions

The kinetics of the reaction of benzyl radicals with [L(1)(H(2)O)RhH{D}](2+) (L(1)=1,4,8,11-tetraazacyclotetradecane) were studied directly by laser-flash photolysis. The rate constants for the two isotopologues, k=(9.3±0.6) × 10(7) M(-1) s(-1) (H) and (6.2±0.3) × 10(7) M(-1) s(-1) (D), lead to a kinetic isotope effect k(H)/k(D)=1.5±0.1. The same value was obtained from the relative yields of PhCH(3) and PhCH(2)D in a reaction of benzyl radicals with a mixture of rhodium hydride and deuteride. Similarly, the reaction of methyl radicals with {[L(1)(H(2)O)RhH](2+) + [L(1)(H(2)O)RhD](2+)} produced a mixture of CH(4) and CH(3)D that yielded k(H)/k(D)=1.42±0.07. The observed small normal isotope effects in both reactions are consistent with reduced sensitivity to isotopic substitution in very fast hydrogen-atom abstraction reactions. These data disprove a literature report claiming much slower kinetics and an inverse kinetic isotope effect for the reaction of methyl radicals with hydrides of L(1)Rh.     

Multicoefficient Gaussian-3 calculation of the rate constant for the OH + CH4 reaction and its 12C/13C kinetic isotope effect with emphasis on the effects of coordinate system and torsional treatment

Rate constants and (12)C/(13)C kinetic isotope effects are calculated by direct dynamics for the OH + CH(4) --> H(2)O + CH(3) reaction. The electronic structure calculations required to generate the implicit potential energy surface were carried out by the high-level multicoefficient Gaussian-3/version-3 (MCG3) method and compared to two other multilevel methods, MC3BB and MC3MPW, and three density functional methods, M06-2X, BB1K, and MPW1K. The rate constants and (12)C/(13)C kinetic isotope effects are shown to depend strongly on the coordinate system used to calculate the frequencies as well as on the method used to account for the torsional anharmonicity of the lowest-frequency vibrational mode of the generalized transition states.     

A computational study of the reaction of ground-state nitrogen atoms with chloromethyl radicals

A computational study of the N(4S) + CH2Cl reaction has been carried out. The first step of the reaction is the formation of an initial intermediate (NCH2Cl), which is relatively stable and does not involve any energy barrier. The two most exothermic products are those resulting from the release of a chlorine atom, H2C=N + Cl and trans-HC=NH + Cl. A kinetic study within the framework of the statistical theories suggests that the kinetically preferred product is also the most exothermic one. This is in contrast with the analogue reaction of nitrogen atoms with CH2F, where the preferred product from both thermodynamic and kinetic points of view is HFCN + H. Therefore, reactions of nitrogen atoms with chloromethyl radicals release chlorine atoms as major products. The rate coefficient for the title reaction is estimated to be about 3.09 x 10(-13) cm3 s(-1) molecule(-1) at 300 K, a value four times smaller than the rate coefficient for its fluorine analogue.     

Quantum chemical and kinetic study of formation of 2-chlorophenoxy radical from 2-chlorophenol: unimolecular decomposition and bimolecular reactions with H, OH, Cl, and O2

This study investigates the kinetic parameters of the formation of the chlorophenoxy radical from the 2-chlorophenol molecule, a key precursor to polychlorinated dibenzo-p-dioxins and dibenzofurans (PCCD/F), in unimolecular and bimolecular reactions in the gas phase. The study develops the reaction potential energy surface for the unimolecular decomposition of 2-chlorophenol. The migration of the phenolic hydrogen to the ortho-C bearing the hydrogen atom produces 2-chlorocyclohexa-2,4-dienone through an activation barrier of 73.6 kcal/mol (0 K). This route holds more importance than the direct fission of Cl or the phenolic H. Reaction rate constants for the bimolecular reactions, 2-chlorophenol + X --> X-H + 2-chlorophenoxy (X = H, OH, Cl, O2) are calculated and compared with the available experimental kinetics for the analogous reactions of X with phenol. OH reaction with 2-chlorophenol produces 2-chlorophenoxy by direct abstraction rather than through addition and subsequent water elimination. The results of the present study will find applications in the construction of detailed kinetic models describing the formation of PCDD/F in the gas phase.     

Theoretical study on the methyl radical with chlorinated methyl radicals CH(3-n)Cln (n = 1, 2, 3) and CCl2

Radical-radical reactions involving chlorinated methyl radicals are particularly important in the mechanism of combustion of chlorinated hydrocarbons. Yet, they are usually difficult to study experimentally. In this paper, four chloride-related radical-radical reactions, i.e., CH3+CH(3-n)Cln (n = 1, 2, 3) and CH3+CCl2, are theoretically studied for the first time by means of the Gaussian-3//B3LYP potential energy surface survey combined with the master equation study over a wide range of temperatures and pressures. Our calculated results show that the three CH3+CH(3-n)Cln reactions can barrierlessly generate the former two kinetically allowed products P1 H(2)C=C(H)(3-n)Cl(n-1)+HCl and P2 CH3CH(3-n)Cl(n-1)+Cl with the very high predominance of P1 over P2. For the CH3 reaction with the biradical CCl2, which inevitably takes place during the CH3+CCl3 reaction and yet has never been studied experimentally or theoretically, H(2)C=CCl2+H and H(2)C=C(H)Cl+Cl are predicted to be the respective major and minor products. The results are compared with the recent laser photolysis/photoionization mass spectroscopy study on the CH3+CH(3-n)Cln (n = 1, 2, 3) reactions. The predicted rate constants and product branching ratios of the CH3+CCl2 reaction await future experimental verification.     

Mechanism of the

Stereochemical studies on [2 + 2] photoaddition of cis-/trans-4-propenylanisole (cis-1 and trans-1) and cis-1-(p-methoxyphenyl)ethylene-2-d(1) (cis-3-d(1)) to C(60) exhibit stereospecificity in favor of the trans-2 cycloadduct in the former case and nonstereoselectivity in the latter. The observed stereoselectivity in favor of the cis-6-d(3) [2 + 2] diastereomer by 12% in the case of the photochemical addition of (E)-1-(p-methoxyphenyl)-2-methyl-prop-1-ene-3,3,3-d(3) (trans-5-d(3)) to C(60) is attributed to a steric kinetic isotope effect (k(H)/k(D) = 0.78). The loss of stereochemistry in the cyclobutane ring excludes a concerted addition and is consistent with a stepwise mechanism. Intermolecular secondary kinetic isotope effects of the [2 + 2] photocycloaddition of 3-d(0) vs 3-d(1), and 3-d(6) as well as 5-d(0) vs 5-d(1), and 5-d(6) to C(60) were also measured. The intermolecular competition due to deuterium substitution of both vinylic hydrogens at the beta-carbon of 3 exhibits a substantial inverse alpha-secondary isotope effect k(H)/k(D) = 0.83 (per deuterium). Substitution with deuterium at both vinylic methyl groups of 5 yields a small inverse k(H)/k(D) = 0. 94. These results are consistent with the formation of an open intermediate in the rate-determining step.     

On the theory of the CO+OH reaction, including H and C kinetic isotope effects

The effect of pressure, temperature, HD isotopes, and C isotopes on the kinetics of the OH+CO reaction are investigated using Rice-Ramsperger-Kassel-Marcus theory. Pressure effects are treated with a step-ladder plus steady-state model and tunneling effects are included. New features include a treatment of the C isotope effect and a proposed nonstatistical effect in the reaction. The latter was prompted by existing kinetic results and molecular-beam data of Simons and co-workers on incomplete intramolecular energy transfer to the highest vibrational frequency mode in HOCO(*). In treating the many kinetic properties two small customary vertical adjustments of the barriers of the two transition states were made. The resulting calculations show reasonable agreement with the experimental data on (1) the pressure and temperature dependence of the HD effect, (2) the pressure-dependent (12)C(13)C isotope effect, (3) the strong non-Arrhenius behavior observed at low temperatures, (4) the high-temperature data, and (5) the pressure dependence of rate constants in various bath gases. The kinetic carbon isotopic effect is usually less than 10 per mil. A striking consequence of the nonstatistical assumption is the removal of a major discrepancy in a plot of the k(OH+CO)k(OD+CO) ratio versus pressure. A prediction is made for the temperature dependence of the OD+CO reaction in the low-pressure limit at low temperatures.