HCNO自由基与羟基在气相中反应机理的理论研究

在CCSD(T)/B3LYP/6-311G(d,p)//B3LYP/6-311G(d,p)+ZPE水平上对反应HCNO+OH进行了计算,建立了反应势能面,对反应中涉及到的6个中间体和12个过渡态都做了详尽的分析.详细阐明了理论上可能得到的7种产物:P1为H2O+CNO,P2为HCO+HNO,P3为HO2+HCN,P4为HONH+CO,P5为H2CO+NO,P6为H2NO+CO和P7为H2O+OCN,以及形成这些产物的各种反应通道.其中最主要通道为由反应物形成反式初始复合物,再连续经过2次1,3-氢迁移最终形成产物HONH+CO,该通道是一条热力学可行的反应通道.并且从反应物、中间体和产物的相对能量来看,此反应是典型的消除型反应.另外,直接的氢提取反应也是比较重要的反应通道.     

HS和HOO抽氢反应机理的研究

在B3LYP/6-311+G(3df,2p)水平上对HS和HOO反应中的所有物种进行了几何构型优化和频率计算,采用QCISD(T)/6-311+G(3df,2p)方法获得了各物种的单点能,构建了HS和HOO反应在单、三重态势能剖面.结果表明,HS与HOO反应体系中存在2种不同的抽氢通道,在单、三重态势能面上生成的产物分别为[1P1(H2O2+1S),1P2(H2S+1 O2)]和[3P1(H2O2+3S),3P2(H2S+3O2)].标题反应主要发生在三重态势能面上,优势通道[R→3 TS2→3P2(H2S+3O2)]的活化能为9.99kJ·mol-1.此结果对认识大气硫迁移转变规律具有实际意义.     

气相中烯丙基负离子与N2O反应机理

采用二阶微扰理论的MP2/6-31G(d,p)方法对气相中烯丙基负离子与N2O的反应机理进行了理论计算研究, 并在相同基组下进一步用CCSD(T)方法进行了单点能的校正. 计算结果表明, 该反应存在三条反应通道, 产物分别为cis-CH2CHCNN-+H2O, trans-CH2CHCNN-+H2O和CH2CCH-+N2+H2O, 其中生成cis-CH2CHCNN-和trans-CH2CHCNN-的两条通道为相互竞争的主反应通道, 计算结果与实验相吻合. 同时利用传统的过渡态理论, 计算了各反应通道在298 K时, 速控步骤的反应速率常数k(T).     

·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生成后就很难再回到初始反应物.     

Criegee中间体RCHOO(R=H,CH_3)与NCO的反应机理

采用CCSD(T)/aug-cc-p VTZ//B3LYP/6-311+G(2df,2p)方法对Criegee中间体RCHOO(R=H,CH_3)与NCO反应的机理进行了研究,利用经典过渡态理论(TST)并结合Eckart校正模型计算了标题反应在298~500 K范围内优势通道的速率常数.结果表明,上述反应包含亲核加成、氧化和抽氢3类机理,其中每类又包括NCO中N和O分别进攻的两种形式.亲核加成反应中O端进攻为优势通道,氧化和抽氢反应则是N端进攻为优势通道;甲基取代使CH_3CHOO反应活性高于CH2OO;anti-CH_3CHOO的加成及氧化反应活性高于syn-CH_3CHOO,而抽氢反应则是syn-CH_3CHOO的活性高于anti-CH_3CHOO.anti-构象对总速率常数的贡献大于syn-构象,且总速率常数具有显著的负温度效应.     

单分子水对H_2O_2+Cl气相反应影响的理论研究

在aug-cc-pVTZ基组下采用CCSD(T)和B3LYP方法,研究了H2O2+Cl反应,并考虑在大气中单个水分子对该反应的影响.结果表明,H2O2+Cl反应只存在一条生成产物为HO2+HCl的通道,其表观活化能为10.21kJ·mol-1.加入一分子水后,H2O2+Cl反应的产物并没有发生改变,但是所得势能面却比裸反应复杂得多,经历了RW1、RW2和RW3三条通道.水分子在通道RW1和RW2中对产物生成能垒的降低起显著的负催化作用,而在通道RW3中则起明显的正催化作用.利用经典过渡态理论(TST)并结合Wigner矫正模型计算了216.7-298.2 K温度范围内标题反应的速率常数.结果显示,298.2 K时通道R1的速率常数为1.60×10-13cm3·molecule-1·s-1,与所测实验值非常接近.此外,尽管通道RW3的速率常数kRW3比对应裸反应的速率常数kR1大了46.6-131倍,但该通道的有效速率常数k'RW3却比kR1小了10-14个数量级,表明在实际大气环境中水分子对H2O2+Cl反应几乎没有影响.     

F+CH3OH碰撞反应机理和反应势能面

在MP2(full)/6-311++g(d,p)水平上详细研究了氟原子与甲醇抽氢反应的多通道反应机理,得到了各条通道中涉及的驻点的构型和振动频率及其能量,给出了两张完整的反应势能面.结果表明,氟原子从C原子上抽氢时有一条明显的最低能量通道,而从氧原子上抽氢时要涉及多条分支通道和多个驻点构型,给出了各分支通道的势能面示意图,结果表明以形成五元环状过渡态通道为优势通道.计算得到经途径1生成CH2OH时反应放热170.62kJ/mol,经分支途径6生成CH3O自由基时反应放热119.41 kJ/mol,此结果与实验值一致.     

1,2-环丙烷乙酰化糖水解机理的理论研究

采用密度泛函B3LYP方法研究了1,2-环丙烷乙酰化糖在氘代氯仿溶液中的水解反应的详细机理。计算结果表明,当H2O分子从不同方向进攻1,2-环丙烷乙酰化糖分子时,会形成不同的反应途径,当H2O从糖分子的六元环平面下方进攻时,反应为一步反应,环丙烷的开环步骤为反应决速步,生成α构型产物。当H2O从糖分子的六元环平面上方进攻时,反应为两步反应,形成产物的氢迁移步骤为决速步,生成β构型产物。H2O从糖分子平面下方进攻的反应途径在热力学及动力学上都更有优势,1,2-环丙烷乙酰化糖的水解反应更有利于生成α构型产物,计算结果与实验结果一致。     

NO2和HSO反应的机理及动力学研究

本文采用CCSD(T)/aug-cc-pVTZ//B3LYP/aug-cc-pVTZ方法构建了NO2 + HSO反应的单、三重态势能面,并对主通道速率常数进行了计算。研究结果表明,该反应在单[R1(HN(O)O + 1SO)、R2(cis-HONO + 1SO)和R3 (trans-HONO + 1SO)]、三重态[R6(HN(O)O + 3SO)、R7(cis-HONO + 3SO)和R8(trans-HONO + 3SO)]均存在3条抽氢反应通道,在单[R4(NO + HS(O)O)和R5(H + SO2 + NO)]、三重态[R9(HS(O)O + NO)和R10(H + SO2 + NO)]均存在两条抽氧通道,其中单重态抽氢通道R2 (cis-HONO + 1SO)是NO2 + HSO反应主通道。利用传统过渡态理论(TST)并结合Wigner校正,计算了上述10条通道在200 ~ 1000 K温度范围内的速率常数。计算结果表明,NO2 + HSO反应主通道在298 K时的速率常数(7.78×10-13cm3?molecule-1?s-1)与实验值(9.6×10-12 cm3?molecule-1?s-1)相吻合。此外,水分子影响主通道R2(cis-HONO + 1SO)经历了NO2 + H2O…HSO和 NO2 + H2O…HSO(HSO + NO2…H2O)两条反应通道,且两条通道的能垒分别比R2升高了49.97和20.43 kcal?mol-1,说明在实际大气环境中水分子对NO2 + HSO反应几乎没有影响。     

CH_4与NO_2反应的微观机理及动力学性质的理论研究

采用量子化学计算方法,研究了CH_4+NO_2反应直接氢抽提反应通道的机理和速率常数.该反应有3条反应通道分别生成CH_3+HNO_2,CH_3+trans-HONO和CH_3+cis-HONO.计算结果表明采用变分过渡态理论加小曲率隧道效应校正计算得到反应速率常数和已有的实验值很吻合.在整个研究温度区间,O原子提取H原子生成CH_3+cis-HONO是反应的主要通道.     

Mechanistic Investigation on the Reaction of O- with CH3CN Using Density Functional Theory

The potential energy profile of the reaction between the atomic oxygen radical anion and acetonitrile has been mapped at the G3MP2B3 level of theory. Geometries of the reactants, products, intermediate complexes, and transition states involved in this reaction have been optimized at the (U)B3LYP/6-31+G(d,p) level, and then their accurate relative energies have been improved using the G3MP2B3 method. The potential energy profile is confirmed via intrinsic reaction coordinate calculations of transition states. Four possible production channels are examined respectively, as H⁺ transfer, H-atom transfer, H²⁺ transfer, and bi-molecular nucleophilic substitution (S_N2) reaction pathways. Based on present calculations, the H²⁺ transfer reaction is major among these four channels, which agrees with previous experimental conclusions     

The O + NH3 reaction: A review

Experimental data for the reaction of O atoms with NH₃ have been reviewed with particular attention to the possible effects of secondary reactions on the deduced rate coefficient. A reaction mechanism was assembled and computer calculations carried out to simulate several sets of experiments. The sensitivity of the calculations to uncertainties in the various rate coefficients was assessed. Transition-state theory calculations were carried out on the rate coefficient k₁ for the O + NH₃ → OH + NH₂ reaction. These studies suggest that the reaction stoichiometry is dependent on temperature, initial reagent ratios, and extent of reaction; that available data are not sufficient to determine whether the initial step is H-atom abstraction (producing OH and NH₂ radicals) or O-atom addition (producing an NH₃O* complex): and that the low temperature values of k₁ (T ? 400 K) are not consistent with values deduced at higher temperatures if the reaction proceeds by H-atom abstraction. Taking all the evidence into account, it is recommended that the expression 1.1 × 10³ T^2.1 exp(?2620/T) L mol^?1 s^?1 be used for k₁ until more definitive experimental measurements are made at low temperatures.     

Ab initio study of the unimolecular decomposition of 2‐butenenitrile: Molecular elimination channels

Ab initio molecular orbital calculations have been performed for the unimolecular decomposition of 2‐butenenitrile (CH₃CH?CHCN), especially for HCN and H₂ molecular elimination channels. Structures and energies of the reactants, products, and relevant species in the individual reaction pathways were determined by MP2 gradient optimization and MP4 CCSD(T) single‐point energy calculations. Direct 1,1 and 1,2 molecular eliminations and H or CN migration followed by elimination channels were identified. Dissociation rates for the individual reaction pathways were calculated from vibrational frequencies at the ab initio transition state geometries by employing Rice–Ramsperger–Kassel–Marcus theory, from which channel branching ratios were determined. It was concluded that the most important reaction channel should be the direct 1,1 three‐center molecular elimination of HCN. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

A computational study of the radical–radical reaction of O(<Superscript>3</Superscript>P) + C<Subscript>2</Subscript>H<Subscript>5</Subscript> with comparisons to gas-phase kinetics and crossed-beam experiments

We present density functional theory (DFT) and complete basis set (CBS) calculations of the prototypical radical–radical reaction of ground–state atomic oxygen [O(³P)] with ethyl (C₂H₅) radicals. The respective reaction mechanisms and dynamics were investigated on the doublet potential energy surfaces using the DFT method and CBS model. In the title reaction, the barrierless addition of O(³P) to C₂H₅ led to the formation of energy-rich intermediates that underwent subsequent isomerization and decomposition to yield various products. The products predicted to be found were: H₂CO + CH₃, CH₃CHO + H, c–CH₂OCH₂ + H, ^1,3CH₃COH + H, ^1,3HCOH + CH₃, CH₂CHOH + H, C₂H₃ + H₂O, and CH₂CH₂ + OH. In particular, unlike previous kinetic results, proposed to proceed only through the direct H-atom abstraction process, two distinctive pathways to the formation of CH₂CH₂ + OH were predicted to be in competition: direct, barrierless H-atom abstraction mechanism versus addition process. The competition was consistent with the recent crossed-beam investigations, and their microscopic dynamic characteristics are discussed at the molecular level.     

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     

Computational study on the kinetics and mechanisms for the reactions of HCO with HONO and HNOH

The kinetics and mechanisms of the HCO reactions with HONO and HNOH have been studied at the G2M level of theory based on the geometric parameters optimized at BH&HLYP/6‐311G(d,p). The rate constants in the temperature range 200–3000 K at different pressures have been predicted by microcanonical RRKM and/or variational transition state theory calculations with Eckart tunneling corrections. For the HCO + HONO reaction, hydrogen abstraction from trans‐HONO and cis‐HONO by HCO produces H₂CO + NO₂, with the latter being dominant. Two other channels involving cis‐HONO by the association/decomposition mechanism via the HC(O)N(O)OH intermediate, which could fragment to give H₂O + CO + NO at high temperatures, were also found to be important. For the HCO + HNOH reaction, three reaction channels were identified: one association reaction giving a stable intermediate, HC(O)N(H)OH (LM2), and two hydrogen abstraction channels producing H₂CO and H₂NOH. The dominant products were predicted to be the formation of LM2 at low temperatures and H₂NOH + CO at middle and high temperatures. © 2004 Wiley Periodicals, Inc. Int J Chem Kinet 36: 178–187 2004     

Quantum chemical investigation of the reaction of O(3 P 2) with certain hydrocarbon radicals

The reaction of ground-state atomic oxygen [O(³ P ₂)] with methyl, ethyl, n-propyl and isopropyl radicals has been studied using the density functional method and the complete basis set model. The energies of the reactants, products, reaction intermediates and various transition states as well as the reaction enthalpies have been computed. The possible product channels and the reaction pathways are identified in each case. In the case of methyl radical the minimum energy reaction pathway leads to the products CO + H₂ + H. In the case of ethyl radical the most facile pathway leads to the products, methanal + CH₃ radical. For propyl radical (n- and iso-), the minimum energy reaction pathways would lead to the channel containing ethanal + methyl radical.     

SIFDT study of the SO+2/H2 H-atom abstraction reaction

The exothermic H-atom abstraction reaction of SO⁺₂ with H₂ has been studied in a selected ion flow drift tube (SIFDT) over a range of center-of-mass energies from thermal (300 K) to about 0.12 eV. The measured rate coefficient at 300 K is 4.2 × 10⁻¹² cm³ s⁻¹ which is very much less than the Langevin capture rate. The increase in rate coefficient with ion kinetic energy gives a linear Arrhenius-type plot with a slope that indicates a barrier of ∼5 kJ mol⁻¹ exists on the potential surface. The H₂SO⁺₂ potential surface is also explored in an ab initio investigation using the G2 procedure. An (SO⁺₂.H₂)¹ transition state between reactants and products is identified, corresponding to the barrier found from experiments.     

Theoretical investigation on the reaction of N2O and CO catalyzed by Fe(C6H6)

The reaction of N₂O with CO, catalyzed by Fe⁺(C₆H₆) and producing N₂ and CO₂, has been investigated at the UB3LYP/6-311+G(d) level. The computation results revealed that the reaction of Fe⁺(C₆H₆), N₂O and CO, is an O-atom abstraction mechanism. For the reaction channels, the geometries and the vibrational frequencies of all species have been calculated and the frequency modes analysis also have been given to elucidate the reaction mechanism. On the basis for geometry optimizations, the thermodynamic data of these reactions channels have been calculated using the statistical theory at 295.15 K and pressure of 0.35 Torr. Using Eyring transition state theory with Wigner correction, the activation thermodynamic data, rate constant and frequency factors for the these reaction channels also have been given. The results showed that CO and N₂O do not react without catalyst and Fe⁺(C₆H₆) can excellently mediate the reaction of N₂O and CO.     

Theoretical study on the hydrogen abstraction reactions from hydrazine derivatives by H atom

The hydrogen abstraction reactions from hydrazine and its methyl derivatives by the H atom have been investigated theoretically by using CBS-QB3//DSD-BLYP-D3(BJ)/Def2-TZVP quantum chemical calculations and transition state theory calculations coupled with various tunneling correction methods. Both the products and transition state energies of the hydrogen abstraction from the amino group were stabilized by the methyl group substitution. The substitution effect on the ^αN site was two times larger than that on the ^βN site. On the other hand, the substitution effect was negligible on the hydrogen abstraction from the methyl group. The overall rate coefficients of N₂H₄ + H reaction calculated by canonical variational transition state theory with the small-curvature tunneling correction agreed well with previously reported values, but those of CH₃NHNH₂/(CH₃)₂NNH₂ + H were slightly lower than a previous experimental value. The product-specific rate coefficients have been proposed for the kinetics modeling of these fuels’ combustion.