Ca+催化的N2O+CO反应机理的理论研究

采用B3LYP/cc-pVTZ理论水平系统研究了Ca+离子催化N2O+CO→N2+CO2反应的微观机理.反应分两步进行:第一步Ca+夺取N2O中的O原子有两条反应通道,其中优势通道为Ca+金属离子与N2O分子中O作用,形成线性分子复合物,活化N2O分子中的N-O键,之后的反应路径为O-N键断裂机理;第二步为CaO+金属...     

气相中O3与HSO自由基间的氢键复合物

气相中O3与HSO自由基之间的相互作用及其反应在大气化学中非常重要. 在DFT-B3LYP/6-311++G**和MP2/6-311++G**水平上求得O3+HSO复合物势能面上的稳定构型, B3LYP方法得到了三种构型(复合物I, II和III), 而MP2方法只能得到一种构型(复合物II). 在复合物I和III中, HSO单元中的1H原子作为质子供体, 与O3分子中的端基O原子作为质子受体相互作用, 形成红移氢键复合物; 而在复合物II中, 虽与复合物I和III中具有相同的质子供体和质子受体, 却形成了蓝移氢键复合物. B3LYP/6-311++G**水平上计算的单体间相互作用能的计算考虑了基组重叠误差(BSSE)和零点振动能(ZPVE)校正, 其值在-3.37到-4.55 kJ·mol-1之间. 采用自然键轨道理论(NBO)对单体间相互作用的本质进行了考查, 并通过分子中原子理论(AIM)分析了三种复合物中氢键的电子密度拓扑性质.     

气相中O3与HSO自由基间的氢键复合物

气相中O3与HSO自由基之间的相互作用及其反应在大气化学中非常重要.在DFT-B3LYP/6-311++G**和MP2/6-311++G**水平上求得O3+HSO复合物势能面上的稳定构型,B3LYP方法得到了三种构型(复合物Ⅰ,Ⅱ和Ⅲ),而MP2方法只能得到一种构犁(复合物Ⅱ).在复合物Ⅰ和Ⅲ中,HSO单元中的1H原子作为质子供体.与O3分子中的端基O原子作为质子受体相互作用,形成红移氢键复合物;而在复合物Ⅱ中,虽与复合物Ⅰ和Ⅲ中具有相间的质子供体和质子受体,却形成了蓝移氢键复合物.B3LYP/6-311++G**水平上计算的单体间相互作用能的计算考虑了基组重甍误差(BSSE)和零点振动能(ZPVE)校正,其值在-3.37到-4.55 kJ·mol-1之间.采用自然键轨道理论(NBO)对单体间相互作用的本质进行了考查,并通过分子中原子理论(AIM)分析了三种复合物中氢键的电子密度拓扑性质.     

复合物C2H6·(H2O)2中分子间相互作用的理论研究

在B3LYP/6-311 G**水平上得到C2H6.(H2O)2复合物势能面上四种稳定构型。在相同基组下经MP2电子相关能和基组叠加误差(BSSE)进行单点能量校正,求得单体间相互作用能的大小。结果发现:四种稳定构型都通过CH…O氢键而形成,相应σ(CH)键都出现了较小的收缩,导致伸缩振动发生蓝移,在最稳定的复合物Comp lex2和Comp lex3中,H2O(A)分子的一个H原子与C2H6的两个H原子相对距离较短,并且具有最大的总相互作用能和两个单体AC分子间相互作用能,这说明三个氢原子间存在着较强的相互作用,并对分子的稳定性起着重要作用.     

SiH2自由基与HNCO反应机理的理论研究

采用密度泛函理论B3LYP方法研究了SiH₂自由基与HNCO的反应机理, 并在B3LYP/6-311＋＋G**水平上对反应物、中间体、过渡态进行了全几何参数优化, 通过频率分析和内禀反应坐标(IRC)确定了中间体和过渡态. 为了得到更精确的能量值, 又用QCISD(T)/6-311＋＋G**方法计算了在B3LYP/6-311＋＋G**水平优化后的各个驻点的相对能量. 计算结果表明SiH₂自由基与HNCO的反应有五条反应通道, 其中顺式反应通道SiH₂＋HNCO→IM3→ TS4→IM5→TS5→IM6→SiH₂NH＋CO反应能垒最低, 为主反应通道.     

SiH2自由基与HNCO反应机理的理论研究

采用密度泛函理论B3LYP方法研究了SiH₂自由基与HNCO的反应机理, 并在B3LYP/6-311＋＋G**水平上对反应物、中间体、过渡态进行了全几何参数优化, 通过频率分析和内禀反应坐标(IRC)确定了中间体和过渡态. 为了得到更精确的能量值, 又用QCISD(T)/6-311＋＋G**方法计算了在B3LYP/6-311＋＋G**水平优化后的各个驻点的相对能量. 计算结果表明SiH₂自由基与HNCO的反应有五条反应通道, 其中顺式反应通道SiH₂＋HNCO→IM3→ TS4→IM5→TS5→IM6→SiH₂NH＋CO反应能垒最低, 为主反应通道.     

H2CCF自由基与HNCO反应机理的理论研究

采用密度泛函理论的B3LYP方法, 在6-311＋＋G(d,p)基组水平上研究了H₂CCF自由基与HNCO的微观反应机理, 优化了反应过程中的反应物、中间体、过渡态和产物, 为了获得更精确的能量信息, 还在QCISD(T)/6-311＋＋G(d,p)基组水平上计算了各物质的能量．振动分析结果和IRC分析结果证实了中间体和过渡态的真实性, 计算所得的成键临界点电荷密度的变化也确认了反应过程．对于H₂CCF自由基与HNCO反应, 我们找到了六条可行的反应通道, 结果分析表明通道H₂CCF＋HNCO→IM3→TS5→H₂CCFH＋NCO控制步骤活化能最低, 是该反应的主要通道, 在此反应过程中有稳定的氢键复合物IM3生成, 还表现出氢原子迁移的反应特征．     

CH2NH与O(3P)反应的量子化学及电子密度拓扑研究

采用B3LYP、MP2(full)和 QCISD 三种方法在6-311G(d, p)和aug-cc-pVDZ基组水平上对三线态O(3P)原子与CH2NH(s)的反应进行了详细的理论研究. 采用B3LYP和MP2(full)方法对反应势能面上的各驻点进行了几何构型优化, 通过振动频率分析证实了过渡态的真实性, 内禀反应坐标(IRC)跟踪验证了过渡态与反应物和产物的连接关系, 用上述三种方法计算得到了各反应通道的反应势垒. 对反应过程中的一些重要点进行了电子密度拓扑分析研究. 研究结果表明, O(3P)原子进攻CH2NH(s)中的N2原子和C1原子生成CH2NHO(t)和OCH2NH(t), CH2NHO(t)中N2上的H5可迁移到C1上异构化为CH3NO(t).     

H2CCF自由基与HNCO反应机理的理论研究

采用密度泛函理论的B3LYP方法, 在6-311＋＋G(d,p)基组水平上研究了H₂CCF自由基与HNCO的微观反应机理, 优化了反应过程中的反应物、中间体、过渡态和产物, 为了获得更精确的能量信息, 还在QCISD(T)/6-311＋＋G(d,p)基组水平上计算了各物质的能量．振动分析结果和IRC分析结果证实了中间体和过渡态的真实性, 计算所得的成键临界点电荷密度的变化也确认了反应过程．对于H₂CCF自由基与HNCO反应, 我们找到了六条可行的反应通道, 结果分析表明通道H₂CCF＋HNCO→IM3→TS5→H₂CCFH＋NCO控制步骤活化能最低, 是该反应的主要通道, 在此反应过程中有稳定的氢键复合物IM3生成, 还表现出氢原子迁移的反应特征．     

N2O…HOCl复合物结构和性质的理论研究

分别在DFT-B3LYP和MP2/6-311++G**水平上求得HOCl + N2O复合物势能面上的六种(S1, S2, S3, S4, S5和S6)和四种(S1, S2, S4和S5)构型. 频率分析表明,其中的S1和S3为过渡态,其它为稳定构型. 在复合物S3, S5 和S6中,HOCl 单体的σ*(5O-6H)作为质子供体,与N2O单体中作为质子受体的3O原子相互作用,形成氢键结构,而在氢键复合物S2中, 质子受体为N2O单体中的端1N原子;复合物S1中,HOCl分子的σ*(5O-4Cl)作为质子供体与N2O分子中的端1N原子(质子受体)相互作用,形成卤键结构,而复合物S4中的卤键结构的质子受体为N2O分子中的端3O原子. 经B3LYP/6-311++G**水平上的计算,考虑了基组重叠误差(BSSE)校正的单体间相互作用能在-1.56 ~ -8.73 kJ·mol-1之间. 采用自然键轨道理论(NBO)对两种单体间相互作用的本质进行了考查,并通过分子中原子理论(AIM)分析了复合物中氢键和卤键键鞍点处的电子密度拓扑性质.     

亚甲基硅烯与乙烯环加成反应机理的理论研究

The mechanism of a cycloaddition reaction between singlet methylidenesilene and ethylene has been investigated with MP2/6-31G^＊ and B3LYP/6-31G^＊ methods, including geometry optimization and vibrational analysis for the involved stationary points on the potential energy surface. Energies of the involved conformers were calculated by CCSD（T）//MP2/6-31G＊ and CCSD（T）//B3LYP/6-31 G＊ methods, respectively. The results show that the dominant reaction pathway of the cycloaddition reaction is that a complex intermediate is firstly formed between the two reactants through a barrier-free exothermic reaction of 13.3 kJ/mol, and the complex is then isomefized to a four-membered ring product P2,1 via a transition state TS2.1 with a barrier of 32.0 kJ/mol.     

亚烷基锡烯与乙烯环加成反应机理的理论研究

The mechanism of a cycloaddition reaction between singlet alkylidenestannylene and ethylene has been investigated with MP2/3-21 G^* and B3LYP/3-21 G* methods, including geometry optimization and vibrational analysis for the involved stationary points on the potential energy surface. Energies for the involved conformations were calculated by CCSD(T)//MP2/3-2 IG^* and CCSD(T)//B3LYP/3-21G^* methods, respectively. The results show that the dominant reaction pathway of the cycloaddition is that an intermediate (INT) is firstly formed between the two reactants through a barrier-free exothermic reaction of 39.7 kJ/mol, and the intermediate then isomerizes to a four-membered ring product (P2.1) via a transition state TS2.1 with a barrier of 66.8 kJ/mol.     

Density Functional Theory and MP2 Calculations of the Transition States and Reaction Paths on Coupling Reaction of Methane through Plasma

The two possible reaction paths of producing ethane on coupling reaction of methane through plasma were theoretically investigated by B3LYP and MP2 methods with 6-311G^* respectively and further compared with the previous results calculated from B3LYP/6-31G^*. The new investigated results consistently confirmed the previous conclusion. And the influences of the calculation methods and basis sets on the calculated results were also discussed.     

Comparative density functional theory study of the binding of ligands to Cu+ and Cu2+: Influence of the coordination and oxidation state

BP86, B3LYP and MP2 methods, generally used to study large systems containing transition metals, were compared for their ability to accuratly evaluate bond dissociation energies of copper complexes. Various [Cu-L]+ and [Cu-L]2+ complexes in which L are small ligands and the higher coordinated complexes, [Cu(NH3)(4)]+ and [Cu(NH3)4]2+ were studied. For monoligated complexes, the BDEs calculated by the three methods differed by 2 to 60 kcal/mol, the larger differences being obtained for [Cu-L]2+ complexes. The BDEs calculated using the B3LYP functional were in general close to the experimental values whereas the BDEs calculated using the BP86 functional were too high and the BDEs calculated using the MP2 were too low. If we rank the whole ligands according to their increased bond strength, the resulting orders obtained with the three methods are different for the [Cu-L]+ complexes, the B3LYP giving the same order as the experimental one. This result indicates that the BDEs of [Cu-L]+ complexes are better modeled using the B3LYP than using the BP86 and MP2 methods. For [Cu-L]2+, B3LYP also gave the most reliable results whereas BP86 gave too large BDEs and MP2 gave too small BDEs. However, symmetries of ground states can be different using DFT and post-Hartree-Fock methods. For [Cu-N2O]2+ the use of the B1LYP provides a better symmetry of the complex than the B3LYP, as has been recently shown in the literature for [Cu-H2O]2+. MP2 led to an incorrect bent structure for [Cu-N2]2+ in contrast to a linear structure obtained with the other methods, including CCSD(T). However, due to the lack of experimental data for [Cu-L]2+ complexes and to contrasted results for the methods, it is not possible to conclude definitely. For the high coordinated complexes [Cu(NH3)4]+ and [Cu(NH3)4]2+, the PBE calculation method was used in addition to the BP86, B3LYP and MP2. The BDE values were very close to each other when there is no change of the oxidation state during the reaction. On the basis of these calculations, the choice of the method was less crucial for high coordinated complexes [Cu(NH3)4]+ and [Cu(NH3)4]2+ so long as the oxidation state remained the same during the reaction. In contrast, when [Cu(NH3)4]2+ is reduced in [Cu(NH3)3]+ and NH3, the BDE calculated using the four methods were markedly different.     

反应类等键反应方法及类反应热动力学参数计算

本文提出了反应类等键反应方法,将通常用于热力学性质计算的等键反应方法推广用于类反应中反应势垒和反应焓变的计算. 对碳氢燃料低温燃烧反应机理中的一类重要反应类：?烷基自由基过氧化氢裂解生成烯烃和HO2自由基的反应势垒和反应焓变进行了计算. 通过对该类16个反应中的5个代表反应分别在不同计算水平HF、DFT、MP2、CCSD(T)的比较计算发现,采用等键反应方法可在较低从头算级别计算得到类反应较高精度的反应势垒,提高了计算的效率和精度. 本文采用反应类等键反应方法在B3LYP/6-311G(d,p)计算水平对该类16个反应进行了反应势垒和反应焓变的计算,并建立了反应势垒和反应能的线性自由能关系：delta V=71.02+0.41?delta E.     

二氯亚甲基锗烯与乙烯环加成反应机理的理论研究

The mechanism of a cycloaddition reaction between singlet dichloromethylene germylene and ethylene has been investigated with B3LYP/6-31G＊ method, including geometry optimization and vibrational analysis for the involved stationary points on the potential energy surface. Energies for the involved conformations were calculated by CCSD（T）//B3LYP/6-31G＊ method. On the basis of the surface energy profile obtained with CCSD（T）// B3LYP/6-31G＊ method for the cycloaddition reaction between singlet dichloromethylene germylene and ethylene, it can be predicted that the dominant reaction pathway is that an intermediate INT1 is firstly formed between the two reactants through a barrier-free exothermic reaction of 61.7 kJ/mol, and the intermediate INT1 then isomerizes to an active four-membered ring product P2.1 via a transition state TS2, an intermediate INT2 and a transition state TS2.1, in which energy barriers are 57.7 and 42.2 kJ/mol, respectively.     

Cl_2 2HI=2HCl I_2反应机理的理论研究

分别在MP2/3-21G!!、CCSD(T)/3-21G!!//MP2/3-21G!!和B3LYP/3-21G!!3种水平上,计算研究了气相反应Cl2 2HI=2HCl I2的机理,求得一系列四中心和三中心的过渡态.通过比较六种反应通道的活化能大小,得到了相同的结论:双分子基元反应Cl2 HI"HCl ICl和ICl HI"I2 HCl的最小活化能小于Cl2、HI和ICl的解离能,从理论上证明了反应Cl2 2HI=2HCl I2将优先以分子与分子作用形式分两步完成.用内禀反应坐标(IRC)验证了MP2/3-21G!!方法计算得到的过渡态.     

FC(O)O自由基与NO2反应的量子化学研究

用量子化学DFT, MP2, G3和G3MP2方法对FC(O)O自由基与NO₂的反应机理进行了理论研究. 优化了反应势能面上各驻点的几何结构, 通过内禀反应坐标(IRC)计算和振动分析, 确认了反应中的过渡态, 并用过渡态理论(TST)计算了相关反应的速率常数.     

Theoretical Study on the Mechanism of CF2 Reaction with CH2O

The insertion reaction mechanism of CF2 with CH2O was investigated at the B3LYP/6-311G（d）//MP2/6-311G（d） level. The geometric conformations at each stationary point in reaction potential surface were fully optimized and the transition states were verified by intrinsic reaction coordinate （IRC） and frequency analysis. The energies of all reactants were calculated with CCSD（T）/6-311G（d）//G2MP2 methods. Results indicated that the P1 reaction route with difuoroaldehyde as product is the dominant reaction pathway, which exhibits nucleophilic character. According to NBO analysis, the starting point of insertion reaction is the interaction between carbene LP（C3） and formaldehyde π（Cl-O2）. Besides, the thermodynamic and dynamic properties of dominated reaction （1） at different temperature were studied with statistic thermodynamic method and Eyring transition state theory adjusted by Wigner means, from which the proper temperature （500- 1200 K） of reaction （1） could be estimated. Finally, the thermo- dynamic and dynamic properties of insertion reaction mechanisms （CF2, CX2 （X = Cl, Br） with CH2O） were compared and discussed.