Zum thermischen unimolekularen Zerfall von Methanol

Full-dimensional potential energy surface for the H + CH3OH reaction. Theoretical kinetics and dynamics study

The dynamics and kinetics of the abstraction reactions of hydrogen atoms with methanol have been studied using quasi-classical trajectory calculations and variational transition state theory with tunnelling corrections, based on a new analytical potential energy surface (PES). The new PES is a valence-bond/molecular mechanics (VB/MM) expression that provides us with the potential energy for any set of Cartesian coordinates. Two reaction channels are considered: hydrogen abstraction from the methyl group (R1) and hydrogen abstraction from the alcohol group (R2), R1 being much more likely to occur in the wide temperature range under study (250-1000 K), as expected from the lower barrier height. Our dynamic calculations at a collision energy of 20 kcal mol^-1 show that the H₂ co-product is produced mainly in its vibrational ground-state and little rotation excitation is found. As for our kinetic results, they agree with those from previous theoretical studies as well as with those from kinetic experimental results (rate constants and kinetic isotopic effects), lending confidence to the analytical PES presented here. Thus, we expect this PES to be a simple yet powerful tool to understand such an important reaction in combustion chemistry at very high temperatures and interstellar chemistry at very low temperatures.

Kinetics of the Reaction of Methyl Radical with Hydroxyl Radical and Methanol Decomposition

The CH3 + OH bimolecular reaction and the dissociation of methanol are studied theoretically at conditions relevant to combustion chemistry. Kinetics for the CH3 + OH barrierless association reaction and for the H + CH2OH and H + CH3O product channels are determined in the high-pressure limit using variable reaction coordinate transition state theory and multireference electronic structure calculations to evaluate the fragment interaction energies. The CH3 + OH --> 3CH2 + H2O abstraction reaction and the H2 + HCOH and H2 + H2CO product channels feature localized dynamical bottlenecks and are treated using variational transition state theory and QCISD(T) energies extrapolated to the complete basis set limit. The 1CH2 + H2O product channel has two dynamical regimes, featuring both an inner saddle point and an outer barrierless region, and it is shown that a microcanonical two-state model is necessary to properly describe the association rate for this reaction over a broad temperature range. Experimental channel energies for the methanol system are reevaluated using the Active Thermochemical Tables (ATcT) approach. Pressure dependent, phenomenological rate coefficients for the CH3 + OH bimolecular reaction and for methanol decomposition are determined via master equation simulations. The predicted results agree well with experimental results, including those from a companion high-temperature shock tube determination for the decomposition of methanol.

High-temperature rate constants for CH3OH + Kr --> products, OH + CH3OH --> products, OH + (CH3)(2)CO --> CH2COCH3 + H2O, and OH + CH3 --> CH) + H2O

The reflected shock tube technique with multipass absorption spectrometric detection of OH radicals at 308 nm (corresponding to a total path length of approximately 4.9 m) has been used to study the dissociation of methanol between 1591 and 2865 K. Rate constants for two product channels [CH3OH + Kr --> CH3 + OH + Kr (1) and CH3OH + Kr --> 1CH2 + H2O + Kr (2)] were determined. During the course of the study, it was necessary to determine several other rate constants that contributed to the profile fits. These include OH + CH3OH --> products, OH + (CH3)2CO --> CH2COCH3 + H2O, and OH + CH3 --> 1,3CH2 + H2O. The derived expressions, in units of cm(3) molecule(-1) s(-1), are k(1) = 9.33 x 10(-9) exp(-30857 K/T) for 1591-2287 K, k(2) = 3.27 x 10(-10) exp(-25946 K/T) for 1734-2287 K, kOH+CH3OH = 2.96 x 10-16T1.4434 exp(-57 K/T) for 210-1710 K, k(OH+(CH3)(2)CO) = (7.3 +/- 0.7) x 10(-12) for 1178-1299 K and k(OH+CH3) = (1.3 +/- 0.2) x 10(-11) for 1000-1200 K. With these values along with other well-established rate constants, a mechanism was used to obtain profile fits that agreed with experiment to within <+/-10%. The values obtained for reactions 1 and 2 are compared with earlier determinations and also with new theoretical calculations that are presented in the preceding article in this issue. These new calculations are in good agreement with the present data for both (1) and (2) and also for OH + CH3 --> products.

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.

Ab initiorate constants from hyperspherical quantum scattering: Application to H+C2H6 and H+CH3OH

The dynamics and kinetics of the abstraction reactions of H atoms with ethane and methanol have been studied using a quantum mechanical procedure. Bonds being broken and formed are treated with explicit hyperspherical quantum dynamics. The ab initio potential energy surfaces for these reactions have been developed from a minimal number of grid points (average of 48 points) and are given by analytical functionals. All the degrees of freedom except the breaking and forming bonds are optimized using the second order perturbation theory method with a correlation consistent polarized valence triple zeta basis set. Single point energies are calculated on the optimized geometries with the coupled cluster theory and the same basis set. The reaction of H with C2H6 is endothermic by 1.5 kcal/mol and has a vibrationally adiabatic barrier of 12 kcal/mol. The reaction of H with CH3OH presents two reactive channels: the methoxy and the hydroxymethyl channels. The former is endothermic by 0.24 kcal/mol and has a vibrationally adiabatic barrier of 13.29 kcal/mol, the latter reaction is exothermic by 7.87 kcal/mol and has a vibrationally adiabatic barrier of 8.56 kcal/mol. We report state-to-state and state-selected cross sections together with state-to-state rate constants for the title reactions. Thermal rate constants for these reactions exhibit large quantum tunneling effects when compared to conventional transition state theory results. For H+CH3OH, it is found that the CH2OH product is the dominant channel, and that the CH3O channel contributes just 2% at 500 K. For both reactions, rate constants are in good agreement with some measurements.