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.     

Br+CH3S(O)CH3反应机理的直接动力学研究

采用直接动力学的方法，对多通道反应体系Br＋CH3S（O）CH3进行了理论研究．在BH＆H-LYP／6-311G（2d，2p）水平下获得了优化几何构型、频率及最小能量路径（MEP），能量信息的进一步确认在MC-QCISD（单点）水平下完成．利用正则变分过渡态理论，结合小曲率隧道效应校正（CVT／SCT）方法计算了该反应的两个可行的反应通道在200K～2000K温度范围内的速率常数．在整个反应区间内，生成HBr的反应通道与生成CHa的反应通道存在着竞争，前者是主反应通道，后者是次反应通道．变分效应和小曲率隧道效应对反应速率常数的计算影响都很小．理论计算得到的两个反应通道的反应速率常数与实验值符合得很好．     

Dual-level direct dynamics studies on the reaction Cl + CHBr(2)Cl

Theoretical investigations are carried out on the multichannel reaction CHBr(2)Cl + Cl by means of direct dynamics methods. The minimum energy path (MEP) is obtained at the BH&H-LYP/6-311G(d,p) level, and energetic information is further refined at the CCSD(T)/6-311+G(2df,2p) (single-point) level. The rate constants for three reaction channels, H-abstraction, Br-abstraction, and Cl-abstraction, are calculated by using the improved canonical variational transition state theory (ICVT) incorporating with the small-curvature tunneling (SCT) correction. The theoretical overall rate constants are in good agreement with the available experimental data and are found to be k=2.58 x 10(-15) T(1.18) exp(-861.17/T) cm(3)molecule(-1)s(-1) over the temperature range 200--2400 K. For the title reaction, H-abstraction reaction channel is the major channel at the lower temperatures, while as the temperature increases, the contribution of Br-abstraction reaction channel should be taken into account. At 2180 K, the rate constants of these two pathways are equal. Cl-abstraction reaction channel is minor channel over the whole temperature region.     

A study of the hydrogen abstraction reactions of C2H radical with CH3CN, C2H5CN, and C3H7CN by dual-level generalized transition state theory

The hydrogen abstraction reactions C2H + CH3CN --> products (R1), C2H + CH3CH2CN --> products (R2), and C2H + CH3CH2CH2CN --> products (R3) have been investigated by dual-level generalized transition state theory. Optimized geometries and frequencies of all the stationary points and extra points along the minimum-energy path (MEP) are performed at the BH&H-LYP and MP2 methods with the 6-311G(d, p) basis set, and the energy profiles are further refined at the MC-QCISD level of theory. The rate constants are evaluated using canonical variational transition state theory (CVT) with a small-curvature tunneling correction (SCT) over a wide temperature range 104-2000 K. The calculated CVT/SCT rate constants are in good agreement with the available experimental values. Our calculations show that for reaction R2, the alpha-hydrogen abstraction channel and beta-hydrogen abstraction channel are competitive over the whole temperature range. For reaction R3, the gamma-hydrogen abstraction channel is preferred at lower temperatures, while the contribution of beta-hydrogen abstraction will become more significant with a temperature increase. The branching ratio to the alpha-hydrogen abstraction channel is found negligible over the whole temperature range.     

Theoretical study and rate constant calculation for reaction of CF(3)CH(2)OH with OH

The reaction mechanism of CF(3)CH(2)OH with OH is investigated theoretically and the rate constants are calculated by direct dynamics method. The potential energy surface (PES) information, which is necessary for dynamics calculation, is obtained at the B3LYP/6-311G (d, p) level. The single-point energy calculations are performed at the MC-QCISD level using the B3LYP geometries. Complexes, with the energies being less than corresponding reactants and products, are found at the entrance and exit channels for methylene-H-abstraction channel, while for the hydroxyl-H-abstraction channel only entrance complex is located. By means of isodesmic reactions, the enthalpies of the formation for the species CF(3)CH(2)OH, CF(3)CHOH, and CF(3)CH(2)O are estimated at the MC-QCISD//B3LYP/6-311G (d, p) level of theory. The rate constants for two kinds of H-abstraction channels are evaluated by canonical variational transition state theory with the small-curvature tunneling correction (CVT/SCT) over a wide range of temperature 200-2000 K. The calculated results are in good agreement with the experimental values in the temperature region 250-430 K. The present results indicate that the two channels are competitive. Below 289 K, hydroxyl-H-abstraction channel has more contribution to the total rate constants than methylene-H-abstraction channel, while above 289 K, methylene-H-abstraction channel becomes more important and then becomes the major reaction channel.     

甲硫醇和氢原子反应的从头算直接动力学

利用双水平直接动力学方法对反应CH3SH+H的微观机理和动力学性质进行了理论研究.对于此反应的三个反应通道,即—SH和—CH3基团上的两个氢提取通道及一个取代通道,在MP2/6-311+G(d,p)水平上优化得到了各稳定点的结构及振动频率,并在G3(MP2)水平上进行了单点能量计算以获得更精确的能量信息;在此基础上运用结合小曲率隧道效应校正的变分过渡态理论(CVT/SCT)计算了各反应通道在220-1000 K温度区间的速率常数.计算结果表明提取—SH基团上H的反应通道R1在整个反应温度区间都是主要通道,而随着温度的升高,低温下的次要反应通道——取代通道R3变得越来越重要,并且在高温下将成为一个竞争的反应通道;提取—CH3基团上H的反应通道(R2)由于具有较高的反应能垒,因而,其对总反应速率常数的贡献可以忽略.计算得到的总反应速率常数与已有的实验值符合得很好,进而我们预测了该反应在220-1000 K温度范围内速率常数的表达式为:k=5.00×10-18T2.39exp(-119.81/T),为将来的实验研究提供参考.     

HNCS与CH2(X2Π)反应微观动力学的理论研究

用量子化学密度泛函理论的UB3LYP/6-311+G**方法和高级电子相关的UQCISD(T)/6-311+G**方法研究了异硫氰酸(HNCS)与乙炔基自由基(C2H(X2Π))反应的微观机理. 采用双水平直接动力学方法IVTST-M, 获取反应的势能面信息, 应用正则变分过渡态理论并考虑小曲率隧道效应, 计算了在250～2500 K温度范围内反应的速率常数. 研究结果表明, HNCS与C2H(X2Π)反应为多通道、多步骤的复杂反应, 共存在三个可能的反应通道, 主反应通道为通过分子间H原子迁移, 生成主要产物NCS+C2H2. 反应速率常数随温度升高而增大, 表现为正温度效应. 速率常数计算中变分效果很小. 在低温区隧道效应对反应速率的贡献较大, 反应为放热反应.     

C2H3与CH3F氢抽提反应机理与动力学性质

采用双水平直接动力学方法对C₂H₃与CH₃F氢抽提反应进行了研究. 在QCISD(T)/6-311++G(d, p)//B3LYP/6-311G(d, p)水平上, 计算的三个反应通道R1、R2和R3的能垒(ΔE^≠)分别为43.2、43.9和44.1 kJ·mol^-1, 反应热为-38.2 kJ·mol^-1. 此外, 利用传统过渡态理论(TST)、正则变分过渡态理论(CVT)和包含小曲率隧道效应(SCT)的CVT, 分别计算了200-3000 K温度范围内反应的速率常数k^TST、k^CVT和k^CVT/SCT. 结果表明: (1) 三个氢抽提反应通道的速率常数随温度的增加而增大, 其中变分效应的影响可以忽略, 隧道效应则在低温段影响显著; (2) R1反应是主反应通道, 但随着温度的升高, R2反应的竞争力增大, 而R3反应对总速率常数的影响很小.     

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     

DFT and ab initio dual-level direct dynamics studies on the reactions of fluorine atom with HOCl and HOBr

The multichannel reactions (1) HOCl + F --> products and (2) HOBr + F --> products have been investigated using the dual-level direct dynamics method. The minimum energy paths (MEPs) are calculated at both the MPW1K/6-311G(d,p) and QCISD/6-311G(d,p) levels, then the single-point energies are further corrected at the QCISD(T)/6-311++G(3df,3pd) level of theory. There are hydrogen-bonded complexes with the energies less than those of the reactants or products located at the entrance or exit channel of both hydrogen abstraction reactions; while for the halogen abstraction channels only one complex exists at the reactant side in the bromine abstraction channel. The rate constants are evaluated by the improved canonical variational transition-state theory (ICVT). The agreement of the rate constants with available experimental values for two reactions at room temperature is good. Theoretical results indicate that for the reaction HOCl + F, hydrogen abstraction channel leading to the formation of HF + ClO will predominate the reaction over the whole temperature range, and the reaction of HOBr + F may proceed mainly through the bromine abstraction channel at the lower temperature while the contribution of hydrogen abstraction channel will become significant as the temperature increases. Because of lack of the kinetic data of these reactions, the present theoretical results are expected to be useful and reasonable to estimate the dynamical properties of these reactions over a wide temperature range where no experimental value is available.     

CH_3自由基与C_2H_5CN高温反应机理的理论研究

用G3(MP2)//B3PW91/6-311G(d,p)双级别方法研究了CH_3自由基与C_2H5_CN的反应机理和动力学性质.计算表明反应存在抽氢、加成-消除和取代3种机理7条反应通道.用CVT方法计算了所有反应通道在1 000K~3 000K温度范围内的速率常数,结果表明计算值与实验值符合得很好.     

CH3S自由基H迁移异构化及脱H2反应的直接动力学研究

采用密度泛函方法(MPW1PW91)在6-311G(d,p)基组水平上研究了CH₃S自由基H迁移反应CH₃S→CH₂SH (R1), 脱H₂反应CH₃S→HCS＋H₂ (R2)以及脱H₂产物HCS异构化反应HCS→CSH (R3)的微观动力学机理. 在QCISD(t)/6- 311＋＋G(d,p)//MPW1PW91/6-311G(d,p)＋ZPE水平上进行了单点能校正. 利用经典过渡态理论(TST)与变分过渡态理论(CVT)分别计算了各反应在200～2000 K温度区间内的速率常数k^TST和k^CVT, 同时获得了经小曲率隧道效应模型(SCT)校正后的速率常数k^CVT/SCT. 结果表明, 反应 R1, R2 和R3的势垒△E^≠分别为160.69, 266.61和241.63 kJ/mol, R1为反应的主通道. 低温下CH₃S比CH₂SH稳定, 高温时CH₂SH比CH₃S更稳定. 另外, 速率常数计算结果显示, 量子力学隧道效应在低温段对速率常数的计算有显著影响, 而变分效应在计算温度段内对速率常数的影响可以忽略.     

CH3SH与NO2的反应机理及动力学

在G3B3, CCSD(T)/6-311＋＋G(d,p)//B3LYP/6-311＋＋G(d,p)水平上详细研究了CH3SH与基态NO2的微观反应机理. 在B3LYP/6-311＋＋G(d,p)水平得到了反应势能面上所有反应物、过渡态和产物的优化构型, 通过振动频率分析和内禀反应坐标(IRC)跟踪验证了过渡态与反应物和产物的连接关系. 在CCSD(T)/6-311＋＋G(d,p)和G3B3水平计算了各物种的能量, 得到了反应势能面. 利用经典过渡态理论(TST)与变分过渡态理论(CVT)并结合小曲率隧道效应模型(SCT), 分别计算了在200～3000 K温度范围内的速率常数kTST, kCVT和kCVT/SCT. 研究结果表明, 该反应体系共存在5个反应通道, 其中N进攻巯基上H原子生成CH3S＋HNO2的通道活化势垒较低, 为主要反应通道. 动力学数据也表明, 该通道在200～3000 K计算温度范围内占绝对优势, 拟合得到的速率常数表达式为k1CVT/SCT＝1.93×10－16T0.21exp(－558.2/T) cm3&;#8226;molecule－1&;#8226;s－1.     

CF3CH2CH3+OH的反应机理及动力学性质的理论研究

采用密度泛函理论BB1K/6-31+G(d,p)计算了反应CF₃CH₂CH₃+OH各反应通道上驻点的稳定结构和振动频率, 并分别在BMC-CCSD, MC-QCISD和G3(MP2)水平上进行了单点能校正. 运用变分过渡态理论, 在BMC-CCSD//BB1K, MC-QCISD//BB1K, G3(MP2)//BB1K以及BB1K水平上计算了各反应通道的速率常数, 讨论了-CH₂和-CH₃基团上H提取通道对总反应的贡献, 并与已有实验和理论结果进行了对比. 计算结果表明, BMC-CCSD水平上的速率常数与实验测量值符合得很好, 进而给出了该水平上反应在200~1000 K温度范围内速率常数k(cm³?molecule^-1?s^-1)的三参数表达式: k=1.90×10^-21T^3.21exp(-292.62/T).     

HNCS与Cl原子的反应机理及电子密度拓扑分析

在QCISD(T)/6-311++G(d,p)和B3LYP/6-311++G(d,p)级别上研究了HNCS与Cl原子的反应机理. 并应用经典过渡态理论和正则变分过渡态理论结合小曲率隧道效应, 计算了200-2500 K温度范围内各反应通道的速率常数. 结果表明, HNCS与Cl原子反应存在3个反应通道. 当温度低于294 K时, 生成HCl+NCS的夺氢反应(a)是优势通道, 温度高于294 K时, 生成HNC(Cl)S的加成反应(c)为主反应通道, Cl进攻N的反应通道(b)因能垒较高而难以进行.     

H原子与CH3NH2抽氢反应的动力学计算

用QC ISD(T)/6-311 G(3DF,3PD)/MP2/6-311G(D,P)方法研究了H原子与CH3NH2的抽氢反应过程。该反应包含两个反应通道:H分别从CH3基团(R1)和NH2(R2)基团上抽氢。R1势垒比R2势垒低3.42kJ/mol,表明R1是主反应通道。在从头算的基础上,用变分过渡态理论(CVT)加小曲率隧道效应(SCT)研究了各反应温度范围为200~4000K内的速率常数,所得结果与实验值符合的很好。动力计算表明,在所研究的温度范围内,变分效应对速率常数的计算影响不大,而在低温范围内,隧道效应起了很重要的作用。     

Reaction Mechanism and Kinetics for HCCO Radical with NO

The mechanism and dynamical properties for the reaction of HCCO radicals with NO were investigated theoretically. The minimum energy paths(MEP) of the reaction were calculated by using the density functional theory(DFT) at the B3LYP/6-311 G^** level, and the energies along the MEP were further refined at the QCISD(T)/6-311 G^** level. It is found that the reaction mechanism of the title reaction involves three channels, producing HCNO CO, HONC CO and HCN CO2 products, respectively. Channel 1 is the most favorable path. The rate constant for channel 1 were calculated over a temperature range of 800-2500 K by using the canonical variational transition-state theory(CVT). The rate constant for the main path is negatively dependent on temperature, which is a characteristic of radical reactions with negative activation energy, and the variational effect for the rate constant calculation is small in the whole temperature range.     

Theoretical Studies on the Abstraction Reaction of Atomic O（^3p） with Si2H6

The hydrogen abstraction reaction of O(^3P) with Si2H6 has been studied theoretially. Two transition states of ^3A″ and ^3A′ symmetries have been located for this abstraction reaction. Geometries have been optimized at the UMP2 leve with 6-311G (d) basis set. G3MP2 has been used for the final single-point energy calculation. The rate constants have been calculated over a wide temperature range of 200-3000K using canonical variational transition-state sheory (CVT) with small curvature tunneling effect(SCT). The calculated CVT/SCT rate constants match well with the experimental value.     

Mechanism and Kinetics for the Reaction of NCS and OH Radicals

The mechanism and dynamical properties for the reaction of NCS and OH radicals have been investigated theoretically. The minimum energy paths (MEP) of the reaction were calculated using the density functional theory(DFT) at the B3LYP/6-311 G^** level, and the energies along the MEP were further refined at the QCISD(T)/6-311 G^** level. As a result, the reaction mechanism of the title reaction involves three channels, producing HCS NO and HNC SO products, respectively. Path Ⅰ and path Ⅱ are competitive, with some advantages for path Ⅰ in kinetics. As for path Ⅲ, it looks difficult to react for its high energy barrier. Moreover, the rate constant have been calculated over the temperature range of 8190-2500K using canonical variational transition-state theory (CVT). It was found that the rate constants for both path Ⅰ and path Ⅱ are negatively dependent on temperature, which is similarwith the experimental results for reactions of NCS with NO and NO2, and the variational effect for the rate constant calculation olavs an important role in whole temperature range.     

Theoretical studies of the reaction of hydroxyl radical with methyl acetate

The mechanisms and the kinetics of the OH (OD) radicals with methyl acetate CH3C(O)OCH3 are investigated theoretically. The dual-level direct dynamics method is employed in the calculation of the rate constants. The optimized geometries and frequencies and the gradients of the stationary points are calculated at the MP2/6-311G(d,p) level. The energetic information of potential energy surfaces is further refined by the multicoefficient correlation method based on QCISD (MC-QCISD) using the MP2/6-311G(d,p) geometries. Four channels are found for the title reaction. The calculated results reveal that there exists an attractive well (reactant complex) in each entrance H-abstraction channel, that is, the H-abstraction reaction makes a stepwise mechanism. The rate constants are calculated by the canonical variational transition-state theory (CVT) with the interpolated single-point energies (ISPE) approach in the temperature range of 200-1200 K. The small-curvature tunneling effect (SCT) approximation is used to evaluate the transmission coefficient. The calculated rate constants are in good agreement with the experimental ones in the measured temperature range. It is shown that the "out-of-plane hydrogen abstraction" from the methoxy end is the dominant channel at the lower temperatures, and the other two H-abstraction channels should be taken into account with the temperatures increasing. The kinetic isotope effects (KIEs) for the three H-abstraction channels and the total reaction are "inverse", and these theoretically calculated KIEs as a function of temperature are expected to be useful for the future laboratory investigation.