Cl_2 2HI=2HCl I_2反应机理的理论研究

分别在MP2/3-21G!!、CCSD(T)/3-21G!!//MP2/3-21G!!和B3LYP/3-21G!!3种水平上,计算研究了气相反应Cl2 2HI=2HCl I2的机理,求得一系列四中心和三中心的过渡态.通过比较六种反应通道的活化能大小,得到了相同的结论:双分子基元反应Cl2 HI"HCl ICl和ICl HI"I2 HCl的最小活化能小于Cl2、HI和ICl的解离能,从理论上证明了反应Cl2 2HI=2HCl I2将优先以分子与分子作用形式分两步完成.用内禀反应坐标(IRC)验证了MP2/3-21G!!方法计算得到的过渡态.     

Cl2+2HI=2HCl+I2反应机理的理论研究

分别在MP2/3-21G**、CCSD(T)/3-21G**//MP2/3-21G**和B3LYP/3-21G**3种水平上, 计算研究了气相反应Cl₂+2HI=2HCl+I₂的机理, 求得一系列四中心和三中心的过渡态. 通过比较六种反应通道的活化能大小, 得到了相同的结论：双分子基元反应Cl₂+HIHCl+ICl和ICl+HII₂+HCl的最小活化能小于Cl₂、HI和ICl的解离能, 从理论上证明了反应Cl₂+2HI=2HCl+I₂将优先以分子与分子作用形式分两步完成. 用内禀反应坐标(IRC)验证了MP2/3-21G**方法计算得到的过渡态.     

DFT studies of the phosphinidene derivative HPNaF and its insertion reaction with R–H(R=F,OH, NH<Subscript>2</Subscript>, CH<Subscript>3</Subscript>)

The formations of the phosphinidene derivative HPNaF and its insertion reactions with R–H (R=F, OH, NH₂, CH₃) have been systematically investigated employing the density functional theory (DFT), such as the B3LYP and MPW1PW91 methods. A comparison with the results of MP2 calculations shows that MPW1PW91 underestimates the barrier heights for the four reactions considered. Similarly, the same is also true for the B3LYP method depending on the selected reactions, but by much less than MPW1PW91, where the barrier heights of the four reactions are 25.2, 85.7, 119.0, and 142.4 kJ/mol at the B3LYP/6-311+G* level of theory, respectively. All the mechanisms of the four reactions are identical to each other, i.e., an intermediate has been located during the insertion reaction. Then, the intermediate could dissociate to substituted phosphinidane(H₂RP) and NaF with a barrier corresponding to their respective dissociation energies. Correspondingly, the reaction energies for the four reactions are −92.2, −68.1, −57.2, and −44.3 kJ/mol at the B3LYP/6-311+G* level of theory, respectively, where both the B3LYP and MPW1PW91 methods underestimate the reaction energies compared with the MP2 results. The linear correlations between the calculated barrier heights and the reaction energies have also been observed. As a result, the relative reactivity among the four insertion reactions should be as follows: H–F > H–OH > H–NH₂ > H–CH₃.     

Mechanisms for the deamination reaction of cytosine with H2O/OH(-) and 2H2O/OH(-): a computational study

Mechanisms for the deamination reaction of cytosine with H 2O/OH (-) and 2H 2O/OH (-) to produce uracil were investigated using ab initio calculations. Optimized geometries of reactants, transition states, intermediates, and products were determined at MP2 and B3LYP using the 6-31G(d) basis set and at B3LYP/6-31+G(d) levels of theory. Single point energies were also determined at MP2/G3MP2Large and G3MP2 levels of theory. Thermodynamic properties (Delta E, Delta H, and Delta G), activation energies, enthalpies, and free energies of activation were calculated for each reaction pathway investigated. Intrinsic reaction coordinate (IRC) analysis was performed to characterize the transition states on the potential energy surface. Seven pathways for the deamination reaction were found. All pathways produce an initial tetrahedral intermediate followed by several conformational changes. The final intermediate for all pathways dissociates to product via a 1-3 proton shift. The activation energy for the rate-determining step, the formation of the tetrahedral intermediate for pathway D, the only pathway that can lead to uracil, is 115.3 kJ mol (-1) at the G3MP2 level of theory, in excellent agreement with the experimental value (117 +/- 4 kJ mol (-1)).     

Computational methods to calculate accurate activation and reaction energies of 1,3-dipolar cycloadditions of 24 1,3-dipoles

Theoretical calculations were performed on the 1,3-dipolar cycloaddition reactions of 24 1,3-dipoles with ethylene and acetylene. The 24 1,3-dipoles are of the formula X≡Y(+)-Z(-) (where X is HC or N, Y is N, and Z is CH(2), NH, or O) or X═Y(+)-Z(-) (where X and Z are CH(2), NH, or O and Y is NH, O, or S). The high-accuracy G3B3 method was employed as the reference. CBS-QB3, CCSD(T)//B3LYP, SCS-MP2//B3LYP, B3LYP, M06-2X, and B97-D methods were benchmarked to assess their accuracies and to determine an accurate method that is practical for large systems. Several basis sets were also evaluated. Compared to the G3B3 method, CBS-QB3 and CCSD(T)/maug-cc-pV(T+d)Z//B3LYP methods give similar results for both activation and reaction enthalpies (mean average deviation, MAD, < 1.5 kcal/mol). SCS-MP2//B3LYP and M06-2X give small errors for the activation enthalpies (MAD < 1.5 kcal/mol), while B3LYP has MAD = 2.3 kcal/mol. SCS-MP2//B3LYP and B3LYP give the reasonable reaction enthalpies (MAD < 5.0 kcal/mol). The B3LYP functional also gives good results for most 1,3-dipoles (MAD = 1.9 kcal/mol for 17 common 1,3-dipoles), but the activation and reaction enthalpies for ozone and sulfur dioxide are difficult to calculate by any of the density functional methods.     

Theoretical Study on Gas—phase Pyrolytic Reactions of N—Ethyl,Isopropyl and N—t—Butyl Substituted 2—Aminopyrazine

Density functional theory (DFT) and ab initio methods were used to study gas‐phase pyrolytic reaction mechanisms of iV‐ethyl, N‐isopropyl and N‐t‐butyl substituted 2‐aminopyrazine at B3LYP/6–31G* and MP2/6–31G*, respectively. Single‐point energies of all optimized molecular geometries were calculated at B3LYP/6–311 + G(2d,p) level. Results show that the pyrolytic reactions were carried out through a unimolecular first‐order mechanism which were caused by the migration of atom H(17) via a six‐member ring transition state. The activation energies which were verified by vibrational analysis and correlated with zero‐point energies along the reaction channel at B3LYP/6–311 + G(2d,p) level were 252.02 kJ. mo^?1 (N‐ethyl substituted), 235.92 kJ‐mol^?1 (N‐t‐isopropyl substituted) and 234.27 kJ‐mol^?1 (N‐t‐butyl substituted), respectively. The results were in good agreement with available experimental data.     

Transition state structure and energetics of the N(2)O + X (X = Cl,Br) reactions

The structural and vibrational properties of the transition state of the N(2)O + X (X = Cl,Br) reactions have been characterized by ab initio methods using density functional theory. We have employed Becke's hybrid functional (B3LYP), and transition state optimizations were performed with 6-31G(d), 6-311G(2d,2p), 6-311+G(3d,2p), and 6-311+G(3df,2p) basis sets. For the chlorine atom reaction the coupled-cluster method (CCSD(T)) with 6-31G(d) basis set was also used. All calculations resulted in transition state structures with a planar cis arrangement of atoms for both reactions. The geometrical parameters of transition states at B3LYP are very similar, and the reaction coordinates involve mainly the breaking of the N-O bond. At CCSD(T)/6-31G(d) level a contribution of the O-Cl forming bond is also observed in the reaction coordinate. In addition, several highly accurate ab initio composite methods of Gaussian-n (G1, G2, G3), their variations (G2(MP2), G3//B3LYP), and complete basis set (CBS-Q, CBS-Q//B3LYP) series of models were applied to compute reaction energetics. All model chemistries predict exothermic reactions. The G3 and G2 methods result in the smallest deviations from experiment, 1.8 and 0 kcal mol(-1), for the enthalpies of reaction for N(2)O reaction with chlorine and bromine, respectively. The G3//B3LYP and G1 methods perform best among the composite methods in predicting energies of the transition state, with a deviation of 1.9 and 3.0 kcal mol(-1), respectively, in the activation energies for the above processes. However, the B3LYP/6-311+G(3df,2p) method gives smaller deviations of 0.4 and -1.0 kcal mol(-1), respectively. The performance of the methodologies applied in predicting transition state energies was analyzed.     

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.     

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.     

Addition reaction of adamantylideneadamantane with Br2 and 2Br2: a computational study

Ab initio calculations were carried out for the reaction of adamantylideneadamantane (Ad=Ad) with Br2 and 2Br2. Geometries of the reactants, transition states, intermediates, and products were optimized at HF and B3LYP levels of theory using the 6-31G(d) basis set. Energies were also obtained using single point calculations at the MP2/6-31G(d)//HF/6-31G(d), MP2/6-31G(d)//B3LYP/6-31G(d), and B3LYP/6-31+G(d)//B3LYP/6-31G(d) levels of theory. Intrinsic reaction coordinate (IRC) calculations were performed to characterize the transition states on the potential energy surface. Only one pathway was found for the reaction of Ad=Ad with one Br2 producing a bromonium/bromide ion pair. Three mechanisms for the reaction of Ad=Ad with 2Br2 were found, leading to three different structural forms of the bromonium/Br3- ion pair. Activation energies, free energies, and enthalpies of activation along with the relative stability of products for each reaction pathway were calculated. The reaction of Ad=Ad with 2Br2 was strongly favored over the reaction with only one Br2. According to B3LYP/6-31G(d) and single point calculations at MP2, the most stable bromonium/Br3- ion pair would form spontaneously. The most stable of the three bromonium/Br3- ion pairs has a structure very similar to the observed X-ray structure. Free energies of activation and relative stabilities of reactants and products in CCl4 and CH2ClCH2Cl were also calculated with PCM using the united atom (UA0) cavity model and, in general, results similar to the gas phase were obtained. An optimized structure for the trans-1,2-dibromo product was also found at all levels of theory both in gas phase and in solution, but no transition state leading to the trans-1,2-dibromo product was obtained.     

Predicting Rate Constants for Nucleophilic Reactions of Amines with Diarylcarbenium Ions Using an ONIOM Method

The rate constants of the nucleophilic reactions between amines and benzhydrylium ions were calculated using first-principles theoretical methods. Solvation models including PCM, CPCM, and COSMORS, as well as different types of atomic radii including UA0, UAKS, UAHF, Bondi, and UFF, and several single-point energy calculation methods (B3LYP, B3P86, B3PW91, BHANDH, PBEPBE, BMK, M06, MP2, and ONIOM method) were exam-ined. By comparing the correlation between experimental rate constants and the calculated values, the ONIOM(CCSD(T)/6-311++G(2df,2p):B3LYP/6-311++G(2df,2p))// B3LYP/6-31G(d)/PCM/UFF) method was found to perform the best. This method was then employed to calculate the rate constants of the reactions between diverse amines and diarylcarbenium ions. The calculated rate constants for 65 reactions of amines with diarylcarbenium ions are in agreement with the experimental values, indicating that it is feasible to predict the rate constant of a reaction between an amine and a diarylcarbenium ion through ab initio calculation.     

Theoretical studies of the reactions of CF3CHCLOCHF2/CF3CHFOCHF2 with OH radical and Cl atom and their product radicals with OH

The mechanisms and dynamics studies of the OH radical and Cl atom with CF(3)CHClOCHF(2) and CF(3)CHFOCHF(2) have been carried out theoretically. The geometries and frequencies of all the stationary points are optimized at the B3LYP/6-311G(d,p) level, and the energy profiles are further refined by interpolated single-point energies (ISPE) method at the G3(MP2) level of theory. For each reaction, two H-abstraction channels are found and four products (CF(3)CHFOCF(2), CF(3)CFOCHF(2), and CF(3)CHClOCF(2), CF(3)CClOCHF(2)) are produced during the above processes. The rate constants for the CF(3)CHClOCHF(2)/CF(3)CHFOCHF(2) + OH/Cl reactions are calculated by canonical variational transition-state theory (CVT) within 200-2000 K, and the small-curvature tunneling is included. The total rate constants calculated from the sum of the individual rate constants and the branching ratios are in good agreement with the experimental data. The Arrhenius expressions for the reactions are obtained. Our calculation shows that the substitution of Cl by F decreases the reactivity of CF(3)CHClOCHF(2) toward OH and Cl. In addition, the mechanisms of subsequent reactions of product radicals and OH radical are further investigated at the G3(MP2)//B3LYP/6-311G(d,p) level, and the main products are predicted in the this article.     

FC(O)O自由基与NO2反应的量子化学研究

用量子化学DFT, MP2, G3和G3MP2方法对FC(O)O自由基与NO₂的反应机理进行了理论研究. 优化了反应势能面上各驻点的几何结构, 通过内禀反应坐标(IRC)计算和振动分析, 确认了反应中的过渡态, 并用过渡态理论(TST)计算了相关反应的速率常数.     

Experimental and theoretical investigation of the decomposition of lithiated hydroxyl side-chain amino acids

Lithium cation complexes with serine (Ser) and threonine (Thr) are collisionally activated with xenon in a guided ion beam tandem mass spectrometer and are observed to exhibit a variety of decomposition pathways in addition to a loss of the intact ligand. Prominent pathways include a loss of H2O, CO2, and aldehydes (XCHO where X=H for Ser and CH3 for Thr). Quantum chemical calculations at the B3LYP/6-311+G(d,p) level are used to explore the reaction mechanisms for these processes in detail. Complete potential energy surfaces for all three processes are elucidated, including all intermediates and transition states. Theoretical molecular parameters for the rate-limiting transition states are then used to analyze the threshold energies in the experimental data, providing experimental measurements of the energies of these transition states. These experimental energies are compared with single-point energies calculated at three different levels, B3LYP, B3P86, and MP2(full), using the 6-311+G(2d,2p) basis set with geometries and zero-point energies calculated at the B3LYP/6-311+G(d,p) level. Good agreement between experiment and theory (especially MP2(full)) suggests that the reaction mechanisms have been reasonably elucidated.     

How computational methods and relativistic effects influence the study of chemical reactions involving Ru‐NO complexes?

Two treatments of relativistic effects, namely effective core potentials (ECP) and all‐electron scalar relativistic effects (DKH2), are used to obtain geometries and chemical reaction energies for a series of ruthenium complexes in B3LYP/def2‐TZVP calculations. Specifically, the reaction energies of reduction ( A ‐ F ), isomerization ( G‐I ), and Cl⁻ negative trans influence in relation to NH₃ ( J ‐ L ) are considered. The ECP and DKH2 approaches provided geometric parameters close to experimental data and the same ordering for energy changes of reactions A ‐ L . From geometries optimized with ECP, the electronic energies are also determined by means of the same ECP and basis set combined with the computational methods: MP2, M06, BP86, and its derivatives, so as B2PLYP, LC‐wPBE, and CCSD(T) (reference method). For reactions A ‐ I , B2PLYP provides the best agreement with CCSD(T) results. Additionally, B3LYP gave the smallest error for the energies of reactions J ‐ L . © 2017 Wiley Periodicals, Inc.     

Electronic structure benchmark calculations of CO2 fixing elementary chemical steps in RuBisCO using the projector-based embedding approach

Ribulose 1,5-bisphosphate carboxylase-oxygenase (RuBisCO) is the main enzyme involved in atmospheric carbon dioxide (CO₂) fixation in the biosphere. This enzyme catalyzes a set of five chemical steps that take place in the same active-site within magnesium (II) coordination sphere. Here, a set of electronic structure benchmark calculations have been carried out on a reaction path proposed by Gready et al. by means of the projector-based embedding approach. Activation and reaction energies for all main steps catalyzed by RuBisCO have been calculated at the MP2, SCS-MP2, CCSD, and CCSD(T)/aug-cc-pVDZ and cc-pVDZ levels of theory. The treatment of the magnesium cation with post-HF methods is explored to determine the nature of its involvement in the mechanism. With the high-level ab initio values as a reference, we tested the performance of a set of density functional theory (DFT) exchange-correlation (xc) functionals in reproducing the reaction energetics of RuBisCO carboxylase activity on a set of model fragments. Different DFT xc-functionals show large variation in activation and reaction energies. Activation and reaction energies computed at the B3LYP level are close to the reference SCS-MP2 results for carboxylation, hydration and protonation reactions. However, for the carbon–carbon bond dissociation reaction, B3LYP and other functionals give results that differ significantly from the ab initio reference values. The results show the applicability of the projector-based embedding approach to metalloenzymes. This technique removes the uncertainty associated with the selection of different DFT xc-functionals and so can overcome some of inherent limitations of DFT calculations, complementing, and potentially adding to modeling of enzyme reaction mechanisms with DFT methods.     

2 + 4] and [4 + 2] cycloadditions of o-thioquinones with 1,3-dienes: a computational study

Cycloadditions of o-thioquinones (o-TQs) with 1,3-dienes could proceed via either a [2 + 4] or a [4 + 2] mechanism. Under kinetic control and with acyclic dienes the reaction affords the spiro cycloadducts 5deriving from the [2 + 4] path as the main products. Under thermodynamic control, or with cyclic dienes, the o-TQs behave as heterodienes to give the benzoxathiin derivatives 4, in most cases with complete regioselectivity. In the present computational study, DFT calculations were performed in order to achieve a deep understanding of both [2 + 4] and [4 + 2] paths. The reactions of three o-TQs with six 1,3-dienes were thoroughly investigated at the B3LYP/TZVP//B3LYP6-31G level, and the two reaction mechanisms were then compared, evidencing that [2 + 4] cycloadditions are kinetically favored, strongly asynchronous, or even unconcerted, while [4 + 2] reactions are thermodynamically favored, quite asynchronous, but undoubtedly concerted. Moreover, the observed regioselectivity was rationalized by mean of the FMO theory and by comparison of the activation energies for different pathways.     

Theoretical investigation of clusters of phosphorus and arsenic: fascination and temptation of high symmetries

We present a theoretical study of the energetic and thermodynamic stability of selected phosphorus and arsenic clusters containing 18 to 168 atoms. For this purpose we employ MP2 as well as DFT functionals BP86 and B3LYP with extended basis sets. All procedures predict the family of one-dimensional polymers X18+12n, each with 2n-1 isomers of virtually identical energy, to be more stable than other structures investigated so far. Furthermore, islands of stability result for ring-shaped clusters X24n with Dnd symmetry for n=4 (only for arsenic), 5, 6, and 7. Phosphorus and arsenic show otherwise a very similar behavior. An investigation of basis set effects shows that a doubly polarized triple zeta valence basis (TZVPP) is both necessary and sufficient. In comparison to the reliable spin component scaled MP2 (SCS-MP2) procedure, DFT methods underestimate and MP2 overestimates the stability of larger clusters; the discrepancy increases with the number of atoms. The addition of a long-range dispersion correction to B3LYP energies does not rectify the shortcomings of DFT in comparison with SCS-MP2.     

F2＋2HI=2HF＋I2反应机理的研究

根据氯离子型层状复合氢氧化物(LDH-Cl)制备过程中溶液浓度变化的监测结果和不同反应进程时产物的EDS、IR、XRD、TEM、TG-DTA表征结果,研究了合成LDH-Cl的共沉淀反应动力学特征及机理.实验结果表明, LDH-Cl的生成符合多核层表面反应动力学模型;反应过程中LDH的晶胞参数c从2.421 nm变为2.399 nm,通道高度h由0.3321 nm减小为0.3228 nm,粒子直径Da由6.40 nm增大为15.16 nm, Dc由7.43 nm增大到10.93 nm,纵横比由0.86增大为1.39; IR和TG-DTA特征变化表明了层板对阴离子作用的强度和层板的结构稳定性随反应进程而提高.