Quasi-Classical Trajectory Investigation of H+SO$ _{2} \rightarrow $OH+SO Reaction on Full-Dimensional Accurate Potential Energy Surface

The reaction H+SO\begin{document}$ _2 $\end{document}\begin{document}$ \rightarrow $\end{document}OH+SO is important in the combustion and atmospheric chemistry, as well as the interstellar medium. It also represents a typical complex-forming reaction with deep complexes, serving as an ideal candidate for testing various kinetics theories and providing interesting reaction dynamical phenomena. In this work, we reported a quasi-classical trajectory study of this reaction on our previously developed accurate full-dimensional potential energy surface. The experimental thermal rate coefficients over the temperature range 1400 K\begin{document}$ \leq $\end{document}\begin{document}$ T $\end{document}\begin{document}$ \leq $\end{document} 2200 K were well reproduced. For the reactant SO\begin{document}$ _2 $\end{document} being sampled at the ground ro-vibrational state, the calculated integral cross sections increased slightly along the collision energy ranging from 31.0 kcal/mol to 40.0 kcal/mol, and then became essentially flat at the collision energy within 40.0\begin{document}$ - $\end{document}55.0 kcal/mol. The product angular distributions are almost symmetric with nearly identical backward-forward double peak structure. The products OH and SO vibrational state distributions were also analyzed.     

O+HCl→OH+Cl反应的立体动力学的准经典轨线研究

运用准经典轨线方法, 基于Peterson从头计算势能面对O+HCl→OH+Cl反应的立体动力学性质进行了研究. 讨论了在31.77和51.04 kJ/mol两种碰撞能情况下极化依赖的微分反应截面(2π/σ)(dσ₀₀/dω_t), (2π/σ)(dσ₂₀/dω_t), (2π/σ)(dσ₂₂₊/dω_t)和(2π/σ)(dσ_21-/dω_t)以及描述k-j′两矢量相关和k-k′-j′三矢量相关的分布函数P(θ_r)和P(φ_r). 计算得到的P(θ_r)分布表明, 产物分子的转动角动量j′具有强烈的取向分布, 并且产物转动角动量的取向效应对散射能的变化比较敏感. 而P(φ_r)的分布表明, 产物分子虽然具有沿着y轴的取向效应, 但是没有明显的定向效应.     

准经典轨线方法对Ca+CD3Ⅰ→CaⅠ+CD3同位素效应的动力学研究

在扩展Lond-Eyring-Polanyi-Sato(LEPS)势能面上,采用准经典轨线方法对反应Ca+CD3I→CaI+CD3进行了动力学计算,并讨论了该反应的同位素效应.在同位素效应作用下,产物CaI的振动态分布向低振动态转移,反应体系的散射截面在低碰撞能和高碰撞能处有较小的变化.同时,受到反应物的质量因子变化的影响,产物转动取向值减少,产物转动取向增强.仅有产物的角分布受同位素效应的影响不明显.     

基态钙原子与二氯甲烷反应的动力学研究

利用准经典轨线方法在扩展的Lond-Eyring-Polanyi-Sato势能面上对Ca+CH2Cl2→CaCl+CH2Cl反应进行了动力学研究. 在低碰撞能Ecol=9.632 kJ/mol下, 计算得到产物CaCl的振动分布峰位于ν=0. 散射截面随着碰撞能的增加而减少. 产物的角分布趋向后向散射. 随着碰撞能的增加, 产物的转动取向趋向强烈, 〈P2(J’,K)〉略大于-0.5. 采用准经典轨线方法计算得到的结果与实验及其它理论计算结果相符合.     

反应物振动激发对反应He+H+2→HeH++H影响的准经典轨线研究

基于PPA势能面(Palmieri, P.; Puzzarini, C.; Aquilanti, V. Mol. Phys., 2000, 98: 1835), 运用准经典轨线(quasi-classical trajectory)方法, 对反应He+H+2 (v=0-2, j=0)→HeH++H的立体动力学性质进行了理论研究. 结果表明, 反应物振动激发对反应的k-j'两矢量相关和k-k'-j'三矢量相关分布都产生了较大影响. 除此之外, 极化微分反应截面对振动量子数v也非常敏感.     

CH3NO异构体及异构化反应机理的量子化学研究

在CCSD(T)/6-311＋＋G(d,p)//B3LYP/6-311＋＋G(d,p)水平上对CH₃NO的异构体及各异构化反应机理进行了详细的理论研究. 包括CH₃NO共找到了23种异构体, 这23种异构体间通过24个过渡态相互转化. 用振动频率分析证实了过渡态的真实性, 内禀反应坐标跟踪(IRC)计算验证了过渡态与各异构体的连接关系. 计算结果表明, 在CH₃NO的各种异构体中, 分子骨架中存在N—C—O结构是最稳定的, 而C—O—N, C—N—O结构能量最高, 不稳定. 异构体间的异构化反应主要通过化学键的旋转, H迁移, C, N, O骨架的重组及包含上述两种机理的混合型机理等4种方式完成.     

Theoretical study of the stereodynamics of the reactions O(D) + H2, D2 and HD

Calculations of the dynamics of the reactions O(¹D) + H₂ → OH + H, O(¹D) + HD → OH + D, O(¹D) + HD → OD + H and O(¹D) + D₂ → OD + D have been performed using the quasi-classical trajectory (QCT) method with symplectic integration. The theoretical calculations were carried out on the ground state 1A′ potential energy surfaces (PES) by Dobbyn and Knowles. The distributions of the dihedral angle P(_r), the angle between k and j′, P(θ_r), and the product vibrational state are presented. The results show that the intermediate geometrical structures and lifetimes of the reactive collisions play a vital role in these reactions.     

硝基甲烷异构化反应势能面的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的异构化反应较为不利.     

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.     

CH3O·2+NO气相反应的密度泛函理论研究

用密度泛函方法分别研究了单态和三态 CH3 O·2 NO CH3 O· NO2 气相反应 .结果表明 ,反应中 NO进攻 CH3 O·2 经过了一个顺反异构化的过程 ,摘取 CH3 O·2 的端基氧 .整个反应是吸热反应 ,理论计算吸热值为 5 0 .93k J/ mol,单态为多通道多步骤反应 ,决定速度步骤的能垒为 1 90 .6 1 k J/ mol.而三态为单通道反应 ,其决定速度步骤的能垒为 1 6 3.31 k J/ mol.三态反应为最佳反应通道 .该反应的研究将为保护臭氧层及大气环境提供重要的理论依据 .     

Quasiclassical Trajectory Study of the Reaction of CD4 with O(1D)

Extensive quasiclassical trajectory calculations for the O(¹D)+CD₄ multichannel reaction were carried out on a new global potential energy surface fit by permutationally invariant polynomials. The product branching ratios, translational energy distributions, and angular distributions of OD+CD₃, D+CD₂OD/CD₃O, and D₂+DCOD/D₂CO product channels were calculated and compared with the available experimental results. Good agreement between theory and experiment has been achieved, indicating small isotope effects for the title reaction. The O(¹D)+CD₄ reaction mainly proceeds through the CD₃OD intermediate via the trapped abstraction mechanism, with initial abstraction of the D atom rather than the direct insertion, followed by decomposition of CD₃OD into various products.     

Cl+HD反应产物的极化与态分布

利用准经典轨线理论,在BW2和G3两个势能面上,研究了Cl+HD反应的动力学.计算结果表明,产物的转动取向对势能面及反应体系的质量因子非常敏感.在BW2势能面上,计算的两个产物的转动取向强于在G3势能面上计算的结果,而无论是在BW2势能面上还是在G3势能面上,DCl产物的取向都强于HCl产物的取向.计算结果还表明,在不同的势能面上反应物的转动激发对反应的影响有着显著的不同.在BW2势能面上,反应物的初始转动激发有利于Cl+HD反应的进行;而在G3势能面上,反应物的初始转动激发消弱了反应的反应性.     

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.     

Product Vibrational State Distributions of F+CH3OH Reaction on Full-Dimensional Accurate Potential Energy Surface

The hydrogen abstraction reaction of methanol with fluorine atoms can produce HF and CH\begin{document}$ _3 $\end{document}O or CH\begin{document}$ _2 $\end{document}OH radicals, which are important in the environment, combustion, radiation, and interstellar chemistry. In this work, the dynamics of this typical reaction is investigated by the quasi-classical trajectory method based on a recently developed globally accurate full-dimensional potential energy surface. Particularly, the vibrational state distributions of the polyatomic products CH\begin{document}$ _3 $\end{document}O and CH\begin{document}$ _2 $\end{document}OH are determined by using the normal mode analysis method. It is found that CH\begin{document}$ _3 $\end{document}O and CH\begin{document}$ _2 $\end{document}OH are dominantly populated in the ground state when the reactants are at the ground ro-vibrational state. The OH stretching mode, torsional mode, H\begin{document}$ _2 $\end{document}CO out-of-plane bending mode and their combination bands in the CH\begin{document}$ _2 $\end{document}OH product can be effectively excited once the OH stretching mode of the reactant CH\begin{document}$ _3 $\end{document}OH is excited to the first vibrationally excited state. Most of the available energy flows into the HF vibrational energy and the translational energy in both channels, while the radical products, CH\begin{document}$ _3 $\end{document}O or CH\begin{document}$ _2 $\end{document}OH, receive a small amount of energy, consistent with experiment, which is an indication of its spectator nature.     

Quasi-classical Trajectory Study of the Intramolecular Isotope Effect in the Reaction O（3p）＋H2/HD

Theoretical studies of the dynamics of the reactions O（3p）＋H2/HD（ν=0, j=0）→OH＋H have been performed with quasi-classical trajectory method （QCT） on an ab initio potential surface for the lowest triplet electronic state of H2O（aA＂）. The QCT-calculated integral cross sections are in good agreement with the earlier time-dependent quantum mechanics results. The state-resolved rotational distributions reveal that the product OH rotational distributions for O＋HD have a preference for populating highly internally excited states compared with the O＋H2 reaction. Distributions of differential cross sections show that directions of scattering are strongly dependent on the choice of quantum state. The polarization dependent generalized differential cross-sections and the distributions were calculated and a pronounced isotopic effect is revealed. The calculated results indicate that the product polarization is very sensitive to the mass factor.     

CH_3SO+HO_2气相反应机理与主通道速率常数

在MPW1PW91/6-311G(d,p)水平上优化了标题反应各驻点物种的几何构型,并在相同水平上通过频率计算和内禀反应坐标(IRC)分析对过渡态结构及连接性进行了验证.采用QCISD(T)/6-311G(d,p)方法对所有驻点及反应路径的部分选择点进行单点能校正,分别构建了CH3SO+HO2反应体系的单、三重态反应势能剖面.研究结果表明,CH3SO+HO2反应体系存在6条反应通道7条路径,优势通道(1)R→3IM→P1(CH3SOH+3O2)发生在三重态势能面上,此通道包含两条路径,其表观活化能分别为12.01和-30.04kJ?mol-1,主路径(2)R→3IM→3TS2→P1(CH3SOH+3O2)是一个无势垒氢迁移过程.利用经典过渡态理论(TST)与变分过渡态理论(CVT)并结合小曲率隧道效应模型(SCT),分别计算了主路径(2)在200～2500K温度范围内的速率常数kTST,kCVT和kCVT/SCT,在此温度区间内的表观反应速率常数三参数表达式为kCVT/SCT=4.08×10-24T3.13exp(8012.2/T)cm3imolecule-1is-1,具有负温度系数效应.速率常数计算结果显示,变分效应在计算温度段内影响较小,而量子力学隧道效应在低温段有显著影响.     

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～3000K温度范围内的速率常数kTST,kCVT和kCVT/SCT.研究结果表明,该反应体系共存在5个反应通道,其中N进攻巯基上H原子生成CH3S+HNO2的通道活化势垒较低,为主要反应通道.动力学数据也表明,该通道在200～3000K计算温度范围内占绝对优势,拟合得到的速率常数表达式为k1CVT/SCT=1.93×10-16T0.21exp(-558.2/T)cm3·molecule-1·s-1.     

A Theoretical Study on the Potential Energy Surface of the C3 + NO Reaction

The reaction of linear form carbon cluster C₃ molecules in their ¹Σ ground state with NO (X²Π) radicals is explored theoretically to investigate the formation of hitherto undetected NCCCO molecules in the interstellar medium via a neutral-neutral reaction. The doublet potential energy surface is worked out by the ab initio MO calculations at the CCSD(T)/cc-pVTZ//B3LYP/6-311G(d,p)+ZPE level of theory. It is shown that the main pathway of the C₃(¹Σ)+NO(X²Π) reaction involves the N-atom of NO attacking the side C-atom of the ¹C₃ molecule first to form the adduct CCCNO, followed by the N-shift to give I6 CCNCO, and then to the main products P₁ (CCN+CO). Alternatively, I6 can be converted via the N-shift again to I9 (CN)CCO, and then it leads to the minor products P₂ (CNC+CO) and P₃ NCCCO. Since the three pathways have zero threshold energy and proceed via strongly bound energized complexes, they should possess the character of negative temperature dependence, and might be quite rapid even at very low temperature. The reaction represents facile neutral-neutral pathways to produce hitherto undetected CCN, CNC and NCCCO molecules in interstellar environments.     

Reaction Mechanism Investigation Using Vibrational Mode Anal—ysis for the Multichannel Reacition of CH3O＋CO

On the basis of the computed results got by the Gaussian 94 package at B3LYP/6-311 G** level,the reaction mechanism of CH3O radical with CO has been investiagted thoroughly via the vibrational model analysis ,And the relationships among the reactants,eight transition states,four intermediates and various products involved this multichannel reation are eluci-dated,The vibrational mode anaysis shows that the reaction mechanism is relialbe.