Ab initio and direct quasiclassical-trajectory study of the F+CH4-->HF+CH3 reaction

We present an electronic structure and dynamics study of the F+CH4-->HF+CH3 reaction. CCSD(T)/aug-cc-pVDZ geometry optimizations, harmonic-frequency, and energy calculations indicate that the potential-energy surface is remarkably isotropic near the transition state. In addition, while the saddle-point F-H-C angle is 180 degrees using MP2 methods, CCSD(T) geometry optimizations predict a bent transition state, with a 153 degrees F-H-C angle. We use these high-quality ab initio data to reparametrize the parameter-model 3 (PM3) semiempirical Hamiltonian so that calculations with the improved Hamiltonian and employing restricted open-shell wave functions agree with the higher accuracy data. Using this specific-reaction-parameter PM3 semiempirical Hamiltonian (SRP-PM3), we investigate the reaction dynamics by propagating quasiclassical trajectories. The results of our calculations using the SRP-PM3 Hamiltonian are compared with experiments and with the estimates of two recently reported potential-energy surfaces. The trajectory calculations using the SRP-PM3 Hamiltonian reproduce quantitatively the measured HF vibrational distributions. The calculations also agree with the experimental HF rotational distributions and capture the essential features of the excitation function. The results of the SRP semiempirical Hamiltonian developed here clearly improve over those using the two prior potential-energy surfaces and suggest that reparametrization of semiempirical Hamiltonians is a promising strategy to develop accurate potential-energy surfaces for reaction dynamics studies of polyatomic systems.     

Ab initio study on the mechanism of reaction HNCO+NH_2

Ab initio UMP2 and UQCISD(T) calculations, with 6-311G** basis sets, were performed for the titled reactions. The results show that the reactions have two product channels: NH2+ HNCO→NH3+NCO (1) and NH2+HNCO-N2H3+CO (2), where reaction (1) is a hydrogen abstraction reaction via an H-bonded complex (HBC), lowering the energy by 32.48 kJ/mol relative to reactants. The calculated QCISD(T)//MP2(full) energy barrier is 29.04 kJ/mol, which is in excellent accordance with the experimental value of 29.09 kJ/mol. In the range of reaction temperature 2300-2700 K, transition theory rate constant for reaction (1) is 1.68 × 1011- 3.29 × 1011 mL · mol-1· s-1, which is close to the experimental one of 5.0 ×1011 mL× mol-1· s-1 or less. However, reaction (2) is a stepwise reaction proceeding via two orientation modes, cis and trans, and the energy barriers for the rate-control step at our best calculations are 92.79 kJ/mol (for cis-mode) and 147.43 kJ/mol (for trans-mode), respectively, which is much higher than     

自由基-分子反应:F+Propene(CH2CHCH3)的从头算研究

The radical-molecule reaction F+propene (CH2CHCH3) was studied in detail by using the Becke's three parameter Lee-Yang-Parr-B3LYP/6-311G(d,p) and coupled cluster with single, double, and triple excitationsCCSD(T)/6-311+G(2d,2p). It is shown that F+propene reaction mainly occurs through complex-formation mechanism: F attacks the double bond of propene leading to the formation of complex 1 and complex 2. As the two radical complexes are metastable, they can quickly dissociate to H+C3H5F, CH3+C2H3F and HF+C3H5. Based on the ab initio calculations, the CH3+C2H3F is the main channel, and the H elimination and HF forming channels also provide some contribution to products. The calculated values are in good agreement with the recently reported experimental results.     

Ab initio and direct quasiclassical-trajectory study of the Cl + CH4-->HCl + CH3 reaction

We present an electronic structure and dynamics study of the Cl + CH(4)--> HCl + CH(3) reaction. We have characterized the stationary points of the ground-state potential-energy surface using various electronic structure methods and basis sets. Our best calculations, CCSD(T) extrapolated to the complete basis-set limit based on geometries and harmonic frequencies obtained at the CCSD(T)/aug-cc-pvtz level, are in agreement with the experimental reaction energy and indirect measurements of the barrier height. Using ab initio information, we have reparametrized a semiempirical Hamiltonian so that the predictions of the improved Hamiltonian agree with the higher-level calculations in various regions of the potential-energy surface. This improved semiempirical Hamiltonian is then used to propagate quasiclassical trajectories and characterize the reaction dynamics. The good agreement of the calculated HCl rotational and angular distributions with the experiment indicates that reparametrizing semiempirical Hamiltonians is a promising approach to derive accurate potential-energy surfaces for polyatomic reactions. However, excessive energy leakage from the initial vibrational energy of the CH(4) molecule to the reaction coordinate in the trajectory calculations calls into question the suitability of the standard quasiclassical-trajectory method to describe energy partitioning in polyatomic reactions.     

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     

Ab initio study of the formation of C3H3+ from the reaction of CH3+ with acetylene

Ab initio molecular orbital theory has been used to study the mechanism of the formation of C₃H₃⁺ from the reaction of CH₃⁺ with acetylene. The highest level geometry optimizations and frequencies were computed at MP2-FC/6-31G**; single point energies of all the critical structures were computed to the MP4-FC/6-31G**//MP2-FC/6-31G** theory level. One of the three alternative transition structures leading to the formation of C₃H₃⁺ gives the cyclopropenyl cation and the other two the propargyl cation. The proportions of C₃H₂D⁺ and C₃HD₂⁺ obtained when CD₃⁺ reacts with acetylene, and the composite nature of the metastable peak observed for the [C₃H₅]⁺→[C₃H₃]⁺ + H₂ fragmentation are explained by assuming a different degree of deuterium scrambling depending on the energy of the system. © 1996 by John Wiley & Sons, Inc.     

Ab initio characterization of (CH3IO3) isomers and the CH3O2 + IO reaction pathways

The geometries, harmonic vibrational frequencies, relative energetics, and enthalpies of formation of (CH(3)IO(3)) isomers and the reaction CH(3)O(2) + IO have been investigated using quantum mechanical methods. Optimization has been performed at the MP2 level of theory, using all electron and effective core potential, ECP, computational techniques. The relative energetics has been studied by single-point calculations at the CCSD(T) level. Methyl iodate, CH(3)OIO(2), is found to be the lowest-energy isomer showing particular stabilization. The two nascent association minima, CH(3)OOOI and CH(3)OOIO, show similar stabilities, and they are considerably higher located than CH(3)OIO(2). Interisomerization barriers have been determined, along with the transition states involved in various pathways of the reaction CH(3)O(2) + IO.     

Ab initio study of the spectroscopy of CH3N and CH3CH2N

Complete active space (CAS) calculations with 6-311++g(3df,3pd) basis sets were performed for a large number of electronic states of the nitrate free radical (CH3N/CH3CH2N) and their positive and negative ions. All calculated states are valence states, and their characters are discussed in detail. To investigate the Jahn-Teller effect on the CH3N radical, Cs symmetry was used for both CH3N and CH3CH2N in calculations. The results (CASPT2 adiabatic excitation energies and CASSI oscillator strengths) suggest that the calculated transitions of CH3N at 32172 and 32139 cm(-1) are attributed to the 2(3)A' ' --> 1(3)A' ' and 1(3)A' --> 1(3)A' ', respectively, which is in accordance with the A3E --> X3A2 emission spectrum at T0 = 31 817 cm(-1). The calculated transitions of CH3CH2N at 334 nm are attributed to the 1(3)A' ' --> 2(3)A' ' and 1(3)A' ' --> 1(3)A', respectively, which is in accordance with the UV absorption spectrum of a series of 11 bands beginning at 335 nm. The vertical and adiabatic ionization energies were obtained to compare with the PES data. These results are in agreement with previous experimental data, which is discussed in detail.     

Theoretical study on the mechanism of the CH2F + NO2 reaction

Despite the importance of the Fluoromethyl radicals in combustion chemistry, very little experimental information on their reactions toward stable molecules is available in the literature. Motivated by recent laboratory characterization about the reaction kinetics of Chloromethyl radicals with NO2, we carried out a detailed potential energy survey on the CH2F + NO2 reaction at the B3LYP/6-311G(d,p) and MC-QCISD (single-point) levels as an attempt toward understanding the CH2F + NO2 reaction mechanism. It is shown that the CH2F radical can react with NO2 to barrierlessly generate adduct a (H2FCNO2), followed by isomerization to b1 (H2FCONO-trans) which can easily interconvert to b2 (H2FCONO-cis). Subsequently, Starting from b (b1, b2), the most feasible pathway is the C--F and N--O1 bonds cleavage along with N--F bond formation of b (b1, b2) leading to P1 (CH2O + FNO), or the direct N--O1 weak-bond fission of b (b1, b2) to give P2 (CH2FO + NO), or the 1,3-H-shift associated with N--O1 bond rupture of b1 to form P3 (CHFO + HNO), all of which may have comparable contribution to the reaction CH2F + NO2. Much less competitively, b2 either take the 1,4-H-shift and O1--N bond cleavage to form product P4 (CHFO + HON) or undergo a concerted H-shift to isomer c2 (HFCONOH), followed by dissociation to P4. Because the rate-determining transition state (TSab1) in the most competitive channels is only 0.3 kcal/mol higher than the reactants in energy, the CH2F + NO2 reaction is expected to be rapid, and may thus be expected to significantly contribute to elimination of nitrogen dioxide pollutants. The similarities and discrepancies among the CH2X + NO2 (X = H, F, and Cl) reactions are discussed in terms of the electronegativity of halogen atom. The present article may assist in future experimental identification of the product distributions for the title reaction, and may be helpful for understanding the halogenated methyl chemistry.     

Ab initio study of the mechanism of the atmospheric reaction: NO2 + O3 --> NO3 + O2

The atmospheric reaction NO2 + O3 --> NO3 + O2 (1) has been investigated theoretically by using the MP2, G2, G2Q, QCISD, QCISD(T), CCSD(T), CASSCF, and CASPT2 methods with various basis sets. The results show that the reaction pathway can be divided in two different parts at the MP2 level of theory. At this level, the mechanism proceeds along two transition states (TS1 and TS2) separated by an intermediate, designated as A. However, when the single-reference higher correlated QCISD methodology has been employed, the minimum A and the transition state TS2 are not found on the hypersurface of potential energy, which confirms a direct reaction mechanism. Single-reference high correlated and multiconfigurational methods consistently predict the barrier height of reaction (1) to be within the range 2.5-6.1 kcal mol(-1), in reasonable agreement with experimental data. The calculated reaction enthalpy is -24.6 kcal mol(-1) and the reaction rate calculated at the highest CASPT2 level, of k = 6.9 x 10(-18) cm(3) molecule(-1) s(-1). Both results can be regarded also as accurate predictions of the methodology employed in this article.     

Ab initio study on the reaction of Sc+ + CH4 → Sc+(SINGLE BOND)CH2 + H2

An ab initio study on the reaction of the ground state (³D) and the excited state (¹D) of Sc⁺ with methane was performed. Reaction channels on the singlet and triplet potential surface (PES) and the reaction mechanism are examined and discussed. Three regions of the potential surface was studied: the molecular complex, the C(SINGLE BOND)H insertion products, and the transition states for the reaction. Comparisons between singlet and triplet PESs show that the excited state (¹D) of Sc⁺ has more reactivity with methane than does the ground state (³D) due to the spin quantum number conservation with the more stable insertion intermediate. © 1997 John Wiley & Sons, Inc.     

Ab initio potential energy surface of CH and reaction dynamics of H + CH+

This work presents a new ground state potential energy surface (PES) for CH. The potential is tested using quasi classical trajectory (QCT) and quantum reactive scattering methods for the H + CH(+) reaction. Cross sections and rate coefficients for all reaction channels up to 300 K are calculated. The abstraction rate coefficients follow the expected slightly decreasing behaviour above 90 K, but have a positive gradient with lower temperatures. The inelastic collision and exchange reaction rate constants are increasing monotonically with temperature. The rate coefficients of the exchange reaction differ significantly between QCT and quantum reactive scattering, due to intrinsic shortcomings of the QCT final state distributions.     

DFT and ab initio study on the reaction mechanism of CH2SH + O2

The mechanism for the CH₂SH + O₂ reaction was investigated by DFT and ab initio chemistry methods. The geometries of all possible stationary points were optimized at the B3LYP/6-311+G(d,p) level, and the single point energy was calculated at the CCSD(T)/cc-pVXZ(X = D and T), G3MP2 and BMC-CCSD levels. The results indicate that the oxidation of CH₂SH by O₂ to form HSCH₂OO is a barrierless process. The most favorable channel is the rearrangement of the initial adduct HSCH₂OO (IM1) to form another intermediate H₂C(S)OOH (IM3) via a five-center transition state, and then the C–O bond fission in IM3 leads to a complex CH₂S. . .HO₂ (MC1), which finally gives out to the major product CH₂S + HO₂. Due to high barriers, other products including cis- and trans-HC(O)SH + HO could be negligible. The direct abstraction channel was also determined to yield CH₂S + HO₂, with the barrier height of 22.3, 18.1 and 15.0 kcal/mol at G3MP2, CCSD(T)/cc-pVTZ and BMC-CCSD levels, respectively, it is not competitive with the addition channel, in which all stationary points are lower than reactant energetically. The other channels to produce cis- and trans-CHSH + HO₂ are also of no importance.     

Ab initio studies on reaction H2C=C(OH)Li+CH3 + (CH3 -)

The possible structures and isomerizations of H₂C=C(OH)Li are studied theoretically by the gradient analytical method at RHF/6-31+G level. According to these results, reactions of H₂C=C(OH)Li with CH₃ ⁺ and CH ₃ ^- are investigated thoroughly. When H₂C=C(OH)Li reacts with CH ₃ ⁺ , H_zC=C(OH)Li firstly changes from structure1 to structure4, and then combines with CH₃ ⁺. In this reaction, the configuration of central carbon is retained. When H₂C=C(OH)Li reacts with CH ₃ ^- , structure1 firstly breaks its C-O bond to give contact ion-pair. Then through transition state16 which is similar to structure2, the attack of CH ₃ ^- from the opposite side of^-OH replaces^-OH group and inverts the configuration of carbenoid carbon atom. All the results show that the ambident reactivity of carbenoid has close relationship with the stability of special structures. Project supported by the National Natural Science Foundation of China (Grant No. 29773025).     

Ab initio方法研究CH~3+OClO反应的可能通道

利用abinitio方法研究了CH~3+OClO反应的三个可能通道,首次应用UMP2(full)/6-31G(d,p)方法得到各反应物、产物、中间物及过渡态的优化构型和谐振频率;然后采用G2MP2理论计算各通道反应焓变和势垒高度。理论计算表明产物通道CH~2O+HOCl是最可能发生的途径,反应放热为443.80kJ·mol^-^1。可能的反应过程为：CH~3和OClO自由基先经无垒过程生成了一个富能中间物,继而通过较低的势垒解离成HOCl+H~2CO。     

Ab initio quantum chemical studies of reaction mechanism for CN with CH2CO

Six product channels have been found in the association reaction of CN + CH₂CO, and a variety of possible complexes and saddle points along the minimum energy reaction paths have been characterized at the UMP2(full)/6‐31G(d) level. The dominant reaction channels are the production of CH₂CN + CO and CH₂NC + CO. The isomerization and dissociation reactions of the major products of CH₂CN and CH₂NC have been investigated using the G2MP2 level. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

用从头算研究Shilov反应的机理

Shilov反应在CH~4活化中占有中心地位,它有氧化加成和σ迁移两种可能的机理。本文用较大基组的从头算研究了这两种机理的反应过程,认为Shilov反应应按氧化加成机理进行。     

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.     

Ab initio/density functional theory and multichannel RRKM study for the ClO + CH2O reaction

The potential energy surface for the CH(2)O + ClO reaction was calculated at the QCISD(T)/6-311G(2d,2p)//B3LYP/6-311G(d,p) level of theory. The rate constants for the lower barrier reaction channels producing HOCl + HCO, H atom, OCH(2)OCl, cis-HC(O)OCl and trans-HC(O)OCl have been calculated by TST and multichannel RRKM theory. Over the temperature range of 200-2000 K, the overall rate constants were k(200-2000K) = 1.19 x 10(-13)T(0.79) exp(-3000.00/T). At 250 K, the calculated overall rate constant was 5.80 x 10(-17) cm(3) molecule(-1) s(-1), which was in good agreement with the experimental upper limit data. The calculated results demonstrated that the formation of HOCl + HCO was the dominant reaction channel and was exothermic by 9.7 kcal/mol with a barrier of 5.0 kcal/mol. When it retrograded to the reactants CH(2)O + ClO, an energy barrier of 14.7 kcal/mol is required. Furthermore, when HOCl decomposed into H + ClO, the energy required was 93.3 kcal/mol. These results suggest that the decomposition in both the forward and backward directions for HOCl would be difficult in the ground electronic state.     

Computational study on the reaction CH2CH2 + F --> CH2CHF + H

Previous ab initio studies on reactions involving radical addition to alkenes showed that such reactions are very sensitive to theoretical levels, and thus are difficult to deal with. This motivates us to theoretically reexamine the title reaction thoroughly, which has been studied only at several low levels of theory. In the present work, the geometry optimizations and energy calculations for all species involved in the title reaction were performed at several high levels of theory. The reaction mechanism of the title reaction is discussed at the CCSD(T)/aug-cc-pVDZ//CCSD/6-31G(d,p) theoretical level. According to our study, the fluorine addition to ethylene occurs via the formation of a prereaction complex with C2v symmetry, which is pointed out for the first time. The prereaction complex evolves into a fluoroethyl radical almost without a barrier, with an exothermicity of 41.49 kcal/mol. The fluoroethyl radical can further decompose into a hydrogen atom and fluoroethylene, with an energy release of 10.33 kcal/mol. Besides the direct departure of the hydrogen atom from the fluoroethyl radical, an indirect decomposition pathway may also be open, which has not been reported before. In addition, the formation of a fluoroethyl radical from a separate fluorine atom and ethylene is described pictorially via the molecular intrinsic characteristic contour (MICC) and the electron density mapped on it. Thereby, strong interpolarization and evident electron transfer between the fluorine atom and ethylene are observed as they approach each other. The transition structure for the fluorine addition to ethylene is clearly shown to be reactant-like. This provides new and intuitional insight into the title reaction.