Thermochemical kinetics of bromine nitrate,bromine nitrite,halogen hydroperoxides,dichlorine pentoxide,peroxycarboxylic acids,and diacyl peroxides

Heats of formation of BrONO₂, BrONO, BrOOH, FOOH, FOOCl, CF₃C(O)OOH, HC(O)OOH, CH₃C(O)OOH, and [CH₃C(O)O]₂ are estimated from bond contributions taken from J. Phys. Chem., 100, 10150 (1996). They agree within ±2 kcal/mol with recent experimental or ab initio data. The resulting BDE(O(SINGLEBOND)O)=36 kcal/mol value in diacetyl peroxide requires the concerted assistance of exothermic C(SINGLEBOND)C(O) weakening in the transition state of its decomposition into free radicals. It also implies the existence of a previously unrecognized 12 kcal/mol nonbonded repulsion in acyl anhydrides. The formation of chloryl chlorate with ΔH_f(O₂ClOClO₂)=50 kcal/mol, a marginally stable species toward dissociation into (ClO₃+OClO), may account for observations made in the [O(³P+OClO] system at low temperatures. © 1998 John Wiley & Sons, Inc. Int J Chem Kinet: 30: 41–45, 1998.     

Reaction of CH3O2 and HO2: Ab initio characterization of dimer structure and vibrational mode analysis for reaction mechanisms

All species involved in the multi‐channel reaction of CH₃O₂ with HO₂ have been investigated using density functional theory (DFT). The molecular geometries for various species are optimized employing the B3LYP method implementing the 6‐311++G** basis set. The relative energies of all species are calculated at the same level theory. The results show that there are two kinds of channels: singlet and triplet. The singlet channel involves four intermediates and six transition states. The triplet channel includes two intermediates and two transition states. There are four kinds of reaction products: CH₃OOH + ¹O₂, CH₃OH + O₃, CH₄ + 2O₂, and CH₃OOH + ³O₂. The vibrational mode analysis is used to elucidate the relationships of the intermediates, the transition states, and the products. The extensive investigation shows that the reaction mechanism is reliable. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

DFT and ab initio study on the reaction mechanism of CH2SH + O2

The mechanism for the CH₂SH + O₂ reaction was investigated by DFT and ab initio chemistry methods. The geometries of all possible stationary points were optimized at the B3LYP/6-311+G(d,p) level, and the single point energy was calculated at the CCSD(T)/cc-pVXZ(X = D and T), G3MP2 and BMC-CCSD levels. The results indicate that the oxidation of CH₂SH by O₂ to form HSCH₂OO is a barrierless process. The most favorable channel is the rearrangement of the initial adduct HSCH₂OO (IM1) to form another intermediate H₂C(S)OOH (IM3) via a five-center transition state, and then the C–O bond fission in IM3 leads to a complex CH₂S. . .HO₂ (MC1), which finally gives out to the major product CH₂S + HO₂. Due to high barriers, other products including cis- and trans-HC(O)SH + HO could be negligible. The direct abstraction channel was also determined to yield CH₂S + HO₂, with the barrier height of 22.3, 18.1 and 15.0 kcal/mol at G3MP2, CCSD(T)/cc-pVTZ and BMC-CCSD levels, respectively, it is not competitive with the addition channel, in which all stationary points are lower than reactant energetically. The other channels to produce cis- and trans-CHSH + HO₂ are also of no importance.     

Insertion reaction of propene into Rh?H bond in HRh(CO)(PH3)2(C3H6) compound: A density functional study

Quantum mechanical calculations at the MP4 (SDQ) level using the BP86‐optimized geometries were carried out to investigate the energies and reaction mechanism for the propene (CH₃ C¹H CH$^{\mathrm{2}}_{\mathrm{2}}$) insertion reaction into the Rh H bond, using the cis‐HRh(CO)(PH₃)₂ compound as a model catalytic species. Since the reaction may occur on the branched carbon 1 or in the normal carbon 2 , which leads to branched and normal Rh(alkyl) compounds, respectively, we investigated these two mechanisms. The results show that the insertion in the branched carbon has an activation energy of 16.2 kcal/mol, and the activation energy for the reaction to take place at the normal carbon is 14.3 kcal/mol. These activation energies, together with the calculated relative energy of the metal–alkyl compounds formed after the insertion considering these two pathways, were used to access the regioselectivity on this reaction. We found a ratio of normal‐ and iso‐products, n:iso, of (96:4), which is in excellent agreement with the experimental regioselectity of (95:5). © 2000 John Wiley & Sons, Inc. Int J Quant Chem 78: 42–51, 2000     

Mechanism of the chain process forming H2 in the photooxidation of formaldehyde

A mechanism for the formation in a chain of H₂, CO, and HCOOH in the photooxidation of formaldehyde is proposed. This mechanism is initiated by the addition of HO₂ to formaldehyde. Hydrogen atoms are produced by the thermal dissociation of the HOCH₂O radical: HOCH₂O → H + HCOOH; ΔH = + 3.2 kcal/mol [5]. Photolysis of CH₂O? O₂? NO mixtures and product analysis were carried out in conjunction with kinetic simulation yielding an estimate for the activation energy of the dissociation reaction : E₅ = 14.9 ± 1.0 kcal/mol. Previous observations of this chain process are considered in view of this mechanism.     

CH3OOH和CH3CH2OOH的表征和活性：光电子能谱的研究

The ionization energies of MHP （CH3OOH） and EHP（CH3CH2OOH） nave been determined by Hel photoelectron spectroscopy （PES） measurement and both Gaussian-2 （G2） calculation and Hartree-Fock （HF） method on the basis of Koopmans theorem at 6.311＋G^＊ basis set level for the first time. The assignment and characterization of PE spectra of MHP and EHP were also supported by the G2 and HF calculations. The first ionization energies of MHP and EHP are 9.87 and 9.65 eV, respectively. Higher solubility of EHP in the atmosphere was attributed to their lower ionization energy values.     

Effects of multidimensional tunneling in the kinetics of hydrogen abstraction reactions of O (3P) with CH3OCHO

Quantum tunneling paths are important in reactions when there is a significant component of hydrogenic motion along the potential energy surface. In this study, variational transition state with multidimensional tunneling corrections are employed in the calculations of the thermal rate constants for hydrogen abstraction from the cis‐CH₃OCHO by O (³P) giving CH₃OCO + OH (R1) and CH₂OCHO + OH (R2). The structures and electronic energies are computed with the M06‐2X method. Benchmark calculations with the CBS_D–T approach give an enthalpy of reaction at 0 K for R1 (−2.8 kcal/mol) and R2 (−2.5 kcal/mol) which are in good agreement with the experiment, i.e. −2.61 and −1.81 kcal/mol. At the low and intermediate values of temperatures, small‐ and large‐curvature tunneling dominate the kinetics of R1, which is the dominant path over the range of temperature from 250 to 1200 K. This study shows the importance of multidimensional tunneling corrections for both R1 and R2, for which the total rate constant at 298 K calculated with the CVT/μOMT method is 8.2 × 10⁻¹⁵ cm³ molecule⁻¹ s⁻¹ which agrees well with experiment value of 9.3 × 10⁻¹⁵ cm³ molecule⁻¹ s⁻¹ (Mori, Bull. Inst. Chem. Res. 1981, 59, 116). © 2018 Wiley Periodicals, Inc.     

Ab initio calculations on the barrier height for the hydrogen addition to ethylene and formaldehyde. The importance of spin projection

Heats of reaction and barrier heights have been computed for H + CH₂CH₂ → C₂H₅, H + CH₂O → CH₃O, and H + CH₂O → CH₂OH using unrestricted Hartree-Fock and Møller–Plesset perturbation theory up to fourth order (with and without spin annihilation), using single-reference configuration interaction, and using multiconfiguration self-consistent field methods with 3-21G, 6-31G(d), 6-31G(d,p), and 6-311G(d,p) basis sets. The barrier height in all three reactions appears to be relatively insensitive to the basis sets, but the heats of reaction are affected by p-type polarization functions on hydrogen. Computation of the harmonic vibrational frequencies and infrared intensities with two sets of polarization functions on heavy atoms [6-31G(2d)] improves the agreement with experiment. The experimental barrier height for H + C₂H₄ (2.04 ± 0.08 kcal/mol) is overestimated by 7?9 kcal/mol at the MP2, MP3, and MP4 levels. MCSCF and CISD calculations lower the barrier height by approximately 4 kcal/mol relative to the MP4 calculations but are still almost 4 kcal/mol too high compared to experiment. Annihilation of the largest spin contaminant lowers the MP4SDTQ computed barrier height by 8?9 kcal/mol. For the hydrogen addition to formaldehyde, the same trends are observed. The overestimation of the barrier height with Møller-Plesset perdicted barrier heights for H + C₂H₄ → C₂H₅, H + CH₂O → CH₃O, and H + CH₂O → CH₂OH at the MP4SDTQ /6-31G(d) after spin annihilation are respectively 1.8, 4.6, and 10.5 kcal/mol.     

UHF-CI studies of energy barriers for the abstraction and exchange reactions in the system H + CH4

UHF and CI calculations, using the direct CI method, and double-zeta plus polarization functions basis sets, have been performed on the more important parts of the energy hypersurface for CH₅. The abstraction H + CH₄ → H₂ + CH₃ and the inversion substitution reaction H′ + CH₄ → CH₃H′ + H have been studied in detail. The predicted barriers for these two reactions are 13.5 and 36.6 kcal/mol, respectively. The abstraction reaction is, in agreement with experiment, found to be almost thermo-neutral with a heat of reaction of 1.5 kcal/mol.     

On the mechanism of thermal oxidation of methane

Using the method of freezing radicals in conjunction with ESR spectroscopic measurements, the kinetics of the thermal oxidation of methane has been studied under atmospheric pressure depending on the temperature, composition of the mixture, and nature of the surface of the reaction vessel. It has been shown that in a reactor treated with boric acid, the intermediates methylhydroperoxide and hydrogen peroxide are responsible for chain branching. It has been established that the leading active centers of the reaction are the HO₂ radicals, while chain branching occurs as a result of the decomposition of peroxy compounds—methylhydroperoxide and hydrogen peroxide. In reactors treated with potassium bromide, the concentrations of radicals and peroxy compounds were found to be lower than the sensitivity of the method of measurement. Computations were performed for the scheme of methane oxidation at 738 K for a reactor treated with boric acid. Satisfactory agreement was found between the experimental and computed kinetic curves of accumulation of main intermediates CH₂O, H₂O₂, CH₃OOH. The influence of their addition on the kinetics of the reaction has been considered. It has been shown that the addition of formaldehyde does not lead to chain branching, however; it contributes to the formation of those peroxy compounds that bring about chain branching. Mathematical modeling confirmed conclusions made on the basis of experimental data concerning the nature of the leading active centers and the products that are responsible for the degenerate chain branching.     

Vinyloxyborane and its isomers. An ab initio study of the C2H5BO potential energy surface,the barrier to 1,3-shifts in β-ketoboranes,and the mechanism of the carbonylation reaction of boranes

Vinyloxyboranes, CH₂?CH? ;O? ;BR₂, are shown by ab initio molecular orbital theory to be more stable than the isomeric β-aldoboranes, R₂B? CH₂? CH?O, by ca. 19 kcal/mol. The MP2/6-31G*/6-31G* + ZPE barrier for the [1,3] boron shift is only 10.9 kcal/mol (R ? Me) relative to the aldoborane. Other C₂H₅BO isomers (β-ketoboranes, boraepoxides and organoboron oxides), which are related to the proposed stages in the carbonylation reaction of boranes, are shown to be plausible intermediates. However, some of the computed barriers for methyl group migrations are unrealistically large, up to ca. 63 kcal/mol.     

Theoretical reinvestigation of the ozonolysis mechanism of allyl alcohol

Allyl alcohol (AA) is the simplest unsaturated alcohol. Ozonolysis is one of the key removal processes for AA in the atmosphere. However, a recent theoretical study suggests that the ozonolysis of AA cannot feasibly occur in atmospheric conditions because of the high barrier (~96 kcal/mol) involved in the primary ozonide (POZ) decomposition. In this work, the ozonolysis mechanism of AA was reinvestigated theoretically. The computed barrier for POZ decomposition is only ~20 kcal/mol. Therefore, the AA ozonolysis can take place in the atmosphere, consistent with the experimental conclusions. Moreover, two new Criegee intermediates (syn- and anti-AA-CI) were found to be produced in this reaction. The wave function analyses indicate that there exists an intermolecular hydrogen bond in syn-AA-CI, which significantly affects its unimolecular decomposition and reactions with H₂O. Compared with the normal reactions of anti-CI-AA, the stabilized syn-AA-CI has two distinct isomerization channels: (i) addition of OH group to the reactive sites of CI forming an ethylene oxide (HOOCH₂OCH₂) and (ii) double H-transfer producing HOOCH₂CHO. Meanwhile, the addition of H₂O in syn-AA-CI also exhibits two different pathways. One is the unique addition-coupled double H-transfer, and the other is the addition-coupled single H-transfer, both leading to the formation of CH₂(OH)CH(OH)OOH.     

Theoretical study on the atmospheric reaction of CH3O2 with OH

A quantum chemical investigation on the reaction mechanism of CH₃O₂ with OH has been performed. Based on B3LYP and QCISD(T) calculations, seven possible singlet pathways and seven possible triplet pathways have been found. On the singlet potential energy surface (PES), the most favorable channel starts with a barrierless addition of O atom to CH₃O₂ leading to CH₃OOOH and then the O? O bond dissociates to give out CH₃O + HO₂. On the triplet PES, the calculations indicate that the dominant products should be ³CH₂O₂ + H₂O with an energy barrier of 29.95 kJ/mol. The results obtained in this work enrich the theoretical information of the title reaction and provide guidance for analogous atmospheric chemistry reactions. © 2015 Wiley Periodicals, Inc.     

Theoretical study of the free‐radical copolymerization of styrene with methyl methacrylate: Comparative study to the styrene–acrylonitrile monomer system

Enthalpic and electronic terminal and penultimate unit effects in the free‐radical copolymerization of styrene (S) with methyl methacrylate (M) were investigated by quantum mechanical calculations at 0 and 298 K. Total energies, zero‐point energies scaled by a 0.96 factor, and thermal enthalpy corrections for all optimized structures were computed at the B3‐LYP/6‐31G(d) level of theory. Differences in enthalpies for elementary propagation reactions at 0 and 298 K did not exceed 0.6 kcal/mol. Enthalpic effects of the replacement of S by M in the penultimate position of the growing radicals in elementary copolymerization propagation reactions (enthalpic penultimate unit effects) were always positive, ranging from 1.2 to 3.3 kcal/mol at 298 K. The values suggested that the elementary propagation reactions involving more S units in the growing polymer chain ends should be slightly thermodynamically preferred. A comparison of these results with those for the S–acrylonitrile monomer system showed that the most crucial feature differentiating enthalpic effects for the two monomer systems is the replacement of M by acrylonitrile in the reaction pair CH₃‐S‐M · + M → CH₃‐S‐M‐M · and CH₃‐M‐M · + M → CH₃‐M‐M‐M ·. © 2004 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 42: 1557–1565, 2004     

Thermochemistry and kinetics of the 2-butanone-4-yl CH3C(=O)CH2CH2• + O2 reaction system

Thermochemistry and kinetic pathways on the 2-butanone-4-yl (CH₃C(=O)CH₂CH₂•) + O₂ reaction system are determined. Standard enthalpies, entropies, and heat capacities are evaluated using the G3MP2B3, G3, G3MP3, CBS-QB3 ab initio methods, and the B3LYP/6-311g(d,p) density functional calculation method. The CH₃C(=O)CH₂CH₂• radical + O₂ association reaction forms a chemically activated peroxy radical with 35 kcal mol⁻¹ excess of energy. The chemically activated adduct can undergo RO−O bond dissociation, rearrangement via intramolecular hydrogen transfer reactions to form hydroperoxide-alkyl radicals, or eliminate HO₂ and OH. The hydroperoxide-alkyl radical intermediates can undergo further reactions forming ketones, cyclic ethers, OH radicals, ketene, formaldehyde, or oxiranes. A relatively new path showing a low barrier and resulting in reactive product sets involves peroxy radical attack on a carbonyl carbon atom in a cyclic transition state structure. It is shown to be important in ketones when the cyclic transition state has five or more central atoms.     

The remarkably invariant interaction energies of lithium first-row compounds with water and with ammonia

Solvation energies of lithium first-row compounds LiX (X ? H, Li, BeH, BH₂, CH₃, NH₂, OH, F) and of the lithium cation with the model solvents, water and ammonia, have been calculated ab inito (MP2/6-31 + G^*//6-31G^* with zero-point vibrational energy corrections at 3-21G//3-21G). The solvation energies are found to be remarkably constant: ?18.0 ± 1.2 and ?21.5 ± 1.3 kcal/mol for the hydrates and ammonia solvates, respectively. This independence on the nature of X is due largely to the ionic character of the LiX compounds (dipole moments 4.7–6.6 debye). The unexpectedly high solvation energies of the lithium molecule (?14.3 and ?17.8 kcal/mol, respectively) are due to the polarizability of Li₂. At the same level, the lithium cation has interaction energies with H₂O and NH₃ of ?34.1 and ?39.7 kcal/mol, respectively. For the hydrates of LiOH and LiF cyclic structures with hydrogen bonds and somewhat increased solvation energies also are described.     

Variational transition state theory rate constants and H/D kinetic isotope effects for CH3 + CH3OCOH reactions

The rate constants and H/D kinetic isotope effect for hydrogen abstraction reactions involving isotopomers of methyl formate by methyl radical are computed employing methods of the variational transition state theory (VTST) with multidimensional tunneling corrections. The energy paths were built with a dual-level method using the moller plesset second-order perturbation theory (MP2) method as the low-level and complete basis set (CBS) extrapolation as the high-level energy method. Benchmark calculations with the CBS_D-T approach give an enthalpy of reaction at 0 K for R1 (−4.5 kcal/mol) and R2 (−4.2 kcal/mol) which are in good agreement with the experiment, that is, −4.0 and − 4.8 kcal/mol. For the reactional paths involving the isotopomers CH₃ + CH₃OCOH → CH₄ + CH₃OCO and CH₃ + CH₃OCOD → CH₃D + CH₃OCO, the value of k^H/k^D (T = 455 K) using the canonical VTST/small-curvature tunneling approximation method is 6.7 in close agreement with experimental value (6.2). © 2019 Wiley Periodicals, Inc.     

Analysis of the intermolecular interaction between CH(3)OCH(3), CF(3)OCH(3), CF(3)OCF(3), and CH(4): high level ab initio calculations

The intermolecular interaction energies of the CH₃OCH₃? CH₄, CF₃OCH₃? CH₄, and CF₃OCF₃? CH₄ systems were calculated by ab initio molecular orbital method with the electron correlation correction at the second order Møller–Plesset perturbation (MP2) method. The interaction energies of 10 orientations of complexes were calculated for each system. The largest interaction energies calculated for the three systems are ?1.06, ?0.70, and ?0.80 kcal/mol, respectively. The inclusion of electron correlation increases the attraction significantly. It gains the attraction ?1.47, ?1.19, and ?1.27 kcal/mol, respectively. The dispersion interaction is found to be the major source of the attraction in these systems. In the CH₃OCH₃? CH₄ system, the electrostatic interaction (?0.34 kcal/mol) increases the attraction substantially, while the electrostatic energies in the other systems are not large. Fluorine substitution of the ether decreases the electrostatic interaction, and therefore, decreases the attraction. In addition the orientation dependence of the interaction energy is decreased by the substitution. © 2002 Wiley Periodicals, Inc. J Comput Chem 23: 1472–1479, 2002     

Clustering and activation in reactions of CoCp with hydrogen and methane

Gas-phase clustering reactions of CoCp⁺ with H₂ and with CH₄ were investigated using temperature-dependent equilibrium experiments. In both systems, the CoCp⁺ ion was found to form strong interactions with two ligands. The first and second H₂ groups cluster to CoCp⁺ with bond energies of 16.2 and 16.8 kcal/mol, respectively, while the first and second CH₄ groups cluster to CoCp⁺ with bond energies of 24.1 and 12.1 kcal/mol, respectively. These bond energies are in good agreement with those determined by density functional theory (DFT). Molecular geometries for the four clusters determined with DFT are also presented. Weak experimental bond energies of 0.9 kcal/mol for the third H₂ and 2.2 kcal/mol for the third CH₄ clustering to CoCp⁺ suggest these ligands occupy the second solvation shell of the ion. In addition to clustering in the methane system, H₂-elimination from CoCp(CH₄)₂⁺ was observed. The mechanism for this reaction was investigated by collision-induced dissociation experiments and DFT, which suggest the predominate H₂-elimination product is (c-C₅H₆)Co⁺---C₂H₅. Theory indicates that dehydrogenation requires the active participation of the Cp ring in the mechanism. Transfer of H and CH₃ groups to the C₅-ring ligand allows the metal center to avoid the high-energy Co(IV) oxidation state required when it forms two covalent bonds in addition to its interaction with a C₅-ring ligand.     

Ab initio study of the potential energy surface and product branching ratios for the reaction of O(D) with CH3CH2F

The potential energy surface of O(¹D) + CH₃CH₂F reaction has been studied using QCISD(T)/6-311++G(d,p)//MP2/6-311G(d,p) method. The calculations reveal an insertion–elimination reaction mechanism of the title reaction. The insertion process has two possibilities: one is the O(¹D) atom inserting into C–F bond of CH₃CH₂F produces one energy-rich intermediate CH₃CH₂OF and another is the O(¹D) atom inserting into one of the C–H bonds of CH₃CH₂F produces two energy-rich intermediates, IM1 and IM2. The three intermediates subsequently decompose to various products. The calculations of the branching ratios of various products formed though the three intermediates have been carried out using RRKM theory at the collision energies of 0, 5, 10, 15, 20, 25 and 30 kcal/mol. CH₃CH₂O is the main decomposition product of CH₃CH₂OF. HF and CH₃ are the main decomposition products for IM1; CH₂OH is the main decomposition product for IM2. Since IM1 is more stable and more likely to form than CH₃CH₂OF and IM2, HF and CH₃ are probably the main products of the O(¹D) + CH₃CH₂F reaction. Our computational results can give insight to reaction mechanism and provide probable explanations for future experiments.