Theoretical Study on the Mechanism of CF2 Reaction with CH2O

The insertion reaction mechanism of CF2 with CH2O was investigated at the B3LYP/6-311G（d）//MP2/6-311G（d） level. The geometric conformations at each stationary point in reaction potential surface were fully optimized and the transition states were verified by intrinsic reaction coordinate （IRC） and frequency analysis. The energies of all reactants were calculated with CCSD（T）/6-311G（d）//G2MP2 methods. Results indicated that the P1 reaction route with difuoroaldehyde as product is the dominant reaction pathway, which exhibits nucleophilic character. According to NBO analysis, the starting point of insertion reaction is the interaction between carbene LP（C3） and formaldehyde π（Cl-O2）. Besides, the thermodynamic and dynamic properties of dominated reaction （1） at different temperature were studied with statistic thermodynamic method and Eyring transition state theory adjusted by Wigner means, from which the proper temperature （500- 1200 K） of reaction （1） could be estimated. Finally, the thermo- dynamic and dynamic properties of insertion reaction mechanisms （CF2, CX2 （X = Cl, Br） with CH2O） were compared and discussed.     

Kinetics of the Reaction of Abstracting Hydrogen from Ethylene by Hydrogen Atom

Hydrogen abstraction reaction, H C2H4 --H2 C2H2 was studied by using A initio SCF method. Ge-ometries were fully optimized at SCF level and energies were computed at STO-3G basis set for reactants and transition state. Vibrational analysis was performed thereupon. Finally, the rate constant calculations were carried out at different temperatures for all range of reaction temperature according to Eyring's sbwlute reaction rate theory. The calculated activation energy is 12. 68 kcal/mol, lower than observed value (H. S kcal/mol) by 1. 82 kcal/mol only. The agreement of the calculated rate constants with the experiments is satisfactory.     

Theoretical Studies on Reaction Mechanisms of HNCS with NH （ X^3∑）

The reaction mechanisms of HNCS with NH（X^3∑ ） were theoretically investigated. The minimum energy paths （MEP） of the reaction were calculated by using the density functional theory（DFT） at the B3LYP/6-311 ＋ ＋ G^＊＊ level. The equilibrium structural parameters, the harmonic vibrational frequencies, the total energies, and the zeropoint energies（ZPE） of all the species were calculated. The single-point energies along the MEP were further refined at the QCISD（T）/6-311 ＋ ＋ G^＊ ＊ level. It was found that the mechanisms of the HNCS ＋ NH（X^3∑） reaction involve two channels producing the HNC ＋ HNS and the N2H2 ＋ CS products. Channel 1 plays a dominant role and the HNC ＋ HNS are the main preduets. The reaction is exothermie.     

邻位吡啶自由基多通道分解反应的动力学和反应机理研究

Utilizing Gaussian94 program package, all species involved in decomposition reactions of o-pyridyl radical were optimized fully at B3LYP/6-311＋＋G^＊＊ level. Intrinsic reaction coordinate calculations were employed to confirm the connections of the transition states and products, and transition states were ascertained by the number of imaginary frequency （0 or 1）. The reaction mechanism was elucidated by the vibrational mode analysis and electronic population analysis, and the reaction rate constants were calculated with transition state theory.     

Ab initio Study of the Potential Energy Surface and Product Branching Ratios for Reaction of O（^1D） with C2H5CI

The potential energy surface of O（^1D）＋C2H5Cl reaction was studied using QCISD（T）/6- 311＋＋G（d,p）//MP2/6-31G（d,p） method. The calculations reveal an insertion-elimination mechanism. The insertion reaction of O（^1D） and C2H5Cl produces two energy-rich intermediates, IM1 and IM2, which subsequently decompose into various products. The calculations of the branching ratios of various products formed through the two intermediates were carried out using RRKM （Rice-Ramsperger-Kassel-Marcus） theory at the collision energies of 0, 20.9, 41.8, 62.7, 83.6, 104.5, and 125.4 kJ/mol. HCl is the main decomposition product for IM1; CH2OH is the main decomposition product for IM2. Since IM1 is more stable than IM2, HCl is probably the main product of the O（^1D）＋C2H5Cl reaction.     

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.     

Theoretical study on the reaction of VS^＋（^3∑^-, 1^Г） with CO

Two possible reaction mechanisms of VS^＋（^3∑^-, 1^Г） with CO in the gas phase have been studied by using B3LYP/TZVP and CCSD（T）/6-311＋G （3df, 3pd） methods： the O/S exchange reaction （VS^＋＋CO→VO^＋＋CS） and the S-transfer reaction （VS^＋ ＋ CO → V^＋ ＋ COS）. The two reactions proceed via two-step and one-step mechanism, respectively. The barriers of the triplet and singlet PESs are 30.6 and 50.9 kcal/mol, respectively, for O/S exchange reaction and 7.3 and 50.2 kcal/mol, respectively, for the S-transfer reaction. The results indicate that the triplet ground state reaction is more favorable, and the S-transfer reaction is more favorable than the O/S exchange reaction, which is in good agreement with the experimental observation.     

DFT Study on Intermolecular Cleavage Reaction of N-（2-Hydroxyphenyl）-phthalamic Acid

The reaction mechanisms of intermolecular cleavage reaction of N-（2-hydroxyphenyl）-phthalamic acid were studied via the density functional theory（DFT）. All geometries of the reactant, transition states, and products were optimized at the B3LYP/6-31G（d, p） level. Vibration analysis was carried out to confirm its identity as transitions＇ structure, and the intrinsic reaction coordinate method（IRC） was used to search the minimum energy path. Two possible reaction channels are reported in this article. The calculated results indicate that O-cyclization reaction channel has the lower activation barrier, and therefore, it occurs more easier than the other.     

Theoretical Study on the Mechanism of the Gas-phase Reaction of Sc＋ with Propargyl Alcohol

In order to elucidate the reaction mechanisms of reaction Sc^＋ with propargyl alcohol (PPA), the triplet potential energy surface for the reactions has been theoretically investigated using a DFT method. The geometries for the reactants, intermediates, transition states and products were completely optimized at B3LYP/DZVP level. The single point energy of each stationary point was calculated at MP4/(6-311＋G** for C, H, O and Lanl2dz for Sc^＋) level. All the transition states were verified by the vibrational analysis and the internal reaction coordinate (IRC) calculations. The present results show that the reaction takes an insertion-elimination mechanism both along the O—H and C—O bond activation branches, but the C—O bond activation is much more favorable in energy than the O—H bond activation. All theoretical results not only support the existing conclusions inferred from early experiment, but also complement the pathway and mechanism for this reaction.     

Thermochemical Properties and Decomposition Kinetics of Ammonium Magnesium Phosphate Monohydrate

Ammonium magnesium phosphate monohydrate NH4MgPO4·H2O was prepared via solid state reaction at room temperature and characterized by XRD, FT-IR and SEM. Thermochemical study was performed by an isoperibol solution calorimeter, non-isothermal measurement was used in a multivariate non-linear regression analysis to determine the kinetic reaction parameters. The results show that the molar enthalpy of reaction above is （28.795 ± 0.182） kJ/mol （298.15 K）, and the standard molar enthalpy of formation of the title complex is （-2185.43 ± 13.80） kJ/mol （298.15 K）. Kinetics analysis shows that the second decomposition of NH4MgPO4·H2O acts as a double-step reaction： an nth-order reaction （Fn） with n=4.28, E1=147.35 kJ/mol, A1=3.63×10^13 s^-1 is followed by a second-order reaction （F2） with E2=212.71 kJ/mol, A2= 1.82 × 10^18 s^-1.     

SiH3+O(3P)反应机理的理论研究

The reaction for SiH3+O(3P) was studied by ab initio method. The geometries of the reactants, intermediates, transition states and products were optimized at MP2/6-311+G(d,p) level. The single point calculations for all the stationary points were carried out at the QCISD(T) /6-311+G(d,p) level using the MP2/6-311+G(d,p) optimized geometries. The results of the theoretical study indicate that the major pathway is the SiH3+O(3P)→IM1→TS3→IM2→TS8→HOSi+H2. The other minor products include the HSiOH+H, H2SiO+H and HSiO+H2. Furthermore, the products HOSi, HSiO and HSiOH(cis) can undergo dissociation into the product SiO. In addition, the calculations provide a possible interpretation for disagreement about the mechanism of the reaction SiH4+O(3P). It suggests that the products HSiOH, H2SiO and SiO observed by Withnall and Andrews are produced from the secondary reaction SiH3+O(3P) and not from the reaction SiH4+O(3P).     

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.     

Thermodynamics of boroxine formation from the aliphatic boronic acid monomers R-B(OH)2 (R = H, H3C, H2N, HO, and F): a computational investigation

Boroxines are the six-membered cyclotrimeric dehydration products of organoboronic acids, 3R–B(OH)2 → R3B3O3 + 3H2O, and in recent years have emerged as a useful class of organoboron molecules with applications in organic synthesis both as reagents and catalysts, as structural components in boronic-acid-derived pharmaceutical agents, and as anion acceptors and electrolyte additives for battery materials [Korich, A. L.; Iovine, P. M. Dalton Trans. 2010, 39, 1423?1431]. Second-order M?ller–Plesset perturbation theory, in conjunction with the Dunning–Woon correlation-consistent cc-pVDZ, aug-cc-pVDZ, cc-pVTZ, and aug-cc-pVTZ basis sets, was used to investigate the structures and relative energies of the endo–exo, anti, and syn conformers of the aliphatic boronic acids R–B(OH)2 (R = H, H3C, H2N, HO, and F), as well as the thermodynamics of their boroxine formation; single-point calculations at the MP2/aug-cc-pVQZ, MP2/aug-cc-pV5Z, and CCSD(T)/aug-cc-pVTZ levels using the MP2/aug-cc-pVTZ optimized geometries were also performed in selected cases. The endo–exo conformer was generally lowest in energy in vacuo, as well as in PCM and CPCM models of aqueous and carbon tetrachloride media. The values of ΔH(298)(0) for boroxine formation via dehydration from the endo–exo conformers of these aliphatic boronic acids ranged from ?2.9 for (H2N)3B3O3 to +12.2 kcal/mol for H3B3O3 at the MP2/aug-cc-pVTZ level in vacuo; for H3B3O3, the corresponding values in PCM/UFF implicit carbon tetrachloride and aqueous media were +11.2 and +9.8 kcal/mol, respectively. On the basis of our calculations, we recommend that ΔHf(298K) for boroxine listed in the JANAF compilation needs to be revised from ?290.0 to approximately ?277.0 kcal/mol.     

CH3CH2+O(3P)反应机理的理论研究

The reaction for CH3CH2+O(3P) was studied by ab initio method. The geometries of the reactants, intermediates, transition states and products were optimized at MP2/6-311+G(d,p) level. The corresponding vibration frequencies were calculated at the same level. The single-point calculations for all the stationary points were carried out at the QCISD(T)/6-311+G(d,p) level using the MP2/6-311+G(d,p) optimized geometries. The results of the theoretical study indicate that the major products are the CH2O＋CH3, CH3CHO+H and CH2CH2+OH in the reaction. For the products CH2O＋CH3 and CH3CHO+H, the major production channels are A1: (R)→IM1→TS3→(A) and B1: (R)→IM1→TS4→(B), respectively. The majority of the products CH2CH2+OH are formed via the direct abstraction channels C1 and C2: (R)→TS1(TS2)→(C). In addition, the results suggest that the barrier heights to form the CO reaction channels are very high, so the CO is not a major product in the reaction.     

C3H2环丙烯基自由基与O(3P)反应机理的理论研究

本文用量子化学密度泛函方法对C3H2 (环丙烯基自由基)与O(3P)反应的机理进行了理论研究。在B3LYP/6-311++G**计算水平上优化了各驻点（过渡态,中间体,产物）的几何结构,在QCISD(T)/6-311++G**水平下计算了各物质的单点能量,在两种水平下计算了298K和600K时的能量。计算结果表明：C3H2 + O(3P) 反应可以生成P1 (C2H +HCO),P2 (C2H2 + CO) 和P3 (HC3O+H)三种产物。生成P1反应通道的能垒最低,即P1为主要产物,与实验的结果一致。产物P1可以通过路径：R→ IM1→ IM2→ P1获得。本文详细地讨论了C3H2 + O(3P) 的反应机理,并从理论上对实验结果进行了验证。研究结果有助于深入理解C3H2 + O(3P)反应机理以及C3H2在大气中的燃烧过程。     

CO oxidation catalyzed by a single gold atom: benchmark calculations and the performance of DFT methods

Quantum chemical calculations were carried out on CO oxidation catalyzed by a single gold atom. To investigate the performance of density functional theory (DFT) methods, 42 DFT functionals have been evaluated and compared with high-level wavefunction based methods. It was found that in order to obtain accurate results the functionals used must treat long range interaction well. The double-hybrid mPW2PLYP and B2PLYP functionals are the two functionals with best overall performance. CAM-B3LYP, a long range corrected hybrid GGA functional, also performs well. On the other hand, the popular B3LYP, PW91, and PBE functionals do not show good performance and the performance of the latter two are even at the bottom of the 42 functionals. Our accurate results calculated at the CCSD(T)/aug-cc-pVTZ//mPW2PLYP/aug-cc-pVTZ level of theory indicate that Au atom is a good catalysis for CO oxidation. The reaction follows the following mechanism where CO and O(2) adsorb on Au atom forming an Au(OCOO) intermediate and subsequently O(2) transfer one oxygen atom to CO to form CO(2) and AuO. Then AuO reacts with CO to form another CO(2) to complete the catalytic cycle. The overall energy barrier at 0 K for the first CO oxidation step (Au + CO + O(2)→ AuO + CO(2)) is just 4.8 kcal mol(-1), and that for the second CO oxidation step (AuO + CO → Au + CO(2)) is just 1.6 kcal mol(-1).     

4-氟苯甲醛、β-萘胺和Meldrum酸-锅反应生成1-(4-氟苯基)-1,2-二氢苯并[f]喹啉-3(4H)-酮的反应机理理论研究

采用密度泛函理论(DFT)研究了4-氟苯甲醛、β-萘胺和Meldrum酸一锅反应生成1-(4-氟苯基)-1,2-二氢苯并[f]喹啉-3(4H)-酮的微观反应机理.在B3LYP/6-311G*基组水平上优化了反应物、过渡态、中间体及产物的几何构型,通过振动分析确认了过渡态的结构,并用内禀反应坐标(IRC)确认反应途径.应用分子中的原子理论(AIM)分析了这些物质的成键特征.采用SCRF(PCM)方法研究了反应体系的溶剂化效应.报道了可能的反应路径,其中Re→TS1→IM1→TS2→IM2→TS3→IM3→TS4→IM5→TS7→IM9→TS13→IM10→TS14→P3具有相对较低的活化能,是反应的主要通道,理论预测的主要产物与实验吻合.     

n(H_2O)(n=1,2)在HO_2+NO→HNO_3反应中的催化机制研究

采用CCSD(T)/aug-cc-p VTZ//B3LYP/6-311+G(2df,2p)方法对n(H_2O)(n=0,1,2)参与HO_2+NO→HNO_3反应的微观机理和速率常数进行了研究.结果表明,由于水分子与HO_2形成的复合物(H_2O…HO_2,HO_2…H_2O)结合NO与水分子形成的复合物(NO…H_2O,ON…H_2O)的反应方式具有较高能垒和较低有效速率,其对HO_2+NO→HNO_3反应的影响远小于双体水(H_2O)2与HO_2(或NO)形成复合物然后再与另一分子反应物NO(或HO_2)的反应方式,因此n(H_2O)(n=1,2)催化HO_2+NO→HNO_3反应主要经历了HO_2…(H_2O)_n(n=1,2)+NO和NO…(H_2O)_n(n=1,2)+HO_22种反应类型.由于HO_2…(H_2O)_n(n=1,2)+NO反应的低能垒和高速率,HO_2…(H_2O)_n(n=1,2)+NO反应优于NO…(H_2O)_n(n=1,2)+HO_2反应.与此同时,由于计算温度范围内HO_2…H_2O+NO反应的有效速率常数比HO_2…(H_2O)2+NO反应对应的有效速率常数大了10~12数量级,可推测(H_2O)_n(n=1,2)催化HO_2+NO→HNO_3反应主要来自于单个水分子.此外,在216.7~298.6 K范围内水分子对HO_2+NO→HNO_3反应起显著的正催化作用,且随温度的升高有明显增大的趋势,在298.2 K时增强因子k'RW1/ktotal达到67.93%,表明在实际大气环境中水蒸气对HO_2+NO→HNO_3反应具有显著影响.     

C2H与HO2双自由基反应的密度泛函理论研究

应用量子化学从头算和密度泛函理论(DFT)对C₂H与HO₂双自由基的单重态反应进行了研究.在UB3LYP/6-311G水平上优化了反应通道上各驻点(反应物、中间体、过渡态和产物)的几何构型.在CCSD(T)/6-311G^{**水平上计算了各物种的单点能,并对总能量进行了零点能校正.研究结果表明,反应物中自由基C₂H的边端C进攻自由基HO₂的边端O是主要的进攻方式.首先形成了中间体1(HCCOOH),由此经过不同的反应通道可以得到主要产物P1,次要产物P2,P3和P5.生成P1的反应热为-814.40kJ/mol.自由基C₂H的中间C进攻自由基HO₂的边端O是次要的进攻方式,可以得到产物P4和P6.根据势能面分析,所有反应均是放热反应.}     

High-level direct-dynamics variational transition state theory calculations including multidimensional tunneling of the thermal rate constants, branching ratios, and kinetic isotope effects of the hydrogen abstraction reactions from methanol by atomic hydrogen

We report a detailed theoretical study of the hydrogen abstraction reaction from methanol by atomic hydrogen. The study includes the analysis of thermal rate constants, branching ratios, and kinetic isotope effects. Specifically, we have performed high-level computations at the MC3BB level together with direct dynamics calculations by canonical variational transition state theory (CVT) with the microcanonically optimized multidimensional tunneling (μOMT) transmission coefficient (CVT/μOMT) to study both the CH(3)OH+H→CH(2)OH+H(2) (R1) reaction and the CH(3)OH+H→CH(3)O+H(2) (R2) reaction. The CVT/μOMT calculations show that reaction R1 dominates in the whole range 298≤T (K)≤2500 and that anharmonic effects on the torsional mode about the C-O bond are important, mainly at high temperatures. The activation energy for the total reaction sum of R1 and R2 reactions changes substantially with temperature and, therefore, the use of straight-line Arrhenius plots is not valid. We recommend the use of new expressions for the total R1 + R2 reaction and for the R1 and R2 individual reactions.