Dyotropic rearrangement of alpha-lactone to beta-lactone: a computational study of small-ring halolactonisation

Transition structures have been optimised using the B3LYP/6-31+G* density functional level method, in vacuum and in implicit (PCM) and explicit (DFT/MM) aqueous solvation, for the degenerate rearrangement of the alpha-lactone derived by the formal addition of Cl(+) to acrylate anion and for the dyotropic rearrangement of this to the beta-lactone. Despite being lower in energy than the alpha-lactone, there is no direct pathway to the beta-lactone from the acrylate chloronium zwitterion, which is the transition structure for the degenerate rearrangement. This may be rationalised by consideration of the unfavorable angle of attack by the carboxylate nucleophile on the beta-position; attack on the alpha-position involves a less unfavorable angle. Formation of the beta-lactone may occur by means of a dyotropic rearrangement of the alpha-lactone. This involves a high energy barrier for the acrylate derived alpha-lactone, but dyotropic rearrangement of the beta,beta-dimethyl substituted alpha-lactone to the corresponding beta-lactone involves a much lower barrier, estimated at about 46 kJ mol(-1) in water, and is predicted to be a facile process.     

DFT法研究3-羟基丙烯醛的双键旋转异构反应机理

利用密度泛函理论(DFT)分别在B3LYP/6-31G**和B3LYP/6-311＋＋G**的计算水平上优化了基态3-羟基丙烯醛分子在双键旋转异构反应过程中的平衡态以及过渡态的几何构型，分析了反应过程中键参数的变化，计算了该反应的内禀反应坐标(IRC)，发现在重排反应途径上存在一个四元环骨架的中间体.通过振动分析对平衡态和过渡态进行了确认，并得到了零点能.计算结果表明，基态3-羟基丙烯醛分子的双键旋转异构反应经过两步完成，第一步反应位垒稍高，第二步反应位垒较低，存在着发生重排反应的可能性.     

DFT Studies on the Isomerization of Butene Double Bond Catalyzed by 1-Butyl-3-methyl-imidazolium in Ionic Liquid

1 INTRODUCTION Butene and its isomers are important petroleum raw materials. Isomerization reaction of butene plays a key role in the course of C4 alkylation and its reaction mechanism has captured the attention of chemists all along[1, 2]. As a green so…     

苯并氧化呋咱稳定性和异构化的DFT和ab initio研究

运用B3LYP/6-31G(d)密度泛函理论（DFT）方法对苯并氧化呋咱、邻二亚硝基苯及其间的异构化反应进行了计算研究。结果表明,苯并氧化呋咱的分子总能量比邻二亚硝基苯的低;由苯并氧化呋咱异构为邻二亚硝基苯的正向反应活化能（Ea+=51.0kJ/mol）,与文献实测值（58.6kJ/mol）较接近,而其逆向反应活化能（Ea-=4.6kJ/mol）很小,从而揭示了苯并氧化呋咱比邻二亚硝基苯更稳定·此外,进行了HF/3-21G、HF/6-31G(d)和MP2/6-31G(d)//6-31G(d)水平下相应的计算,发现B3LYP-DFT的结果较abinitio为优。谐振动频率的B3LYP/6-31G(d)计算还支持了邻二亚硝基苯为苯并氧化呋咱“自-自”互变重排反应的中间体。     

Dissociative and associative mechanisms of cope rearrangements of fluorinated 1,5-hexadienes and 2,2'-bis-methylenecyclopentanes

The effects of fluorine substitution on the Cope rearrangements of 1,5-hexadiene and 2,2'-bis-methylenecyclopentane have been examined computationally using (U)B3LYP/6-31+G(d,p) and CASPT2/6-31G(d) methodology. The calculations indicate that fluorine substituents at the hexadiene termini generally stabilize the transition states relative to the ground states of the chair conformations and destabilize pathways that occur via boat conformations, in accord with the experimental observations of Dolbier. With meso-2,2'-bis(difluoromethylene)cyclopentane, this destabilization is sufficient to favor a relatively dissociative concerted transition state resembling weakly interacting allyl radicals over a normal aromatic concerted boat transition state. This preference is due partly to an increase in the activation enthalpy for the concerted rearrangement coupled with the more favorable entropy for dissociation. In octafluoro- and decafluorohexadienes, the situation is reversed, and reaction through a cyclohexane-1,4-diyl is favored. Even in the octafluoro system with no radical stabilizing substituents at C(2) and C(5), the preference of fluorine for sp(3) centers causes reaction via the cyclohexane-1,4-diyl. In establishing methodology for this study, the conformations of 1,2-difluoroethane and 1,1,2,2-tetrafluoroethane were also examined thoroughly by the B3LYP method using three basis sets.     

2,2-二氟-1-溴-1-锂乙烯的结构及重排反应的理论研究

用量子化学中的密度泛函DFT方法,在B3LYP／6-31G(d,p)水平上研究了2,2-二氟-1-溴-1-锂乙烯F2C—CLiBr的结构．结果表明,F2C=CLiBr有2种平衡结构,其中只有1种是稳定的．对稳定的平衡结构,找到了其可能的重排反应过渡态,根据计算得到的重排反应势垒,解释了氟原子只从溴原子对位发生迁移的原因．     

Theoretical investigation on the GaCl3-catalyzed ring-closing metathesis reaction of N-2,3-butadienyl-2-propynyl-1-amine: three-membered ring versus four-membered ring mechanism

The gallium chloride (GaCl(3))-catalyzed ring-closing metathesis reaction mechanism of N-2,3-butadienyl-2-propynyl-1-amine has been studied at the Becke three-parameter hybrid functional combined with Lee-Yang-Parr correlation functional (B3LYP)/6-31G(d), B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p)//B3LYP/ 6-31G(d) and the second-order M?ller-Plesset perturbation (MP2)/6-311++G(d,p)//B3LYP/6-31+G(d,p) levels. It was found that the final metathesis product can be yielded via a three-membered or four-membered ring mechanism. The three-membered ring pathway is favorable due to its low energy barrier at the rate determining step. The whole reaction is stepwise and strongly exothermic.     

A computational study of conformational interconversions in 1,4-dithiacyclohexane (1,4-dithiane)

Ab initio molecular orbital theory with the 6-31G(d), 6-31G(d,p), 6-31+G(d), 6-31+G(d,p), 6-31+G(2d,p), 6-311G(d), 6-311G(d,p), and 6-311+G(2d,p) basis sets and density functional theory (BLYP, B3LYP, B3P86, B3PW91) have been used to locate transition states involved in the conformational interconversions of 1,4-dithiacyclohexane (1,4-dithiane) and to calculate the geometry optimized structures, relative energies, enthalpies, entropies, and free energies of the chair and twist conformers. In the chair and 1,4-twist conformers the C-Hax and C-Heq bond lengths are equal at each carbon, which suggest an absence of stereoelectronic hyperconjugative interactions involving carbon-hydrogen bonds. The 1,4-boat transition state structure was 9.53 to 10.5 kcal/mol higher in energy than the chair conformer and 4.75 to 5.82 kcal/mol higher in energy than the 1,4-twist conformer. Intrinsic reaction coordinate (IRC) calculations showed that the 1,4-boat transition state structure was the energy maximum in the interconversion of the enantiomers of the 1,4-twist conformer. The energy difference between the chair conformer and the 1,4-twist conformer was 4.85 kcal/mol and the chair-1,4-twist free energy difference (deltaG degrees (c-t)) was 4.93 kcal/mol at 298.15 K. Intrinsic reaction coordinate (IRC) calculations connected the transition state between the chair conformer and the 1,4-twist conformer. This transition state is 11.7 kcal/mol higher in energy than the chair conformer. The effects of basis sets on the 1,4-dithiane calculations and the relative energies of saturated and unsaturated six-membered dithianes and dioxanes are also discussed.     

7-Ethyl-2,3,5,6,8-pentahydroxy-1,4-naphthoquinone (echinochrome A): a DFT study of the mechanism of the antioxidant action. 4. Migration of protons and the Na<Superscript>+</Superscript> cation in echinochrome a monosodium salts

According to B3LYP/6-31G(d,p) and B3LYP/6-31G(d) calculations, intramolecular migrations of H⁺ and Na⁺ cations in the molecule of 7-ethyl-2,3,5,6,8-pentahydroxy-1,4-naphthoquinone (echinochrome A) monosodium salt in the gas phase and in clusters containing five water molecules can proceed asynchronously through the chain of O-H...O hydrogen bonds by the jump-over mechanism.     

DFT study on mechanism of the classical Biginelli reaction

The condensation of benzaldehyde, urea, and ethyl acetoacetate according to the procedure described by Biginelli was investigated at the B3LYP/6-31G（d）, B3LYP/6-31 ＋G（d,p）, and B3LYP/6-311＋G（3df,2p）//B3LYP/6-31＋G（d,p） levels to explore the reaction mechanism. According to the mechanism proposed by Kappe, structures of five intermediates were optimized and four transition states were found. The calculation results proved that the mechanism proposed by Kappe is right.     

Reaction of ketenyl radical with acetylene: a promising route for cyclopropenyl radical

The reaction of the ketenyl radical (HCCO) with acetylene (C(2)H(2)) is very relevant to the oxygen-acetylene flames and fuel-rich combustion process for nitrogen-containing compounds. Unfortunately, except for several rate constant measurements, the mechanism is completely unknown for this reaction. In this paper, detailed theoretical investigations are performed for the HCCO + C(2)H(2) reaction at the G3B3 level using the B3LYP/6-31G(d), B3LYP/6-311++G(d,p), and QCISD/6-31G(d) geometries. The exclusive fragmentation channel is the formation of the cyclopropenyl radical (c-C(3)H(3)) and carbon monoxide (CO) via the chainlike OCCHCHCH and three-membered ring OC-cCHCHCH intermediates. Thus, the mass spectroscopic peak of C(3)H(3)(+) in a previous experiment can be explained. The calculated overall reaction barrier is 4.4, 4.4, and 5.3 kcal/mol at the G3B3//B3LYP/6-31G(d), G3B3//B3LYP/6-311++G(d,p), and G3B3//QCISD/6-31G(d) levels, respectively. The title reaction may provide an effective route for generating the long-sought cyclopropenyl radical in the laboratory, which has been the long-standing subject of numerous theoretical studies as the simplest cyclic conjugate radical, and its bulky derivatives were already known. Future experimental investigations for the HCCO + C(2)H(2) reaction are greatly desired to test the predicted fragmentation channel. The implication of the present study in combustion and interstellar processes is discussed.     

Solvent and stereoelectronic effects on the solvolysis rates of oxaspirocyclopropanated 1-norbornyl triflates and related bridgehead derivatives

The study of the stereochemical outcome of the solvolysis of oxaspirocyclopropanated 1-norbornyl triflates is highly interesting since these reactions do not lead to the usual retention or fragmentation products but only synthetically interesting rearranged products are enantiospecifically formed. There is no correlation between the experimental solvolysis rates (ln k) and the B3LYP/6-31G(d)-computed ionization energies (Delta E) of the corresponding bridgehead hydrocarbons in gas phase. However, this work demonstrates the existence of a fair linear correlation between the experimental reaction rates and the PCM//B3LYP/6-31G(d)-computed free ionization energies in solution (Delta G). This theoretically relevant result reveals that the reason for the lack of linearity in gas phase is not the rearrangement of the intermediate carbocations but unspecific solvent effects on the solvolysis rates, accounted for by the PCM model.     

Formation mechanism of furfuryl sulfides from o-furfuryl dithiocarbonates: density functional theory study for aromatic [3,3]-sigmatropic rearrangement

Density functional theory (DFT) calculations at the B3LYP/6-31G(d) and B3LYP/6-31G+(d) levels demonstrated that O-furfuryl S-alkyl dithiocarbonate (1) undergoes aromatic [3,3]-sigmatropic rearrangement to the energetically unfavorable S-(2-methylene-2,3-dihydrofuran-3-yl) S-alkyl dithiocarbonate (2'), which then rearranges to furfuryl alkyl sulfide (3) with COS extrusion to regain the aromaticity lost in the first step.     

过氧烷基自由基分子内氢迁移反应类速率常数的计算

过氧烷基自由基分子内氢迁移是低温燃烧反应中的一类重要基元反应. 本文用等键反应方法计算了该类反应的动力学参数. 所有反应物、过渡态、产物的几何结构均在B3LYP/6-311+G(d,p)水平下优化得到. 本文提出了用过渡态反应中心几何结构守恒作为反应类判据, 并将该分子内氢迁移反应分为四类, 包括(1,3)、(1,4)、(1,5)、(1,n) (n=6, 7, 8)氢迁移类. 分别将这4 类反应类中最小反应体系作为类反应的主反应, 并分别在B3LYP/6-311+G(d,p)低水平和CBS-QB3 高水平下得到其近似能垒和精确能垒. 其余氢迁移反应作为目标反应, 在B3LYP/6-311+G(d,p)低水下计算得到其近似能垒, 再采用等键反应方法校正得到目标反应的精确反应势垒和精确速率常数. 研究表明, 采用等键反应方法只需在低水平用从头算计算就可以得到大分子反应体系的高精度能垒和速率常数值, 且本文按等键反应本质的分类方法更能揭示反应类的本质, 并对反应类的定义给出了客观标准. 本文的研究为碳氢化合物低温燃烧模拟中重要的过氧烷基分子内氢迁移反应提供了准确的动力学参数.     

An ab initio study on thermal rearrangement reactions of 1-silylprop-2-en-1-ol H3SiCH(OH)CH=CH2

The thermal rearrangement reactions of 1-silylprop-2-en-1-ol H3SiCH(OH)CH=CH2 were studied by ab initio calculations at the G2(MP2) and G3 levels. The reaction mechanisms were revealed through ab initio molecular orbital theory. On the basis of the MP2(full)/6-31G(d) optimized geometries, harmonic vibrational frequencies of various stationary points were calculated. The reaction paths were investigated and confirmed by intrinsic reaction coordinate (IRC) calculations. The results show that the thermal rearrangements of H3SiCH(OH)CH=CH2 happen in two ways. One is via the Brook rearrangement reactions (reaction A), and the silyl group migrates from carbon atom to oxygen atom passing through a double three-membered ring transition state, forming allyloxysilane CH2=CHCH2OSiH3. In the other, the reactant undergoes a dyotropic rearrangement; the hydroxyl group migrates from carbon atom to silicon atom coupled with a simultaneous migration of a hydrogen atom from silicon atom to carbon atom, forming allylsilanol CH2=CHCH2SiH2OH (reaction B). The barriers for reactions A and B were computed to be 343.5 and 203.7 kJ/mol, respectively, at the G3 level. The changes of the thermodynamic functions, entropy (DeltaS), entropy (DeltaS(doubledagger)) for the transition state, enthalpy (DeltaH), and free energy (DeltaG) were calculated by using the MP2(full)/6-31G(d) optimized geometries, and harmonic vibrational frequencies of reactants, transition states, and products with statistical mechanical methods, and equilibrium constant K(T) and reaction rate constant k(T) in canonical variational transition-state theory (CVT) with centrifugal-dominant small-curvature tunneling approximation (SCT) were calculated over a temperature range 400-1300 K. The conventional transition-state theory (TST) rate constants were also calculated for the purposes of comparison. The influences of the vinyl group attached to the center carbon of the alpha-silyl alcohols on reactions were discussed.     

过充电保护添加剂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.     

Substituent effects on the rearrangements of cyclohexyl to cyclopentyl radicals involving avermectin-related radicals

The rearrangement of a substituted cyclohexyl radical to a cyclopentylmethyl radical on the skeleton of avermectin B1 has been investigated using density functional (UB3LYP/6-31G(d)) and G3MP2B3 computational methods. The rearrangement is preferred when highly radical stabilizing groups are present at the 2- and 3-positions of the cyclohexyl radical. A substituent on the 3-position of the cyclohexyl radical enables ring-cleavage of the cyclohexyl radical, while a radical stabilizing substituent on the 2-position of the cyclohexyl radical stabilizes the final cyclopentylmethyl radical, enabling the overall rearrangement and reversing the normal thermodynamic preference for the hexenyl radical ring closure.     

CH_3SH+CN微观反应机理的理论研究

采用BMC-CCSD//B3LYP/6-311G(d,p)方法对CH3SH+CN反应机理进行了详细的理论研究.反应中涉及的各稳定点的构型、振动频率和零点能在B3LYP/6-311G(d,p)水平下计算得到,计算结果表明,该反应存在两种反应机理,5条可能的反应通道.SN2机理由于能垒太高,与直接氢抽提机理相比可以忽略.该反应的最可行通道为CN中的C原子进攻SH中的H原子经由一个前期和一个后期分子络合物生成产物CH3S和HCN.计算得到的反应焓变与已有实验值非常吻合.     

Computational study of the deamination reaction of cytosine with H2O and OH-

The mechanism for the deamination reaction of cytosine with H(2)O and OH(-) to produce uracil was investigated using ab initio calculations. Optimized geometries of reactants, transition states, intermediates, and products were determined at RHF/6-31G(d), MP2/6-31G(d), and B3LYP/6-31G(d) levels and for anions at the B3LYP/6-31+G(d) level. Single-point energies were also determined at B3LYP/6-31+G(d), MP2/GTMP2Large, and G3MP2 levels of theory. Thermodynamic properties (DeltaE, DeltaH, and DeltaG), activation energies, enthalpies, and free energies of activation were calculated for each reaction pathway that was investigated. Intrinsic reaction coordinate analysis was performed to characterize the transition states on the potential energy surface. Two pathways for deamination with H(2)O were found, a five-step mechanism (pathway A) and a two-step mechanism (pathway B). The activation energy for the rate-determining steps, the formation of the tetrahedral intermediate for pathway A and the formation of the uracil tautomer for pathway B, are 221.3 and 260.3 kJ/mol, respectively, at the G3MP2 level of theory. The deamination reaction by either pathway is therefore unlikely because of the high barriers that are involved. Two pathways for deamination with OH(-) were also found, and both of them are five-step mechanisms. Pathways C and D produce an initial tetrahedral intermediate by adding H(2)O to deprotonated cytosine which then undergoes three conformational changes. The final intermediate dissociates to product via a 1-3 proton shift. Deamination with OH(-), through pathway C, resulted in the lowest activation energy, 148.0 kJ/mol, at the G3MP2 level of theory.