Ab initio chemical kinetics for the NH2 + HNOx reactions,part II: Kinetics and mechanism for NH2 + HONO

The kinetics and mechanism for the reaction of NH₂ with HONO have been investigated by ab initio calculations with rate constant prediction. The potential energy surface of this reaction has been computed by single‐point calculations at the CCSD(T)/6‐311+G(3df, 2p) level based on geometries optimized at the CCSD/6‐311++G(d, p) level. The reaction producing the primary products, NH₃ + NO₂, takes place via precomplexes, H₂N???c‐HONO or H₂N???t‐HONO with binding energies, 5.0 or 5.9 kcal/mol, respectively. The rate constants for the major reaction channels in the temperature range of 300–3000 K are predicted by variational transition state theory or Rice–Ramsperger–Kassel–Marcus theory depending on the mechanism involved. The total rate constant can be represented by k_total = 1.69 × 10^?20 × T^2.34 exp(1612/T) cm³ molecule^?1 s^?1 at T = 300–650 K and 8.04 × 10^?22 × T^3.36 exp(2303/T) cm³ molecule^?1 s^?1 at T = 650–3000 K. The branching ratios of the major channels are predicted: k₁ + k₃ producing NH₃ + NO₂ accounts for 1.00–0.98 in the temperature range 300–3000 K and k₂ producing OH + H₂NNO accounts for 0.02 at T > 2500 K. The predicted rate constant for the reverse reaction, NH₃ + NO₂ → NH₂ + HONO represented by 8.00 × 10^?26 × T^4.25 exp(?11,560/T) cm³ molecule^?1 s^?1, is in good agreement with the experimental data. © 2009 Wiley Periodicals, Inc. Int J Chem Kinet 41: 678–688, 2009     

A quantum mechanical approach to the kinetics of the hydrogen abstraction reaction H2O2 + •OH → HO2 + H2O

The kinetics of the hydrogen abstraction from H₂O₂ by ^?OH has been modeled with MP2/6‐31G*//MP2/6‐31G*, MP2‐SAC//MP2/6‐31G*, MP2/6‐31+G**//MP2/6‐31+G**, MP2‐SAC// MP2/6‐31+G**, MP4(SDTQ)/6‐311G**//MP2/6‐31G*, CCSD(T)/6‐31G*//CCSD(T)/6‐31G*, CCSD(T)/6‐31G**//CCSD(T)/6‐31G**, CCSD(T)/6‐311++G**//MP2/6‐31G* in the gas phase. MD simulations have been used to generate initial geometries for the stationary points along the potential energy surface for hydrogen abstraction from H₂O₂. The effective fragment potential (EFP) has been used to optimize the relevant structures in solution. Furthermore, the IEFPCM model has been used for the supermolecules generated via MD calculations. IEFPCM/MP2/6‐31G* and IEFPCM/CCSD(T)/6‐31G* calculations have also been performed for structures without explicit water molecules. Experimentally, the rate constant for hydrogen abstraction by ^?OH drops from 1.75 × 10^?12 cm³ molecule^?1 s^?1 in the gas phase to 4.48 × 10^?14 cm³ molecule^?1 s^?1 in solution. The same trend has been reproduced best with MP4 (SDTQ)/6‐311G**//MP2/6‐31G* in the gas phase (0.415 × 10^?12 cm³ molecule^?1 s^?1) and with EFP (UHF/6‐31G*) in solution (3.23 × 10^?14 cm³ molecule^?1 s^?1). © 2005 Wiley Periodicals, Inc. Int J Chem Kinet 37: 502–514, 2005     

Formic acid catalyzed gas-phase reaction of H2O with SO3 and the reverse reaction: a theoretical study

The formic acid catalyzed gas‐phase reaction between H₂O and SO₃ and its reverse reaction are respectively investigated by means of quantum chemical calculations at the CCSD(T)//B3LYP/cc‐pv(T+d)z and CCSD(T)//MP2/aug‐cc‐pv(T+d)z levels of theory. Remarkably, the activation energy relative to the reactants for the reaction of H₂O with SO₃ is lowered through formic acid catalysis from 15.97 kcal mol^?1 to ?15.12 and ?14.83 kcal mol^?1 for the formed H₂O ??? SO₃ complex plus HCOOH and the formed H₂O ??? HCOOH complex plus SO₃, respectively, at the CCSD(T)//MP2/aug‐cc‐pv(T+d)z level. For the reverse reaction, the energy barrier for decomposition of sulfuric acid is reduced to ?3.07 kcal mol^?1 from 35.82 kcal mol^?1 with the aid of formic acid. The results show that formic acid plays a strong catalytic role in facilitating the formation and decomposition of sulfuric acid. The rate constant of the SO₃+H₂O reaction with formic acid is 10⁵ times greater than that of the corresponding reaction with water dimer. The calculated rate constant for the HCOOH+H₂SO₄ reaction is about 10^?13 cm³ molecule^?1 s^?1 in the temperature range 200–280 K. The results of the present investigation show that formic acid plays a crucial role in the cycle between SO₃ and H₂SO₄ in atmospheric chemistry.     

Ab initio chemical kinetics for the NH2 + HNOx Reactions,Part I: Kinetics and Mechanism for NH2 + HNO

The kinetics and mechanism for the reaction of NH₂ with HNO have been investigated by ab initio calculations with rate constant prediction. The potential energy surface of this reaction has been computed by single‐point calculations at the CCSD(T)/6‐311+G(3df, 2p) level based on geometries optimized at the CCSD/6‐311++G(d, p) level. The major products of this reaction were found to be NH₃ + NO formed by H‐abstraction via a long‐lived H₂N???HNO complex and the H₂NN(H)O radical intermediate formed by association with 26.9 kcal/mol binding energy. The rate constants for formation of primary products in the temperature range of 300–3000 K were predicted by variational transition state or RRKM theories. The predicted total rate constants at the 760 Torr Ar pressure can be represented by k_total = 3.83 × 10^?20 × T^+2.47exp(1450/T) at T = 300–600 K; 2.58 × 10^?22 × T^+3.15 exp(1831/T) cm³ molecule^?1 s^?1 at T = 600?3000 K. The branching ratios of major channels at 760 Torr Ar pressure are predicted: k₁ + k₃ + k₄ producing NH₃ + NO accounts for 0.59–0.90 at T = 300–3000 K peaking around 1000 K, k₂ accounts for 0.41–0.03 at T = 300–600 K decreasing with temperature, and k₅ accounts for 0.07–0.27 at T > 600 K increasing gradually with temperature. The NH₃ + NO formation rate constant was found to be a factor of 3–10 smaller than that of the isoelectronic reaction CH₃ + HNO producing CH₄ + NO, which has been shown to take place by barrierless H‐abstraction without involving a hydrogen‐bonding complex as in the NH₂ case. © 2009 Wiley Periodicals, Inc. Int J Chem Kinet 41: 677–677, 2009     

Ab initio chemical kinetics for the NH2 + HNOx reactions,part III: Kinetics and mechanism for NH2 + HONO2

The kinetics and mechanism for the reaction of NH₂ with HONO₂ have been investigated by ab initio calculations with rate constant prediction. The potential energy surface of this reaction has been computed by single‐point calculations at the CCSD(T)/6‐311+G(3df, 2p) level based on geometries optimized at the B3LYP/6‐311+G(3df, 2p) level. The reaction producing the primary products, NH₃ + NO₃, takes place via a precursor complex, H₂N…HONO₂ with an 8.4‐kcal/mol binding energy. The rate constants for major product channels in the temperature range 200–3000 K are predicted by variational transition state or variational Rice–Ramsperger–Kassel–Marcus theory. The results show that the reaction has a noticeable pressure dependence at T < 900 K. The total rate constants at 760 Torr Ar‐pressure can be represented by k_total = 1.71 × 10^?3 × T^?3.85 exp(?96/T) cm³ molecule^?1 s^?1 at T = 200–550 K, 5.11 × 10^?23 × T^+3.22 exp(70/T) cm³ molecule^?1 s^?1 at T = 550–3000 K. The branching ratios of primary channels at 760 Torr Ar‐pressure are predicted: k₁ producing NH₃ + NO₃ accounts for 1.00–0.99 in the temperature range of 200–3000 K and k₂ + k₃ producing H₂NO + HONO accounts for less than 0.01 when temperature is more than 2600 K. The reverse reaction, NH₃ + NO₃ → NH₂ + HONO₂ shows relatively weak pressure dependence at P < 100 Torr and T < 600 K due to its precursor complex, NH₃…O₃N with a lower binding energy of 1.8 kcal/mol. The predicted rate constants can be represented by k_?1 = 6.70 × 10^?24 × T^+3.58 exp(?850/T) cm³ molecule^?1 s^?1 at T = 200–3000 K and 760 Torr N₂ pressure, where the predicted rate at T = 298 K, 2.8 × 10^?16 cm³ molecule^?1 s^?1 is in good agreement with the experimental data. The NH₃ + NO₃ formation rate constant was found to be a factor of 4 smaller than that of the reaction OH + HONO₂ producing the H₂O + NO₃ because of the lower barrier for the transition state for the OH + HONO₂. © 2009 Wiley Periodicals, Inc. Int J Chem Kinet 42: 69–78, 2010     

Ab initio study of the OH + CH2O reaction: The effect of the OH··OCH2 complex on the H‐abstraction kinetics

Kinetics for the reaction of OH radical with CH₂O has been studied by single‐point calculations at the CCSD(T)/6‐311+G(3df, 2p) level based on the geometries optimized at the B3LYP/6‐311+G(3df, 2p) and CCSD/6‐311++G(d,p) levels. The rate constant for the reaction has been computed in the temperature range 200–3000 K by variational transition state theory including the significant effect of the multiple reflections above the OH··OCH₂ complex. The predicted results can be represented by the expressions k₁ = 2.45 × 10^‐21 T^2.98 exp (1750/T) cm³ mol^?1 s^?1 (200–400 K) and 3.22 × 10^‐18 T^2.11 exp(849/T) cm³ mol^?1 s^?1 (400–3000 K) for the H‐abstraction process and k₂ = 1.05 × 10^‐17 T^1.63 exp(?2156/T) cm³ mol^?1 s^?1 in the temperature range of 200–3000 K for the HO‐addition process producing the OCH₂OH radical. The predicted total rate constants (k₁ + k₂) can reproduce closely the recommended kinetic data for OH + CH₂O over the entire range of temperature studied. © 2006 Wiley Periodicals, Inc. Int J Chem Kinet 38: 322–326, 2006     

Direct dynamics study on the mechanism and the kinetics of the reaction of CH3NH2 with OH

A theoretical study of the mechanism and the kinetics for the hydrogen abstraction reaction of methylamine by OH radical has been presented at the CCSD(T)/6‐311 ++G(2d,2p)//CCSD/6‐31G(d) level of theory. Our theoretical calculations suggest a stepwise mechanism involving the formation of a prereactant complex in the entrance channel and a preproduct complex in the exit channel, for the two hydrogen abstraction channels involving the methyl and amine groups. For clarity, the diagram of potential for the reaction is given. The calculated standard reaction enthalpies are ?98.48 and ?76.50 kJ mol^?1 and barrier heights are 0.36 and 25.25 kJ mol^?1, respectively. The rate constants are evaluated by means of the improved canonical variational transition state theory with small‐curvature tunneling correction (ICVT/SCT) in the temperature range of 299–3000 K. The calculated results show that the rate constants at experimentally measured temperatures are in good agreement with the experimental values. It is shown that the calculated rate constants exhibit a non‐Arrhenius behavior. Moreover, the variational effect is obvious in the calculated temperature range. The dominant product channel is to form CH₂NH₂ and H₂O via hydrogen abstraction from the CH₃ group of CH₃NH₂ by OH in the calculated temperature range. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

A Theoretical Study of the Kinetics of the Hydrogen Atom Abstraction Reactions from Cyclopropane by H,O (3P), and Cl (2P3/2) Atoms and OH Radicals

The rate constants of the H‐abstraction reactions from cyclopropane by H, O (³P), Cl (²P_3/2), and OH radicals have been calculated over the temperature range of 250?2500 K using two different levels of theory. Calculations of optimized geometrical parameters and vibrational frequencies are performed using the MP2 method combined with the cc‐pVTZ basis set and the 6–311++G(d,p) basis set. Single‐point energy calculations have been carried out with the highly correlated ab initio coupled cluster method in the space of single, double, and triple (perturbatively) electron excitations CCSD(T) using either the cc‐pVTZ, aug‐cc‐pVTZ, and aug‐cc‐pVQZ basis sets or the 6–311++G(3df,3pd) basis set. The CCSD(T) calculated potential energies have been extrapolated to the complete basis limit (CBS) limit. The Full Configuration Interaction (FCI) energies have been also estimated using the continued‐fraction approximation as proposed by Goodson (J. Chem. Phys., 2002, 116, 6948–6956). Canonical transition‐state theory combined with an Eckart tunneling correction has been used to predict the rate constants as a function of temperature using two kinetic models (direct abstraction or complex mechanism) at two levels of theory (CCSD(T)‐cf/CBS//MP2/cc‐pVTZ and CCSD(T)‐cf/6–311++G(3df,3pd)//MP2/6–311++G(d,p)). The calculated kinetic parameters are in reasonable agreement with their literature counterparts for all reactions. In the light of these trends, the use of the Pople‐style basis sets for studying the reactivity of other systems such as larger cycloalkanes or halogenated cycloalkanes is recommended because the 6–311++G(3df,3pd) basis set is less time consuming than the aug‐cc‐pVQZ basis set. Based on our calculations performed at the CCSD(T)‐cf/CBS//MP2/cc‐pVTZ level of theory, the standard enthalpy of formation at 298 K for the cyclopropyl radical has been reassessed and its value is (290.5 ± 1.6) kJ mol^?1.     

A theoretical study of the thermodynamics and kinetics of small organosulfur compounds

The performance of a large set of ab initio procedures in predicting geometries, thermochemical and kinetic data of small sulfur compounds is assessed. Geometries and thermochemical data for H₂S, (CH₃)₂S, H₂S₂, (CH₃)₂S₂ and H₂C=S are studied using the HF method, density functional theory methods (B3LYP, BHandHLYP, MPW1PW91 and BMK), post-HF methods [MP2, MP3, MP4, CCSD, CCSD(T) and QCISD] and composite techniques (G3, G3B3, CBS-QB3 and W1U). A set of five reactions involving these small organosulfur compounds is studied and the influence of the level of theory on transition state geometries, reaction barriers and rate coefficients is assessed. Independent of the level of theory used, accurate geometries are obtained with the 6-311G(2d,d,p) and cc-pVTZ basis sets, both reproducing experimental bond lengths and bond angles within 2 pm and 0.5°. Besides composite methods, the BMK/cc-pVTZ method is the only studied method that succeeds to predict standard enthalpies of formation within 10 kJ mol^?1 of the experimental data. The best agreement with experimental rate coefficients is obtained with the BHandHLYP/cc-pVTZ method, closely followed by the composite methods and the BMK/cc-pVTZ method. All these methods succeed to reproduce the experimental rate coefficients within a factor 4. To obtain an accurate prediction of both thermochemical and kinetic data for organosulfur compounds, the commonly used composite methods G3B3 and CBS-QB3 and the BMK/cc-pVTZ method prove to be valuable tools.     

Computational study on the mechanisms and kinetics of the CH<Subscript>2</Subscript>=CHCH<Subscript>2</Subscript>I with OH reaction

The potential energy surface for the reaction of OH with CH₂═CHCH₂I has been studied at the CCSD(T)//M06-2X/6-311++G(d,p) level of theory. Three different reaction entrances were revealed, namely, terminal-C addition, central-C addition, and H-abstraction, leading to CH₂OHCHCH₂I (IM1), CH₂CHOHCH₂I (IM2), and H₂O?+?C₃H₄I, respectively. Several conceivable decomposition and isomerization channels were also examined for IM1 and IM2. The total and individual rate constants were calculated by using multichannel RRKM and TST theories over a wide range of temperatures (200–3000 K) and pressures(10^?14–10¹⁴ Torr).     

Mechanism and rate constants for ammonia photochemical oxidation with O2 on the singlet and triplet potential energy surfaces

The mechanism and thermodynamic of NH₃ + O₂ reaction on the singlet and triplet potential energy surfaces (PES), were carried out using the RMP2 and CCSD (T)//RMP2 theoretical approaches in connection with the 6-311++G(d, p) basis set. Three pre-reactive complexes, ¹C1, ¹C2, and ³C1 on the singlet and triplet PES were formed between ammonia and molecular oxygen. With variety of pre-reactive complexes, six types of products are obtained, of which two types are found to be thermodynamically stable. The mechanistic properties of all products channels are discussed. Results show that production of HONO + H₂ and HN(OH)₂ are the main reaction channels in thermodynamic viewpoint with the Gibbs free energy of ?G° = ?34.681 and ?27.153 kcal/mol, respectively. Rate constants of the title reaction over the temperature range of (200–1000 K) show kinetic products are different from thermodynamic products.     

Theoretical study on the kinetics and branching ratios of the gas phase reactions of 1, 1-Dichlorodimethylether (DCDME) with Cl atom

Quantum mechanical calculations are carried out on the reactions of CH₃OCHCl₂ (DCDME) with Cl atom by means of DFT and couple cluster methods. The geometries of the reactants, products, and transition states involved in the reaction pathways are optimized at BHandHLYP level of theory using 6-311G(d,p) basis set. Transition states are searched on the potential energy surface involved during the reaction channels, and each of the transition states is characterized by the presence of only one imaginary frequency. The existence of transition states on the corresponding potential energy surface is ascertained by performing intrinsic reaction coordinate calculation. Single point energy calculations are performed at CCSD(T) level using the same basis set. The hydrogen abstraction rate constant for the title reaction is calculated at 298 K and atmospheric pressure using the canonical transition state theory including tunneling correction. The calculated value for rate constant as 1.204 × 10^?12 cm³ molecule^?1 s^?1 is found to be in very good agreement with the recent experimental data. The percentage contributions of both reaction channels are also reported at 298 K.     

Theoretical Studies of the Reactions CFxH3−xCOOR+Cl and CF3COOCH3+OH

The mechanism and kinetics of the reactions of CF₃COOCH₂CH₃, CF₂HCOOCH₃, and CF₃COOCH₃ with Cl and OH radicals are studied using the B3LYP, MP2, BHandHLYP, and M06‐2X methods with the 6‐311G(d,p) basis set. The study is further refined by using the CCSD(T) and QCISD(T)/6‐311++G(d,p) methods. Seven hydrogen‐abstraction channels are found. All the rate constants, computed by a dual‐level direct method with a small‐curvature tunneling correction, are in good agreement with the experimental data. The tunneling effect is found to be important for the calculated rate constants in the low‐temperature range. For the reaction of CF₃COOCH₂CH₃+Cl, H‐abstraction from the CH₂ group is found to be the dominant reaction channel. The standard enthalpies of formation for the species are also calculated. The Arrhenius expressions are fitted within 200–1000 K as k_T(1)=8.4×10^?20T ^2.63exp(381.28/T), k_T(2)=2.95×10^?21T ^3.13exp(?103.21/T), k_T(3)=1.25×10^?23T ^3.37exp(791.98/T), and k_T(4)=4.53×10^?22T ^3.07exp(465.00/T).     

Accurate binding energies of hydrogen,halogen, and dihydrogen bonded complexes and cation enhanced binding strengths

Interaction energy (E_int) values of a variety of hydrogen, halogen, and dihydrogen bonded complexes in the weak, medium, and strong regimes have been computed using W1BD, MP2, M06L density functional theory, and hybrid methods MP4//MP2, MP4//M06L, and CCSD(T)//MP2. W1BD level E_int and CCSD(T) results reported in the literature show very good agreement (mean absolute deviation = 0.19 kcal/mol). MP2 underestimates E_int while M06L shows accurate behavior for all except halogen and charge‐assisted hydrogen bonds. MP4//MP2, MP4//M06L, and CCSD(T)//MP2 yield E_int very close to those obtained from W1BD. The high accuracy energy data at MP4/MP2 is used to study the effect of a cation (Li⁺, NH₄⁺) on the E_int. The cation enhances electron donation from the donor to noncovalent bonding region leading to substantial enhancement in E_int (～141–566% for Li⁺ and ～105–539% for NH₄⁺) and promotes a noncovalent bond in the weak regime to medium regime and that in the medium regime to strong regime. © 2014 Wiley Periodicals, Inc.     

Conformational properties of 1-methyl-1-germacyclohexane: low-temperature NMR and quantum chemical calculations

The conformational preference of the methyl group of 1-methyl-1-germacyclohexane was studied experimentally in solution (low-temperature ¹³C NMR) and by quantum chemical calculations (CCSD(T), MP2 and DFT methods). The NMR experiment resulted in an axial/equatorial ratio of 44/56 mol% at 114 K corresponding to an A value (A = G _ax ? G _eq) of 0.06 kcal mol^?1. An average value for ΔG _e→a ^#  = 5.0 ± 0.1 kcal mol^?1 was obtained for the temperature range 106–134 K. The experimental results are very well reproduced by the calculations. CCSD(T)/CBS calculations + thermal corrections resulted in an A value of 0.02 kcal mol^?1, whereas a ΔE value of ?0.01 kcal mol^?1 at 0 K was obtained.     

Direct dynamics study on mechanism and kinetics of the biradical self‐reaction of HOO

A theoretical study of the mechanism and the kinetics for the hydrogen abstraction reaction of the biradical hydroperoxy radical has been presented at the CCSD(T)/6‐311++G(3d,2p)//CCSD/6‐31+G(d,p) level of theory. Our theoretical calculations suppose a stepwise mechanism involving the formation of a postreactant complex in the triplet and singlet entrance channels. Four transition states of the six‐membered chain complexes (³TS1 and ¹TS1) and six‐membered ring complexes (³TS2 and ¹TS2) are located at the high dual level CCSD(T)/6‐311++G(3d,2p)//CCSD/6‐31+G(d,p) method. The rate constants of Path 1 ～ Path 4 at the CCSD(T)/6‐311++G(3d,2p)//CCSD/6‐31+G (d,p) level are calculated by means of the conventional transition state theory (TST) and canonical variational TST without and with small‐curvature tunneling (SCT) correction within the temperature range of 200–2,500 K. The calculated results show that the triplet channel is the dominating reaction channel and Path 2 is found to be the most favorable pathway. The rate constants of Path 2 are in good agreement with the experimental values at the experimentally measured temperatures. Moreover, the variational effect is not obvious in the low temperature range but is not neglectable in the high temperature range. The SCT plays an important role particularly in the low temperature range. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Computational mechanistic study of methanol and molecular oxygen reaction on the triplet and singlet potential energy surfaces

The reaction of CH₃OH with the O₂ on the triplet and singlet potential energy surfaces (PES) was carried out using the B3LYP, MP2, and CCSD(T)//B3LYP theoretical approaches in connection with the 6-311++G(3df–3pd) basis set. Three pre-reactive complexes, ¹C1, ¹C2, and ³C1, on the singlet and triplet PES were formed between methanol and molecular oxygen. From a variety of the complexes, seven types of products are obtained, of which four types are found to be thermodynamically stable. Results reveal that there exists one intersystem crossing between triplet and singlet PES. For P4 adduct that is the main and kinetically the most favorable product, the rate constants are calculated in the temperature range of 200–1,000 K in the reliable pathway.     

An ab initio chemical kinetic study on the reactions of H,OH, and Cl with HOClO3

The mechanisms for reactions of H, HO, and Cl with HOClO₃, important elementary processes in the early stages of the ammonium perchlorate (AP) combustion reaction, have been investigated at the CCSD(T)/6‐311+G(3df,2p)//PW91PW91/6‐311+G(3df) level of theory. The rate constants for the low‐energy channels have been calculated by statistical theory. For the reaction of H and HOClO₃, the main channels are the production of H₂ + ClO₄ (k_1a) and HO + HOClO₂ (k_1b); k_1a and k_1b can be represented as 1.07 × 10^?17 T^1.97 exp(?7484/T) and 6.08 × 10^?17T^1.96 exp(?7729/T) cm³ molecule^?1 s^?1, respectively. For the HO + HOClO₃ reaction, the main pathway is the H₂O + ClO₄ (k_2a) production process, with the predicted rate constant k_2a = 1.24 × 10 ^?8 T^?2.99 exp(1664/T) for 300–500 K and k_2a = 1.27×10^?19 T^2.12 exp(?1474/T) for 500–3000 K. For the Cl + HOClO₃ reaction, the formation of HOCl + ClO₃ (k_3a) and HCl + ClO₄ (k_3b) is dominant, with k_3a = 1.33×10^?12 T^0.67 exp(?9658/T) and k_3b = 1.75×10¹⁶ T^1.63 exp(?11156/T) cm³ molecules^?1 in the range of 300–3000 K. In addition, the heats of formation of ClO₃ and HOClO₃ have been predicted based on several isodesmic and/or isogyric reactions with Δ_fH^o₀ (ClO₃) = 47.0 ± 1.0 and Δ_fH^o₀ (HOClO₃) = 5.5 ± 1.5 kcal/mol, respectively. These data may be used for kinetic simulation of the AP decomposition and combustion reaction. © 2010 Wiley Periodicals, Inc. Int J Chem Kinet 42: 253–261, 2010     

CCl2与CH2O插入反应机理及热力学与动力学特性的理论研究

采用密度泛函B3LYP/6-311G*和高级电子相关耦合簇[CCSD(T)/6-311G*]方法计算研究了CCl₂与CH₂O的插入反应机理, 全参数优化了反应势能面各驻点的几何构型, 用内禀反应坐标(IRC)和频率分析方法, 对过渡态进行了验证. 研究结果表明: 反应(1)是单重态二氯卡宾与甲醛插入反应的主反应通道. 该反应由两步组成: (i)两反应物首先经一无能垒的放热反应, 放出9.73 kJ•mol^－1的热量, 生成一中间体IM1, (ii)中间体IM1经一过渡态TS1, 发生H的转移, 生成产物P1, 其势垒为47.32 kJ•mol^－1. 用RRKM-TST理论计算了300～1900 K温度范围内反应(1)的压力效应. 用经Wigner校正的Eyring过渡态理论研究了不同温度下该反应的热力学和动力学性质. 从热力学和动力学角度综合分析, 在高压限101325 Pa下, 该反应进行的适宜温度范围为400～1800 K, 如此, 反应既有较大的自发趋势和平衡常数, 又具有较快的反应速率.     

Features of the potential energy surface for the reaction of HO2 radical with acetone

The mechanism of the reaction of acetone with HO₂ has been studied by quantum chemical computations. Different stationary points on the potential energy surface (PES) of the reaction have been characterized. These stationary points are the reactants, products, molecular complexes, and transition states. Three pathways have been studied: two H‐abstraction channels and one HO₂‐addition channel. The MP2 level of theory with the 6‐311G(d,p) basis set was employed for geometry optimization. The electronic energies was obtained at the PMP2, PMP4, and CCSD(T) level of theory with the 6‐311G(d,p) basis set on the computed geometries. The addition pathway is clearly the more favorable, contrary to the acetone + OH system. The pre‐reactive hydrogen‐bonded complexes have been characterized and show a large red shift between the O? H stretching frequency in the HO₂ radical and the one in the HO₂ fragment of intermolecular complexes. Our addition rate constant k₊ at T = 298 K (3.49 × 10^?16 cm³ s^?1) is consistent with previous experimental results (giving an upper limit of the rate constant of 6 × 10^?16 cm³ s^?1 at 298 K). © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007