Kinetic study of the reaction H2O2+H→H2O+OH by ab initio and density functional theory calculations

Theoretical investigations on the kinetics of the elementary reaction H₂O₂+H→H₂O+OH were performed using the transition state theory (TST). Ab initio (MP2//CASSCF) and density functional theory (B3LYP) methods were used with large basis set to predict the kinetic parameters; the classical barrier height and the pre-exponential factor. The ZPE and BSSE corrected value of the classical barrier height was predicted to be 4.1 kcal mol⁻¹ for MP2//CASSCF and 4.3 kcal mol⁻¹ for B3LYP calculations. The experimental value fitted from Arrhenius expressions ranges from 3.6 to 3.9 kcal mol⁻¹. Thermal rate constants of the title reaction, based on the ab initio and DFT calculations, was evaluated for temperature ranging from 200 to 2500 K assuming a direct reaction mechanism. The modeled ab initio-TST and DFT–TST rate constants calculated without tunneling were found to be in reasonable agreement with the observed ones indicating that the contribution of the tunneling effect to the reaction was predicted to be unimportant at ambient temperature.     

Ab initio molecular orbital study on the gas phase SN2 reaction F + CH3Cl → CH3F + Cl

Ab-initio molecular orbital (MO) and direct ab initio dynamics calculations have been applied to the gas phase S_N2 reaction F⁻ + CH₃Cl → CH₃F + Cl⁻. Several basis sets were examined in order to select the most convenient and best fitted basis set to that of high-quality calculations. The Hartree–Fock (HF) 3−21+G(d) calculation reasonably represents a potential energy surface calculated at the MP2/6−311++G(2df,2pd) level. A direct ab initio dynamics calculation at the HF/3−21+G(d) level was carried out for the S_N2 reaction. A full dimensional ab initio potential energy surface including all degrees of freedom was used in the dynamics calculation. Total energies and gradients were calculated at each time step. Two initial configurations at time zero were examined in the direct dynamics calculations: one is a near collinear collision, and the other is a side-attack collision. It was found that in the near collinear collision almost all total available energy is partitioned into two modes: the relative translational mode between the products (40%) and the C − F stretching mode (60%). The other internal modes of CH₃F were still in the ground state. The lifetimes of the early- and late-complexes F⁻ … CH₃Cl and FCH₃ … Cl⁻ are significantly short enough to dissociate directly to the products. On the other hand, in the side-attack collision, the relative translation energy was about 20% of total available energy.     

Direct dynamics studies on the hydrogen abstraction reactions CF3O+CH4(CD4)→CF3OH(CF3OD)+CH3(CD3)

The dynamics properties of the hydrogen abstraction reaction CF₃O+CH₄→CF₃OH+CH₃ are studied by dual-level direct dynamics method. Optimization calculations are preformed by B3LYP and MP2 with the 6-311G(d,p) basis set, and the single-point calculations are done at the multi-coefficient correction method based on quadratic configuration interaction with single and double excitations (MC-QCISD) method. The rate constants are evaluated by canonical variational transition-state theory with a small-curvature tunneling correction over a wide range of temperature 200–2000 K. The agreement between theoretical and experimental rate constants is good in the measured temperature range. The calculated results show that the variational effect is small and almost neglected over the whole temperature range, whereas, the tunneling correction plays a role in the lower temperature range. The kinetic isotope effect for the reaction is ‘normal’. The value of k_H/k_D is 2.38 at room temperature and it decreases with the temperature increasing.     

CH2+O2反应的反应机理

The mechanisms of the CH2＋ O2→ H2O＋ CO and CH2＋ O2→ H2＋ CO2 reactions have been studied by performing ab initio CAS(8,8)/6-31G(d,p) calculations, and five intermediates(IMn) and eight transitions(TSn) have been located along the reaction paths. The predicted path for the CH2＋ O2→ H2O＋ CO is: CH2＋ O2→ TS1→ IM1→ TS2→ IM2→ TS3→ IM3→ TS4→ IM4a→ TS5→ H2O＋ CO. For the CH2＋ O2→ H2＋ CO2 reaction, there are two paths: (i) CH2＋ O2→ TS1→ IM1→ TS2→ IM2→ TS3→ IM3→ TS6→ H2＋ CO2 and (ii) CH2＋ O2→ TS1→ IM1→ TS2→ IM2→ TS3→ IM3→ TS4→ IM4a→ TS7→ IM4b→ TS8→ H2＋ CO2, with the latter path more favorable energetically.     

Ab initio molecular orbital study of potential energy surface for the H2NO(B1)→NO(Π)+H2 reaction

The mechanism of the H₂NO(²B₁)→NO(²Π)+H₂ reaction has been examined using ab initio molecular orbital methods. Ground-state and first-excited-state potential surfaces were plotted at the FOCI/cc-pVTZ level of theory as functions of two appropriate internal degrees of freedom. A conical intersection was found on the C_s pathway that is symmetric with respect to the plane perpendicular to the molecular plane of C_2v H₂NO(²B₁). It is therefore considered that trajectories that start from H₂NO(²B₁) towards the product region detour around the conical intersection, pass through the neighborhood of the transition state that is located at the saddle point on the C_s pathway, and finally reach the products, NO(²Π)+H₂. Thus we can explain the mechanism of the H₂NO(²B₁)→NO(²Π)+H₂ reaction, which has remained unclear to date.     

CH3NO2和CH3自由基吸氢反应途径和变分速率常数计算

采用MP2(full)/6-311G(d, p)从头算方法,优化了硝基甲烷和甲基自由基吸氢反应的过渡态结构,经QCISD(T)方法进行能量校正,得出该反应的正逆向反应的活化位垒分别是58.21 kJ•mol－1和67.17 kJ•mol－1.沿IRC分析指出该反应是氢转移协同反应,而且在反应途径上存在一个引导反应进行的振动模式,这一反应模式引导反应进行的区间在反应坐标S的-0.9～1.0(amu)1/2bohr之间;在温度为800～2600 K范围内,运用改进的变分过渡态理论(ICVT),计算了该反应的速率常数,并与实验类比所得的速率常数随温度的变化趋势进行了比较.     

Theoretical kinetic study of two competing paths in the H2O2 + O(P) → HO2 + OH reaction

Ab initio-TST calculations were carried out to study the kinetics of the title reaction. The H atom and the OH abstraction paths leading to the same products HO₂ and OH have been considered. The ZPE and BSSE corrected classical barrier heights were predicted to be 7.4 and 17.3 kcal/mol, respectively. Calculated thermal rate constants over the temperature range 300–5000 K showed that the H-abstraction path was the most likely to occur for temperatures below 2500 K which confirms the result found in a previous study [Y. Tarchouna, M. Bahri, N. Jaïdane, Z. Ben Lakdar, J. Mol. Struct. (Theochem), 189 (2003) 664]. The contribution of OH abstraction path to the reaction was predicted to be important for high temperatures.     

Ab initio and density function theory studies on NLi6 molecule: compared with CLi6 and OLi6 species

A wide variety of geometrical structures of NLi₆ molecule were studied using HF ab initio and BLYP-DFT techniques. Three stationary configurations which take D_4h, D_3d and D_2d configurations were obtained. Their equilibrium geometries and fundamental frequencies were calculated at HF and BLYP-DFT levels. Among the three stable states, the global minimum prefers D_3d configuration, which is different from those of CLi₆ and OLi₆. The D_3d isomer of NLi₆ is 3.43 and 28.45 kcal/mol lower in energy than the D_4h and D_2d ones in the DFT calculation, respectively. All calculations were performed with 6-31G* polarized split-valence basis set.     

Ab initio study of the two competitive channels HO2 + H → H2 + O2 and HO2 + H → H2O + O

Saddle point geometries and barrier heights have been calculated for the H abstraction reaction HO₂(²A″)+H(²S) → H₂(¹Σ⁺_g)+O₂(³Σ⁻_g) and the concerted H approach-O removing reaction HO₂ (²A″)+H(²S) → H₂O(¹A₁)+O(³P) by using SDCI wavefunctions with a valence double-zeta plus polarization basis set. The saddle points are found to be of C_s symmetry and the barrier heights are respectively 5.3 and 19.8 kcal by including size consistent correction. Moreoever kinetic parameters have been evaluated within the framework of the TST theory. So activation energies and the rate constants are estimated to be respectively 2.3 kcal and 0.4×10⁹ ℓ mol⁻¹ s⁻¹ for the first reaction, 20.0 kcal and 5.4.10⁻⁵ ℓ mol⁻¹ s⁻¹ for the second. Comparison of these results with experimental determinations shows that hydrogen abstraction on HO₂ is an efficient mechanism for the formation of H₂ + O₂, while the concerted mechanism envisaged for the formation of H₂O + O is highly unlikely.     

Ab initio and equilibrium bond angles. Structures of HNO and H2O2

The possibility of calculating accurate ab initio bond angles is examined using a sample of 29 molecules (35 independent angles) containing only first row atoms and whose equilibrium structures are known. Three different correlated methods are compared: MP2, CCSD(T), and DFT, using the hybrid functional B3LYP. The convergence of Dunning's correlation consistent polarized valence basis sets, cc-pVnZ is also studied. It is found that the CCSD(T) method is consistently the most accurate; the DFT/B3LYP being slightly less reliable than MP2. It is shown that when convergence of the basis set is achieved (which is dependent on the kind of bonding) and when the effect of diffuse functions on electronegative atoms is taken into account, a high accuracy may be obtained: 0.03° for the median of absolute deviations or 0.07° for the mean absolute deviation. It does not exclude the possibility that the ab initio method may fail in some particular case, for instance when a large amplitude motion is involved. The MP2/cc-pVQZ method gives a mean absolute deviation of 0.22° to be compared with the 0.07° of the CCSD(T) method. To obtain these results, it was necessary to reanalyze the structure of a few molecules, particularly, a new and more accurate structure is proposed for nitroxyl, HNO and hydrogen peroxide, H₂O₂.     

2D and 3D quantum scattering calculations on the CH4 + OH → CH3 + H2O reaction

The effect on the thermal rate constant and the differential cross-sections of varying the dimensionality of quantum scattering calculations of a polyatomic reaction is investigated. The rotating bond approximation (RBA; 3D) and a rotating line approximation (RLA; 2D) are used for the CH₄ + OH → CH₃ + H₂O reaction. It is found that the RBA and RLA results are in close agreement when an adiabatic treatment is used for the degree of freedom which is treated explicitly in the RBA but not in the RLA.     

Ab initio and density functional theory calculations of molecular structure and vibrational spectrum of ethyl azidoacetate

Here we report ab initio and density functional results for molecular properties of ethyl azidoacetate (N₃CH₂COOC₂H₅) and for the corresponding singly ionized structure (N₃CH₂COOC₂H₅⁺). Ab initio ionization energies based on Koopmans’ theorem are in excellent agreement with the experimental data from ultraviolet photoelectron spectroscopy. DFT adiabatic energy differences between neutral and ionized structures are very sensitive to electronic correlation effects and are not in very good agreement with experiment. The results for the structure and vibrational frequencies are compared with the experimental data of related molecular structures.     

On the NH3+HCO3H→HCOOH+H3NO mechanism: a density functional theory study

The B1LYP, B3LYP and MPW1PW91 density functional theory methods combined with the 6-311G(2d, 2p) basis set were used to carry out a density functional theory study of the NH₃+HCO₃H→HCOOH+H₃NO reaction. The purpose of this work is to study the reaction mechanism from the viewpoint of bond order transformations throughout the course of the reaction, and propose the reasons for the apparent differences in activation barriers.     

Ab initio studies of the intramolecular amide hydrolysis in N-methylmaleamic acids

The intramolecular amide hydrolysis reactions of N-methylmaleamic acids (NMMA) are studied at the MP2/6-31G**//RHF/4-31G level of theory as model reactions of peptide bond cleavage by a proteolytic enzyme. In contrast to the previously reported results for a bimolecular reaction model of peptide hydrolysis, the unimolecular reactions studied here proceed via the concerted pathway in which the C–O bond formation and the release of methylamine occur simultaneously in preference to the stepwise one. The determination of an intrinsic reaction coordinate shows that the reaction is facilitated by the intramolecular proton transfer from the undissociated carboxyl group to the nitrogen of the leaving amine group. Mainly because of the increase in activation energy, methyl substitution at the 2-position retards the hydrolysis reaction rate by a factor of 14 compared to the reaction of the unsubstituted molecule. In contrast, additional methyl substitution at 3-position leads to 35-fold increase in the reaction rate. These variations of reactivity are caused by the charge redistributions in the amide group induced by methyl substitutions and the resulting changes in electrophilicity of the aminocarbonyl carbon.     

Theoretical molecular structures for partially bonded complexes of trimethylamine with SO2 and SO3: ab initio and density functional theory calculations

The self-consistent reaction field (SCRF) method based on Onsager's reaction field theory is applied to investigate the effect of polar media on molecular structures of complexes of trimethylamime (TMA) with SO_x (x=2,3). The calculated SCRF N–S bond lengths at the MPW1PW91/6-311+G(3df) level are in satisfactory agreement with the experimental N–S bond lengths for the TMA–SO_x upon crystallization. The results are enough to demonstrate the usefulness of the reaction field theory in providing qualitative understanding of the medium effect on the partially bonded system such as TMA–SO_x.     

High level of ab initio and density functional theory computational study of structural and energetic properties of tetrafluorhydrazine rotamers

The HF, MP2, MP3, MP4, and QCISD ab initio methods were compared with local, hybrid, and gradient-corrected density functional theory (DFT) methods for computing structures and energies of N2F4 rotamers. In all DFT calculations 6-311 + G(2d) basis set was used. The generated structures energies of trans- and gauche-N₂F₄ rotamers, and their dissociation energies to nitrogen difluoride were compared with experimental data. Suitable hybrid and gradient-corrected DFT methods for determining structures and energies for these and similar molecular systems were discussed.     

Conformation of N,N′-bis (2,6-dimethylmorpholino), N″-dichloroacetyl phosphoric triamide: CHCl2C(O)NHP(O)[2,6(CH3)2(NC4H6O)]2. NMR and ab initio studies

Using phosphorus pentachloride as a substrate, a new carbacyclamidophosphate, N,N″-bis (2,6-dimethylmorpholino), N″-dichloroacetyl phosphoric triamide (1) has been synthesized and characterized by ¹H, ³¹P and ¹³C NMR, IR spectroscopy and elemental analysis. Due to the presence of methyl disubstituted morpholine rings and the dichloroacetamide group, several conformers can be considered for this molecule. The ³¹P{¹H} NMR spectra for the isomeric mixture of synthesized compound showed four signals with the ratio 67.1; 19.0; 12.2; 1.7, which indicates four independent conformers. The ¹H NMR spectra confirmed these results. The conformational space and the molecular geometry of the molecule in the gaseous phase have been studied using the B3LYP method of approximation, with 6-31G and 6-311++G^** basis sets.     

HNCO+OH->H2O+NCO的反应机理

采用从头算分子轨道法 (UHF/6 31G 水平 ,并用MP4加以相关能校正 )研究了HNCO OHH2 O NCO反应机理 .同时用Morokuma数值法获得了反应途径即内禀反应坐标 (IRC) .沿着IRC ,运用反应途径哈密顿理论 ,获得反应途径动态学信息 .在此基础上 ,根据过渡态理论和相应隧道效应校正 ,计算了在不同温度下的反应速率常数 ,得到了和实验相一致的结果 .计算结果表明 ,此反应是一步直接型的抽提H反应 .     

Thermochemistry of the higher chlorine oxides ClOx (x=3, 4) and Cl2Ox (x=3–7)†

Heats of formation for ClO₃, ClO₄, Cl₂O₃, Cl₂O₄, Cl₂O₅, Cl₂O₆ and Cl₂O₇ molecules are determined at the B3LYP, B3PW91, mPW1PW91 and B1LYP levels of the density functional theory employing a series of extended basis sets, and using Gaussian-3 model chemistries. Modified Gaussian-3 calculations, which employ accurate B3LYP/6-311+G(3d2f) molecular geometries and vibrational frequencies, were also performed. Heats of formation were calculated from both total atomization energies and isodesmic reaction schemes. The latter method in conjunction with Gaussian-3 models leads to the most reliable results. The best values at 298 K for ClO₃, ClO₄, Cl₂O₃ and Cl₂O₄ as derived from an average of G3//B3LYP and G3//B3LYP/6-311+G(3d2f) calculations are 43.1, 54.8, 31.7 and 37.4 kcal mol⁻¹. From calculations carried out at the G3(MP2)//B3LYP and G3(MP2)//B3LYP/6-311+G(3d2f) levels, heats of formation for Cl₂O₅, Cl₂O₆ and Cl₂O₇ are predicted to be 53.2, 52.2 and 61.5 kcal mol⁻¹. All best values are reproduced within 1 kcal mol⁻¹ by using mPW1PW91/6-311+G(3d2f) isodesmic energies. Enthalpy changes for relevant Cl–O bond fission reactions are reported. Comparisons with previous thermodynamics data are made.     

An ab initio study of the Cl(P)+C2H6→C2H5+HCl abstraction reaction

Electronic energies, geometries, and harmonic vibration frequencies for the reactants, products, and transition state for the Cl(³P)+C₂H₆→C₂H₅+HCl abstraction reaction were evaluated at the HF and MP2 levels using several correlation consistent polarized-valence basis sets. Single-point calculations at PMP2, MP4, QCISD(T), and CCSD(T) levels were also carried out. The values of the forward activation energies obtained at the MP4/cc-pVTZ, QCISD(T)/cc-pVTZ, and CCSD(T)/cc-pVTZ levels using the MP2/cc-pVTZ structures are equal to −0.1, −0.4, and −0.3 kcal/mol, respectively. The experimental value is equal to 0.3±0.2 kcal/mol. We found that the MP2/aug-cc-pVTZ adiabatic vibration energy for the reaction (−2.4 kcal/mol) agrees well with the experimental value −(2.2–2.6) kcal/mol. Rate constants calculated with the zeroth-order interpolated variational transition state (IVTST-0) method are in good agreement with experiment. In general, the theoretical rate constants differ from experiment by, at most, a factor of 2.6.