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键断裂反应可以忽略.     

CH3CHF自由基与HNCO反应机理的理论

采用密度泛函理论的B3LYP方法, 在6- 311++G(d, p)基组水平上研究了CH3CHF自由基与HNCO的微观反应机理, 优化了反应过程中的反应物、中间体、过渡态和产物, 在QCISD(T)/6- 311++G(d, p)水平上计算体系在反应通道各驻点的能量. 振动分析结果和IRC分析结果证实了中间体和过渡态的真实性, 计算所得的成键临界点电荷密度变化也确认了该反应过程, 并找到了七条反应通道. 其中生成氟代烷基酰亚胺稳定分子的通道活化能垒最低, 在该反应体系中是与氢迁移平行竞争较易发生的一条反应通道.     

BH2+与H2O反应机理的量子拓扑研究

采用密度泛函方法B3LYP/6-311+G(d, p)和耦合簇方法CCSD/6-311+G(d, p)研究了BH2+与H2O的气相离子-分子反应机理. 优化得到了反应途径中各驻点的几何构型, 并采用内禀反应坐标法进行追踪. 从量子拓扑学的角度, 讨论了在反应过程中各化学键的变化. 反应(I)经历了一个四元环过渡态, 找到了这个反应的能量过渡态和两个结构过渡态.     

CH3F与C2H3反应机理的理论研究

采用密度泛函B3LYP/6-311G(d,p)方法对CH3F与C2H3的反应体系进行了理论研究,获得了反应的势能面信息及可能的微观机理.在QCISD(T)/6-311++G(d,p)水平上精确计算了各反应物种的单点能.结果表明,除抽提氢反应外,标题反应还存在抽提氟(R1)、消氟化氢(R2)、消氢(R3)和自由基形成(R4)四类反应.在QCISD(T)/6-311++G(d,p)//B3LYP/6-311G(d,p)水平上,R1,R2,R3和R4反应的能垒分别是163.9,152.2,209.8和224.2kJ·mol-1,相应反应能为-56.6,-164.3,-2.7和-156.0kJ·mol-1,所有反应均放热,为热力学允许的反应.     

烷烃与氢过氧自由基氢提取反应类反应能垒与速率常数的精确计算

氢过氧自由基从烷烃分子中提取氢的反应是碳氢燃料中低温燃烧化学中非常重要的一类反应。本文用等键反应方法计算了这一类反应的动力学参数。所有反应物、过渡态、产物的几何结构均在HF/6-31+G(d)水平下优化得到。以反应中的过渡态反应中心的几何结构守恒为判据,该反应类可用等键反应处理。本文选取了乙烷和氢过氧自由基的氢提取反应为参考反应,其它反应作为目标反应,用等键反应方法对目标反应在HF/6-31+G(d)水平的近似能垒和反应速率常数进行了校正。为了验证方法的可靠性,选取C5以下的烷烃分子体系,对等键反应方法校正结果和高精度CCSD(T)/CBS直接计算结果进行了比较,最大绝对误差为5.58k J?mol~(-1),因此,采用等键反应方法只需用低水平HF从头算方法就可以再现高精度CCSD(T)/CBS计算结果,从而解决了该反应类中大分子体系的能垒的精确计算。本文的研究为碳氢化合物中低温燃烧模拟中重要的烷烃与氢过氧自由基氢提取反应提供了准确的动力学参数。     

喹唑啉酮互变异构的理论计算

采用密度泛函理论B3LYP/6-311+G(d,p)方法,计算并考察了喹唑啉酮进行结构互变的质子迁移过程的两种可能途径:(a)分子内质子迁移,(b)水助质子迁移.结果表明,途经b所需要的能垒小,氢键在降低反应能垒方面起重要作用.     

H2ClCS单分子分解反应机理的理论研究

用量子化学B3LYP/6 - 311+G(d,p)方法优化了H2ClCS单分子分解反应驻点物种的几何构型,并在相同水平上通过频率计算和内禀反应坐标(IRC)分析对过渡态结构及连接性进行了验证.用QCISD(T)/6-311++G(d,p)方法计算各物种的单点能,并对总能量进行了零点能校正.利用经典过渡态理论(TST)与...     

n‐Butane Monomolecular Cracking on Acidic Zeolites: A Density Functional Study

采用5T簇模型,利用密度泛函理论在B3LYP/6-311+G(3df,2p)//B3LYP/6-31G(d)水平下研究正丁烷在酸性分子筛上的单分子催化裂解反应。本文重点详细研究了正丁烷在分子筛表面不同C位的脱氢反应。在B3LYP/6-311+G(3df,2p)//B3LYP/6-31G(d)水平下计算所得第一和第二位C-C键裂解的活化能垒分别为 238、217 kJ/mol。而第一第二序位脱氢反应能垒分别为296、242 kJ/mol。正丁烷不同序位脱氢反应的活化能垒相差54 kJ/mol。从计算结果可以看出,正丁烷在分子筛上催化裂解脱氢反应优先发生在第二位C原子上。此外,本文还讨论了簇模型结构与酸性的关系,结果显示改变封端Si-H键的键长的方法可以用来模拟分子筛酸性变化。最后研究了分子筛酸性变化与正丁烷催化裂解反应能垒的关系。     

NO+HCCCO反应势能面的理论研究

在CCSD(T)/6-311G(d,p)//B3LYP/6-311G(d,p)+ZPVE水平下, 对反应NO＋HCCCO进行了研究. 建立了反应势能面, 揭示了该反应的反应机理, 通过O迁移、C—C键或N—O键断裂等多步反应, 得到4种产物, 其中, 最主要产物为P1(HCCO+NCO).     

烷基过氧化氢中氢提取反应类大分子体系的反应能垒与速率常数的精确计算

自由基与分子反应是一类具有负活化能的非基元反应,通常认为是通过反应复合物的两步过程,在大气化学和碳氢燃料燃烧机理中广泛存在,且在理论计算和实验上广泛研究.以碳氢燃料燃烧机理中重要反应类羟基自由基提取烷基过氧化氢α位氢的反应为研究对象,通过量化计算揭示其反应规律,计算得到其精确动力学参数.在所研究反应类中,定义第一步反应复合物的生成反应的标准摩尔吉布斯自由能变化等于零时所对应的温度为其转折温度T_c,并表明了当T >> T_c时可采用稳态近似法处理该类反应体系,得到总包反应速率常数.所有反应涉及的物种几何结构优化和频率分析均在BHandHLYP/6-311G（d,p）水平下得到,并在所研究反应类中选取了5个代表反应,通过CCSD（T）/CBS单点能计算,得到其最高转折温度为195.17 K,远远低于碳氢燃料燃烧模拟通常关注温度范围的最低温度650 K,表明用稳态近似法处理该类负活化能反应体系是合理的.计算还表明,该类反应的过渡态反应中心几何结构守恒,因此可将等键反应方法引入类反应,通过对低水平从头算得到的反应能垒进行校正,以得到高精度的结果.为了验证等键反应方法的可靠性,选取5个反应作为研究对象,将低水平BHandHLYP/6-311G（d,p）的校正结果和高水平CCSD（T）/CBS直接计算的结果进行比较,反应能垒最大绝对偏差由校正前的19.99 kJ·mol^－1降到校正后的1.47 kJ·mol^－1,表明用等键反应方法,只需在低水平从头算水平下就可以得到高水平的计算结果,从而可解决大分子体系精确动力学参数缺乏的问题.利用等键反应方法计算了20个反应的反应能垒,并结合过渡态理论计算得到了总包反应的速率常数,并揭示了该类反应只在低温段呈现负活化能关系.     

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

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

Comprehensive density functional theory study on serine and related ions in gas phase: conformations, gas phase basicities, and acidities

Density functional theory (DFT) calculations have been performed to investigate the gas-phase conformations of serine and its three related ions (serineH(+), serine(-), and serine(2-)). The full ensemble of possible conformations, 324 conformations for serine, 108 for serineH(+), 162 for serine(-) and 54 for serine(2-), were first surveyed at B3LYP/6-31G level, and then the obtained unique conformations were further refined at B3LYP/6-311+G level. From full optimizations, 74 unique conformations for seine, 14 for serineH(+), 11 for serine(-), and 4 for serine(2-) were located, and their relative energies were also determined at B3LYP/6-311+G level. Atoms in molecules (AIM) analysis was carried out to establish rigorous definition of hydrogen bonds. Six types of intramolecular H-bonds in conformers of serine, six types in serineH(+), three types in serine(-), and two types in serine(2-) were identified within the framework of AIM theory and their relative strengths were determined based on topological properties at bond critical points (BCPs) of H-bonds. The intramolecular H-bonds were demonstrated to play an important role in deciding the relative stability of conformations of amino acids and the related ions. The enthalpies and Gibbs free energies of protonation and deprotonation reactions of serine and its related ions were calculated at B3LYP/6-311+G//B3LYP/6-31G, and B3LYP/6-311+G//B3LYP/6-311+G level. The calculated results are both in excellent agreement with the experimental data. We demonstrate in this study that B3LYP is an efficient and accurate method to predict the thermochemical and structural parameters of amino acids and the related ions.     

Isomerization of B_6,B_6~- and B_6~+ clusters

The interconversions between isomers with the same spin multiplicity of neutral B6 and charged B6-and B6+ clusters have been investigated at the B3LYP/6-311+G level of theory,including determination of the minimum energy pathways with transition states connecting the corresponding reactants and products.In dynamic calculations,26 isomers were optimized,including 11 novel isomers.In order to further refine the energies,single-point B3LYP/6-311+G(3df) calculations were carried out on the corresponding B3LYP/6...     

Theoretical studies on the intramolecular hydrogen bond and tautomerism of 8-mercaptoquinoline in the gaseous phase and in solution using modern DFT methods

A theoretical quantum chemical study of the intramolecular hydrogen bonding interactions in 8-mercaptoquinoline has been carried out. Special attention has been paid to the rotation of S-H bond and intramolecular proton-transfer reactions. Therewith, the B3LYP/6-311++G(d,p), B3LYP/6-31+G(2d,2p), MPW1K/6-311++G(d,p), MPW1K/6-31+G(2d,2p), BH&HLYP/6-311++G(d,p), and G96LYP/6-311++G(d,p) methods have been used. By means of the Onsager and PCM reaction field methods, the effects of solvent on hydrogen-bond energies, conformational equilibria, rotational barriers, and tautomerism in aqueous solution have been studied. These simulations were done at the MPW1K/6-311++G(d,p) and B3LYP/6-311++G(d,p) levels. Natural-bond orbital analysis has been performed to study the intramolecular hydrogen bond (IHB) in the gaseous phase and in aqueous medium. The stability of forms under consideration in solution does not coincide with that in the gaseous phase, underlining a great importance of the electrostatic influence of solvent. Double-proton transfer in the prototropic tautomerization of 8-mercaptoquinoline, one water molecule complex in the gaseous phase and in solution, has been systematically studied. The double-proton transfer occurs concertedly and synchronously. The water-assisted tautomerization is kinetically less, but thermodynamically more favorable, compared to that of the single-proton transfer. As in the case with single-proton transfer, for water-assisted reaction, the tautomerization energies and barrier heights decrease with the increase in dielectric constant, which implies faster and more complete tautomerization of 8-mercaptoquinoline in a polar solvent.     

Theoretical Studies on Energetic Property and Stability of Pyrazine and Pyridine Derivatives

Density functional calculations at the B3LYP level with 6‐311G** and aug‐cc‐pVDZ basis sets were performed to predict the heats of formation (HOFs) for two pyrazine derivatives and eight pyridine derivatives. In the isodesmic reactions designed for the computation of heats of formation (HOFs), pyrazine and pyridine were chosen as reference compounds. The N‐oxidations for the ring nitrogen of pyrazine and pyridine derivatives decrease the HOF values when N‐oxide oxygen is neighboring with amino groups, but increase when it neighbors with nitro groups. Thermal stability was evaluated via bond dissociation energies (BDE) at the UB3LYP/6‐311G** level. As a whole, the homolysis of C–NO₂ bonds is the main step for bond dissociation of the title compounds. The BDE values of title compounds are influenced by intramolecular hydrogen bonds. The hydrogen bond effects associated with the length of the H···O bonds were analyzed by the electron density at the critical points and natural bond orbital.     

吡咯-HCN体系在气相及溶液中相互作用的理论研究

用量子化学B3LYP方法在6-311G(d, p)水平上优化了吡咯-HCN氢键复合物，通过振动频率分析确定了两个吡咯-HCN体系稳定构型.为了得到更加精确的氢键作用能，采用相关一致基组aug-cc-pVDZ以及Boys 和Bernardi的CP(counterpoise)校正方法消除基组重叠误差后得到C－H…π和N－H…N型复合物的氢键相互作用能.为了确定B3LYP方法计算的相互作用能的可靠性，在MP2/aug-cc-pVDZ水平计算了复合物的氢键相互作用能，结果分别为－25.10和－19.30 kJ·mol－1.采用自然键轨道(NBO)分析考察了吡咯与HCN分子间轨道相互作用.以自洽场理论(SCRF)中的Onsager模型研究了不同极性溶剂对吡咯-氰化氢体系N－H…N型氢键几何构型，频率位移，电荷分布以及相对能量的影响.研究发现，当溶液的介电常数在1.5～30.0范围时，溶液作用十分显著，而当介电常数超过30.0以后，溶液作用已经达到了极限.     

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.     

Theoretical study of the reactions of CF3OCHF2 with the hydroxyl radical and the chlorine atom.

Theoretical investigations are carried out on the reaction mechanism of the reactions of CF3OCHF2 (HFOC-125) with the OH radials and Cl atoms, as well as the heats of formation of CF3OCHF2 and CF3OCF2. The electronic structure information on the potential energy surface for each reaction is obtained at the B3LYP/6-311G(d,p) level, and energetic information is further refined by G3(MP2) theory. The direct dynamics calculation of the hydrogen abstraction reactions are also performed at the G3(MP2)//B3LYP/ 6-311G(d,p) level. The classical energy profile is corrected by interpolated single-point energies (ISPE) approach, incorporating the small-curvature tunnelling effect (SCT) calculated by the variational transition-state theory (VTST). The rate constants are in good agreement with the experimental data and are found to be k1 = 4.95 x 10(-30)T(5.40)exp(-347/T) and k2 = 1.86 x 10(-23)T(3.43)exp(-1579/T) cm3 molecule(-1)s(-1) over the temperature range 200-2000 K. The rate constants at 298 K for these two reactions are 3.38 x 10(-16) and 2.80 x 10(-17) cm3 molecule(-1)s(-1), respectively. Using group-balanced isodesmic reactions as working chemical reactions, the standard enthalpies of formation for CF3OCHF2 and CF3OCF2 are -312.3 +/- 1.0 and -257.3 +/- 1.0 kcalmol(-1), respectively, evaluated by G3(MP2) theory based on the B3LYP/6-311G(d,p) geometries. The theoretical studies provide rate constants of the title reactions and the enthalpies of the species, which are important parameters in determining the atmospheric lifetime and the feasible pathways for the loss of HFOC-125.     

Thermochemistry of the higher chlorine oxides ClOx (x=3, 4) and Cl2Ox (x=3–7)†

Heats of formation for ClO₃, ClO₄, Cl₂O₃, Cl₂O₄, Cl₂O₅, Cl₂O₆ and Cl₂O₇ molecules are determined at the B3LYP, B3PW91, mPW1PW91 and B1LYP levels of the density functional theory employing a series of extended basis sets, and using Gaussian-3 model chemistries. Modified Gaussian-3 calculations, which employ accurate B3LYP/6-311+G(3d2f) molecular geometries and vibrational frequencies, were also performed. Heats of formation were calculated from both total atomization energies and isodesmic reaction schemes. The latter method in conjunction with Gaussian-3 models leads to the most reliable results. The best values at 298 K for ClO₃, ClO₄, Cl₂O₃ and Cl₂O₄ as derived from an average of G3//B3LYP and G3//B3LYP/6-311+G(3d2f) calculations are 43.1, 54.8, 31.7 and 37.4 kcal mol⁻¹. From calculations carried out at the G3(MP2)//B3LYP and G3(MP2)//B3LYP/6-311+G(3d2f) levels, heats of formation for Cl₂O₅, Cl₂O₆ and Cl₂O₇ are predicted to be 53.2, 52.2 and 61.5 kcal mol⁻¹. All best values are reproduced within 1 kcal mol⁻¹ by using mPW1PW91/6-311+G(3d2f) isodesmic energies. Enthalpy changes for relevant Cl–O bond fission reactions are reported. Comparisons with previous thermodynamics data are made.