Potential energy surface and reaction mechanism for the ion-molecule reaction: CH4 + CH4+ → CH3 + CH5+

The potential energy surface for the CH₄+CH₄⁺ reaction system has been calculated with the ab initio method. A stable complex, responsible for the complex mechanism, has been found but is hard to reach. Each of the two direct mechanisms, hydrogen transfer and proton transfer, has been shown to consist of a combination of electron transfer and hydrogen atom transfer processes.     

Quasi-classical trajectory study of the dynamics of the Cl + CH4→ HCl + CH3 reaction

We present an on-the-fly classical trajectory study of the Cl + CH(4)→ HCl + CH(3) reaction using a specific reaction parameter (SRP) AM1 Hamiltonian that was previously optimized for the Cl + ethane reaction [S. J. Greaves et al., J. Phys Chem A, 2008, 112, 9387]. The SRP-AM1 Hamiltonian is shown to be a good model for the potential energy surface of the title reaction. Calculated differential cross sections, obtained from trajectories propagated with the SRP-AM1 Hamiltonian compare favourably with experimental results for this system. Analysis of the vibrational modes of the methyl radical shows different scattering distributions for ground and vibrationally excited products.     

Mechanistic and kinetic study of CF3CH═CH2 + OH reaction

The potential energy surfaces of the CF(3)CH═CH(2) + OH reaction have been investigated at the BMC-CCSD level based on the geometric parameters optimized at the MP2/6-311++G(d,p) level. Various possible H (or F)-abstraction and addition/elimination pathways are considered. Temperature- and pressure-dependent rate constants have been determined using Rice-Ramsperger-Kassel-Marcus theory with tunneling correction. It is shown that IM1 (CF(3)CHCH(2)OH) and IM2 (CF(3)CHOHCH(2)) formed by collisional stabilization are major products at 100 Torr pressure of Ar and in the temperature range of T < 700 K (at P = 700 Torr with N(2) as bath gas, T ≤ 900 K), whereas CH(2)═CHOH and CF(3) produced by the addition/elimination pathway are the dominant end products at 700-2000 K. The production of CF(3)CHCH and CF(3)CCH(2) produced by hydrogen abstractions become important at T ≥ 2000 K. The calculated results are in good agreement with available experimental data. The present theoretical study is helpful for the understanding the characteristics of the reaction of CF(3)CH═CH(2) + OH.     

MCSCF and VTST studies of reaction CH + H_2→CH_3

The reaction path of CH+H2→CH3 is traced with Fukui's IRC theory and MCSCF/6-31G method. On this basis. the dynamical properties along the reaction path and variational transition state theory, (VTST. including CVT, ICVT,μVT and US methods) rate constants are investigated. The results show that the recrossing effect is small, but the curvature effect on rate constant is considerable and it must be taken into account in theoretical study.     

CF3CH3 --> HF + CF2CH2: a non-RRKM reaction?

The experimental shock tube data recently reported by Kiefer et al. [J. Phys. Chem. A 2004, 108, 2443-2450] for the title reaction at temperatures between 1600 and 2400 K have been compared to master equation simulations using three models: (a) standard RRKM theory, (b) RRKM theory modified by local random matrix theory, which introduces dynamical corrections arising from slow intramolecular vibrational energy randomization, and (c) an ad hoc empirical non-RRKM model. Only the third model provides a good fit of the Kiefer et al. unimolecular reaction rate data. In separate simulations, all three models accurately reproduce the experimental 300 K chemical activation data of Marcoux and Setser [J. Phys. Chem. 1978, 82, 97-108] when the energy transfer parameters are freely varied to fit the data. When experimental energy transfer parameters for a geometrical isomer (1,1,2-trifluoroethane) are used, the standard RRKM model fits the chemical activation data better than the other models, but if energy transfer in the 1,1,1-trifluoroethane is significantly reduced in comparison to the 1,1,2 isomer, then the empirical ad hoc non-RRKM model also gives a good fit. While the ad hoc empirical non-RRKM model can be made to fit the data, it is not based on theory, and we argue that it is physically unrealistic. We also show that the master equation simulations can mimic the Kiefer et al. vibrational relaxation data, which was the first shock tube observation of double-exponential relaxation. We conclude that, until more data on the trifluoroethanes become available, the current evidence is insufficient to decide with confidence whether non-RRKM effects are important in this reaction, or whether the Kiefer et al. data can be explained in some other way.     

RISM-SCF study of the free-energy profile of the Menshutkin-type reaction NH3+CH3Cl→NH3CH3++Cl− in aqueous solution

The free-energy profile for the Menshutkin-type reaction NH₃ + CH₃Cl → NH₃CH₃ ⁺ + Cl⁻ in aqueous solution is studied using the RISM-SCF method. The effect of electron correlation on the free-energy profile is estimated by the RISM-MP2 method at the HF optimized geometries along the reaction coordinate. Solvation was found to have a large influence on the vibrational frequencies at the reactant, transition state and product; these vibrational frequencies are utilized to calculate the zero-point energy correction of the free-energy profile. The computed barrier height and reaction exothermicity are in reasonable agreement with those of experiment and previous calculations. The change of solvation structure along the reaction path is represented by radial distribution functions between solute-solvent atomic sites. The mechanisms of the reaction are discussed from the view points of solute electronic and solvation structures. Received: 26 June 1998/Accepted: 28 August 1998 / Published online: 2 November 1998     

Crossed-beams studies of the dynamics of the H-atom abstraction reaction, O(3P) + CH4 → OH + CH3, at hyperthermal collision energies

The H-atom abstraction reaction, O((3)P) + CH(4) → OH + CH(3), has been studied at a hyperthermal collision energy of 64 kcal mol(-1) by two crossed-molecular-beams techniques. The OH products were detected with a rotatable mass spectrometer employing electron-impact ionization, and the CH(3) products were detected with the combination of resonance-enhanced multiphoton ionization (REMPI) and time-sliced ion velocity-map imaging. The OH products are mainly formed through a stripping mechanism, in which the reagent O atom approaches the CH(4) molecule at large impact parameters and the OH product is scattered in the forward direction: roughly the same direction as the reagent O atoms. Most of the available energy is partitioned into product translation. The dominance of the stripping mechanism is a unique feature of such H-atom abstraction reactions at hyperthermal collision energies. In the hyperthermal reaction of O((3)P) with CH(4), the H-atom abstraction reaction pathway accounts for 70% of the reactive collisions, while the H-atom elimination pathway to produce OCH(3) + H accounts for the other 30%.     

Semiempirical calculation of the CH4 + CH3 · → CH3 · + CH4 radical-reaction potential surfaces in a basis of frontier orbitals

A semiempirical approach is suggested to describe potential-energy surfaces (PESs) of some radical reactions in the zero-differential overlap (ZDO) approximation. An incomplete basis set is used including only frontier (single-filled) radical molecular orbitals (MOs) depending on geometrical parameters. All possible configurations are taken into account. The parameter selection techniques are analyzed. The approach is applied to PES calculation of the CH₄ + CH₃ · CH₃ · + CH₄ reaction. Some points of the PES are verified by a nonempirical method using the perturbation theory and taking into account the correlation energy. The relaxation energies are calculated. The one-center parameters are determined nonempirically from the CH₃ + CH₃ ·, and CH₃ ^– energetics. The two-center parameters are found by modeling the CH₄ CH₃ · + H and C₂H₆ 2CH₃ reactions in the same single-orbital approximation. The energy parameters of the reactions considered are overestimated by 10%, whereas the geometrical parameters are under-estimated by 15%. Further, a comparative analysis of the Hartree-Fock solutions and those including correlation interactions (CIs) is given. The variations in the spin and charge densities on the reaction centers are considered.Institute of Chemical Kinetics and Combustion, Academy of Sciences of the USSR, Siberian Branch, Novosibirsk. Translated from Teoreticheskaya i Éksperimental'naya Khimiya, Vol. 27, No. 4, pp. 499–506, July–August, 1991. Original article submitted February 10, 1989.     

Ab initio potential energy surface and quantum dynamics for the H + CH4 → H2 + CH3 reaction

A new full-dimensional potential energy surface for the title reaction has been constructed using the modified Shepard interpolation scheme. Energies and derivatives were calculated using the UCCSD(T) method with aug-cc-pVTZ and 6-311++G(3df,2pd) basis sets, respectively. A total number of 30,000 data points were selected from a huge number of molecular configurations sampled by trajectory method. Quantum dynamical calculations showed that the potential energy surface is well converged for the number of data points for collision energy up to 2.5 eV. Total reaction probabilities and integral cross sections were calculated on the present surface, as well as on the ZBB3 and EG-2008 surfaces for the title reaction. Satisfactory agreements were achieved between the present and the ZBB3 potential energy surfaces, indicating we are approaching the final stage to obtain a global potential energy surface of quantitative accuracy for this benchmark polyatomic system. Our calculations also showed that the EG-2008 surface is less accurate than the present and ZBB3 surfaces, particularly in high energy region.     

A photoionization study of the ion-neutral complexes CH3CH +CH 3 · CH2CH3] and CH3CH2CH+CH 3 · CH3 in the gas phase: Formation,H-transfer and C—C bond formation between partners,and channeling of energy into dissociation

Photoionization mass spectrometry was used to investigate the dynamics of ion-neutral complex-mediated dissociations of the n-pentane ion (1). Reinterpretation of previous data demonstrates that a fraction of ions 1 isomerizes to the 2-methylbutane ion (2) through the complex CH₃CH⁺CH ₃ ^· CH₂CH₃ (3), but not through CH₃CH⁺CH₂CH ₃ ^· CH₃ (4). The appearance energy for C₃Hin ₇ ⁺ formation from 1 is 66 kJ mol^?1 below that expected for the formation of n-C₃H ₇ ⁺ and just above that expected for formation of i-C₃H ₇ ⁺ . This demonstrates that the H shift that isomerizes C₃H ₇ ⁺ is synchronized with bond cleavage at the threshold for dissociation to that product. It is suggested that ions that contain n-alkyl chains generally dissociate directly to more stable rearranged carbenium ions. Ethane elimination from 3 is estimated to be about seven times more frequent than is C-C bond formation between the partners in that complex to form 2, which demonstrates a substantial preference in 3 for H abstraction over C-C bond formation. In 1 → CH₃CH⁺CH₂CH₃ + CH₃ by direct cleavage of the C1–C2 bond, the fragments part rapidly enough to prevent any reaction between them. However, 1 → 2 → 4 → C4H ₈ ⁺ + CH₄ occurs in this same energy range. Thus some of the potential energy made available by the isomerization of n-C₄H₉ in 1 is specifically channeled into the coordinate for dissociation. In contrast, analogous formation of 3 by 1 → 3 is predominantly followed by reaction between the electrostatically bound partners.     

Bimolecular reaction rates from ring polymer molecular dynamics: application to H + CH4 → H2 + CH3

In a recent paper, we have developed an efficient implementation of the ring polymer molecular dynamics (RPMD) method for calculating bimolecular chemical reaction rates in the gas phase, and illustrated it with applications to some benchmark atom-diatom reactions. In this paper, we show that the same methodology can readily be used to treat more complex polyatomic reactions in their full dimensionality, such as the hydrogen abstraction reaction from methane, H + CH(4) → H(2) + CH(3). The present calculations were carried out using a modified and recalibrated version of the Jordan-Gilbert potential energy surface. The thermal rate coefficients obtained between 200 and 2000 K are presented and compared with previous results for the same potential energy surface. Throughout the temperature range that is available for comparison, the RPMD approximation gives better agreement with accurate quantum mechanical (multiconfigurational time-dependent Hartree) calculations than do either the centroid density version of quantum transition state theory (QTST) or the quantum instanton (QI) model. The RPMD rate coefficients are within a factor of 2 of the exact quantum mechanical rate coefficients at temperatures in the deep tunneling regime. These results indicate that our previous assessment of the accuracy of the RPMD approximation for atom-diatom reactions remains valid for more complex polyatomic reactions. They also suggest that the sensitivity of the QTST and QI rate coefficients to the choice of the transition state dividing surface becomes more of an issue as the dimensionality of the reaction increases.     

Quantum dynamics of heavy light heavy reactions: Application to (F + CH4 → FCH3 + H) reaction

We apply the time-dependent theory to the collinear exchange reaction F + CH₄ → FCH₃ + H. We have performed detailed calculations on two-dimensional potential surfaces representing the ground electronic potentials of the collinear F + CH₄ → FCH₃ + H reaction, at incident energies. Transmission coefficients range from zero to unity, depending upon the incident energy. Normal modes of vibrations are displayed along the reaction path.     

A transition state view on reactive scattering: initial state-selected reaction probabilities for the H + CH4 → H2 + CH3 reaction studied in full dimensionality

Initial state-selected reaction probabilities for the H+CH(4)→H(2)+CH(3) reaction are computed for vanishing total angular momentum by full-dimensional calculations employing the multiconfigurational time-dependent Hartree approach. An ensemble of wave packets completely describing reactivity for total energies up to 0.58 eV is constructed in the transition state region by diagonalization of the thermal flux operator. These wave packets are then propagated into the reactant asymptotic region to obtain the initial state-selected reaction probabilities. Reaction probabilities for reactants in all rotational states of the vibrational 1A(1), 1F(2), and 1E levels of methane are presented. Vibrational excitation is found to decrease reactivity when reaction probabilities at equivalent total energies are compared but to increase reaction probabilities when the comparison is done at the basis of equivalent collision energies. Only a fraction of the initial vibrational energy can be utilized to promote the reaction. The effect of rotational excitation on the reactivity differs depending on the initial vibrational state of methane. For the 1A(1) and 1F(2) vibrational states of methane, rotational excitation decreases the reaction probability even when comparing reaction probabilities at equivalent collision energies. In contrast, rotational energy is even more efficient than translational energy in increasing the reaction probability when the reaction starts from the 1E vibrational state of methane. All findings can be explained employing a transition state based interpretation of the reaction process.     

Ten-Dimensional Quantum Dynamics Study of H+CH3D→H2+CH2D Reaction

The mode selectivity of the H+CH3D→H2+CH2D reaction was studied using a recently developed ten-dimensional time-dependent wave packet method.The reac-tion dynam...     

Reduced dimensionality spin-orbit dynamics of CH3 + HCl ⇌ CH4 + Cl on ab initio surfaces

A reduced dimensionality quantum scattering method is extended to the study of spin-orbit nonadiabatic transitions in the CH(3) + HCl ? CH(4) + Cl((2)P(J)) reaction. Three two-dimensional potential energy surfaces are developed by fitting a 29 parameter double-Morse function to CCSD(T)/IB//MP2/cc-pV(T+d)Z-dk ab initio data; interaction between surfaces is described by geometry-dependent spin-orbit coupling functions fit to MCSCF/cc-pV(T+d)Z-dk ab initio data. Spectator modes are treated adiabatically via inclusion of curvilinear projected frequencies. The total scattering wave function is expanded in a vibronic basis set and close-coupled equations are solved via R-matrix propagation. Ground state thermal rate constants for forward and reverse reactions agree well with experiment. Multi-surface reaction probabilities, integral cross sections, and initial-state selected branching ratios all highlight the importance of vibrational energy in mediating nonadiabatic transition. Electronically excited state dynamics are seen to play a small but significant role as consistent with experimental conclusions.     

F-+CH3Cl→CH3F+Cl-反应过程中的分子形貌变化

双分子亲核(SN2)反应是重要的基本有机反应之一,其中电子从亲核基团向离去基团的转移发挥着关键作用.利用从头计算方法CCSD(T)/aug-cc-pVDZ和我们发展的分子形貌理论,对反应F-+CH3Cl→CH3F+Cl-进行了研究,给出了反应过程中分子形状和电子转移的动态变化图像.结果表明,沿内禀反应坐标,从反应开始到生成反应前复合物,亲核试剂F-的分子内禀特征轮廓在缓慢收缩,而其上的电子密度在缓慢增大.此后,F的轮廓迅速膨胀,电子密度急剧下降,尤其是从过渡态到产物复合物的过程中.而在反应过程中,离去基团Cl的轮廓一直在收缩,其上的电子密度一直在增大.对反应过程中电子所受到作用势的研究表明,随着反应的进行,电子在F与C间受到的作用势逐渐降低,而在C与Cl间受到的作用势逐渐升高,清楚地展现反应过程中F与C间化学键生成和C与Cl间化学键断裂的动态过程.     

Ab initio study of the reaction path of the reaction HNCO+CN

HNCO is a convenient photolytic source of NCO and NH radicals for laboratory kinetics studies of elementary reaction[1] and plays an important role in the combustion and atmosphere chemistry. It can re- move deleterious compounds rapidly from exhausted ga…     

CH3CH2+O(3P)反应机理的理论研究

The reaction for CH3CH2+O(3P) was studied by ab initio method. The geometries of the reactants, intermediates, transition states and products were optimized at MP2/6-311+G(d,p) level. The corresponding vibration frequencies were calculated at the same level. The single-point calculations for all the stationary points were carried out at the QCISD(T)/6-311+G(d,p) level using the MP2/6-311+G(d,p) optimized geometries. The results of the theoretical study indicate that the major products are the CH2O＋CH3, CH3CHO+H and CH2CH2+OH in the reaction. For the products CH2O＋CH3 and CH3CHO+H, the major production channels are A1: (R)→IM1→TS3→(A) and B1: (R)→IM1→TS4→(B), respectively. The majority of the products CH2CH2+OH are formed via the direct abstraction channels C1 and C2: (R)→TS1(TS2)→(C). In addition, the results suggest that the barrier heights to form the CO reaction channels are very high, so the CO is not a major product in the reaction.     

Electronic state dependence of the ion-molecule reaction CH(3)CN(+) + CH(3)CN → CH(4)CN(+) + CH(2)CN: threshold electron-secondary ion coincidence (TESICO) and direct ab initio molecular dynamics study

The ion-molecule reaction, CH(3)CN(+) + CH(3)CN → CH(3)CNH(+) + CH(2)CN, has been investigated using the threshold electron-secondary ion coincidence (TESICO) technique. Relative reaction cross sections for two microscopic reaction mechanisms, i.e., proton transfer (PT) from the acetonitrile ion CH(3)CN(+) to neutral acetonitrile CH(3)CN and hydrogen atom abstraction (HA) by CH(3)CN(+) from CH(3)CN, have been determined for two low-lying electronic states, (2)E and (2)A(1) of the CH(3)CN(+) primary ion. The cross section for PT of the (2)A(1) state was smaller than that of the (2)E state, whereas that of HA are almost the same in the two states. Ab initio calculations showed that the dissociation of the C-H(+) bond of CH(3)CN(+) is easier in the (2)E state than that in the (2)A(1) state. The direct ab initio molecular dynamics (MD) calculations showed that two mechanisms, direct proton transfer and complex formation, contribute the reaction dynamics.