CH自由基和NO~2反应研究: I. 反应的热力学计算

探讨了CH自由基与NO~2反应的可能路径,通过计算确定了反应物,产物和稳定中间体的电子状态和平衡构型,并运用Gaussian-3方法和MRCISD方法对可能的反应路径进行了热力学计算。在多数情况下与实验值符合较好。对于个别与理论计算差别较大的实验值进行了评述。     

CH自由基和NO2反应研究Ⅰ.反应的热力学计算

探讨了CH自由基与NO2反应的可能路径,通过计算确定了反应物,产物和稳定中间体的电子状态和平衡构型,并运用Gaussian-3方法和MRCISD方法对可能的反应路径进行了热力学计算.在多数情况下与实验值符合较好.对于个别与理论计算差别较大的实验值进行了评述.     

对两个SN2反应的从头计算研究

对S_N2反应势能面的研究是探讨有机和生物大分子反应机理的中心课题.对气相反应,精确的从头计算可得到与实验值相近的能垒值和实验无法确定的过渡态和反应路径等细     

碳前驱体CH3ArCH2NH2热解反应的热力学和动力学DFT研究

在实验研究基础上,通过量子化学理论计算对碳前驱体CH3ArCH2NH2的热裂解机理作了进一步的研究.利用Gaussian98程序包中AM1方法及DFTUB3LYP/3-21G*方法,对化合物5种可能热裂解路径的热力学和动力学计算结果表明,CH3ArCH2NH2热裂解的主反应路径为生成自由基CH3ArCH2*和NH2*,其主反应路径AM1计算的活化能Ea=230.78kJ/mol,DFT计算的活化能Ea=321.18kJ/mol;比较键焓计算的数据与相应的实验数据,发现DFT计算结果与实验结果吻合得较好;通过分析优化的反应物及产物自由基的部分结构参数,了解了理论支持主反应的原因;计算的产物自由基的空间构型表明主反应路径生成的产物自由基相互间若进行稠环缩合反应,将获得分子平面取向性很好的稠环芳烃产物.     

碳前驱体CH3ArCH2NH2热解反应的热力学和动力学DFT研究

在实验研究基础上 ,通过量子化学理论计算对碳前驱体 CH3 Ar CH2 NH2 的热裂解机理作了进一步的研究 .利用 Gaussian98程序包中 AM1方法及 DFT UB3 LYP/3 -2 1 G*方法 ,对化合物 5种可能热裂解路径的热力学和动力学计算结果表明 ,CH3 Ar CH2 NH2 热裂解的主反应路径为生成自由基 CH3 Ar CH2 ·和 NH2 · ,其主反应路径 AM1计算的活化能 Ea=2 3 0 .78k J/mol,DFT计算的活化能 Ea=3 2 1 .1 8k J/mol;比较键焓计算的数据与相应的实验数据 ,发现 DFT计算结果与实验结果吻合得较好 ;通过分析优化的反应物及产物自由基的部分结构参数 ,了解了理论支持主反应的原因 ;计算的产物自由基的空间构型表明主反应路径生成的产物自由基相互间若进行稠环缩合反应 ,将获得分子平面取向性很好的稠环芳烃产物     

烷基自由基β位裂解反应类反应势垒与速率常数的精确计算

提出反应类等键方法并用于高温燃烧机理中一类重要反应——烷基自由基β位裂解反应的反应势垒和速率常数的精确校正计算. 通过10种不同从头算水平对类反应中5个代表反应的反应势垒的计算发现, 用反应类等键反应方法和直接从头算方法获得的5 个代表反应的反应势垒最大绝对偏差的平均值分别为5.32 和16.16 kJ·mol^-1, 表明反应类等键反应方法计算的反应势垒对不同水平从头算方法的依赖性小, 可在较低从头算水平计算得到精确的反应势垒, 解决大分子体系反应势垒的精确计算问题. 此外应用反应类等键反应方法在BHandHLYP/cc-pVDZ 从头算水平计算了3 个代表反应的速率常数, 并与文献报道的实验值进行了比较, 其在500-2000 K温度区间内计算速率常数与实验速率常数中较大值与较小值的比值k_max/k_min的平均值为1.67, 最大值也仅有2.49. 表明应用反应类等键反应方法在较低从头算水平即可对同类反应的速率常数进行精确计算.最后在BHandHLYP/cc-pVDZ从头算水平用反应类等键反应方法计算了13个烷基自由基β位裂解反应的速率常数.     

CH自由基和NO~2反应研究: II.反应的动力学计算

采用MP2(FULL)/6-31G(d)方法从动力学计算上探讨了CH自由基与NO~2反应的可能途径,找到了反应物,中间体及产物之间的能量通道和过渡态,报道了它们的构型、电子态及能量。并通过频率分析和IRC方法对所有的过渡态进行了验证。在此基础上求出了各步反应的活化能。在以前热力学研究的基础上,对于可能的反应通道进一步作了动力学分析,找到了反应主产物通道的分支比,与实验得到的分支比基本吻合。     

氯原子与甲酸反应的反应路径和动力学研究

用从头算UMP2方法采用6-311+G基组详细研究了Cl原子与甲酸反应的最佳反应 路径，对所有驻点都进行了G2单点能校正，结果表明：在Cl与HC0OH反应的过程中 ，Cl原子主要夺取trans-HCOOH中C端的H原子，产生cis—HOCO基团，反应的势垒为 16．7kJ／mol．考虑Cl原子存在下的二级反应，cis—HOCO将通过通道Rla1和Rla2 生成trans-HOCOCl再分解为最终产物CO2.cis-HOCO和trans-HOCO自由基与Cl原子的 反应为放热和无势垒过程．而HOCO和HOCOCl两个中间体的cis—到trans-转化势垒 分别为15．1和33．5kJ/mol．由于构型的原因，HOCOCl分解成CO2和HCl必须通过 trans-构型来进行．这步反应的势垒较高，为131．0kJ/mol．298K时，采用经典过 渡态理论计算得到的引发步骤的H-迁移反应的速率常数与实验值符合很好．也计算 了甲酸二聚体与Cl原子的H-迁移反应的速率常数．     

CH自由基和NO~2反应研究: II.反应的动力学计算

采用MP2(FULL)/6-31G(d)方法从动力学计算上探讨了CH自由基与NO~2反应的可能途径,找到了反应物,中间体及产物之间的能量通道和过渡态,报道了它们的构型、电子态及能量。并通过频率分析和IRC方法对所有的过渡态进行了验证。在此基础上求出了各步反应的活化能。在以前热力学研究的基础上,对于可能的反应通道进一步作了动力学分析,找到了反应主产物通道的分支比,与实验得到的分支比基本吻合。     

吡啶光氯化反应过渡态和反应途径的量子化学研究

用量子化学B3LYP方法在3---21G*水平上优化吡啶光氯化反应加成取代反应机理生成邻、间、对位氯代吡啶不同反应途径的过渡态并对反应热和活化能进行了计算,对邻位反应途径进行了IRC反应路解析,计算结果表明邻位反应途径过渡态的能量最低,为－704.830027a.u.,生成2－氯吡啶所需的活化能最低,为114.60kJ/mol。光氯化反应主要产物为2－氯吡啶,与实验结果一致。IRC反应路径显示在反应过程中C(2)---H(7)键的断裂和C(2)-----Cl(8)键的生成是协同但不是同步的。     

CH自由基和NO2反应研究Ⅱ.反应的动力学计算

采用MP2(FULL)/6-31G(d)方法从动力学计算上探讨了CH自由基与NO2反应的可能途径,找到了反应物,中间体及产物之间的能量通道和过渡态,报道了它们的构型、电子态及能量.并通过频率分析和IRC方法对所有的过渡态进行了验证.在此基础上求出了各步反应的活化能.在以前热力学研究的基础上,对于可能的反应通道进一步作了动力学分析,找到了反应主产物通道的分支比,与实验得到的分支比基本吻合.     

基于量子化学计算的BH4-水解制氢反应机理研究

采用密度泛函理论(DFT)和二阶微扰理论(MP2)对BH-4与两分子水反应制氢的微观机理进行了研究.在B3PW91/6-311++G(2df,2p)和MP2(full)/6-311++G(2df,2p)水平上优化了反应体系中所有反应物、中间体、过渡态和产物的几何构型,通过振动频率分析和内禀反应坐标(IRC)跟踪验证了过渡态的正确性.计算结果表明,反应生成氢气的过程是分四步进行,并存在两种可能的反应路径:路径1中,两分子的水分别在第一步、第三步参加反应;而路径2中,两分子的水分别在第一步、第二步参加反应;其中路径2的势垒比路径1低应为主反应路径.此外,两种反应路径都是强放热过程,且MP2法所得焓变值与所报道的实验值更接近.     

HNCO+HCO→NCO+CH2O氢转移反应的从头算及动力学研究

在UMP2(Full)/6-311G(d,p)计算水平上,优化了标题反应的反应物、过渡态、产物的几何结构,沿最小能量途径讨论了异氰酸(HNCO)和甲酰自由基(HCO)发生氢转移反应位能面上驻点的结构以及相互作用分子结构变化.指出该反应是一个N-H键断裂和C-H键生成的协同反应.进一步采用UQCISD(T,Full)方法对反应途径上的驻点进行了单点能量校正,得出该反应的计算位垒是91.47 kJ/mol,与实验值108.92 kJ/mol接近在500～2500K实验温度范围内,运用变分过渡态理论(CVT)计算得到的速率常数与实验观测值进行了比较     

一类有机体系分子内电子转移反应的理论研究

对1984年Miller和Closs设计的有机分子内电子转移反应体系的反应速率常数进行了直接理论计算研究,获得与实验值接近的结果,理论计算结果进一步证明了体系中反转区的存在.     

电子转移反应O2+O2·-→O2·-+O2的从头算研究

利用abinitio方法,在UHF,UMP2及不同基组3-21G,6-31G^*,6-311+G^*和UMP2(full)/6-311+G^*水平上,研究了O~2/O~2^.^-自交换电子转移反应。优化了电子转移前后反应物和产物的结构,研究了体系能量的变化,计算了自交换电子转移反应的内重组能。对UHF方法和UMP2方法的计算结果进行了比较,并与实验结果进行了对照。结果表明UHF方法由于没有考虑组态相互作用,计算结果存在较大偏差,UMP2(full)/6-311+G^*水平上的计算结果与实验值吻合较好。在UMP2(full)/6-311+G^*水平上计算了气相自交换电子转移反应速率常数。在优化了电子转移复合物结构的基础上考虑了溶剂效应的影响,计算了水溶液中的溶剂重组能。研究结果表明O~2/O~2^.^-体系电子转移反应的活化能主要来源于溶剂重组能的贡献。最后计算了该反应在水溶液中的反应速率常数。理论计算结果与实验值吻合得很好。     

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

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

CH_3CN+CH_3→CH_2CN+CH_4氢消除反应的直接动力学研究

采用直接动力学方法,对乙腈与甲基的反应进行了理论研究.在BHandHLYP/6-311G(d,p)和MP2/6-311G(d,p)水平下获得,稳定点的几何结构、振动频率及最小能量路径(MEP),在G3(MP2)和MC-QCISD水平下对能量信息进一步确认.利用正则变分过渡态理论,结合小曲率隧道效应校正(CVT/SCT)方法计算了该反应在220K～2000K的速率常数,与实验值符合得很好.     

电子转移反应O~2+O~2^.^-→O~2^.^-+O~2的从头算研究

利用abinitio方法,在UHF,UMP2及不同基组3-21G,6-31G^*,6-311+G^*和UMP2(full)/6-311+G^*水平上,研究了O~2/O~2^.^-自交换电子转移反应。优化了电子转移前后反应物和产物的结构,研究了体系能量的变化,计算了自交换电子转移反应的内重组能。对UHF方法和UMP2方法的计算结果进行了比较,并与实验结果进行了对照。结果表明UHF方法由于没有考虑组态相互作用,计算结果存在较大偏差,UMP2(full)/6-311+G^*水平上的计算结果与实验值吻合较好。在UMP2(full)/6-311+G^*水平上计算了气相自交换电子转移反应速率常数。在优化了电子转移复合物结构的基础上考虑了溶剂效应的影响,计算了水溶液中的溶剂重组能。研究结果表明O~2/O~2^.^-体系电子转移反应的活化能主要来源于溶剂重组能的贡献。最后计算了该反应在水溶液中的反应速率常数。理论计算结果与实验值吻合得很好。     

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

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

H_20 CN→HCN OH的IRC解析及其反应动力学研究

本文用能量梯度法,在UHF/3-21G水平上,优化了反应H_2O CN→HCN OH的反应物、产物及其过渡态的几何构型,得到了该反应的活化能为32.6kJ/mol,与实验所得的测量值相一致;同时用Morokuma的数值分析方法,得到了该反应的内禀反应坐标(IRC);沿着IRC,对反应过程中体系构型的变化进行了分析,并计算IRC运动与垂直于IRC简正振动之间的偶合常数,讨论振动频率的变化,并结合偶合常数进行分子动态学分析;用传统过渡态理论、变分过渡态理论及相关的隧道校正等方法计算该反应的速率常数,结果与实验值基本吻合(如T=763K时,K_(计算值)~(CVT/SCSAG)=309×10~(10),K_(实验值)=(5.1±0.6)×10~(10),单位cm~3·mol~(-1)·s~(-1)).