Mechanism and Rate Constants of the CH3+ CH2CO Reaction: A Theoretical Study

The mechanism of the reaction of ketene with methyl radical has been studied by ab initio CCSD(T)‐F12/cc‐pVQZ‐f12//B2PLYPD3/6‐311G** calculations of the potential energy surface. Temperature‐ and pressure‐dependent reaction rate constants have been computed using the Rice–Ramsperger–Kassel–Marcus (RRKM)–Master Equation and transition state theory methods. Three main channels have been shown to dominate the reaction; the formation of the collisionally stabilized CH₃COCH₂ radical and the production of the C₂H₅ + CO and HCCO + CH₄ bimolecular products. Relative contributions of the CH₃COCH₂, C₂H₅ + CO, and HCCO + CH₄ channels strongly depend on the reaction conditions; the formation of thermalized CH₃COCH₂ is favored at low temperatures and high pressures, HCCO + CH₄ is dominant at high temperatures, whereas the yield of C₂H₅ + CO peaks at intermediate temperatures around 1000 K. The C₂H₅ + CO channel is favored by a decrease in pressure but remains the second most important reaction pathway after HCCO + CH₄ under typical flame conditions. The calculated rate constants at different pressures are proposed for kinetic modeling of ketene reactions in combustion in the form of modified Arrhenius expressions. Only rate constant to form CH₃COCH₂ depends on pressure, whereas those to produce C₂H₅ + CO and HCCO + CH₄ appeared to be pressure independent.     

Theoretical study on the gas‐phase reaction mechanism between rhodium monoxide and methane for methanol production

The gas‐phase reaction mechanism between methane and rhodium monoxide for the formation of methanol, syngas, formaldehyde, water, and methyl radical have been studied in detail on the doublet and quartet state potential energy surfaces at the CCSD(T)/6‐311+G(2d, 2p), SDD//B3LYP/6‐311+G(2d, 2p), SDD level. Over the 300–1100 K temperature range, the branching ratio for the Rh(⁴F) + CH₃OH channel is 97.5–100%, whereas the branching ratio for the D‐CH₂ORh + H₂ channel is 0.0–2.5%, and the branching ratio for the D‐CH₂ORh + H₂ channel is so small to be ruled out. The minimum energy reaction pathway for the main product methanol formation involving two spin inversions prefers to both start and terminate on the ground quartet state, where the ground doublet intermediate CH₃RhOH is energetically preferred, and its formation rate constant over the 300–1100 K temperature range is fitted by k_CH3RhOH = 7.03 × 10⁶ exp(?69.484/RT) dm³ mol^?1 s^?1. On the other hand, the main products shall be Rh + CH₃OH in the reactions of RhO + CH₄, CH₂ORh + H₂, Rh + CO +2H₂, and RhCH₂ + H₂O, whereas the main products shall be CH₂ORh + H₂ in the reaction of Rh + CH₃OH. Meanwhile, the doublet intermediates H₂RhOCH₂ and CH₃RhOH are predicted to be energetically favored in the reactions of Rh + CH₃OH and CH₂ORh + H₂ and in the reaction of RhCH₂ + H₂O, respectively. © 2009 Wiley Periodicals, Inc. J Comput Chem 2010     

Direct ab initio dynamics calculations on the rate constants for the hydrogen-abstraction reaction of C2H5F with O (3P)

The hydrogen-abstraction reaction C₂H₅F+O → C₂H₄F+OH has been studied by a dual-level direct dynamics method. For the reaction, three reaction channels, one for α-abstraction and two for β-abstraction, have been identified. The potential-energy surface information is obtained at the MP2(full)/6-311G(d,p) and PMP2(full)/6-311G(3df,3pd) (single-point) levels. By canonical variational transition-state theory, rate constants for each reaction channel are calculated with a small-curvature tunneling correction. The total rate constant is calculated from the sum of the individual rate constants and the temperature dependence of the branching ratios is obtained over a wide range of temperatures from 300 to 5,000 K. The agreement of the rate constants with experiment is good in the experimental temperature range from 1,000 to 1,250 K. The calculated results indicate that at low temperatures α-abstraction is most likely to be the major reaction channel, while β-abstraction channels will significantly contribute to the whole reaction rate as the temperature increases. Received: 23 January 2002 / Accepted: 23 June 2002 / Published online: 20 September 2002     

A Theoretical Investigation of the Decomposition Reactions of Ethyl Formate in the S0 State

This study revisits the stability of the possible conformations and the decomposition reactions of ethyl formate in the S₀ state using the (U)MP2, MP4SDTQ, CCSD(T), and (U)B3LYP methods with various basis sets. The transition states of the decomposition channels to HCOOH + C₂H₄, CO + CH₃CH₂OH, CH₂O + CH₃CHO, HCOH + CH₃CHO, C₂H₆ + CO₂, and H₂ + CH₂CHOCHO are determined. The microcanonical rate constants derived from the RRKM theory are calculated for each of the decomposition reactions. The high‐pressure limit rate constants are calculated for the decomposition channels to HCOOH + C₂H₄, CO + CH₃CH₂OH, and CH₂O + CH₃CHO.     

Mechanism and rate constants of the CH2 + CH2CO reactions in triplet and singlet states: A theoretical study

Ab initio and density functional CCSD(T)-F12/cc-pVQZ-f12//B2PLYPD3/6-311G** calculations have been performed to unravel the reaction mechanism of triplet and singlet methylene CH₂ with ketene CH₂CO. The computed potential energy diagrams and molecular properties have been then utilized in Rice–Ramsperger–Kassel–Marcus-Master Equation (RRKM-ME) calculations of the reaction rate constants and product branching ratios combined with the use of nonadiabatic transition state theory for spin-forbidden triplet-singlet isomerization. The results indicate that the most important channels of the reaction of ketene with triplet methylene lead to the formation of the HCCO + CH₃ and C₂H₄ + CO products, where the former channel is preferable at higher temperatures from 1000 K and above. In the C₂H₄ + CO product pair, the ethylene molecule can be formed either adiabatically in the triplet electronic state or via triplet-singlet intersystem crossing in the singlet electronic state occurring in the vicinity of the CH₂COCH₂ intermediate or along the pathway of CO elimination from the initial CH₂CH₂CO complex. The predominant products of the reaction of ketene with singlet methylene have been shown to be C₂H₄ + CO. The formation of these products mostly proceeds via a well-skipping mechanism but at high pressures may to some extent involve collisional stabilization of the CH₃CHCO and cyclic CH₂COCH₂ intermediates followed by their thermal unimolecular decomposition. The calculated rate constants at different pressures from 0.01 to 100 atm have been fitted by the modified Arrhenius expressions in the temperature range of 300–3000 K, which are proposed for kinetic modeling of ketene reactions in combustion. © 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.     

Probing the pyrolysis of ethyl formate in the dilute gas phase by synchrotron radiation and theory

The thermal decomposition of the atmospheric constituent ethyl formate was studied by coupling flash pyrolysis with imaging photoelectron photoion coincidence (iPEPICO) spectroscopy using synchrotron vacuum ultraviolet (VUV) radiation at the Swiss Light Source (SLS). iPEPICO allows photoion mass-selected threshold photoelectron spectra (ms-TPES) to be obtained for pyrolysis products. By threshold photoionization and ion imaging, parent ions of neutral pyrolysis products and dissociative photoionization products could be distinguished, and multiple spectral carriers could be identified in several ms-TPES. The TPES and mass-selected TPES for ethyl formate are reported for the first time and appear to correspond to ionization of the lowest energy conformer having a cis (eclipsed) configuration of the O = C (H)– O – C (H₂)–CH₃ and trans (staggered) configuration of the O= C (H)– O – C (H₂)– C H₃ dihedral angles. We observed the following ethyl formate pyrolysis products: CH₃CH₂OH, CH₃CHO, C₂H₆, C₂H₄, HC(O)OH, CH₂O, CO₂, and CO, with HC(O)OH and C₂H₄ pyrolyzing further, forming CO + H₂O and C₂H₂ + H₂. The reaction paths and energetics leading to these products, together with the products of two homolytic bond cleavage reactions, CH₃CH₂O + CHO and CH₃CH₂ + HC(O)O, were studied computationally at the M06-2X-GD3/aug-cc-pVTZ and SVECV-f12 levels of theory, complemented by further theoretical methods for comparison. The calculated reaction pathways were used to derive Arrhenius parameters for each reaction. The reaction rate constants and branching ratios are discussed in terms of the residence time and newly suggest carbon monoxide as a competitive primary fragmentation product at high temperatures.     

Reaction of propionitrile with methylidyne: A theoretical study

The results of a CCSD(T)-F12/cc-pVTZ-f12//ωB97XD/cc-pVTZ quantum-chemical study of the potential energy surface (PES) for the reaction of propionitrile with methylidyne are combined with Rice-Ramsperger-Kassel-Marcus (RRKM) calculations of the reaction rate constants and product branching ratios in the deep space conditions corresponding to the zero-pressure limit at various collision energies. The most energetically favorable reaction pathways have been identified. The reaction outcome has been shown to strongly depend on the branching in the entrance reaction channel, between CH additions and insertions into various C-H and C-C bonds. For instance, CH addition to the N atom predominantly leads to 3H-pyrrole + H (p9), with CH₂NC + C₂H₄ (p2) also being a significant product. CH addition to the triple C≡N bond mostly results in 2-methylene-2H-azirine + CH₃ (p13), whereas CH insertions into C-H bonds in the CH₃ and CH₂ groups of propionitrile form CH₂CN + C₂H₄ (p1) and CH₂CHCN + CH₃ (p7) respectively. Less likely CH insertions into single C-C bonds yield CH₃CHCHCN + H (p5) and CH₂CHCH₂CN + H (p8). The results indicate that the methylidyne + propionitrile reaction may represent a critical step toward the formation of heterocyclic N-containing molecules in the interstellar medium and in planetary atmospheres.     

The √t law in solid-phase reactions. Analysis of the hypothesis on rate constant distribution

The reaction C₂H₅ + O₂ → C₂H₅O₂ in glassy methanol-d₄ and the H-atom abstraction by CH₃, C₂H₅, and n-C₄H₉ radicals in C₂H₅OH + C₂D₅OH and CD₃CH₂OH + C₂D₅OH glassy mixtures have been studied by electron spin resonance. The analysis of the dependence of the reaction rates on the concentration of O₂ (oxidation) and C₂H₅OH, CD₃CH₂OH (H-atom abstraction) has shown that the √t law is not conditioned by the existence of regions characterized by different rate constants.     

CH3SO3裂解反应的机理和热力学性质

在G3XMP2//B3LYP/6-311+G(3df,2p)水平上对CH₃SO₃裂解反应的机理进行了研究, 获得了6 条通道(10 条路径), 并构建了其势能剖面. 同时采用单分子反应理论计算了各个通道在温度200-3000 K区间的速率常数. 研究结果表明, 在计算温度范围内, CH₃SO₃裂解反应的主产物为P1(CH₃+SO₃), 产物P2(CH₃O+SO₂)和P3(HCHO+HOSO)仅在温度大于3000 K时对总产物有贡献, 而产物P4(CHSO₂+H₂O), P5(CH₂SO₃+H)和P6(CHSO₃+H₂)贡献相对较少. 将裂解反应总的速率常数拟合为k_total=1.40×10¹²T^0.15exp(7831.58/T). 此外, 根据统计热力学原理, 预测了所有物种的生成焓(D_fH^Θ_{298 K}, D_fH_{0 K}), 熵(S^Θ_{298 K})和热容(C_p, 298-2000 K), 计算的结果与实验值较接近.     

Ab initio and RRKM calculations for the reaction channels of O(1D) + CH3CHF2

Ab initio and Rice–Ramsperger–Kassel–Marcus theories are carried out to study the potential energy surface and the energy‐dependent rate constants and branching ratios of the products for O(¹D) + CH₃CHF₂ reaction. Optimized geometries and vibrational frequencies have been obtained by MP2/6‐311G(d,p) method. The main products of the title reaction are CH₃CFO + HF, CH₂CFOH + HF, and CH₃ + CF₂OH at lower collision energy; and CH₃ + CF₂OH, CH₃CF₂ + OH are the main products at higher collision energy. CHF₂ + CH₂OH are the main products in the whole range of collision energy. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

The DFT study on mechanisms for NCO with C2H5 reaction

A detailed investigation has been performed at the QCISD(T)/6‐311++G(d,p)//B3LYP/6‐311+G(d,p) level for the reaction of NCO with C₂H₅ by constructing singlet and triplet potential energy surfaces (PES). The results show that the title reaction is more favorable on the singlet PES than on the triplet PES. On the singlet PES, the initial addition processes are barrierless and release lots of energy. The dominant channel occurs via the fragmentations of the initial adduct C₂H₅NCO and C₂H₅OCN to form C₂H₄ + HNCO and HOCN, respectively. With higher barrier heights, other products such as CH₄ + HNC + CO, CH₃CHNH + CO, CH₃CH + HNCO, and CH₃CN + H₂ + CO are less competitive. On the triplet PES, the entrance reactions surpass significant barriers; therefore, it could be negligible at the normal atmospheric condition. However, the most feasible channel on the triplet PES is the direct hydrogen abstraction channel to form CH₂CH₂ + HNCO. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

A Theoretical Investigation of the Decomposition Reactions of Ethyl Formate in the S1 and T1 States

This study investigates the decomposition reactions of ethyl formate in the S₁ and T₁ states. The dissociation channels leading to HCOOH + C₂H₄, CH₃CH₂O + HCO, CH₃CH₂OCO + H, and CH₃CH₂ + HCO₂ were studied. The major reactions of ethyl formate in the S₁ and T₁ states are isomerization to the biradical CH₂CH₂OC(OH)H and dissociation to CH₃CH₂O + HCO. All the stationary and intersection points were optimized at the CAS(10,8) level of theory with the 6‐31G(d,p) and 6‐311G(2df,2pd) basis sets. Single‐point CASPT3 energy was calculated for each of the stationary and intersection points. Microcanonical rate constants were also calculated for each of the reactions by using the RRKM theory.     

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     

An experimental and modeling study of ethanol oxidation kinetics in an atmospheric pressure flow reactor

Experimental profiles of stable species concentrations and temperature are reported for the flow reactor oxidation of ethanol at atmospheric pressure, initial temperatures near 1100 K and equivalence ratios of 0.61–1.24. Acetaldehyde, ethene, and methane appear in roughly equal concentrations as major intermediate species under these conditions. A detailed chemical mechanism is validated by comparison with the experimental species profiles. The importance of including all three isomeric forms of the C₂H₅O radical in such a mechanism is demonstrated. The primary source of ethene in ethanol oxidation is verified to be the decomposition of the C₂H₄OH radical. The agreement between the model and experiment at 1100 K is optimized when the branching ratio of the reactions of C₂H₅OH with OH and H is defined by (30% C₂H₄OH + 50% CH₃CHOH + 20% CH₃CH₂O) + XH. As in methanol oxidation, HO₂ chemistry is very important, while the H + O₂ chain branching reaction plays only a minor role until late in fuel decay, even at temperatures above 1100 K.     

In cation interactions with some organics: ab initio molecular orbital and density functional theory

Ab initio molecular orbital and density functional theory studies were undertaken to investigate the structural and energetic characteristics of complexes of In⁺ with several different organic molecules for the first time. HF, MP2, QCISD, and CCD levels of theory in ab initio MO as well as B3LYP, B3PW91 hybrid functionals in density functional theory were used. A valence TZ+P basis set with relativistic effective core potentials was used for the In atom while the 6-311++G(3d, 2p) basis set was utilized for all other atoms. Both closed-shell (H₂O, CH₄, CH₃OH, and C₆H₆) and open-shell (CH₃ and C₂H₃) molecules were considered for complexation with In⁺. In⁺ affinities of 21.5, 24.8, 28.6, 18.4, and 23.0 kcal/mol were obtained with the B3PW91 hybrid functional for H₂O, CH₃OH, C₆H₆, CH₃, and C₂H₃, respectively. The large values for the calculated affinities indicate the validity of our recent experimental detection of In⁺ ion attachment to some organic molecules.     

On the kinetic mechanism of reactions of hydroxyl radical with CHF2CH3 − n F n (n = 1–3)

The potential energy surfaces of the reactions CHF₂CH_3 − n F _n (n = 1–3) + OH were investigated by MPWB1K and BMC-CCSD (single-point) methods. Furthermore, with the aid of canonical variational transition state theory including the small-curvature tunneling correction, the rate constants of the title reactions were calculated over a wide temperature range of 220–1,500 K. Agreement between the CVT/SCT rate constants and the experimental values is good. Our results show that the order of rate constants is CHF₂CH₂F + OH > CHF₂CHF₂ + OH > CHF₂CF₃ + OH. For reaction CHF₂CH₂F + OH, the 1-H-abstraction channel dominates the reaction at the whole temperature, while 2-H-abstraction channel appears to be competitive with the increase of temperature.     

A detailed chemical kinetic model for high temperature ethanol oxidation

A detailed chemical kinetic model for ethanol oxidation has been developed and validated against a variety of experimental data sets. Laminar flame speed data (obtained from a constant volume bomb and counterflow twin‐flame), ignition delay data behind a reflected shock wave, and ethanol oxidation product profiles from a jet‐stirred and turbulent flow reactor were used in this computational study. Good agreement was found in modeling of the data sets obtained from the five different experimental systems. The computational results show that high temperature ethanol oxidation exhibits strong sensitivity to the fall‐off kinetics of ethanol decomposition, branching ratio selection for C₂H₅OH + OH ↔ Products, and reactions involving the hydroperoxyl (HO₂) radical. The multichanneled ethanol decomposition process is analyzed by RRKM/Master Equation theory, and the results are compared with those obtained from earlier studies. The ten‐parameter Troe form is used to define the C₂H₅OH(+M) ↔ CH₃ + CH₂OH(+M) rate expression as k^∞ = 5.94E23 T^−1.68 exp(−45880 K/T) (s⁻¹) k^o = 2.88E85 T^−18.9 exp(−55317 K/T) (cm³/mol/sec) F_cent = 0.5 exp(−T/200 K) + 0.5 exp(−T/890 K) + exp(−4600 K/T) and the C₂H₅OH(+M) ↔ C₂H₄ + H₂O(+M) rate expression as k^∞ = 2.79E13 T^0.09 exp(−33284 K/T) (s⁻¹) k^o = 2.57E83 T^−18.85 exp(−43509 K/T) (cm³/mol/sec) F _cent = 0.3 exp(−T/350 K) + 0.7 exp(−T/800 K) + exp(−3800 K/T) with an applied energy transfer per collision value of <ΔE_down> = 500 cm⁻¹. An empirical branching ratio estimation procedure is presented which determines the temperature dependent branching ratios of the three distinct sites of hydrogen abstraction from ethanol. The calculated branching ratios for C₂H₅OH + OH, C₂H₅OH + O, C₂H₅OH + H, and C₂H₅OH + CH₃ are compared to experimental data. © 1999 John Wiley & Sons, Inc. Int J Chem Kinet 31: 183–220, 1999     

Theoretical study of the reaction of ethynyl radical with acetonitrile

A detailed computational study has been performed on the mechanism and kinetics of the C₂H + CH₃CN reaction. The geometries were optimized at the BHandHLYP/6–311G(d, p) level. The single-point energies were calculated using the BMC-CCSD, MC-QCISD and QCISD(T)/6–311+G(2df, 2pd) methods. Five mechanisms were investigated, namely, direct hydrogen abstraction, C-addition/elimination, N-addition/elimination, C₂H–to–CN substitution and H-migration. The kinetics of the title reaction were studied using TST and multichannel RRKM methodologies over a wide range of temperatures (150–3,000 K) and pressures (10^?4–10⁴ torr). The total rate constants show positive temperature dependence and pressure independence. At lower temperatures, the C-addition step is the most feasible channel to produce CH₃ and HCCCN. At higher temperatures, the direct hydrogen abstraction path is the dominant channel leading to C₂H₂ and CH₂CN. The calculated overall rate constants are in good agreement with the experimental data.     

Study of the ethane oxidation reaction by the kinetic tracer method

The kinetics of ethane oxidation was studied at 320, 340, 353 and 380°C, mixture composition 2 C₂H₆ + 1 O₂, and total pressure 609 torr. It was found that at 320°C CH₂O and CH₃CHO were branching agents. A series of experiments was conducted on 2C₂H₆ + O₂ oxidation in the presence of 0.7% ¹⁴C-labeled ethylene. The ethylene oxide was found to form only from C₂H₄, formaldehyde formed from C₂H₄ and C₂H₆; and CH₃CHO, C₂H₅OH, and CH₃OH formed only from ethane. The formation rates of C₂H₄, C₂H₄O, and CH₂O were calculated by the kinetic tracer method. At 320°C the fraction of oxygen-containing products formed from C₂H₄ was 16–18%, and at 353 and 380°C it was 30–40%.