High-pressure falloff curves and specific rate constants for the reactions atomic hydrogen + molecular oxygen .dblharw. perhydroxyl .dblharw. hydroxyl + atomic oxygen

Molecular Dynamics of Combustion Reactions in Supercritical Carbon Dioxide. 6. Computational Kinetics of Reactions between Hydrogen Atom and Oxygen Molecule H + O2 ⇌ HO + O and H + O2 ⇌ HO2

Reactions of the hydrogen atom and the oxygen molecule are among the most important ones in the hydrogen and hydrocarbon oxidation mechanisms, including combustion in a supercritical CO₂ (sCO₂) environment, known as oxy-combustion or the Allam cycle. Development of these energy technologies requires understanding of chemical kinetics of H + O₂ ⇌ HO + O and H + O₂ ⇌ HO₂ in high pressures and concentrations of CO₂. Here, we combine quantum treatment of the reaction system by the transition state theory with classical molecular dynamics simulation and the multistate empirical valence bonding method to treat environmental effects. Potential of mean force in the sCO₂ solvent at various temperatures 1000-2000 K and pressures 100-400 atm was obtained. The reaction rate for H + O₂ ⇌ HO + O was found to be pressure-independent and described by the extended Arrhenius equation 4.23 × 10^-7 T^-0.73 exp(-21 855.2 cal/mol/RT) cm³/molecule/s, while the reaction rate H + O₂ ⇌ HO₂ is pressure-dependent and can be expressed as 5.22 × 10^-2 T^-2.86 exp(-7247.4 cal/mol/RT) cm³/molecule/s at 300 atm.

Consistent Conditional Moment Closure Modelling of a Lifted Turbulent Jet Flame Using the Presumedβ-PDF Approach

A liftedH2/N2turbulent jet flame issuing into a vitiated coflow is investigated using the conditional moment closure. The conditional velocity (CV) and the conditional scalar dissipation rate (CSDR) submodels are chosen such that they are fully consistent with the moments of the presumedβprobability density function (PDF). The CV is modelled using the PDF-gradient diffusion model. Two CSDR submodels based on the double integration of the homogeneous and inhomogeneous mixture fraction PDF transport equations are implemented. The effect of CSDR modelling is investigated over a range of coflow temperatures (Tc) and the stabilisation mechanism is determined from the analysis of the transport budgets and the history of radical build-up ahead of the stabilisation height. For allTc, the balance between chemistry, axial convection, and micromixing, and the absence of axial diffusion upstream of the stabilisation height indicate that the flame is stabilized by autoignition. This conclusion is confirmed from the rapid build-up ofHO2ahead ofH,O, andOH. The inhomogeneous CSDR modelling yields higher dissipation levels at the most reactive mixture fraction, which results in longer ignition delays and larger liftoff heights. The effect of the spurious sources arising from homogeneous modelling is found to be small but nonnegligible, mostly notably within the flame zone.

Study of the H+O2 reaction by means of quantum mechanical and statistical approaches: the dynamics on two different potential energy surfaces

The possible existence of a complex-forming pathway for the H+O(2) reaction has been investigated by means of both quantum mechanical and statistical techniques. Reaction probabilities, integral cross sections, and differential cross sections have been obtained with a statistical quantum method and the mean potential phase space theory. The statistical predictions are compared to exact results calculated by means of time dependent wave packet methods and a previously reported time independent exact quantum mechanical approach using the double many-body expansion (DMBE IV) potential energy surface (PES) [Pastrana et al., J. Phys. Chem. 94, 8073 (1990)] and the recently developed surface (denoted XXZLG) by Xu et al. [J. Chem. Phys. 122, 244305 (2005)]. The statistical approaches are found to reproduce only some of the exact total reaction probabilities for low total angular momenta obtained with the DMBE IV PES and some of the cross sections calculated at energy values close to the reaction threshold for the XXZLG surface. Serious discrepancies with the exact integral cross sections at higher energy put into question the possible statistical nature of the title reaction. However, at a collision energy of 1.6 eV, statistical rotationally resolved cross sections managed to reproduce the experimental cross sections for the H+O(2)(v=0,j=1)-->OH(v(')=1,j('))+O process reasonably well.

Reaction of O2 with the Hydrogen Atom in Water up to 350 °C

The reaction of the H* atom with O2, giving the hydroperoxyl HO2* radical, has been investigated in pressurized water up to 350 degrees C using pulse radiolysis and deep-UV transient absorption spectroscopy. The reaction rate behavior is highly non-Arrhenius, with near diffusion-limited behavior at room temperature, increasing to a near constant limiting value of approximately 5 x 10(10) M(-1) s(-1) above 250 degrees C. The high-temperature rate constant is in near-perfect agreement with experimental extrapolations and ab initio calculations of the gas-phase high-pressure limiting rate. As part of the study, reaction of the OH* radical with H2 has been reevaluated at 350 degrees C, giving a rate constant of (6.0 +/- 0.5) x 10(8) M(-1) s(-1). The mechanism of the H* atom reaction with the HO2* radical is also investigated and discussed.

Quantum Mechanical Rate Constants for H + O2 ↔ O + OH and H + O2 → HO2 Reactions

Canonical rate constants for both the forward and reverse H + O(2) <--> O + OH reactions were calculated using a quantum wave packet-based statistical model on the DMBE IV potential energy surface of Varandas and co-workers. For these bimolecular reactions, the results show reasonably good agreement with available experimental and theoretical data up to 1500 K. In addition, the capture rate for the H + O(2) --> HO(2) addition reaction at the high-pressure limit was obtained on the same potential using a time-independent quantum capture method. Excellent agreement with experimental and quasi-classical trajectory results was obtained except for at very low temperatures, where a reaction threshold was found and attributed to the centrifugal barrier of the orbital motion.

A Quasiclassical Trajectory Study of the Reaction H + O2 ⇔ OH + O with the O2 Reagent Vibrationally Excited

Quasiclassical trajectories have been computed on the Melius-Blint (MB) Potential Energy Surface (PES) and on the Double Many-Body Expansion (DMBE) IV PES of Pastrana et al. describing the H + O(2) <==> OH + O reaction with the nonrotating (J = 0) O(2) reagent vibrationally excited to levels v = 6, 7, 8, 9, and 10 at four temperatures: 1000, 2000, 3000, and 4000 K. The vibrational energy levels were selected by using a semiclassical Einstein-Brillouin-Keller (EBK) quantization procedure while the relative translational energy was sampled from a Boltzmann weighted distribution. The rate coefficient for the formation of the OH + O products is seen to increase monotonically with quantum number and nearly monotonically with temperature. On the MB PES, at T = 1000 K, the total rate coefficient increases by a factor of 5.2 as the initial vibrational quantum number of the O(2) diatom increases from v = 6 to v = 10. For T = 2000 K, this factor drops to 3.3, to 2.9 for T = 3000 K, and to 2.5 for T = 4000 K. On the DMBE IV PES, at T = 1000 K the total rate coefficient increases by a factor of 4.1 as the initial vibrational quantum number of the O(2) diatom increases from v = 6 to v = 10. For T = 2000 K, this factor drops to 3.5, to 2.1 for T = 3000 K, and to 2.0 for T = 4000 K. The less-direct group (defined below) of trajectories is sensitive to the initial O(2) vibrational excitation in several different temperature ranges, apparently retaining the effect of reagent vibrational excitation. The more-direct group (defined below) of trajectories does not exhibit this behavior. Reagent vibrational excitation does not increase the total rate coefficients for the title reaction more than the increase due to a simple temperature increase. The less-direct and more-direct groups of trajectories differ in their contribution to the rate coefficient for the title reaction. In particular, at T = 4000 K, the two PESs used in this work differ dramatically in the roles of the less-direct and more-direct trajectories. The behavior of the more-direct and less-direct groups of trajectories can be understood in terms of the efficiency of intramolecular vibrational energy transfer. This work utilizes the recently introduced PES Library, POTLIB 2001, which made the comparisons between the two PESs discussed in this work possible in a very straightforward way.