HCCO与CH(2Π)双自由基反应微观动力学的理论研究

用量子化学密度泛函理论的UB3LYP/6-311+G**方法和高级电子相关的UQCISD(T)/6-311+G**方法研究了HCCO与CH(2Π)自由基反应的微观机理. 采用双水平直接动力学方法IVTST-M和正则变分过渡态理论研究了在1 000～2 500 K温度范围内反应的速率常数. 结果表明, HCCO与CH(2Π)双自由基反应过程中存在3个反应通道, 生成产物为C2H2+CO. 通道2为主要反应路径, 通道1也占一定的比例. 在所研究的温度范围内, 速率常数计算的变分效果均较小, 反应为放热反应.     

氢抽提反应CHBr2+HBr→CH2Br2+Br的直接动力学研究

采用直接动力学方法,对CHBr2+HBr→CH2Br2+Br反应通道进行了理论研究,在B3LYP/6-311+G(d,p)水平下获得了优化几何构型、频率以及最小能量路径,更精确的单点能在B3LYP/6-311++G(3df,2pd)水平下完成.利用正则变分过渡态理论,结合小曲率隧道效应校正方法计算了反应通道在220 K～2 000 K温度范围内的速率常数.在整个反应区间,隧道效应对反应的影响比较大;变分效应在低温时有一定的影响,在高温区间的影响很小可以忽略.计算得到的速率常数和已有实验值很好地吻合.     

N_2O与氧负离子微观反应机理的理论研究

采用G3B3//B3LYP理论水平对反应O-+N2O的双重电子态势能面反应机理进行了详细的理论研究.该反应涉及的各个稳定点的构型、振动频率是在B3LYP/6-311++G(d,p)理论水平下计算的.计算结果表明,得到的反应焓变与已有实验值相吻合,该反应主反应通道是O-+N2O→NO+NO-,而生成O2-+N2的反应通道是次反应通道.     

HCCO与CH（2∏）双自由基反应微观动力学的理论研究

用量子化学密度泛函理论的UB3LY/6-311 G**方法和高级电子相关的UQCISD(T)／6-311 G**方法研究了HCCO与CH(^2∏)自由基反应的微观机理．采用双水平直接动力学方法IVTST-M和正则变分过渡态理论研究了在l000～2500K温度范围内反应的速率常数．结果表明，HCCO与CH(^2∏)双自由基反应过程中存在3个反应通道。生成产物为C2H2 CO．通道2为主要反应路径，通道1也占一定的比例．在所研究的温度范围内，速率常数计算的变分效果均较小，反应为放热反应．     

SiH2自由基与HNCO反应机理的理论研究

采用密度泛函理论B3LYP方法研究了SiH₂自由基与HNCO的反应机理, 并在B3LYP/6-311＋＋G**水平上对反应物、中间体、过渡态进行了全几何参数优化, 通过频率分析和内禀反应坐标(IRC)确定了中间体和过渡态. 为了得到更精确的能量值, 又用QCISD(T)/6-311＋＋G**方法计算了在B3LYP/6-311＋＋G**水平优化后的各个驻点的相对能量. 计算结果表明SiH₂自由基与HNCO的反应有五条反应通道, 其中顺式反应通道SiH₂＋HNCO→IM3→ TS4→IM5→TS5→IM6→SiH₂NH＋CO反应能垒最低, 为主反应通道.     

SiH2自由基与HNCO反应机理的理论研究

采用密度泛函理论B3LYP方法研究了SiH₂自由基与HNCO的反应机理, 并在B3LYP/6-311＋＋G**水平上对反应物、中间体、过渡态进行了全几何参数优化, 通过频率分析和内禀反应坐标(IRC)确定了中间体和过渡态. 为了得到更精确的能量值, 又用QCISD(T)/6-311＋＋G**方法计算了在B3LYP/6-311＋＋G**水平优化后的各个驻点的相对能量. 计算结果表明SiH₂自由基与HNCO的反应有五条反应通道, 其中顺式反应通道SiH₂＋HNCO→IM3→ TS4→IM5→TS5→IM6→SiH₂NH＋CO反应能垒最低, 为主反应通道.     

O(3P)与卤代甲基反应的理论研究--Ⅰ.O(3P)+CH2Cl的反应机理

用量子化学从头算法对氧原子与CH2 Cl自由基的反应进行了研究 ,采用G2MP2方法计算了势能面上各驻点的构型参数、振动频率和能量 .给出了O( 3 P)与CH2 Cl反应的明确机理 .反应首先形成富能中间体OCH2 Cl,而后经各种复杂的解离或异构化途径生成产物 .计算的各个通道的反应热与实验结果相符 ,预测H CHClO和Cl CH2 O是反应的主要通道 .根据从头算的结果 ,用过渡态理论计算了反应的总速率常数 .反应速率常数与压力无关 ,在低温下有弱的负温度效应 .计算值与实验值符合很好 .     

Cl_2 2HI=2HCl I_2反应机理的理论研究

分别在MP2/3-21G!!、CCSD(T)/3-21G!!//MP2/3-21G!!和B3LYP/3-21G!!3种水平上,计算研究了气相反应Cl2 2HI=2HCl I2的机理,求得一系列四中心和三中心的过渡态.通过比较六种反应通道的活化能大小,得到了相同的结论:双分子基元反应Cl2 HI"HCl ICl和ICl HI"I2 HCl的最小活化能小于Cl2、HI和ICl的解离能,从理论上证明了反应Cl2 2HI=2HCl I2将优先以分子与分子作用形式分两步完成.用内禀反应坐标(IRC)验证了MP2/3-21G!!方法计算得到的过渡态.     

硅烯与异硫氰酸反应机理的量子化学研究

采用密度泛函理论B3LYP方法在B3LYP/6-311＋＋G**水平上对反应物、中间体、过渡态进行了全几何参数优化, 通过频率分析和IRC方法确认了中间体和过渡态. 又用QCISD(T)/6-311＋＋G**//B3LYP/6-311＋＋G**方法计算了各个驻点的单点能, 计算结果表明单重态的硅烯与异硫氰酸的反应有抽提硫、插入、抽提亚氨基的路径. 而经由三元环中间体的抽提硫反应SiH2＋HNCS→IM1→TS2→IM3→TS3→IM4→SiH2S＋HNC(P1), 反应能垒最低, 为主反应通道, 硫代硅甲醛和异氰氢酸为主产物. 硅烯直接抽提硫、插入N—H键和经由三元环中间体的亚氨基抽提反应为竞争反应通道, 在室温下可以发生, 应为次反应通道.     

CH_3SH+CN微观反应机理的理论研究

采用BMC-CCSD//B3LYP/6-311G(d,p)方法对CH3SH+CN反应机理进行了详细的理论研究.反应中涉及的各稳定点的构型、振动频率和零点能在B3LYP/6-311G(d,p)水平下计算得到,计算结果表明,该反应存在两种反应机理,5条可能的反应通道.SN2机理由于能垒太高,与直接氢抽提机理相比可以忽略.该反应的最可行通道为CN中的C原子进攻SH中的H原子经由一个前期和一个后期分子络合物生成产物CH3S和HCN.计算得到的反应焓变与已有实验值非常吻合.     

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.     

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.     

Theoretical study and rate constant calculation for the reactions of SH (SD) with Cl2, Br2, and BrCl

The mechanisms of the SH (SD) radicals with Cl2 (R1), Br2 (R2), and BrCl (R3) are investigated theoretically, and the rate constants are calculated using a dual-level direct dynamics method. The optimized geometries and frequencies of the stationary points are calculated at the MP2/6-311G(d,p) and MPW1K/6-311G(d,p) levels. Higher-level energies are obtained at the approximate QCISD(T)/6-311++G(3df, 2pd) level using the MP2 geometries as well as by the multicoefficient correlation method based on QCISD (MC-QCISD) using the MPW1K geometries. Complexes with energies less than those of the reactants or products are located at the entrance or the exit channels of these reactions, which indicate that the reactions may proceed via an indirect mechanism. The enthalpies of formation for the species XSH/XSD (X = Cl and Br) are evaluated using hydrogenation working reactions method. By canonical variational transition-state theory (CVT), the rate constants of SH and SD radicals with Cl2, Br2, and BrCl are calculated over a wide temperature range of 200-2000 K at the a-QCISD(T)/6-311++G(3df, 2pd)//MP2/6-311G(d, p) level. Good agreement between the calculated and experimental rate constants is obtained in the measured temperature range. Our calculations show that for SH (SD) + BrCl reaction bromine abstraction (R3a or R3a') leading to the formation of BrSH (BrSD) + Cl in a barrierless process dominants the reaction with the branching ratios for channels 3a and 3a' of 99% at 298 K, which is quite different from the experimental result of k3a'/k3' = 54 +/- 10%. Negative activation energies are found at the higher level for the SH + Br2 and SH + BrCl (Br-abstraction) reactions; as a result, the rate constants show a slightly negative temperature dependence, which is consistent with the determination in the literature. The kinetic isotope effects for the three reactions are "inverse". The values of kH/kD are 0.88, 0.91, and 0.69 at room temperature, respectively, and they increase as the temperature increases.     

A dual-level ab initio and hybrid density functional theory dynamics study on the unimolecular decomposition reaction C2H5O-->CH2O + CH3

We present a direct ab initio and hybrid density functional theory dynamics study of the thermal rate constants of the unimolecular decomposition reaction of C2H5O-->CH2O + CH3 at a high-pressure limit. MPW1K/6-31+G(d,p), MP2/6-31+G(d,p), and MP2(full)/6-31G(d) methods were employed to optimize the geometries of all stationary points and to calculate the minimum energy path (MEP). The energies of all the stationary points were refined at a series of multicoefficient and multilevel methods. Among all methods, the QCISD(T)/aug-cc-pVTZ energies are in good agreement with the available experimental data. The rate constants were evaluated based on the energetics from the QCISD(T)/aug-cc-pVTZ//MPW1K/6-31+G(d,p) level of theory using both microcanonical variational transition state theory (microVT) and RRKM theory with the Eckart tunneling correction in the temperature range of 300-2500 K. The calculated rate constants at the QCISD(T)/aug-cc-pVTZ/MPW1K/6-31+G(d,p) level of theory are in good consistent with experimental data. The fitted three-parameter Arrhenius expression from the microVT/Eckart rate constants in the temperature range 200-2500 K is k = 2.52 x 10(12)T(0.41)e(-8894.0/T) s(-1). The falloff curves of pressure-dependent rate constants are performed using master-equation method within the temperature range of 391-471 K. The calculated results are in good agreement with the available experimental data.     

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

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

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.     

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内的速率常数,所得结果与实验值符合的很好。动力计算表明,在所研究的温度范围内,变分效应对速率常数的计算影响不大,而在低温范围内,隧道效应起了很重要的作用。     

Theoretical study for the reaction of CH3OCl with Cl atom

A direct dynamics method is employed to study the kinetics of the multiple channel reaction CH(3)OCl + Cl. The potential energy surface (PES) information is explored from ab initio calculations. Two reaction channels, Cl- and H-abstractions, have been identified. The optimized geometries and frequencies of the stationary points and the minimum-energy paths (MEPs) are calculated at the MP2 level of theory using the 6-311G(d, p) and cc-pVTZ basis sets, respectively. The single-point energies along the MEPs are further refined at the G3(MP2)//MP2/6-311G(d, p), G3//MP2/6-311G(d, p), as well as by the multicoefficient correlation method based on QCISD (MC-QCISD) using the MP2/cc-pVTZ geometries. The enthalpies of formation for the species CH(3)OCl and CH(2)OCl are calculated via isodesmic reactions. The rate constants of the two reaction channels are evaluated by using the variational transition-state theory over a wide range of temperature, 220-2200 K. The calculated rate constants exhibit the slightly negative temperature dependence and show good agreement with the available experimental data at room temperature at the G3(MP2)//MP2/6-311G(d, p) level. The present calculations indicate that the two channels are competitive at low temperatures while H-abstraction plays a more important role with the increase of temperature. The calculated k(1a)/k(1) ratio of 0.5 at 298 K is in general agreement with the experimental one, 0.8 +/- 0.2. The high rate constant for CH(3)OCl + Cl shows that removal by reaction with Cl atom is a potentially important loss process for CH(3)OCl in the polar stratosphere.     

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.     

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.