Quantum Scattering Calculations of the H + O2 → O + OH Reaction

Impact of Geometric Phase on Dynamics of Complex-Forming Reactions: H + O2 → OH + O

Reaction dynamics on the ground electronic state might be significantly influenced by conical intersections (CIs) via the geometric phase (GP), as demonstrated for activated reactions (i.e., the H + H2 exchange reaction). However, there have been few investigations of GP effects in complex-forming reactions. Here, we report a full quantum dynamical study of an important reaction in combustion (H + O2 → OH + O), which serves as a proving ground for studying GP effects therein. The results reveal significant differences in reaction probabilities and differential cross sections (DCSs) obtained with and without GP, underscoring its strong impact. However, the GP effects are less pronounced for the reaction integral cross sections, apparently due to the integral of the DCS over the scattering angle. Further analysis indicated that the cross section has roughly the same contributions from the two topologically distinct paths around the CI, namely, the direct and looping paths.

3D time-dependent wave-packet approach in hyperspherical coordinates for the H + O2 reaction on the CHIPR and DMBE IV potential energy surfaces

3D wavepacket quantum dynamics methodology ICS calculation of H + O₂ reaction on the CHIPR and DMBE IV PESs by J-shifting scheme.

Quantum dynamics study on the CHIPR potential energy surface for the hydroperoxyl radical: the reactions O + OH⇋O2 + H

Quantum scattering calculations of the O((3)P)+OH((2)Π)⇌O2((3)Σg (-))+H((2)S) reactions are presented using the combined-hyperbolic-inverse-power-representation potential energy surface [A. J. C. Varandas, J. Chem. Phys. 138, 134117 (2013)], which employs a realistic, ab initio-based, description of both the valence and long-range interactions. The calculations have been performed with the ABC time-independent quantum reactive scattering computer program based on hyperspherical coordinates. The reactivity of both arrangements has been investigated, with particular attention paid to the effects of vibrational excitation. By using the J-shifting approximation, rate constants are also reported for both the title reactions.

State-to-state reaction probabilities for the H+O2(v,j)→O+OH(v′,j′) reaction on three potential energy surfaces

We report state-to-state and total reaction probabilities for J=0 and total reaction probabilities for J=2 and 4 for the title reaction, both for ground-state and initially rovibrationally excited reactants. The results for three different potential energy surfaces are compared and contrasted. The potential energy surfaces employed are the DMBE IV surface by Pastrana et al. [J. Phys. Chem. 94, 8073 (1990)], the surface by Troe and Ushakov (TU) [J. Chem. Phys. 115, 3621 (2001)], and the new XXZLG ab initio surface by Xu et al. [J. Chem. Phys. 122, 244305 (2005)]. Our results show that the total reaction probabilities from both the TU and XXZLG surfaces are much smaller in magnitude for collision energies above 1.2 eV compared to the DMBE IV surface. The three surfaces also show different behavior with regards to the effect of initial state excitation. The reactivity is increased on the XXZLG and the TU surfaces and decreased on the DMBE IV surface. Vibrational and rotational product state distributions for the XXZLG and the DMBE IV surface show different behaviors for both types of distributions. Our results show that for energies above 1.25 eV the dynamics on the DMBE IV surface are not statistical. However, there is also evidence that the dynamics on the XXZLG surface are not purely statistical for energies above the onset of the first excited product vibrational state v'=1. The magnitude of the total reaction probability is decreased for J>0 for the DMBE IV and the XXZLG surfaces for ground-state reactants. However, for initially rovibrationally excited reactants, the total reaction probability does not decrease as expected for both surfaces. As a result the total cross section averaged over all Boltzmann accessible rotational states may well be larger than the cross section reported in the literature for j=1.

Cross sections and rate constants for the C(P3)+OH(XΠ2)→CO(XΣ+1)+H(S2) reaction using a quasiclassical trajectory method

First quasiclassical trajectory calculations have been carried out for the C(3P)+OH(X 2Pi)-->CO(X 1Sigma+)+H(2S) reaction using a recent ab initio potential energy surface for the ground electronic state, X 2A', of HCO/COH. Total and state-specific integral cross sections have been determined for a wide range of collision energies (0.001-1 eV). Then, thermal and state-specific rate constants have been calculated in the 1-500 K temperature range. The thermal rate constant varies from 1.78x10(-10) cm3 s-1 at 1 K down to 5.96x10(-11) cm3 s-1 at 500 K with a maximum value of 3.39x10(-10) cm3 s-1 obtained at 7 K. Cross sections and rate constants are found to be almost independent of the rovibrational state of OH.

Global analytical potential energy surfaces for HO2(X2A") based on high-level ab initio calculations

Two global analytical potential energy surfaces for the HO2(X2A") system have been developed by fitting approximately 15,000 ab initio points at the icMRCI+Qaug-cc-pVQZ level of theory, using the reproducing kernel Hilbert space method. One analytical potential is designed to give a very accurate representation of the low energy range that determines the vibrational spectrum, while the other attempts to provide a fast and uniformly accurate potential function for reaction dynamics. The quality of the fitted potential functions is confirmed by good agreement of the (J=0) HO2 vibrational spectrum and (J=0) quantum reaction probability for the H+O2(ji=0,nui=0) reaction with those obtained using the spline fitted potential. Quasiclassical trajectory calculations carried out on the new potential energy surface provided the reaction probability with a zero impact parameter (b=0), which is in reasonably good agreement with the J=0 quantum results.

A new ab initio potential-energy surface of HO2(X2A") and quantum studies of HO2 vibrational spectrum and rate constants for the H + O2 <--> O + OH reactions

A new global potential-energy surface for the ground electronic state of HO(2)(X(2)A(")) has been developed by three-dimensional cubic spline interpolation of more than 15 000 ab initio points, which were calculated at the multireference configuration-interaction level with Davidson correction using the augmented correlation-consistent polarized valence quadruple zeta basis set. Low-lying vibrational states were obtained in this new potential using the Lanczos method and assigned. The calculated vibrational frequencies are in much better agreement with the available experimental band origins than those obtained from a previous potential. In addition, rate constants for the H+O(2) <--> O + OH reactions were obtained using a wave-packet-based statistical model. Reasonably good agreement with experimental data was obtained. These results demonstrate the accuracy of the potential.

Quantum Dynamics of the H + O2 → O + OH Reaction on an Accurate ab Initio Potential Energy Surface

We report exact time-dependent and time-independent quantum mechanical studies of the title reaction on an accurate ab initio potential energy surface of Xu et al. (J. Chem. Phys. 2005, 122, 24305). The J = 0 reaction probabilities for several reactant states show sharp resonance structures superimposed on relatively low backgrounds, and they are remarkably different from existing quantum results on an earlier potential energy surface (DMBE-IV). The new findings reported here suggest that our current understanding of this important reaction might require significant revision.

Theories of reactive scattering

This paper is an overview of the theory of reactive scattering, with emphasis on fully quantum mechanical theories that have been developed to describe simple chemical reactions, especially atom-diatom reactions. We also describe related quasiclassical trajectory applications, and in all of this review the emphasis is on methods and applications concerned with state-resolved reaction dynamics. The review first provides an overview of the development of the theory, including a discussion of computational methods based on coupled channel calculations, variational methods, and wave packet methods. Choices of coordinates, including the use of hyperspherical coordinates are discussed, as are basis set and discrete variational representations. The review also summarizes a number of applications that have been performed, especially the two most comprehensively studied systems, H+H2 and F+H2, along with brief discussions of a large number of other systems, including other hydrogen atom transfer reactions, insertion reactions, electronically nonadiabatic reactions, and reactions involving four or more atoms. For each reaction we describe the method used and important new physical insight extracted from the results.

Quantum mechanical investigations of the N(S4)+O2(XΣg−3)→NO(XΠ2)+O(P3) reaction

The reaction between energetic nitrogen atoms and oxygen molecules has received important attention in connection with nitric oxide chemistry in the lower thermosphere. We report time-independent quantum mechanical calculations of the N(4S)+O2-->NO+O reaction employing the X 2A' and a 4A' electronic potential energy surfaces of Sayos et al. [J. Chem. Phys. 117, 670 (2002)]. We confirm the production of highly vibrationally excited NO molecules, consistent with previous semiclassical and more recent time-dependent quantum wave packet studies. Calculations are carried out for total angular momentum quantum number J=0 and cross sections and rate coefficients are extracted using the J-shifting approximation. The results are in good agreement with available experimental and theoretical data.

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.