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.     

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.     

CH3CH2S自由基H迁移异构化及裂解反应的理论研究

采用密度泛函方法(B3LYP)在6-311+G(d,p)基组水平上研究了CH3CH2S自由基H迁移异构化以及裂解反应的微观动力学机理. 在QCISD(T)/6-311++G(d,p)//B3LYP/6-311+G(d,p)+ZPE水平上进行了单点能校正. 利用经典过渡态理论(TST)与变分过渡态理论(CVT)分别计算了在200~2000 K温度区间内的速率常数kTST和kCVT, 同时获得了经小曲率隧道效应模型(SCT)校正后的速率常数kCVT/SCT. 研究结果表明, CH3CH2S自由基1,2-H迁移、1,3-H迁移、C—C键断裂和β-C—H键断裂反应的势垒ΔE≠分别为149.74, 144.34, 168.79和198.29 kJ/mol. 当温度低于800 K时, 主要发生1,2-H迁移反应, 高于1800 K时, 主要表现为C—C键断裂反应, 在1300—1800 K范围内, 1,3-H迁移反应是优势通道, 在计算的整个温度段内, β-C—H键断裂反应可以忽略.     

H3PO→H2POH异构化反应的直接动力学研究

在QCISD(T)/6-311C++G(2df,2pd)//QCISD/6-311C++G(d,p)+ZPE水平上,对H3PO的异构化反应H3PO→(1)H2POH(trans)→(2)H2POH(cis)进行了计算研究.结果表明,H原子由P原子向O原子迁移反应(1)的能垒为250.0kJ/mol,是反应速率控制步骤,而O_H键绕P_O键旋转的构型转化反应(2)的能垒只为12.3kJ/mol.利用经典过渡态理论(TST)与变分过渡态理论(CVT)分别计算了反应(1)在200～2000K温度区间内的速率常数kTST和kCVT,获得了经小曲率隧道效应(SCT)及Eckart模型校正后的速率常数kTST/Eckart和kCVT/SCT.对只涉及H原子迁移的反应(1),量子力学隧道效应的影响在低温段非常明显,而变分效应对反应速率常数的影响很小.     

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,具有负温度系数效应.速率常数计算结果显示,变分效应在计算温度段内影响较小,而量子力学隧道效应在低温段有显著影响.     

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

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

CH4-nFn(n=1～3)与CH3氢抽提反应微观动力学的理论研究

采用量子化学的QCISD(T)/6-311 G(d,p)//BHandHLYP/6-311G(d,p)方法研究了氟代甲烷CH4-nFn(n=1~3)与CH3自由基氢抽提反应的微观动力学性质.并利用Polyrate程序分别计算了3个反应在200~3000K范围内的速率常数.计算结果表明,R1a,R2a和R3三个反应路径的反应能量分别为-12.7,-9.5和11.8kJ/mol,相应的能垒依次为67.0,62.2和67.5kJ/mol.在437K时,kCVT/SCT分别为6.72×10-19,8.01×10-18和8.82×10-20cm3/(molecule.s).计算结果还表明,在低温段反应的量子隧道效应显著,在计算温度范围内变分效应对反应速率常数的影响可以忽略.     

HNCS与C_2H(X~2Π)反应微观动力学的理论研究

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

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. 反应速率常数随温度升高而增大, 表现为正温度效应. 速率常数计算中变分效果很小. 在低温区隧道效应对反应速率的贡献较大, 反应为放热反应.     

HOSO+NO反应的机理及动力学

在B3LYP/6-311++G(d,p)水平上研究了HOSO+NO的反应机理. 优化得到了反应势能面上各驻点的几何构型, 通过内禀反应坐标(IRC)确认了反应物、中间体、过渡态和产物的相关性. 在CCSD(T)/6-311++G(d,p)水平上对计算得到的构型进行了能量校正. 应用经典过渡态理论(TST)与变分过渡态理论(CVT), 并结合小曲率隧道(SCT)效应模型校正的方法计算了标题反应在200-3000 K温度范围内的速率常数k^TST、k^CVT和k^CVT/SCT. 计算结果表明: HOSO+NO反应在单重态和三重态条件下均可发生, 其中单重态反应为主反应通道, HNO+SO₂为主产物. 并利用电子密度拓扑分析方法研究主反应通道反应过程中的化学键变化.     

Direct dynamic study on the hydrogen abstraction reaction CH3CN + OH-->CH2CN + H2O

The reaction of acetonitrile with hydroxyl has been studied using the direct ab initio dynamics methods. The geometries, vibrational frequencies of the stationary points, as well as the minimum energy paths were computed at the BHandHLYP and MP2 levels of theory with the 6-311G(d, p) basis set. The energies were further refined at the PMP4/6-311+G(2df, 2pd) and QCISD(T)/6-311+G(2df, 2pd) levels of theory based on the structures optimized at BHandHLYP/6-311G(d, p) and MP2/6-311G(d, p) levels of theory. The Polyrate 8.2 program was employed to predict the thermal rate constants using the canonical variational transition state theory incorporating a small-curvature tunneling correction. The computed rate constants are in good agreement with the available experimental data.     

TDDFT investigation on nucleic acid bases: comparison with experiments and standard approach

A comprehensive theoretical study of electronic transitions of canonical nucleic acid bases, namely guanine, adenine, cytosine, uracil, and thymine, was performed. Ground state geometries were optimized at the MP2/6-311G(d,p) level. The nature of respective potential energy surfaces was determined using the harmonic vibrational frequency analysis. The MP2 optimized geometries were used to compute electronic vertical singlet transition energies at the time-dependent density functional theory (TDDFT) level using the B3LYP functional. The 6-311++G(d,p), 6-311(2+,2+)G(d,p), 6-311(3+,3+)G(df,pd), and 6-311(5+,5+)G(df,pd) basis sets were used for the transition energy calculations. Computed transition energies were found in good agreement with the corresponding experimental data. However, in higher transitions, the Rydberg contaminations were also obtained. The existence of pisigma* type Rydberg transition was found near the lowest singlet pipi* state of all bases, which may be responsible for the ultrafast deactivation process in nucleic acid bases.     

Excited state proton transfer in guanine in the gas phase and in water solution: a theoretical study

Theoretical investigations were performed to study the phenomena of ground and electronic excited state proton transfer in the isolated and monohydrated forms of guanine. Ground and transition state geometries were optimized at both the B3LYP/6-311++G(d,p) and HF/6-311G(d,p) levels. The geometries of tautomers including those of transition states corresponding to the proton transfer from the keto to the enol form of guanine were also optimized in the lowest singlet pipi* excited state using the configuration interaction singles (CIS) method and the 6-311G(d,p) basis set. The time-dependent density function theory method augmented with the B3LYP functional (TD-B3LYP) and the 6-311++G(d,p) basis set was used to compute vertical transition energies using the B3LYP/6-311++G(d,p) geometries. The TD-B3LYP/6-311++G(d,p) calculations were also performed using the CIS/6-311G(d,p) geometries to predict the adiabatic transition energies of different tautomers and the excited state proton transfer barrier heights of guanine tautomerization. The effect of the bulk aqueous environment was considered using the polarizable continuum model (PCM). The harmonic vibrational frequency calculations were performed to ascertain the nature of potential energy surfaces. The excited state geometries including that of transition states were found to be largely nonplanar. The nonplanar fragment was mostly localized in the six-membered ring. Geometries of the hydrated transition states in the ground and lowest singlet pipi* excited states were found to be zwitterionic in which the water molecule is in the form of hydronium cation (H3O(+)) and guanine is in the anionic form, except for the N9H form in the excited state where water molecule is in the hydroxyl anionic form (OH(-)) and the guanine is in the cationic form. It was found that proton transfer is characterized by a high barrier height both in the gas phase and in the bulk water solution. The explicit inclusion of a water molecule in the proton transfer reaction path reduces the barrier height drastically. The excited state barrier height was generally found to be increased as compared to that in the ground state. On the basis of the current theoretical calculation it appears that the singlet electronic excitation of guanine may not facilitate the excited state proton transfer corresponding to the tautomerization of the keto to the enol form.     

乙醇-水分子团簇C2H5OH(H2O)n (n=1-9)稳定结构的量子化学研究

采用密度泛函理论B3LYP方法, 在B3LYP/6-311++G(2d,2p)//B3LYP/6-311++G(d,p)基组水平上对乙醇-水分子团簇(C₂H₅OH(H₂O)_n (n=1-9))的各种性质进行研究, 如: 优化的几何构型、结构参数、氢键、结合能、平均氢键强度、自然键轨道(NBO)电荷分布、团簇的生长规律等. 结果表明, 从二维(2-D)环状结构到三维(3-D)笼状结构的过渡出现在n=5的乙醇-水分子团簇中. 此外, 利用团簇结合能的二阶差分、形成能、能隙等性质, 发现在n=6时乙醇-水分子团簇的最低能量结构稳定性较好, 可能为幻数结构. 最后, 为了进一步探讨氢键本质, 将C₂H₅OH(H₂O)_n (n=2-9)最低能量结构的各种性质与纯水分子团簇(H₂O)_n (n=3-10)比较, 结果表明前者与后者中的水分子之间氢键相似.     

水助甲酰胺-H2O2氧化乙烯制环氧乙烷反应机理的理论研究

采用MP2方法研究了甲酰胺-H2O2氧化乙烯制取环氧乙烷的反应机理.优化得到了反应物、过渡态、中间体及产物的几何构型并计算了反应势垒.研究结果表明:没有水参与时,反应需要通过四元环过渡态完成,反应势垒很高,在常温下难以进行;有水参与时,在水分子的协助下,反应可以通过六元环过渡态完成,反应势垒较低,常温下反应容易进行.     

Theoretical study on the hydrogen abstraction reaction of CH3CH2F (HFC-161) with Cl atom

A direct dynamics method is employed to study the hydrogen abstraction reaction of CH₃CH₂F+Cl. Three distinct transition states are located, one for -H abstraction and two for β-H abstraction. The potential energy surface (PES) information is obtained at the QCISD(T)/6-311+G(3df,2p)//MP2/6-311G(d,p), CCSD(T)/6-311+G(3df,2p)//MP2/6-311G(d,p) and G2//MP2/6-311G(d,p) level. Based on the QCISD(T)/6-311+G(3df,2p)//MP2/6-311G(d,p) results, the rate constants of the three reaction channels are evaluated by using the canonical variational transition state theory (CVT) with small-curvature tunneling (SCT) contributions over the temperature range of 220–2800 K. The calculated results indicate that -H abstraction dominates the total reaction almost over the whole temperature range.     

A computational study of conformational interconversions in 1,4-dithiacyclohexane (1,4-dithiane)

Ab initio molecular orbital theory with the 6-31G(d), 6-31G(d,p), 6-31+G(d), 6-31+G(d,p), 6-31+G(2d,p), 6-311G(d), 6-311G(d,p), and 6-311+G(2d,p) basis sets and density functional theory (BLYP, B3LYP, B3P86, B3PW91) have been used to locate transition states involved in the conformational interconversions of 1,4-dithiacyclohexane (1,4-dithiane) and to calculate the geometry optimized structures, relative energies, enthalpies, entropies, and free energies of the chair and twist conformers. In the chair and 1,4-twist conformers the C-Hax and C-Heq bond lengths are equal at each carbon, which suggest an absence of stereoelectronic hyperconjugative interactions involving carbon-hydrogen bonds. The 1,4-boat transition state structure was 9.53 to 10.5 kcal/mol higher in energy than the chair conformer and 4.75 to 5.82 kcal/mol higher in energy than the 1,4-twist conformer. Intrinsic reaction coordinate (IRC) calculations showed that the 1,4-boat transition state structure was the energy maximum in the interconversion of the enantiomers of the 1,4-twist conformer. The energy difference between the chair conformer and the 1,4-twist conformer was 4.85 kcal/mol and the chair-1,4-twist free energy difference (deltaG degrees (c-t)) was 4.93 kcal/mol at 298.15 K. Intrinsic reaction coordinate (IRC) calculations connected the transition state between the chair conformer and the 1,4-twist conformer. This transition state is 11.7 kcal/mol higher in energy than the chair conformer. The effects of basis sets on the 1,4-dithiane calculations and the relative energies of saturated and unsaturated six-membered dithianes and dioxanes are also discussed.     

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.     

Modern valence-bond description of chemical reaction mechanisms: the 1,3-dipolar addition of diazomethane to ethene

The electronic mechanism for the gas-phase concerted 1,3-dipolar cycloaddition of diazomethane (CH2N2) to ethene (C2H4) is described through spin-coupled (SC) calculations at a sequence of geometries along the intrinsic reaction coordinate obtained at the MP2/6-31G(d) level of theory. It is shown that the bonding rearrangements occurring during the course of this reaction follow a heterolytic pattern, characterized by the movement of three well-identifiable orbital pairs, which are initially responsible for the pi bond in ethene and the C-N pi bond and one of the N-N pi bonds in diazomethane and are retained throughout the entire reaction path from reactants to product. Taken together with our previous SC study of the electronic mechanism of the 1,3-dipolar cycloaddition of fulminic acid (HCNO) to ethyne (C2H2) (Theor. Chim. Acc. 1998, 100, 222), the results of the present work suggest strongly that most gas-phase concerted 1,3-dipolar cycloaddition reactions can be expected to follow a heterolytic mechanism of this type, which does not involve an aromatic transition state. The more conventional aspects of the gas-phase concerted 1,3-dipolar cycloaddition of diazomethane to ethene, including optimized transition structure geometry, electronic activation energy, activation barrier corrected for zero-point energies, standard enthalpy, entropy and Gibbs free energy of activation, have been calculated at the HF/6-31G(d), B3LYP/6-31G(d), MP2/6-31G(d), MP2/6-31G(d,p), QCISD/6-31G(d) and CCD/6-31G(d) levels of theory. We also report the CCD/6-311++G(2d, 2p)//CCD/6-31G(d), MP4(SDTQ)/6-311++G(2d,2p)//CCD/6-31G(d) and CCSD(T)/6-311++G(2d, 2p)//CCD/6-31G(d) electronic activation energies.     

类硅烯H2C＝SiLiBr与RH (R＝F, OH, NH2)的插入反应

采用DFT B3LYP和QCISD方法研究了类硅烯H2C＝SiLiBr与RH (R＝F, OH, NH2)的插入反应. 在B3LYP/6- 311＋G(d,p)水平上优化了反应势能面上的驻点构型. 结果表明, H2C＝SiLiBr与HF, H2O或NH3发生插入反应的机理相同. QCISD/6-311＋＋G(d,p)//B3LYP/6-311＋G(d,p)计算的三个反应的势垒分别为148.62, 164.42和165.07 kJ•mol－1, 反应热分别为－69.63, －43.02和－28.27 kJ•mol－1. 相同条件下发生插入反应时, 反应活性都是H—F＞H—OH＞H—NH2.