Direct ab initio dynamics calculation of the reaction rates of CH3OCl with OH

A direct ab initio dynamics method was carried out for the reaction CH3OCl + OH --> products. Three abstraction channels from chlorine atom, in-plane hydrogen, and out-of-plane hydrogen atoms at the CH3 group have been found. The optimized geometries and frequencies of the stationary points and the minimum-energy paths (MEPs) were calculated at the MP2/6-311G(d,p) level. To improve the reaction enthalpy and potential barrier, single-point calculations were made at three higher levels of theory, the approximate QCISD(T)/6-311++G(3df,2pd), G3, and G3(MP2) levels. Furthermore, the rate constants for three abstraction channels were evaluated using canonical variational transition state theory (CVT) with the small-curvature tunneling correction (SCT) over a wide temperature range of 220-2000 K at above three higher theory levels, respectively. The calculated rate constants as well as branching rates are in reasonable agreement with the experimental values in the temperature region 250-341 K. The present results indicate H-abstraction especially from out-of-plane hydrogen is the main reaction pathway, while Cl-abstraction is much less competitive.     

Direct dynamics study on the hydrogen abstraction reaction CH2O + HO2 --> CHO + H2O2

We present a direct ab initio dynamics study on the hydrogen abstraction reaction CH2O + HO2 --> CHO + H2O2, which is predicted to have four possible reaction channels caused by different attacking orientations of HO2 radical to CH2O. The structures and frequencies at the stationary points and the points along the minimum energy paths (MEPs) of the four reaction channels are calculated at the B3LYP/cc-pVTZ level of theory. Energetic information of stationary points and the points along the MEPs is further refined by means of some single-point multilevel energy calculations (HL). The rate constants of these channels are calculated using the improved canonical variational transition-state theory with the small-curvature tunneling correction (ICVT/SCT) method. The calculated results show that, in the whole temperature range, the more favorable reaction channels are Channels 1 and 3. The total ICVT/SCT rate constants of the four channels at the HL//B3LYP/cc-pVTZ level of theory are in good agreement with the available experiment data over the measured temperature ranges, and the corresponding three-parameter expression is k(ICVT/SCT) = 3.13 x 10(-20) T(2.70) exp(-11.52/RT) cm3 mole(-1) s(-1) in the temperature range of 250-3000 K. Additionally, the flexibility of the dihedral angle of H2O2 is also discussed to explain the different experimental values.     

Thermal decomposition of ethanol. 4. Ab initio chemical kinetics for reactions of H atoms with CH3CH2O and CH3CHOH radicals

The potential energy surfaces of H-atom reactions with CH(3)CH(2)O and CH(3)CHOH, two major radicals in the decomposition and oxidation of ethanol, have been studied at the CCSD(T)/6-311+G(3df,2p) level of theory with geometric optimization carried out at the BH&HLYP/6-311+G(3df,2p) level. The direct hydrogen abstraction channels and the indirect association/decomposition channels from the chemically activated ethanol molecule have been considered for both reactions. The rate constants for both reactions have been calculated at 100-3000 K and 10(-4) Torr to 10(3) atm Ar pressure by microcanonical VTST/RRKM theory with master equation solution for all accessible product channels. The results show that the major product channel of the CH(3)CH(2)O + H reaction is CH(3) + CH(2)OH under atmospheric pressure conditions. Only at high pressure and low temperature, the rate constant for CH(3)CH(2)OH formation by collisonal deactivation becomes dominant. For CH(3)CHOH + H, there are three major product channels; at high temperatures, CH(3)+CH(2)OH production predominates at low pressures (P < 100 Torr), while the formation of CH(3)CH(2)OH by collisional deactivation becomes competitive at high pressures and low temperatures (T < 500 K). At high temperatures, the direct hydrogen abstraction reaction producing CH(2)CHOH + H(2) becomes dominant. Rate constants for all accessible product channels in both systems have been predicted and tabulated for modeling applications. The predicted value for CH(3)CHOH + H at 295 K and 1 Torr pressure agrees closely with available experimental data. For practical modeling applications, the rate constants for the thermal unimolecular decomposition of ethanol giving key accessible products have been predicted; those for the two major product channels taking place by dehydration and C-C breaking agree closely with available literature data.     

C2H3与CH3F氢抽提反应机理与动力学性质

采用双水平直接动力学方法对C₂H₃与CH₃F氢抽提反应进行了研究. 在QCISD(T)/6-311++G(d, p)//B3LYP/6-311G(d, p)水平上, 计算的三个反应通道R1、R2和R3的能垒(ΔE^≠)分别为43.2、43.9和44.1 kJ·mol^-1, 反应热为-38.2 kJ·mol^-1. 此外, 利用传统过渡态理论(TST)、正则变分过渡态理论(CVT)和包含小曲率隧道效应(SCT)的CVT, 分别计算了200-3000 K温度范围内反应的速率常数k^TST、k^CVT和k^CVT/SCT. 结果表明: (1) 三个氢抽提反应通道的速率常数随温度的增加而增大, 其中变分效应的影响可以忽略, 隧道效应则在低温段影响显著; (2) R1反应是主反应通道, 但随着温度的升高, R2反应的竞争力增大, 而R3反应对总速率常数的影响很小.     

A systematic computational study of the reactions of HO2 with RO2: the HO2 + CH2ClO2, CHCl2O2, and CCl3O2 reactions

Potential energy surfaces for the reactions of HO(2) with CH(2)ClO(2), CHCl(2)O(2), and CCl(3)O(2) have been calculated using coupled cluster theory and density functional theory (B3LYP). It is revealed that all the reactions take place on both singlet and triplet surfaces. Potential wells exist in the entrance channels for both surfaces. The reaction mechanism on the triplet surface is simple, including hydrogen abstraction and S(N)2-type displacement. The reaction mechanism on the singlet surface is more complicated. Interestingly, the corresponding transition states prefer to be 4-, 5-, or 7-member-ring structures. For the HO(2) + CH(2)ClO(2) reaction, there are two major product channels, viz., the formation of CH(2)ClOOH + O(2) via hydrogen abstraction on the triplet surface and the formation of CHClO + OH + HO(2) via a 5-member-ring transition state. Meanwhile, two O(3)-forming channels, namely, CH(2)O + HCl + O(3) and CH(2)ClOH + O(3) might be competitive at elevated temperatures. The HO(2) + CHCl(2)O(2) reaction has a mechanism similar to that of the HO(2) + CH(2)ClO(2) reaction. For the HO(2) + CCl(3)O(2) reaction, the formation of CCl(3)O(2)H + O(2) is the dominant channel. The Cl-substitution effect on the geometries, barriers, and heats of reaction is discussed. In addition, the unimolecular decomposition of the excited ROOH (e.g., CH(2)ClOOH, CHCl(2)OOH, and CCl(3)OOH) molecules has been investigated. The implication of the present mechanisms in atmospheric chemistry is discussed in comparison with the experimental measurements.     

Pyrolysis mechanisms of thiophene and methylthiophene in asphaltenes

The pyrolysis mechanisms of thiophene in asphaltenes have been investigated theoretically using density functional and ab initio quantum chemical techniques. All of the possible reaction pathways were explored using B3LYP, MP2, and CBS-QB3 models. A comparison of the calculated heats of reaction with the available experimental values indicates that the CBS-QB3 level of theory is quantitatively reliable for calculating the energetic reaction paths of the title reactions. The pyrolysis process is initiated via four different types of hydrogen migrations. According to the reaction barrier heights, the dominant 1,2-H shift mechanism involves two competitive product channels, namely, C(2)H(2) + CH(2)CS and CS + CH(3)CCH. The minor channels include the formation of CS + CH(2)CCH(2), H(2)S + C(4)H(2), HCS + CH(2)CCH, CS + CH(2)CHCH, H + C(4)H(3)S, and HS + C(4)H(3). The methyl substitution effect was investigated with the pyrolysis of 2-methylthiophene and 3-methylthiophene. The energetics of such systems were very similar to that for unsubstituted thiophene, suggesting that thiophene alkylation may not play a significant role in the pyrolysis of asphaltene compounds.     

CH3OCF2CF2OCH3+Cl的反应机理及动力学性质

利用密度泛函理论直接动力学方法研究了反应CH₃OCF₂CF₂OCH₃+Cl的微观机理和动力学性质. 在BB1K/6-31+G(d,p)水平上获得了反应的势能面信息, 计算中考虑了反应物CH₃OCF₂CF₂OCH₃两个稳定构象(SC1和SC2)的氢提取通道和取代反应通道. 利用改进的正则变分过渡态理论结合小曲率隧道效应(ICVT/SCT)计算了各氢提取通道的速率常数, 进而根据Boltzmann配分函数得到总包反应速率常数(k_T)以及每个构象对总反应的贡献. 结果表明296 K温度下计算的k_T(ICVT/SCT)值与已有实验值符合得很好. 由于缺乏其他温度速率常数的实验数据, 我们预测了该反应在200-2000 K温度区间内反应速率常数的三参数表达式: k_T=0.40×10^-14T^1.05exp(-206.16/T).     

Direct dynamics study on hydrogen abstraction reaction of CF3CF2CH2OH with OH radical

The dual-level direct dynamics method has been employed to investigate the H-abstraction reaction of CF(3)CF(2)CH(2)OH with OH radical, which is predicted to have two classes of possible reaction channels caused by different positions of hydrogen atom attack. The minimum-energy path is calculated at the B3LYP/6-311G(d,p) level, and the energetic information is further refined by the MC-QCISD method. To compare the structures, the other method MPW1K/6-311G(d,p) is also applied to this system. Hydrogen-bonded complexes are presented in the reactant and product sides of the three channels, indicating that each reaction may proceed via an indirect mechanism. The rate constants for each reaction channel are evaluated by canonical variational transition-state theory (CVT) with the small-curvature tunneling correction (SCT) over a wide range of temperatures from 200 to 2000 K. The calculated CVT/SCT rate constants are found to be in good agreement with the available experimental values. The result shows that the variational effect is small, and in the lower-temperature range, the SCT effect is important for each reaction. It is shown that hydrogen abstracted from the -CH(2)- position is the major channel, while H-abstraction from the -OH position may be neglected with the temperature increasing.     

Mechanistic Study of Ozone-Water Reaction and Their Reaction Products with Gas Phase Elemental Mercury:A Computational Study

The interaction of oxygen of water and central oxygen of ozone produces stable H₂O‐O₃ complex with no barrier. With decomposition of this complex through H‐abstraction by O₃ and O‐abstraction by H₂O, four possible product channels have been found. The reaction of mercury and the products of water‐ozone reaction have been studied. All geometrical and AIM parameters of intermediate, transition states, and the products of reactions are calculated and thermodynamic parameters are obtained. The negative value of free energy show that channels Hg+H₂OO, Hg+H₂O₂ and Hg+H₂O₄ in hydrogen tetroxide form (HTO) may be the main reaction channels.     

Bond-forming reactions of dications with molecules: a computational and experimental study of the mechanisms for the formation of HCF2+ from CF3(2+) and H2

The QCISD and QCISD(T) quantum chemical methods have been used to characterize the energetics of various possible mechanisms for the formation of HCF2+ from the bond-forming reaction of CF3(2+) with H2. The stationary points on four different pathways leading to the product combinations HCF2+ + H+ + F and HCF2+ + HF+ have been calculated. All four pathways begin with the formation of a collision complex [H2-CF3]2+, followed by an internal hydrogen atom migration to give HC(FH)F2(2+). In two of the mechanisms, immediate charge separation of HC(FH)F2(2+) via loss of either HF+ or a proton, followed by loss of an F atom, yields the experimentally observed bond-forming product HCF2+. For the other two mechanisms, internal hydrogen rearrangement of HC(FH)F2(2+) to give C(FH)2F(2+), followed by charge separation, yields the product CF2H+. This product can then overcome a 2.04 eV barrier to rearrange to the HCF2+ isomer, which is 1.80 eV more stable. All four calculated mechanisms are in agreement with the isotope effects and collision energy dependencies of the product ion cross sections that have been previously observed experimentally following collisions between CF3(2+) and H2/D2. We find that in this open-shell system, CCSD(T) and QCISD(T) T1-diagnostic values of up to 0.04 are acceptable. A series of angularly resolved crossed-beam scattering experiments on collisions of CF3(2+) with D2 have also been performed. These experiments show two distinct channels leading to the formation of DCF2+. One channel appears to correspond to the pathway leading to the ground state 1DCF2+ + D+ + F product asymptote and the other to the 3DCF2+ + D+ + F product asymptote, which is 5.76 eV higher in energy. The experimental kinetic energy releases for these channels, 7.55 and 1.55 eV respectively, have been determined from the velocities of the DCF2+ product ion and are in agreement with the reaction mechanisms calculated quantum chemically. We suggest that both of these observed experimental channels are governed by the reaction mechanism we calculate in which charge separation occurs first by loss of a proton, without further hydrogen atom rearrangement, followed by loss of an F atom to give the final products 1DCF2+ + D+ + F or 3DCF2+ + D+ + F.     

Multichannel RRKM study on the mechanism and kinetics of the CH3Cl + OH reaction

The kinetics and mechanism of the reaction of OH with CH₃Cl have been theoretically studied. The potential energy surface for each possible pathway has been investigated by the G2MP2 method. The rate constants for channels leading to several products have been calculated by multichannel‐Rice‐Ramsperger‐Kassel‐Marcus (RRKM) theory over a temperature range 200–2000 K. The results show the major channel is hydrogen abstraction mechanism. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

Direct dynamics study on the hydrogen abstraction reactions N2H4 + R-->N2H3 + RH (R=NH2, CH3)

We present a direct ab initio dynamics study on the hydrogen abstraction reactions N(2)H(4)+R-->N(2)H(3)+RH (R=NH(2),CH(3)), which are predicted to have six possible reaction channels for NH(2) abstraction and four for CH(3) abstraction caused by the different N(2)H(4) isomers and various attacking orientations of foreign radical to N(2)H(4). The structures and frequencies at the stationary points and the points along the minimum energy paths (MEPs) of all reaction channels are obtained at the UMP2(full)6-31+G(d,p) level of theory. Energetic information of stationary points and the points along the MEPs is further refined by means of MC-QCISD method. The rate constants of these channels are calculated using the improved canonical variational transition-state theory with the small-curvature tunneling correction (ICVT/SCT) method. The calculated results show that the favorable reaction channels are channels (n1) and (n4) as well as (c1) and (c3) (refer to Scheme 1) in the whole temperature range. The total ICVT/SCT rate constants of all channels for the two reactions at the MC-QCISDUMP2(full)6-31+G(d,p) level are both in good agreement with the available experimental data, and corresponding three-parameter expressions of k(ICVTSCT) in 220-3000 K are fitted as 6.46 x 10(-15)(T298)(3.60) exp(-386T) cm(3) mol(-1) s(-1) for NH(2) abstraction and 1.04 x 10(-14)(T298)(4.00) exp(-2037T) cm(3) mol(-1) s(-1) for CH(3) abstraction. Additionally, the long range interaction between the H atom of X-H bond in foreign radicals and the lone pair on the nonreactive N atom of the transition states is further discussed to explain the various transition-state numbers of the two similar hydrogen abstraction reactions.     

Kinetics of the hydrogen abstraction reactions of 1,1‐ and 1,2‐difluoroethane with hydroxyl radical: an ab initio study

The hydrogen abstraction reactions of 1,1‐ and 1,2‐difluoroethane with the OH radical have been investigated by the ab initio molecular orbital theory. The geometries of the reactants, products, and transition states have been optimized at the (U)MP2=full level of theory in conjunction with 6‐311G(d,p) basis functions. Single‐point (U)MP2=full with larger basis set, such as 6‐311G(3d,2p), and QCISD(T)=full/6‐311G(d,p) calculations have also been carried out to observe the effects of basis sets utilized and higher order electron correlation. Three and four reaction channels have been identified for 1,1‐ and 1,2‐difluoroethane, respectively. In the case of 1,1‐difluoroethane, hydrogen abstraction from the α‐carbon has been found to be easier than that from the β‐carbon. The barriers of the four reaction channels for 1,2‐difluoroethane are close to each other. Weak hydrogen bonding interactions have been observed between hydroxyl hydrogen and a fluorine atom in the transition states. Rate constants for the reactions of 1,1‐ and 1,2‐difluoroethane with the OH radical have been calculated using the standard transition state theory and found to be in good agreement with the experimental results. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 1305–1318, 2000     

Reactivity of C10H7+ and C10D7+ with H2 AND D2

We have investigated, both theoretically and experimentally, the reactions of naphthylium C10H7+ and d-naphthylium C10D7+ ions with H2 and D2. Cross sections as functions of the collision energy have been measured for a variety of reaction channels. Theoretical calculations have been carried out at the density functional theory level which utilizes the hybrid functional B3LYP and the split-valence 6-31G* basis set. The key features of the potential energy surfaces and the relevant thermochemical parameters have been calculated and they provide insights on the reaction mechanisms. The bimolecular reactivity of C10H7+ with H2 is dominated by the production of naphthalene cation C10H8+. The reaction is not a direct atom-abstraction process, but instead it proceeds via the formation of a stable intermediate complex C10H9+ of sigma type geometry, with a significant mobility of hydrogen along the ring. This mobility allows the scrambling of the hydrogen atoms and causes the successive statistical fragmentation of the complex into a variety of product channels. Elimination of one H(D) atom appears to be favored over elimination of one H2 or HD molecule. Alternatively, the intermediate complex can be stabilized either by collision with a third body or by emission of a photon.     

A theoretical study on the mechanism and thermodynamics of ozone-water gas phase reaction

Ozone water reaction including a complex was studied at the MP2/6-311++G(d,p) and CCSD/6-311++G(2df,2p)//MP2/6-311++G(d,p) levels of theory. The interaction between water oxygen and central oxygen of ozone produces stable H₂O-O₃ complex with no barrier. With decomposition of this complex through H-abstraction by O₃ and O-abstraction by H₂O, three possible product channels were found. Intrinsic reaction coordinate, topological analyses of atom in molecule, and vibrational frequency calculation have been used to confirm the preferred mechanism. Thermodynamic data at T = 298.15 K and atmospheric pressure have been calculated. The results show that the production of hydrogen peroxide is the main reaction channel with ΔG = ?21.112 kJ mol^-1.     

Rate Coeffcients and Branching Ratio for Multi-Channel Hydrogen Abstractions from CH3OH by F

The hydrogen abstraction reaction F+CH₃OH has two possible reaction pathways: HF+CH₃O and HF+CH₂OH. Despite the absence of intrinsic barriers for both channels, the former has a branching ratio comparable to the latter, which is far from the statistical limit of 0.25 (one out of four available H atoms). Furthermore, the measured branching ratio of the two abstraction channels spans a large range and is not quantitatively reproduced by previous theoretical predictions based on the transition-state theory with the stationary point information calculated at the levels of M?ller-Plesset perturbation theory and G2. This work reports a theoretical investigation on the kinetics and the associated branching ratio of the two competing channels of the title reaction using a quasi-classical trajectory approach on an accurate full-dimensional potential energy surface (PES) fitted by the permutation invariant polynomial-neural network approach to ca. 1.21x10⁵ points calculated at the explicitly correlated (F12a) version of coupled cluster singles doubles and perturbative triples (CCSD(T)) level with the aug-cc-pVDZ basis set. The calculated room temperature rate coeffcient and branching ratio of the HF+CH₃O channel are in good agreement with the available experimental data. Furthermore, our theory predicts that rate coeffcients have a slightly negative temperature dependence, consistent with barrierless nature of the reaction.     

Ab initio direct dynamics studies on the reactions of chlorine atom with CH3-nFnCH2OH (n = 1-3)

The hydrogen abstraction reactions of Cl atom with a series of fluorinated alcohols, i.e., CH(3-n)F(n)CH(2)OH + Cl (n = 1-3) (R1-R3) have been studied systematically by ab initio direct dynamics method and the canonical variational transition state theory (CVT). The potential energy surface information is calculated at the MP2/6-311G(d,p) level. Energies along the minimum energy paths are improved by a series of single-point calculations at the higher modified GAUSSIAN-2 (G2M) level of theory. Theoretical analysis shows that three kinds of hydrogen atoms can be abstracted from the reactants CH(2)FCH(2)OH and CHF(2)CH(2)OH, and for CF(3)CH(2)OH, two possible pathways are found. The rate constants for each reaction channel are evaluated by CVT with the small-curvature tunneling correction (SCT) over a wide range of temperature from 200 to 2000 K. The calculated CVT/SCT rate constants are in good agreement with the available experimental values for the reactions CHF(2)CH(2)OH + Cl and CF(3)CH(2)OH + Cl. However, for the reaction CH(2)FCH(2)OH + Cl, there is negative temperature dependence below 500 K, which is different from the experimental fitted. It is shown that in the low temperature ranges, the three reactions all proceed predominantly via H-abstraction from the methylene positions, and with the increase of the temperature the H-abstraction channels from the fluorinated-methyl positions should be taken into account, while the H-abstraction channels from the hydroxyl groups are negligible over the whole temperature ranges. Also, the reactivity decreases substantially with fluorine substitution at the methyl position of alcohol.     

Density Functional Theoretical Study on the Reaction Mechanism of HOCO and Its Radical

In this paper, the reaction mechanism of HOCO with its radical has been investigated deeply by density functional theory(DFT), while the geometries and harmonic vibration frequencies of the reactants, intermediates, transition states and products have been calculated at the B3LYP/6-311++G** level. The CCSD(T)/cc-p VQZ method was used to further calculate the single-point energy of each stationary point along the reaction channel. The result shows that channels(b) and(d) via carbon-carbon intermediates are competing with the ones(e),(f),(i) and(l) through carbon-oxygen intermediates, and the six channels are dominant for the title reaction to produce HCOOH and CO2. The result also indicates that the hydrogen abstraction from the HOCO radical can be performed quickly by the self-interaction of HOCO in the absence of other radicals or atoms, which indirectly proves for the first time that the ground-state cis-HOCO radical is not decomposed by the tunneling effect. In addition, the charges of natural population analysis(NPA) and the variation of chemical bonds have been analyzed by the Natural Bond Orbital(NBO) method along with the important reaction path.     

Theoretical Studies on the Reactions of 1,1,2,2,3,3,4‐Heptafluorocyclopentane with Hydroxyl and Hydrogen Free Radicals

1,1,2,2,3,3,4‐Heptafluorocyclopentane (F7A) has considerable potential to be a new halon replacement due to its environmental friendliness and low‐toxicity. However, the reaction processes of F7A with hydroxyl and hydrogen free radicals, which are of great importance for investigating its fire suppression mechanisms, are still unclear. In this paper, ab inito and density functional theory are used to deduce the possible reaction pathways for the reactions of F7A with hydroxyl and hydrogen free radicals at the CCSD/cc‐pVDZ//B3LYP/6‐311++G (d,p) level of theory. Two distinct reaction pathways including ten elementary reaction channels for F7A with hydroxyl free radical, and five distinct reaction pathways including twenty elementary reaction channels for F7A with hydrogen free radical are investigated. The geometries, vibrational frequencies and reaction energy barriers are also determined. Based on the calculated results, the possible reaction mechanisms are proposed and discussed. The most feasible reaction channel for F7A with hydroxyl free radical is that leads to CH(OH)CH₂(CF₂)₃+·F, and the most feasible reaction channel for F7A with hydrogen free radical is that leads to (CF₂)₃CH₂CH·+HF. The study is helpful to further study its fire suppression mechanisms and promote it to be a new generation of halon replacement.     

Dual-level direct dynamics studies for the reactions of OH radical with bromine-substituted ethanes

The dynamic properties of the multichannel hydrogen abstraction reactions of CH(3)CH(2)Br + OH --> products and CH(3)CHBr(2) + OH --> products are studied by dual-level direct dynamics method. For each reaction, three reaction channels, one for alpha-hydrogen abstraction and two for beta-hydrogen abstractions, have been identified. The minimum energy paths (MEPs) of both the reactions are calculated at the Becke's half-and-half (BH&H)-Lee-Yang-Parr (LYP)/6-311G(d, p) level and the energy profiles along the MEPs are further refined with interpolated single-point energies (ISPE) method at the G2M(RCC5)//BH&H-LYP level. There are complexes with energies less than those of the reactants or products located at the entrance or exit channels, which indicates that the reactions may proceed via an indirect mechanism. By canonical variational transition-state theory (CVT) the rate constants are calculated incorporating the small-curvature tunneling (SCT) correction in the temperature range of 220-2000 K. The agreement of the rate constants with available experimental values for two reactions is good in the measured temperature range. The calculated results show that alpha-hydrogen abstraction channel is the major reaction pathway in the lower temperature for two reactions, while the contribution of beta-hydrogen abstraction will increase with the increase in temperature.