CH2CO＋CH2的多通道反应的理论研究

利用密度泛函(DFT)和自然键轨道理论(NBO)及高级电子耦合簇[CCSD(T)]和电子密度拓扑(AIM)方法, 对单重态和三重态CH₂与CH₂CO反应的微观机理进行了研究. 在B3LYP/6-311＋G(d,p)水平上优化了反应通道各驻点的几何构型. 在CCSD(T)/6-311＋G(d,p)水平上计算了各物种的单点能量, 并对总能量进行了校正. 计算表明, 单重态CH₂与CH₂CO的C—H键可发生插入反应, 与C＝C、C＝O可发生加成反应, 存在三条反应通道, 产物为CO和C₂H₄, 从能量变化和反应速控步骤能垒两方面考虑, 反应II更容易发生. 对反应通道中的关键点进行了自然键轨道及电子密度拓扑分析. 三重态CH₂与CH₂CO的反应存在三条反应通道, 一条是与C-H键的插入反应, 另一条是三重态CH₂与C＝C发生加成反应, 产物为CO和三重态C₂H₄, 通道II势垒较低, 更容易发生. 最后一条涉及双自由基的反应活化能最大, 最难发生.     

二氯卡宾与甲醇及甲硫醇插入与氢抽提反应的量子化学计算

利用密度泛函和自然键轨道理论及电子密度拓扑分析方法，对单、三重态CCl2与CH3MH（M=O，S）中C—H和M—H键的插入反应及抽提氢反应进行了研究．在B3LYP／6-311G（d，p）水平上优化了势能面上构型，并以频率分析和内禀反应坐标法进行了确认．计算了各物种的CCSD（T）／6-211G（d，p）能量．结果表明，主反应通道主要发生在单重态势能面中，单重态CCl2既可以与C-H及M—H键发生插入反应，存在四条主反应通道，分别生成P1[CH3OCHCl2，反应Ⅰ（1）]，P3[Cl2HCCH2OH，反应Ⅰ（2）]和PS[CH3SCHCl2，反应Ⅱ（1）]，P7[Cl2HCCH2SH，反应Ⅱ（2）]，也可以与CH3MH发生抽提氢反应，分别生成P4[CH2O＋CH2Cl2，反应Ⅰ（3）]和P8[CH2S＋CH2Cl2，反应Ⅱ（3）]．同时，存在三重态CCl2与CH3SH插入生成^3P4[CH3SCHCl＋Cl]的反应通道．进一步对反应通道上的关键点进行了自然键轨道和电子密度拓扑分析．     

密度泛函理论研究SiC与C2H4反应的产物通道

在B3LYP/6-311++G(d,p)水平上研究了SiC与乙烯的单重态和三重态反应机理,优化得到了反应物、过渡态、中间体和产物的几何构型;通过振动分析对过渡态和中间体构型进行了确认。在CCSD(T)/cc-pVTZ水平上对计算得到的构型进行了能量校正。计算结果表明,SiC+C2H4反应在单重态和三重态条件下均可发生,其中单重态反应为主反应通道,1P5为主产物。     

CH_3CH_2O与HO_2抽氢反应与主通道速率常数的理论研究

在B3LYP//6-311++G(2df,2p)水平完成了对CH_3CH_2O与HO_2反应各驻点物种的几何构型优化,并在相同水平上对相关物种进行了频率分析和内禀反应坐标(IRC)计算.为得到较准确的单点能信息,同时采用CCSD(T)/cc-pVTZ方法对反应途径中各驻点进行单点能校正.标题反应的单、三重态势能剖面图揭示,反应在单、三重态势能面上均有3条抽氢通道.其中单重态通道分别生成1CH3CHO+H_2O_2(R1),CH_3CH_2OH+1 O_2(R2)和1CH_2CH_2O+H_2O_2(R3),三重态通道分别生成3CH3CHO+H_2O_2(R4),CH_3CH_2OH+3O_2(R5)和3CH_2CH_2O+H_2O_2(R6);势垒高度揭示三重态在动力学和热力学上比单重态更具优势.200K~1 200K区间内标题反应速率常数计算结果表明,除通道R5的速率常数几乎不随温度变化外,所有其他通道的速率常数均随温度升高而增大(即速率呈现正温度系数效应),同时发现抽氢通道R5的分支比始终大于94%,因而是绝对优势通道.     

CH_3SH和CH_3SLi分子间的锂键和氢键共存型相互作用

在CH3SLi+CH3SH势能面上求得锂键和氢键共存型复合物的两种稳定构型.频率分析表明,与单体相比复合物中S(5)—Li(6)键伸缩振动频率发生红移,而C(8)—H(10)键伸缩振动频率发生蓝移.经B3LYP/6-311++G**,MP2/6-311++G**及MP2/AUG-CC-PVDZ水平计算的含基组重叠误差(BSSE)校正的复合物Ⅰ中相互作用能分别为-58.99,-57.87和-62.89kJ·mol-1.采用自然键轨道(NBO)理论,分析了复合物中单体轨道间的电荷转移,电子密度重排及其与相关键键长变化的本质等.采用分子中的原子(AIM)理论分析了复合物中氢键和锂键的电子密度拓扑性质.在极化连续模型(PCM)下,考察了溶剂化效应.结果表明,所考察的水、二甲亚砜、乙醇和乙醚等四种溶剂均使单体间的相互作用能增大,且溶剂对复合物中的锂键结构及其振动频率具有显著的影响,而对复合物中的氢键的振动频率影响不大.     

CH_2=CClF→(hv) ·CH=CClF+H光解反应的理论研究

采用UHF,CIS和CASSCF方法,在aug-cc-pvdz基组水平上对CH2=CClF?h?v→?CH=CClF+H的光解反应通道及其后续反应作了研究.计算表明:分子吸收一个光子后,在第一电子激发态(S1)经过一个过渡态解离与Cl原子同侧的C—H键,这与用CIS方法计算垂直激发得到的π→σ*C-H跃迁及其对Frank-Condon点的计算中分子的单占轨道和键电荷密度变化所预测的结果是一致的.光解产物?CH=CClF(基态)还可再发生反应,经过渡态解离C—Cl键或是C—F键.     

H CH_2CO反应机理的G2计算

分别在UQCISD/6-311G(d,p)和G2理论计算水平上,对CH2CO和H反应可能存在的四条反应通道进行了研究,详细分析了每个通道的反应机理;通过振动分析的虚频数和内禀反应坐标(IRC)计算,确认了反应涉及的每一个过渡态.通过反应位能剖面的比较,发现经过一个中间体生成CH3 CO的一条途径是主反应通道,该通道是个放热反应,总焓变为-146.07kJ·mol-1,速控步骤的位垒为55.09kJ·mol-1.理论计算结果较好地解释了实验观察到的主要产物和副产物并存的现象.     

S_N2反应X_l~-+CH_3X_r→X_lCH_3+X_l~-(X_l=X_r=F,Cl,Br,I)的价键方法研究

应用最近发展的价键组态相互作用 ( VBCI)方法计算了 SN2反应 X-l +CH3 Xr→ Xl CH3 +X-r ( Xl=Xr=F,Cl,Br,I)的反应能垒和价键相关参数 .计算结果表明 ,VBCI能垒与采用分子轨道理论的 CCSD( T)方法计算的能垒相一致 .讨论了 SN2反应的反应参数 .     

HO_2+NO气相反应机理及主通道速率常数的理论研究

采用CCSD(T)/aug-cc-p VTZ//MP2/aug-cc-p VTZ方法构建了HO_2+NO反应的单、三重态反应势能剖面.研究结果表明,标题反应在单、三重态均存在3条反应通道,分别经抽氧、加成和抽氢反应生成产物NO_2+OH、HNO_3和HNO+O_2.其中抽氧、加成反应的优势通道发生在单重态势能面上(通道R1和R2),而抽氢的优势通道则发生在三重态势能面上(通道R6).利用经典过渡态理论(TST)并结合Wigner矫正模型在200K~1500K温度范围内计算了优势通道R1,R2和R6的速率常数kTST/W.结果显示,在计算温度范围内抽氧、加成反应的通道R1和R2始终存在竞争,且抽氢反应通道R6在高温区500K~1500K竞争作用也逐渐显现,但抽氧通道R1分支比始终大于81%,具有绝对优势.     

CH2=CCIFhv→·CH=CCIF+H光解反应的理论研究

采用UHF,CIS和CASSCF方法,在aug-cc-pvdz基组水平上对CH2=CCIFhv→·CH=CClF+H的光解反应通道及其后续反应作了研究.计算表明:分子吸收一个光子后,在第一电子激发态(S1)经过一个过渡态解离与Cl原子同侧的C-H键,这与用CIS方法计算垂直激发得到的π→σ*C-H跃迁及其对Frank-Condon点的计算中分子的单占轨道和键电荷密度变化所预测的结果是一致的.光解产物·CH=CCIF(基态)还可再发生反应,经过渡态解离C-Cl键或是C-F键.     

OH-与CH2ClF反应的阴离子产物通道

理论研究了羟基负离子(OH^-)与氟氯代甲烷(CH₂ClF)反应的阴离子产物通道. 分别在B3LYP/6-31+G(d,p)和B3LYP/6-311++G(2d,p)水平上得到反应势能面上各关键物种的优化构型, 进而计算得到谐振频率和零点能. 基于CCSD(T)/6-311+G(3df,3dp)水平的相对能量, 描述了由质子转移和双分子亲核取代(S_N2)过程生成各阴离子产物的途径. 各阴离子产物途径势垒的计算结果表明质子转移过程是实验中的主要产物通道, 与以往实验测量的结论相符. 此外, 计算还显示双分子亲核取代过程得到了非典型的阴离子产物, 其中动力学效应可能会导致F^-的生成.     

CH3CH2,CH3CHCl,CH3CCl2与NO2反应微观动力学理论研究

运用量子化学密度泛函理论UB3LYP/6-311+G*和高级电子相关校正的偶合簇(CCSD(T)/6-311+G*)方法,对CH3CH2,CH3CHCl和CH3CCl2自由基与NO2反应的机理和动力学进行了理论研究,得到了体系的势能面信息和可能的反应机理.根据计算得到的各反应热力学参数及反应能垒,采用传统过渡态理论计算了各反应在温度T=298 K和T=700 K时的速率常数.研究结果表明,该类反应均通过1个中间体和1个过渡态生成产物,产物分别为CH3CHO+HNO,CH3CHO+ClNO和CH3CClO+ClNO.     

H-Atom abstraction from CH(3)NHNH(2) by NO(2): CCSD(T)/6-311++G(3df,2p)//MPWB1K/6-31+G(d,p) and CCSD(T)/6-311+G(2df,p)//CCSD/6-31+G(d,p) calculations

Stationary points of paths for H atom abstraction from CH(3)NHNH(2) (monomethylhydrazine) by NO(2) were characterized via CCSD(T)/6-311++G(3df,2p)//MPWB1K/6-31+G(d,p) and CCSD(T)/6-311+G(2df,p)//CCSD/6-31+G(d,p) calculations. Five transition states connecting CH(3)NHNH(2)-NO(2) complexes to a manifold that includes CH(3)NHNH-HONO, CH(3)NNH(2)-HONO, CH(3)NNH(2)-HNO(2), and CH(3)NHNH-HNO(2) complexes were identified. Transition states that connect CH(3)NHNH-HONO, CH(3)NNH(2)-HONO, CH(3)NNH(2)-HNO(2), and CH(3)NHNH-HNO(2) complexes to each other via H atom exchange and/or hindered internal rotation were also identified. The high point in the minimum energy path from the CH(3)NHNH(2) + NO(2) reactant asymptote to the manifold of HONO-containing product states is a transition state 8.6 kcal/mol above the reactant asymptote. From a kinetics standpoint, this value is considerably higher than the 5.9 kcal/mol value that was estimated for it based on theoretical results for H atom abstraction from NH(3) by NO(2).     

亚甲基自由基(3CH2)与So反应机理的理论研究

白洪涛,黄旭日,于广涛,李吉来,于健康,孙家钟. 亚甲基自由基(3CH2)与SO反应机理的理论研究[J]. 化学学报, 2006, 64(2): 139-144.     

Insight into global reaction mechanism of [C2, H4, O] system from ab initio calculations by the scaled hypersphere search method

A detailed computational study is performed on the singlet potential energy surface (PES) for possible isomerization and dissociation reactions of CH(3)CHO at the DFT (B3LYP/6-311++G(d,p)) and CCSD(T)/cc-pVTZ//B3LYP/6-311++G(d,p) levels. The pathways around the equilibrium structures can be discovered by the scaled hypersphere search (SHS) method, which enables us to make a global analysis of the PES for a given chemical composition. Fourteen isomers inclusive of 11 single-molecules and three "non-stabilized" oxygen-based ylides, 5 energetically favored complexes, and 79 interconversion transition states have been found on the singlet PES. Four lowest lying isomers with thermodynamic stability are also kinetically stable with respect to metastable intermediates. It was revealed that vinyl alcohols, which could be generated by the tautomerization of acetaldehyde, could undergo dissociation to form acetylene and water. In addition, recombination channels between some fragments, such as H(2)CO + (1)CH(2) and (1)CHOH + (1)CH(2), are energetically accessible via collision complex or oxygen-based ylides. Most of available unimolecular decompositions are found to be responsible for favorable hydrogen abstraction processes.     

Substituent effects on singlet-triplet gaps and mechanisms of 1,2-rearrangements of 1,3-oxazol-2-ylidenes to 1,3-oxazoles

Electronic structures, partial atomic charges, singlet-triplet gaps (Delta E ST), substituent effects, and mechanisms of 1,2-rearrangements of 1,3-oxazol-2-ylidene ( 5) and 4,5-dimethyl- ( 6), 4,5-difluoro- ( 7), 4,5-dichloro- ( 8), 4,5-dibromo- ( 9), and 3-methyl-1,3-oxazol-2-ylidene ( 10) to the corresponding 1,3-oxazoles have been studied using complete-basis-set methods (CBS-QB3, CBS-Q, CBS-4M), second-order M?ller-Plesset perturbation method (MP2), hybrid density functionals (B3LYP, B3PW91), coupled-cluster theory with single and double excitations (CCSD) and CCSD plus perturbative triple excitations [CCSD(T)], and the quadratic configuration interaction method including single and double excitations (QCISD) and QCISD plus perturbative triple excitations [QCISD(T)]. The 6-311G(d,p), 6-31+G(d,p), 6-311+G(d,p), and correlation-consistent polarized valence double-xi (cc-pVDZ) basis sets were employed. The carbenes have singlet ground states, and the CBS-QB3 and CBS-Q methods predict Delta E ST values for 5- 8 and 10 of 79.9, 79.8, 74.7, 77.0, and 82.0 kcal/mol, respectively. CCSD(T), QCISD(T), B3LYP, and B3PW91 predict smaller Delta E ST values than CBS-QB3 and CBS-Q, with the hybrid density functionals predicting the smallest values. The concerted unimolecular exothermic out-of-plane 1,2-rearrangements of singlet 1,3-oxazol-2-ylidenes to their respective 1,3-oxazoles proceed via cyclic three-center transition states. The CBS-predicted barriers to the 1,2-rearrangements of singlet carbenes 5- 9 to their respective 1,3-oxazoles are 41.4, 40.4, 37.8, 40.4, and 40.5 kcal/mol, respectively. During the 1,2-rearrangements of singlet 1,3-oxazol-2-ylidenes 5- 9, there is a decrease in electron density at oxygen, N3 (the migration origin), and C5 and an increase in electron density at C2 (the migration terminus), C4, and the partially positive migrating hydrogen.     

A systematic computational study on the reactions of HO2 with RO2: The HO2 + CH3O2(CD3O2) and HO2 + CH2FO2 reactions

A systematic theoretical study of the reactions of HO2 with RO2 has been carried out. The major concern of the present work is to gain insight into the reaction mechanism and then to explain experimental observations and to predict new product channels for this class of reactions of importance in the atmosphere. In this paper, the reaction mechanisms for two reactions, namely, HO2 + CH3O2 and HO2 + CH2FO2, are reported. Both singlet and triplet potential energy surfaces are investigated. The complexity of the present system makes it impossible to use a single ab initio method to map out all the reaction paths. Various ab initio methods including MP2, CISD, QCISD(T), CCSD(T), CASSCF, and density function theory (B3LYP) have been employed with the basis sets ranging from 6-31G(d) to an extrapolated complete basis set (CBS) limit. It has been established that the CCSD(T)/cc-pVDZ//B3LYP/6-311G(d,p) scheme represents the most feasible method for our systematic study. For the HO2 + CH3O2 reaction, the production of CH3OOH is determined to be the dominant channel. For the HO2 + CH2FO2 reaction, both CH2FOOH and CHFO are major products, whereas the formation of CHFO is dominant in the overall reaction. The computational findings give a fair explanation for the experimental observation of the products.     

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.     

Theoretical Studies on Structures and Stabilities of C4H2+ Isomers

The structures, energies, stabilities and spectroscopies of doublet C₄H₂⁺ cations were explored at the DFT/B3LYP/6-311G(d,p), CCSD(T)/6-311+G(2df,2pd)(single-point), and G3B3 levels. Ten minimum isomers including the chainlike, three-member-ring, and four-member-ring structures are interconverted by means of 15 interconversion transition states. The potential energy surface was investigated. At the CCSD(T)/6-311+G(2df,2pd) and G3B3 levels, the global minimum isomer was found to be a linear HCCCCH. The structures of the stable isomer and its relevant transition state are further optimized at the QCISD/6-311G(d,p) level. The bonding nature and structure of isomer HCCCCH were analyzed.     

A theoretical study of the H-abstraction reactions from HOI by moist air radiolytic products (H, OH, and O (3P)) and iodine atoms (2P(3/2))

The rate constants of the reactions of HOI molecules with H, OH, O ((3)P), and I ((2)P(3/2)) atoms have been estimated over the temperature range 300-2500 K using four different levels of theory. Geometry optimizations and vibrational frequency calculations are performed using MP2 methods combined with two basis sets (cc-pVTZ and 6-311G(d,p)). Single-point energy calculations are performed with the highly correlated ab initio coupled cluster method in the space of single, double, and triple (pertubatively) electron excitations CCSD(T) using the cc-pVTZ, cc-pVQZ, 6-311+G(3df,2p), and 6-311++G(3df,3pd) basis sets. Reaction enthalpies at 0 K were calculated at the CCSD(T)/cc-pVnZ//MP2/cc-pVTZ (n = T and Q), CCSD(T)/6-311+G(3df,2p)//MP2/6-311G(d,p), and CCSD(T)/6-311++G(3df,3pd)//MP2/6-311G(d,p) levels of theory and compared to the experimental values taken from the literature. Canonical transition-state theory with an Eckart tunneling correction is used to predict the rate constants as a function of temperature. The computational procedure has been used to predict rate constants for H-abstraction elementary reactions because there are actually no literature data to which the calculated rate constants can be directly compared. The final objective is to implement kinetics of gaseous reactions in the ASTEC (accident source term evaluation code) program to improve speciation of fission products, which can be transported along the reactor coolant system (RCS) of a pressurized water reactor (PWR) in the case of a severe accident.