Ab initio chemical kinetics for the OH+HNCN reaction

The kinetics and mechanism of the reaction of the cyanomidyl radical (HNCN) with the hydroxyl radical (OH) have been investigated by ab initio calculations with rate constants prediction. The single and triplet potential energy surfaces of this reaction have been calculated by single-point calculations at the CCSD(T)/6-311+G(3df,2p) level based on geometries optimized at the B3LYP/6-311+G(3df,2p) and CCSD/6-311++G(d,p) levels. The rate constants for various product channels in the temperature range of 300-3000 K are predicted by variational transition-state and Rice-Ramsperger-Kassel-Marcus (RRKM) theories. The predicted total rate constants can be represented by the expressions ktotal=2.66 x 10(+2)xT-4.50 exp(-239/T) in which T=300-1000 K and 1.38x10(-20)xT2.78 exp(1578/T) cm3 molecule(-1) s(-1) where T=1000-3000 K. The branching ratios of primary channels are predicted: k1 for forming singlet HON(H)CN accounts for 0.32-0.28, and k4 for forming singlet HONCNH accounts for 0.68-0.17 in the temperature range of 300-800 K. k2+k7 for producing H2O+NCN accounts for 0.55-0.99 in the high-temperature range of 800-3000 K. The branching ratios of k3 for producing HCN+HNO, k6 for producing H2N+NCO, k8 for forming 3HN(OH)CN, k9 for producing CNOH+3NH, and k5+k10 for producing NH2+NCO are negligible. The rate constants for key individual product channels are provided in a table for different temperature and pressure conditions.     

Theoretical Study of Reaction Mechanism of 1-Propenyl Radical with NO

The reaction system of 1-propenyl radical with NO is an ideal model for studying the intermolecular and intramolecular reactions of complex organic free radicals containing C=C double bonds. On the basis of the full optimization of all species with the Gaussian 98 package at the B3LYP/6-311++G** level, the reaction mechanism was elucidated extensively using the vibrational mode analysis. There are seven reaction pathways and five sets of small molecule end products: CH2O+CH3CN, CH2CHCN+H2O, CH3CHO+HCN, CH3CHO+HNC, and CH3CCH+HNO. The channel of C3H5￠+NO→ IM1→TS1→IM2→TS2→IM3→TS3→CH3CHO+HCN is thermodynamically most favorable.     

CH3SS与OH自由基反应机理的理论研究

应用密度泛函理论(DFT)对CH3SS与OH自由基单重态反应机理进行了研究.在B3PW91/6-311+G(d,p)水平上优化了反应通道上各驻点(反应物、中间体、过渡态和产物)的几何构型,用内禀反应坐标(IRC)计算和频率分析方法对过渡态进行了验证.在QCISD(T)/6-311++G(d,p)水平上计算了各物种的单点能,并对总能量进行了零点能校正.研究结果表明,CH3SS与OH反应为多通道反应,有5条可能的反应通道.反应物首先通过不同的S—O键相互作用形成具有竞争反应机理的中间体IM1和IM2.再经过氢迁移、脱氢和裂解等机理得到主要产物P1(CH2SS+H2O),次要产物P2(CH2S+HSOH),P3(CH3SH+1SO)和P4(CH2SSO+H2),其中最低反应通道的势垒为174.6kJ.mol-1.     

Direct ab initio dynamics calculation of the reaction rates of CH3OCl with OH

A direct ab initio dynamics method was carried out for the reaction CH3OCl + OH --> products. Three abstraction channels from chlorine atom, in-plane hydrogen, and out-of-plane hydrogen atoms at the CH3 group have been found. The optimized geometries and frequencies of the stationary points and the minimum-energy paths (MEPs) were calculated at the MP2/6-311G(d,p) level. To improve the reaction enthalpy and potential barrier, single-point calculations were made at three higher levels of theory, the approximate QCISD(T)/6-311++G(3df,2pd), G3, and G3(MP2) levels. Furthermore, the rate constants for three abstraction channels were evaluated using canonical variational transition state theory (CVT) with the small-curvature tunneling correction (SCT) over a wide temperature range of 220-2000 K at above three higher theory levels, respectively. The calculated rate constants as well as branching rates are in reasonable agreement with the experimental values in the temperature region 250-341 K. The present results indicate H-abstraction especially from out-of-plane hydrogen is the main reaction pathway, while Cl-abstraction is much less competitive.     

硝基甲烷异构化反应势能面的ab initio研究

在B3LYP/6-311+ +G(2d,2p)水平上,优化得到硝基甲烷CH3NO2的10种异构体和23个异构化反应过渡态,并用G2MP2方法进行了单点能计算.根据计算得到的G2MP2相对能量,探讨了CH3NO2势能面上异构化反应的微观机理.研究表明,反应初始阶段的CH3NO2异构化过程具有较高的能垒,其中CH3NO2的两个主要异构化反应通道,即CH3NO2→CH3ONO和CH3NO2→CH2N(O)OH的活化能分别为270.3和267.8 kJ/mol,均高于CH3NO2的C-N键离解能.因而,从动力学角度考虑, CH3NO2的异构化反应较为不利.     

CH3 CH2 +N(4S)反应机理的理论研究

The reaction for CH3CH2+N(4S) 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 CH2CH2+3NH and H2CN＋CH3, and the minor products are the CH3CHN+H in the reaction. The majority of the products CH2CH2+3NH are formed via a direct hydrogen abstraction channel. The products H2CN＋CH3 are produced via an addition/dissociation channel. The products CH3CHN+H are produced via an addition/dissociation channel.     

·OH自由基与CH3CN反应机理及动力学

在CBS-QB3水平上研究了CH3CN 和·OH反应的势能面, 其中包括两个中间体和9个反应过渡态. 分别给出了各主要物质的稳定构型、相对能量及各反应路径的能垒. 根据计算的CBS-QB3势能面, 探讨了CH3CN+·OH反应机理. 计算结果表明, 生成产物P1(·CH2CN+H2O)的反应路径在整个反应体系中占主要地位. 运用过渡态理论对产物通道P1(·CH2CN+H2O)的速率常数k1(cm3·molecule-1·s-1)进行了计算. 预测了k1(cm3·molecule-1·s-1)在250-3000 K温度范围内的速率常数表达式为k1(250-3000 K)=2.06×10^-20T^3.045exp(-780.00/T). 通过与已有的实验值进行对比得出, 在实验所测定的250-320 K 范围内, 计算得到的k1的数值与已有的实验值比较吻合. 由初始反应物生成产物P1 (·CH2CN+H2O)只需要克服一个14.2 kJ·mol-1的能垒. 而产物·CH2CN+H2O生成后要重新回到初始反应物CH3CN+·OH, 则需要克服一个高达111.2 kJ·mol-1的能垒,这就表明一旦产物P1生成后就很难再回到初始反应物.     

CH3NO2+X(X=H,OH, CH3, CH2[^3B1], O[^3P])→CH2NO2+XH吸氢反应过渡态 结 构和反应位垒的DFT计算研究

采用DFT(B3LYP)方法,分别在6-311g(d,p),6-311++g(d,p)和自洽相关基组cc-pVIZ水平上优化了基态硝基甲烷和自由基H,OH,CH3,CH2[^3B1]以及O[^3P]等发生吸氢反应时的过渡态结构,并计算了反应的位垒。研究表明,对同一反应,不同基组下优化得到的过渡态几何结构基本一致;反应位垒数值的大小也基本接近,经校正,硝基甲烷同自由基反应位垒的理论计算值同实验结果基本吻合。     

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.     

On the mechanism of gas-phase reaction of C1-C3 aliphatic thiols + OH radicals

A theoretical study on the mechanism of the OH + aliphatic thiols reactions is presented. Optimum geometries and frequencies have been computed at the BHandHLYP/6-311++G(2d,2p) level of theory for all stationary points. Energy values have been improved by single-point calculations at the above geometries using CCSD(T)/6-311++G(d,p). Twelve possible channels have been modeled, three of them including the possible influence of molecular oxygen, and three of them involving excess of OH. The only channels that have been found to significantly contribute to the overall reaction in the troposphere are the hydrogen abstractions from the -SH group and from the alkyl groups. Our analysis supports a stepwise mechanism involving the formation of a short-lived, weakly bonded adduct in the entrance channel, for the abstraction paths. The results proposed in the present work seem to provide a viable explanation for diverse findings previously reported from experimental investigations.     

3-甲基环状乙撑磷酸二酯醇解的反应途径的理论研究

采用密度泛函理论和MP2方法研究了3-甲基环状乙撑磷酸二酯(MEP)与甲醇的反应途径:(Ⅰ)CH3O-+MEP;(Ⅱ)CH3OH+MEP;(Ⅲ)CH3O-+HMEP(MEP的质子化形式);(Ⅳ)CH3OH+HMEP.在B3LYP/6-31++G(d,p)水平上优化了四条反应途径的反应物、中间体、过渡态及产物的几何构型,并在同水平上进行了自然电荷分析,然后在MP2/6-311++G(3df,2p)水平上计算了各驻点的单点能.采用极化连续介质模型(PCM)研究了各途径在苯、甲醇和水溶液中的溶剂化效应.计算结果表明,溶剂效应使途径(Ⅰ)的自由能垒降低,而使途径(Ⅱ)和(Ⅳ)的决速步骤的自由能垒升高.在气相和苯溶剂中途径(Ⅳ)是反应的优势途径,在甲醇和水溶剂中途径(Ⅰ)则成为最优.研究结果进一步表明实验条件下途径(Ⅱ)与(Ⅳ)对总醇解反应的贡献相当.     

Mechanistic and kinetic study of the CH3CO + O2 reaction

Potential-energy surface of the CH3CO + O2 reaction has been calculated by ab initio quantum chemistry methods. The geometries were optimized using the second-order Moller-Plesset theory (MP2) with the 6-311G(d,p) basis set and the coupled-cluster theory with single and double excitations (CCSD) with the correlation consistent polarized valence double zeta (cc-pVDZ) basis set. The relative energies were calculated using the Gaussian-3 second-order Moller-Plesset theory with the CCSD/cc-pVDZ geometries. Multireference self-consistent-field and MP2 methods were also employed using the 6-311G(d,p) and 6-311++G(3df,2p) basis sets. Both addition/elimination and direct abstraction mechanisms have been investigated. It was revealed that acetylperoxy radical [CH3C(O)OO] is the initial adduct and the formation of OH and alpha-lactone [CH2CO2(1A')] is the only energetically accessible decomposition channel. The other channels, e.g., abstraction, HO2 + CH2CO, O + CH3CO2, CO + CH3O2, and CO2 + CH3O, are negligible. Multichannel Rice-Ramsperger-Kassel-Marcus theory and transition state theory (E-resolved) were employed to calculate the overall and individual rate coefficients and the temperature and pressure dependences. Fairly good agreement between theory and experiments has been obtained without any adjustable parameters. It was concluded that at pressures below 3 Torr, OH and CH2CO2(1A') are the major nascent products of the oxidation of acetyl radicals, although CH2CO2(1A') might either undergo unimolecular decomposition to form the final products of CH2O + CO or react with OH and Cl to generate H2O and HCl. The acetylperoxy radicals formed by collisional stabilization are the major products at the elevated pressures. In atmosphere, the yield of acetylperoxy is nearly unity and the contribution of OH is only marginal.     

Oxygen Atom Exchange Mechanism in Reaction of OH Radical with AsO

IntroductionNitrous acid,HONO,has been extensivelystudied by means of experimental[1,2 ] andtheoretical methods due to its importance inatomspheric chemistry[3 ] . Those studies includedits experimental spectroscopy[1,2 ] and potentialenergy surface with the aid of density functionaltheory[3 ] . Its phosphorus analogue,HOPO,hasbeen studied by virtue of theoreticalcomputations[4] ,and detected in the gas by infraredlaser spectroscopy in2 0 0 0 by Bell and coworkers[5]and matrix isolation in…     

Mechanistic and kinetic study of the O + CH2OH reaction

The mechanism for the O + CH2OH reaction was investigated by various ab initio quantum chemistry methods. For the chemical activation mechanism, that is, the addition/elimination path, the couple-cluster methods including CCSD and CCSD(T) were employed with the cc-pVXZ (X = D, T, Q, 5) basis sets. For the abstraction channels, multireference methods including CASSCF, CASPT2, and MRCISD were used with the cc-pVDZ and cc-pVTZ basis sets. It has been shown that the production of H + HCOOH is the major channel in the chemical activation mechanism. The minor channels include HCO + H2O and OH + CH2O. The hydrogen abstraction by an O atom from the CH2OH radical produces either OH + CH2O or OH + HCOH. Moreover, the two abstraction reactions are essentially barrierless processes. The rate constants for the association of O with CH2OH have been calculated using the flexible transition state theory. A weak negative temperature dependence of the rate constants is found in the range 250-1000 K. Furthermore, it is estimated that the abstraction processes also play an important role in the O + CH2OH reaction. Additionally, the falloff behavior for the OCH2OH --> H + HCOOH reaction has been investigated. The present theoretical results are compared to the experimental measurements to understand the mechanism and kinetic behavior of the O + CH2OH reaction and the unimolecular reaction of the OCH2OH radical.     

Theoretical Investigation of CH3CF2O2+HOO Reaction

The reactions of CH3CF2O2 with HOO are important chemical cyclic processes of photochemical contamination. In this paper, the reaction pathways and reaction mechanism of CH3CF2O2+HOO are investigated extensively with the Gaussian 98 package at the B3LYP/6-311++G** basis sets. The use of vibrational mode analysis and electron population analysis to reveal the reaction mechanism is firstly reported. The study shows that CH3CF2CO2+HOO→IM1→TS1→CH3CF2O2H+O2 channel is the energetically most favorable, CH3CF2CO2H and O2 are the principal products, and the formation of CH3OH and CF2O is also possible.     

CH2=CHCl与O(3P)反应的理论研究

用量子化学密度泛函理论和QCISD(Quadratic configuration interaction calculation)方法,对0(^3P)与CH2CHCl的反应进行了理论研究．在UB3LYP／6—311 G(d,p),UB3LYP／6—31 (3df,3pd)计算水平上,优化了反应物、产物、中间体和过渡态的几何构型,并在UQCISD(T)／6—311 G(2df,2pO)水平上计算了单点能量．为了确证过渡态的真实性,在UB3LYP／6—311 G(3df,3pd)水平上进行了内禀坐标(IRC)计算和频率分析,并确定了反应机理．研究结果表明,反应主要产物为CH2CHO和Cl．     

The OH* + CH3SH reaction: support for an addition-elimination mechanism from ab initio calculations

Several intermediates for the CH(3)SH + OH(*) --> CH(3)S(*) + H(2)O reaction were identified using MP2(full) 6-311+g(2df,p) ab initio calculations. An adduct, CH(3)S(H)OH(*), I, with electronic energy 13.63 kJ mol(-1) lower than the reactants, and a transition state, II(double dagger), located 5.14 kJ mol(-1) above I, are identified as the entrance channel for an addition-elimination reaction mechanism. After adding zero-point and thermal energies, DeltaH(r,298) ( degrees )(reactants --> I) = -4.85 kJ mol(-1) and DeltaH(298) (double dagger)(I --> II(double dagger)) = +0.10 kJ mol(-1), which indicates that the potential energy surface is broad and flat near the transition state. The calculated imaginary vibrational frequency of the transition state, 62i cm(-1), is also consistent with an addition-elimination mechanism. These calculations are consistent with experimental observations of the OH(*) + CH(3)SH reaction that favored an addition-elimination mechanism rather than direct hydrogen atom abstraction. An alternative reaction, CH(3)SH + OH(*) --> CH(3)SOH + H(*), with DeltaH(r,298) ( degrees ) = +56.94 kJ mol(-1) was also studied, leading to a determination of DeltaH(f,298) ( degrees )(CH(3)SOH) = -149.8 kJ mol(-1).     

Reaction mechanism of HCN+ + C2H4: a theoretical study

The complex doublet potential energy surface for the ion-molecule reaction of HCN(+) with C(2)H(4) is investigated at the B3LYP/6-311G(d,p) and CCSD(T)/6-311++G(3df,2pd) (single-point) levels. The initial association between HCN(+) and C(2)H(4) forms three energy-rich addition intermediates, 1 (HCNCH(2)CH(2)(+)), 2 (HC-cNCH(2)CH(2)(+)), and 3 (N-cCHCH(2)CH(2)(+)), which are predicted to undergo subsequent isomerization and decomposition steps. A total of nine kinds of dissociation products, including P(1) (HCN + C(2)H(4)(+)), P(2) (HCNCHCH(2)(+) + H), P(3) (NCCH(2) + CH(3)(+)), P(4) (CN + C(2)H(5)(+)), P(5) (NCCHCH(2)(+) + H(2)), P(6) (HNCCHCH(2)(+) + H), P(7) (c-CHCCH(2)N(+) + H(2)), P(8) (c-NHCCH(2)C(+) + H(2)), and P(9) (HNCCCH(+) + H(2) + H), are obtained. Among the nine products, P(1) is the most abundant product. P(2) is the second feasible product but is much less competitive than P(1). P(3), P(4), P(5), and P(6) may have the lowest yields observed. Other products, P(7), P(8), and P(9), may become feasible at high temperature. Because the intermediates and transition states involved in the most favorable pathway all lie below the reactant, the HCN(+) + C(2)H(4) reaction is expected to be rapid, which is confirmed by experiment. The present calculation results may provide a useful guide for understanding the mechanism of HCN(+) toward other unsaturated hydrocarbons.     

Theoretical mechanistic study on the radical-molecule reaction of CH2OH with NO2

The complex singlet potential energy surface for the reaction of CH2OH with NO2, including 14 minimum isomers and 28 transition states, is explored theoretically at the B3LYP/6-311G(d,p) and Gaussian-3 (single-point) levels. The initial association between CH2OH and NO2 is found to be the carbon-to-nitrogen approach forming an adduct HOCH2NO2 (1) with no barrier, followed by C-N bond rupture along with a concerted H-shift leading to product P1 (CH2O + trans-HONO), which is the most abundant. Much less competitively, 1 can undergo the C-O bond formation along with C-N bond rupture to isomer HOCH2ONO (2), which will take subsequent cis-trans conversion and dissociation to P2 (HOCHO + HNO), P3 (CH2O + HNO2), and P4 (CH2O + cis-HONO) with comparable yields. The obtained species CH2O in primary product P1 is in good agreement with kinetic detection in experiment. Because the intermediate and transition state involved in the most favorable pathway all lie blow the reactants, the CH2OH + NO2 reaction is expected to be rapid, as is confirmed by experiment. These calculations indicate that the title reaction proceeds mostly through singlet pathways; less go through triplet pathways. In addition, a mechanistic comparison is made with the reactions CH3 + NO2 and CH3O + NO2. The present results can lead us to deeply understand the mechanism of the title reaction and may be helpful for understanding NO2-combustion chemistry.     

Ab initio study of the F + CH3NHNH2 reaction mechanism

The F + CH(3)NHNH(2) reaction mechanism is studied based on ab initio quantum chemistry methods as follows: the minimum energy paths (MEPs) are computed at the UMP2/6-311++G(d,p) level; the geometries, harmonic vibrational frequencies, and energies of all stationary points are predicted at the same level of theory; further, the energies of stationary points and the points along the MEPs are refined by UCCSD(T)/6-311++g(3df,2p). The ab initio study shows that, when the F atom approaches CH(3)NHNH(2), the heavy atoms, namely N and C atoms, are the favorable combining points. For the two N atoms, two prereaction complexes with C(s) symmetry are generated and there exists seven possible subsequent reaction routes, of which routes 1, 2, 5, and 7 are the main channels. Routes 1, 2, and 5 are associated with HF elimination, with H from the amino group or imido group, and route 7 involves the N-N bond break. Routes 3 and 6 with relation to HF elimination with H from methyl, and route 4 involved the C-N bond break, are all energetically disfavored. For the C atom, the attack of F results in the break of the C-N bond and the products are CH(3)F + NHNH(2). This route is very competitive.