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

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

异硫氰酸与甲亚胺环加成反应机理的理论研究

对异硫氰酸与甲亚胺形成四元或六元环产物的环加成反应进行了理论研究,结果表明,一分子异硫氰酸与一分子甲亚胺形成四元环(1:1)产物P1的反应(1)为经过一个两性离子中间体的分步反应,其中第二步为速控步骤,其活化热垒为107.86kJ/mol.此外,反应(1)的中间体还可与另一甲亚胺或异硫氰酸分子继续反应形成两个不同的六元环(1:2或2:1)产物P2或P3;这两个反应均为协同反应,其活化势垒分别为15.88和21.82kJ/mol.这些结果与当异硫氰酸酯与亚胺发生环加成反应时,只有类似于P2和P3的两种六元环产物生成的实验事实一致。     

Double Bond Migration Mechanism in Allyl Systems Involving the Hydroxide Ion. 1. Gas-Phase and Born–Onsager Models

The reaction profile of the 1,3-prototropic rearrangement of propene involving the hydroxide ion was studied by the RHF/6-31+G*, MP2/6-31+G*, and B3LYP/6-31+G* ab initio methods within the framework of the gas-phase and Born–Onsager models (the latter including solvent effects). Propene isomerization in the presence of the hydroxide ion in the gas phase may occur with participation of a base proton with the intermediate formation of a water complex of the allyl ion. The transition state energy of this transformation is lower than the total energy of the starting hydroxide ion and propene and much lower than the sum of the energies of the isolated propenide ion and water molecule. An activation barrier arises when the solvent effect is included in calculation within the framework of the Born–Onsager model; the intermediate complex is much less stable than the complex considered in the gas-phase model. As in the latter, the mechanism of multiple bond migration is energetically preferable to the mechanism involving proton transfer to the reaction medium.     

Modern valence-bond description of chemical reaction mechanisms: the 1,3-dipolar addition of diazomethane to ethene

The electronic mechanism for the gas-phase concerted 1,3-dipolar cycloaddition of diazomethane (CH2N2) to ethene (C2H4) is described through spin-coupled (SC) calculations at a sequence of geometries along the intrinsic reaction coordinate obtained at the MP2/6-31G(d) level of theory. It is shown that the bonding rearrangements occurring during the course of this reaction follow a heterolytic pattern, characterized by the movement of three well-identifiable orbital pairs, which are initially responsible for the pi bond in ethene and the C-N pi bond and one of the N-N pi bonds in diazomethane and are retained throughout the entire reaction path from reactants to product. Taken together with our previous SC study of the electronic mechanism of the 1,3-dipolar cycloaddition of fulminic acid (HCNO) to ethyne (C2H2) (Theor. Chim. Acc. 1998, 100, 222), the results of the present work suggest strongly that most gas-phase concerted 1,3-dipolar cycloaddition reactions can be expected to follow a heterolytic mechanism of this type, which does not involve an aromatic transition state. The more conventional aspects of the gas-phase concerted 1,3-dipolar cycloaddition of diazomethane to ethene, including optimized transition structure geometry, electronic activation energy, activation barrier corrected for zero-point energies, standard enthalpy, entropy and Gibbs free energy of activation, have been calculated at the HF/6-31G(d), B3LYP/6-31G(d), MP2/6-31G(d), MP2/6-31G(d,p), QCISD/6-31G(d) and CCD/6-31G(d) levels of theory. We also report the CCD/6-311++G(2d, 2p)//CCD/6-31G(d), MP4(SDTQ)/6-311++G(2d,2p)//CCD/6-31G(d) and CCSD(T)/6-311++G(2d, 2p)//CCD/6-31G(d) electronic activation energies.     

环丙基硅烯的理论研究: 环丙基硅烯C3H5SiH的重排反应及其机理

本文用限制的Hartree-Fock解析梯度方法在3-21G和6-31G^*水平上对环丙基硅烯的重排反应及其机理进行了从头算研究。以6-31G^*优化构型作了二级微扰计算, 并计算了各构型的频率。在此基础上得到了重排反应的热焓△H, 自由能△G和平衡常数K, 用Eyring过渡态理论计算了反应的速度常数k(T), 应用Woodward-Hoffmann规则讨论了环丙基硅烯重排反应过程中端基的旋转机理。结果表明, 环丙基硅烯经过113.4kJ/mol的势垒扩环重排为硅杂环丁烯为自发反应; 而其1,2-氢迁移重排反应热垒为190.0kJ/mol, 是非自发反应, 难于进行,不能与扩环重排相竞争。另外, 扩环重排反应可分为多步过程, 每步的端基旋转均可用Woodward-Hoffmann规则说明。     

Theoretical studies on the gas-phase pyrolysis of 2-phenoxycarboxylic acids: an ONIOM approach

Various ONIOM combinations-ONIOM(HF/6-31G*: PM3), ONIOM(B3LYP/6-31G*: PM3), ONIOM(MP2/6-31G*: PM3), and ONIOM(MP2/6-31G*: HF/3-21G)--were applied to investigate thermal decomposition mechanisms of four 2-phenoxycarboxylic acids (2-phenoxyacetic acid, 2-phenoxypropionic acid, 2-phenoxybutyric acid, and 2-phenoxyisobutyric acid) in the gas phase. All the transition states and intermediates of the reaction paths were optimized. The reaction pathway of four reactants yielding the phenol, CO, and the corresponding carbonyl compound was characterized on the potential energy surface and found to proceed stepwise. The first step corresponds to the elimination of phenol and the formation of alpha-lactone intermediate through a five-membered ring transition state, and the second step is the cycloreversion process of alpha-lactone intermediate to form CO and the corresponding carbonyl compound. The reaction pathway of latter three compounds to produce the carboxylic acid and phenol via a four-membered cyclic transition structure was also examined theoretically. Comparison with experiment indicates that the activation parameters for the fist reaction channel are accurately predicted at the ONIOM(MP2/6-31G*: HF/3-21G) level of theory.     

Theoretical study of geometrical effect on the deoxygenation of epoxide by singlet carbenes

The mechanism of the abstraction of oxygen from epoxide by carbenes has been investigated using the B3LYP/6-31G* and MP2/6-311++G(d,p)//B3LYP/6-31G* levels of theory. According to our model calculations, the reactivity of carbene decreases in the order: four-membered ring carbene > six-membered ring carbene > three-membered ring carbene > five-membered ring carbene vinylidene. Moreover, the present work shows that the singlet–triplet splitting of a carbene can be used as a guide for the prediction of its reactivity.     

Theoretical studies on the conformations of selenamides

Ab initio HF/6-31+G*, MP2/6-31+G*, B3LYP/6-31+G* level calculations have been performed on HSe-NH₂ to estimate the Se-N rotational barriers and N-inversion barriers. Two conformers have been found withsyn andanti arrangement of the NH₂ hydrogens with respect to Se-H bond. The N inversion barriers in selenamide are 1.65, 2.47, 1.93 kcal/mol and the Se-N rotational barriers are 6.58, 6.56 and 6.12 kcal/mol respectively at HF/6-31+G*, MP2/6-31+G* and B3LYP/6-31+G* levels respectively. The n_N →Σ *_Se-H negative hyperconjugation is found to be responsible for the higher rotational barriers.     

Adiabatic Potential Energy Surface of the 1,2,3-Trifluorobenzene Cation-Radical

Ab initio study of the adiabatic potential energy surface (PES) for 1,2,3-trifluorobenzene cation-radical was carried out by the ROHF, UHF, MCSCF, MP2, and DFT (B3LYP) methods with the 6-31G* basis set. The PES in question is a pseudorotation surface at all calculation levels. The pseudorotation barrier height does not exceed 3.2 kcal/mol, suggesting a possibility for its manifestations in experimental EPR spectra.     

Synthesis,Characterization and Theoretical Studies on a Thiocarbamide Derivative Containing Schiff Base Group

A thiocarbamide derivative containing Schiff base groups,1,5-bis[4-(dimethylamino)benzylidene]thiocarbonohydrazide,has been synthesized and characterized by elemental analysis,IR,1H NMR,UV and X-ray single-crystal diffraction.Density function theory(DFT) calculations at the B3LYP/6-31G* and PBE0/6-31G* levels for optimized geometries and electronic transition spectra have been performed.Comparative studies show that both B3LYP/6-31G* and PBE0/6-31G* methods can well reproduce the molecular structure,and the latter is more reliable than the former to simulate electronic spectra.NPA calculational results at the B3LYP/6-31G* level indicate the title compound to be a potential multidentate ligand to the metallic ions.Based on the vibrational analysis,thermodynamic properties at different temperatures have been obtained.     

Conformational study of the structure of 12-crown-4-alkali metal cation complexes

A conformational search was performed for the 12-crown-4 (12c4)-alkali metal cation complexes using two different methods, one of them is the CONFLEX method, whereby eight conformations were predicted. Computations were performed for the eight predicted conformations at the HF/6-31+G*, MP2/6-31+G*//HF/6-31+G*, B3LYP/6-31+G*, MP2/6-31+G*//B3LYP/6-31+G*, and MP2/6-31+G* levels. The calculated energies predict a C4 conformation for the 12c4-Na+, -K+, -Rb+, and -Cs+ complexes and a C(s) conformation for the 12c4-Li+ complex to be the lowest energy conformations. For most of the conformations considered, the relative energies, with respect to the C4 conformation, at the MP2/6-31+G*//B3LYP/6-31+G* are overestimated, compared to those at the MP2/6-31+G* level, the highest level of theory considerd in this report, by 0.2 kcal/mol. Larger relative energy differences are attributed to larger differences between the B3LYP and MP2 optimized geomtries. Binding enthalpies (BEs) were calculated at the above-mentioned levels for the eight conformations. The agreement between the calculated and experimental BEs is discussed.     

Double Bond Migration Mechanism in Allyl Systems Involving the Hydroxide Ion. 2. Polarizable Continuum Model (PCM)

Isomerization processes of a double bond site in propene and methylthiopropene molecules with the hydroxide ion were studied in the framework of the RHF/6-31+G*, MP2/6-31+G*, and B3LYP/6-31+G* (density functional) ab initio methods. The solvent effect was taken into account using PCM in its IEFPCM and SCIPCM versions. It is shown that to construct a reaction profile for propene rearrangement, it suffices to perform geometry optimization of stationary points within the Born–Onsager model with further refinement of the energy using IEFPCM. The reaction profiles obtained display that the multiple bond migration mechanism involving the hydroxide ion proton is energetically preferable to the two-stage mechanism forming a solvated carbanion for the propene molecule and for the methylthiopropene molecule that forms a much more stable carbanion.     

Ab initio and DFT calculations on the initial step in thiamin catalysis

The diphosphate ester (ThDP) of thiamin (vitamin B₁) is an important cofactor of enzymes within the carbohydrate metabolism. The initial reaction step shared by all ThDP-dependent enzymes is the deprotonation of the C2–H of the thiazolium ring. The replacement of the 4′-amino group by a hydroxyl one in the pyrimidine ring leads to the oxy-ThDP analogue which is known as an antagonist in thiamin catalysis.

Ab initio and DFT calculations on the MP2/6-31G* and B3LYP/6-31G* level were performed to study the proton relay function in thiamin and oxythiamin systems. Both MP2 and B3LYP calculations show significant differences of the reaction coordinate of the ylide formation in the systems. Tautomers, protonated and deprotonated species of both systems show different trends regarding their stability. The influence of correlation effects on the results is discussed by comparison with the HF-SCF/6-31G* data. Frequency calculations on the B3LYP/6-31G* level were performed to characterize the minima and transition state structures, respectively.     

Structure and conformational stability of CH2CHCH2X (X=F, Cl and Br) molecules—post Hartree–Fock and density functional theory methods

The molecular structure and conformational stability of CH₂CHCH₂X (X=F, Cl and Br) molecules were studied using ab initio and density functional theory (DFT) methods. The molecular geometries of 3-fluoropropene were optimized employing BLYP and B3LYP levels of theory of DFT method implementing 6-311+G(d,p) basis set. The MP2/6-31G^*, BLYP and B3LYP levels of theory of ab initio and DFT methods were used to optimize the 3-chloropropene and 3-bromopropene molecules. The structural and physical parameters of the molecules are discussed with the available experimental values. The rotational potential energy surface of the above molecules were obtained at MP2/6-31G^* and B3LYP/6-311+G(d,p) levels of theory. The Fourier decomposition of the rotational potentials were analyzed. The HF/6-31G^* and MP2/6-31G^* levels of theory have predicted the cis conformer as the minimum energy structure for 3-fluoropropene, which is in agreement with the experimental values, whereas the BLYP/6-311+G(d,p) and B3LYP/6-311+G(d,p) levels of theory reverses the order of conformation. The ΔE values calculated for 3-chloropropene at MP2/6-31G^*, BLYP/6-311+G(d,p) and B3LYP/6-311+G(d,p) levels of theory show that the gauche form is more stable than the cis form, which is in agreement with the experimental value. The same levels of theory have also predicted that the gauche form is stable than cis for 3-bromopropene molecule. The maximum hardness principle has been able to predict the stable conformer of 3-fluoropropene at HF/6-31G^* level of theory, but the same level of theory reverses the conformational stability of 3-chloropropene and 3-bromopropene molecules and MP2/6-31G^* level of theory predicted the stable conformer correctly.     

Comparative theoretical study of transition structures,barrier heights,and reaction energies for the intramolecular tautomerization in acetaldehyde/vinyl alcohol and acetaldimine/vinylamine systems

The transition structures associated with the possible intramolecular tautomerization for acetaldehyde/vinyl alcohol and acetaldimine/vinylamine systems as models of keto/enol and imine/enamine interconversion processes, respectively, were characterized. The relative stabilities of the tautomers and the associated barrier heights were calculated. Ab initio analytical gradients and second derivatives at the HF level of theory and 3-21G, 6-31G, 6-31G**, 6-31++G**, and 6-311++G** basis-set, DFT (BP86/6-311++G** and BLYP/6-311++G**), and semiempirical (AM1 and PM3) procedures were used to identify the stationary points. Correlation effects were estimated using the perturbational approach at MP2/6-31G**, MP2/6-311++G**, and MP2/6-311++G (3df,2p) levels. The geometry, electronic structure, harmonic vibrational frequencies, and transition vector associated with the transition structures as well as the relative stabilities of different isomers and barrier heights were analyzed. The dependence of these properties upon theoretical methods is analyzed and discussed. The transition structures are four-membered rings and the corresponding transition vectors are associated to collective fluctuations. The 1,3 intramolecular hydrogen migration is much more advanced than are the hybridization changes on donor and acceptor centers at the transition structure. The corresponding barrier heights can be related to the change of bond orders and acid/base properties of these centers. A comparison of the results obtained with different methods renders that the nature of the transition structure seems to be a rather robust entity. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 66 : 9–24, 1998     

Ab initio studies on mechanisms of cycloaddition reaction between dichloro-vinylidene and ethylene

Mechanisms of cycloaddition reaction between singlet dichloro-vinylidene (R1) and ethylene (R2) have been investigated with MP2 and B3LYP /6-31G* methods, including geometry optimization, vibrational analysis, and energy for the involved stationary points on the potential energy surface. CCSD(T)/6-31G* single point calculations are also applied to the geometries from both methods. CCSD(T) relative energies for the stationary points predicted by MP2 and B3LYP are in a very good agreement. Three reaction pathways are located. Along the first one, one intermediate (INT1) is firstly generated, which then rearranges into a three-membered ring compound (P1) via a small barrier of 5.4 kJ/mol; the other two paths share the other intermediate (INT2), which isomerizes into two four-membered ring compounds (P2 and P3) via a chlorine and a hydrogen transfer, respectively. The barriers of the latter two paths are significantly higher by approximately 25 kJ/mol than that of the former (27.2 and 28.8 vs 5.4 kJ/mol), the major reaction is therefore the formation of P1.     

D-苯丙氨酸与Cu＋(1S0, 3d10)气相反应理论研究

用量子化学密度泛函(DFT)方法研究D-苯丙氨酸与一价基态金属阳离子Cu^＋在气相中反应的机理. 在B3LYP/6-31G*水平上, 优化了反应包含的4个反应通道的反应物、中间体、过渡态和产物的几何构型, 并采用B3LYP/DZVP, B3LYP/[6-311＋G**(C,H,O)＋Lanl2dz(Cu)], B3LYP/6-311＋G**, MP2/6-311＋G** 等方法对各驻点进行了单点能计算. 通过对计算结果的分析, 获得了其单重态反应势能面的一般轮廓、各驻点几何构型优化参数, 明确了其反应机理.     

Reaction of hypochlorous acid with imidazole: formation of 2-chloro- and 2-oxoimidazoles

Reaction of hypochlorous acid (HOCl) with imidazole (Im) taken as a model for the 5-membered ring of guanine, leading to the products 2-chloro- and 2-oxo-imidazoles was investigated at the B3LYP/6-31+G* and B3LYP/AUG-cc-pVDZ levels of density functional theory. For all cases, single point energy calculations were performed at the MP2/AUG-cc-pVDZ level of theory using the geometries optimised at the B3LYP/AUG-cc-pVDZ level. Intrinsic reaction coordinate calculations were performed to ensure genuineness of all the calculated transition states. Effect of aqueous media was investigated by solvating all the species involved in the reactions using the polarizable continuum model. It is found that 2-chloroimidazole (2-ClIm) can be formed following three different reaction schemes while 2-oxoimidazole (2-oxoIm) can be formed following two different reaction schemes. The calculated barrier energies show that formation of 2-oxoIm would be less favored than that of 2-ClIm, which explains the experimental observations on relative yields of 8-chlorodeoxyguanosine and 8-oxodeoxyguanosine.     

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

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.     

Theoretical Study on the Tautomerization Reaction of Formamide to Formamidic Acid Encapsulated in a Nano-sized Molecular Container

The hybrid calculations with ONIOM(B3LYP/6-31G*:AM1) method were carried out on the tautomerization reaction of formamide to formamidic acid in the microcontainer-encapsulated state. The free-state tautomerization process was also investigated with B3LYP/6-31G**//B3LYP/6-31G* method for the purpose of comparison. Bare tautomerization, H2O-assisted(single-H2O or multiple-H2O) and self-assisted mechanisms were all taken into consideration for the encapsulated state. The results show that only bare tautomerization and single-H2O catalysis mechanisms are possible to the encapsulated for-mamide tautomerization owing to the container's size confinement effect. Geometrical changes in the complexed container and guest molecules are discussed to deeply understand the complex's structural properties. The bare tautomerization barrier in the encapsulated state increases by 23.826 kJ/mol, ac-counting for 12% of the corresponding total energy barrier in the free state, and the increased values for the single-H2O catalysis process are 12.958 kJ/mol, accounting for 16% of the corresponding total energy barrier, respectively. This finding suggests that the encapsulation can make the tautomerization process slightly difficult.