Theoretical calculations of energetics, structures, and rate constants for the H+CH3OH hydrogen abstraction reactions

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.

Theoretical studies of hydrogen abstraction from H2X and CH3XH (X = O, S) by trichloromethyl radicals

We have theoretically investigated the hydrogen abstraction reactions of H₂O, H₂S, CH₃OH, and CH₃SH by the CCl₃ radical, which is of interest in atmospheric chemistry research. In this study, mechanistic and kinetic analyses for the title reactions have been performed at the W1 and the CCSD(T)/cc-pVTZ//M06-2X/cc-pVTZ level. Standard Gibbs free energies for all reaction channels at the W1 level with ZPE corrections were also calculated. Intrinsic Reaction Coordinate (IRC) calculations were performed to verify the connectivity of all the transition states with the reactants and products. All rate coefficients are computed by conventional transition state theory (CTST) with the zero-curvature tunneling (ZCT) and also Wigner's tunneling (W) correction in the temperature range from 200 K to 2000 K. The rate coefficients for each reaction channel are also evaluated by canonical variational transition state theory (CVT) with the small-curvature tunneling correction method (SCT) in the same temperature range. Three-parameter Arrhenius expressions have been obtained by fitting to the computed rate coefficients of all abstraction channels between 200 and 2000 K. The branching ratios for these reactions have been determined. This study provides the first theoretical and kinetic determination of the CCl₃ rate coefficient for reactions with H₂O, H₂S, CH₃OH, and CH₃SH over a large temperature range.

Application of one-dimensional semiclassical transition state theory to the CH3OH + H ⇌ CH2OH/CH3O + H2 reactions

The rate constants of the two branches of H-abstractions from CH₃OH by the H-atom and the corresponding reactions in the reverse direction are calculated using the one-dimensional semiclassical transition state theory (1D SCTST). In this method, only the reaction mode vibration of the transition state (TS) is treated anharmonically, while the remaining internal degrees of freedom are treated as they would have been in a standard TS theory calculation. A total of eight ab initio single-point energy calculations are performed in addition to the computational cost of a standard TS theory calculation. This allows a second-order Richardson extrapolation method to be employed to improve the numerical estimation of the third- and fourth-order derivatives, which in turn are used in the calculation of the anharmonic constant. Hindered-rotor (HR) vibrations are identified in the equilibrium states of CH₃OH and CH₂OH, and the TSs of the reactions. The partition function of the HRs are calculated using both a simple harmonic oscillator model and a more sophisticated one-dimensional torsional eigenvalue summation (1D TES) method. The 1D TES method can be easily adapted in 1D SCTST computation. The resulting 1D SCTST with 1D TES rate constants show good agreement to previous theoretical and experimental works. The effects of the HR on rate constants for different reactions are also investigated.This article is part of the theme issue 'Modern theoretical chemistry'.

A novel assessment of the role of the methyl radical and water formation channel in the CH3OH + H reaction

A number of experimental and theoretical papers accounted almost exclusively for two channels in the reaction of atomic hydrogen with methanol: H-abstraction from the methyl (R1) and hydroxyl (R2) functional groups. Recently, several astrochemical studies claimed the importance of another channel for this reaction, which is crucial for kinetic simulations related to the abundance of molecular constituents in planetary atmospheres: methyl radical and water formation (R3 channel). Here, motivated by the lack of and uncertainties about the experimental and theoretical kinetic rate constants for the third channel, we developed first-principles Car-Parrinello molecular dynamics thermalized at two significant temperatures - 300 and 2500 K. Furthermore, the kinetic rate constant of all three channels was calculated using a high-level deformed-transition state theory (d-TST) at a benchmark electronic structure level. d-TST is shown to be suitable for describing the overall rate constant for the CH₃OH + H reaction (an archetype of the moderate tunnelling regime) with the precision required for practical applications. Considering the experimental ratios at 1000 K, k_R1/k_R2 ≈ 0.84 and k_R1/k_R3 ≈ 15-40, we provided a better estimate when compared with previous theoretical work: 7.47 and 637, respectively. The combination of these procedures explicitly demonstrates the role of the third channel in a significant range of temperatures and indicates its importance considering the thermodynamic control to estimate methyl radical and water formation. We expect that these results can help to shed new light on the fundamental kinetic rate equations for the CH₃OH + H reaction.

Signatures of a quantum diffusion limited hydrogen atom tunneling reaction

We are studying the details of hydrogen atom (H atom) quantum diffusion in parahydrogen quantum solids in an effort to better understand H atom transport and reactivity under these conditions.

Thermochemical and Kinetics of CH3SH + H Reactions: The Sensitivity of Coupling the Low and High-Level Methodologies

The reaction system formed by the methanethiol molecule (CH₃SH) and a hydrogen atom was studied via three elementary reactions, two hydrogen abstractions and the C-S bond cleavage (CH₃SH + H → CH₃S + H₂ (R1); → CH₂SH + H₂ (R2); → CH₃ + H₂S (R3)). The stable structures were optimized with various methodologies of the density functional theory and the MP2 method. Two minimum energy paths for each elementary reaction were built using the BB1K and MP2 methodologies, and the electronic properties on the reactants, products, and saddle points were improved with coupled cluster theory with single, double, and connected triple excitations (CCSD(T)) calculations. The sensitivity of coupling the low and high-level methods to calculate the thermochemical and rate constants were analyzed. The thermal rate constants were obtained by means of the improved canonical variational theory (ICVT) and the tunneling corrections were included with the small curvature tunneling (SCT) approach. Our results are in agreement with the previous experimental measurements and the calculated branching ratio for R1:R2:R3 is equal to 0.96:0:0.04, with k_R1 = 9.64 × 10^-13 cm³ molecule^-1 s^-1 at 298 K.

Thermochemical and kinetics studies of the CH3SH+S (3P) hydrogen abstraction and insertion reactions

Sulfur-containing molecules have a significant impact on atmosphere and biosphere. In this work we studied, from the point of view of electronic structure and chemical kinetics methods, the elementary reactions between a methanethiol molecule and a sulfur atom leading to hydrogen abstraction C-S bond cleavage (CH(3)SH+S; R1:→ CH(3)S+SH; R2: → CH(2)SH+SH; R3:→ CH(3)+HS(2)). The geometrical structures of the reactants, products, and saddle points for the three reaction paths were optimized using the BB1K method with the aug-cc-pV(T+d)Z basis set. The thermochemical properties were improved using single point coupled-cluster (CCSD(T)) calculations on the BB1K geometries followed by extrapolation to the complete basis set (CBS) limit. This methodology was previously applied and has given accurate values of thermochemical and kinetics properties when compared to benchmark calculations and experimental data. For each reaction, the thermal rate constants were calculated using the improved canonical variational theory (ICVT) including the zero-curvature (ICVT/ZCT) and small-curvature (ICVT/SCT) tunneling corrections. For comparison, the overall ICVT/SCT reaction rate constant at 300 K obtained with single-point CCSD(T)/CBS calculations for the CH(3)SH+S reaction is approximately 1400 times lower than the isovalent CH(3)SH+O reaction, obtained with CVT/SCT. The reaction path involving the hydrogen abstraction from the thiol group is the most important reactive path in all temperatures.