Predictive theory for hydrogen atom-hydrocarbon radical association kinetics

Procedures for accurately predicting the kinetics of hydrogen atom associations with hydrocarbon radicals are described and applied to a series of reactions. The approach is based on CASPT2/cc-pvdz evaluations of the orientation-dependent interaction energies within variable reaction coordinate transition state theory. One-dimensional corrections to the interaction energies are estimated from CAS+1+2/aug-cc-pvtz evaluations for the H + CH3 reaction, and a dynamical correction factor of 0.9 is applied. This corrected CASPT2 approach yields results that are within 10% of those obtained with the full CAS+1+2/aug-cc-pvtz potential for the H + CH3, H + C2H5, H + C2H3, and H + C2H reactions. New predictions are made for the H + iso-C3H7, H + tert-C4H9, H + C6H5, and H + C10H7 reactions. For the H + CH3 and H + C2H3 reactions, where the experimental values appear to be the most well-determined, theory and experiment essentially agree to within their error bars. For the other reactions, the agreement is reasonably satisfactory given the often large dispersion in the experimental results. For the reactions with saturated alkyl radicals, the theory predicts that each additional CH3 group increases the steric factor by approximately a factor of 2. In contrast, for the unsaturated radicals, the H + C6H5 and H + C10H7 high-pressure association rate coefficients are nearly identical to that for H + C2H3.     

Identification and chemistry of C4H3 and C4H5 isomers in fuel-rich flames

Quantitative identification of isomers of hydrocarbon radicals in flames is critical to understanding soot formation. Isomers of C4H3 and C4H5 in flames fueled by allene, propyne, cyclopentene, or benzene are identified by comparison of the observed photoionization efficiencies with theoretical simulations based on calculated ionization energies and Franck-Condon factors. The experiments combine molecular-beam mass spectrometry (MBMS) with photoionization by tunable vacuum-ultraviolet synchrotron radiation. The theoretical simulations employ the rovibrational properties obtained with B3LYP/6-311++G(d,p) density functional theory and electronic energies obtained from QCISD(T) ab initio calculations extrapolated to the complete basis set limit. For C4H3, the comparisons reveal the presence of the resonantly stabilized CH2CCCH isomer (i-C4H3). For C4H5, contributions from the CH2CHCCH2 (i-C4H5) and some combination of the CH3CCCH2 and CH3CHCCH isomers are evident. Quantitative concentration estimates for these species are made for allene, cyclopentene, and benzene flames. Because of low Franck-Condon factors, sensitivity to n-isomers of both C4H3 and C4H5 is limited. Adiabatic ionization energies, as obtained from fits of the theoretical predictions to the experimental photoionization efficiency curves, are within the error bars of the QCISD(T) calculations. For i-C4H3 and i-C4H5, these fitted adiabatic ionization energies are (8.06 +/- 0.05) eV and (7.60 +/- 0.05) eV, respectively. The good agreement between the fitted and theoretical ionization thresholds suggests that the corresponding theoretically predicted radical heats of formation (119.1, 76.3, 78.7, and 79.1 kcal/mol at 0 K for i-C4H3, i-C4H5, CH3CCCH2, and CH3CHCCH, respectively) are also quite accurate.     

Energetics, transition states, and intrinsic reaction coordinates for reactions associated with O(3P) processing of hydrocarbon materials

Electronic structure calculations based on multiconfiguration wave functions are used to investigate a set of archetypal reactions relevant to O(3P) processing of hydrocarbon molecules and surfaces. These include O(3P) reactions with methane and ethane to give OH plus methyl or ethyl radicals, O(3P) + ethane to give CH3O + CH3, and secondary reactions of the OH product radical with ethane and the ethyl radical. Geometry optimization is carried out with CASSCF/cc-pVTZ for all reactions, and with CASPT2/cc-pVTZ for O(3P) + methane/ethane. Single-point energy corrections are applied with CASPT2, CASPT3, and MRCI + Q with the cc-pVTZ and cc-pVQZ basis sets, and the energies extrapolated to the complete basis set limit (CBL). Where comparison of computed barriers and energies of reaction with experiment is possible, the agreement is good to excellent. The best agreement (within experimental error) is found for MRCI + Q/CBL applied to O(3P) + methane. For the other reactions, CASPT2/CBL and MRCI + Q/CBL predictions differ from experiment by 1-5 kcal/mol for 0 K enthalpies of reaction, and are within 1 kcal/mol of the best-estimate experimental range of 0 K barriers for O(3P) + ethane and OH + ethane. The accuracy of MRCI + Q/CBL is limited mainly by the quality of the active space. CASPT2/CBL barriers are consistently lower than MRCI + Q/CBL barriers with identical reference spaces.     

Crossed molecular beam studies of the reactions of allyl radicals, C3H5(X2A2), with methylacetylene (CH3CCH(X1A1)), allene (H2CCCH2(X1A1)), and their isotopomers

The chemical dynamics of the reaction of allyl radicals, C(3)H(5)(X(2)A(2)), with two C(3)H(4) isomers, methylacetylene (CH(3)CCH(X(1)A(1))) and allene (H(2)CCCH(2)(X(1)A(1))) together with their (partially) deuterated counterparts, were unraveled under single-collision conditions at collision energies of about 125 kJ mol(-1) utilizing a crossed molecular beam setup. The experiments indicate that the reactions are indirect via complex formation and proceed via an addition of the allyl radical with its terminal carbon atom to the terminal carbon atom of the allene and of methylacetylene (alpha-carbon atom) to form the intermediates H(2)CCHCH(2)CH(2)CCH(2) and H(2)CCHCH(2)CHCCH(3), respectively. The lifetimes of these intermediates are similar to their rotational periods but too short for a complete energy randomization to occur. Experiments with D4-allene and D4-methylacetylene verify explicitly that the allyl group stays intact: no hydrogen emission was observed but only the release of deuterium atoms from the perdeuterated reactants. Further isotopic substitution experiments with D3-methylacetylene combined with the nonstatistical nature of the reaction suggest that the intermediates decompose via hydrogen atom elimination to 1,3,5-hexatriene, H(2)CCHCH(2)CHCCH(2), and 1-hexen-4-yne, H(2)CCHCH(2)CCCH(3), respectively, via tight exit transition states located about 10-15 kJ mol(-1) above the separated products. The overall reactions were found to be endoergic by 98 +/- 4 kJ mol(-1) and have characteristic threshold energies to reaction between 105 and 110 kJ mol(-1). Implications of these findings to combustion and interstellar chemistry are discussed.     

Predictive theory for the combination kinetics of two alkyl radicals

An ab initio transition state theory based procedure for accurately predicting the combination kinetics of two alkyl radicals is described. This procedure employs direct evaluations of the orientation dependent interaction energies at the CASPT2/cc-pvdz level within variable reaction coordinate transition state theory (VRC-TST). One-dimensional corrections to these energies are obtained from CAS+1+2/aug-cc-pvtz calculations for CH3 + CH3 along its combination reaction path. Direct CAS+1+2/aug-cc-pvtz calculations demonstrate that, at least for the purpose of predicting the kinetics, the corrected CASPT2/cc-pvdz potential energy surface is an accurate approximation to the CAS+1+2/aug-cc-pvtz surface. Furthermore, direct trajectory simulations, performed at the B3LYP/6-31G* level, indicate that there is little local recrossing of the optimal VRC transition state dividing surface. The corrected CASPT2/cc-pvdz potential is employed in obtaining direct VRC-TST kinetic predictions for the self and cross combinations of methyl, ethyl, iso-propyl, and tert-butyl radicals. Comparisons with experiment suggest that the present dynamically corrected VRC-TST approach provides quantitatively accurate predictions for the capture rate. Each additional methyl substituent adjacent to a radical site is found to reduce the rate coefficient by about a factor of two. In each instance, the rate coefficients are predicted to decrease quite substantially with increasing temperature, with the more sterically hindered reactants having a more rapid decrease. The simple geometric mean rule, relating the capture rate for the cross reaction to those for the self-reactions, is in remarkably good agreement with the more detailed predictions. With suitable generalizations the present approach should be applicable to a wide array of radical-radical combination reactions.     

Association rate constants for reactions between resonance-stabilized radicals: C3H3 + C3H3, C3H3 + C3H5, and C3H5 + C3H5

Reactions between resonance-stabilized radicals play an important role in combustion chemistry. The theoretical prediction of rate coefficients and product distributions for such reactions is complicated by the fact that the initial complex-formation steps and some dissociation steps are barrierless. In this paper direct variable reaction coordinate transition state theory (VRC-TST) is used to predict accurately the association rate constants for the self and cross reactions of propargyl and allyl radicals. For each reaction, a set of multifaceted dividing surfaces is used to account for the multiple possible addition channels. Because of their resonant nature the geometric relaxation of the radicals is important. Here, the effect of this relaxation is explicitly calculated with the UB3LYP/cc-pvdz method for each mutual orientation encountered in the configurational integrals over the transition state dividing surfaces. The final energies are obtained from CASPT2/cc-pvdz calculations with all pi-orbitals in the active space. Evaluations along the minimum energy path suggest that basis set corrections are negligible. The VRC-TST approach was also used to calculate the association rate constant and the corresponding number of states for the C(6)H(5) + H --> C(6)H(6) exit channel of the C(3)H(3) + C(3)H(3) reaction, which is also barrierless. For this reaction, the interaction energies were evaluated with the CASPT2(2e,2o)/cc-pvdz method and a 1-D correction is included on the basis of CAS+1+2+QC/aug-cc-pvtz calculations for the CH(3) + H reference system. For the C(3)H(3) + C(3)H(3) reaction, the VRC-TST results for the energy and angular momentum resolved numbers of states in the entrance channels and in the C(6)H(5) + H exit channel are incorporated in a master equation simulation to determine the temperature and pressure dependence of the phenomenological rate coefficients. The rate constants for the C(3)H(3) + C(3)H(3) and C(3)H(5) + C(3)H(5) self-reactions compare favorably with the available experimental data. To our knowledge there are no experimental rate data for the C(3)H(3) + C(3)H(5) reaction.     

Identification of C5Hx isomers in fuel-rich flames by photoionization mass spectrometry and electronic structure calculations

The isomeric composition of C(5)H(x) (x = 2-6, 8) flame species is analyzed for rich flames fueled by allene, propyne, cyclopentene, or benzene. Different isomers are identified by their known ionization energies and/or by comparison of the observed photoionization efficiencies with theoretical simulations based on calculated ionization energies and Franck-Condon factors. The experiments combine flame-sampling molecular-beam mass spectrometry with photoionization by tunable vacuum-UV synchrotron radiation. The theoretical simulations employ the rovibrational properties obtained with B3LYP/6-311++G(d,p) density functional theory and electronic energies obtained from QCISD(T) electronic structure calculations extrapolated to the complete basis set limit. For C(5)H(3), the comparison reveals the presence of both the H(2)CCCCCH (i-C(5)H(3)) and the HCCCHCCH (n-C(5)H(3)) isomer. The simulations also suggest a modest amount of cyclo-CCHCHCCH-, which is consistent with a minor signal for C(5)H(2) that is apparently due to cyclo-CCHCCCH-. For C(5)H(4), contributions from the CH(2)CCCCH(2) (1,2,3,4-pentatetraene), CH(2)CCHCCH, and CH(3)CCCCH (1,3-pentadiyne) isomers are evident, as is some contribution from CHCCH(2)CCH (1,4-pentadiyne) in the cyclopentene and benzene flames. Signal at m/z = 65 originates mainly from the cyclopentadienyl radical. For C(5)H(6), contributions from cyclopentadiene, CH(3)CCCHCH(2), CH(3)CHCHCCH, and CH(2)CHCH(2)CCH are observed. No signal is observed for C(5)H(7) species. Cyclopentene, CH(2)CHCHCHCH(3) (1,3-pentadiene), CH(3)CCCH(2)CH(3) (2-pentyne), and CH(2)CHCH(2)CHCH(2) (1,4-pentadiene) contribute to the signal at m/z = 68. Newly derived ionization energies for i- and n-C(5)H(3) (8.20 +/- 0.05 and 8.31 +/- 0.05 eV, respectively), CH(2)CCHCCH (9.22 +/- 0.05 eV), and CH(2)CHCH(2)CCH (9.95 +/- 0.05 eV) are within the error bars of the QCISD(T) calculations. The combustion chemistry of the observed C(5)H(x) intermediates and the impact on flame chemistry models are discussed.     

Kinetics of the reaction of methyl radical with hydroxyl radical and methanol decomposition

The CH3 + OH bimolecular reaction and the dissociation of methanol are studied theoretically at conditions relevant to combustion chemistry. Kinetics for the CH3 + OH barrierless association reaction and for the H + CH2OH and H + CH3O product channels are determined in the high-pressure limit using variable reaction coordinate transition state theory and multireference electronic structure calculations to evaluate the fragment interaction energies. The CH3 + OH --> 3CH2 + H2O abstraction reaction and the H2 + HCOH and H2 + H2CO product channels feature localized dynamical bottlenecks and are treated using variational transition state theory and QCISD(T) energies extrapolated to the complete basis set limit. The 1CH2 + H2O product channel has two dynamical regimes, featuring both an inner saddle point and an outer barrierless region, and it is shown that a microcanonical two-state model is necessary to properly describe the association rate for this reaction over a broad temperature range. Experimental channel energies for the methanol system are reevaluated using the Active Thermochemical Tables (ATcT) approach. Pressure dependent, phenomenological rate coefficients for the CH3 + OH bimolecular reaction and for methanol decomposition are determined via master equation simulations. The predicted results agree well with experimental results, including those from a companion high-temperature shock tube determination for the decomposition of methanol.     

Thermochemistry and accurate quantum reaction rate calculations for H2/HD/D2 + CH3

Accurate quantum-mechanical results for thermodynamic data, cumulative reaction probabilities (for J = 0), thermal rate constants, and kinetic isotope effects for the three isotopic reactions H2 + CH3 --> CH4 + H, HD + CH3 --> CH4 + D, and D2 + CH3 --> CH(3)D + D are presented. The calculations are performed using flux correlation functions and the multiconfigurational time-dependent Hartree (MCTDH) method to propagate wave packets employing a Shephard interpolated potential energy surface based on high-level ab initio calculations. The calculated exothermicity for the H2 + CH3 --> CH4 + H reaction agrees to within 0.2 kcal/mol with experimentally deduced values. For the H2 + CH3 --> CH4 + H and D2 + CH3 --> CH(3)D + D reactions, experimental rate constants from several groups are available. In comparing to these, we typically find agreement to within a factor of 2 or better. The kinetic isotope effect for the rate of the H2 + CH3 --> CH4 + H reaction compared to those for the HD + CH3 --> CH4 + D and D2 + CH3 --> CH(3)D + D reactions agree with experimental results to within 25% for all data points. Transition state theory is found to predict the kinetic isotope effect accurately when the mass of the transferred atom is unchanged. On the other hand, if the mass of the transferred atom differs between the isotopic reactions, transition state theory fails in the low-temperature regime (T < 400 K), due to the neglect of the tunneling effect.     

Approach to the atmospheric chemistry of methyl nitrate and methylperoxy nitrite. Chemical mechanisms of their formation and decomposition reactions in the gas phase

Potential energy surfaces, minimum energy reaction paths, minima, transition states, reaction barriers, and conical intersections for the most important atmospheric reactions of methyl nitrate (CH(3)ONO(2)) and methylperoxy nitrite (C(3)HOONO) on the electronic ground state have been studied (i) with the second-order multiconfigurational perturbation theory (CASPT2) by computation of numerical energy gradients for stationary points and (ii) with the density functional theory (DFT). The proposed mechanism explains the conversion of unreactive alkyl peroxy radicals into alkoxy radicals: CH(3)O(2) + NO <=> CH(3)OONO <=> CH(3)O + NO(2) left arrow over right arrow CH(3)ONO(2). Additionally, several discrepancies found in the comparison of the results obtained from the two employed approaches are analyzed. CASPT2 predicts that all dissociation reactions into radicals occur without an extra exit energy barrier. In contrast, DFT finds transition states for the dissociations of cis- and trans-methylperoxy nitrite into CH(3)O + NO(2). Furthermore, multiconfigurational methods [CASPT2 and complete active space SCF (CAS-SCF)] predict the isomerization of CH3ONO2 to CH3OONO to occur in a two-step mechanism: (i) CH(3)ONO(2) --> CH(3)O + NO(2); and (ii) CH(3)O + NO(2) --> CH(3)OONO. The reason for this has to do with the coupling of the ground electronic state with the first excited state. Therefore, it is demonstrated that DFT methods based on single determinantal wave functions give an incorrect picture of the aforementioned reaction mechanisms.     

A study of the hydrogen abstraction reactions of C2H radical with CH3CN, C2H5CN, and C3H7CN by dual-level generalized transition state theory

The hydrogen abstraction reactions C2H + CH3CN --> products (R1), C2H + CH3CH2CN --> products (R2), and C2H + CH3CH2CH2CN --> products (R3) have been investigated by dual-level generalized transition state theory. Optimized geometries and frequencies of all the stationary points and extra points along the minimum-energy path (MEP) are performed at the BH&H-LYP and MP2 methods with the 6-311G(d, p) basis set, and the energy profiles are further refined at the MC-QCISD level of theory. The rate constants are evaluated using canonical variational transition state theory (CVT) with a small-curvature tunneling correction (SCT) over a wide temperature range 104-2000 K. The calculated CVT/SCT rate constants are in good agreement with the available experimental values. Our calculations show that for reaction R2, the alpha-hydrogen abstraction channel and beta-hydrogen abstraction channel are competitive over the whole temperature range. For reaction R3, the gamma-hydrogen abstraction channel is preferred at lower temperatures, while the contribution of beta-hydrogen abstraction will become more significant with a temperature increase. The branching ratio to the alpha-hydrogen abstraction channel is found negligible over the whole temperature range.     

Ab initio/Rice-Ramsperger-Kassel-Marcus study of the singlet C4H4 potential energy surface and of the reactions of C2(X1 Sigmag+) with C4H4(X1A1g) and C(1D) with C3H4 (allene and methylacetylene)

Ab initio modified Gaussian-2 G2M(RCC,MP2) calculations have been performed for various isomers and transition states on the singlet C4H4 potential energy surface. The computed relative energies and molecular parameters have then been used to calculate energy-dependent rate constants for different isomerization and dissociation processes in the C4H4 system employing Rice-Ramsperger-Kassel-Marcus theory and to predict branching ratios of possible products of the C2(1Sigmag+)+C2H4, C(1D)+H2CCCH2, and C(1D)+H3CCCH reactions under single-collision conditions. The results show that C2 adds to the double C=C bond of ethylene without a barrier to form carbenecyclopropane, which then isomerizes to butatriene by a formal C2 "insertion" into the C-C bond of the C2H4 fragment. Butatriene can rearrange to the other isomers of C4H4, including allenylcarbene, methylenecyclopropene, vinylacetylene, methylpropargylene, cyclobutadiene, tetrahedrane, methylcyclopropenylidene, and bicyclobutene. The major decomposition products of the chemically activated C4H4 molecule formed in the C2(1Sigmag+)+C2H4 reaction are calculated to be acetylene+vinylidene (48.6% at Ecol = 0) and 1-buten-3-yne-2-yl radical [i-C4H3(X2A'), H2C=C=C=CH*]+H (41.3%). As the collision energy increases from 0 to 10 kcal/mol, the relative yield of i-C4H3+H grows to 52.6% and that of C2H2+CCH2 decreases to 35.5%. For the C(1D)+allene reaction, the most important products are also i-C4H3+H (55.2%) and C2H2+CCH2 (30.1%), but for C(1D)+methylacetylene, which accesses a different region of the C4H4 singlet potential energy surface, the calculated product branching ratios differ significantly: 65%-69% for i-C4H3+H, 18%-14% for C2H2+CCH2, and approximately 8% for diacetylene+H2.     

Origin of the SN2 benzylic effect

The S N2 identity exchange reactions of the fluoride ion with benzyl fluoride and 10 para-substituted derivatives (RC6H 4CH 2F, R = CH3, OH, OCH 3, NH2, F, Cl, CCH, CN, COF, and NO2) have been investigated by both rigorous ab initio methods and carefully calibrated density functional theory. Groundbreaking focal-point computations were executed for the C6H5CH 2F + F (-) and C 6H 5CH2Cl + Cl (-) SN2 reactions at the highest possible levels of electronic structure theory, employing complete basis set (CBS) extrapolations of aug-cc-pV XZ (X = 2-5) Hartree-Fock and MP2 energies, and including higher-order electron correlation via CCSD/aug-cc-pVQZ and CCSD(T)/aug-cc-pVTZ coupled cluster wave functions. Strong linear dependences are found between the computed electrostatic potential at the reaction-center carbon atom and the effective SN2 activation energies within the series of para-substituted benzyl fluorides. An activation strain energy decomposition indicates that the SN2 reactivity of these benzylic compounds is governed by the intrinsic electrostatic interaction between the reacting fragments. The delocalization of nucleophilic charge into the aromatic ring in the SN2 transition states is quite limited and should not be considered the origin of benzylic acceleration of SN2 reactions. Our rigorous focal-point computations validate the benzylic effect by establishing SN2 barriers for (F (-), Cl (-)) identity exchange in (C6H5CH2F, C6H 5CH2Cl) that are lower than those of (CH3F, CH3Cl) by (3.8, 1.6) kcal mol (-1), in order.     

Ab initio methods for reactive potential surfaces

Case studies of ten reactions using a variety of standard electronic structure methods are presented. These case studies are used to illustrate the usefulness and shortcomings of these standard methods for various classes of reactions. Limited comparisons with experiment are made. The reactions studied include four radical-radical combinations, H + CH(3)--> CH(4), CH(3) + CH(3)--> C(2)H(6), H + HCO --> H(2)CO and CH(3) + HCO --> CH(3)CHO, three abstraction reactions, H + HO(2)--> H(2) + O(2), H + HCO --> H(2) + CO and CH(3) + HCO --> CH(4) + CO, a radical-molecule addition, H + HCCH --> C(2)H(3), and two molecular decompositions, H(2)CO --> H(2) + CO and CH(3)CHO --> CH(4) + CO. The electronic structure methods used are DFT, MP2, CCSD(T), QCISD(T), CASSCF, CASPT2, and CAS+1+2+QC.     

Dual-level direct dynamics studies for the reactions of dimethyl ether with hydrogen atom and methyl radical

The dual-level direct dynamics approach is employed to study the dynamics of the CH(3)OCH(3) + H (R1) and CH(3)OCH(3) + CH(3) (R2) reactions. Low-level calculations of the potential energy surface are carried out at the MP2/6-311+G(d,p) level of theory. High-level energetic information is obtained at the QCISD(T) level of theory with the 6-311+G(3df,3pd) basis set. The dynamics calculations are performed using variational transition state theory (VTST) with the interpolated single-point energies (ISPE) method, and small-curvature tunneling (SCT) is included. It is shown that the reaction of CH(3)OCH(3) with H (R1) may proceed much easier and with a lower barrier height than the reaction with CH(3) radical (R2). The calculated rate constants and activation energies are in good agreement with the experimental values. The calculated rate constants are fitted to k(R1) = 1.16 x 10(-19) T(3) exp(-1922/T) and k(R2) = 1.66 x 10(-28) T(5) exp(-3086/T) cm(3) mol(-1) s(-1) over a temperature range 207-2100 K. Furthermore, a small variational effect and large tunneling effect in the lower temperature range are found for the two reactions.     

Computational study of the ion-molecule reactions involving fluxional cations: CH4+ + H2--> CH5+ + H and isotope effect

Potential energy surfaces for the reactions of CH4+ with H2, HD, and D2 have been calculated using high-level ab initio methods, including coupled cluster theory, complete active space self-consistent field, and multireference configuration interaction. The energies are extrapolated to the complete basis set limit using the basis sets aug-cc-pVXZ (X = D, T, Q, 5, 6). The CH4+ + H2 reaction produces CH5+ and H exclusively. Three types of reaction mechanisms have been found, namely, complex-forming abstraction, scrambling, and S(N)2 displacement. The abstraction occurs via a very minor barrier and it is dominant. The other two mechanisms are negligible because of the significant barriers involved. Quantum phase space theory and variational transition state theory are used to calculate the rate coefficients as a function of temperatures in the range of 5-1000 K. The theoretical rate coefficients are compared with the available experimental data and the discrepancy is discussed. The significance of isotope effect, tunneling effect, and nuclear spin effect is investigated. The title reaction is predicted to be slightly exothermic with DeltaHr = -12.7 +/- 5.2 kJ/mol at 0 K.     

Criegee中间体气相反应热力学的G2理论计算

Criegee intermediate is believed to play an important role in the atmospheric chemistry. Because of its short life and the difficulty in experimental study, we carried out ah initio calculations on the thermochemistry of the Criegee involving reactions in this study. Thermochemistry data of reaction enthalpies and Gibbs free energies for four different stable structures of the Criegee intermediates (singlet CH2OO ①1 A1 in C2v, triplet CH2OO ②3B1 in C2v, singlet CH2OO ③1A' in Cs and triplet CH2OO ④ in C1 symmetry) involved in some of the gas-phase reactions were calculated at the standard Gaussian-2 [G2(MP2) and G2] and a modified G2, G2(fu1)[10],levels of theory. Relative energies among those Criegees and formic acid were compared. Chemical reactions include the formation of Criegees, re-arrangement from Criegee to formic acid, dissociations (producing CH2(3B1)+O2, CH2(1A1)+O2, CO2+H2, CO2+2H, CO+H2O, OH+HCO) and the reactions between Criegee and NO/H2O. Standard equilibrium constants for some reactions were investigated and may be obtained for all of the rest reactions involved in this study by the standard Gibbs free energies. It is shown that the formation of Criegee ①-④ by ethylene and ozone, the re-arrangement from any Criegee to formic acid, the dissociation in producing CO2+O2and CO+H2O and the reactions between any Criegee and NO/H2O are all favourable thermodynamicaly. The dissociation in forming CO2+2H and OH+HCO is less favourable. While the dissociation in forming carbene (either in 3B1 or 1A1 state) is not allowed by ΔrGm? values. Standard enthalpies of formation at 298 K for the four Criegees were predicted at the G2(ful) level of theory. Each value is the average value from ten of the above reactions and they are -4.3, 74.8,98.9 and 244.6 kJ mol-1 at the G2(ful) level for Criegee ① to Criegee ④, respectively. In addition, tile standard enthalpy of formation at 298 K for HOCH2OOH is further predicted to be -315.6 kJ mol-1 at the G2(MP2) level.     

Mechanistic and kinetic study of the O + CH2OH reaction

The mechanism for the O + CH2OH reaction was investigated by various ab initio quantum chemistry methods. For the chemical activation mechanism, that is, the addition/elimination path, the couple-cluster methods including CCSD and CCSD(T) were employed with the cc-pVXZ (X = D, T, Q, 5) basis sets. For the abstraction channels, multireference methods including CASSCF, CASPT2, and MRCISD were used with the cc-pVDZ and cc-pVTZ basis sets. It has been shown that the production of H + HCOOH is the major channel in the chemical activation mechanism. The minor channels include HCO + H2O and OH + CH2O. The hydrogen abstraction by an O atom from the CH2OH radical produces either OH + CH2O or OH + HCOH. Moreover, the two abstraction reactions are essentially barrierless processes. The rate constants for the association of O with CH2OH have been calculated using the flexible transition state theory. A weak negative temperature dependence of the rate constants is found in the range 250-1000 K. Furthermore, it is estimated that the abstraction processes also play an important role in the O + CH2OH reaction. Additionally, the falloff behavior for the OCH2OH --> H + HCOOH reaction has been investigated. The present theoretical results are compared to the experimental measurements to understand the mechanism and kinetic behavior of the O + CH2OH reaction and the unimolecular reaction of the OCH2OH radical.     

Experimental and theoretical studies of the reactions Y (a2D) + H2CO and Y (a2D) + CH3CHO

The reactions of ground state Y (a(2)D) with H(2)CO and CH(3)CHO were studied at a range of collision energies in crossed molecular beams. For reaction with H(2)CO, three product channels were observed: formation of YH(2) + CO, YCO + H(2), and YHCO + H. Reaction with CH(3)CHO led to three analogous product channels involving formation of HYCH(3) + CO, YCH(2)CO + H(2), and YCH(3)CO + H. The calculated CCSD(T) energetics and DFT geometries for key intermediates in both reactions, together with RRKM theory, are used to calculate a priori the branching ratios between various product channels. These calculated values are compared to those obtained experimentally.     

Ab initio rate constants from hyperspherical quantum scattering: application to H+C2H6 and H+CH3OH

The dynamics and kinetics of the abstraction reactions of H atoms with ethane and methanol have been studied using a quantum mechanical procedure. Bonds being broken and formed are treated with explicit hyperspherical quantum dynamics. The ab initio potential energy surfaces for these reactions have been developed from a minimal number of grid points (average of 48 points) and are given by analytical functionals. All the degrees of freedom except the breaking and forming bonds are optimized using the second order perturbation theory method with a correlation consistent polarized valence triple zeta basis set. Single point energies are calculated on the optimized geometries with the coupled cluster theory and the same basis set. The reaction of H with C2H6 is endothermic by 1.5 kcal/mol and has a vibrationally adiabatic barrier of 12 kcal/mol. The reaction of H with CH3OH presents two reactive channels: the methoxy and the hydroxymethyl channels. The former is endothermic by 0.24 kcal/mol and has a vibrationally adiabatic barrier of 13.29 kcal/mol, the latter reaction is exothermic by 7.87 kcal/mol and has a vibrationally adiabatic barrier of 8.56 kcal/mol. We report state-to-state and state-selected cross sections together with state-to-state rate constants for the title reactions. Thermal rate constants for these reactions exhibit large quantum tunneling effects when compared to conventional transition state theory results. For H+CH3OH, it is found that the CH2OH product is the dominant channel, and that the CH3O channel contributes just 2% at 500 K. For both reactions, rate constants are in good agreement with some measurements.