Reflected shock tube studies of high-temperature rate constants for OH + NO2 --> HO2 + NO and OH + HO2 --> H2O + O2

The motivation for the present study comes from the preceding paper where it is suggested that accepted rate constants for OH + NO2 --> NO + HO2 are high by approximately 2. This conclusion was based on a reevaluation of heats of formation for HO2, OH, NO, and NO2 using the Active Thermochemical Table (ATcT) approach. The present experiments were performed in C2H5I/NO2 mixtures, using the reflected shock tube technique and OH-radical electronic absorption detection (at 308 nm) and using a multipass optical system. Time-dependent profile decays were fitted with a 23-step mechanism, but only OH + NO2, OH + HO2, both HO2 and NO2 dissociations, and the atom molecule reactions, O + NO2 and O + C2H4, contributed to the decay profile. Since all of the reactions except the first two are known with good accuracy, the profiles were fitted by varying only OH + NO2 and OH + HO2. The new ATcT approach was used to evaluate equilibrium constants so that back reactions were accurately taken into account. The combined rate constant from the present work and earlier work by Glaenzer and Troe (GT) is k(OH+NO2) = 2.25 x 10(-11) exp(-3831 K/T) cm3 molecule(-1) s(-1), which is a factor of 2 lower than the extrapolated direct value from Howard but agrees well with NO + HO2 --> OH + NO2 transformed with the updated equilibrium constants. Also, the rate constant for OH + HO2 suitable for combustion modeling applications over the T range (1200-1700 K) is (5 +/- 3) x 10(-11) cm3 molecule(-1) s(-1). Finally, simulating previous experimental results of GT using our updated mechanism, we suggest a constant rate for k(HO2+NO2) = (2.2 +/- 0.7) x 10(-11) cm3 molecule(-1) s(-1) over the T range 1350-1760 K.     

Study of the kinetics of the reactions of NO3 with HO2 and OH

The Absolute rate constants for the gas-phase reactions of NO₃ with HO₂ and OH have been determined using the discharge flow laser magnetic resonance method (DF-LMR). Since OH was found to be produced in the reaction of HO₂ with NO₃, C₂F₃Cl was used to scavenge it. The overall rate constant, k₁, for the reaction, HO₂ + NO₃ → products, was measured to be k₁=(3.0 ± 0.7)×10^?12 cm³ molecule^?1 s^?1 at (297 ± 2) K and P=(1.4 – 1.9) torr. This result is in reasonable agreement with the previous studies. Direct detection of HO₂ and OH radicals and the use of three sources of NO₃ enabled us to confirm the existence of the channel producing OH:HO₂+NO₃→OH+NO₂+O₂ (1a); the other possible channel is HO₂+NO₃→HNO₃+O₂ (1b). From our measurements and the computer simulations, the branching ratio, k_1a/(k_1a + k_1b), was estimated to be (1.0). The rate coefficient for the reaction of OH with NO₃ was determined to be (2.1 ± 1.0) × 10^?11 cm³ molecule^?1 s^?1. © 1993 John Wiley & Sons, Inc.     

Semiclassical transition state theory/master equation kinetics of HO + CO: Performance evaluation

Previously, master equation (ME) simulations using semiclassical transition state theory (SCTST) and high-accuracy extrapolated ab initio thermochemistry (HEAT) predicted rate constants in excellent agreement with published experimental data over a wide range of pressure and temperatures ≳250 K, but the agreement was not as good at lower temperatures. Possible reasons for this reduced performance are investigated by (a) critically evaluating the published experimental data and by investigating; (b) three distinct ME treatments of angular momentum, including one that is exact at the zero- and infinite-pressure limits; (c) a hindered-rotor model for HOCO that implicitly includes the cis- and trans-conformers; (d) possible empirical adjustments of the thermochemistry; (e) possible empirical adjustments to an imaginary frequency controlling tunneling; (f) including or neglecting the prereaction complex PRC1; and (g) its possible bimolecular reactions. Improvements include better approximations to factors in SCTST and using the Hill and van Vleck treatment of angular momentum coupling. Evaluation of literature data does not reveal any specific shortcomings, but the stated uncertainties may be underestimated. All ME treatments give excellent fits to experimental data at T ≥ 250 K, but the discrepancy at T < 250 K persists. Note that each ME model requires individual empirical energy transfer parameters. Thermochemical adjustments were unable to match the experimental H/D kinetic isotope effects. Adjusting an imaginary frequency can achieve good fits, but the adjustments are unacceptably large. Whether PRC1 and its possible bimolecular reactions are included had little effect. We conclude that none of the adjustments is an improvement over the unadjusted theory. Note that only one set of experimental data exists in the regime of the discrepancy with theory, and data for DO + CO are scanty.     

Kinetics and mechanism of the glyoxal + HO2 reaction: conversion of HO2 to OH by carbonyls

The kinetics of the glyoxal + HO(2) reaction have been investigated using computational chemistry and statistical reaction rate theory techniques, with consideration of a novel pathway that results in the conversion of HO(2) to OH. Glyoxal is shown to react with HO(2) to form an α-hydroxyperoxy radical with additional α-carbonyl functionality. Intramolecular H atom abstraction from the carbonyl moiety proceeds with a relatively low barrier, facilitating decomposition to OH + CO + HC(O)OH (formic acid). Time-dependent master equation simulations demonstrate that direct reaction to form OH is relatively slow at ambient temperature. The major reaction product is predicted to be collisionally deactivated HC(OH)(OO)CHO, which predominantly dissociates to reform the reactants under low-NO(x) conditions. The mechanism described here for the conversion of OH to HO(2) is available to a diverse range of carbonyls, including methylglyoxal, glycolaldehyde, hydroxyacetone, and glyoxylic acid, and energy surfaces are reported for the reaction of these species with HO(2).     

Kinetics of the reaction HO2 + NO2(+M) = HO2NO2 using molecular modulation spectrometry

Rate constants for the reaction HO₂ + NO₂(+ M) = HO₂NO₂(+ M) have been obtained from direct observations of the HO₂ radical using the technique of molecular modulation ultraviolet spectrometry. HO₂ was generated by periodic photolysis of Cl₂ in the presence of excess H₂ and O₂, and k₁ was determined from the measured concentrations and lifetime of HO₂ with NO₂ present. k₁ increased with pressure in the range of 40–600 Torr, and a simple energy transfer model gave the following limiting second- and third-order rate constants at 283 K: k₁^∞ = 1.5 ± 0.5 × 10^?12 cm³/molec·sec and k₁^III = 2.5 ± 0.5 × 10^?31 cm⁶/molec·sec. The ultraviolet absorption spectrum of peroxynitric acid was also recorded in the range of 195–265 nm; it showed a broad feature with a maximum at 200 nm, σ_max = 4.4 × 10^?18 cm².     

Direct determination of the rate constant of the reaction NO + HO2 → NO2 + OH with the LMR

The rate of the reaction was determined in an isothermal discharge flow reactor with a combined ESR–LMR detection under pseudo-first-order conditions in HO₂. The rate constant was identical in experiments with two different HO₂ sources: F + H₂O₂ and H + O₂ + M. The absolute rate constant at T = 293 K was measured as In the range 2 ≤ p mbar ≤ 17 no pressure dependence for k₁ was found.     

Reaction kinetics of CO + HO(2) --> products: ab initio transition state theory study with master equation modeling

The kinetics of the reaction CO + HO2* --> CO2 + *OH was studied using a combination of ab initio electronic structure theory, transition state theory, and master equation modeling. The potential energy surface was examined with the CCSD(T) and CASPT2 methods. The classical energy barriers were found to be about 18 and 19 kcal/mol for CO + HO2* addition following the trans and cis paths, respectively. For the cis path, rate constant calculations were carried out with canonical transition state theory. For the trans path, master equation modeling was also employed to examine the pressure dependence. Special attention was paid to the hindered internal rotations of the HOOC*O adduct and transition states. The theoretical analysis shows that the overall rate coefficient is independent of pressure up to 500 atm for temperature ranging from 300 to 2500 K. On the basis of this analysis, we recommend the following rate expression for reaction R1 k(cm(3)/mol x s) = 1.57 x 10(5) T(2.18)e(-9030/T) for 300 < or = T < or = 2500 K with the uncertainty factor equal to 8, 2, and 1.7 at temperatures of 300, 1000, and 2000 K, respectively.     

Experimental and master equation study of the kinetics of OH + C2H2: temperature dependence of the limiting high pressure and pressure dependent rate coefficients

The kinetics of the reaction OH + C2H2 have been studied using laser flash photolysis at 248 nm to generate OH radicals and laser-induced fluorescence to monitor OH removal. An attempt was made to use the rate coefficients OH (v = 1,2) + C2H2 to determine the limiting high-pressure rate coefficient, k(1a)(infinity), over the temperature range of 195-823 K. This method is usually applicable if the reaction samples the potential energy well of the adduct, HOC2H2, and if intramolecular vibrational relaxation is fast. In the present case, however, the rate coefficients for loss of the vibrationally excited states by reaction with C2H2 also contain a substantial contribution from nonreactive vibrational relaxation, which occurs via a mechanism that does not sample the adduct potential energy well but involves, at least at low temperatures, collisions that access a shallower, longer range van der Waals well. The data were analyzed using a composite mechanism that incorporates both reactive and nonreactive energy transfer mechanisms, which allows the determination of k(1a)(infinity)(T) for OH + C2H2 with satisfactory accuracy over the temperature range 195-823 K. The kinetics of the reaction OH (v = 0) + C2H2 were also studied in He over the range of conditions: 210-373 K and 5-760 Torr. A one-dimensional master equation (ME) analysis of the experimental data provided a further determination of k(1a)(infinity)(T) and also (down) for He. Combining the two sets of results gives a consistent dataset for k(1a)(infinity) and the Arrhenius parameters A1ainfinity = 7.3 x 10(-12) cm(3) molecule(-1) s(-1) and E(1a)(infinity) = 5.3 kJ mol(-1), with (down) = 150(T/300 K) cm(-1). Additional experiments were conducted at room temperature in N(2) and SF(6) by laser flash photolysis with cavity ring down spectroscopy, and ME calculations were then optimized for the pressure falloff in N(2) by varying the average downward energy transfer parameter ((down)). The output from the best fit ME was parametrized using a modified Troe expression to provide rate data for use in atmospheric modeling.     

Methyl vinyl ketone+OH and methacrolein+OH oxidation reactions: a master equation analysis of the pressure- and temperature-dependent rate constants

High-level electronic structure calculations and master equation analyses were carried out to obtain the pressure- and temperature-dependent rate constants of the methyl vinyl ketone+OH and methacrolein+OH reactions. The balance between the OH addition reactions at the high-pressure limit, the OH addition reactions in the fall-off region, and the pressure-independent hydrogen abstractions involved in these multiwell and multichannel systems, has been shown to be crucial to understand the pressure and temperature dependence of each global reaction. In particular, the fall-off region of the OH addition reactions contributes to the inverse temperature dependence of the rate constants in the Arrhenius plots, leading to pressure-dependent negative activation energies. The pressure dependence of the methyl vinyl ketone+OH reaction is clearly more important than in the case of the methacrolein+OH reaction owing to the weight of the hydrogen abstraction process in this second system. Comparison of the theoretical rate constants and the experimental measurements shows quite good agreement.     

The rate constant of the reaction HO2 + COCO2 + OH

The high-temperature oxidation of formaldehyde in the presence of carbon monoxide was investigated to determine the rate constant of the reaction HO₂ + CO ? CO₂ + OH (10). In the temperature range of 878–952°K from the initial parts of the kinetic curves of the HO₂ radicals and CO₂ accumulation at small extents of the reaction, when the quantity of the reacted formaldehyde does not exceed 10%, it was determined that the rate constant k₁₀ is A computer program was used to solve the system of differential equations which correspond to the high-temperature oxidation of formaldehyde in the presence of carbon monoxide. The computation confirmed the experimental results. Also discussed are existing experimental data related to the reaction of HO₂ with CO.     

HO2 formation from the OH + benzene reaction in the presence of O2

In this study we investigated the secondary formation of HO(2) following the benzene + OH reaction in N(2) with variable O(2) content at atmospheric pressure and room temperature in the absence of NO. After pulsed formation of OH, HO(x) (= OH + HO(2)) and OH decay curves were measured by means of a laser-induced fluorescence technique (LIF). In synthetic air the total HO(2) yield was determined to be 0.69 ± 0.10 by comparison to results obtained with CO as a reference compound. HO(2) is expected to be a direct product of the reaction of the intermediately formed OH-benzene adduct with O(2). The HO(2) yield is slightly greater than the currently recommended yield of the proposed HO(2) co-product phenol (～53%). This hints towards other, minor HO(2) forming channels in the absence of NO, e.g. the formation of epoxide species that was proposed in the literature. For other test compounds upper limits of HO(2) yields of 0.10 (isoprene) and 0.05 (cyclohexane) were obtained, respectively. In further experiments at low O(2) concentrations (0.06-0.14% in N(2)) rate constants of (2.4 ± 1.1) × 10(-16) cm(3) s(-1) and (5.6 ± 1.1) × 10(-12) cm(3) s(-1) were estimated for the OH-benzene adduct reactions with O(2) and O(3), respectively. The rate constant of the unimolecular dissociation of the adduct back to benzene + OH was determined to be (3.9 ± 1.3) s(-1). The HO(2) yield at low O(2) was similar to that found in synthetic air, independent of O(2) and O(3) concentrations indicating comparable HO(2) yields for the adduct + O(2) and adduct + O(3) reactions.     

Temperature dependence of the HO2 + ClO reaction. 2. Reaction kinetics using the discharge-flow resonance-fluorescence technique

The total rate coefficient, k3, for the reaction HO2 + ClO --> products has been determined over the temperature range of 220-336 K at a total pressure of approximately 1.5 Torr of helium using the discharge-flow resonance-fluorescence technique. Pseudo-first-order conditions were used with both ClO and HO2 as excess reagents using four different combinations of precursor molecules. HO2 molecules were formed by using either the termolecular association of H atoms in an excess of O2 or via the reaction of F atoms with an excess of H(2)O(2). ClO molecules were formed by using the reaction of Cl atoms with an excess of O3 or via the reaction of Cl atoms with Cl(2)O. Neither HO2 nor ClO were directly observed during the course of the experiments, but these species were converted to OH or Cl radicals, respectively, via reaction with NO prior to their observation. OH fluorescence was observed at 308 nm, whereas Cl fluorescence was observed at approximately 138 nm. Numerical simulations show that under the experimental conditions used secondary reactions did not interfere with the measurements; however, some HO2 was lost on conversion to OH for experiments in excess HO2. These results were corrected to compensate for the simulated loss. At 296 K, the rate coefficient was determined to be (6.4 +/- 1.6) x 10(-12) cm3 molecule(-1) s(-1). The temperature dependence expressed in Arrhenius form is (1.75 +/- 0.52) x 10-12 exp[(368 +/- 78)/T] cm3 molecule(-1) s(-1). The Arrhenius expression is derived from a fit weighted by the reciprocal of the measurement errors of the individual data points. The uncertainties are cited at the level of two standard deviations and contain contributions from statistical errors from the data analysis in addition to estimates of the systematic experimental errors and possible errors from the applied model correction.     

A reinvestigation of the kinetics and the mechanism of the CH3C(O)O2 + HO2 reaction using both experimental and theoretical approaches

The kinetics and the mechanism of the reaction CH(3)C(O)O(2)+ HO(2) were reinvestigated at room temperature using two complementary approaches: one experimental, using flash photolysis/UV absorption technique and one theoretical, with quantum chemistry calculations performed using the density functional theory (DFT) method with the three-parameter hybrid functional B3LYP associated with the 6-31G(d,p) basis set. According to a recent paper reported by Hasson et al., [J. Phys. Chem., 2004, 108, 5979-5989] this reaction may proceed by three different channels: CH(3)C(O)O(2)+ HO(2)--> CH(3)C(O)OOH + O(2) (1a); CH(3)C(O)O(2)+ HO(2)--> CH(3)C(O)OH + O(3) (1b); CH(3)C(O)O(2)+ HO(2)--> CH(3)C(O)O + OH + O(2) (1c). In experiments, CH(3)C(O)O(2) and HO(2) radicals were generated using Cl-initiated oxidation of acetaldehyde and methanol, respectively, in the presence of oxygen. The addition of amounts of benzene in the system, forming hydroxycyclohexadienyl radicals in the presence of OH, allowed us to answer that channel (1c) is <10%. The rate constant k(1) of reaction (1) has been finally measured at (1.50 +/- 0.08) x 10(-11) cm(3) molecule(-1) s(-1) at 298 K, after having considered the combination of all the possible values for the branching ratios k(1a)/k(1,)k(1b)/k(1,)k(1c)/k(1) and has been compared to previous measurements. The branching ratio k(1b)/k(1), determined by measuring ozone in situ, was found to be equal to (20 +/- 1)%, a value consistent with the previous values reported in the literature. DFT calculations show that channel (1c) is also of minor importance: it was deduced unambiguously that the formation of CH(3)C(O)OOH + O(2) (X (3)Sigma(-)(g)) is the dominant product channel, followed by the second channel (1b) leading to CH(3)C(O)OH and singlet O(3) and, much less importantly, channel (1c) which corresponds to OH formation. These conclusions give a reliable explanation of the experimental observations of this work. In conclusion, the present study demonstrates that the CH(3)C(O)O(2)+ HO(2) is still predominantly a radical chain termination reaction in the tropospheric ozone chain formation processes.     

Quantum chemical and master equation study of OH + CH2O → H2O + CHO reaction rates in supercritical CO2 environment

We investigated the reaction rates of OH + CH₂O → H₂O + CHO at CO₂ pressures of up to 1000 atm with and without CO₂ molecule included in a reactive complex. Both mechanisms begin with formation of the hydrogen-bonded prereactive complexes. Our ab initio calculations indicate a possibility of catalytic effect, predicting an activation barrier that one order of magnitude lower when the CO₂ molecule is involved. To verify this effect, we use the Rice–Ramsperger–Kassel–Marcus theory and solve unimolecular master equations in the steady-state approximation. We assume the equilibrium between prereactive complexes and reactants and compare the bimolecular reaction rates for the two mechanisms. The catalyzed reaction mechanism is found to be faster at higher CO₂ pressures and lower temperatures, when prereactive complexes have nonnegligible concentration. Therefore, this catalytic effect may be important for this reactive process in room temperature supercritical CO₂ solvent, but is unlikely to play a role during oxy-combustion.     

Temperature dependence of the reactions of HO2 with NO and NO2

Mixtures of N₂O, H₂, O₂, and trace amounts of NO and NO₂ were photolyzed at 213.9 nm, at 245°–328°K, and at about 1 atm total pressure (mostly H₂). HO₂ radicals are produced from the photolysis and they react as follows: Reaction (1b) is unimportant under all of our reaction conditions. Reaction (1a) was studied in competition with reaction (3) from which it was found that k_1a/k₃^1/2 = 6.4 × 10^?6 exp { z?(1400 ± 500)/RT} cm^3/2/sec^1/2. If k₃ is taken to be 3.3 × 10^?12 cm³/sec independent of temperature, k_1a = 1.2 × 10^?11 exp {?(1400 ± 500)/RT} cm³/sec. Reaction (2a) is negligible compared to reaction (2b) under all of our reaction conditions. The ratio k_2b/k₁ = 0.61 ± 0.15 at 245°K. Using the Arrhenius expression for k_1a given above leads to k_2b = 4.2 × 10^?13 cm³/sec, which is assumed to be independent of temperature. The intermediate HO₂NO₂ is unstable and induces the dark oxidation of NO through reaction (?2b), which was found to have a rate coefficient k_?2b = 6 × 10¹⁷ exp {?26,000/RT} sec^?1 based on the value of k_1a given above. The intermediate can also decompose via Reaction (10b) is at least partially heterogeneous.     

Determination of the temperature and pressure dependence of the reaction OH + C2H4 from 200-400 K using experimental and master equation analyses

The pressure and temperature dependence for the reaction of OH + C(2)H(4) was studied over the range of conditions: 200-400 K and 5-600 Torr by laser flash photolysis, laser-induced fluorescence (FP-LIF). Additional experiments were conducted at room temperature by laser flash photolysis, cavity ring-down spectroscopy to facilitate determination of the high pressure limit. One-dimensional master equation calculations were conducted to test the temperature and pressure dependence of the reaction in He and in N(2). The energetics of the reaction and geometries of intermediate species were calculated by ab initio calculations (DFT-BH&HLYP/6-311+G(3df,2p) and CBS-APNO level along DFT-IRC, respectively. An investigation into the importance of a pre-reaction van der Waals complex on the kinetics over the pressure range of the troposphere is discussed. The high pressure rate coefficient was extracted by fitting the master equation calculations to the data and yields k(infinity) = 5.01 x 10(-12) exp(148/T) cm(3) molecule(-1) s(-1). The master equation calculations were then optimized for the pressure fall-off in He and N(2) by varying the average downward energy transfer parameter (DeltaE(down)) for the different collision partners and finally fitted to a Troe expression to determine k(o) and F(cent) for use in atmospheric modeling.     

A quantum mechanical approach to the kinetics of the hydrogen abstraction reaction H2O2 + •OH → HO2 + H2O

The kinetics of the hydrogen abstraction from H₂O₂ by ^?OH has been modeled with MP2/6‐31G*//MP2/6‐31G*, MP2‐SAC//MP2/6‐31G*, MP2/6‐31+G**//MP2/6‐31+G**, MP2‐SAC// MP2/6‐31+G**, MP4(SDTQ)/6‐311G**//MP2/6‐31G*, CCSD(T)/6‐31G*//CCSD(T)/6‐31G*, CCSD(T)/6‐31G**//CCSD(T)/6‐31G**, CCSD(T)/6‐311++G**//MP2/6‐31G* in the gas phase. MD simulations have been used to generate initial geometries for the stationary points along the potential energy surface for hydrogen abstraction from H₂O₂. The effective fragment potential (EFP) has been used to optimize the relevant structures in solution. Furthermore, the IEFPCM model has been used for the supermolecules generated via MD calculations. IEFPCM/MP2/6‐31G* and IEFPCM/CCSD(T)/6‐31G* calculations have also been performed for structures without explicit water molecules. Experimentally, the rate constant for hydrogen abstraction by ^?OH drops from 1.75 × 10^?12 cm³ molecule^?1 s^?1 in the gas phase to 4.48 × 10^?14 cm³ molecule^?1 s^?1 in solution. The same trend has been reproduced best with MP4 (SDTQ)/6‐311G**//MP2/6‐31G* in the gas phase (0.415 × 10^?12 cm³ molecule^?1 s^?1) and with EFP (UHF/6‐31G*) in solution (3.23 × 10^?14 cm³ molecule^?1 s^?1). © 2005 Wiley Periodicals, Inc. Int J Chem Kinet 37: 502–514, 2005     

Reactions of the HO2 radical with OH,H, Fe2+ and Cu2+ at elevated temperatures

The temperature dependence of the rate constant for the reactions of HO₂ with OH, H, Fe²⁺ and Cu²⁺ has been determined using pulse radiolysis technique. The following rate constants, k (dm³ mol⁻¹ s⁻¹) at 20°C and activation energies, E_a (kJ mol⁻¹) have been found. The reaction with OH was studied in the temperature range 20–296°C (k=7.0×10⁹, E_a=7.4) and the reaction with H in the temperature range 5–149°C (k=8.5×10⁹, E_a=17.5). The reaction with Fe²⁺ was studied in the temperature range 16–118°C (k=7.9×10⁵, E_a=36.8) and the reaction with Cu²⁺ in the temperature range 17–211°C (k=1.1×10⁸, E_a=14.9).     

Dynamics and kinetics of the S + HO2 reaction: A theoretical study

We report a quasi‐classical trajectory study of the S + HO₂ reaction using a previously reported global potential energy surface for the ground electronic state of HSO₂. Zero‐point energy leakage is approximately accounted for by using the vibrational energy quantum mechanical threshold method. Calculations are carried out both for specific ro‐vibrational states of the reactants and thermalized ones, with rate constants being reported as a function of temperature. The results suggest that the title reaction is capture type, with OH and SO showing as the most favorable products. The internal energy distribution of such products and the reaction mechanism are also investigated. © 2008 Wiley Periodicals, Inc. Int J Chem Kinet 40: 533–540, 2008     

Rate constant for the reaction Cl + HO2NO2 → products

The total rate constant for the reaction of Cl atoms with HO₂NO₂ was found to be less than 1.0 × 10^?13 cm³ s^?1 at 296 K by the discharge flow/resonance fluorescence technique. The reaction was also studied by the discharge flow/mass spectrometric technique. k_1a + k_1b was measured to be (3.4 ± 1.4) × 10^?14 cm³ s^?1 at 296 K. The reaction is too slow to be of any importance in stratospheric chemistry.