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.     

Canonical statistical adiabatic channel model calculations of the H + CH3 → CH4 recombination reaction on an ab initio potential energy surface. The role of the transitional modes

The canonical formalism of the statistical adiabatic channel model is used to calculate limiting high pressure rate constants for the H + CH₃ → CH₄ recombination reaction on a recently reported analytic potential energy surface based on ab initio calculations. An effective adiabatic channel potential which incorporates the G_?? matrix element of the twofold degenerate H₃C? H transitional bending mode, quartic anharmonicity, and state selected mode coupling effects is implemented. The rate constants calculated over the temperature range 200–1000 K are in very good agreement with recent canonical variational transition state theory calculations performed on the same surface. The comparison with experimental results is also discussed.     

Theoretical study of the reaction H + ClCH3 → HCl + CH3

The reaction H + ClCH₃ has theoretically studied in a LEPS potential energy surface with a single-particle approximation for the methyl group. The LEPS adjustable parameters were selected to reach a good agreement with experimental values of activation energy and exothermicity. A wide set of quasi-classical trajectories for that system has been calculated within a energy range covering the significative values of relative velocities at temperatures between 300 and 1000 K. Calculated reactive cross sections increase with translational energy and with the initial vibrational level, but they are not influenced by rotational excitation of the reactants. Microscopic and total reaction rate constants have been obtained within the temperature range and agree quite well with available experimental results. Final energy distribution shows that most of the exoergicity is consumed in increasing the relative velocity of the products, while HCl molecules remain in their vibrational ground state.     

Direct dynamics study of the hydrogen abstraction reaction of CF3CH2Cl + Cl → CF3CHCl + HCl

The hydrogen abstraction reaction of Cl atoms with CF₃CH₂Cl (HCFC‐133a) is investigated by using density function theory and ab initio approach, and the rate constants are calculated by using the dual‐level direct dynamics method. Optimized geometries and frequencies of reactants, transition state, and products are computed at the B3LYP/6‐311+G(2d,2p) level. To refine the energetic information along the minimum energy path, single‐point energy calculations are carried out at the G3(MP2) level of theory. The interpolated single‐point energy method is employed to correct the energy profiles for the title reaction. The rate constants are evaluated by using the canonical variational transition state theory with a small‐curvature tunneling correction over a wide range of temperature, 200–2000 K. The variational effect for the reaction is moderate at low temperatures and very small at high temperatures. However, the tunneling correction has an important contribution in the lower temperature range. The agreement between calculated rate constants and available experimental values is good at lower temperatures but diverges significantly at higher temperatures. © 2012 Wiley Periodicals, Inc. Int J Chem Kinet 44: 661–667, 2012     

Variational transition‐state theory study of the atmospheric reaction OH + O3 → HO2 + O2

We report variational transition‐state theory calculations for the OH + O₃→ HO₂ + O₂ reaction based on the recently reported double many‐body expansion potential energy surface for ground‐state HO₄ [Chem Phys Lett 2000, 331, 474]. The barrier height of 1.884 kcal mol^?1 is comparable to the value of 1.77–2.0 kcal mol^?1 suggested by experimental measurements, both much smaller than the value of 2.16–5.11 kcal mol^?1 predicted by previous ab initio calculations. The calculated rate constant shows good agreement with available experimental results and a previous theoretical dynamics prediction, thus implying that the previous ab initio calculations will significantly underestimate the rate constant. Variational and tunneling effects are found to be negligible over the temperature range 100–2000 K. The O₁? O₂ bond is shown to be spectator like during the reactive process, which confirms a previous theoretical dynamics prediction. © 2007 Wiley Periodicals, Inc. 39: 148–153, 2007     

Direct dynamic study on the hydrogen abstraction reaction of H2CO with NCO

A direct ab initio dynamics method is used to investigate the hydrogen‐abstraction reaction of H₂CO with NCO. The potential energy surface information is obtained at the MP2/6‐311G(d,p) level. More accurate single‐point energy is refined at the G3(MP2)//MP2/6‐311G(d,p) level. Furthermore, the rate constants of reaction H₂CO + NCO are evaluated by using the canonical variational transition state theory with small‐curvature tunneling contributions over a wide temperature range of 200–2000 K. The calculated reaction enthalpy and rate constants are in good agreement with the available experimental values. © 2009 Wiley Periodicals, Inc. Int J Chem Kinet 41: 394–400, 2009     

Reaction dynamics of Ca(4s3d 1D2)+CH4-->CaH(X 2Sigma+)+CH3: reaction pathway and energy disposal for the CaH product

The reaction pathway for Ca(4s3d 1D2)+CH4-->CaH(X 2Sigma+)+CH3 has been investigated by using a pump-probe technique in combination with potential-energy surface (PES) calculations. The nascent product distributions of CaH have been characterized with Boltzmann rotational temperatures of 1013+/-102 and 834+/-70 K for the v=0 and 1 levels, respectively, and a Boltzmann vibrational temperature of 1313+/-173 K. The rotational and vibrational energy partitions in CaH have been estimated to be 461+/-45 and 252+/-15 cm(-1), respectively. According to the PES calculations, the pathway favors an insertion mechanism. Ca(3 1D2) approaches CH4 along an attractive potential surface in a C2v (or Cs) symmetry and then the collision complex undergoes nonadiabatic transition to the reactive ground-state surface. An Arrhenius plot shows a potential-energy requirement of 2695+/-149 cm(-1), which accounts for the endothermicity of 2930 cm(-1) for the reaction scheme. The Ca-C bond distance in the transition state structure is short enough to allow for tight orbital overlap between CaH and CH3. The strong coupling between the moieties renders the energy transfer sufficient from CaH into the CH3 radical. As compared to the Ca(4 1P1) reaction, the dissociation lifetime of the intermediate complex with less excess energy is prolonged so as to cause much less vibrational energy disposal into CaH.     

Quasiclassical trajectory study of the reaction NH(X~3∑~-)+H→N(~4S)+H_2

The dynamics of the NH + H→N+H2 reaction has been investigated by means of the 3D quasiclassical trajectory approach by using the LEPS potential energy surface.The calculated rate coefficient is in good agreement with the experimental value.The reaction was found to occur via a direct channel.The product H2 has a cold excitation of rotational state,but has a reverse distribution of the vibrational state with a peak at v=1.Based on the potential energy surface and the trajectory analysis,the reaction mechanism has been explained successfully.     

An analytic hindered rotor model for calculating microcanonical variational unimolecular rate constants from reaction path properties

A model is proposed for performing microcanonical variational transition state theory calculations which incorporates ideas from vibrator and flexible variational transition state models. Vibrational frequencies, moments of inertia, and potential energy for the variational transition state are found by reaction path following as for the vibrator model. However, the transitional modes are treated as hindered rotors using an analytic potential and an analytic density of states, which are fit to barriers for hindered rotation determined from reaction path following. The model proposed here differs from the flexible transition state model in that the density of states for the transitional modes is analytic and transitional modes and external rotational angular momenta are uncoupled. For the H + CH₃ ? CH₄ system, rate constants calculated with this new model are only 6–23% smaller on average from those of the flexible transition state model for values of total angular momentum which correspond to average rotational temperatures of 0–2000 K. Harmonic frequencies calculated for the transitional modes from the hindered rotor Hamiltonian are in good agreement with the exact values found by a reaction path analysis. © 1994 John Wiley & Sons, Inc.     

Ab initio quantum chemical studies of the reactions of CF3CFHO2 with HO2

The potential energy surface (PES) for the CF₃CFHO₂+HO₂ reaction has been theoretically investigated using the DFT [B3LYP/6‐311G(d,p)] and B3LYP/6‐311++G(3df,3pd)//B3LYP/6‐311G(d,p) levels of theory. Both singlet and triplet PESs are investigated. The reaction mechanism on the triplet surface is simple. It is revealed that the formation of CF₃CFHOOH+³O₂ is the dominant channel on the triplet surface. On the basis of the ab initio data, the total rate constants for the reaction CF₃CFHO₂+HO₂ in the T = 210–500 K range have been computed using conventional transition state theory with Wigner's tunneling correction and have been fitted by a rate constant expression as k = 1.04 ×10^?12(cm³ molecule^?1 s^?1) exp (700.33/T). Calculated transition state rate constants with Wigner's tunneling correction for the reaction CF₃CFHO₂+HO₂ are in good agreement with the available experimental values. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

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     

CH3NHNH2 + OH reaction: Mechanism and dynamics studies

A direct dynamics study was carried out for the multichannel reaction of CH₃NHNH₂ with OH radical. Two stable Conformers (I, II) of CH₃NHNH₂ are identified by the rotation of the ? CH₃ group. For each conformer, five hydrogen‐abstraction channels are found. The reaction mechanisms of product radicals (CH₃NNH₂ and CH₃NHNH) with OH radical are also investigated theoretically. The electronic structure information on the potential energy surface is obtained at the B3LYP/6‐311G(d,p) level and the energetics along the reaction path is refined by the BMC‐CCSD method. Hydrogen‐bonded complexes are presented at both the reactant and product sides of the five channels, indicating that the reaction may proceed via an indirect mechanism. The influence of the basis set superposition error (BSSE) on the energies of all the complexes is discussed by means of the CBS‐QB3 method. The rate constants of CH₃NHNH₂ + OH are calculated using canonical variational transition‐state theory with the small‐curvature tunneling correction (CVT/SCT) in the temperature range of 200–1000 K. Slightly negative temperature dependence of rate constant is found in the temperature range from 200 to 345 K. The agreement between the theoretical and experimental results is good. It is shown that for Conformer I, hydrogen‐abstraction from ? NH? position is the primary pathway at low temperature; the hydrogen‐abstraction from ? NH₂ is a competitive pathway as the temperature increases. A similar case can be concluded for Conformer II. The overall rate constant is evaluated by considering the weight factors of each conformer from the Boltzmann distribution function, and the three‐term Arrhenius expressions are fitted to be k_T = 1.6 × 10^?24T^4.03exp (1411.5/T) cm³ molecule^?1 s^?1 between 200–1000 K. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2009     

反式2-甲基-2-丁烯酸甲酯与臭氧反应机理的计算研究

采用G3B3方法构建反式2-甲基-2-丁烯酸甲酯与O3反应体系以及后续Criegee自由基有、无水分子参与下异构化反应的势能面剖面.结果表明,反式2-甲基-2-丁烯酸甲酯与O3反应首先生成一个稳定的五元环中间体,此中间体按断键位置不同后续裂解反应存在两条路径,分别生成产物P1(CH3CHOO+CH3OC(O)C(CH3)O)和P2(CH3CHO+CH3OC(O)C(CH3)OO).利用经典过渡态理论(TST)并结合Wigner矫正模型计算了200-1200 K温度区间内标题反应的速率常数kTST/W.计算结果显示,294 K时,该反应速率常数为7.55×10-18cm3molecule-1s-1,与Bernard等对类似反应所测实验值非常接近.生成的Criegee自由基(CH3CHOO和CH3OC(O)C(CH3)OO)可分别与水分子发生α-加成及β-氢迁移反应,其中Criegee自由基与水的α-加成反应较其与水的β-氢迁移反应具有优势.另外与无水分子参与CH3CHOO和CH3OC(O)C(CH3)OO异构化反应相比,水分子的参与使得异构化反应较为容易进行.     

An experimental and theoretical study on the kinetics of the reaction between 4‐hydroxy‐3‐hexanone CH3CH2C(O)CH(OH)CH2CH3 and OH radicals

Experimental and theoretical rate coefficients are determined for the first time for the reaction of 4‐hydroxy‐3‐hexanone (CH₃CH₂C(O)CH(OH)CH₂CH₃) with OH radicals as a function of temperature. Experimental studies were carried out using two techniques. Absolute rate coefficients were measured using a cryogenically cooled cell coupled to the pulsed laser photolysis‐laser‐induced fluorescence technique with temperature and pressure ranges of 280‐365 K and 5‐80 Torr, respectively. Relative values of the studied reaction were measured under atmospheric pressure in the range of 298‐354 K by using a simulation chamber coupled to a FT‐IR spectrometer. In addition, the reaction of 4H3H with OH radicals was studied theoretically by using the density functional theory method over the range of 278‐350 K. Results show that H‐atom abstraction occurs more favorably from the C–H bound adjacent to the hydroxyl group with small barrier height. Theoretical rate coefficients are in good agreement with the experimental data. A slight negative temperature dependence was observed in both theoretical and experimental works. Overall, the results are deliberated in terms of structure–reactivity relationship and atmospheric implications.     

Measurements and modeling of DO2 formation in the reactions of C2D5 and C3D7 radicals with O2

Time-resolved production of HO2 and DO2 from the reactions of nondeuterated and deuterated ethyl and propyl radicals with O2 are measured as a function of temperature and pressure in the "transition region" between 623 and 748 K using the technique of laser photolysis/long path frequency modulation spectroscopy. Experimental measurements, using both pulsed-photolytic Cl-atom-initiated oxidation of ethane and propane and direct photolysis of ethyl, n-propyl, and isopropyl iodides, are compared to kinetic models based on the results of time-dependent master equation calculations with ab initio characterization of stationary points. The formation of DO2 and HO2 from the subsequent reaction of the alkyl radicals with O2 is followed by infrared frequency modulation spectroscopy. The concentration of I atoms is simultaneously monitored by direct absorption of a second laser probe on the spin-orbit transition. The kinetic models accurately describe the time scale and amplitude of the DO2 and HO2 formation resulting from C2D5 + O2, n-C3D7 + O2, i-C3D7 + O2, and i-C3H7 + O2. Overall, a very good level of agreement is found between theory and experiments over a wide range of temperatures, pressures, and O2 concentrations. Good agreement is also found between previous literature studies and the theory presented in this work except in the case of the high-temperature rate coefficients for the reaction of i-C3H7 + O2 to form propene. A reinvestigation of the high-temperature kinetics of the i-C3H7 + O2 reaction appears warranted. The results from the present work suggest that the theory for formation of HO2 from the reactions of ethyl and both isomeric forms of propyl radicals with O2 are very well established at this time. It is hoped that these reactions can now form the groundwork for the study and interpretation of larger and more complex R + O2 systems.     

Influence of reagent rotation on (H-, D2) and (D-, H2) collisions: a quantum mechanical study

Time-independent quantum mechanical (TIQM) approach (helicity basis truncated at k = 2) has been used for computing differential and integral cross sections for the exchange reaction H- + D2 (v = 0, j = 0-4) --> HD + D- and D- + H2 (v = 0, j = 0-3) --> HD + H- in three dimensions on an accurate ab initio potential energy surface. It is shown that the j-weighted differential reaction cross section values are in good agreement with the experimental results reported by Zimmer and Linder at four different relative translational energies (Etrans = 0.55, 0.93, 1.16 and 1.48 eV) for (H-, D2) and at one relative translational energy (Etrans = 0.6 eV) by Haufler et al. for both (H-, D2) and (D-, H2) collisions. The j-weighted integral reaction cross section values are in good agreement with the crossed beam measurements by Zimmer and Linder in the Etrans range 0.5-1.5 eV and close to the guided ion beam results by Haufler et al. for (H-, D2) in the range 0.8-1.2 eV. Time-dependent quantum mechanical (TDQM) results obtained using centrifugal sudden approximation are reported in the form of integral reaction cross section values as a function of Etrans in the range 0.3-3.0 eV for both reactions in three dimensions on the same potential energy surface. The TDQM reaction cross section values decline more sharply than the TIQM results with increase in the initial rotational quantum number (j) for the D2 molecules in their ground vibrational state (v = 0) for (H-, D2) collisions. The computed j-weighted reaction cross section values are in good agreement with the experimental results reported by Zimmer and Linder for (H-, D2) collisions and guided ion beam results by Haufler et al. for both (H-, D2) and (D-, H2) collisions for energies below the threshold for electron detachment channel.     

The Reaction X + Cl2→XCl + Cl (X = Mu,H, D). II. Comparison of experimental data with theoretical results derived from a new potential energy surface

We consider experimental implications for the Mu + Cl₂, H + Cl₂, and D + Cl₂ reactions of the extended London—Eyring—Polanyi—Sato (LEPS) potential energy surface derived from experimental data in paper I. In the present calculations, it is necessary to make additional implicit and explicit assumptions concerning the three-dimensional (3D) nature of the potential surface, since the inversion procedure of paper I yields information only on the collinear (1D) part of the surface. We have performed accurate 1D quantum calculations of reaction probabilities, which are then transformed into 3D by an information theoretic 1D → 3D transformation incorporating a constraint to allow for angular momentum transfer effects in light+heavy—heavy atom reactions. This procedure implicitly accounts for the 3D nature of the potential surface. The calculated vibrational and vibrotational product distributions are in good agreement with those determined in thermal chemiluminescence experiments. The Sato parameters for the 1D surface also define a full 3D surface. This is used as an approximation to the true surface, and its properties are explored in 3D quasiclassical trajectory calculations. Comparison is made for the H and D reactions with available chemiluminescence, molecular beam and kinetic experimental data for differential and total reaction cross sections, energy disposal, rate coefficients and Arrhenius parameters. Some kinetic isotope effects in the Mu, H, and D reactions are discussed using vibrationally adiabatic theory. Comparison is also made with results from other calculations in the literature for the H + Cl₂ and D + Cl₂ reactions.     

Theoretical study of the dynamics of Cl + O3 reaction I. Ab initio potential energy surface and quasiclassical trajectory results

We present a global full dimensional potential energy surface (PES) for the Cl + O(3)→ ClO + O(2) reaction, which is an elementary step in a catalytic cycle that leads to the destruction of ozone in the stratosphere. The PES is constructed by interpolation of quantum chemistry data using the method developed by Collins and co-workers. Ab initio data points (energy, gradients and Hessian matrix elements) have been calculated at the UQCISD/aug-cc-pVDZ (unrestricted quadratic configuration interaction with single and double excitations) level of theory. The ab initio calculations predict a markedly non-coplanar (dihedral angle of 80°) transition state for the reaction, located very early in the reactant valley and slightly below the energy of the reactants as long as the spin-orbit splitting is neglected. Quasiclassical trajectory (QCT) calculations have been carried out at several collision energies to investigate the reaction dynamics. The QCT excitation function shows no threshold, displays a minimum at a collision energy of 2.5 kcal mol(-1), and then increases monotonically at larger collision energies. This behaviour is consistent with a barrierless reaction dominated by an oxygen-abstraction mechanism. The calculated product vibrational distributions (strongly inverted for ClO) and rate constants are compared with experimental determinations. Differential cross sections (DCS) summed over all final states are found to be in fairly good agreement with those derived from crossed molecular beam experiments.     

A density functional theory study of the oxidative addition of methyl iodide to square planar [Rh(acac)(P(OPh)3)2] complex and simplified model systems

The transition state for the oxidative addition reaction [Rh(acac)(P(OPh)₃)₂] + CH₃I, as well as two simplified models viz. [Rh(acac)(P(OCH₃)₃)₂] and [Rh(acac)(P(OH)₃)₂], are calculated with the density functional theory (DFT) at the PW91/TZP level of theory. The full experimental model, as well as the simplified model systems, gives a good account of the experimental Rh-ligand bond lengths of both the rhodium(I) and rhodium(III) β-diketonatobis(triphenylphosphite) complexes. The relative stability of the four possible rhodium(III) reaction products is the same for all the models, with trans-[Rh(acac)(P(OPh)₃)₂(CH₃)(I)] (in agreement with experimental data) as the most stable reaction product. The best agreement between the theoretical and experimental activation parameters was obtained for the full experimental system.     

Effects of multidimensional tunneling in the kinetics of hydrogen abstraction reactions of O (3P) with CH3OCHO

Quantum tunneling paths are important in reactions when there is a significant component of hydrogenic motion along the potential energy surface. In this study, variational transition state with multidimensional tunneling corrections are employed in the calculations of the thermal rate constants for hydrogen abstraction from the cis‐CH₃OCHO by O (³P) giving CH₃OCO + OH (R1) and CH₂OCHO + OH (R2). The structures and electronic energies are computed with the M06‐2X method. Benchmark calculations with the CBS_D–T approach give an enthalpy of reaction at 0 K for R1 (−2.8 kcal/mol) and R2 (−2.5 kcal/mol) which are in good agreement with the experiment, i.e. −2.61 and −1.81 kcal/mol. At the low and intermediate values of temperatures, small‐ and large‐curvature tunneling dominate the kinetics of R1, which is the dominant path over the range of temperature from 250 to 1200 K. This study shows the importance of multidimensional tunneling corrections for both R1 and R2, for which the total rate constant at 298 K calculated with the CVT/μOMT method is 8.2 × 10⁻¹⁵ cm³ molecule⁻¹ s⁻¹ which agrees well with experiment value of 9.3 × 10⁻¹⁵ cm³ molecule⁻¹ s⁻¹ (Mori, Bull. Inst. Chem. Res. 1981, 59, 116). © 2018 Wiley Periodicals, Inc.