3-羟基-3-甲基-2-丁酮和苯甲酰甲酸甲酯热分解反应的 理论研究

用从头算的方法在6-31G水平上研究了3-羟基-3-甲基-2-丁酮(1)和苯甲酰甲酸甲酯(2)热分解反应的机理。结果是:前一反应是经历五元环过渡态到达氢键中间体,它接着直接分解成乙醛的异构体和丙酮,最后乙醛的异构体异构化成乙醛;后一反应经历六元环过渡态形成中间体1(INT1),中间体1(INT1)直接分解成中间体2(INT2)和甲醛,中间体2(INT2)经过第二个过渡态分解成苯甲醛的异构体和一氧化碳,最后苯甲醛异构体异构化成苯甲醛。其中氢迁过程是反应的速控步骤。在MP~2/6-31G//HF/6-31G+ZPE水平上,对应于这两个反应速控步骤的活化位垒分别是251.42kJ/moL和247.94kJ/mol。采用传统过渡态理论计算了两反应的热反应速率常数,理论的计算结果与实验值吻合较好。     

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

用密度泛函理论方法,在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以单键相结合,则分子就没有芳香性。     

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

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

H CH_2CO反应机理的G2计算

分别在UQCISD/6-311G(d,p)和G2理论计算水平上,对CH2CO和H反应可能存在的四条反应通道进行了研究,详细分析了每个通道的反应机理;通过振动分析的虚频数和内禀反应坐标(IRC)计算,确认了反应涉及的每一个过渡态.通过反应位能剖面的比较,发现经过一个中间体生成CH3 CO的一条途径是主反应通道,该通道是个放热反应,总焓变为-146.07kJ·mol-1,速控步骤的位垒为55.09kJ·mol-1.理论计算结果较好地解释了实验观察到的主要产物和副产物并存的现象.     

二氯乙烯锗烯与甲硫醛环加成的反应机理

利用MP2/6-311+G*方法计算了单线态二氯乙烯锗烯与甲硫醛的各种反应机理. 计算结果表明两者之间的环加成反应具有很好的选择性, 优势反应路径分为三步: 首先两种反应物经过无能垒的放热反应形成中间体INT, 然后INT经历过渡态TS3异构化为P31, 最后P31继续与甲硫醛反应形成二环杂环化合物P33. 其中第一步反应放热103.4 kJ·mol-1, 后两步反应能垒分别为0.7 和32.3 kJ·mol-1.     

β-羟基丙醛基态与激发态氢迁移反应的从头算研究

本文用从头算RHF和UHF方法在3-21G基组上研究了β-羟基丙醛基态和激发态分解为甲醛和乙烯醇的反应机理。优化得到了各反应途径的过渡态和中间体, 其结果为: 基态β-羟基丙醛经过一个六元环过渡态和一个氢键中间体形成产物, 反应属于氢迁移和断键的协同过程; 激发三态β-羟基丙醛的分解途径首先经过一个氢迁移六元环过渡态形成双自由基中间体, 然后该中间体的分解包括两条相互竞争的途径, 它们各自经过一个断碳碳键的过渡态和一个氢键激-基态配合物中间体而形成两类产物, 一类为甲醛的基态和乙烯醇的激发态, 另一类为甲醛的激发态和乙烯醇的基态。激发态反应的两条通道均属于先氢迁移后断键分解的分步过程, 且反应的第二步为速控步骤。计算结果表明, 激发态反应活化位垒都比基态的低。     

反式共轭类碳烯CH2=CH-CH=C:LiF氢迁移反应的DFT研究

应用密度泛函理论DFT方法,在B3LYP/6-31G*水平上研究了反式共轭类碳烯CH2=CH-CH=C∶LiF的氢迁移反应.计算优化了反应过程中的所有反应物、中间体、过渡态和产物的几何构型,通过频率振动分析确定中间体和过渡态.结果表明,在β位上的H原子迁移过程中,经历一个带有三元环结构的中间体和两个带有三元环结构的过渡态.第一步反应势垒较大.     

硅杂环丙烯异构化反应的理论研究

用量子化学理论方法研究了硅杂环丙烯单重态的异构化反应。结果表明:该异构化反应的过渡态为三元环结构,该反应为氢迁移反应;反应的势垒高度为276.67kJ·mol^-^1[MP2/6-31G(d)]。通过内禀反应坐标(IRC)计算,获得了沿反应途径的势能剖面。     

单重态CCl2与O3反应机理的理论研究

用量子化学从头算方法,研究了单重态CCl_2与O_3反应的机理.在HF/6-31G(d)水平上用梯度解析技术全参数优化上述反应的反应物、中间体、过渡态和产物构型,MP2/6-31G(d)//HF/6-31G(d)方法计算能量.给出了有关化合物的结构数据.结果表明:CCl_2与O_3首先生成富能中间体CCl_2O_3,然后中间体裂解生成CCl_2O和O_2.该反应为强放热反应,放出的热量为516.88kJ·mol~(-1)[MP2/6-31G(d)//HF/6-31G(d)].通过内禀反应坐标(IRC)计算,获得了沿反应途径的势能剖面.     

HNCO+OH——NH~2+CO~2反应理论研究

用从头算UHF/6-31G基组研究了异氰酸和羟基生成氨基和二氧化碳即HNCO+OH--NH~2+CO~2的反应机理.优化得到了反应途径上的过渡态和中间体,并通过振动分析对过渡态和中间体进行了确认.在UMP4/6-31G水平上计算了它们的能量,同时对零点能进行了较正.计算结果表明:此反应是多步反应,先后通过3个过渡态(TS1,TS2,TS3),2个内旋转位垒(TSI,TSII),4个中间体(IM1,IM2,IM3,IM4),其中,IM3--TS2这一步为整个反应的决速步骤,速控步的活化能为202.388kJ/mol.与异氰酸和羟基作用的另一反应通道(即HNCO+OH--H~2O+NCO)的活化能(69.038kJ/mol)比较,可看出所研究反应通道为次要反应通道,这与实验结果是一致的。     

Molecular decomposition of polyolefins: Kinetic parameters and geometry of the transition state

The experimental data on the molecular decomposition of olefins ( 1 -pentens) of various structures to two olefins in a gas phase were analyzed by means of the parabola intersection method. The enthalpies and kinetic parameters characterizing such decomposition were calculated for eighteen reactions. Decomposition of olefins representing two-centered concerted re action was found to be characterized by a very high classic potential barrier of thermoneutral reaction (197.4 kJ/mol). The kinetic parameters (activation energy and velocity constant) of fifteen reverse reactions of formation of olefins from two alkenes were calculated using the parabola-intersection method. The factors affect ing the activation energy of the reactions of olefin decomposition and formation are discussed. Quantum chemical calculations of the transition state energy and geometry for six reactions of olefin decomposition were performed.     

二水醋酸镉与邻氨基苯甲酸的分步固相反应的热化学研究

用溶解量热法,以一定比例的盐酸(0.2mol·L     

The catalytic power of pyruvate decarboxylase. A stochastic model for the molecular evolution of enzymes.

Pyruvate decarboxylase (PDC) catalyzes the decarboxylation of pyruvate anion by a factor of around 10(12), compared with the non-enzymic decarboxylation by thiamine, under standard state conditions of 1 mM pyruvate and thiamine diphosphate (TDP), pH 6.2. Free-energy diagrams constructed on the basis of earlier measurements for the enzymic and non-enzymic reactions give some information on catalysis by PDC. PDC stabilizes the reactant state preceding TDP addition to pyruvate by 76 kJ mol-1 and the transition state for the addition by 83 kJ mol-1. PDC stabilizes the reactant state preceding decarboxylation (presumably alpha-lactyl-TDP) by 27 kJ mol-1 and the transition state for decarboxylation by 68 kJ mol-1. In addition, the free-energy diagrams reveal a leveling of reactant-state free energies in the enzymic reaction compared with the non-enzymic reaction, in that the former are nearly equal to each other. The enzyme-bound transition-state energies are similarly leveled. The energetic leveling of reactant states has been noted by Albery, Knowles and their coworkers in many enzymic reactions and termed 'matched internal thermodynamics.' They showed that the result would arise naturally (and inevitably) in the 'evolution to perfection' of enzymes, when the evolutionary process was treated by a deterministic model. The critical assumption of this model was the validity of a Marcus-type or Br?nsted-type linear free-energy relationship between rate and equilibrium constants for reactions occurring wholly within enzyme complexes. Here a completely stochastic simulation of molecular evolution, with no deterministic assumptions, is shown to reproduce both 'matched internal thermodynamics' and the 'matched internal kinetics' or leveling of transition-state energies noted here. The Albery-Knowles result is thus more general than might have been supposed.     

大气中一元酸催化亚硫酸分解反应的从头算及动力学研究

采用从头算的理论计算方法,在CCSD（T）/aug-cc-pVDZ//MP2/aug-cc-pVDZ水平下,对甲酸、乙酸、丙酸和硝酸催化亚硫酸分解成二氧化硫和水的微观反应机理进行了理论研究.结果表明,这4种一元酸均可使亚硫酸分解反应能垒显著降低,降低幅度由大到小的顺序为丙酸 > 乙酸 > 甲酸 > 硝酸.其中,丙酸甚至使反应能垒从裸反应时的99.84 kJ/mol降至27.24 kJ/mol.在此基础上,计算了200~320 K温度范围内4种一元酸催化亚硫酸分解反应的速率常数,并结合它们在大气中的实际浓度计算了有效速率常数.结果表明,在实际大气环境中,乙酸对亚硫酸分解反应的催化效果最为显著.当乙酸存在时,亚硫酸在室温下的大气寿命仅为0.02 s.     

Mechanistic investigation of the oxygen-atom-transfer reactivity of dioxo-molybdenum(VI) complexes

The oxygen-atom-transfer (OAT) reactivity of [LiPrMoO2(OPh)] (1, LiPr=hydrotris(3-isopropylpyrazol-1-yl)borate) with the tertiary phosphines PEt3 and PPh2Me in acetonitrile was investigated. The first step, [LiPrMoO2(OPh)]+PR3-->[LiPrMoO(OPh)(OPR3)], follows a second-order rate law with an associative transition state (PEt3, DeltaH not equal=48.4 (+/-1.9) kJ mol-1, DeltaS not equal=-149.2 (+/-6.4) J mol-1 K-1, DeltaG not equal=92.9 kJ mol-1; PPh2Me, DeltaH not equal=73.4 (+/-3.7) kJ mol-1, DeltaS not equal=-71.9 (+/-2.3) J mol-1 K-1, DeltaG not equal=94.8 kJ mol-1). With PMe3 as a model substrate, the geometry and the free energy of the transition state (TS) for the formation of the phosphine oxide-coordinated intermediate were calculated. The latter, 95 kJ mol-1, is in good agreement with the experimental values. An unexpectedly large O-P-C angle calculated for the TS suggests that there is significant O-nucleophilic attack on the P--C sigma* in addition to the expected nucleophilic attack of the P on the Mo==O pi*. The second step of the reaction, that is, the exchange of the coordinated phosphine oxide with acetonitrile, [LiPrMoO(OPh)(OPR3)]+MeCN-->[LiPrMoO(OPh)(MeCN)]+OPR3, follows a first-order rate law in MeCN. A dissociative interchange (Id) mechanism, with activation parameters of DeltaH not equal=93.5 (+/-0.9) kJ mol-1, DeltaS not equal=18.2 (+/-3.3) J mol-1 K-1, DeltaG not equal=88.1 kJ mol-1 and DeltaH not equal=97.9 (+/-3.4) kJ mol-1, DeltaS not equal=47.3 (+/-11.8) J mol-1 K-1, DeltaG not equal=83.8 kJ mol-1, for [LiPrMoO(OPh)(OPEt3)] (2 a) and [LiPrMoO(OPh)(OPPh2Me)] (2 b), respectively, is consistent with the experimental data. Although gas-phase calculations indicate that the Mo--OPMe3 bond is stronger than the Mo--NCMe bond, solvation provides the driving force for the release of the phosphine oxide and formation of [LiPrMoO(OPh)(MeCN)] (3).     

A quantum chemistry study of the Cl atom reaction with formaldehyde

The elementary vapor-phase reaction between Cl atoms and HCHO has been studied by ab initio methods. Calculations at the MP2, MP3, MP4(SDTQ), CCSD, CCSD(T), and MRD-CI levels of theory show that the reaction is characterized by a low electronic barrier; excluding the effects of spin-orbit splitting in Cl, our best estimate at the MRD-CI/aug-cc-pVTZ//RHF-RCCSD(T)/aug-cc-pVTZ level of theory predicts a Born-Oppenheimer barrier height of 0.7 kJ mol-1. The energies of the lowest six electronic states as resulting from MRD-CI calculations are presented at discrete points along the reaction path, and two avoided crossings are found in the transition state region. The spin-orbit splitting in Cl is also calculated along the reaction path; it is not negligible in the transition state region and is found to increase the barrier by only 1.4 kJ mol-1 at the RCCSD(T)/aug-cc-pVTZ transition state geometry. The minimum energy path of the reaction connects an energetically weakly stabilized adduct on the flat potential surface on the reactant side and an energetically strongly stabilized postreaction adduct. The reaction rate coefficient and the kinetic isotope effects were calculated using improved canonical variational theory with small curvature tunneling (ICVT/SCT), and the results were compared to experimental data. The experimental reaction rate coefficient is reproduced within its uncertainty limits by variational transition state theory with interpolated single-point energy corrections (ISPE) at the MP4(SDTQ) level of theory and by conventional transition state theory with interpolated optimized energies (IOE) at the MRD-CI//RCCSD(T) level of theory and interpolated optimized geometries at the RCCSD(T) level of theory on an MP2/aug-cc-pVTZ potential energy surface when employing scaled vibrational frequencies.     

Cl+F2→ClF+F和F+ClF′→ClF+F′反应机理的密度泛函理论研究

用密度泛函理论(DFT)B3LYP方法, 在6-311G基组下,计算研究了反应Cl+F2→ClF+F和对称反应F+ClF′→ClF+F′的机理。求得前者的过渡态为三角形,活化能为15.57 kJ*mol-1;后者的过渡态为线形和三角形,活化能分别为11.52和196.25 kJ*mol-1。结果均经过振动分析和IRC计算验证。     

Calculations on the competition between association and reaction for C3H(+) + H2

A potential energy surface has been calculated for the competing associative and reactive ion-molecule processes involving the reactants C3H(+) + H2. Our ab initio results show that the linear ion C3H+ and H2 can directly access the deep potential well of the propargyl ion H2CCCH+, which is calculated to lie 390 kJ mol-1 below the zero-point energy of the reactants. Isomerization between the propargyl ion and the lower energy, cyclic C3H3+ ion, calculated to lie 501 kJ mol-1 below the zero-point energy of reactants, can subsequently occur via two pathways. One of these pathways involves a transition state lying 22 kJ mol-1 below the energy of the reactants while the other, which occurs at much lower energies, involves two transition states and an intermediate. The dissociation of c-C3H3+ into c-C3H2(+) + H is calculated to occur directly, without any intermediate potential energy maximum, but the energy of the products lies 7.3 kJ mol-1 above the energy of the reactants. Using the minimum energy potential pathway and properties of the stationary point structures determined via ab initio methods, we have calculated both the association rate coefficient to produce C3H3+ as a function of density and the branching ratio between the propargyl and cyclic structures of the ion. Our results are in good agreement with some experimental results and in conflict with others. Specifically, we agree with the 1:1 branching ratio measured for the propargyl and cyclic isomers of C3H3+ at 80 and 300 K and we agree with the rate coefficient for radiative association measured at 80 K. We cannot reproduce reported measurements that the reactive channel (C3H2(+) + H) is the dominant channel at 80 K and at low gas densities, or that the association channel at high densities saturates at an effective rate coefficient well below the Langevin value -2x10(-11) cm3 s-1 at 300K and 1x10(-10) cm3 s-1 at 80K.     

HgH4 and HgH6: further candidates for high-valent mercury compounds

Mercury tetrahydride (D4h) is calculated to have similar bond lengths and vibrational frequencies as the already known HgH2 and to lie energetically 200 kJ mol-1 above HgH2 + H2, in a local well, about 40 kJ mol-1 below a transition state.