O（^3P）与CH2=CHF反应的理论研究

用量子化学密度泛函理论（DFT）方法,对O（^3P）与CH2=CHF的反应进行了理论研究.在UB3LYP/6-311＋＋G（d,p）,UB3LYP/6-31＋＋G（3df,3pd）计算水平上,优化了反应物、产物、中间体和过渡态的几何构型,并在UQC ISD（T）/6-311＋＋G（2df,2pd）水平上计算了单点能量.为了确证过渡态的真实性,在UB3LYP/6-311＋＋G（d,p）水平上进行了内禀反应坐标（IRC）计算和频率分析,并确定了反应机理,研究结果表明,反应主要产物为CH2CFO和H.     

O(1D)与CF3Cl反应的理论研究

用量子化学密度泛函理论（DFT）方法，对O（^1D）与CF3CI的反应进行了研究．在B3LYP／6-311＋G（d），B3LYP／6-311＋G（2df，2pd）计算水平上，优化了反应物、产物、中间体和过渡态的几何构型。并在QCISD（T）／6-311＋＋G（d，p）水平上计算了单点能量，为了确证过渡态的真实性，在B3LYP／6-311＋G（d）水平上进行了振动分析．研究结果表明，O（^1D）与CF3Cl反应中存在4个反应通道，其中生成产物：P1（CF3O＋Cl）和P2（CF3＋CIO）的通道应为主反应通道．     

反应CH3CN+OH→CH2CN+H2O的动力学

采用Gaussian 98从头算程序包,从理论上对反应CH3CN+OH→CH2CN+H2O的动力学进行了研究.反应物、产物和过渡态的几何优化分别在UMP2/6-311+G(d,p)和UB3LYP/6-311+G(d,p)水平下进行.用MP4(SDTQ)/6-311+G(d,p)和CCSD(T)/6-311G(3df,2p)方法对UMP2/6-311+G(d,p)优化的反应物、产物和过渡态进行了单点能校正,并计算出了反应在300～2 100 K温度范围的速率常数和谐振频率.     

Cl原子与CH3OCH2Br反应机理的量子化学研究

对C1与CH3OCH2Br的反应进行了量子化学研究．在MP2／6-311G（d，p）水平上优化了反应势能面上各驻点的几何结构，并在QCISD（T）／6-311G（d，p）水平上计算了各驻点的能量．通过内禀反应坐标（IRC）计算，对反应过渡态进行了确认，并确定了反应的机理．     

CH(X2Ⅱ)与NH3反应机理的量子化学研究

对CH(X2Ⅱ)自由基与NH3的反应进行了量子化学研究.分别在B3LYP/6-311++G(d,p),B3LYP/6-311++G(3df,3pd),MP2/6-311++G(d,p)和MP2/6-311++G(3df,3pd)水平优化了反应势能面上各驻点的几何结构,并在QCISD(T)/6-311++G(3df,3pd)水平上计算了各驻点的能量.通过IRC计算确认了过渡态,确定了反应的机理.     

CO与HO2或NO2反应机理的计算研究

用密度泛函DFT（B3LYP）方法和MP2从头算方法,在6-311G（d,p）,6-311＋＋G（d,p）基组水平上研究了CO＋HO2-C02＋OH,CO＋NO2-CO2＋NO2个反应的过渡态结构和反应机理,沿IRC分析指出反应均为协同完成的一步反应．同时计算了2个反应的活化位垒和反应焓变,并与实验估计值进行了比较．对CO与H02的反应采用MP2／6—311＋＋G（d,p）方法,计算所得反应位垒为100．98kJ／mol,与实验估计值93．60kJ／mol基本接近,对CO与N02的反应采用B3LYP／6—311＋＋G（d,p）方法,计算所得反应位垒为115．99kJ／mol,比实验估计值141．51kJ／mol低25．52kJ／mol．计算同时表明2个反应均为放热反应．     

OXO(X=Cl,Br)自由基与CO和NO自由基反应机理的理论研究

对CO和NO自由基与OXO(X=Cl,Br)自由基反应机理进行了理论研究.在MP2/6-311 +G(d,P)水平下优化了反应物、产物、中间体和过渡态的几何构型,并对所有平衡构型进行了简谐振动频率分析.在相同的计算水平下以过渡态为出发点,利用内禀反应坐标(IRC)理论得到了反应的最小能量路径(MEP),能量信息的进一步校正在较高的MCQCISD//MP2/6-311+G(d,p)水平下完成.研究结果表明,该反应体系可发生氧抽提反应,其中NO与OClO的反应是最具有动力学优势的.     

O(3P)与CH2=CHF反应的理论研究

用量子化学密度泛函理论(DFT)方法,对O(3P)与CH2=CHF的反应进行了理论研究.在UB3LYP/6-311 G(d,p),UB3LYP/6-31 G(3df,3pd)计算水平上,优化了反应物、产物、中间体和过渡态的几何构型,并在UQCISD(T)/6-311 G(2df,2pd)水平上计算了单点能量.为了确证过渡态的真实性,在UB3LYP/6-311 G(d,p)水平上进行了内禀反应坐标(IRC)计算和频率分析,并确定了反应机理,研究结果表明,反应主要产物为CH2CFO和H.     

O(3 P)与C2H2反应机理的量子化学研究

用量子化学计算方法对O(^3P)与C2H2的反应进行了研究．在HF／6—311 G(d，p)，HF／6—311 G(3df，3pd)，MP2／6—311 G(d，p)和MP2／6—311 G(3df，3pd)计算水平上，优化了反应势能面上各驻点的几何结构．通过内禀反应坐标(IRC)计算和振动分析．对反应过渡态进行了确认，并确定了反应机理．     

二氯卡宾与异硫氰酸反应机理的量子化学研究

用密度泛函理论B3LYP方法研究了单、三重态CCl2与HNCS的反应机理.在B3LYP/6-311＋＋G＊＊水平上对反应物、中间体和过渡态进行了全几何参数优化;通过频率分析和内禀反应坐标（IRC）方法确认了中间体和过渡态;并用CCSD（T）/6-311＋＋G＊＊//B3LYP/6-311＋＋G＊＊方法计算了各个驻点的单点能.结果表明,对于单、三重态CCl2与HNCS反应,共有4条可行的反应通道,且通道CCl2（S）＋HNCS→IM1→TS1→IM2→HCN＋CCl2S（P1）反应能垒最低,为该反应的主反应通道;单重态抽提亚氨基通道反应能垒也比较低,为主反应通道的竞争通道.可为进一步实验研究提供参考.     

从头算和RRKM理论研究O（1D）+trans-CH2FCH2F反应

运用MP2/6-311G(d,p)方法优化得到了O(1D)+trans-CH2FCH2F反应势能面上稳定点的优化几何构型和振动频率。精算的能量是在QCISD(T)/6-311++G(d,p)水平上得到的。RRKM(Rice-Ramsperger-Kassel-Marcus)理论用来计算碰撞能在0-350 kcal/mol范围内通过中间体的微正则速率常数。计算得出,当碰撞能低于140 kcal/mol时标题反应的主要产物是HF和CH2F;碰撞能在140-350 kcal/mol时反应的主要产物是HF,CH2F和OH.在整个碰撞能范围内,HF和CH2F是主要产物。     

甲酸分子内质子转移过程的分子轮廓的动态变化

在MP2/6-311+G(d,p)水平下,获得甲酸分子内质子转移反应的过渡态,并在MP2/6-31G(d)水平下,研究了该反应的反应路径(IRC).沿着IRC反应路径选取了9个特征点,使用MELD精密从头计算中的CISD/3-21+G(d,p)方法,结合我们自编的分子形貌程序,绘制了甲酸的分子形貌以及其分子界面的特征轮廓,定义和计算了分子形貌轮廓参数[包括分子形貌轮廓(R)以及电子密度(ED)],比较了各原子端轮廓值与孤立原子半径的相对大小,反映了甲酸分子形貌轮廓的动态变化,以及该反应过程中电子密度的变化.生动形象地展现了甲酸分子内质子转移过程中的空间和前沿电子密度的动态变化信息,清晰地再现了电子的转移过程.为该反应,提供了一种微观认识.     

脯氨酸-水二聚体分子间相互作用的理论研究

在B3LYP/6-311＋G（d,p）水平下优化得到了脯氨酸-水二聚体势能面上的4种低能优化构型。利用MP2/6-311＋G（d,p）方法计算脯氨酸-水二聚物平衡结构的偶极矩,原子电荷布居和分子间相互作用能。经BSSE校正后的相互作用能依次增强的顺序为Ⅳ〈Ⅱ〈Ⅰ〈Ⅲ,最稳定结构Ⅲ的相互作用能为-33.74 KJ/mol。     

二氯卡宾消耗臭氧的反应机理研究

采用密度泛函理论对单重态势能面上的二氯卡宾1 CCl2与臭氧1 O3反应的微观机理进行了研究.在B3LYP/6-311G(d,p)水平上全参数优化反应物、中间体、过渡态和产物的几何构型,同时使用内禀反应坐标(IRC)在同一水平上对过渡态与中间体之间的联系进行了确认,并且在QCISD/6-311G(d,p)水平上计算了各驻点的单点能.研究得到1 CCl2与1 O3反应的8种产物通道,即P1(Cl2CO+1 O2)、P2(CO2+ClClO)、P3(CO3+Cl2)、P4(CO+ClOOCl)、P5(CO2+ClOCl)、P6(v-CO2+ClOCl)、P7(ClCClO+1 O2)和P8(COO+ClOCl),其中通道P1(Cl2CO+1 O2)是最主要的产物通道,通道P2(CO2+ClClO)、P3(CO3+Cl2)、P4(CO+ClOOCl)对反应体系的产物贡献依次减小,而通道P5(CO2+ClOCl)、P6(v-CO2+ClOCl)、P7(ClCClO+1 O2)和P8(COO+ClOCl)很难生成.     

全金属团簇Be32-芳香性的理论研究

将芳香性概念扩展到全金属团簇Be32-,运用密度泛函理论（B3LYP,B3PW91）,从头算方法（MP2）,对全金属团簇Be32-,NaBe3-和Na2Be3的稳定结构、振动频率与电子总能量（考虑了零点能ZPE）作理论计算．计算的结果显示,NaBe3-和NaBe32-团簇中包含一个正三角形的阴离子Be32-．并根据芳香性的平面、电子结构、核独立化位移（NICS）以及它们的分子轨道几个标准进行分析．分析的结果指出,2个离域化的π电子遵循4n＋2电子计算规则,并且呈现出π芳香性．     

水合锌（Ⅱ）团簇Zn2+（H2O）n（n=10~12）的理论化学研究

应用高水平量子化学方法 MP2/6-311++G(2d,2p)//M062X/6-311+G(d,p)对气相中的Zn2+(H2O)n(n=10~12)团簇进行优化,获得其多种稳定结构,计算并比较不同Zn2+的配位数下的结合能及电荷分布.结果表明:当n=10~12时,Zn2+的配位数可以为4、5或6,其中最稳定结构的Zn2+配位数为5或6.随着水分子数的增加,体系最稳定结构的氢键数和离子与内层水分子中氧的距离RZn-O都在不断增加.     

铍氯类锗烯H_2GeClBeCl与RH(R=OH,NH_2,CH_3)的插入反应

用DFT B3LYP和QCISD方法研究了铍氯类锗烯H2GeClBeCl与RH(R=OH,NH2,CH3)的插入反应.在B3LYP/6-311+G(d,p)水平上优化了反应过程中所有驻点的构型并用QCISD/6-311+G(d,p)方法计算了单点能量,并考察了溶剂化效应对反应的影响.结果表明,在插入反应势能面上有一个过渡态(TS)和一个中间体(IM)连接反应物和产物.计算的反应势垒分别为177.62(R=OH),186.30(R=NH2),214.90(R=CH3)kJ/mol,表明在相同反应条件下,反应活性大小为H-OH,H-NH2,H-CH3,随着溶剂极性的增大,反应越来越容易进行.     

密度泛函理论研究叶绿素a水合物的结构、结合能和电离能

用M06-2X/6-311G(d,p)密度泛函理论对叶绿素a及其两种水合物的结构进行优化,然后用M06-2X/6-311++G(2d,2p)方法计算结合能和电离能。优化结果表明水合分子会使叶绿素a分子结构发生改变。结合能计算结果表明叶绿素a的二水合物比一水合物更为稳定。电离能计算结果表明随着叶绿素a水合分子数的增加,电离能逐渐减小,叶绿素a的水合有利于光合作用的发生。     

几种醇类阳离子的密度泛函理论研究

用B3LYP方法及6—31G（d,P）、6—31＋G（d,p）、6—311G（d,P）和6—311＋G（d,P）基组,对若干醇类阳离子CH3OH＋、CH3CH2OH＋和CH3CH2CH2OH＋做了理论研究,优化了它们的电子基态的结构,计算了对应分子的垂直电离势（VIP）和绝热电离势（AIP）．结果表明,这些离子的对称性与对应分子相同;自然布居分析计算表明,这些离子的正电荷主要分布在各H原子上;B3LYP／6—3ll＋G（d,P）级别上计算的各分子的VIP值和实验值符合得很好．     

NCO自由基与H反应的机理和动力学

研究采用密度泛函理论中的B3LYP方法,在6-311G(d,p)基组水平下,优化反应物、产物、中间异构体和过渡态分子,得到可靠的几何构型和频率.通过内禀反应坐标计算,确认了四条反应通道,并且以通道中中间体、过渡态相对能量为基础,讨论了各反应通道的优先顺序.在反应中,H和NCO的初始连接有两种方式,H分别进攻NCO两端的N和O,进行无势垒加合,得到两个低能中间体HNCO和HOCN.HNCO中∠NCO从172.9°减小到69.0°,生成一个相对能量较高但N—C键较长的中间体,该中间体发生N—C键断裂生成主要产物P1(CO+NH).从HNCO和HOCN出发均可生成次要产物P2(CN+OH),而从HNCO出发还有另一条生成产物P3(CH+NO)的路径,反应沿该路径进行的可能性较小.