Br2+Cl2=2BrCl反应机理的理论和实验研究

用(U)MP2方法, 取6-311G*基组, 研究了反应Br2+Cl2=2BrCl的机理, 求得四中心和三中心的过渡态, 通过比较反应通道的活化能的大小, 得到如下结论: 双分子基元反应的最小活化能小于Cl2和Br2的离解能, 在没有光引发的条件下, 标题反应将以分子与分子作用形式完成; 若有光引发, Br2或Cl2先解离成原子, 再经过Br原子与Cl2反应或Cl原子与Br2反应, 能较快完成标题反应. 分别测定了光照和避光两种条件下的反应体系在412 nm处吸光度的变化, 证实了理论研究的结果.     

Ab Initio Study of the Dehydrogcnation of Methylenimine

The ab initio calculations were performed by the intrinsic reaction coordinated(IRC) method for the reaction paths of the dehydrogenations of methylenimine. We determined the geometries and energies of the transition states and obtained the activation energies, activation entropies and the statistical A factors on the RHF/4 -31G singlet potential energy surface. The mode selective study reveals the intrinsic reaction coordi-nates(URC) of 1 ,1 and 1 ,2 dehydrogenations of methylenimine to be related with the H2C= group symmetrical stretch mode and the =NH group stretch mode, respectively. The vibration coupling between the normal coordinates occurs in either of these two reaction paths.     

A hybrid elastic band string algorithm for studies of enzymatic reactions

A common challenge in theoretical biophysics is the identification of a minimum energy path (MEP) for the rearrangement of a group of atoms from one stable configuration to another. The structure with maximum energy along the MEP approximates the transition state for the process and the energy profile itself permits estimation of the transition rates. In this work we describe a computationally efficient algorithm for the identification of minimum energy paths in complicated biosystems. The algorithm is a hybrid of the nudged elastic band (NEB) and string methods. It has been implemented in the pDynamo simulation program and tested by examining elementary steps in the reaction mechanisms of three enzymes: citrate synthase, RasGAP, and lactate dehydrogenase. Good agreement is found for the energies and geometries of the species along the reaction profiles calculated using the new algorithm and previous versions of the NEB and string techniques, and also those obtained by the common method of adiabatic exploration of the potential energy surface as a function of predefined reaction coordinates. Precisely refined structures of the saddle points along the paths may be subsequently obtained with the climbing image variant of the NEB algorithm. Directions in which the utility of the methods that we have implemented can be further improved are discussed.     

Transition structures, energetics, and secondary kinetic isotope effects for cope rearrangements of cis-1,2-divinylcyclobutane and cis-1,2-divinylcyclopropane: a DFT study

The minimum energy reaction paths and secondary kinetic isotope effects (KIE) for the Cope rearrangements of cis-1,2-divinylcyclobutane and cis-1,2-divinylcyclopropane obtained by (U)B3LYP calculations are reported. Both reactions proceed through endo-boatlike reaction paths, and have aromatic transition states. The predicted activation energies are in agreement with the experimental data. The reaction paths of the rearrangements are intervened by enantiomerization saddle points of the products (and the reactant in the case of divinylcyclobutane). The calculated KIEs are similar in the two systems, and consistent with the geometries of the transition structures. There is computational evidence that the isotope effect associated with the conversion of a pure sp(2) C-H bond into a pure sp(3) one might be the same in all molecules. The predicted KIEs agree with experiment for divinylcyclopropane, but not for divinylcyclobutane.     

A theoretical study of the alkylation reaction of toluene with methanol catalyzed by acidic mordenite.

A theoretical study of the alkylation reaction of toluene with methanol catalyzed by the acidic Mordenite (Si/Al = 23) is reported. Cluster DFT as well as periodical structure DFT calculations have been performed. Full reaction energy diagrams of the elementary reaction steps that lead to the formation of the three xylene isomers are given. The use of periodical structure calculations allows one to account for zeolite framework electrostatic contributions and steric constraints that take place in zeolitic catalysts. Especially the steric constraint energy contribution has a significant effect on the energies and bond formation paths. The activation energy barrier of p-xylene formation is found to be approximately 20 kJ/mol lower than the corresponding values for the formation of its isomers. Computed host-guest binding energies according to the DFT method need a correction due to the absence of the dispersive interaction with the zeolite wall. Apparent activation energies obtained with this correction are in good agreement with experimental data.     

A DFT study of SiH(4) activation by Cp(2)LnH

A theoretical study of SiH(4) activation by Cp(2)LnH complexes for the entire series of lanthanides has been carried out at the DFT-B3PW91 level of theory. The reaction paths corresponding to H/H exchange and silylation, formation of Cp(2)Ln(SiH(3)), have been computed. They both occur via a single-step sigma-bond metathesis mechanism. For the athermal H/H exchange reaction, the calculated activation barrier averages 1.8 kcal.mol(-)(1) relative to the precursor adduct Cp(2)LnH(eta(2)-SiH(4)) for all lanthanide elements. The silylation path is slightly exogenic (DeltaE approximately -6.5 kcal.mol(-1)) with an activation barrier averaging 5.2 kcal.mol(-1) relative to the precursor adduct where SiH(4) is bonded by two Si-H bonds. Both pathways are therefore thermally accessible. The H/H exchange path is calculated to be kinetically more favorable whereas the silylation reaction is thermodynamically preferred. The reactivity of this familly of lanthanide complexes with SiH(4) contrasts strongly with that obtained previously with CH(4). The considerably lower activation barrier for silylation relative to methylation is attributed to the ability of Si to become hypervalent.     

Haptic quantum chemistry

We present an implementation designed to physically experience quantum mechanical forces between reactants in chemical reactions. This allows one to screen the profile of potential energy surfaces for the study of reaction mechanisms. For this, we have developed a interface between the user and a virtual laboratory by means of a force‐feedback haptic device. Potential energy surfaces of chemical reactions can be explored efficiently by rendering in the haptic device the gradients calculated with first‐principles methods. The underlying potential energy surface is accurately fitted on the fly by the interpolating moving least‐squares (IMLS) scheme to a grid of quantum chemical electronic energies (and geometric gradients). In addition, we introduce a new IMLS‐based method to locate minimum‐energy paths between two points on a potential energy surface. © 2009 Wiley Periodicals, Inc. J Comput Chem 2009     

CuCl2催化（2-甲基辛烷-2，3-二烯-4-基）磷酸乙酯氯代环化反应机理理论研究

采用密度泛函理论（DFT）中的B3LYP方法对CuCl2催化的（2-甲基辛烷-2,3-二烯-4-基）磷酸乙酯氯代环化反应机理进行了理论研究．在6-31＋G（d）基组水平上对反应机理中所有反应物、过渡态、中间体和产物进行了优化,通过能量和振动频率分析以及IRC计算证实了中间体和过渡态的合理性．在相同基组水平上应用自然键轨道（NBO）理论和分子中的原子（AIM）理论分析了复合物的成键特征和轨道间相互作用．反应物R和催化剂CuCl2可通过IA和IB两条可行反应通道生成中间体IM9,控制步骤活化能分别是129．61和142．10kJ／m01．中间体IM9到产物P也有两条反应路径PA和PB,控制步骤活化能分别是179．55和9．83kJ／m01．整个反应机理中IA—PB和IB—PB反应通道可能同时发生,反应控制步骤活化能最低反应通道为IA→PB．     

A hybrid potential reaction path and free energy study of the chorismate mutase reaction.

We present a combination of two techniques--QM/MM statistical simulation methods and QM/MM internal energy minimizations--to get a deeper insight into the reaction catalyzed by the enzyme chorismate mutase. Structures, internal energies and free energies, taken from the paths of the reaction in solution and in the enzyme have been analyzed in order to estimate the relative importance of the reorganization and preorganization effects. The results we obtain for this reaction are in good agreement with experiment and show that chorismate mutase achieves its catalytic efficiency in two ways; first, it preferentially binds the active conformer of the substrate and, second, it reduces the free energy of activation for the reaction relative to that in solution by providing an environment which stabilizes the transition state.     

Theoretical Study on the Reaction Mechanism of SiCl_4 with H in the Gas Phase

1 INTRODUCTION Semiconductor silicon materials are vital for mi- croelectronic and information industry. Silicon has many advantages, for example, rich resource, out- standing quality and sophisticated processing tech- nology. So it has been widely used in semiconduc- tor industry. One of the key techniques of mo- dern microelectronic industry is epitaxial growth of single crystal thin film on single crystal silicon and its ba- cking materials. In the chemical vapour deposition of Si, g…     

Parallel iterative reaction path optimization in ab initio quantum mechanical/molecular mechanical modeling of enzyme reactions

The determination of reaction paths for enzyme systems remains a great challenge for current computational methods. In this paper we present an efficient method for the determination of minimum energy reaction paths with the ab initio quantum mechanical/molecular mechanical approach. Our method is based on an adaptation of the path optimization procedure by Ayala and Schlegel for small molecules in gas phase, the iterative quantum mechanical/molecular mechanical (QM/MM) optimization method developed earlier in our laboratory and the introduction of a new metric defining the distance between different structures in the configuration space. In this method we represent the reaction path by a discrete set of structures. For each structure we partition the atoms into a core set that usually includes the QM subsystem and an environment set that usually includes the MM subsystem. These two sets are optimized iteratively: the core set is optimized to approximate the reaction path while the environment set is optimized to the corresponding energy minimum. In the optimization of the core set of atoms for the reaction path, we introduce a new metric to define the distances between the points on the reaction path, which excludes the soft degrees of freedom from the environment set and includes extra weights on coordinates describing chemical changes. Because the reaction path is represented by discrete structures and the optimization for each can be performed individually with very limited coupling, our method can be executed in a natural and efficient parallelization, with each processor handling one of the structures. We demonstrate the applicability and efficiency of our method by testing it on two systems previously studied by our group, triosephosphate isomerase and 4-oxalocrotonate tautomerase. In both cases the minimum energy paths for both enzymes agree with the previously reported paths.     

Theoretical studies of the NTO unimolecular decomposition

This work studies 39 decomposition paths among 18 intermediates and 14 transition states. Three types of intra-molecular proton migration and the direct scission of C–NO₂ were regarded as the initial steps of the unimolecular decomposition of NTO. The activation energies of the radicalization C–NO₂ homolysis step are 79.158, 79.781 and 80.652 kcal mol⁻¹. The activation energies of the ionization C–NO⁻¹₂ scission step are 262.488, 263.138 and 272.278 kcal mol⁻¹. The bottle neck activation energies of the C–NO₂H cleavage are 54.936, 63.257 and 71.247 kcal mol⁻¹. Two paths have the smallest bottle neck activation energy. Both of them have two proton migration steps and one internal rotation step prior to C–NO₂H cleavage. At lower temperatures, energy accumulated slowly. When the energy is high enough and reaction time is long enough for structure transformation, these two mechanisms should be the most probable decomposition paths. At high temperatures, the shortest (four steps) mechanism which goes through radicalization C–NO₂ scission should be the dominant path. There are five tautomers found in this study. Four of them are intra-molecular proton migration tautomers. The other one is an internal rotational tautomer. Their energy barriers for structure transfer are lower than any of the activation energies of the decomposition reactions. It may be regarded as one explanation of the insensitive property of NTO.     

The Role of Liquid Water in Crystalline Hydrate Dehydration: Copper sulphate pentahydrate

Simultaneous Differential Thermal Analysis/Thermogravimetric experiments carried out on one large single crystal, several small single crystals and powdered crystals of pentahydrate copper sulphate have been used to demonstrate the role that retained liquid water plays in maintaining crystal morphology during dehydration. Measured activation energies for stepwise dehydration in the system show the presence of solution-based transformations provide lower energy paths for the dehydration steps and stress relieving mechanisms. Skeletal anhydrous crystals from large-sized pentahydrate copper sulphate have the same morphology as the starting crystal on complete dehydration at controlled heating rates as long as a solution phase is maintained within the crystal during decomposition. The athermal activation energies, in this work, are in agreement with those obtained by isothermal methods as long as coincident reaction paths for the two techniques are maintained. The literature has been reviewed in the light of this work and a three-stage process is presented to rationalise the conflicting information obtained by workers using a variety of different experimental techniques.This revised version was published online in November 2005 with corrections to the Cover Date.     

Activation energy for a model ferrous-ferric half reaction from transition path sampling

Activation parameters for the model oxidation half reaction of the classical aqueous ferrous ion are compared for different molecular simulation techniques. In particular, activation free energies are obtained from umbrella integration and Marcus theory based thermodynamic integration, which rely on the diabatic gap as the reaction coordinate. The latter method also assumes linear response, and both methods obtain the activation entropy and the activation energy from the temperature dependence of the activation free energy. In contrast, transition path sampling does not require knowledge of the reaction coordinate and directly yields the activation energy [C. Dellago and P. G. Bolhuis, Mol. Simul. 30, 795 (2004)]. Benchmark activation energies from transition path sampling agree within statistical uncertainty with activation energies obtained from standard techniques requiring knowledge of the reaction coordinate. In addition, it is found that the activation energy for this model system is significantly smaller than the activation free energy for the Marcus model, approximately half the value, implying an equally large entropy contribution.     

Reaction Mechanism and Kinetics for HCCO Radical with NO

The mechanism and dynamical properties for the reaction of HCCO radicals with NO were investigated theoretically. The minimum energy paths(MEP) of the reaction were calculated by using the density functional theory(DFT) at the B3LYP/6-311 G^** level, and the energies along the MEP were further refined at the QCISD(T)/6-311 G^** level. It is found that the reaction mechanism of the title reaction involves three channels, producing HCNO CO, HONC CO and HCN CO2 products, respectively. Channel 1 is the most favorable path. The rate constant for channel 1 were calculated over a temperature range of 800-2500 K by using the canonical variational transition-state theory(CVT). The rate constant for the main path is negatively dependent on temperature, which is a characteristic of radical reactions with negative activation energy, and the variational effect for the rate constant calculation is small in the whole temperature range.     

化学反应的内禀反应坐标法VI. 甲胺脱氢反应的反应路径解析

本文用内禀反应坐标(IRC)法对甲胺的1,1-和1,2-脱氢反应途径进行了微观动力学解析。在RHF/4-31G水平上确定了反应的过渡态结构和势能变化曲线,得到了活化能、反应热、频率因子和活化熵等物理量。模式选择研究表明:1,1-和1,2-脱氢反应过程的反应坐标(IRC)分别与反应物中频率为1057和1348cm^-^1的两个变形振动模式相联结。在每条反应途径上都存在正则坐标间的振动耦合。本文的计算为激光诱导和模式选择化学反应的研究提供了重要信息和理论指导。     

硅氮烷添加剂的水解稳定性对硅橡胶热稳定性的影响

采用一种新的方法研究了硅氮烷水解反应的动力学规律 ,根据实验数据测出硅氮烷水解速率对硅氮烷的浓度符合一级动力学关系 .计算出了几种硅氮化合物与水蒸气反应的表观活化能 ,结果表明硅氮烷的结构与其水解反应的表观活化能有密切的关系 .硅原子或氮原子上带有较大空间位阻的基团后 ,其水解稳定性提高 .其中六苯基环三硅氮烷 4的水解表观活化能为 2 14kJ mol,而苯基硅氮聚合物 5的水解表观活化能更达到 2 91 3kJ mol.添加到硅橡胶中的硅氮化合物水解表观活化能越大 ,即水解稳定性越高 ,其改进硅橡胶热稳定性的效果越好 .将 4和 5添加到硅橡胶生胶中 ,35 0℃下老化 2 4h的热失重分别为 0 96 %和 0 6 % .     

The effect of external forces on the initial dissociation of RDX (1,3,5‐trinitro‐1,3,5‐triazine): A mechanochemical study

Experimental and theoretical studies have proposed different initiation reactions for the decomposition of hexahydro‐1,3,5‐trinitro‐1,3,5‐triazine (RDX). Three primary reactions are considered to start RDX decomposition: homolytic N? N bond fission, HONO elimination, and concerted fission of C? N bonds. The focus of this article is to study the effect of external forces on the energy barrier and reaction energies of all three mechanisms. We used the Nudged Elastic Band method along with ab initio Density Functional Theory within the framework of a generalized force‐modified potential energy surface (G‐FMPES) to calculate the minimum energy paths at different compressive (corresponding to pressure between approximately 6 and 294 MPa) and expansive force values (between 10 and 264 pN). For all three reactions, the application of an expansive force increases the exothermicity and lowers the energy barriers to different extents, while a compressive force decreases the exothermicity and raises the energy barrier to different extents.