Substitution Reactions in the Pyrolysis of Acetone Revealed through a Modeling, Experiment, Theory Paradigm.

The development of high-fidelity mechanisms for chemically reactive systems is a challenging process that requires the compilation of rate descriptions for a large and somewhat ill-defined set of reactions. The present unified combination of modeling, experiment, and theory provides a paradigm for improving such mechanism development efforts. Here we combine broadband rotational spectroscopy with detailed chemical modeling based on rate constants obtained from automated ab initio transition state theory-based master equation calculations and high-level thermochemical parametrizations. Broadband rotational spectroscopy offers quantitative and isomer-specific detection by which branching ratios of polar reaction products may be obtained. Using this technique, we observe and characterize products arising from H atom substitution reactions in the flash pyrolysis of acetone (CH3C(O)CH3) at a nominal temperature of 1800 K. The major product observed is ketene (CH2CO). Minor products identified include acetaldehyde (CH3CHO), propyne (CH3CCH), propene (CH2CHCH3), and water (HDO). Literature mechanisms for the pyrolysis of acetone do not adequately describe the minor products. The inclusion of a variety of substitution reactions, with rate constants and thermochemistry obtained from automated ab initio kinetics predictions and Active Thermochemical Tables analyses, demonstrates an important role for such processes. The pathway to acetaldehyde is shown to be a direct result of substitution of acetone's methyl group by a free H atom, while propene formation arises from OH substitution in the enol form of acetone by a free H atom.

Reaction Profiles and Kinetics for Radical-Radical Hydrogen Abstraction via Multireference Coupled Cluster Theory.

Radical-radical abstractions in hydrocarbon oxidation chemistry are disproportionation reactions that are generally exothermic with little or no barrier yet are underappreciated and poorly studied. Such challenging multireference electronic structure problems are tackled here using the recently developed state-specific multireference coupled cluster methods Mk-MRCCSD and Mk-MRCCSD(T), as well as the companion perturbation theory Mk-MRPT2 and the established MRCISD, MRCISD+Q, and CASPT2 approaches. Reaction paths are investigated for five prototypes involving radical-radical hydrogen abstraction: H + BeH → H2+ Be, H + NH2 → H2 + NH, CH3 + C2H5 → CH4 + C2H4, H + C2H5 → H2 + C2H4, and H + HCO → H2 + CO. Full configuration interaction (FCI) benchmark computations for the H + BeH, H + NH2, and H + HCO reactions prove that Mk-MRCCSD(T) provides superior accuracy for the interaction energies in the entrance channel, with mean absolute errors less than 0.3 kcal mol-1 and percentage deviations less than 10% over the fragment separations of relevance to kinetics. To facilitate combustion studies, energetics for the CH3 + C2H5, H + C2H5, and H + HCO reactions were computed at each level of theory with correlation-consistent basis sets (cc-pVXZ, X = T, Q, 5) and extrapolated to the complete basis set (CBS) limit. These CBS energies were coupled with CASPT2 projected vibrational frequencies along a minimum energy path to obtain rate constants for these three reactions. The rigorous Mk-MRCCSD(T)/CBS results demonstrate unequivocally that these three reactions proceed with no barrier in the entrance channel, contrary to some earlier predictions. Mk-MRCCSD(T) also reveals that the economical CASPT2 method performs well for large interfragment separations but may deteriorate substantially at shorter distances.

Anharmonic Rovibrational Partition Functions at High Temperatures: Tests of Reduced-Dimensional Models for Systems with up to Three Fluxional Modes.

Curvilinear coordinate Monte Carlo phase space integration and a series of full-dimensional fitted potential energy surfaces are used to study the effectiveness of reduced-dimensional models for predicting rovibrational anharmonicity at high temperatures. Fully coupled and fully anharmonic, but classical, rovibrational partition functions Q are computed for 14 species with two or three fluxional modes (inversions or torsions) and as many as 30 degrees of freedom. These results are converted to semiclassical anharmonicity correction factors f and are analyzed alongside results obtained previously for 22 species with up to two fluxional modes. As expected, fluxional species show considerable variation in f at high temperatures; f is as small as 0.2 for acetone and is as large as 9 for methylene glycol at 2500 K. This set of full-dimensional results is used to test the accuracy of reduced-dimensional models where fluxional modes are treated as coupled to one another but as separable from the remaining nonfluxional modes. For most systems, we find that such an approximation is accurate at high temperatures, with average errors in Q of just ∼25%. For some systems, however, larger errors are found, and these are attributed to strong coupling of the fluxional modes to one or more nonfluxional modes. In particular, we identify strong coupling to low-frequency bends for some systems, and we show that by comparing curvilinear and rectilinear harmonic frequencies for the fluxional modes, we can estimate the effect of this coupling on rovibrational anharmonicity. We also quantify the accuracy of the more severe but common assumption of treating fluxional modes as separable from one another, that is, as sets of uncoupled one-dimensional inversions and torsions. This approach can work well for methyl and alkyl rotors, but it is shown to have errors as large as a factor of 7 at high temperatures for more complex systems. Finally, we note that while the present analysis focuses on the treatment of fluxional modes, the collective anharmonicity correction associated with the more numerous nonfluxional modes, although simpler to describe, comprises a significant fraction of the overall anharmonicity.

Active Thermochemical Tables: The Partition Function of Hydroxymethyl (CH2OH) Revisited.

The best currently available set of temperature-dependent nonrigid rotor anharmonic oscillator (NRRAO) thermochemical and thermophysical properties of hydroxymethyl radical is presented. The underlying partition function relies on a critically evaluated complement of accurate experimental and theoretical data and is constructed using a two-pronged strategy that combines contributions from large amplitude motions obtained from direct counts, with contributions from the other internal modes of motion obtained from analytic NRRAO expressions. The contributions from the two strongly coupled large-amplitude motions of CH2OH, OH torsion and CH2 wag, are based on energy levels obtained by solving the appropriate two-dimensional projection of a fully dimensional potential energy surface that was recently obtained at the CCSD(T)/cc-pVTZ level of theory. The contributions of the remaining seven, more rigid, vibrational modes and of the external rotations are captured by NRRAO corrections to the standard rigid rotor harmonic oscillator (RRHO) treatment, which include corrections for vibrational anharmonicities, rotation-vibration interaction, Coriolis effects, and low temperature. The basic spectroscopic constants needed for the construction of the initial RRHO partition function rely on experimental ground-state rotational constants and the best available experimental fundamentals, additionally complemented by fundamentals obtained from the variational solution of the full-dimensional potential energy surface using a recently developed two-component multilayer Lanczos algorithm. The higher-order spectroscopic constants necessary for the NRRAO corrections are extracted from a second-order variational perturbation treatment (VPT2) of the same potential energy surface. The Lanczos solutions of the fully dimensional surface are validated against available experimental data, and the VPT2 results and the solutions of the reduced dimensionality surface are validated both against the Lanczos solutions and available experiments. The NRRAO thermophysical and thermochemical properties, given both in tabular form and as seven- and nine-coefficient NASA polynomials, are compared to previous results. In addition, the latest ATcT values for the enthalpy of formation of CH2OH at 298.15 K (0 K), -16.75 ± 0.27 kJ/mol (-10.45 ± 0.27 kJ/mol), and of other related CH nO m species ( n = 0-4, m = 0,1) are reported, together with a plethora of related bond dissociation enthalpies (BDEs), such as the C-H, O-H, and C-O bond dissociation enthalpies of methanol, 402.16 ± 0.26 kJ/mol (395.61 ± 0.26 kJ/mol), 440.34 ± 0.26 kJ/mol (434.86 ± 0.26 kJ/mol), and 384.85 ± 0.15 kJ/mol (377.14 ± 0.15 kJ/mol), respectively, and analogous BDEs for hydroxymethyl, 343.67 ± 0.37 kJ/mol (339.16 ± 0.37 kJ/mol), 125.54 ± 0.28 kJ/mol (121.11 ± 0.28 kJ/mol), and 445.86 ± 0.29 kJ/mol (438.76 ± 0.29 kJ/mol), respectively. The reasons governing the alternation between strong and weak sequential H atom BDEs of methanol are also discussed.

Nascent energy distribution of the Criegee intermediate CH2OO from direct dynamics calculations of primary ozonide dissociation.

Ozonolysis produces chemically activated carbonyl oxides (Criegee intermediates, CIs) that are either stabilized or decompose directly. This branching has an important impact on atmospheric chemistry. Prior theoretical studies have employed statistical models for energy partitioning to the CI arising from dissociation of the initially formed primary ozonide (POZ). Here, we used direct dynamics simulations to explore this partitioning for decomposition of c-C2H4O3, the POZ in ethylene ozonolysis. A priori estimates for the overall stabilization probability were then obtained by coupling the direct dynamics results with master equation simulations. Trajectories were initiated at the concerted cycloreversion transition state, as well as the second transition state of a stepwise dissociation pathway, both leading to a CI (H2COO) and formaldehyde (H2CO). The resulting CI energy distributions were incorporated in master equation simulations of CI decomposition to obtain channel-specific stabilized CI (sCI) yields. Master equation simulations of POZ formation and decomposition, based on new high-level electronic structure calculations, were used to predict yields for the different POZ decomposition channels. A non-negligible contribution of stepwise POZ dissociation was found, and new mechanistic aspects of this pathway were elucidated. By combining the trajectory-based channel-specific sCI yields with the channel branching fractions, an overall sCI yield of (48 ± 5)% was obtained. Non-statistical energy release was shown to measurably affect sCI formation, with statistical models predicting significantly lower overall sCI yields (∼30%). Within the range of experimental literature values (35%-54%), our trajectory-based calculations favor those clustered at the upper end of the spectrum.

Anharmonic Rovibrational Partition Functions for Fluxional Species at High Temperatures via Monte Carlo Phase Space Integrals.

Monte Carlo phase space integration (MCPSI) is used to compute full dimensional and fully anharmonic, but classical, rovibrational partition functions for 22 small- and medium-sized molecules and radicals. Several of the species considered here feature multiple minima and low-frequency nonlocal motions, and efficiently sampling these systems is facilitated using curvilinear (stretch, bend, and torsion) coordinates. The curvilinear coordinate MCPSI method is demonstrated to be applicable to the treatment of fluxional species with complex rovibrational structures and as many as 21 fully coupled rovibrational degrees of freedom. Trends in the computed anharmonicity corrections are discussed. For many systems, rovibrational anharmonicities at elevated temperatures are shown to vary consistently with the number of degrees of freedom and with temperature once rovibrational coupling and torsional anharmonicity are accounted for. Larger corrections are found for systems with complex vibrational structures, such as systems with multiple large-amplitude modes and/or multiple minima.

Time-Resolved Kinetic Chirped-Pulse Rotational Spectroscopy in a Room-Temperature Flow Reactor.

Chirped-pulse Fourier transform millimeter-wave spectroscopy is a potentially powerful tool for studying chemical reaction dynamics and kinetics. Branching ratios of multiple reaction products and intermediates can be measured with unprecedented chemical specificity; molecular isomers, conformers, and vibrational states have distinct rotational spectra. Here we demonstrate chirped-pulse spectroscopy of vinyl cyanide photoproducts in a flow tube reactor at ambient temperature of 295 K and pressures of 1-10 µbar. This in situ and time-resolved experiment illustrates the utility of this novel approach to investigating chemical reaction dynamics and kinetics. Following 193 nm photodissociation of CH2CHCN, we observe rotational relaxation of energized HCN, HNC, and HCCCN photoproducts with 10 µs time resolution and sample the vibrational population distribution of HCCCN. The experimental branching ratio HCN/HCCCN is compared with a model based on RRKM theory using high-level ab initio calculations, which were in turn validated by comparisons to Active Thermochemical Tables enthalpies.

Ab Initio Computations and Active Thermochemical Tables Hand in Hand: Heats of Formation of Core Combustion Species.

The fidelity of combustion simulations is strongly dependent on the accuracy of the underlying thermochemical properties for the core combustion species that arise as intermediates and products in the chemical conversion of most fuels. High level theoretical evaluations are coupled with a wide-ranging implementation of the Active Thermochemical Tables (ATcT) approach to obtain well-validated high fidelity predictions for the 0 K heat of formation for a large set of core combustion species. In particular, high level ab initio electronic structure based predictions are obtained for a set of 348 C, N, O, and H containing species, which corresponds to essentially all core combustion species with 34 or fewer electrons. The theoretical analyses incorporate various high level corrections to base CCSD(T)/cc-pVnZ analyses (n = T or Q) using H2, CH4, H2O, and NH3 as references. Corrections for the complete-basis-set limit, higher-order excitations, anharmonic zero-point energy, core-valence, relativistic, and diagonal Born-Oppenheimer effects are ordered in decreasing importance. Independent ATcT values are presented for a subset of 150 species. The accuracy of the theoretical predictions is explored through (i) examination of the magnitude of the various corrections, (ii) comparisons with other high level calculations, and (iii) through comparison with the ATcT values. The estimated 2σ uncertainties of the three methods devised here, ANL0, ANL0-F12, and ANL1, are in the range of ±1.0-1.5 kJ/mol for single-reference and moderately multireference species, for which the calculated higher order excitations are 5 kJ/mol or less. In addition to providing valuable references for combustion simulations, the subsequent inclusion of the current theoretical results into the ATcT thermochemical network is expected to significantly improve the thermochemical knowledge base for less-well studied species.

Accurate Anharmonic Zero-Point Energies for Some Combustion-Related Species from Diffusion Monte Carlo.

Full-dimensional analytic potential energy surfaces based on CCSD(T)/cc-pVTZ calculations have been determined for 48 small combustion-related molecules. The analytic surfaces have been used in Diffusion Monte Carlo calculations of the anharmonic zero-point energies. The resulting anharmonicity corrections are compared to vibrational perturbation theory results based both on the same level of electronic structure theory and on lower-level electronic structure methods (B3LYP and MP2).

Comment on "A novel and facile decay path of Criegee intermediates by intramolecular insertion reactions via roaming transition states" [J. Chem. Phys. 142, 124312 (2015)].

We discuss the recent report of a roaming type mechanism for the decomposition of the Criegee intermediate. We show that the predicted barrier height for this new pathway is too low by ∼30 kcal/mol owing to an inconsistent use of spin-restricted and spin-unrestricted calculations. As a result, this new pathway is not expected to compete significantly with the well-known dioxirane pathways for the decomposition of the Criegee intermediate.

Temperature and Pressure-Dependent Rate Coefficients for the Reaction of Vinyl Radical with Molecular Oxygen.

State-of-the-art calculations of the C2H3O2 potential energy surface are presented. A new method is described for computing the interaction potential for R + O2 reactions. The method, which combines accurate determination of the quartet potential along the doublet minimum energy path with multireference calculations of the doublet/quartet splitting, decreases the uncertainty in the doublet potential and thence the rate constants by more than a factor of 2. The temperature- and pressure-dependent rate coefficients are computed using variable reaction coordinate transition-state theory, variational transition-state theory, and conventional transition-state theory, as implemented in a new RRKM/ME code. The main bimolecular product channels are CH2O + HCO at lower temperatures and CH2CHO + O at higher temperatures. Above 10 atm, the collisional stabilization of CH2CHOO directly competes with these two product channels. CH2CHOO decomposes primarily to CH2O + HCO. The next two most significant bimolecular products are OCHCHO + H and (3)CHCHO + OH, and not C2H2 + HO2. C2H3 + O2 will be predominantly chain branching above 1700 K. Uncertainty analysis is presented for the two most important transition states. The uncertainties in these two barrier heights result in a significant uncertainty in the temperature at which CH2CHO + O overtakes all other product channels.

Thermal Dissociation and Roaming Isomerization of Nitromethane: Experiment and Theory.

The thermal decomposition of nitromethane provides a classic example of the competition between roaming mediated isomerization and simple bond fission. A recent theoretical analysis suggests that as the pressure is increased from 2 to 200 Torr the product distribution undergoes a sharp transition from roaming dominated to bond-fission dominated. Laser schlieren densitometry is used to explore the variation in the effect of roaming on the density gradients for CH3NO2 decomposition in a shock tube for pressures of 30, 60, and 120 Torr at temperatures ranging from 1200 to 1860 K. A complementary theoretical analysis provides a novel exploration of the effects of roaming on the thermal decomposition kinetics. The analysis focuses on the roaming dynamics in a reduced dimensional space consisting of the rigid-body motions of the CH3 and NO2 radicals. A high-level reduced-dimensionality potential energy surface is developed from fits to large-scale multireference ab initio calculations. Rigid body trajectory simulations coupled with master equation kinetics calculations provide high-level a priori predictions for the thermal branching between roaming and dissociation. A statistical model provides a qualitative/semiquantitative interpretation of the results. Modeling efforts explore the relation between the predicted roaming branching and the observed gradients. Overall, the experiments are found to be fairly consistent with the theoretically proposed branching ratio, but they are also consistent with a no-roaming scenario and the underlying reasons are discussed. The theoretical predictions are also compared with prior theoretical predictions, with a related statistical model, and with the extant experimental data for the decomposition of CH3NO2, and for the reaction of CH3 with NO2.

Resolving Some Paradoxes in the Thermal Decomposition Mechanism of Acetaldehyde.

The mechanism for the thermal decomposition of acetaldehyde has been revisited with an analysis of literature kinetics experiments using theoretical kinetics. The present modeling study was motivated by recent observations, with very sensitive diagnostics, of some unexpected products in high temperature microtubular reactor experiments on the thermal decomposition of CH3CHO and its deuterated analogs, CH3CDO, CD3CHO, and CD3CDO. The observations of these products prompted the authors of these studies to suggest that the enol tautomer, CH2CHOH (vinyl alcohol), is a primary intermediate in the thermal decomposition of acetaldehyde. The present modeling efforts on acetaldehyde decomposition incorporate a master equation reanalysis of the CH3CHO potential energy surface (PES). The lowest-energy process on this PES is an isomerization of CH3CHO to CH2CHOH. However, the subsequent product channels for CH2CHOH are substantially higher in energy, and the only unimolecular process that can be thermally accessed is a reisomerization to CH3CHO. The incorporation of these new theoretical kinetics predictions into models for selected literature experiments on CH3CHO thermal decomposition confirms our earlier experiment and theory-based conclusions that the dominant decomposition process in CH3CHO at high temperatures is C-C bond fission with a minor contribution (∼10-20%) from the roaming mechanism to form CH4 and CO. The present modeling efforts also incorporate a master-equation analysis of the H + CH2CHOH potential energy surface. This bimolecular reaction is the primary mechanism for removal of CH2CHOH, which can accumulate to minor amounts at high temperatures, T > 1000 K, in most lab-scale experiments that use large initial concentrations of CH3CHO. Our modeling efforts indicate that the observation of ketene, water, and acetylene in the recent microtubular experiments are primarily due to bimolecular reactions of CH3CHO and CH2CHOH with H-atoms and have no bearing on the unimolecular decomposition mechanism of CH3CHO. The present simulations also indicate that experiments using these microtubular reactors when interpreted with the aid of high-level theoretical calculations and kinetics modeling can offer insights into the chemistry of elusive intermediates in the high-temperature pyrolysis of organic molecules.

Predictive a priori pressure-dependent kinetics.

The ability to predict the pressure dependence of chemical reaction rates would be a great boon to kinetic modeling of processes such as combustion and atmospheric chemistry. This pressure dependence is intimately related to the rate of collision-induced transitions in energy E and angular momentum J. We present a scheme for predicting this pressure dependence based on coupling trajectory-based determinations of moments of the E,J-resolved collisional transfer rates with the two-dimensional master equation. This completely a priori procedure provides a means for proceeding beyond the empiricism of prior work. The requisite microcanonical dissociation rates are obtained from ab initio transition state theory. Predictions for the CH4 = CH3 + H and C2H3 = C2H2 + H reaction systems are in excellent agreement with experiment.

Electronic states of the quasilinear molecule propargylene (HCCCH) from negative ion photoelectron spectroscopy.

We use gas-phase negative ion photoelectron spectroscopy to study the quasilinear carbene propargylene, HCCCH, and its isotopologue DCCCD. Photodetachment from HCCCH­ affords the X̃(3B) ground state of HCCCH and its ã(1A), b̃ (1B), d̃(1A2), and B̃(3A2) excited states. Extended, negatively anharmonic vibrational progressions in the X̃(3B) ground state and the open-shell singlet b̃ (1B) state arise from the change in geometry between the anion and the neutral states and complicate the assignment of the origin peak. The geometry change arising from electron photodetachment results in excitation of the ν4 symmetric CCH bending mode, with a measured fundamental frequency of 363 ± 57 cm(­1) in the X̃(3B) state. Our calculated harmonic frequency for this mode is 359 cm(­1). The Franck­Condon envelope of this progression cannot be reproduced within the harmonic approximation. The spectra of the ã(1A), d̃(1A2), and B̃(3A2) states are each characterized by a short vibrational progression and a prominent origin peak, establishing that the geometries of the anion and these neutral states are similar. Through comparison of the HCCCH­ and DCCCD­ photoelectron spectra, we measure the electron affinity of HCCCH to be 1.156 ± (0.095)(0.010) eV, with a singlet­triplet splitting between the X̃(3B) and the ã(1A) states of ΔEST = 0.500 ± (0.01)(0.10) eV (11.5 ± (0.2)(2.3) kcal/mol). Experimental term energies of the higher excited states are T0 [b̃(1B)] = 0.94 ± (0.20)(0.22) eV, T0 [d̃(1A2)] = 3.30 ± (0.02)(0.10) eV, T0 [B̃(3A2)] = 3.58 ± (0.02)(0.10) eV. The photoelectron angular distributions show significant π character in all the frontier molecular orbitals, with additional σ character in orbitals that create the X̃(3B) and b̃(1B) states upon electron detachment. These results are consistent with a quasilinear, nonplanar, doubly allylic structure of X̃(3B) HCCCH with both diradical and carbene character.

Predictive theory for the addition and insertion kinetics of 1CH2 reacting with unsaturated hydrocarbons.

The reactions of singlet methylene, (1)CH2, with unsaturated hydrocarbons are of considerable significance to the formation and growth of polycyclic aromatic hydrocarbons (PAHs). In this work, we employ high level ab initio transition state theory to predict the high pressure rate coefficient for singlet methylene reacting with acetylene (C2H2), ethylene (C2H4), propyne (CH3CCH), propene (CH3CHCH2), allene (CH2CCH2), 1,3-butadiene (CH2CHCHCH2), 2-butyne (CH3CCCH3), and benzene (C6H6). Both addition and insertion channels are found to contribute significantly to the kinetics, with the insertion kinetics of increasing importance for larger hydrocarbons due to the increasing number of CH bonds and increasingly attractive interactions. We treat the addition kinetics with direct CASPT2 based variable-reaction-coordinate transition state theory. One-dimensional corrections to the CASPT2 interaction energies are obtained from geometry relaxation calculations and CCSD(T)/CBS evaluations. The insertion kinetics is treated with traditional variational TST methods employing CCSD(T)/CBS energies obtained along the CASPT2/cc-pVTZ distinguished coordinate reaction paths. The overall rate constant and branching fractions are obtained from a multiple transition state model that accounts for the physical distinction between tight inner and loose outer transition states. The predicted rate constants, which cover the range from 200 to 2000 K, are found to be in excellent agreement with the available experimental data, with a maximum observed discrepancy of about 40%.

Rate constant and branching fraction for the NH2 + NO2 reaction.

The NH2 + NO2 reaction has been studied experimentally and theoretically. On the basis of laser photolysis/LIF experiments, the total rate constant was determined over the temperature range 295-625 K as k1,exp(T) = 9.5 × 10(-7)(T/K)(-2.05) exp(-404 K/T) cm(3) molecule(-1) s(-1). This value is in the upper range of data reported for this temperature range. The reactions on the NH2 + NO2 potential energy surface were studied using high level ab initio transition state theory (TST) based master equation methods, yielding a rate constant of k1,theory(T) = 7.5 × 10(-12)(T/K)(-0.172) exp(687 K/T) cm(3) molecule(-1) s(-1), in good agreement with the experimental value in the overlapping temperature range. The two entrance channel adducts H2NNO2 and H2NONO lead to formation of N2O + H2O (R1a) and H2NO + NO (R1b), respectively. The pathways through H2NNO2 and H2NONO are essentially unconnected, even though roaming may facilitate a small flux between the adducts. High- and low-pressure limit rate coefficients for the various product channels of NH2 + NO2 are determined from the ab initio TST-based master equation calculations for the temperature range 300-2000 K. The theoretical predictions are in good agreement with the measured overall rate constant but tend to overestimate the branching ratio defined as ß = k1a/(k1a + k1b) at lower temperatures. Modest adjustments of the attractive potentials for the reaction yield values of k1a = 4.3 × 10(-6)(T/K)(-2.191) exp(-229 K/T) cm(3) molecule(-1) s(-1) and k1b = 1.5 × 10(-12)(T/K)(0.032) exp(761 K/T) cm(3) molecule(-1) s(-1), in good agreement with experiment, and we recommend these rate coefficients for use in modeling.

Unconventional Peroxy Chemistry in Alcohol Oxidation: The Water Elimination Pathway.

Predictive simulation for designing efficient engines requires detailed modeling of combustion chemistry, for which the possibility of unknown pathways is a continual concern. Here, we characterize a low-lying water elimination pathway from key hydroperoxyalkyl (QOOH) radicals derived from alcohols. The corresponding saddle-point structure involves the interaction of radical and zwitterionic electronic states. This interaction presents extreme difficulties for electronic structure characterizations, but we demonstrate that these properties of this saddle point can be well captured by M06-2X and CCSD(T) methods. Experimental evidence for the existence and relevance of this pathway is shown in recently reported data on the low-temperature oxidation of isopentanol and isobutanol. In these systems, water elimination is a major pathway, and is likely ubiquitous in low-temperature alcohol oxidation. These findings will substantially alter current alcohol oxidation mechanisms. Moreover, the methods described will be useful for the more general phenomenon of interacting radical and zwitterionic states.

Separability of tight and roaming pathways to molecular decomposition.

Recent studies have questioned the separability of the tight and roaming mechanisms to molecular decomposition. We explore this issue for a variety of reactions including MgH(2) → Mg + H(2), NCN → CNN, H(2)CO → H(2) + CO, CH(3)CHO → CH(4) + CO, and HNNOH → N(2) + H(2)O. Our analysis focuses on the role of second-order saddle points in defining global dividing surfaces that encompass both tight and roaming first-order saddle points. The second-order saddle points define an energetic criterion for separability of the two mechanisms. Furthermore, plots of the differential contribution to the reactive flux along paths connecting the first- and second-order saddle points provide a dynamic criterion for separability. The minimum in the differential reactive flux in the neighborhood of the second-order saddle point plays the role of a mechanism divider, with the presence of a strong minimum indicating that the roaming and tight mechanisms are dynamically distinct. We show that the mechanism divider is often, but not always, associated with a second-order saddle point. For the formaldehyde and acetaldehyde reactions, we find that the minimum energy geometry on a conical intersection is associated with the mechanism divider for the tight and roaming processes. For HNNOH, we again find that the roaming and tight processes are dynamically separable but we find no intrinsic feature of the potential energy surface associated with the mechanism divider. Overall, our calculations suggest that roaming and tight mechanisms are generally separable over broad ranges of energy covering most kinetically relevant regimes.