Search

Recent high-level theoretical calculations predict a mild temperature dependence for HO 2 + HO 2 inconsistent with state-of-the-art experimental determinations that upheld the stronger temperature dependence observed in early experiments. Via MultiScale Informatics analysis of the theoretical and experimental data, we identified an alternative interpretation of the raw experimental data that uses HO 2 + HO 2 rate constants nearly identical to theoretical predictions─implying that the theoretical and experimental data are actually consistent, at least when considering the raw data from experimental studies. Similar analyses of typical signals from low-temperature experiments indicate that an HOOOOH intermediate─identified by recent theory but absent from earlier interpretations─yields modest effects that are smaller than, but may have contributed to, the scatter in data among different experiments. More generally, the findings demonstrate that modern chemical theories and experiments have progressed to a point where meaningful comparison requires joint consideration of their data simultaneously.

Rovibrationally excited ephemeral complexes AB ** , formed from the association of two molecules A + B, are generally considered to undergo collisions only with an inert bath gas M that transfer energy-inducing termolecular association reactions A + B (+M) → AB (+M). Recent studies have demonstrated that reactive collisions of AB ** with a third molecule C-inducing chemically termolecular reactions A + B + C → products-can also be significant in combustion and planetary atmospheres. Previous studies on systems with reactive collisions have primarily focused on limited ranges of reactive collider mole fraction, X C , and pressure, P , specific to the chosen application. Yet, it remains to be established how such systems, and the rate constants of their emergent phenomenological reactions, behave over the wide X C and P ranges of potential interest-a gap in the present understanding that has impeded the development of broadly applicable rate laws and general treatment of such systems in kinetic modeling. Here, we present results from master equation calculations for HO 2 ** formed from H + O 2 and its reactions with H to advance understanding and explore representations of systems with reactive colliders across wide ranges of X C and P . With regard to understanding, we demonstrate that reactive collisions can both (1) increase the overall rate of conversion of reactants to products and (2) alter the branching ratio among final products. With regard to representations in kinetic models, we find that rate constants of all emergent phenomenological reactions-termolecular association A + B (+M), chemically termolecular A + B + C, and bimolecular AB + C-exhibit a rich X C and P dependence. We also present analyses to explore the existence of a unique phenomenological representation (or lack thereof) and assess ways for the distinct effects of reactive collisions to be represented in kinetic models.

Nearly all studies of and available data for pressure-dependent reactions focus on pure bath gases. Of the comparatively fewer studies on bath gas mixtures, important to combustion and planetary atmospheres, nearly all focus on single-channel reactions. The present study explores, and seeks reliable representations of, bath gas mixture effects on multichannel reactions. Analytical and numerical solutions of the master equation here reveal several unique manifestations of mixture effects for multichannel reactions, including behavior completely opposite to trends observed for single-channel reactions. The most common way of evaluating mixture rate constants from data for pure components, the classic linear mixture rule, is found to yield errors exceeding a factor of ∼10. A new linear mixture rule based on the reduced pressure, instead of the absolute pressure, is found to be accurate within ∼30% for rate constants (and ∼50% for the branching ratio). A new nonlinear mixture rule that additionally incorporates analytically derived activity coefficients is found to be accurate within ∼10% for rate constants and branching ratios. These new mixture rules are therefore recommended for use in fundamental and applied chemical kinetics investigations of reacting mixtures, including reacting flow codes and experimental interpretations of third-body efficiencies.

Low-temperature propane oxidation was studied at P = 4 Torr and T = 530, 600, and 670 K by time-resolved multiplexed photoionization mass spectrometry (MPIMS), which probes the reactants, intermediates, and products with isomeric selectivity using tunable synchrotron vacuum UV ionizing radiation. The oxidation is initiated by pulsed laser photolysis of oxalyl chloride, (COCl)2, at 248 nm, which rapidly generates a ∼1:1 mixture of 1-propyl (n-propyl) and 2-propyl (i-propyl) radicals via the fast Cl + propane reaction. At all three temperatures, the major stable product species is propene, formed in the propyl + O2 reactions by direct HO2 elimination from both n- and i-propyl peroxy radicals. The experimentally derived propene yields relative to the initial concentration of Cl atoms are (20 ± 4)% at 530 K, (55 ± 11)% at 600 K, and (86 ± 17)% at 670 K at a reaction time of 20 ms. The lower yield of propene at low temperature reflects substantial formation of propyl peroxy radicals, which do not completely decompose on the experimental time scale. In addition, C3H6O isomers methyloxirane, oxetane, acetone, and propanal are detected as minor products. Our measured yields of oxetane and methyloxirane, which are coproducts of OH radicals, suggest a revision of the OH formation pathways in models of low-temperature propane oxidation. The experimental results are modeled and interpreted using a multiscale informatics approach, presented in detail in a separate publication (Burke, M. P.; Goldsmith, C. F.; Klippenstein, S. J.; Welz, O.; Huang H.; Antonov I. O.; Savee J. D.; Osborn D. L.; Zádor, J.; Taatjes, C. A.; Sheps, L. Multiscale Informatics for Low-Temperature Propane Oxidation: Further Complexities in Studies of Complex Reactions. J. Phys. Chem A. 2015, DOI: 10.1021/acs.jpca.5b01003). The model predicts the time profiles and yields of the experimentally observed primary products well, and shows satisfactory agreement for products formed mostly via secondary radical-radical reactions.

The present paper describes further development of the multiscale informatics approach to kinetic model formulation of Burke et al. (Burke, M. P.; Klippenstein, S. J.; Harding, L. B. Proc. Combust. Inst. 2013, 34, 547-555) that directly incorporates elementary kinetic theories as a means to provide reliable, physics-based extrapolation of kinetic models to unexplored conditions. Here, we extend and generalize the multiscale informatics strategy to treat systems of considerable complexity-involving multiwell reactions, potentially missing reactions, nonstatistical product branching ratios, and non-Boltzmann (i.e., nonthermal) reactant distributions. The methodology is demonstrated here for a subsystem of low-temperature propane oxidation, as a representative system for low-temperature fuel oxidation. A multiscale model is assembled and informed by a wide variety of targets that include ab initio calculations of molecular properties, rate constant measurements of isolated reactions, and complex systems measurements. Active model parameters are chosen to accommodate both "parametric" and "structural" uncertainties. Theoretical parameters (e.g., barrier heights) are included as active model parameters to account for parametric uncertainties in the theoretical treatment; experimental parameters (e.g., initial temperatures) are included to account for parametric uncertainties in the physical models of the experiments. RMG software is used to assess potential structural uncertainties due to missing reactions. Additionally, branching ratios among product channels are included as active model parameters to account for structural uncertainties related to difficulties in modeling sequences of multiple chemically activated steps. The approach is demonstrated here for interpreting time-resolved measurements of OH, HO2, n-propyl, i-propyl, propene, oxetane, and methyloxirane from photolysis-initiated low-temperature oxidation of propane at pressures from 4 to 60 Torr and temperatures from 300 to 700 K. In particular, the multiscale informed model provides a consistent quantitative explanation of both ab initio calculations and time-resolved species measurements. The present results show that interpretations of OH measurements are significantly more complicated than previously thought-in addition to barrier heights for key transition states considered previously, OH profiles also depend on additional theoretical parameters for R + O2 reactions, secondary reactions, QOOH + O2 reactions, and treatment of non-Boltzmann reaction sequences. Extraction of physically rigorous information from those measurements may require more sophisticated treatment of all of those model aspects, as well as additional experimental data under more conditions, to discriminate among possible interpretations and ensure model reliability.

We consider an alternative formulation of the master equation for complex-forming chemical reactions with multiple wells and bimolecular products. Within this formulation the dynamical phase space consists of only the microscopic populations of the various isomers making up the reactive complex, while the bimolecular reactants and products are treated equally as sources and sinks. This reformulation yields compact expressions for the phenomenological rate coefficients describing all chemical processes, i.e., internal isomerization reactions, bimolecular-to-bimolecular reactions, isomer-to-bimolecular reactions, and bimolecular-to-isomer reactions. The applicability of the detailed balance condition is discussed and confirmed. We also consider the situation where some of the chemical eigenvalues approach the energy relaxation time scale and show how to modify the phenomenological rate coefficients so that they retain their validity.