鸟嘌呤受羟基自由基损伤反应机理的量子化学研究

用密度泛函理论(DFT)研究羟基自由基与鸟嘌呤分子加成反应的过渡态, 并进行内禀反应坐标(IRC)反应路径解析, 结果表明, 羟基自由基加成到鸟嘌呤碳碳双键上. 利用B3LYP/6-31++G**对反应物、反应物络合物、过渡态以及产物络合物等反应通道上各个能量驻点的能量进行了计算, 得到反应活化能E_a=28.0867 kJ/mol. AIM计算结果显示, 过渡态结构中鸟嘌呤分子碳碳双键结构被削弱, 羟基自由基氧原子与鸟嘌呤分子碳碳双键中的C4原子具有较强的相互作用, 双键中剩余的π电子离域到了环体系中.     

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)计算得到的速率常数与实验观测值进行了比较     

尿酸质子转移反应的理论研究I——分子内质子转移研究

使用量子化学中的Hartree-Fock方法和密度泛函理论中的B3LYP方法，分别在3－21G^*和6－31G（d)水平上，计算了尿酸分子从三羰基异构体向三羟基异构体的转化。结果表明，转化过程经历了单羟基和双羟基异构体2种中间物和3种过渡态时的分子内质子转移（IPT），转移中的H原邻近的N，O和C原子形成了具有四元环结构的过渡态。随着IPT的进行，N－H键逐渐被削弱和断裂，O－H键则逐渐生成。3个反应的活化能分别为190.3kJ/mol,181.4kJ/mol和249.9kJ/mol(B3LYP/6-31G(d))。较高的活化能表明在室温下，无催化剂的IPT难以进行。     

OH+C_2H_2N←C_2H_3+NO→OH+C_2H_2N反应机理的密度泛函理论研究

应用密度泛函理论研究了反应通道(a)C2H3+NO→CH3+NCO和(b)C2H3+NO→OH+C2H2N的反应机理.在B3LYP/6-31G(d)水平上优化了反应物、中间体、过滤态、产物的几何构型,通过频率分析确定了11个中间体和10个过渡态.所有的反应物、中间体、过渡态、产物都在CCSD/6-311++G(d,p)水平上进行了单点能较正.并讨论了反应的异构化过程.计算结果表明10是能量最低的中间体,比反应物的能量低308 479kJ/mol;过渡态1/3,2/5,3/4,4/8比反应物的能量高,其中3/4是能量最高的过渡态,比反应物的能量高91 894kJ/mol.通道(a)和(b)的理论放热值分别为111 059和96 619kJ/mol.     

单嘧磺隆除草剂水解机理的理论研究

用量子化学中密度泛函B3LYP方法在6-31G*水平上, 对单嘧磺隆的水解机理进行了理论研究. 优化了反应过程中反应物、中间体、过渡态及产物的几何构型, 并对各个过渡态和中间体进行了振动分析确认其真实性. 结果表明: 在碱性条件下, 水解反应主要经历一个五元环的过渡态TS1, 该反应速控步骤的能垒为25.7 kJ/mol; 在中性条件下, 水解反应涉及到一个三元环的过渡态和两个四元环的过渡态, 反应的速控步骤为M3'＋H₂O→TS2', 能垒为236.5 kJ/mol. 从能量上看, 碱性条件更有利于水解反应的进行, 与实验结果吻合. 同时, 我们还考虑了溶剂效应对水解反应的影响.     

OH+ C2H2N←C2H3 + NO→CH3 + NCO反应机理的密度泛函理论研究

应用密度泛函理论研究了反应通道(a)C2H3 NO→CH3 NCO和(b)C2H3 NO→OH C2H2N的反应机理．在B3LYP／6-31G(d)水平上优化了反应物、中间体、过滤态、产物的几何构型，通过频率分析确定了11个中间体和10个过渡态．所有的反应物、中间体、过渡态、产物都在CCSD／6-311 G(d，P)水平上进行了单点能较正．并讨论了反应的异构化过程．计算结果表明10是能量最低的中间体，比反应物的能量低308．479kJ／mol；过渡态1／3，2／5，3／4，4／8比反应物的能量高，其中3／4是能量最高的过渡态，比反应物的能量高91．894kJ／mol．通道(a)和(b)的理论放热值分别为111．059和96．619kJ／mol．     

4-三甲基硅-1-丁烯-3-炔合成反应机理的密度泛函法研究

借助于G03W程序的B3LYP方法,对Pd催化三甲基硅乙炔与溴乙烯生成4-三甲基硅-1-丁烯-3-炔历程中的反应物、中间体、过渡态以及产物的分子几何构型、能量等进行了优化,获得了反应势能面上各反应驻点的几何构型及其活化能,并确定了过渡态的结构,同时利用过渡态理论研究了其反应速率,探讨了反应机理.计算结果表明,Sonogashira反应历程中的最大吉布斯自由能出现在炔基加成阶段,ΔG298为150.72kJ/mol,即炔基加成为整个反应的决速步.     

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

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

巴比妥酸的质子迁移异构化反应及其芳香性的密度泛函研究

用密度泛函理论方法,在B3LYP/6-31G水平上研究了巴比妥酸的醇式和酮式结构的分子内质子迁移异构化反应,对反应势能面的研究发现,该化合物可能存在6种中间体和12种过渡态,对反应物、产物和过渡态进行了几何构型优化,在同一水平下计算了单点能量,并用频率振动模式和内稟坐标(IRC)确证了过渡态的存在,反应的活化能最小为107.50kJ·mol-1,最大为330.84kJ·mol-1。在GIAO-HF/6-31G//B3LYP/6-31G水平对反应物、产物,中间体和过渡态的非相关性化学位移(NICS)进行了计算,利用NICS值对反应过程中芳香性进行分析,发现如果分子中的C5与C6存在双键或者形成过渡态的四员环状结构,分子就有不同程度的芳香性,如果有C6上的羟基氢迁移到C5上,C5-C6以单键相结合,则分子就没有芳香性。     

2-磷酸甘油酸脱水反应机理的理论研究

采用密度泛函理论B3LYP方法研究了蛋白质和水环境下2-磷酸甘油酸脱水生成磷酸烯醇式丙酮酸的反应机理. 优化得到了反应物、过渡态及产物的几何构型并计算了反应势垒. 研究结果表明: 没有H₂O参与时, 反应需要通过四元环过渡态完成, 反应势垒高达287.7 kJ/mol, 常温下难以进行|有H₂O参与时, 反应可以通过六元环过渡态完成, 反应势垒大为降低|Mg^2＋的参与可使反应势垒进一步降低|蛋白质环境下两个Mg^2＋和一个H₂O的共同作用可使反应势垒降低至91.2 kJ/mol, 从而使反应在常温下容易进行.     

(4-溴甲基双环[4.4.1]-1,3,5,7,9-十一碳五烯基-3-)甲醇环氧化反应机理的理论研究

A reaction mechanism of epoxidation reaction of the 4-(bromomethyl)bicyclo4.4.1undeca-1,3,5,7,9-pentaen-3-ylmethanol has been studied by using the density functional theory(DFT) method at B3LYP level with 6 31G* basis set. The geometric structures of reactant, product and transition state have been optimized. The transition state is found by the QST2 method and characterized by the vibration frequency analysis. The intrinsic reaction coordinate(IRC) for this reaction is traced and confirms the reaction mechanism. The changes and the nature of related chemical bonds along the IRC path have been analyzed by the theory of electronic charge density. The result shows that the elimination reaction and ring closing reaction are synergistic, and the hydrogen in the elimination reaction derives from hydroxyl. In addition, the activation energy of the reaction is 139.2 kJ/mol.     

Theoretical Studies on the Reaction Mechanism of 1-Chloroethane with Hydroxyl Radical

The reaction mechanism of 1-chloroethane with hydroxyl radical has been inves- tigated by using density functional theory (DFT) B3LYP/6-31G (d, p) method. All bond dissociation enthalpies were computed at the same theoretical level. It was found that hydrogen abstraction pathway is the most favorable. There are two hydrogen abstraction pathways with activation barriers of 0.630 and 4.988 kJ/mol, respectively, while chlorine abstraction pathway was not found. It was observed that activation energies have a more reasonable correlation with the reaction enthalpy changes (△Hr) than with bond dissociation enthalpies (BDE).     

F原子与瞬态自由基CH_2SH反应的理论研究

用量子化学从头算和密度泛函理论（DFT）对F原子与自由基CH_2SH在势能面上 的反应进行了研究。在B3LYP/6-311G水平上计算出了各物种的优化构型、振动频率 和零点振动能(ZPVE);各物种的总能量由B3LYP/6-311 + G(2df, pd)//B3LYP/6- 311G计算,另外对反应物和产物还计算了其G3能量。结果表明：首先F通过与C或S 结合的两种途径与CH_2SH相配位,再通过H（4）原子转移形成甲基,然后甲基再旋 转,甲基中H（4）原子最终与F结合,反应产物为HF和CH_2S。反应为放热反应,分 别为ΔH_r = -370.7 kJ/mol (DFT)和-396.94 kJ/mol (G3)。此外依据计算出的反 应热,可得自由基·CH_2SH的生成热Δ_fH°_(298.15) = 146.44 kJ/mol (DFT), 而Δ_fH°_0 = 167.36 kJ/mol (G3)。它们与以前的实验和理论值是一致的。     

DFT法研究离子液中EMIM＋催化丁烯双键异构反应机理

利用密度泛函方法(DFT)分别在B3LYP/6-31G**和B3LYP/6-311＋＋G**的计算水平上优化了离子液体中1-乙基-3-甲基咪唑阳离子(EMIM＋)催化丁烯双键异构反应过程中的反应物、产物以及过渡态的几何构型,分析了反应过程中键参数的变化.通过振动分析对平衡态和过渡态进行了验证,并得到了零点能.通过计算内禀反应坐标(IRC),确认了对应于过渡态的反应物和产物. 计算结果表明,EMIM＋催化丁烯双键异构可以基元反应的方式一步完成,1-丁烯异构化为2-丁烯的活化能约为192 kJ•mol－1, 逆反应活化能约为208 kJ•mol－1, 可在室温或高于室温条件下进行.     

Transition state for the gas-phase reaction of uranium hexafluoride with water

Density functional theory and small-core, relativistic pseudopotentials were used to look for symmetric and asymmetric transition states of the gas-phase hydrolysis reaction of uranium hexafluoride, UF 6, with water. At the B3LYP/6-31G(d,p)/SDD level, an asymmetric transition state leading to the formation of a uranium hydroxyl fluoride, U(OH)F 5, and hydrogen fluoride was found with an energy barrier of +77.3 kJ/mol and an enthalpy of reaction of +63.0 kJ/mol (both including zero-point energy corrections). Addition of diffuse functions to all atoms except uranium led to only minor changes in the structures and relative energies of the reacting complex and transition state. However, a significant change in the structure of the product complex was found, significantly reducing the enthalpy of reaction to +31.9 kJ/mol. Similar structures and values were found for PBE0 and MP2 calculations with this larger basis set, supporting the B3LYP results. No symmetric transition state leading to the direct formation of uranium oxide tetrafluoride, UOF 4, was found, indicating that the reaction under ambient conditions likely includes several more steps than the mechanisms commonly mentioned. The transition state presented here appears to be the first published transition state for the important gas-phase reaction of UF 6 with water.     

Quantum mechanical study of the ring-closing reaction of the hexatriene radical cation

The ring-closing reaction of hexatriene radical cation 1(*)(+) to 1,3-cyclohexadiene radical cation 2(*)(+) was studied computationally at the B3LYP/6-31G* and QCISD(T)/6-311G*//QCISD/6-31G* levels of theory. Both, concerted and stepwise mechanisms were initially considered for this reaction. Upon evaluation at the B3LYP level of theory, three of the possible pathways-a concerted C(2)-symmetric via transition structure 3(*)(+) and stepwise C(1)-symmetric pathways involving three-membered ring intermediate 5(*)(+) and four-membered ring intermediate 6(*)(+)-were rejected due to high-energy stationary points along the reaction pathway. The two remaining pathways were found to be of competing energy. The first proceeds through the asymmetric, concerted transition structure 4(*)(+) with an activation barrier E(a) = 16.2 kcal/mol and an overall exothermicity of -23.8 kcal/mol. The second pathway, beginning from the cis,cis,trans rotamer of 1(*)(+), proceeds by a stepwise pathway to the cyclohexadiene product with an overall exothermicity of -18.6 kcal/mol. The activation energy for the rate-determining step in this process, the formation of the intermediate bicyclo[3.1.0]hex-2-ene via transition structure 9(*)(+), was found to be 20.4 kcal/mol. More rigorous calculations of a smaller subsection of the potential energy hypersurface at the QCISD(T)//QCISD level confirmed these findings and emphasized the importance of conformational control of the reactant.     

Biphenalenylidene: the forgotten bistricyclic aromatic ene. A theoretical study

The Lawesson reagent and P(2)S(5) mediated reductive coupling of phenalenone (6) gives LPAH peropyrene (5) in 47% and 54% yields. The mechanism of the reaction involves the formation of phenalenethione (10), Z- and E-1,1'-biphenalenylidene (3), and 9 as intermediates. The electrocyclization reaction of Z-3 to 9, followed by aromatization, gives 5. The results of an ab initio and DFT study of 3 and 2,2'-biphenalenylidene (12) are reported. E-3 and Z-3 have a diradical character with twist angles of 44.8 degrees and 57.8 degrees (at UB3LYP/6-311G**). Delta E(++)(Tot) = 10.2 kJ/mol and Delta G(++)(298) = 10.6 kJ/mol for E-3 <==> Z-3 diastereomerization. These unusually low energy barriers are due to the ground-state diradical destabilization and the aromatic stabilization of the transition state TS-3. Triplet Z-3 is higher in energy than singlets E-3 and Z-3 by 10.4 and 3.1 kJ/mol. In the concealed non-Kekulé 12, singlet 12 is more stable than the triplet by 1.3 kJ/mol. Singlet 12 is more stable than singlet E-3 by 2.0 kJ/mol, and orthogonal singlet TS-12 is lower in energy than singlet TS-3 by 6.0 kJ/mol. The energy barriers for the hexatriene-cyclohexadiene electrocyclization Z-3 --> 9 are Delta E(++)(Tot) = 94.8 and Delta G(++)(298) = 98.3 kJ/mol (at (U)B3LYP/6-31G). The reaction occurs thermally in a conrotatory mode.     

氧原子和甲基自由基反应机理的理论研究

用分子轨道从头算和密度泛函理论（DFT）中的B3LYP方法以及适中基组6-311+G(2df,2p)对氧原子与甲基CH3反应进行了系统的研究。计算给出了反应通道上各驻点物种的构型参数、振动频率和能量。结果表明： CH2OH比CH3O稳定,能量约低26.63 kJ/mol,且生成氢和甲醛为其最主要反应通道。     

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

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.     

过充电保护添加剂1,4-二甲氧基苯的反应机理

采用密度泛函理论(B3LYP/6-311+G(d,p))和MP2/6-311+G(d,p)方法,研究锂离子电池过充电保护添加剂1,4-二甲氧基苯(p-DMOB)的作用机理.计算结果表明,在过充时,p-DMOB优先于溶剂分子(乙基甲基碳酸酯、二甲基碳酸酯、碳酸乙酯)发生氧化反应.用B3LYP和MP2计算所得的p-DMOB理论氧化电位接近,分别为4.12和4.05V(vsLi/Li+).p-DMOB氧化时首先失去一个电子,生成p-DMOB+·正离子自由基,用B3LYP和MP2方法计算所得的相应能量变化分别为701.24和728.27kJ·mol-1.失去电子后苯环的共轭性受到破坏,随后p-DMOB+·苯环上的C―H键发生断裂,失去H+并形成p-DMOB·自由基.用B3LYP和MP2方法计算所得的相应能量变化分别为1349.78和1810.99kJ·mol-1.p-DMOB·自由基很不稳定,会在电极表面发生聚合反应形成聚合物膜,用B3LYP和MP2方法计算所得的相应能量变化分别为-553.37和-1331.20kJ·mol-1.