Search

A transient carbon-centered hydroperoxyalkyl intermediate (•QOOH) in the oxidation of cyclopentane is identified by IR action spectroscopy with time-resolved unimolecular decay to hydroxyl (OH) radical products that are detected by UV laser-induced fluorescence. Two nearly degenerate •QOOH isomers, β- and γ-QOOH, are generated by H atom abstraction of the cyclopentyl hydroperoxide precursor. Fundamental and first overtone OH stretch transitions and combination bands of •QOOH are observed and compared with anharmonic frequencies computed by second-order vibrational perturbation theory. An OH stretch transition is also observed for a conformer arising from torsion about a low-energy CCOO barrier. Definitive identification of the β-QOOH isomer relies on its significantly lower transition state (TS) barrier to OH products, which results in rapid unimolecular decay and near unity branching to OH products. A benchmarking approach is utilized to compute high-accuracy stationary point energies, most importantly TS barriers, for cyclopentane oxidation (C 5 H 9 O 2 ), building on higher level reference calculations for ethane oxidation (C 2 H 5 O 2 ). The experimental OH product appearance rates are compared with computed statistical microcanonical rates using RRKM theory, including heavy-atom tunneling, thereby validating the computed TS barrier. The results are extended to thermal unimolecular decay rate constants at temperatures and pressures relevant to cyclopentane combustion via master-equation modeling. The various torsional and ring puckering states of the wells and transition states are explicitly considered in these calculations.

Infrared (IR) action spectroscopy is utilized to characterize carbon-centered hydroperoxy-cyclohexyl radicals (·QOOH) transiently formed in cyclohexane oxidation. The oxidation pathway leads to three nearly degenerate ·QOOH isomers, β-, γ-, and δ-QOOH, which are generated in the laboratory by H-atom abstraction from the corresponding ring sites of the cyclohexyl hydroperoxide (CHHP) precursor. The IR spectral features of jet-cooled and stabilized ·QOOH radicals are observed from 3590 to 7010 cm-1 (∼10-20 kcal mol-1) at energies in the vicinity of the transition state (TS) barrier leading to OH radicals that are detected by ultraviolet laser-induced fluorescence. The experimental approach affords selective detection of β-QOOH, arising from its significantly lower TS barrier to OH products compared to γ and δ isomers, which results in rapid unimolecular decay and near unity branching to OH products. The observed IR spectrum of β-QOOH includes fundamental and overtone OH stretch transitions, overtone CH stretch transitions, and combination bands involving OH or CH stretch with lower frequency modes. The assignment of β-QOOH spectral features is guided by anharmonic frequencies and intensities computed using second-order vibrational perturbation theory. The overtone OH stretch (2νOH) of β-QOOH is shifted only a few wavenumbers from that observed for the CHHP precursor, yet they are readily distinguished by their prompt vs slow dissociation rates to OH products., (© 2024 Author(s). Published under an exclusive license by AIP Publishing.)

Recent work in combustion and atmospheric chemistry has revealed cases in which diastereomers must be distinguished to accurately model a reacting flow. This paper presents an open-source framework for introducing such stereoisomer resolution into a kinetic mechanism. We detail our definitions and algorithms for labeling and enumerating the stereoisomers of a molecule and then generalize our system to describe the transition state (TS) of a reaction. This allows for the stereospecific enumeration of reactants and products while accounting for "fleeting" stereochemistry that is unique to the TS. We also present the AutoMech Chemical Identifier (AMChI), an InChI-like string identifier that accounts for stereocenters omitted by InChI. This identifier is extended to describe the TSs of reactions, providing a universal lookup key for specific reaction channels. The final piece of our methodology is an analytic formula to remove redundancy from a stereoresolved mechanism when its enantiomers exist as a racemic mixture, making it as compact as possible while fully accounting for the differences between diastereomers. In applying our methodology to two subsets of the NUIGMech1.1 mechanism, we find that our approach reduces the extra species added for large-fuel oxidation from 2231 (133%, full expansion) to 694 (41%, nonredundant expansion). We also find that for pyrolysis more than a quarter of the species in the expanded mechanism cannot be properly described by an InChI string, requiring an AMChI string to communicate their identity. Finally, we find that roughly one-quarter of the large-fuel oxidation reactions and one-third of the pyrolysis reactions include fleeting TS stereochemistry, which may have relevant effects on their kinetics.

Machine learning (ML) provides a great opportunity for the construction of models with improved accuracy in classical molecular dynamics (MD). However, the accuracy of a ML trained model is limited by the quality and quantity of the training data. Generating large sets of accurate ab initio training data can require significant computational resources. Furthermore, inconsistent or incompatible data with different accuracies obtained using different methods may lead to biased or unreliable ML models that do not accurately represent the underlying physics. Recently, transfer learning showed its potential for avoiding these problems as well as for improving the accuracy, efficiency, and generalization of ML models using multifidelity data. In this work, ab initio trained ML-based MD (aML-MD) models are developed through transfer learning using DFT and multireference data from multiple sources with varying accuracy within the Deep Potential MD framework. The accuracy of the force field is demonstrated by calculating rate constants for the H + HO 2 → H 2 + 3 O 2 reaction using quasi-classical trajectories. We show that the aML-MD model with transfer learning can accurately predict the rate constants while reducing the computational cost by more than five times compared to the use of more expensive quantum chemistry training data sets. Hence, the aML-MD model with transfer learning shows great potential in using multifidelity data to reduce the computational cost involved in generating the training set for these potentials.

The oxidation of cycloalkanes is important in the combustion of transportation fuels and in atmospheric secondary organic aerosol formation. A transient carbon-centered radical intermediate (•QOOH) in the oxidation of cyclohexane is identified through its infrared fingerprint and time- and energy-resolved unimolecular dissociation dynamics to hydroxyl (OH) radical and bicyclic ether products. Although the cyclohexyl ring structure leads to three nearly degenerate •QOOH isomers (β-, γ-, and δ-QOOH), their transition state (TS) barriers to OH products are predicted to differ considerably. Selective characterization of the β-QOOH isomer is achieved at excitation energies associated with the lowest TS barrier, resulting in rapid unimolecular decay to OH products that are detected. A benchmarking approach is employed for the calculation of high-accuracy stationary point energies, in particular TS barriers, for cyclohexane oxidation (C 6 H 11 O 2 ), building on higher-level reference calculations for the smaller ethane oxidation (C 2 H 5 O 2 ) system. The isomer-specific characterization of β-QOOH is validated by comparison of experimental OH product appearance rates with computed statistical microcanonical rates, including significant heavy-atom tunneling, at energies in the vicinity of the TS barrier. Master-equation modeling is utilized to extend the results to thermal unimolecular decay rate constants at temperatures and pressures relevant to cyclohexane combustion., Competing Interests: Competing interests statement:D.E.H. participated in a Faraday Discussion on Unimolecular Reactions and published general discussions with several co-authors within the last 4 y. The authors declare no other potential conflicts of interest.

The formation of low-volatility organic compounds in the ozonolysis of α-pinene, the dominant atmospheric monoterpene, provides an important route to aerosol formation. In this work, we consider a previously unexplored set of pathways for the formation of highly oxygenated molecules in α-pinene ozonolysis. Pioneering, direct experimental observations of Lester and co-workers have demonstrated a significant production of hydroxycarbonyl products in the dissociation of Criegee intermediates. Theoretical analyses indicate that this production arises from OH roaming-induced pathways during the OO fission of the vinylhydroperoxides (VHPs), which in turn come from internal H transfers in the Criegee intermediates. Ab initio kinetics computations are used here to explore the OH roaming-induced channels that arise from the ozonolysis of α-pinene. For computational reasons, the calculations consider a surrogate for α-pinene, where two spectator methyl groups are replaced with H atoms. Multireference electronic structure calculations are used to illustrate a variety of energetically accessible OH roaming pathways for the four VHPs arising from the ozonolysis of this α-pinene surrogate. Ab initio transition-state theory-based master equation calculations indicate that for the dissociation of stabilized VHPs, these OH roaming pathways are kinetically significant with a branching that generally increases from ∼20% at room temperature up to ∼70% at lower temperatures representative of the troposphere. For one of the VHPs, this branching already exceeds 60% at room temperature. For the overall ozonolysis process, these branching ratios would be greatly reduced by a limited branching to the stabilized VHP, although there would also be some modest roaming fraction for the nonthermal VHP dissociation process. The strong exothermicities of the roaming-induced isomerizations/additions and abstractions suggest new routes to fission of the cyclobutane rings. Such ring fissions would facilitate further autoxidation reactions, thereby providing a new route for producing highly oxygenated nonvolatile precursors to aerosol formation.

Methyl formate (MF; CH 3 OCHO) is the smallest representative of esters, which are common components of biodiesel. The present study characterizes the thermal dissociation kinetics of the radicals formed by H atom abstraction from MF─CH 3 OCO and CH 2 OCHO─through a combination of modeling, experiment, and theory. For the experimental effort, excimer laser photolysis of Cl 2 was used as a source of Cl atoms to initiate reactions with MF in the gas phase. Time-resolved species profiles of MF, Cl 2 , HCl, CO 2 , CH 3 , CH 3 Cl, CH 2 O, and CH 2 ClOCHO were measured and quantified using photoionization mass spectrometry at temperatures of 400-750 K and 10 Torr. The experimental data were simulated using a kinetic model, which was informed by ab initio-based theoretical kinetics calculations and included chlorine chemistry and secondary reactions of radical decomposition products. We calculated the rate coefficients for the H-abstraction reactions Cl + MF → HCl + CH 3 OCO (R1a) and Cl + MF → HCl + CH 2 OCHO (R1b): k 1a,theory = 6.71 × 10 -15 · T 1.14 ·exp(-606/ T ) cm 3 /molecule·s; k 1b,theory = 4.67 × 10 -18 · T 2.21 ·exp(-245/ T ) cm 3 /molecule·s over T = 200-2000 K. Electronic structure calculations indicate that the barriers to CH 3 OCO and CH 2 OCHO dissociation are 13.7 and 31.6 kcal/mol and lead to CH 3 + CO 2 (R3) and CH 2 O + HCO (R5), respectively. The master equation-based theoretical rate coefficients are k 3,theory ( P = ∞) = 2.94 × 10 9 · T 1.21 ·exp(-6209/ T ) s -1 and k 5,theory ( P = ∞) = 8.45 × 10 8 · T 1.39 ·exp(-15132/ T ) s -1 over T = 300-1500 K. The calculated branching fractions into R1a and R1b and the rate coefficient for R5 were validated by modeling of the experimental species time profiles and found to be in excellent agreement with theory. Additionally, we found that the bimolecular reactions CH 2 OCHO + Cl, CH 2 OCHO + Cl 2 , and CH 3 + Cl 2 were critical to accurately model the experimental data and constrain the kinetics of MF-radicals. Inclusion of the kinetic parameters determined in this study showed a significant impact on combustion simulations of larger methyl esters, which are considered as biodiesel surrogates.

Quantitative descriptions of non-adiabatic transition rates at intermediate temperatures are challenging due to the simultaneous importance of quantum and anharmonic effects. In this paper, the interplay between quantum effects-for motion across or along the seam of crossing-and anharmonicity in the seam potential is considered within the weak coupling limit. The well-known expression for quantized 1-D motion across the seam (i.e., tunneling) in the linear terms approximation is derived in the thermal domain using the Lagrangian formalism, which is then applied to the case when tunneling is distributed along the seam of crossing (treating motion along the seam classically). For high-frequency quantum modes, a vibrationally adiabatic (VA) approach is developed that introduces to the non-adiabatic rate constant a factor associated with high-frequency wavefunction overlap; this approach treats the high-frequency motion along the seam quantum mechanically. To test these methodologies, the reaction N2O ↔ N2 + O(3P) was chosen. CCSD(T)-F12b/cc-pVTZ-F12 explorations of the 3A'-1A' seam of N2O revealed that seam anharmonicity has a strong effect on the rate constant (a factor of ∼20 at 2000 K). Several quantum effects were found to be significant at intermediate/lower temperatures, including the quantum N-N vibration that was coupled with seam anharmonicity using the VA approach. Finally, a 1-D approximation to non-adiabatic instanton theory is presented to estimate the validity limit of the linear terms model at low temperatures (∼250 K for N2O). We recommend that the assumptions built into many statistical theories for non-adiabatic reactions-harmonic behavior, classical motion, linear terms, and weak coupling-should be verified on a case-by-case basis., (© 2023 Author(s). Published under an exclusive license by AIP Publishing.)

Methyl-ethyl-substituted Criegee intermediate (MECI) is a four-carbon carbonyl oxide that is formed in the ozonolysis of some asymmetric alkenes. MECI is structurally similar to the isoprene-derived methyl vinyl ketone oxide (MVK-oxide) but lacks resonance stabilization, making it a promising candidate to help us unravel the effects of size, structure, and resonance stabilization that influence the reactivity of atmospherically important, highly functionalized Criegee intermediates. We present experimental and theoretical results from the first bimolecular study of MECI in its reaction with SO 2 , a reaction that shows significant sensitivity to the Criegee intermediate structure. Using multiplexed photoionization mass spectrometry, we obtain a rate coefficient of (1.3 ± 0.3) × 10 -10 cm 3 s -1 (95% confidence limits, 298 K, 10 Torr) and demonstrate the formation of SO 3 under our experimental conditions. Through high-level theory, we explore the effect of Criegee intermediate structure on the minimum energy pathways for their reactions with SO 2 and obtain modified Arrhenius fits to our predictions for the reaction of both syn and anti conformers of MECI with SO 2 ( k syn = 4.42 × 10 11 T -7.80 exp(-1401/T) cm 3 s -1 and k anti = 1.26 × 10 11 T -7.55 exp(-1397/T) cm 3 s -1 ). Our experimental and theoretical rate coefficients (which are in reasonable agreement at 298 K) show that the reaction of MECI with SO 2 is significantly faster than MVK-oxide + SO 2 , demonstrating the substantial effect of resonance stabilization on Criegee intermediate reactivity.

Recent modeling studies of ammonia oxidation, which are motivated by the prospective role of ammonia as a zero-carbon fuel, have indicated significant discrepancies among the existing literature mechanisms. In this study, high-level theoretical kinetics predictions have been obtained for reactions on the NH 2 O potential energy surface, including the NH 2 + O, HNO + H, and NH + OH reactions. These reactions have previously been highlighted as important reactions in NH 3 oxidation with high sensitivity and high uncertainty. The potential energy surface is explored with coupled cluster calculations, including large basis sets and high-level corrections to yield high-accuracy (∼0.2 kcal/mol 2σ uncertainty) estimates of the stationary point energies. Variational transition state theory is used to predict the microcanonical rate constants, which are then incorporated in master equation treatments of the temperature- and pressure-dependent kinetics. For radical-radical channels, the microcanonical rates are obtained from variable reaction coordinate transition state theory implementing directly evaluated multireference electronic energies. The analysis yields predictions for the total rate constants as well as the branching ratios. We find that the NO + H 2 channel contributes 10% of the total NH 2 + O flux at combustion temperatures, although this channel is not included in modern NH 3 oxidation mechanisms. Modeling is used to illustrate the ramifications of these rate predictions on the kinetics of NH 3 oxidation and NO x formation. The present results for NH 2 + O are important for predicting the chain branching and formation of NO in the oxidation of NH 3 and thermal DeNO x .

Alkene ozonolysis generates short-lived Criegee intermediates that are a significant source of hydroxyl (OH) radicals. This study demonstrates that roaming of the separating OH radicals can yield alternate hydroxycarbonyl products, thereby reducing the OH yield. Specifically, hydroxybutanone has been detected as a stable product arising from roaming in the unimolecular decay of the methyl-ethyl-substituted Criegee intermediate (MECI) under thermal flow cell conditions. The dynamical features of this novel multistage dissociation plus a roaming unimolecular decay process have also been examined with ab initio kinetics calculations. Experimentally, hydroxybutanone isomers are distinguished from the isomeric MECI by their higher ionization threshold and distinctive photoionization spectra. Moreover, the exponential rise of the hydroxybutanone kinetic time profile matches that for the unimolecular decay of MECI. A weaker methyl vinyl ketone (MVK) photoionization signal is also attributed to OH roaming. Complementary multireference electronic structure calculations have been utilized to map the unimolecular decay pathways for MECI, starting with 1,4 H atom transfer from a methyl or methylene group to the terminal oxygen, followed by roaming of the separating OH and butanonyl radicals in the long-range region of the potential. Roaming via reorientation and the addition of OH to the vinyl group of butanonyl is shown to yield hydroxybutanone, and subsequent C-O elongation and H-transfer can lead to MVK. A comprehensive theoretical kinetic analysis has been conducted to evaluate rate constants and branching yields (ca. 10-11%) for thermal unimolecular decay of MECI to conventional and roaming products under laboratory and atmospheric conditions, consistent with the estimated experimental yield (ca. 7%).

The mechanism for hydrocarbon ring growth in sooting environments is still the subject of considerable debate. The reaction of phenyl radical (C 6 H 5 ) with propargyl radical (H 2 CCCH) provides an important prototype for radical-radical ring-growth pathways. We studied this reaction experimentally over the temperature range of 300-1000 K and pressure range of 4-10 Torr using time-resolved multiplexed photoionization mass spectrometry. We detect both the C 9 H 8 and C 9 H 7 + H product channels and report experimental isomer-resolved product branching fractions for the C 9 H 8 product. We compare these experiments to theoretical kinetics predictions from a recently published study augmented by new calculations. These ab initio transition state theory-based master equation calculations employ high-quality potential energy surfaces, conventional transition state theory for the tight transition states, and direct CASPT2-based variable reaction coordinate transition state theory (VRC-TST) for the barrierless channels. At 300 K only the direct adducts from radical-radical addition are observed, with good agreement between experimental and theoretical branching fractions, supporting the VRC-TST calculations of the barrierless entrance channel. As the temperature is increased to 1000 K we observe two additional isomers, including indene, a two-ring polycyclic aromatic hydrocarbon, and a small amount of bimolecular products C 9 H 7 + H. Our calculated branching fractions for the phenyl + propargyl reaction predict significantly less indene than observed experimentally. We present further calculations and experimental evidence that the most likely cause of this discrepancy is the contribution of H atom reactions, both H + indenyl (C 9 H 7 ) recombination to indene and H-assisted isomerization that converts less stable C 9 H 8 isomers into indene. Especially at low pressures typical of laboratory investigations, H-atom-assisted isomerization needs to be considered. Regardless, the experimental observation of indene demonstrates that the title reaction leads, either directly or indirectly, to the formation of the second ring in polycyclic aromatic hydrocarbons.

It is generally challenging to obtain high-accuracy predictions for the heat of formation for species with more than a handful of heavy atoms, such as those of importance in standard combustion mechanisms. To this end, we construct the CBH-ANL approach and illustrate that, for a set of 194 alkane oxidation species, it can be used to produce Δ H f (0 K) values with 2σ uncertainties of 0.2-0.5 kcal mol -1 . This set includes the alkanes, hydroperoxides, and alkyl, peroxy, and hydroperoxyalkyl radicals for 17 representative hydrocarbon fuels containing up to 10 heavy atoms with various degrees of branching in the alkane backbone. The CBH-ANL approach, automated in the QTC and AutoMech software suites, builds balanced chemical equations for the calculation of Δ H f (0 K), in which the reference species may be up to five heavy atoms. The high-level ANL0 and ANL1 reference Δ H f (0 K) values are further refined for even the largest of these reference species with a novel laddering approach. We perform a comprehensive quantification of the uncertainties for both the individual reference species (the largest of which is 0.15 kcal mol -1 ) and the propagation of those uncertainties when used in the calculation of Δ H f (0 K) for the 194 target species. We examine the sensitivity of the predicted Δ H f (0 K) values to (i) electronic energies from various methods, including ωB97X-D/cc-pVTZ, B2PLYP-D3/cc-pVTZ, CCSD(T)-F12b/cc-pVDZ-F12//B2PLYP-D3/cc-pVTZ, and CCSD(T)-F12b/cc-pVTZ-F12//B2PLYP-D3/cc-pVTZ; (ii) the zero-point vibrational energies (ZPVEs), where we consider harmonic ZPVEs as well as two scaling-based estimates of the anharmonic ZPVEs, all implemented for both ωB97X-D/cc-pVTZ and B2PLYP-D3/cc-pVTZ calculations; (iii) the particular CBH-ANL scheme employed; and (iv) the procedure for choosing the reference conformer for the analyses. The discussion concludes with a summary of the estimated overall uncertainty in the predictions and a validation of the predictions for the alkane subset.

The kinetics of peroxy radical (RO 2 ) reactions have been of long-standing interest in atmospheric and combustion chemistry. Nevertheless, the lack of kinetic studies at higher temperatures for their reactions with other radicals such as OH has precluded the inclusion of this class of reactions in detailed kinetics models developed for combustion applications. In this work, guided by the limited room-temperature experimental studies on selected alkyl-peroxy radicals and literature theoretical kinetics on the prototypical CH 3 O 2 + OH system, we have performed parametric studies on the effect of uncertainties in the rate coefficients and branching ratios to potential product channels for RO 2 + OH reactions at higher temperatures. Literature kinetics models were used to simulate autoignition delays, laminar flame speeds, and speciation profiles in flow and stirred reactors for a variety of common combustion-relevant fuels. Inclusion of RO 2 + OH reactions was found to retard autoignition in fuel-lean (φ = 0.5) mixtures of ethane and dimethyl ether in air. The observed effects were noticeably more pronounced in ozone-enriched combustion of ethane and dimethyl ether. The simulations also examined the influence of ozone doping levels, pressures, and equivalence ratios for both ethane and dimethyl ether oxidation. Sensitivity and flux analyses revealed that the RO 2 + OH reaction is a significant sink of RO 2 radicals at the early stage of autoignition, affecting fuel oxidation through RO 2 ↔ QOOH, RO 2 ↔ alkene + HO 2 , or RO 2 + HO 2 ↔ ROOH + O 2 . Additionally, the kinetic stability of the trioxide formed from RO 2 + OH reactions was investigated using master equation analyses. Last, we discuss other bimolecular reactions that are missing in literature kinetics models but are relevant to hydrocarbon oxidation initiated by external radical sources (plasma-enhanced, ozone-enriched combustion, etc.). The present simulations provide a strong motivation for better characterizing the bimolecular kinetics of peroxy radicals.

Hydroperoxyalkyl radicals (˙QOOH) are transient intermediates in the atmospheric oxidation of volatile organic compounds and combustion of hydrocarbon fuels in low temperature (<1000 K) environments. The carbon-centered ˙QOOH radicals are a critical juncture in the oxidation mechanism, but have generally eluded direct experimental observation of their structure, stability, and dissociation dynamics. Recently, this laboratory demonstrated that a prototypical ˙QOOH radical [˙CH 2 (CH 3 ) 2 COOH] can be synthesized by an alternative route, stabilized in a pulsed supersonic expansion, and characterized by its infrared (IR) spectroscopic signature and unimolecular dissociation rate to OH radical and cyclic ether products. The present study focuses on a partially deuterated ˙QOOD analog ˙CH 2 (CH 3 ) 2 COOD, generated in the laboratory by H-atom abstraction from partially deuterated tert -butyl hydroperoxide, (CH 3 ) 3 COOD. IR spectral features associated with jet-cooled and isolated ˙QOOD radicals are observed in the vicinity of the transition state (TS) barrier leading to OD radical and cyclic ether products. The overtone OD stretch (2 ν OD ) of ˙QOOD is identified by IR action spectroscopy with UV laser-induced fluorescence detection of OD products. Direct time-domain measurement of the unimolecular dissociation rate for ˙QOOD (2 ν OD ) extends prior rate measurements for ˙QOOH. Partial deuteration results in a small increase in the TS barrier predicted by high level electronic structure calculations due to changes in zero-point energies; the imaginary frequency is unchanged. Comparison of the unimolecular decay rates obtained experimentally with those predicted theoretically for both ˙QOOH and ˙QOOD confirm that unimolecular decay is enhanced by heavy-atom tunneling involving simultaneous O-O bond elongation and C-C-O angle contraction along the reaction pathway.

One hundred years ago, at an earlier Faraday Discussion meeting, Lindemann presented a mechanism that provides the foundation for contemplating the pressure dependence of unimolecular reactions. Since that time, our ability to model and predict the kinetics of such reactions has grown in leaps and bounds through the synergy of ever more sophisticated experimental studies with increasingly high level theoretical analyses. This review begins with a brief historical overview of the progress from the Lindemann mechanism to the master equation, which now provides the cornerstone for most theoretical analyses. The current status of ab initio transition state theory based implementations of the master equation is then reviewed, beginning with a discussion of the energy resolved chemical conversion rates, followed by a review of the pressure dependence of the thermal kinetics. The latter discussion focusses on the collisional energy and angular momentum transfer rates as well as the complexities of non-thermal effects and of multiple-well multiple-channel reactions. The synergy with recent state-of-the-art experiments is used to motivate these discussions. Attempts to automate the calculations are also briefly reviewed. Throughout the discussion, we provide our perspective on current understanding and continuing challenges and opportunities. One key challenge relates to the coupling of sequences of reactions, which frequently leads to deviations from the classic presumption of thermal reactants.

A novel allylic 1,6 hydrogen-atom-transfer mechanism is established through infrared activation of the 2-butenal oxide Criegee intermediate, resulting in very rapid unimolecular decay to hydroxyl (OH) radical products. A new precursor, Z / E -1,3-diiodobut-1-ene, is synthesized and photolyzed in the presence of oxygen to generate a new four-carbon Criegee intermediate with extended conjugation across the vinyl and carbonyl oxide groups that facilitates rapid allylic 1,6 H-atom transfer. A low-energy reaction pathway involving isomerization of 2-butenal oxide from a lower-energy ( tZZ ) conformer to a higher-energy ( cZZ ) conformer followed by 1,6 hydrogen transfer via a seven-membered ring transition state is predicted theoretically and shown experimentally to yield OH products. The low-lying ( tZZ ) conformer of 2-butenal oxide is identified based on computed anharmonic frequencies and intensities of its conformers. Experimental IR action spectra recorded in the fundamental CH stretch region with OH product detection by UV laser-induced fluorescence reveal a distinctive IR transition of the low-lying ( tZZ ) conformer at 2996 cm -1 that results in rapid unimolecular decay to OH products. Statistical RRKM calculations involving a combination of conformational isomerization and unimolecular decay via 1,6 H-transfer yield an effective decay rate k eff ( E ) on the order of 10 8 s -1 at ca. 3000 cm -1 in good accord with the experiment. Unimolecular decay proceeds with significant enhancement due to quantum mechanical tunneling. A rapid thermal decay rate of ca. 10 6 s -1 is predicted by master-equation modeling of 2-butenal oxide at 298 K, 1 bar. This novel unimolecular decay pathway is expected to increase the nonphotolytic production of OH radicals upon alkene ozonolysis in the troposphere.

The reactivity of carbonyl oxides has previously been shown to exhibit strong conformer and substituent dependencies. Through a combination of synchrotron-multiplexed photoionization mass spectrometry experiments (298 K and 4 Torr) and high-level theory [CCSD(T)-F12/cc-pVTZ-F12//B2PLYP-D3/cc-pVTZ with an added CCSDT(Q) correction], we explore the conformer dependence of the reaction of acetaldehyde oxide (CH 3 CHOO) with dimethylamine (DMA). The experimental data support the theoretically predicted 1,2-insertion mechanism and the formation of an amine-functionalized hydroperoxide reaction product. Tunable-vacuum ultraviolet photoionization probing of anti - or anti - + syn -CH 3 CHOO reveals a strong conformer dependence of the title reaction. The rate coefficient of DMA with anti -CH 3 CHOO is predicted to exceed that for the reaction with syn -CH 3 CHOO by a factor of ∼34,000, which is attributed to submerged barrier ( syn ) versus barrierless ( anti ) mechanisms for energetically downhill reactions.

Infrared (IR) action spectroscopy is utilized to characterize a prototypical carbon-centered hydroperoxyalkyl radical (•QOOH) transiently formed in the oxidation of volatile organic compounds. The •QOOH radical formed in isobutane oxidation, 2-hydroperoxy-2-methylprop-1-yl, •CH 2 (CH 3 ) 2 COOH, is generated in the laboratory by H-atom abstraction from tert-butyl hydroperoxide (TBHP). IR spectral features of jet-cooled and stabilized •QOOH radicals are observed from 2950 to 7050 cm -1 at energies that lie below and above the transition state barrier leading to OH radical and cyclic ether products. The observed •QOOH features include overtone OH and CH stretch transitions, combination bands involving OH or CH stretch and a lower frequency mode, and fundamental OH and CH stretch transitions. Most features arise from a single vibrational transition with band contours well simulated at a rotational temperature of 10 K. In each case, the OH products resulting from unimolecular decay of vibrationally activated •QOOH are detected by UV laser-induced fluorescence. Assignments of observed •QOOH IR transitions are guided by anharmonic frequencies computed using second order vibrational perturbation theory, a 2 + 1 model that focuses on the coupling of the OH stretch with two low-frequency torsions, as well as recently predicted statistical •QOOH unimolecular decay rates that include heavy-atom tunneling. Most of the observed vibrational transitions of •QOOH are readily distinguished from those of the TBHP precursor. The distinctive IR transitions of •QOOH, including the strong fundamental OH stretch, provide a general means for detection of •QOOH under controlled laboratory and real-world conditions.

Diastereomers have historically been ignored when building kinetic mechanisms for combustion. Low-temperature oxidation kinetics, which continues to gain interest in both combustion and atmospheric communities, may be affected by the inclusion of diastereomers in radical chain-branching pathways. In this work, key intermediates and transition states lacking stereochemical specification in an existing diethyl ether low-temperature oxidation mechanism were replaced with their diastereomeric counterparts. Rate coefficients for reactions involving diastereomers were computed with ab initio transition state theory master equation calculations. The presence of diastereomers increased rate coefficients by factors of 1.2-1.6 across various temperatures and pressures. Ignition delay simulations incorporating these revised rate coefficients indicate that the diastereomers enhanced the overall reactivity of the mechanism by almost 15% and increased the peak ketohydroperoxide concentration by 30% in the negative temperature coefficient region at combustion-relevant pressures. These results provide an illustrative indication of the important role of stereomeric effects in oxidation kinetics.

A prototypical hydroperoxyalkyl radical (•QOOH) intermediate, transiently formed in the oxidation of volatile organic compounds, was directly observed through its infrared fingerprint and energy-dependent unimolecular decay to hydroxyl radical and cyclic ether products. Direct time-domain measurements of •QOOH unimolecular dissociation rates over a wide range of energies were found to be in accord with those predicted theoretically using state-of-the-art electronic structure characterizations of the transition state barrier region. Unimolecular decay was enhanced by substantial heavy-atom tunneling involving O-O elongation and C-C-O angle contraction along the reaction pathway. Master equation modeling yielded a fully a priori prediction of the pressure-dependent thermal unimolecular dissociation rates for the •QOOH intermediate-again increased by heavy-atom tunneling-which are required for global models of atmospheric and combustion chemistry., (Copyright © 2021, American Association for the Advancement of Science.)

The important, and often dominant, role of tunneling in low temperature kinetics has resulted in numerous theoretical explorations into the methodology for predicting it. Nevertheless, there are still key aspects of the derivations that are lacking, particularly for nonseparable systems in the low temperature regime, and further explorations of the physical factors affecting the tunneling rate are warranted. In this work we obtain a closed-form rate expression for the tunneling rate constant that is a direct analog of the rigid-rotor-harmonic-oscillator expression. This expression introduces a novel "entanglement factor" that modulates the reaction rate. Furthermore, we are able to extend this expression, which is valid for nonseparable systems at low temperatures, to properly account for the conservation of angular momentum. In contrast, previous calculations have considered only vibrational transverse modes and so effectively employ a decoupled rotational partition function for the orientational modes. We also suggest a simple theoretical model to describe the tunneling effects in the vicinity of the crossover temperature (the temperature where tunneling becomes the dominating mechanism). This model allows one to naturally classify, interpret, and predict experimental data. Among other things, it quantitatively explains in simple terms the so-called "quantum bobsled" effect, also known as the negative centrifugal effect, which is related to curvature of the reaction path. Taken together, the expressions obtained here allow one to predict the thermal and E -resolved rate constants over broad ranges of temperatures and energies.

Methacrolein oxide (MACR-oxide) is a four-carbon, resonance-stabilized Criegee intermediate produced from isoprene ozonolysis, yet its reactivity is not well understood. This study identifies the functionalized hydroperoxide species, 1-hydroperoxy-2-methylallyl formate (HPMAF), generated from the reaction of MACR-oxide with formic acid using multiplexed photoionization mass spectrometry (MPIMS, 298 K = 25 °C, 10 torr = 13.3 hPa). Electronic structure calculations indicate the reaction proceeds via an energetically favorable 1,4-addition mechanism. The formation of HPMAF is observed by the rapid appearance of a fragment ion at m / z 99, consistent with the proposed mechanism and characteristic loss of HO 2 upon photoionization of functional hydroperoxides. The identification of HPMAF is confirmed by comparison of the appearance energy of the fragment ion with theoretical predictions of its photoionization threshold. The results are compared to analogous studies on the reaction of formic acid with methyl vinyl ketone oxide (MVK-oxide), the other four-carbon Criegee intermediate in isoprene ozonolysis.

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 (CH 3 C(O)CH 3 ) at a nominal temperature of 1800 K. The major product observed is ketene (CH 2 CO). Minor products identified include acetaldehyde (CH 3 CHO), propyne (CH 3 CCH), propene (CH 2 CHCH 3 ), 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.

Isoprene is the most abundant non-methane hydrocarbon emitted into the Earth's atmosphere. Ozonolysis is an important atmospheric sink for isoprene, which generates reactive carbonyl oxide species (R1R2C[double bond, length as m-dash]O+O-) known as Criegee intermediates. This study focuses on characterizing the catalyzed isomerization and adduct formation pathways for the reaction between formic acid and methyl vinyl ketone oxide (MVK-oxide), a four-carbon unsaturated Criegee intermediate generated from isoprene ozonolysis. syn-MVK-oxide undergoes intramolecular 1,4 H-atom transfer to form a substituted vinyl hydroperoxide intermediate, 2-hydroperoxybuta-1,3-diene (HPBD), which subsequently decomposes to hydroxyl and vinoxylic radical products. Here, we report direct observation of HPBD generated by formic acid catalyzed isomerization of MVK-oxide under thermal conditions (298 K, 10 torr) using multiplexed photoionization mass spectrometry. The acid catalyzed isomerization of MVK-oxide proceeds by a double hydrogen-bonded interaction followed by a concerted H-atom transfer via submerged barriers to produce HPBD and regenerate formic acid. The analogous isomerization pathway catalyzed with deuterated formic acid (D2-formic acid) enables migration of a D atom to yield partially deuterated HPBD (DPBD), which is identified by its distinct mass (m/z 87) and photoionization threshold. In addition, bimolecular reaction of MVK-oxide with D2-formic acid forms a functionalized hydroperoxide adduct, which is the dominant product channel, and is compared to a previous bimolecular reaction study with normal formic acid. Complementary high-level theoretical calculations are performed to further investigate the reaction pathways and kinetics.

Isoprene has the highest emission into Earth's atmosphere of any nonmethane hydrocarbon. Atmospheric processing of alkenes, including isoprene, via ozonolysis leads to the formation of zwitterionic reactive intermediates, known as Criegee intermediates (CIs). Direct studies have revealed that reactions involving simple CIs can significantly impact the tropospheric oxidizing capacity, enhance particulate formation, and degrade local air quality. Methyl vinyl ketone oxide (MVK-oxide) is a four-carbon, asymmetric, resonance-stabilized CI, produced with 21 to 23% yield from isoprene ozonolysis, yet its reactivity has not been directly studied. We present direct kinetic measurements of MVK-oxide reactions with key atmospheric species using absorption spectroscopy. Direct UV-Vis absorption spectra from two independent flow cell experiments overlap with the molecular beam UV-Vis-depletion spectra reported recently [M. F. Vansco, B. Marchetti, M. I. Lester, J. Chem. Phys. 149, 44309 (2018)] but suggest different conformer distributions under jet-cooled and thermal conditions. Comparison of the experimental lifetime herein with theory indicates only the syn -conformers are observed; anti -conformers are calculated to be removed much more rapidly via unimolecular decay. We observe experimentally and predict theoretically fast reaction of syn -MVK-oxide with SO 2 and formic acid, similar to smaller alkyl-substituted CIs, and by contrast, slow removal in the presence of water. We determine products through complementary multiplexed photoionization mass spectrometry, observing SO 3 and identifying organic hydroperoxide formation from reaction with SO 2 and formic acid, respectively. The tropospheric implications of these reactions are evaluated using a global chemistry and transport model., Competing Interests: Competing interest statement: R.L.C., C.A.T., and A. R. Ravishankara are amongst numerous coauthors in the General Discussion associated with the 2017 Faraday Discussion of Atmospheric Chemistry in the Anthropocene [Faraday Discuss. 200, 353–378 (2017)].

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 → H 2 + Be, H + NH 2 → H 2 + NH, CH 3 + C 2 H 5 → CH 4 + C 2 H 4 , H + C 2 H 5 → H 2 + C 2 H 4 , and H + HCO → H 2 + CO. Full configuration interaction (FCI) benchmark computations for the H + BeH, H + NH 2 , 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 CH 3 + C 2 H 5 , H + C 2 H 5 , and H + HCO reactions were computed at each level of theory with correlation-consistent basis sets (cc-pV X Z, 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.

The infrared (IR) action spectrum of the doubly substituted methyl-ethyl Criegee intermediate (MECI) is observed in the CH stretch overtone region with detection of OH products. The MECI exhibits four conformers, all of which undergo unimolecular decay via a 1,4 H-atom transfer mechanism, followed by the rapid release of OH products. Conformers with different orientations of the carbonyl oxide group with respect to the methyl and ethyl substituents (i.e., anti and syn) decay via distinct transition state barriers (16.1 kcal mol -1 and 15.4 kcal mol -1 , respectively). The observed IR action spectrum is in good agreement with the predicted anharmonic IR absorption spectrum, but exhibits significant congestion, which is attributed to couplings between spectroscopic bright states and nearby dark states. Energy-dependent OH appearance rates are measured upon IR excitation of the strongest features in the IR action spectrum and are found to be on the order of 10 6 -10 7 s -1 . The experimental rates are in good agreement with computed Rice-Ramsperger-Kassel-Marcus rates for the unimolecular decay of MECI at these energies, which incorporate quantum mechanical tunneling and sophisticated hindered rotor treatments, as well as high-level theoretical calculations of the TS barrier heights, rovibrational properties, and torsional barriers associated with the MECI conformers. Master equation modeling is used to predict thermal rates for the unimolecular decay of anti- and syn-MECI of 473 s -1 and 660 s -1 , respectively. Comparison with other previously studied Criegee intermediate systems provides insights into substituent effects on unimolecular decay under both energy-dependent and thermal conditions.

The 193-nm photolysis of CH 2 CHCN illustrates the capability of chirped-pulse Fourier transform millimeter-wave spectroscopy to characterize transition states. We investigate the HCN, HNC photofragments in highly excited vibrational states using both frequency and intensity information. Measured relative intensities of J = 1-0 rotational transition lines yield vibrational-level population distributions (VPD). These VPDs encode the properties of the parent molecule transition state at which the fragment molecule was born. A Poisson distribution formalism, based on the generalized Franck-Condon principle, is proposed as a framework for extracting information about the transition-state structure from the observed VPD. We employ the isotopologue CH 2 CDCN to disentangle the unimolecular 3-center DCN elimination mechanism from other pathways to HCN. Our experimental results reveal a previously unknown transition state that we tentatively associate with the HCN eliminated via a secondary, bimolecular reaction., Competing Interests: The authors declare no competing interest.

Ozonolysis of isoprene, one of the most abundant volatile organic compounds in the earth's atmosphere, generates the four-carbon unsaturated methacrolein oxide (MACR-oxide) Criegee intermediate. The first laboratory synthesis and direct detection of MACR-oxide is achieved through reaction of photolytically generated, resonance-stabilized iodoalkene radicals with oxygen. MACR-oxide is characterized on its first π* ← π electronic transition using a ground-state depletion method. MACR-oxide exhibits a broad UV-visible spectrum peaked at 380 nm with weak oscillatory structure at long wavelengths ascribed to vibrational resonances. Complementary theory predicts two strong π* ← π transitions arising from extended conjugation across MACR-oxide with overlapping contributions from its four conformers. Electronic promotion to the 1 1 ππ* state agrees well with experiment, and results in nonadiabatic coupling and prompt release of O 1 D products observed as anisotropic velocity-map images. This UV-visible detection scheme will enable study of its unimolecular and bimolecular reactions under thermal conditions of relevance to the atmosphere.

Quasiclassical trajectories are used to compute nonthermal rate constants, k * , for abstraction reactions involving highly-excited methane CH 4 * and the radicals H, O, OH, and O 2 . Several temperatures and internal energies of methane, E vib , are considered, and significant nonthermal rate enhancements for large E vib are found. Specifically, when CH 4 * is internally excited close to its dissociation threshold (E vib ≈ D 0 = 104 kcal/mol), its reactivity with H, O, and OH is shown to be collision-rate-limited and to approach that of comparably-sized radicals, such as CH 3 , with k * > 10 -10 cm 3 molecule -1 s -1 . Rate constants this large are more typically associated with barrierless reactions, and at 1000 K, this represents a nonthermal rate enhancement, k * /k, of more than two orders of magnitude relative to thermal rate constants k. We show that large nonthermal rate constants persist even after significant internal cooling, with k * /k > 10 down to E vib ≈ D 0 /4. The competition between collisional cooling and nonthermal reactivity is studied using a simple model, and nonthermal reactions are shown to account for up to 35%-50% of the fate of the products of H + CH 3 = CH 4 * under conditions of practical relevance to combustion. Finally, the accuracy of an effective temperature model for estimating k * from k is quantified.

A priori rate predictions for gas phase reactions have undergone a gradual but dramatic transformation, with current predictions often rivaling the accuracy of the best available experimental data. The utility of such kinetic predictions would be greatly magnified if they could more readily be implemented for large numbers of systems. Here, we report the development of a new computational environment, namely, EStokTP, that reduces the human effort involved in the rate prediction for single channel reactions essentially to the specification of the methodology to be employed. The code can also be used to obtain all the necessary master equation building blocks for more complex reactions. In general, the prediction of rate constants involves two steps, with the first consisting of a set of electronic structure calculations and the second in the application of some form of kinetic solver, such as a transition state theory (TST)-based master equation solver. EStokTP provides a fully integrated treatment of both steps through calls to external codes to perform first the electronic structure and then the master equation calculations. It focuses on generating, extracting, and organizing the necessary structural properties from a sequence of calls to electronic structure codes, with robust automatic failure recovery options to limit human intervention. The code implements one or multidimensional hindered rotor treatments of internal torsional modes (with automated projection from the Hessian and with optional vibrationally adiabatic corrections), Eckart and multidimensional tunneling models (such as small curvature theory), and variational treatments (based on intrinsic reaction coordinate following). This focus on a robust implementation of high-level TST methods allows the code to be used in high accuracy studies of large sets of reactions, as illustrated here through sample studies of a few hundred reactions. At present, the following reaction types are implemented in EStokTP: abstraction, addition, isomerization, and beta-decomposition. Preliminary protocols for treating barrierless reactions and multiple-well and/or multiple-channel potential energy surfaces are also illustrated.

The reaction of ozone with isoprene, one of the most abundant volatile organic compounds in the atmosphere, produces three distinct carbonyl oxide species (RR'COO) known as Criegee intermediates: formaldehyde oxide (CH 2 OO), methyl vinyl ketone oxide (MVK-OO), and methacrolein oxide (MACR-OO). The nature of the substituents (R,R' = H, CH 3 , CH═CH 2 ) and conformations of the Criegee intermediates control their subsequent chemistry in the atmosphere. In particular, unimolecular decay of MVK-OO is predicted to be the major source of hydroxyl radicals (OH) in isoprene ozonolysis. This study reports the initial laboratory synthesis and direct detection of MVK-OO through reaction of a photolytically generated, resonance-stabilized monoiodoalkene radical with O 2 . MVK-OO is characterized utilizing infrared (IR) action spectroscopy, in which IR activation of MVK-OO with two quanta of CH stretch at ca. 6000 cm -1 is coupled with ultraviolet detection of the resultant OH products. MVK-OO is identified by comparison of the experimentally observed IR spectral features with theoretically predicted IR absorption spectra. For syn-MVK-OO, the rate of appearance of OH products agrees with the unimolecular decay rate predicted using statistical theory with tunneling. This validates the hydrogen atom transfer mechanism and computed transition-state barrier (18.0 kcal mol -1 ) leading to OH products. Theoretical calculations reveal an additional roaming pathway between the separating radical fragments, which results in other products. Master equation modeling yields a thermal unimolecular decay rate for syn-MVK-OO of 33 s -1 (298 K, 1 atm). For anti-MVK-OO, theoretical exploration of several unimolecular decay pathways predicts that isomerization to dioxole is the most likely initial step to products.

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-C 2 H 4 O 3 , 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 (H 2 COO) and formaldehyde (H 2 CO). 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.

Alkylated aromatics constitute a significant fraction of the components commonly found in commercial fuels. Toluene is typically considered as a reference fuel. Together with n-heptane and iso-octane, it allows for realistic emulations of the behavior of real fuels by the means of surrogate mixture formulations. Moreover, it is a key precursor for the formation of poly-aromatic hydrocarbons, which are of relevance to understanding soot growth and oxidation mechanisms. In this study the POLIMI kinetic model is first updated based on the literature and on recent kinetic modelling studies of toluene pyrolysis and oxidation. Then, important reaction pathways are investigated by means of high-level theoretical methods, thereby advancing the present knowledge on toluene oxidation. H-Abstraction reactions by OH, HO2, O and O2, and the reactivity on the multi well benzyl-oxygen (C6H5CH2 + O2) potential energy surface (PES) were investigated using electronic structure calculations, transition state theory in its conventional, variational, and variable reaction coordinate forms (VRC-TST), and master equation calculations. Exploration of the effect on POLIMI model performance of literature rate constants and of the present calculations provides valuable guidelines for implementation of the new rate parameters in existing toluene kinetic models.

Alkene ozonolysis is a primary oxidation pathway for anthropogenic and biogenic alkenes emitted into the troposphere. It is also an important source of atmospheric hydroxyl (OH) radicals, often called the atmosphere's detergent. Alkene ozonolysis takes place through a highly exothermic reaction pathway with multiple intermediates and barriers prior to releasing the OH radical products. This Account focuses on a key reaction intermediate with a carbonyl oxide functional group (-COO), known as the Criegee intermediate, which is formed along with a carbonyl coproduct in alkene ozonolysis reactions. Under atmospheric conditions, the initially energized Criegee intermediates may promptly decay to OH products or be collisionally stabilized prior to thermal decay to OH radicals and other products. Alternatively, the stabilized Criegee intermediates may undergo bimolecular reactions with atmospheric species, including water vapor and sulfur dioxide, which can lead to nucleation and growth of aerosols. The dimethyl-substituted Criegee intermediate, (CH 3 ) 2 COO, is utilized in this Account to showcase recent efforts to experimentally measure and theoretically predict the rates for prompt and thermal unimolecular decay processes of prototypical Criegee intermediates under laboratory and atmospheric conditions. The experimental laboratory studies utilize an alternative synthesis method to efficiently generate Criegee intermediates via the reaction of iodoalkyl radicals with O 2 . Infrared excitation is then used to prepare the (CH 3 ) 2 COO Criegee intermediates at specific energies in the vicinity of the transition state barrier or significantly below the barrier for 1,4-hydrogen transfer that leads to OH products. The rate of unimolecular decay is revealed through direct time-domain measurements of the appearance of OH products utilizing ultraviolet laser-induced fluorescence detection under collision-free conditions. Complementary high-level theoretical calculations are carried out to evaluate the transition state barrier and the energy-dependent unimolecular decay rates for (CH 3 ) 2 COO using Rice-Ramsperger-Kassel-Marcus (RRKM) theory, which are in excellent accord with the experimental measurements. Quantum mechanical tunneling through the barrier, incorporated through Eckart and semiclassical transition state theory models, is shown to make a significant contribution to the unimolecular decay rates at energies in the vicinity of and much below the barrier. Master equation modeling is used to extend the energy-dependent unimolecular rates to thermal decay rates of (CH 3 ) 2 COO under tropospheric conditions (high pressure limit), which agree well with recent laboratory measurements [ Smith et al. J. Phys. Chem. A 2016 , 120 , 4789 and Chhantyal-Pun et al. J. Phys. Chem. A 2017 , 121 , 4 - 15 ]. Again, tunneling is shown to enhance the thermal decay rate by orders of magnitude. The experimentally validated unimolecular rates are also utilized in modeling the prompt and thermal unimolecular decay of chemically activated (CH 3 ) 2 COO formed upon ozonolysis of 2,3-dimethyl-2-butene under atmospheric conditions [ Drozd et al. J. Phys. Chem. A 2017 , 121 , 6036 - 6045 ]. Future challenges lie in extension of these spectroscopic and dynamical methods to Criegee intermediates derived from more complex ozonolysis reactions involving biogenic alkenes.

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.

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 CH 2 CHCN, 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.

Ozonolysis of alkenes, an important nonphotolytic source of hydroxyl (OH) radicals in the atmosphere, proceeds through unimolecular decay of Criegee intermediates. Here, we report a large kinetic isotope effect associated with the rate-limiting hydrogen-transfer step that releases OH radicals for a prototypical Criegee intermediate, CH 3 CHOO. IR excitation of selectively deuterated syn -CD 3 CHOO is shown to result in deuterium atom transfer and release OD radical products. Vibrational activation of syn -CD 3 CHOO is coupled with direct time-resolved detection of OD products to measure a 10-fold slower rate of unimolecular decay upon deuteration in the vicinity of the transition state barrier, which is confirmed by microcanonical statistical theory that incorporates quantum mechanical tunneling. The corresponding kinetic isotope effect of ∼10 is attributed primarily to the decreased probability of D-atom vs. H-atom transfer arising from tunneling. Master equation modeling is utilized to compute the thermal unimolecular decay rates for selectively and fully deuterated syn methyl-substituted Criegee intermediates under atmospheric conditions. At 298 K (1 atm), tunneling is predicted to enhance the thermal decay rate of syn -CH 3 CHOO compared with the deuterated species, giving rise to a significant kinetic isotope effect of ∼50., Competing Interests: The authors declare no conflict of interest.

Termolecular association reactions involve ephemeral collision complexes-formed from the collision of two molecules-that collide with a third and chemically inert 'bath gas' molecule that simply transfers energy to/from the complex. These collision complexes are generally not thought to react chemically on collision with a third molecule in the gas-phase systems of combustion and planetary atmospheres. Such 'chemically termolecular' reactions, in which all three molecules are involved in bond making and/or breaking, were hypothesized long ago in studies establishing radical chain branching mechanisms, but were later concluded to be unimportant. Here, with data from ab initio master equation and kinetic-transport simulations, we reveal that reactions of H + O 2 collision complexes with other radicals constitute major kinetic pathways under common combustion situations. These reactions are also found to influence flame propagation speeds, a common measure of global reactivity. Analogous chemically termolecular reactions mediated by ephemeral collision complexes are probably of significance in various combustion and planetary environments.

The oxidation chemistry of the simplest conjugated hydrocarbon, 1,3-butadiene, can provide a first step in understanding the role of polyunsaturated hydrocarbons in combustion and, in particular, an understanding of their contribution toward soot formation. On the basis of our previous work on propene and the butene isomers (1-, 2-, and isobutene), it was found that the reaction kinetics of Ḣ-atom addition to the C═C double bond plays a significant role in fuel consumption kinetics and influences the predictions of high-temperature ignition delay times, product species concentrations, and flame speed measurements. In this study, the rate constants and thermodynamic properties for Ḣ-atom addition to 1,3-butadiene and related reactions on the Ċ 4 H 7 potential energy surface have been calculated using two different series of quantum chemical methods and two different kinetic codes. Excellent agreement is obtained between the two different kinetics codes. The calculated results including zero-point energies, single-point energies, rate constants, barrier heights, and thermochemistry are systematically compared among the two quantum chemical methods. 1-Methylallyl (Ċ 4 H 7 1-3) and 3-buten-1-yl (Ċ 4 H 7 1-4) radicals and C 2 H 4 + Ċ 2 H 3 are found to be the most important channels and reactivity-promoting products, respectively. We calculated that terminal addition is dominant (>80%) compared to internal Ḣ-atom addition at all temperatures in the range 298-2000 K. However, this dominance decreases with increasing temperature. The calculated rate constants for the bimolecular reaction C 4 H 6 + Ḣ → products and C 2 H 4 + Ċ 2 H 3 → products are in excellent agreement with both experimental and theoretical results from the literature. For selected C 4 species, the calculated thermochemical values are also in good agreement with literature data. In addition, the rate constants for H atom abstraction by Ḣ atoms have also been calculated, and it is found that abstraction from the central carbon atoms is the dominant channel (>70%) at temperatures in the range of 298-2000 K. Finally, by incorporating our calculated rate constants for both Ḣ atom addition and abstraction into our recently developed 1,3-butadiene model, we show that laminar flame speed predictions are significantly improved, emphasizing the value of this study.

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 H 2 , CH 4 , H 2 O, and NH 3 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.

The NCN diradical is an important intermediate of prompt nitric oxide formation in flames. The mechanism of intersystem crossing (ISC) in the NCN molecule formed via pyrolysis or photolysis of NCN 3 is of relevance to the interpretation of experiments that utilize NCN 3 as a precursor for laboratory studies of NCN kinetics. This mechanism has been investigated by means of multi-reference configuration interaction calculations. From the potential energy surfaces for NCN 3 dissociation, it was inferred that both thermal and photo-chemical decomposition initially lead to NCN in its lowest singlet state, ã 1 Δ g , with a possible contribution from the b̃ 1 Σ g + state at low photolysis wavelengths. Direct formation of the triplet ground state X̃  3 Σ g - is also feasible for the photolytic pathway. An analysis of surface crossings between ã or b̃ and the triplet ground state X̃  3 Σ g - in the absence and presence of a helium atom revealed an ISC channel NCN1(ã)→ 3 NCN(X̃) via a strongly bent structure. However, its barrier of 38 kcal mol -1 relative to the singlet minimum turned out to be much too high to explain the fast ISC observed in experiments. A rigid-bender model including Renner-Teller interactions was used to examine the occurrence of mixed-multiplicity rovibrational states-so-called gateway states-that could enhance collision-induced ISC. The results of this study indicate that a gateway mechanism is probably not operative in the case of the ã/X̃ pair of states in NCN.