Effective rate constants for the surface reaction between solid methanol and deuterium atoms at 10 K

The surface reactions of CH3OH, CH2DOH, and CHD2OH with cold D atoms at 10 K were investigated using an atomic beam source and FTIR. Methyl-deuterated isotopologues CH2DOH, CHD2OH, and CD3OH were produced by exposure of amorphous solid CH3OH to D atoms at 10 K, and the pseudo-first-order rates for the reactions CH3OH + D --> CH2OH + HD, CH2DOH + D --> CHDOH + HD, and CHD2OH + D --> CD2OH + HD were estimated. The ratios of the reaction rates of the second and third reactions to the first reaction were 0.69 +/- 0.11 and 0.52 +/- 0.14, respectively. The difference in reaction rates is thought to be due to a secondary kinetic isotope effect on the H-abstraction reaction from the methyl side by D atoms.     

Imaging the reaction dynamics of OH + CD4. 3. Isotope effects

Isotope effects in the reaction of hydroxyl radical with methane are investigated in a crossed-beam experiment. By exploiting different combinations of OH/OD + CH4/CD4/CHD3, a total of eight isotopically variant reactions are examined to decipher the dynamical consequence of both primary and secondary isotope effects. The most prominent observation is from isotopic substitution of the transferred atom (primary isotope effect), which yields more stretch-excited water product for the D atom case than H atom. The opposite is found, however, for excitations in bending and combination modes. The secondary isotope effects, from either substituted hydroxyl radical or the three H atoms that are not being abstracted, are relatively minor and manifest themselves in spreading the vibrational distribution of the water coproducts.     

Quasiclassical trajectory study of the C(1D) + H2 reaction and isotopomeric variants: kinetic isotope effect and CD/CH branching ratio

The recently proposed ab initio single-sheeted double many-body expansion potential energy for the methylene molecule has been used to perform quasiclassical trajectory (QCT) calculations for the title reaction. Thermal and initial state-specific (v = 0, j = 0) rate constants for the C((1)D) + H(2)/HD/D(2) reactions have been obtained over a wide range of temperatures. Cross sections for the reaction C((1)D) + H(2) and its deuterated isotopes have also been calculated, as well as the CD/CH branching ratios for the C((1)D) + HD reaction. It is found that the CD + H product channel in the C((1)D) + HD reaction is preferred relative to the CH + D channel. The estimated rate constants are predicted to be in the order k(H2) > k(HD) > k(D2) and the calculated cross sections and rate constants compared with available theoretical and experimental data.     

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.     

Tunneling chemical reactions D + H2 --> DH + H and D + DH --> D2 + H in solid D2-H2 and HD-H2 mixtures: an electron-spin-resonance study

Tunneling chemical reactions D + H2 --> DH + H and D + DH --> D2 + H in solid HD-H2 and D2-H2 mixtures were studied in the temperature range between 4 and 8 K. These reactions were initiated by UV photolysis of DI molecules doped in these solids for 30 s and followed by measuring the time course of electron-spin-resonance (ESR) intensities of D and H atoms. ESR intensity of D atoms produced by the photolysis decreases but that of H atoms increases with time. Time course of the D and H intensities has the fast and slow processes. The fast process, which finishes within approximately 300 s after the photolysis, is assigned to the reaction of D atom with one of its nearest-neighboring H2 molecules, D(H2)n(HD)(12-n) --> H(H2)(n-1)(HD)(13-n) or D(H2)n(D2)(12-n) --> H(HD)(H2)(n-1)(D2)(12-n) for 12 > or = n > or = 1. Rate constant for the D + H2 reaction between neighboring D atom-H2 molecule pair is determined to be (7.5 +/- 0.7) x 10(-3) s(-1) in solid HD-H2 and (1.3+/-0.3) x 10(-2) s(-1) in D2-H2 at 4.1 K, which is very close to that calculated based on the theory of chemical reaction in gas phase by Hancock et al. [J. Chem. Phys. 91, 3492 (1989)] and Takayanagi and Sato [J. Chem. Phys. 92, 2862 (1990)]. This rate constant was found to be independent of temperature up to 7 K within experimental error of +/-30%. The slow process is assigned to the reaction of D atom produced in a cage fully surrounded by HD or D2 molecules, D(HD)12 or D(D2)12. This D atom undergoes the D + DH reaction with one of its nearest-neighboring HD molecules in solid HD-H2 or diffuses to the neighbor of H2 molecules to allow the D + H2 reaction in solid HD-H2 and D2-H2. The former is the main channel in solid HD-H2 below 6 K where D atoms diffuse very slowly, whereas the latter dominates over the former above 6 K. Rate for the reactions in the slow process is independent of temperature below 6 K but increases with the increase in temperature above 6 K. We found that the increase is due to the increase in hopping rate of D atoms to the neighbor of H2 molecules. Rate constant for the D + DH reaction was found to be independent of temperature up to 7 K as well.     

Isotope effects on the unimolecular dissociation of ionized 3-methyl-2-butanol: reactions via a long-lived C-H-C hydrogen-bridged ion-neutral complex

A pronounced isotope effect causes metastable CD3CHOHCH(CH3)2+* ions to expell C3H6D2 in preference to C3H7D in a ratio of approximately 33:1; a number of related compounds show similar effects. High-level ab initio calculations suggest that the reactant alcohol molecular ion possesses an extraordinarily long alpha-carbon-carbon bond and that the reaction proceeds via the formation of an intermediate hydrogen-bridged complex of propane and ionized vinyl alcohol, in which the bridging hydrogen atom is almost midway between the two carbon termini. The isotopic preference reflects the difference between the zero-point vibrational energies of the isotopically different product pairs rather than kinetic isotope effects on the hydrogen atom transfer reactions that precede dissociation.     

Reactions of C(1D) with H2 and its deuterated isotopomers, a wave packet study

Using a Chebyshev wave packet method, total and state-resolved reaction probabilities (J=0) were calculated for the reactions of C(1D) with various hydrogen isotopomers (H2, D2, and HD, nu i=0, j i=0) on a recent ab initio potential energy surface. For all the isotopic variants, it was found that the initial state specified reaction probabilities have no energy threshold and are strongly oscillatory, indicative of the involvement of long-lived resonances in this barrierless reaction. The J=0 product vibrational and rotational distributions for all three isotopic reactions, and the CH/CD branching ratio for the C+HD reaction, show strong dependence on the collision energy, further underscoring the important role played by the resonances. The generally decaying vibrational distributions and highly excited rotational distributions, which corroborate an insertion mechanism, and the dominance of the CD+H channel in the C+HD reaction are consistent with existing experimental observations. Initial state specified integral cross sections and rate constants were estimated using a capture model. The estimated rate constants were found to be close and in the order kHD>kH2>kD2. Finally, a method to calculate branching ratio in the C+HD reaction is proposed.     

钒促铑基催化剂上合成气反应中的H2／D2同位素效应

铑基催化剂可催化合成气生成甲烷、乙醇等,CO中C—O键断裂是该反应的关键步骤。Ⅴ助剂可显著提高Rh的催化活性,其作用本质被认为是促进了CO的直接解离。本文考察了Ⅴ助剂对Rh催化性能的影响及Rh/SiO_2、Rh-V(1:1)/SiO_2上合成气反应的H_2/D_2同位素效应,根据实验结果及键级守恒-Morse势法的计算结果初步探讨了Ⅴ助剂的作用本质。     

Deuterium enrichment of interstellar methanol explained by atom tunneling

We calculate, down to low temperature and for different isotopes, the reaction rate constants for the hydrogen abstraction reaction H + H(3)COH → H(2) + CH(2)OH/CH(3)O. These explain the known abundances of deuterated forms of methanol in interstellar clouds, where CH(2)DOH can be almost as abundant as CH(3)OH. For abstraction from both the C- and the O-end of methanol, the barrier-crossing motion involves the movement of light hydrogen atoms. Consequently, tunneling plays a dominant role already at relatively high temperature. Our implementation of harmonic quantum transition state theory with on the fly calculation of forces and energies accounts for these tunneling effects. The results are in good agreement with previous semiclassical and quantum dynamics calculations (down to 200 K) and experimental studies (down to 295 K). Here we extend the rate calculations down to lower temperature: 30 K for abstraction from the C-end of methanol and 80 K for abstraction from the OH-group. At all temperatures, abstraction from the C-end is preferred over abstraction from the O-end, more strongly so at lower temperature. Furthermore, the tunneling behavior strongly affects the kinetic isotope effects (KIEs). D + H(3)COH → HD + CH(2)OH has a lower vibrationally adiabatic barrier than H + H(3)COH → H(2) + CH(2)OH, giving rise to an inverse KIE (k(H)/k(D) < 1) at high temperature, in accordance with previous experiments and calculations. However, since tunneling is more facile for the light H atom, abstraction by H is favored over abstraction by D below ~135 K, with a KIE k(H)/k(D) of 11.2 at 30 K. The H + D(3)COD → HD + CD(2)OD reaction is calculated to be much slower than the D + H(3)COH → HD + CH(2)OH, in agreement with low-temperature solid-state experiments, which suggests the preference for H (as opposed to D) abstraction from the C-end of methanol to be the mechanism by which interstellar methanol is deuterium-enriched.     

Quantum reactive scattering of H + hydrocarbon reactions

A practical quantum-dynamical method is described for predicting accurate rate constants for general chemical reactions. The ab initio potential energy surfaces for these reactions can be built from a minimal number of grid points (average of 50 points) and expressed in terms of analytical functionals. All the degrees of freedom except the breaking and forming bonds are optimised using the MP2 method with a cc-pVTZ basis set. Single point energies are calculated on the optimised geometries at the CCSD(T) level of theory with the same basis set. The dynamics of these reactions occur on effective reduced dimensionality hyper-surfaces accounting for the zero-point energy of the optimised degrees of freedom. Bonds being broken and formed are treated with explicit hyperspherical time independent quantum dynamics. Application of the method to the H + CH(4)--> H(2)+ CH(3), H + C(2)H(6)--> H(2)+ C(2)H(5), H + C(3)H(8)--> H(2)+n-C(3)H(7)/H(2)+i-C(3)H(7) and H + CH(3)OH --> H(2)+ CH(3)O/H(2)+ CH(2)OH reactions illustrate the potential of the approach in predicting rate constants, kinetic isotope effects and branching ratios. All studied reactions exhibit large quantum tunneling in the rate constants at lower temperatures. These quantum calculations compare well with the experimental results.     

H(D)-atom yields in the quenching of Xe(6s[3/2]1) by methane, ethane, ethene, ethyne, and their deuterated isotopologues

The yields for the production of H(D) atoms in the reactions of Xe(6s[3/2]1) with simple hydrocarbons and their deuterated variants were determined. Xe(6s[32](1)) was produced by two-photon laser excitation of Xe(6p[1/2]0) followed by concomitant amplified spontaneous emission. H(D) atoms are detected using a vacuum-ultraviolet laser-induced fluorescence (LIF) technique. The H(D)-atom yields were evaluated from the LIF intensities and the overall rate constants for the quenching, which were determined from the temporal profile measurements of the resonance fluorescence from Xe(6s[3/2](1)). HD isotope effects were observed not only in the overall rate constants but also in the H(D)-atom yields. The yields for CH4, C2H4, and C2H2 were determined to be 0.89, 1.43, 1.03, respectively, while those for CD4, C2D4, and C2D2 were found to be smaller; 0.63, 0.86, and 0.79, respectively. The HD yield ratio for CH2D2 was 1.76. The presence of the isotope effects both in the rate constants and the yields suggests that electronic-to-electronic energy transfer processes and abstractive processes are competing.     

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.     

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.     

Mechanistic and kinetic study of the O + CH2OH reaction

The mechanism for the O + CH2OH reaction was investigated by various ab initio quantum chemistry methods. For the chemical activation mechanism, that is, the addition/elimination path, the couple-cluster methods including CCSD and CCSD(T) were employed with the cc-pVXZ (X = D, T, Q, 5) basis sets. For the abstraction channels, multireference methods including CASSCF, CASPT2, and MRCISD were used with the cc-pVDZ and cc-pVTZ basis sets. It has been shown that the production of H + HCOOH is the major channel in the chemical activation mechanism. The minor channels include HCO + H2O and OH + CH2O. The hydrogen abstraction by an O atom from the CH2OH radical produces either OH + CH2O or OH + HCOH. Moreover, the two abstraction reactions are essentially barrierless processes. The rate constants for the association of O with CH2OH have been calculated using the flexible transition state theory. A weak negative temperature dependence of the rate constants is found in the range 250-1000 K. Furthermore, it is estimated that the abstraction processes also play an important role in the O + CH2OH reaction. Additionally, the falloff behavior for the OCH2OH --> H + HCOOH reaction has been investigated. The present theoretical results are compared to the experimental measurements to understand the mechanism and kinetic behavior of the O + CH2OH reaction and the unimolecular reaction of the OCH2OH radical.     

Kinetics of the reaction of the heaviest hydrogen atom with H2, the 4Heμ + H2 → 4HeμH + H reaction: experiments, accurate quantal calculations, and variational transition state theory, including kinetic isotope effects for a factor of 36.1 in isotopic mass

The neutral muonic helium atom (4)Heμ, in which one of the electrons of He is replaced by a negative muon, may be effectively regarded as the heaviest isotope of the hydrogen atom, with a mass of 4.115 amu. We report details of the first muon spin rotation (μSR) measurements of the chemical reaction rate constant of (4)Heμ with molecular hydrogen, (4)Heμ + H(2) → (4)HeμH + H, at temperatures of 295.5, 405, and 500 K, as well as a μSR measurement of the hyperfine coupling constant of muonic He at high pressures. The experimental rate constants, k(Heμ), are compared with the predictions of accurate quantum mechanical (QM) dynamics calculations carried out on a well converged Born-Huang (BH) potential energy surface, based on complete configuration interaction calculations and including a Born-Oppenheimer diagonal correction. At the two highest measured temperatures the agreement between the quantum theory and experiment is good to excellent, well within experimental uncertainties that include an estimate of possible systematic error, but at 295.5 K the quantum calculations for k(Heμ) are below the experimental value by 2.1 times the experimental uncertainty estimates. Possible reasons for this discrepancy are discussed. Variational transition state theory calculations with multidimensional tunneling have also been carried out for k(Heμ) on the BH surface, and they agree with the accurate QM rate constants to within 30% over a wider temperature range of 200-1000 K. Comparisons between theory and experiment are also presented for the rate constants for both the D + H(2) and Mu + H(2) reactions in a novel study of kinetic isotope effects for the H + H(2) reactions over a factor of 36.1 in isotopic mass of the atomic reactant.     

Deuterium kinetic isotope effects on the gas-phase reactions of C2H with H2(D2) and CH4(CD4)

Kinetics of the ethynyl (C(2)H) radical reactions with H(2), D(2), CH(4) and CD(4) was studied over the temperature range of 295-396 K by a pulsed laser photolysis/chemiluminescence technique. The C(2)H radicals were generated by ArF excimer-laser photolysis of C(2)H(2) or CF(3)C(2)H and were monitored by the chemiluminescence of CH(A(2)Δ) produced by their reaction with O(2) or O((3)P). The measured absolute rate constants for H(2) and CH(4) agreed well with the available literature data. The primary kinetic isotope effects (KIEs) were determined to be k(H(2))/k(D(2)) = 2.48 ± 0.14 and k(CH(4))/k(CD(4)) = 2.45 ± 0.16 at room temperature. Both of the KIEs increased as the temperature was lowered. The KIEs were analyzed by using the variational transition state theory with semiclassical small-curvature tunneling corrections. With anharmonic corrections on the loose transitional vibrational modes of the transition states, the theoretical predictions satisfactorily reproduced the experimental KIEs for both C(2)H + H(2)(D(2)) and C(2)H + CH(4)(CD(4)) reactions.     

Direct ab initio dynamics study on the rate constants and kinetic isotope effect for the reactions of H atoms with GeDn(CH3)4-n (n = 1-4)

Direct ab initio dynamic calculations are performed on the reactions of atomic hydrogen with GeD(n)(CH(3))(4-n) (n = 1-4) over the temperature range 200-2000 K at the PMP4SDTQ/6-311 +G(3df,2p)//MP2/6-31 +G(d) (for n = 2-4) and G2//MP2/6-31 +G(d) (for n = 1) levels. The corresponding k(H)/k(D) ratios are then calculated in order to determine the kinetic isotope effect for the four reactions. For the simplest GeD(4) +H reaction, the only one that has available experimental data, the calculated canonical variational transition state theory incorporates small-curvature tunneling correction (CVT/SCT) thermal rate constants, and the k(H)/k(D) values are in good agreement with the experimental values within the experimental temperature range 293-550 K. For the four GeD(n)(CH(3))(4-4) (n = 1-4) reactions, the variational effect is small over the whole temperature range, whereas the small-curvature effect is important in the lower temperature range. Finally, the overall rate constants are fitted to the three-parameter expression over the whole temperature range 200-2000 K as 5.8 x 10(8)T(1.68)exp(-929/T), 1.7 x 10(8)T(1.80)exp(-691/T), 2.58 x 10(8)T(1.71)exp(-706/T), and 1.0 x 10(7)T(2.08)exp(-544/T) cm(3) mol(-1) s(-1) for the n = 4, 3, 2, and 1 reactions. Our work may represent the first theoretical study of the kinetic isotope effect for the H-attack on the G-H bonding.     

Deuterium kinetic isotope effects in gas-phase S(N)2 and E2 reactions: comparison of experiment and theory

The competition between bimolecular nucleophilc substitution and base-induced elimination is investigated through kinetic isotope effect measurements for gas-phase reactions of RCl + ClO- (R = methyl, ethyl, isopropyl, and tert-butyl) utilizing a FA-SIFT instrument. The overall reaction rate constants and the kinetic isotope effect for the reaction of C2H5Cl + ClO- are compared to computational results. [Hu, W. P.; Truhlar, D. G. J. Am. Chem. Soc. 1996, 118, 860.] Experimental results show that as the degree of substitution in the neutral reactant increases the E2 channel becomes dominant. The systematic change in the overall kinetic isotope effects indicates that, for the reaction of ClO- with C2H5Cl, both the SN2 and E2 pathways do occur, as predicted by computation; however the experimental reaction rate constants and KIE deviate strongly from the computational result.     

CO2 reforming of CH4 on Ni(111): a density functional theory calculation

CO(2) reforming of CH(4) on Ni(111) was investigated by using density functional theory. On the basis of thermodynamic analyses, the first step is CH(4) sequential dissociation into surface CH (CH(4) --> CH(3) --> CH(2) --> CH) and hydrogen, and CO(2) dissociation into surface CO and O (CO(2) --> CO + O). The second step is CH oxygenation into CHO (CH + O --> CHO), which is more favored than dissociation into C and hydrogen (CH --> C + H). The third step is the dissociation of CHO into surface CO and H (CHO --> CO + H). This can explain the enhanced selectivity toward the formation of CO and H(2) on Ni catalysts. It is found that surface carbon formation by the Bouduard back reaction (2CO = C((ads)) + CO(2)) is more favored than by CH(4) sequential dehydrogenation. The major problem of CO(2) reforming of CH(4) is the very strong CO adsorption on Ni(111), which results in the accumulation of CO on the surface and hinders the subsequent reactions and promotes carbon deposition. Therefore, promoting CO desorption should maintain the reactivity and stability of Ni catalysts. The computed energy barriers of the most favorable elementary reaction identify the CH(4) activation into CH(3) and H as the rate-determining step of CO(2) reforming of CH(4) on Ni(111), in agreement with the isotopic experimental results.