Direct Ab Initio Dynamics Study on the Reaction of CH3CHF2 (HFC‐152a) with the Cl Atom

A direct ab initio dynamics method is used to investigate the hydrogen‐abstraction reaction CH₃CHF₂+Cl. One transition state is located for α‐H abstraction, and two are identified for β‐H abstraction. The potential‐energy surface (PES) is obtained at the G3(MP2)//MP2/6‐311G(d, p) level. Furthermore, the rate constants of the three channels are evaluated by using canonical variational transition‐state theory (CVT) with small‐curvature tunneling (SCT) contributions over a wide temperature range of 200–2500 K. The dynamic calculations show that the reaction proceeds mainly by α‐H abstraction over the whole temperature range. The calculated rate constants and branching ratios are both in good agreement with the available experimental values.     

Kinetics for the hydrogen‐abstraction of CH4 with NO2

The detailed mechanism of the NO₂+CH₄ reaction has been computationally investigated at the M06‐2X/MG3S, B3LYP/6‐311G(2d,d,p), and MP2/6‐311+G(2df,p) levels. The direct dynamics calculations were preformed using canonical transition state theory with tunneling correction and scaled generalized normal‐mode frequencies including anharmonic torsion. The calculated results indicate that the NO₂+CH₄ reaction proceeds by three distinct channels simultaneously, leading to the formation of trans‐HONO (1a), cis‐HONO (1b), and HNO₂ (1c), and each channel involves the formation of intermediate having lower energy than the final product. The anti‐Hammond behavior observed in channel 1a is well analyzed. Proper treatment of anharmonic torsions about the C···H···O (or N) axis in the transition structures greatly improves the accuracy of kinetics predictions. The activation energy for each channel increases substantially with temperature, but is not strictly a linear function of temperature. Therefore, the thermal rate constants are fitted to the four‐parameter expression recommended for this case over the wide temperature range 400–4000 K. © 2012 Wiley Periodicals, Inc.     

On the kinetic mechanism of reactions of hydroxyl radical with CHF2CH3 − n F n (n = 1–3)

The potential energy surfaces of the reactions CHF₂CH_3 − n F _n (n = 1–3) + OH were investigated by MPWB1K and BMC-CCSD (single-point) methods. Furthermore, with the aid of canonical variational transition state theory including the small-curvature tunneling correction, the rate constants of the title reactions were calculated over a wide temperature range of 220–1,500 K. Agreement between the CVT/SCT rate constants and the experimental values is good. Our results show that the order of rate constants is CHF₂CH₂F + OH > CHF₂CHF₂ + OH > CHF₂CF₃ + OH. For reaction CHF₂CH₂F + OH, the 1-H-abstraction channel dominates the reaction at the whole temperature, while 2-H-abstraction channel appears to be competitive with the increase of temperature.     

Variational transition state theory rate constants and H/D kinetic isotope effects for CH3 + CH3OCOH reactions

The rate constants and H/D kinetic isotope effect for hydrogen abstraction reactions involving isotopomers of methyl formate by methyl radical are computed employing methods of the variational transition state theory (VTST) with multidimensional tunneling corrections. The energy paths were built with a dual-level method using the moller plesset second-order perturbation theory (MP2) method as the low-level and complete basis set (CBS) extrapolation as the high-level energy method. Benchmark calculations with the CBS_D-T approach give an enthalpy of reaction at 0 K for R1 (−4.5 kcal/mol) and R2 (−4.2 kcal/mol) which are in good agreement with the experiment, that is, −4.0 and − 4.8 kcal/mol. For the reactional paths involving the isotopomers CH₃ + CH₃OCOH → CH₄ + CH₃OCO and CH₃ + CH₃OCOD → CH₃D + CH₃OCO, the value of k^H/k^D (T = 455 K) using the canonical VTST/small-curvature tunneling approximation method is 6.7 in close agreement with experimental value (6.2). © 2019 Wiley Periodicals, Inc.     

Direct dynamics study on the mechanism and the kinetics of the reaction of CH3NH2 with OH

A theoretical study of the mechanism and the kinetics for the hydrogen abstraction reaction of methylamine by OH radical has been presented at the CCSD(T)/6‐311 ++G(2d,2p)//CCSD/6‐31G(d) level of theory. Our theoretical calculations suggest a stepwise mechanism involving the formation of a prereactant complex in the entrance channel and a preproduct complex in the exit channel, for the two hydrogen abstraction channels involving the methyl and amine groups. For clarity, the diagram of potential for the reaction is given. The calculated standard reaction enthalpies are ?98.48 and ?76.50 kJ mol^?1 and barrier heights are 0.36 and 25.25 kJ mol^?1, respectively. The rate constants are evaluated by means of the improved canonical variational transition state theory with small‐curvature tunneling correction (ICVT/SCT) in the temperature range of 299–3000 K. The calculated results show that the rate constants at experimentally measured temperatures are in good agreement with the experimental values. It is shown that the calculated rate constants exhibit a non‐Arrhenius behavior. Moreover, the variational effect is obvious in the calculated temperature range. The dominant product channel is to form CH₂NH₂ and H₂O via hydrogen abstraction from the CH₃ group of CH₃NH₂ by OH in the calculated temperature range. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

A dual‐level direct dynamics study on the hydrogen abstraction reaction of oxygen atom with methylhydrazine

The mechanism of the multichannel reaction CH₃NHNH₂ (SC₁ and SC₂) + O → products is investigated theoretically using ab initio and density functional theory, and dynamics properties are explored by a dual‐level direct dynamics method. The calculation of the potential energy surface is carried out at the BMC‐CCSD//MPW1K/6‐311G(d,p) level. Using canonical variational transition state theory with a small‐curvature tunneling correction, the rate constants of each channel are evaluated over a wide temperature range of 200–2000 K on the basis of obtained electronic structures and energy information. The total rate constants are calculated from the sum of the individual rate constants taking into account the Boltzmann distribution of two conformers. The reactivity of the H atom located in different groups is compared. © 2013 Wiley Periodicals, Inc.     

A Computational Study of the Kinetics and Mechanism for the C2H3 + CH3OH Reaction

The mechanism for the C₂H₃ + CH₃OH reaction has been investigated by the Gaussian‐4 (G4) method based on the geometric parameters of the stationary points optimized at the B3LYP/6–31G(2df, p) level of theory. Four transition states have been identified for the production of C₂H₄ + CH₃O (TSR/P1), C₂H₄ + CH₂OH (TSR/P2), C₂H₃OH + CH₃ (TSR/P3), and C₂H₃OCH₃ + H (TSR/P4) with the corresponding barriers 8.48, 9.25, 37.62, and 34.95 kcal/mol at the G4 level of theory, respectively. The rate constants and branching ratios for the two lower energy H‐abstraction reactions were calculated using canonical variational transition state theory with the Eckart tunneling correction at the temperature range 300–2500 K. The predicted rate constants have been compared with existing literature data, and the uncertainty has been discussed. The branching ratio calculation suggests that the channel producing CH₃O is dominant up to about 1070 K, above which the channel producing CH₂OH becomes very competitive.     

Theoretical studies on the reactions CH3SCH3 with OH,CF3, and CH3 radicals

The multiple‐channel reactions OH + CH₃SCH₃ → products, CF₃ + CH₃SCH₃ → products, and CH₃ + CH₃SCH₃ → products are investigated by direct dynamics method. The optimized geometries, frequencies, and minimum energy path are all obtained at the MP2/6‐31+G(d,p) level, and energetic information is further refined by the MC‐QCISD (single‐point) method. The rate constants for eight reaction channels are calculated by the improved canonical variational transition state theory with small‐curvature tunneling contribution over the temperature range 200–3000 K. The total rate constants are in good agreement with the available experimental data and the three‐parameter expressions k₁ = 4.73 × 10^?16T^1.89 exp(?662.45/T), k₂ = 1.02 × 10^?32T^6.04 exp(933.36/T), k₃ = 3.98 × 10^?35T^6.60 exp(660.58/T) (in unit of cm³ molecule^?1 s^?1) over the temperature range of 200–3000 K are given. Our calculations indicate that hydrogen abstraction channels are the major channels and the others are minor channels over the whole temperature range. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

Ligand Effects on the Mechanisms of Thermal Bond Activation in the Gas‐Phase Reactions NiX+/CH4→Ni(CH3)+/HX (X=H,CH3, OH,F). Short Communication

The thermal ion‐molecule reactions NiX⁺+CH₄→Ni(CH₃)⁺+HX (X=H, CH₃, OH, F) have been studied by mass spectrometric methods, and the experimental data are complemented by density functional theory (DFT)‐based computations. With regard to mechanistic aspects, a rather coherent picture emerges such that, for none of the systems studied, oxidative addition/reductive elimination pathways are involved. Rather, the energetically most favored variant corresponds to a σ‐complex‐assisted metathesis (σ‐CAM). For X=H and CH₃, the ligand exchange follows a ‘two‐state reactivity (TSR)’ scenario such that, in the course of the thermal reaction, a twofold spin inversion, i.e., triplet→singlet→triplet, is involved. This TSR feature bypasses the energetically high‐lying transition state of the adiabatic ground‐state triplet surface. In contrast, for X=F, the exothermic ligand exchange proceeds adiabatically on the triplet ground state, and some arguments are proposed to account for the different behavior of NiX⁺/Ni(CH₃)⁺ (X=H, CH₃) vs. NiF⁺. While the couple Ni(OH)⁺/CH₄ does not undergo a thermal ligand switch, the DFT computations suggest a potential‐energy surface that is mechanistically comparable to the NiF⁺/CH₄ system. Obviously, the ligands X act as a mechanistic distributor to switch between single vs. two‐state reactivity patterns.     

On the kinetic mechanism of the hydrogen abstraction reactions of the hydroxyl radical with CH3CF2Cl and CH3CFCl2: A dual level direct dynamics study

By means of the dual‐level direct dynamics method, the mechanisms of the reactions, CH₃CF₂Cl + OH → products (R1) and CH₃CFCl₂ + OH → products (R2), are studied over a wide temperature range 200–2000 K. The optimized geometries and frequencies of the stationary points are calculated at the MP2/6‐311G(d,p) level, and then the energy profiles of the reactions are refined with the interpolated single‐point energy method at the G3(MP2) level. The canonical variational transition‐state theory with the small‐curvature tunneling (SCT) correction method is used to calculate the rate constants. For the title reactions, three reaction channels are identified and the H‐abstraction channel is the major pathway. The results indicate that F substitution has a significant (reductive) effect on hydrochlorofluorocarbon reactivity. Also, for all H‐abstraction reaction channels the variational effect is small and the SCT effect is only important in the lower temperature range on the rate constants calculation. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Transition state theory thermal rate constants and RRKM‐based branching ratios for the N(2D) + CH4 reaction based on multi‐sState and multi‐reference Ab Initio calculations of interest for the titan's chemistry

Multireference single and double configuration interaction (MRCI) calculations including Davidson (+Q) or Pople (+P) corrections have been conducted in this work for the reactants, products, and extrema of the doublet ground state potential energy surface involved in the N( ² D) + CH₄ reaction. Such highly correlated ab initio calculations are then compared with previous PMP4, CCSD(T), W1, and DFT/B3LYP studies. Large relative differences are observed in particular for the transition state in the entrance channel resolving the disagreement between previous ab initio calculations. We confirm the existence of a small but positive potential barrier (3.86 ± 0.84 kJ mol^?1 (MR‐AQCC) and 3.89 kJ mol^?1 (MRCI+P)) in the entrance channel of the title reaction. The correlation is seen to change significantly the energetic position of the two minima and five saddle points of this system together with the dissociation channels but not their relative order. The influence of the electronic correlation into the energetic of the system is clearly demonstrated by the thermal rate constant evaluation and it temperature dependance by means of the transition state theory. Indeed, only MRCI values are able to reproduce the experimental rate constant of the title reaction and its behavior with temperature. Similarly, product branching ratios, evaluated by means of unimolecular RRKM theory, confirm the NH production of Umemoto et al., whereas previous works based on less accurate ab initio calculations failed. We confirm the previous findings that the N( ² D) + CH₄ reaction proceeds via an insertion–dissociation mechanism and that the dominant product channels are CH₂NH + H and CH₃ + NH. © 2012 Wiley Periodicals, Inc.     

Theoretical study and rate constants calculation for the reactions of SiH3 radical with SiH3CH3 and SiH2(CH3)2

The multiple‐channel reactions SiH₃ + SiH₃CH₃ → products and SiH₃ + SiH₂(CH₃)₂ → products are investigated by direct dynamics method. The minimum energy path (MEP) is calculated at the MP2/6‐31+G(d,p) level, and energetic information is further refined by the MC‐QCISD method. The rate constants for individual reaction channels are calculated by the improved canonical variational transition state theory (ICVT) with small‐curvature tunneling (SCT) correction over the temperature range of 200–2400 K. The theoretical three‐parameter expression k₁(T) = 2.39 × 10⁻²³T^4.01exp(−2768.72/T) and k₂(T) = 9.67 × 10⁻²⁷T^4.92exp(−2165.15/T) (in unit of cm³ molecule⁻¹ s⁻¹) are given. Our calculations indicate that hydrogen abstraction channel from SiH group is the major channel because of the smaller barrier height among eight channels considered. © 2009 Wiley Periodicals, Inc. J Comput Chem 2010     

Theoretical study of the mechanism of CH2CO + CN reaction

The potential energy surface information of the CH₂CO + CN reaction is obtained at the B3LYP/6‐311+G(d,p) level. To gain further mechanistic knowledge, higher‐level single‐point calculations for the stationary points are performed at the QCISD(T)/6‐311++G(d,p) level. The CH₂CO + CN reaction proceeds through four possible mechanisms: direct hydrogen abstraction, olefinic carbon addition–elimination, carbonyl carbon addition–elimination, and side oxygen addition–elimination. Our calculations demonstrate that R→IM1→TS3→P3: CH₂CN + CO is the energetically favorable channel; however, channel R→IM2→TS4→P4: CH₂NC + CO is considerably competitive, especially as the temperature increases (R, IM, TS, and P represent reactant, intermediate, transition state, and product, respectively). The present study may be helpful in probing the mechanism of the CH₂CO + CN reaction. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Mechanism and Rate Constants of the CH3+ CH2CO Reaction: A Theoretical Study

The mechanism of the reaction of ketene with methyl radical has been studied by ab initio CCSD(T)‐F12/cc‐pVQZ‐f12//B2PLYPD3/6‐311G** calculations of the potential energy surface. Temperature‐ and pressure‐dependent reaction rate constants have been computed using the Rice–Ramsperger–Kassel–Marcus (RRKM)–Master Equation and transition state theory methods. Three main channels have been shown to dominate the reaction; the formation of the collisionally stabilized CH₃COCH₂ radical and the production of the C₂H₅ + CO and HCCO + CH₄ bimolecular products. Relative contributions of the CH₃COCH₂, C₂H₅ + CO, and HCCO + CH₄ channels strongly depend on the reaction conditions; the formation of thermalized CH₃COCH₂ is favored at low temperatures and high pressures, HCCO + CH₄ is dominant at high temperatures, whereas the yield of C₂H₅ + CO peaks at intermediate temperatures around 1000 K. The C₂H₅ + CO channel is favored by a decrease in pressure but remains the second most important reaction pathway after HCCO + CH₄ under typical flame conditions. The calculated rate constants at different pressures are proposed for kinetic modeling of ketene reactions in combustion in the form of modified Arrhenius expressions. Only rate constant to form CH₃COCH₂ depends on pressure, whereas those to produce C₂H₅ + CO and HCCO + CH₄ appeared to be pressure independent.     

Sensitivity of the H + CH3 → CH4 recombination rate constant to the shape of the CH stretching potential

Using a potential-energy surface obtained in part from ab initio calculations, the H + CH₃ → CH₄ bimolecular rate constant at T = 300 K is determined from a Monte Carlo classical trajectory study. Representing the CH stretching potential with a standard Morse function instead ofthe ab initio curve increases the calculated rate constant by an order of magnitude. The experimental recombination rate constant is intermediate of the rate constants calculated with the Morse and ab initio stretching potentials.Two properties of the H + CH₃ α CH₄ potential-energy surface which significantly affect the recombination rate constant are the shape of the CH stretching potential and the attenuation of the H₃CH bending frequencies. Ab initio calculations with a hierarchy of basis sets and treatment of electron correlation indicate the latter is properly described [13]. The exact shape of the CH stretching potential is not delineated by the ab initio calculations, since the ab initio calculations are not converged for bond lengths of 2.0–3.0 Å [12]. However, the form of this stretching potential deduced from the highest-level ab initio calculations, and fit analytically by eq. (2), is significantly different from a Morse function. The experimental recombination rate constant is intermediate of the rate constants calculated with the Morse and ab initio CH stretching potentials. This indicates that the actual CH potential energy curve lies between the Morse and ab initio curves. This is consistent with the finding that potential energy curves for diatomics are not well described by a Morse function [12].     

Mechanism and rate constants of the CH2 + CH2CO reactions in triplet and singlet states: A theoretical study

Ab initio and density functional CCSD(T)-F12/cc-pVQZ-f12//B2PLYPD3/6-311G** calculations have been performed to unravel the reaction mechanism of triplet and singlet methylene CH₂ with ketene CH₂CO. The computed potential energy diagrams and molecular properties have been then utilized in Rice–Ramsperger–Kassel–Marcus-Master Equation (RRKM-ME) calculations of the reaction rate constants and product branching ratios combined with the use of nonadiabatic transition state theory for spin-forbidden triplet-singlet isomerization. The results indicate that the most important channels of the reaction of ketene with triplet methylene lead to the formation of the HCCO + CH₃ and C₂H₄ + CO products, where the former channel is preferable at higher temperatures from 1000 K and above. In the C₂H₄ + CO product pair, the ethylene molecule can be formed either adiabatically in the triplet electronic state or via triplet-singlet intersystem crossing in the singlet electronic state occurring in the vicinity of the CH₂COCH₂ intermediate or along the pathway of CO elimination from the initial CH₂CH₂CO complex. The predominant products of the reaction of ketene with singlet methylene have been shown to be C₂H₄ + CO. The formation of these products mostly proceeds via a well-skipping mechanism but at high pressures may to some extent involve collisional stabilization of the CH₃CHCO and cyclic CH₂COCH₂ intermediates followed by their thermal unimolecular decomposition. The calculated rate constants at different pressures from 0.01 to 100 atm have been fitted by the modified Arrhenius expressions in the temperature range of 300–3000 K, which are proposed for kinetic modeling of ketene reactions in combustion. © 2018 Wiley Periodicals, Inc.     

Dynamics of the CH3 + OH reaction

We study dynamics of the CH₃ + OH reaction over the temperature range of 300–2500 K using a quasiclassical method for the potential energy composed of explicit forms of short‐range and long‐range interactions. The explicit potential energy used in the study gives minimum energy paths on potential energy surfaces showing barrier heights, channel energies, and van der Waals well, which are consistent with ab initio calculations. Approximately, 20% of CH₃ + OH collisions undergo OH dissociation in a direct‐mode mechanism on a subpicosecond scale (<50 fs) with the rate coefficient as high as ～10^?10 cm³ molecule^?1 s^?1. Less than 10% leads to the formation of excited intermediates CH₃OH? with excess vibrational energies in CO and OH bonds. CH₃OH? stabilizes to CH₃OH, redissociates back to reactants, or forms one of various products after intramolecular energy redistribution via bond dissociation and formation on the time scale of 50–200 fs. The principal product is ¹CH₂ (k being ～10^?11), whereas ks for CH₂OH, CH₂O, and CH₃O are ～10^?12. The minor products are HCOH and CH₄ (k～10^?13). The total rate coefficient for CH₃ + OH → CH₃OH? → products is ～10^?11 and is weakly dependent on temperature. © 2011 Wiley Periodicals, Inc. Int J Chem Kinet 43: 455–466, 2011     

Theoretical Study of the Rate Coefficients for CH3NHNH2 + NO2 and Related Reactions

The kinetics and mechanisms of H atom abstraction reactions from CH₃NHNH₂ by NO₂ (R1) and related reactions have been investigated theoretically by using ωB97X‐D and CCSD(T)‐F12 quantum chemical calculations and the steady‐state unimolecular master equation analysis based on Rice–Ramsperger–Kassel–Marcus (RRKM) theory. For reaction (R1), both dissociation and isomerization steps between intermediate complexes were found to be important for the distribution of the dissociated bimolecular products. The dominant products of (R1) were found to be cis‐CH₃NHNH and HONO at lower temperature. The branching ratios for CH₃NNH₂ formation paths increased with increasing temperature. On the same reaction potential energy surface, six reactions including isomerization reactions between CH₃NNH₂ and cis‐/trans‐CH₃NHNH catalyzed by HONO were suggested to compete with the reverse reaction of (R1). The temperature‐ and pressure‐dependent rate expressions are proposed for kinetic modeling.     

Mechanism and kinetics properties for the reaction: Chloroethane with atomic O (P)

In this paper, we present direct dynamics calculations for the multiple-channel reaction of CH₃CH₂Cl with atomic O (³P) in a wide temperature range (200–3000 K), based on canonical variational transition state theory including small curvature corrections. Four distinct saddle points, one for α-abstraction and three for β-abstraction, have been located for this reaction. The potential energy surface information has been calculated at the MP2/6-311G(d,p) level. The energies along the minimum energy path have been further improved by single-point energy calculations at the G3MP2 level. In the β-abstraction channel, Jahn–Teller effect has been found. Changes of geometries, generalized normal-mode vibrational frequencies, and potential energies along the reaction paths for all channels have been discussed and compared. The calculated total rate constants match the available experimental values reasonable well over the measured temperature range. The results show the variational effect can be negligible and the small curvature tunneling contribution plays an important role for the calculation of the rate constant. At low temperature α-abstraction may be the major reaction channel, while β-abstraction will have more contribution to the whole reaction rate as the temperature increase.     

Ab initio direct dynamics studies on the reaction of H atom with CH3CH2Cl

The multiple channel reaction H + CH(3)CH(2)Cl --> products has been studied by the ab initio direct dynamics method. The potential energy surface information is calculated at the MP2/6-311G(d,p) level of theory. The energies along the minimum energy path are further improved by single-point energy calculations at the PMP4(SDTQ)/6-311+G(3df,2p) level of theory. For the reaction, four reaction channels (one chlorine abstraction, one alpha-hydrogen abstraction, and two beta-hydrogen abstractions) have been identified. The rate constants for each reaction channel are calculated by using canonical variational transition state theory incorporating the small-curvature tunneling correction in the temperature range 298-5000 K. The total rate constants, which are calculated from the sum of the individual rate constants, are in good agreement with the experimental data. The calculated temperature dependence of the branching fractions indicates that for the title reaction, H-abstraction reaction is the major reaction channel in the whole temperature range 298-5000 K.