Multiple-Well, multiple-path unimolecular reaction systems. II. 2-methylhexyl free radicals

Exploring the OH-initiated reactions of styrene in the atmosphere and the role of van der Waals complex

Reacting with OH provides a major sink for styrene in the atmosphere, with three possible pathways including OH-addition, H-abstraction and addition-dissociation reactions. However, the total rate coefficients of styrene + OH were measured as 1.2-6.2 × 10^-11 cm³ molecule^-1 s^-1 under atmospheric conditions, varying by a maximum factor of 5. On the other hand, only one theoretical work reported this rate coefficient as 19.1 × 10^-11 cm³ molecule^-1 s^-1, which exhibits up to 16 times that measured in laboratory studies. In the present study, the reaction kinetics of styrene + OH was extensively studied with high-level quantum chemical methods combined with RRKM/master equation simulations. In particular, we carried out theoretical treatments for the formation of pre-reaction Van der Waals complexes of styrene + OH, and examined their influence on the reaction kinetics. The total rate coefficient for styrene + OH is calculated to be 1.7 × 10^-11 cm³ molecule^-1 s^-1 at 300 K, 1 atm. The main products are addβ (88.2%), add5 (6.9%), addα (1.9%) and add3 (1.7%). Using our computed rate coefficient and the global atmospheric hydroxyl radical concentration (2 × 10⁶ radicals per cm³), the lifetime of styrene in the atmosphere is estimated at 8.0 h. The degradation of styrene might be negligible for the formation of ozone in the atmosphere based upon the photochemical ozone creation potentials calculation. The computed product yields indicate that addβ via subsequent reactions could significantly produce formaldehyde and benzaldehyde that were observed in previous experimental studies on styrene oxidation, and contribute to the formation of secondary organic aerosols.

Energy Dependence of Ensemble-Averaged Energy Transfer Moments and Its Effect on Competing Decomposition Reactions

We carried out a direct dynamics study on the internal-energy dependence of the ensemble-averaged energy transfer moments of the isobutyl radical in collisions with N₂ bath gas. We find a linear dependence of the downward moment ⟨ΔE_d⟩ and the root-mean-square moment ⟨ΔE2⟩ on the initial internal energy, but the upward moment ⟨ΔE_u⟩ is found to be independent of the molecule's internal energy. We improved the exponential-down relaxation model by including a linear dependence of ⟨ΔE_d⟩ on the initial energy, and we used the improved treatment in the 1D master equation for isobutyl radical decomposition reactions and for a model of competitive reactions with a larger difference in barrier heights. We calculated phenomenological rate constants and branching ratios from chemically significant eigenmodes of the master equation and showed that the energy dependence of ⟨ΔE_d⟩ has a greater influence on channels with higher barriers in competitive reactions. Rate constants and branching ratios from master equation calculations indicate that for a given temperature and pressure, there is a constant ⟨ΔE_d⟩ that can reproduce results obtained with an E-dependent ⟨ΔE_d⟩. But a constant ⟨ΔE_d⟩ cannot do this for all temperatures and pressures, with larger differences when the barriers for the competing channels differ more. We conclude that when the branching ratio of competitive reactions is sensitive to pressure, including the energy dependence of ⟨ΔE_d⟩ in master equation simulations can make a significant difference in the results.

Two-Dimensional Master Equation Modeling of Some Multichannel Unimolecular Reactions

Multichannel thermal decomposition reactions of n-propyl radicals, 1-pentyl radicals, and toluene are investigated by solving a two-dimensional master equation formulated as a function of total energy (E) and angular momentum (J). The primary aim of this study is to elucidate the role of angular momentum in the kinetics of multichannel unimolecular reactions. The collisional transition processes of the reactants colliding with argon are characterized based on the classical trajectory calculations and implemented in the master equation. The rate constants calculated by using the two-dimensional master equation are compared with those of one-dimensional master equations. The consequence of the explicit treatment of angular momentum depends on the J dependence of the microscopic rate constants and is particularly emphasized in the thermal decomposition of toluene, for which the C-H and C-C bond fission channels are considered. The centrifugal effect is insignificant in the energetically favored C-H bond fission but is substantial in the energetically higher C-C bond fission, which causes rotational channel switching of the microscopic rate constants. The proper treatment of the J-dependent channel coupling effect, weak collisional transfer of J, and initial-J-dependent collisional energy transfer are found to be essential for predicting the branching fractions at low pressures.

Low-Pressure Limit of Competitive Unimolecular Reactions

Barker and Ortiz found unusual falloff effects in the flux coefficients of the competitive unimolecular reactions of 2-methylhexyl radicals, and they concluded that this might have important effects on the rate constants of reactions with higher thresholds. To study this effect, we carried out master equation calculations of the same reaction system to learn whether this effect shows up in measurable rate constants, and the answer is yes. We also studied specially designed mechanisms to reveal that the various reactive pathways connecting the reagents can have a large effect on the rate constants, causing them to be quite different than if the reactions proceeded independently, and that reactions with significantly higher barriers may nevertheless have larger rate constants. This provides a new perspective for interpreting and predicting the kinetics of competitive unimolecular reactions.

Atmospheric oxidation mechanism and kinetics of isoprene initiated by chlorine radicals: A computational study

The reaction with chlorine radicals (·Cl) has been considered to be one of indispensable sinks for isoprene. However, the mechanism of ·Cl initiated isoprene reaction was not fully understood. Herein, the reaction of isoprene with ·Cl, and ensuing reactions of the resulting isoprene relevant radicals were investigated by combined quantum chemistry calculations and kinetics modeling. The results indicate that ·Cl addition to two terminal C-atoms of two double bonds of isoprene, forming IM_1-1 and IM_1-4, are more favorable than H-abstractions from isoprene. Interestingly, the predicted reaction rate constant for the direct H-abstraction pathway is much lower than that of the indirect one, clarifying a direct H-abstraction mechanism for previously experimental observation. The IM_1-1 and IM_1-4 have distinct fate in their subsequent transformation. The reaction of IM_1-1 ends after the one-time O₂ addition. However, IM_1-4 can follow auto-oxidation mechanism with two times O₂ addition to finally form highly oxidized multi-functional molecules (HOMs), C₅H₇ClO₃ and ·OH. More importantly, the estimated contribution of ·Cl on HOMs (monomer only) formation from isoprene is lower than that of ·OH in addition pathway, implying overall HOMs yield from atmospheric isoprene oxidation could be overestimated if the role of ·Cl in transforming isoprene is ignored.

Atmospheric Oxidation of Piperazine Initiated by ·Cl: Unexpected High Nitrosamine Yield

Chlorine radicals (·Cl) initiated amine oxidation plays an important role for the formation of carcinogenic nitrosamine in the atmosphere. Piperazine (PZ) is considered as a potential atmospheric pollutant since it is an alternative solvent to monoethanolamine (MEA), a benchmark solvent in a leading CO₂ capture technology. Here, we employed quantum chemical methods and kinetics modeling to investigate ·Cl-initiated atmospheric oxidation of PZ, particularly concerning the potential of PZ to form nitrosamine compared to MEA. Results showed that the ·Cl-initiated PZ reaction exclusively leads to N-center radicals (PZ-N) that mainly react with NO to produce nitrosamine in their further reaction with O₂/NO. Together with the PZ + ·OH reaction, the PZ-N yield from PZ oxidation is still lower than that of the corresponding MEA reactions. However, the nitrosamine yield of PZ is higher than the reported value for MEA when [NO] is <5 ppb, a concentration commonly encountered in a polluted urban atmosphere. The unexpected high nitrosamine yield from PZ compared to MEA results from a more favorable reaction of N-center radicals with NO compared to O₂. These findings show that the yield of N-center radicals cannot directly be used as a metric for the yield of the corresponding carcinogenic nitrosamine.

β-Bond Scission and the Yields of H and CH3 in the Decomposition of Isobutyl Radicals

The relative rates of C-C and C-H β-scission reactions of isobutyl radicals (2-methylprop-1-yl, C₄H₉) were investigated with shock tube experiments at temperatures of (950 to 1250) K and pressures of (200 to 400) kPa. We produced isobutyl radicals from the decomposition of dilute mixtures of isopentylbenzene and observed the stable decomposition products, propene and isobutene. These alkenes are characteristic of C-C and C-H bond scission, respectively. Propene was the main product, approximately 30 times more abundant than isobutene, indicating that C-C β-scission is the primary pathway. Uncertainty in the ratio of [isobutene]/[propene] from isobutyl decomposition is mainly due to a small amount of side chemistry, which we account for using a kinetics model based on JetSurF 2.0. Our data are well-described after adding chemistry specific to our system and adjusting some rate constants. We compare our data to other commonly used kinetics models: JetSurF 2.0, AramcoMech 2.0, and multiple models from Lawrence Livermore National Laboratory (LLNL). With the kinetics model, we have determined an upper limit of 3.0% on the branching fraction for C-H β-scission in the isobutyl radical for the temperatures and pressures of our experiments. While this agrees with previous high quality experimental results, many combustion kinetics models assume C-H branching values above this upper limit, possibly leading to large systematic inaccuracies in model predictions. Some kinetics models additionally assume contributions from 1,2-H shift reactions-which for isobutyl would produce the same products as C-H β-scission-and our upper limit includes possible involvement of such reactions. We suggest kinetics models should be updated to better reflect current experimental measurements.

Atmospheric oxidation of halogenated aromatics: comparative analysis of reaction mechanisms and reaction kinetics

Atmospheric transport is the major route for global distribution of semi-volatile compounds such as halogenated aromatics as well as their major exposure route for humans. Their major atmospheric removal process is oxidation by hydroxyl radicals. There is very little information on the reaction mechanism or reaction-path dynamics of atmospheric degradation of halogenated benzenes. Furthermore, the measured reaction rate constants are missing for the range of environmentally relevant temperatures, i.e. 230-330 K. A series of recent theoretical studies have provided those valuable missing information for fluorobenzene, chlorobenzene, hexafluorobenzene and hexachlorobenzene. Their comparative analysis has provided additional and more general insight into the mechanism of those important tropospheric degradation processes as well as into the mobility, transport and atmospheric fate of halogenated aromatic systems. It was demonstrated for the first time that the addition of hydroxyl radicals to monohalogenated as well as to perhalogenated benzenes proceeds indirectly, via a prereaction complex and its formation and dynamics have been characterized including the respective transition-state. However, in fluorobenzene and chlorobenzene reactions hydroxyl radical hydrogen is pointing approximately to the center of the aromatic ring while in the case of hexafluorobenzene and hexachlorobenzene, unexpectedly, the oxygen is directed towards the center of the aromatic ring. The reliable rate constants are now available for all environmentally relevant temperatures for the tropospheric oxidation of fluorobenzene, chlorobenzene, hexafluorobenzene and hexachlorobenzene while pentachlorophenol, a well-known organic micropollutant, seems to be a major stable product of tropospheric oxidation of hexachlorobenzene. Their calculated tropospheric lifetimes show that fluorobenzene and chlorobenzene are easily removed from the atmosphere and do not have long-range transport potential while hexafluorobenzene seems to be a potential POP chemical and hexachlorobenzene is clearly a typical persistent organic pollutant.

Atmospheric chemical reaction mechanism and kinetics of 1,2-bis(2,4,6-tribromophenoxy)ethane initiated by OH radical: a computational study

The mechanism and kinetics of OH-initiated oxidation of BTBPE, an alternative of PBDEs, were investigated.

Variational transition state theory: theoretical framework and recent developments

This article reviews the fundamentals of variational transition state theory (VTST), its recent theoretical development, and some modern applications.

Atmospheric chemical reactions of alternatives of polybrominated diphenyl ethers initiated by OH: A case study on triphenyl phosphate

Many studies have been performed to evaluate the environmental risk caused by alternative flame retardants (AFRs) of polybrominated diphenyl ethers due to their ubiquitous occurrence in the environment. However, as an indispensable component of the environmental risk assessment, the information on atmospheric fate of AFRs is limited although some AFRs have been frequently and highly detected in the atmosphere. Here, a combined quantum chemical method and kinetics modeling were used to investigate atmospheric transformation mechanism and kinetics of AFRs initiated by OH in the presence of O2, taking triphenyl phosphate (TPhP) as a case. Results show that the pathway involving initial OH addition to phenyl of TPhP to form TPhP-OH adduct, and subsequent reaction of the TPhP-OH adduct with O2 to finally form phenol phosphate, is the most favorable for the titled reaction. The calculated overall reaction rate constant is 1.6×10(-12)cm(3) molecule(-1)s(-1), translating 7.6days atmospheric lifetime of TPhP. This clarifies that gaseous TPhP has atmospheric persistence. In addition, it was found that ice surface, as a case of ubiquitous water in the atmosphere, has little effect on the kinetics of the rate-determining step in the OH-initiated TPhP reaction.

Atmospheric oxidation of hexachlorobenzene: New global source of pentachlorophenol

Hexachlorobenzene is highly persistent, bioaccumulative, toxic and globally distributed, a model persistent organic pollutant. The major atmospheric removal process for hexachlorobenzene is its oxidation by hydroxyl radicals. Unfortunately, there is no information on the reaction mechanism of this important atmospheric process and the respective degradation rates were measured in a narrow temperature range not of environmental relevance. Thus, the geometries and energies of all stationary points significant for the atmospheric oxidation of hexachlorobenzene are optimized using MP2/6-311G(d,p) method. Furthermore, the single point energies were calculated with G3 method on the optimized minima and transition-states. It was demonstrated for the first time that the addition of hydroxyl radicals to hexachlorobenzene proceeds indirectly, via a prereaction complex. In the prereaction complex the hydroxyl radical is almost perpendicular to the aromatic ring while oxygen is pointing to its center. In contrast, in the transition state it is nearly parallel with the aromatic ring. The reliable rate constants are calculated for the first time for the atmospheric oxidation of hexachlorobenzene for all environmentally relevant temperatures. It was also demonstrated for the first time that pentachlorophenol is the major stable product in the addition of hydroxyl radicals to hexachlorobenzene and that atmosphere seems to be a new global secondary source of pentachlorophenol.

Quantum Chemical Study on ·Cl-Initiated Atmospheric Degradation of Monoethanolamine

Recent findings on the formation of ·Cl in continental urban areas necessitate the consideration of ·Cl initiated degradation when assessing the fate of volatile organic pollutants. Monoethanolamine (MEA) is considered as a potential atmospheric pollutant since it is a benchmark and widely utilized solvent in a leading CO2 capture technology. Especially, ·Cl may have specific interactions with the N atom of MEA, which could make the MEA + ·Cl reaction have different pathways and products from those of the MEA + ·OH reaction. Hence, ·Cl initiated reactions with MEA were investigated by a quantum chemical method [CCSD(T)/aug-cc-pVTZ//MP2/6-31+G(3df,2p)] and kinetics modeling. Results show that the overall rate constant for ·Cl initiated H-abstraction of MEA is 5 times faster than that initiated by ·OH, and the tropospheric lifetimes of MEA will be overestimated by 6-46% when assuming that [·Cl]/[·OH] = 1-10% if the role of ·Cl is ignored. The MEA + ·Cl reaction exclusively produces MEA-N that finally transforms into several products including mutagenic nitramine and carcinogenic nitrosamine via further reactions with O2/NOx, and the contribution of ·Cl to their formation is about 25-250% of that of ·OH. Thus, it is necessary to consider ·Cl initiated tropospheric degradation of MEA for its risk assessment.

The pyrolysis of 2-methylfuran: a quantum chemical, statistical rate theory and kinetic modelling study

Due to the rapidly growing interest in the use of biomass derived furanic compounds as potential platform chemicals and fossil fuel replacements, there is a simultaneous need to understand the pyrolysis and combustion properties of such molecules. To this end, the potential energy surfaces for the pyrolysis relevant reactions of the biofuel candidate 2-methylfuran have been characterized using quantum chemical methods (CBS-QB3, CBS-APNO and G3). Canonical transition state theory is employed to determine the high-pressure limiting kinetics, k(T), of elementary reactions. Rice-Ramsperger-Kassel-Marcus theory with an energy grained master equation is used to compute pressure-dependent rate constants, k(T,p), and product branching fractions for the multiple-well, multiple-channel reaction pathways which typify the pyrolysis reactions of the title species. The unimolecular decomposition of 2-methylfuran is shown to proceed via hydrogen atom transfer reactions through singlet carbene intermediates which readily undergo ring opening to form collisionally stabilised acyclic C5H6O isomers before further decomposition to C1-C4 species. Rate constants for abstraction by the hydrogen atom and methyl radical are reported, with abstraction from the alkyl side chain calculated to dominate. The fate of the primary abstraction product, 2-furanylmethyl radical, is shown to be thermal decomposition to the n-butadienyl radical and carbon monoxide through a series of ring opening and hydrogen atom transfer reactions. The dominant bimolecular products of hydrogen atom addition reactions are found to be furan and methyl radical, 1-butene-1-yl radical and carbon monoxide and vinyl ketene and methyl radical. A kinetic mechanism is assembled with computer simulations in good agreement with shock tube speciation profiles taken from the literature. The kinetic mechanism developed herein can be used in future chemical kinetic modelling studies on the pyrolysis and oxidation of 2-methylfuran, or the larger molecular structures for which it is a known pyrolysis/combustion intermediate (e.g. cellulose, coals, 2,5-dimethylfuran).

Mechanisms and reaction-path dynamics of hydroxyl radical reactions with aromatic hydrocarbons: the case of chlorobenzene

All geometries and energies significant for the first step of tropospheric degradation of chlorobenzene were characterized using the MP2/6-31+G(d,p) and G3 methods. A pre-reaction complex for the addition of OH radical to chlorobenzene was found and the associated transition state was determined for the first time. The reaction path for the association of OH radical and chlorobenzene into the pre-reaction complex was extrapolated from the selected low frequency normal mode of pre-reaction complex. The reaction rate constant for addition of OH radical to chlorobenzene was determined for the temperature range 230-330K, using RRKM theory and G3 energies. The calculated rate constants are in agreement with the experimental results. Regioselectivity was also determined for the title reaction from the ratio of respective reaction rates and our results are in very good agreement with the experimental results, which show the dominance of the ortho and para channels as well as a negligible contribution by the ipso channel.

Atmospheric fate of methacrolein. 2. Formation of lactone and implications for organic aerosol production

We investigate the oxidation of methacryloylperoxy nitrate (MPAN) and methacrylicperoxy acid (MPAA) by the hydroxyl radical (OH) theoretically, using both density functional theory [B3LYP] and explicitly correlated coupled cluster theory [CCSD(T)-F12]. These two compounds are produced following the abstraction of a hydrogen atom from methacrolein (MACR) by the OH radical. We use a RRKM master equation analysis to estimate that the oxidation of MPAN leads to formation of hydroxymethyl-methyl-α-lactone (HMML) in high yield. HMML production follows a low potential energy path from both MPAN and MPAA following addition of OH (via elimination of the NO(3) and OH from MPAN and MPAA, respectively). We suggest that the subsequent heterogeneous phase chemistry of HMML may be the route to formation of 2-methylglyceric acid, a common component of organic aerosol produced in the oxidation of methacrolein. Oxidation of acrolein, a photo-oxidation product from 1,3-butadiene, is found to follow a similar route generating hydroxymethyl-α-lactone (HML).

Quantum chemical and master equation studies of the methyl vinyl carbonyl oxides formed in isoprene ozonolysis

Methyl vinyl carbonyl oxide is an important intermediate in the reaction of isoprene and ozone and may be responsible for most of the (*)OH formed in isoprene ozonolysis. We use CBS-QB3 calculations and RRKM/master equation simulations to characterize all the pathways leading to the formation of this species, all the interconversions among its four possible conformers, and all of its irreversible isomerizations. Our calculations, like previous studies, predict (*)OH yields consistent with experiment if thermalized syn-methyl carbonyl oxides form (*)OH quantitatively. Natural bond order analysis reveals that the vinyl group weakens the C=O bond of the carbonyl oxide, making rotation about this bond accessible to this chemically activated intermediate. The vinyl group also allows one conformer of the carbonyl oxide to undergo electrocyclization to form a dioxole, a species not previously considered in the literature. Dioxole formation, which has a CBS-QB3 reaction barrier of 13.9 kcal/mol, is predicted to be favored over vinyl hydroperoxide formation, dioxirane formation, and collisional stabilization. Our calculations also predict that two dioxole derivatives, 1,2-epoxy-3-butanone and 3-oxobutanal, should be major products of isoprene ozonolysis.

Master equation simulations of competing unimolecular and bimolecular reactions: application to OH production in the reaction of acetyl radical with O2

Master equation calculations were carried out to simulate the production of hydroxyl free radicals initiated by the reaction of acetyl free radicals (CH3(C=O).) with molecular oxygen. In particular, the competition between the unimolecular reactions and bimolecular reactions of vibrationally excited intermediates was modeled by using a single master equation. The vibrationally excited intermediates (isomers of acetylperoxyl radicals) result from the initial reaction of acetyl free radical with O2. The bimolecular reactions were modeled using a novel pseudo-first-order microcanonical rate constant approach. Stationary points on the multi-well, multi-channel potential energy surface (PES) were calculated at the DFT(B3LYP)/6-311G(2df,p) level of theory. Some additional calculations were carried out at the CASPT2(7,5)/6-31G(d) level of theory to investigate barrierless reactions and other features of the PES. The master equation simulations are in excellent agreement with the experimental OH yields measured in N2 or He buffer gas near 300 K, but they do not explain a recent report that the OH yields are independent of pressure in nearly pure O2 buffer gas.

Evaluation of the Kinetic Data for Intramolecular 1,x-Hydrogen Shifts in Alkyl Radicals and Structure/Reactivity Predictions from the Carbocyclic Model for the Transition State

Experimental and computational kinetic data for the intramolecular 1,x-hydrogen shift in alkyl radicals are compiled in Arrhenius format for x = 2-5. Significant experimental disparity remains, especially for x = 2 and 3. Experimental data for radicals with tert centers or bearing spectator substituents are lacking for all x, and none exist for x = 6. The common use of the strain energy of the unsubstituted (x+1)-carbocycle to coarsely model the activation energy for the 1,x-shift is extended to explore more subtle differences in progressively methyl-substituted systems by use of molecular mechanics estimates of differences in strain between radicals and carbocycles. For x = 5 and 6, a sterically driven increase in E is predicted for shifts in the tert --> tert class that apparently runs counter to the behavior of bimolecular hydrogen transfers. In contrast, a sterically driven decrease in E is predicted to result from spectator methyl groups for the prim --> prim reaction class for all x. There is no experimental basis to test these predictions; fragmentary computational evidence lends some support to the second but is ambiguous concerning the first. Possible deficiencies in the use of carbocycles as transition state models are discussed.

Test of a chemical timing method for measuring absolute vibrational relaxation rate constants for S1 p-difluorobenzene

A chemical timing (CT) method for measuring absolute rate constants for collisional vibrational relaxation has been tested for the 5(1) state of S(1) p-difluorobenzene (pDFB) where an alternative method exists to provide benchmark values. The CT method was originally developed to treat vibrational energy transfer (VET) in large molecules excited to high vibrational levels where the intramolecular vibrational redistribution (IVR) resulting from large vibrational state densities completely eliminates vibrational structure in the emission spectrum. Here we apply the same method to a low-lying state (5(1) with epsilon(vib) = 818 cm(-1)) located in the low-density region of the vibrational manifold where IVR plays no role. For high vibrational levels, the chemical timing method involves addition of high O(2) pressures (kTorr) to a low-pressure pDFB sample, introducing vibrational structure in the fluorescence spectrum. Response of this spectrum to vibrational relaxation by Ar is then examined. For levels such as 5(1), the fully structured fluorescence spectrum allows the rate constant for single-collision VET into the surrounding vibrational field to be measured directly without the presence of O(2). The measurements of 5(1) VET have been repeated with various O(2) pressures (kTorr) for comparison with the O(2)-free benchmark. In the presence of O(2), the rate constant for VET by Ar is (4.0 +/- 0.5) x 10(6) Torr(-1) s(-1) and independent of high O(2) pressure variations. The rate constant as found by the standard O(2)-free method is (3.6 +/- 0.4) x 10(6) Torr(-1) s(-1). This comparison suggests that the chemical timing method is capable of providing a reasonably accurate measure of the VET rate constant for high vibrational levels provided that details of the kinetics are known.

Quantum Chemical and Master Equation Simulations of the Oxidation and Isomerization of Vinoxy Radicals

The vinoxy radical, a common intermediate in gas-phase alkene ozonolysis, reacts with O2 to form a chemically activated alpha-oxoperoxy species. We report CBS-QB3 energetics for O2 addition to the parent (*CH2CHO, 1a), 1-methylvinoxy (*CH2COCH3, 1b), and 2-methylvinoxy (CH3*CHCHO, 1c) radicals. CBS-QB3 predictions for peroxy radical formation agree with experimental data, while the G2 method systematically overestimates peroxy radical stability. RRKM/master equation simulations based on CBS-QB3 data are used to estimate the competition between prompt isomerization and thermalization for the peroxy radicals derived from 1a, 1b, and 1c. The lowest energy isomerization pathway for radicals 4a and 4c (derived from 1a and 1c, respectively) is a 1,4-shift of the acyl hydrogen requiring 19-20 kcal/mol. The resulting hydroperoxyacyl radical decomposes quantitatively to form *OH. The lowest energy isomerization pathway for radical 4b (derived from 1b) is a 1,5-shift of a methyl hydrogen requiring 26 kcal/mol. About 25% of 4a, but only approximately 5% of 4c, isomerizes promptly at 1 atm pressure. Isomerization of 4b is negligible at all pressures studied.

Product Study of the Reaction of CH3 with OH Radicals at Low Pressures and Temperatures of 300 and 612 K†

The product distribution for the title reaction was studied using our time-of-flight mass spectrometer (TOFMS) connected to a tubular flow reactor. The methyl and hydroxyl radicals were produced by an excimer laser pulse (lambda = 193 nm) photolyzing acetone and nitrous oxide in the presence of excess water or hydrogen. Helium was used as the bath gas; the total density was held constant at 1.2 x 10(17) cm(-3). At 300 K the observations were consistent with singlet methylene ((1)CH(2)) and water as the main product channel with a small contribution of methanol. In contrast, at about 610 K three channels-formaldehyde isomers and methanol in addition to (1)CH(2) + H(2)O-are formed with similar yields. When acetone-d(6) was used, the production of both CHDO and CD(2)O was observed, indicating that two different formaldehyde-producing channels are operating simultaneously. These experimental results are compared with RRKM and master equation calculations on the basis of the properties of the methanol potential energy surface from a recent ab initio study.

Pressure dependence for the CO quantum yield in the photolysis of acetone at 248 nm: A combined experimental and theoretical study

The quantum yield of CO in the laser pulse photolysis of acetone at 248 nm and at 298 K in the pressure range 20-900 mbar (N2) has been measured directly using quantitative infrared diode laser absorption of CO. It is found that the quantum yield of CO shows a significant dependence on total pressure with Phi(CO) decreasing with pressure from around 0.45 at 20 mbar to approximately 0.25 at 900 mbar. From a combination of ab initio quantum chemical calculations on the molecular properties of the acetyl (CH3CO) radical and its unimolecular fragmentation as well as the application of statistical (RRKM) and dynamical calculations we show that CO production results from prompt secondary fragmentation (via(2a)) of the internally excited primary CH3CO* photolysis product with an excess energy of approximately 62.8 kJ mol(-1). Hence, our findings are consistent with a consecutive photochemically induced decomposition model, viz. step (1): CH3COCH3+hv--> CH3CO*+ CH3, step (2a): CH3CO*--> CH3+ CO or step (2b) CH3CO*-(+M)--> CH3CO. Formation of CO via a direct and/or concerted channel CH3COCH3+hv--> 2CH(3)+ CO (1') is considered to be unimportant.