Kinetic study of the reactions of OH and OD with HBr and DBr

Convergent ab initio analysis of the multi-channel HOBr + H reaction

High-level potential energy surfaces for three reactions of hypobromous acid with atomic hydrogen were computed at the CCSDTQ/CBS//CCSDT(Q)/complete basis set level of theory. Focal point analysis was utilized to extrapolate energies and gradients for energetics and optimizations, respectively. The H attack at Br and subsequent Br-O cleavage were found to proceed barrierlessly. The slightly submerged transition state lies -0.2 kcal mol-1 lower in energy than the reactants and produces OH and HBr. The two other studied reaction paths are the radical substitution to produce H2O and Br with a 4.0 kcal mol-1 barrier and the abstraction at hydrogen to produce BrO and H2 with an 11.2 kcal mol-1 barrier. The final product energies lie -37.2, -67.9, and -7.3 kcal mol-1 lower in energy than reactants, HOBr + H, for the sets of products OH + HBr, H2O + Br, and H2 + BrO, respectively. Additive corrections computed for the final energetics, particularly the zero-point vibrational energies and spin-orbit corrections, significantly impacted the final stationary point energies, with corrections up to 6.2 kcal mol-1.

Experimental and Theoretical Study of the Kinetics of the CH3 + HBr → CH4 + Br Reaction and the Temperature Dependence of the Activation Energy of CH4 + Br → CH3 + HBr

The rate coefficient of the reaction of CH3 with HBr was measured and calculated in the temperature range 225-960 K. The results of the measurements performed in a flow apparatus with mass spectrometric detection agree very well with the quasiclassical trajectory calculations performed on a previously developed potential energy surface. The experimental rate coefficients are described well with a double-exponential fit, k1(exp) = [1.44 × 10-12 exp(219/T) + 6.18 × 10-11 exp(-3730/T)] cm3 molecule-1 s-1. The individual rate coefficients below 500 K accord with the available experimental data as does the slightly negative activation energy in this temperature range, -1.82 kJ/mol. At higher temperatures, the activation energy was found to switch sign and it rises up to about an order of magnitude larger positive value than that below 500 K, and the rate coefficient is about 50% larger at 960 K than that around room temperature. The rate coefficients calculated with the quasiclassical trajectory method display the same tendencies and are within about 8% of the experimental data between 960 and 300 K and within 25% below that temperature. The significant variation of the magnitude of the activation energy can be reconciled with the tabulated heats of formation only if the activation energy of the reverse CH4 + Br reaction also significantly increases with the temperature.

The reaction between the bromine atom and the water trimer: high level theoretical studies

The Br + (H2O)3 → HBr + (H2O)2OH reaction has been investigated using the CCSD(T) method with the basis sets as large as cc-pVQZ(-PP). The Br + (H2O)3 reaction is also compared with related Br + H2O/(H2O)2 and F/Cl + (H2O)3 reactions.

Using quantum dynamics to study the effect of energy efficiency on the reactivity of the OH + DBr reaction

We report a time-dependent, full dimensional, wave-packet calculation for the reaction of OH + DBr to examine the effect of the energy efficiency on the reactivity. This study shows that the vibrational excitations of the OH and DBr enhance the reaction. However, the rotational excitations of OH and DBr both hinder the reaction. As a result, the vibrational energies of both the OH and DBr reactants are more efficient at promoting the reactivity than the translational energy, while the rotational energies of OH and DBr are less effective than the translational energy. By analyzing the state population of the vibrational and rotational states along the reaction pathway, we also developed an approach in order to explain the enhancement of the vibrational excitation and the hindrance of the rotational excitation of the reaction. We found that the initial-state selected vibrational excited states of OH and DBr are the dominant components, respectively, for surmounting the barrier. However, the initial-state selected rotational excited states of OH and DBr are no longer the dominant states for surmounting the transition state owing to their population changes in the van der Waals well. This quantitative analysis demonstrates the potential well in the entrance valley plays an important role in the energy efficiency with regards to the reactivity.

Rate Constant of the Reaction of OH Radicals with HBr over the Temperature Range 235-960 K

The kinetics of the reaction of hydroxyl radicals with HBr, important in atmospheric and combustion chemistry, has been studied in a discharge flow reactor combined with an electron impact ionization quadrupole mass spectrometer in the temperature range 235-960 K. The rate constant of the reaction OH + HBr → H₂O + Br (1) was determined using both a relative rate method (using the reaction of OH with Br₂ as a reference) and absolute measurements, monitoring the kinetics of OH consumption under pseudo-first-order conditions in excess of HBr. The observed U-shaped temperature dependence of k₁ is well represented by the sum of two exponential functions: k₁ = 2.53 × 10^-11 exp(-364/T) + 2.79 × 10^-13 exp(784/T) cm³ molecule^-1 s^-1 (with an estimated conservative uncertainty of 15% at all temperatures). This expression for k₁, recommended for T = 240-960 K, combined with that from previous low temperature studies, k₁ = 1.06 × 10^-11 (T/298)^-0.9 cm³ molecule^-1 s^-1 at T = 23-240 K, allows to describe the temperature behavior of the rate constant over an extended temperature range 23-960 K. The current direct measurements of k₁ at temperatures above 460 K, the only ones to date, provide an experimental dataset for use in combustion and volcanic plume modeling and an experimental basis to test theoretical calculations.

Quantum dynamics study of kinetic isotope effects of OD with HBr and DBr

Comparison of kinetic isotope effects between quantum dynamics calculations and experiments shows that they agree well with each other both qualitatively and quantitatively.

Kinetics and Products of the Reaction of OH Radicals with ClNO from 220 to 940 K

The kinetics and products of the reaction of OH radicals with ClNO have been studied in a flow reactor coupled with an electron impact ionization mass spectrometer at nearly 2 Torr total pressure of helium and over a wide temperature range, T = 220-940 K. The rate constant of the reaction OH + ClNO → products was determined under pseudo-first order conditions, monitoring the kinetics of OH consumption in excess of ClNO: k₁ = 1.48 × 10^-18 × T^2.12 exp(146/T) cm³ molecule^-1 s^-1 (uncertainty of 15%). HOCl, Cl, and HONO were observed as the reaction products. As a result of quantitative detection of HOCl and Cl, the partial rate constants of the HOCl + NO and Cl + HONO forming reaction pathways were determined in the temperature range 220-940 K: k_1a = 3.64 × 10^-18 × T^1.99 exp(-114/T) and k_1b = 4.71 × 10^-18 × T^1.74 exp(246/T) cm³ molecule^-1 s^-1 (uncertainty of 20%). The dynamics of the title reaction and, in particular, non-Arrhenius behavior observed for both k_1a and k_1b in a wide temperature range, seems to be an interesting topic for theoretical research.

Quantum dynamics study of energy requirement on reactivity for the HBr + OH reaction with a negative-energy barrier

A time-dependent, quantum reaction dynamics approach in full dimensional, six degrees of freedom was carried out to study the energy requirement on reactivity for the HBr + OH reaction with an early, negative energy barrier. The calculation shows both the HBr and OH vibrational excitations enhance the reactivity. However, even this reaction has a negative energy barrier, the calculation shows not all forms of energy are equally effective in promoting the reactivity. On the basis of equal amount of total energy, the vibrational energies of both the HBr and OH are more effective in enhancing the reactivity than the translational energy, whereas the rotational excitations of both the HBr and OH hinder the reactivity. The rate constants were also calculated for the temperature range between 5 to 500 K. The quantal rate constants have a better slope agreement with the experimental data than quasi-classical trajectory results.

Stereodirectional Origin of anti-Arrhenius Kinetics for a Tetraatomic Hydrogen Exchange Reaction: Born-Oppenheimer Molecular Dynamics for OH + HBr

Among four-atom processes, the reaction OH + HBr → H2O + Br is one of the most studied experimentally: its kinetics has manifested an unusual anti-Arrhenius behavior, namely, a marked decrease of the rate constant as the temperature increases, which has intrigued theoreticians for a long time. Recently, salient features of the potential energy surface have been characterized and most kinetic aspects can be considered as satisfactorily reproduced by classical trajectory simulations. Motivation of the work reported in this paper is the investigation of the stereodirectional dynamics of this reaction as the prominent reason for the peculiar kinetics: we started in a previous Letter ( J. Phys. Chem. Lett. 2015 , 6 , 1553 - 1558 ) a first-principles Born-Oppenheimer "canonical" molecular dynamics approach. Trajectories are step-by-step generated on a potential energy surface quantum mechanically calculated on-the-fly and are thermostatically equilibrated to correspond to a specific temperature. Here, refinements of the method permitted a major increase of the number of trajectories and the consideration of four temperatures -50, +200, +350, and +500 K, for which the sampling of initial conditions allowed us to characterize the stereodynamical effect. The role is documented of the adjustment of the reactants' mutual orientation to encounter the entrance into the "cone of acceptance" for reactivity. The aperture angle of this cone is dictated by a range of directions of approach compatible with the formation of the specific HOH angle of the product water molecule; and consistently the adjustment is progressively less effective the higher the kinetic energy. Qualitatively, this emerging picture corroborates experiments on this reaction, involving collisions of aligned and oriented molecular beams, and covering a range of energies higher than the thermal ones. The extraction of thermal rate constants from this molecular dynamics approach is discussed and the systematic sampling of the canonical ensemble is indicated as needed for quantitative comparison with the kinetic experiments.

Stereodynamical Origin of Anti-Arrhenius Kinetics: Negative Activation Energy and Roaming for a Four-Atom Reaction

The OH + HBr → H2O + Br reaction, prototypical of halogen-atom liberating processes relevant to mechanisms for atmospheric ozone destruction, attracted frequent attention of experimental chemical kinetics: the nature of the unusual reactivity drop from low to high temperatures eluded a variety of theoretical efforts, ranking this one among the most studied four-atom reactions. Here, inspired by oriented molecular-beams experiments, we develop a first-principles stereodynamical approach. Thermalized sets of trajectories, evolving on a multidimensional potential energy surface quantum mechanically generated on-the-fly, provide a map of most visited regions at each temperature. Visualizations of rearrangements of bonds along trajectories and of the role of specific angles of reactants' mutual approach elucidate the mechanistic change from the low kinetic energy regime (where incident reactants reorient to find the propitious alignment leading to reaction) to high temperature (where speed hinders adjustment of directionality and roaming delays reactivity).

Energy disposal and thermal rate constants for the OH + HBr and OH + DBr reactions: quasiclassical trajectory calculations on an accurate potential energy surface

We report reaction cross sections, energy disposal, and rate constants for the OH + HBr → Br + H2O and OH + DBr → Br + HDO reactions from quasiclassical trajectory calculations using an ab initio potential energy surface [ de Oliveira-Filho , A. G. S. ; Ornellas , F. R. ; Bowman , J. M. J. Phys. Chem. Lett. 2014 , 5 , 706 - 712 ]. Comparison with available experiments are made and generally show good agreement.

Reaction of limonene with F2: rate coefficient and products

The kinetics of the reaction of limonene (C10H16) with F2 has been studied using a low pressure (P = 1 Torr) and a high pressure turbulent (P = 100 Torr) flow reactor coupled with an electron impact ionization and chemical ionization mass spectrometers, respectively: F2 + Limonene → products (1). The rate constant of the title reaction was determined under pseudo-first-order conditions by monitoring either limonene or F2 decay in excess of F2 or C10H16, respectively. The reaction rate constant, k1 = (1.15 ± 0.25) × 10(-12) exp(160 ± 70)/T) was determined over the temperature range 278-360 K, independent of pressure between 1 (He) and 100 (N2) Torr. F atom and HF were found to be formed in reaction 1 , with the yields of 0.60 ± 0.13 and 0.39 ± 0.09, respectively, independent of temperature in the range 296-355 K.

Quasiclassical Trajectory Calculations of the Rate Constant of the OH + HBr → Br + H2O Reaction Using a Full-Dimensional Ab Initio Potential Energy Surface Over the Temperature Range 5 to 500 K

We report a permutationally invariant, ab initio potential energy surface (PES) for the OH + HBr → Br + H2O reaction. The PES is a fit to roughly 26 000 spin-free UCCSD(T)/cc-pVDZ-F12a energies and has no classical barrier to reaction. It is used in quasiclassical trajectory calculations with a focus on the thermal rate constant, k(T), over the temperature range 5 to 500 K. Comparisons with available experimental data over the temperature range 23 to 416 K are made using three approaches to treat the OH rotational and associated electronic partition function. All display an inverse temperature dependence of k(T) below roughly 160 K and a nearly constant temperature dependence above 160 K, in agreement with experiment. The calculated rate constant with no treatment of spin-orbit coupling is overall in the best agreement with experiment, being (probably fortuitously) within 20% of it.

Rotationally inelastic scattering of OH (Π3∕22, v=0, J=3∕2, f) by HBr (Σ1, v=0, J<4)

Relative state-to-state cross sections of OH molecules in the (2)Pi(32), v=0, J=32, M(J)=32, f state have been determined for transitions up to (2)Pi(32), v=0, J=112, f and (2)Pi(12), v=0, J=72, e states by collisions with HBr molecules ((1)Sigma, v=0, J<4) at 750 cm(-1) collision energy. In order to investigate features of the anisotropy of the OH-HBr potential energy surface, the steric asymmetries, which account for the effect of the OH orientation with respect to the collision partner, have been measured. A comparison with other systems previously studied shows strong similarities with the OH-HCl system.

Temperature Dependence and Kinetic Isotope Effects for the OH + HBr Reaction and H/D Isotopic Variants at Low Temperatures (53−135 K) Measured Using a Pulsed Supersonic Laval Nozzle Flow Reactor

The reactions of OH + HBr and all isotopic variants have been measured in a pulsed supersonic Laval nozzle flow reactor between 53 and 135 K, using a pulsed DC discharge to create the radical species and laser induced fluorescence on the A 2sigma <-- X 2pi (v' = 1 <-- v'' = 0) transition. All reactions are found to possess an inverse temperature dependence, in accord with previous work, and are fit to the form k = A(T/298)(-n), with k1 (OH + HBr) = (10.84 +/- 0.31) x 10(-12) (T/298)(-0.67+/-0.02) cm3/s, k2 (OD + HBr) = (6.43 +/- 2.60) x 10(-12) (T/298)(-1.19+/-0.26) cm3/s, k3 (OH + DBr) = (5.89 +/- 1.93) x 10(-12) (T/298)(-0.76+/-0.22) cm3/s, and k4 (OD + DBr) = (4.71 +/- 1.56) x 10(-12) (T/298)(-1.09+/-0.21) cm3/s. A global fit of k vs T over the temperature range 23-360 K, including the new OH + HBr data, yields kT = (1.06 +/- 0.02) x 10(-11) (T/298)(-0.90+/-0.11) cm3/s, and (0.96 +/- 0.02) x 10(-11) (T/298)(-0.90+/-0.03) exp((-2.88+/-1.82 K)/T) cm3/s, in accord with previous fits. In addition, the primary and secondary kinetic isotope effects are found to be independent of temperature within experimental error over the range investigated and take on the value of (kH/kD)(AVG) = 1.64 for the primary effect and (kH/kD)(AVG) = 0.87 for the secondary effect. These results are discussed within the context of current experimental and theoretical work.