ABSTRACT:The HO 2 + HO 2 → H 2 O 2 + O 2 chemical reaction is studied using statistical rate theory in conjunction with high level ab initio electronic structure calculations. A new theoretical rate coefficient is generated that is appropriate for both high and low temperature regimes. The transition state region for the ground triplet potential energy surface is characterized using the CASPT2/CBS/aug-cc-pVTZ method with 14 active electrons and 10 active orbitals. The reaction is found to proceed through an intermediate complex bound by approximately 9.79 kcal/mol. There is no potential barrier in the entrance channel, although the free energy barrier was determined using a large Monte Carlo sampling of the HO 2 orientations. The inner (tight) transition state lies below the entrance threshold. It is found that this inner transition state exhibits two saddle points corresponding to torsional conformations of the complex. A unified treatment based on vibrational adiabatic theory is presented that permits the reaction to occur on an equal footing for any value of the torsional angle. The quantum tunneling is also reformulated based on this new approach. The rate coefficient obtained is in good agreement with low temperature experimental results but is significantly lower than the results of shock tube experiments for high temperatures.

We present a model of the surface kinetics of the dehydrogenation reaction of methanol on the Pd(111), Pt(111), and Ni(111) metal surfaces. The mechanism consists of 10 reversible dehydrogenation reactions that lead to the final products of CO and H2. The rate coefficients for each step are calculated using ab initio transition state theory that employs a new approach to obtain the symmetry factors. The potential energies and frequencies of the reagents and transition states are computed using plane wave DFT with the PW91 exchange correlation functional. The mechanism is investigated for low coverages using a global sensitivity analysis that monitors the response of a target function of the kinetics to the value of the rate coefficients. On Pd(111) and Ni(111), the reaction COH → CO + H is found to be rate limiting, and overall rates are highly dependent upon the decomposition time of the COH intermediate. Reactions at branches in the reaction network are also particularly important in the kinetics. A stochastic atom-following approach to pathway analysis is used to elucidate both the pathway probabilities in the kinetics and the dependence of the pathways on the values of the key rate coefficients of the mechanisms. On Pd(111) and Ni(111) there exists significant competition between the pathway containing the slow step and faster pathways that bypass the slow step. A discussion is given of the dependence of the model target’s probability density function on the chemical pathways.

We study the role of individual reaction rates on engine performance, with an emphasis on the contribution of quantum tunneling. It is demonstrated that the effect of quantum tunneling corrections for the reaction HO2 + HO2 = H2O2 + O2 can have a noticeable impact on the performance of a high-fidelity model of a compression-ignition (e.g., diesel) engine, and that an accurate prediction of ignition delay time for the engine model requires an accurate estimation of the tunneling correction for this reaction. The three-dimensional model includes detailed descriptions of the chemistry of a surrogate for a biodiesel fuel, as well as all the features of the engine, such as the liquid fuel spray and turbulence. This study is part of a larger investigation of how the features of the dynamics and potential energy surfaces of key reactions, as well as their reaction rate uncertainties, affect engine performance, and results in these directions are also presented here.

A model for the combustion of butanol is studied using a recently developed theoretical method for the systematic improvement of the kinetic mechanism. The butanol mechanism includes 1446 reactions, and we demonstrate that it is straightforward and computationally feasible to implement a full global sensitivity analysis incorporating all the reactions. In addition, we extend our previous analysis of ignition-delay targets to include species targets. The combination of species and ignition targets leads to multitarget global sensitivity analysis, which allows for a more complete mechanism validation procedure than we previously implemented. The inclusion of species sensitivity analysis allows for a direct comparison between reaction pathway analysis and global sensitivity analysis.