O+OH-->O(2)+H: A key reaction for interstellar chemistry. New theoretical results and comparison with experiment

Accurate Full-Dimensional Global Diabatic Potential Energy Matrix for the Two Lowest-Lying Electronic States of the H + O2 ↔ HO + O Reaction

A new and more accurate diabatic potential energy matrix (DPEM) is developed for the two lowest-lying electronic states of HO₂, covering both the strong interaction region and reaction asymptotes. The ab initio calculations were performed at the Davidson corrected multireference configuration interaction level with the augmented correlation-consistent polarized valence quintuple-zeta basis set (MRCI+Q/AV5Z). The accuracy of the electronic structure calculations is validated by excellent agreement with the experimental HO₂ equilibrium geometry, fundamental vibrational frequencies, and H + O₂ ↔ OH + O reaction energy. Through the combination of an electronic angular momentum-method and a configuration interaction vector-based method, the mixing angle between the first two ²A″ states of HO₂ was successfully determined. Elements of the 2×2 DPEM were fit to neural networks with a proper account of the complete nuclear permutation inversion symmetry of HO₂. The DPEM correctly predicted the properties of conical intersection seams at linear and T-shape geometries, thus providing a reliable platform for studying both the spectroscopy of HO₂ and the nonadiabatic dynamics for the H + O₂ ↔ OH + O reaction.

Dynamical calculations of O(3P) + OH(2Π) reaction on the CHIPR potential energy surface using the fully coupled time-dependent wave-packet approach in hyperspherical coordinates

We have carried out quantum dynamics calculations for the O + OH → H + O₂ reaction on the CHIPR [A. J. C. Varandas, J. Chem. Phys., 2013, 138, 134117] potential energy surface (PES) for ground state HO₂ using the fully coupled 3D time-dependent wavepacket formalism [S. Adhikari and A. J. C. Varandas, Comput. Phys. Commun., 2013, 184, 270] in hyperspherical coordinates. Reaction probabilities for J > 0 are calculated for different initial rotational states of the OH radical (v = 0; j = 0, 1). State-to-state as well as total integral cross sections and rate-coefficients are evaluated and compared with previous theoretical calculations and available experimental studies. Using the rate constant for the forward (hereinafter considered to be H + O₂ → O + OH) and backward (O + OH → H + O₂) reactions of this reactive system, the equilibrium constant of the reversible process [H + O₂ ⇌ O + OH] is calculated as a function of temperature and compared with previous experimental measurements.

Quasi-Classical Trajectory Study of the CN + NH3 Reaction Based on a Global Potential Energy Surface

Based on a combination of valence-bond and molecular mechanics functions which were fitted to high-level ab initio calculations, we constructed an analytical full-dimensional potential energy surface, named PES-2020, for the hydrogen abstraction title reaction for the first time. This surface is symmetrical with respect to the permutation of the three hydrogens in ammonia, it presents numerical gradients and it improves the description presented by previous theoretical studies. In order to analyze its quality and accuracy, stringent tests were performed, exhaustive kinetics and dynamics studies were carried out using quasi-classical trajectory calculations, and the results were compared with the available experimental evidence. Firstly, the properties (geometry, vibrational frequency and energy) of all stationary points were found to reasonably reproduce the ab initio information used as input; due to the complicated topology with deep wells in the entrance and exit channels and a “submerged” transition state, the description of the intermediate complexes was poorer, although it was adequate to reasonably simulate the kinetics and dynamics of the title reaction. Secondly, in the kinetics study, the rate constants simulated the experimental data in the wide temperature range of 25–700 K, improving the description presented by previous theoretical studies. In addition, while previous studies failed in the description of the kinetic isotope effects, our results reproduced the experimental information. Finally, in the dynamics study, we analyzed the role of the vibrational and rotational excitation of the CN(v,j) reactant and product angular scattering distribution. We found that vibrational excitation by one quantum slightly increased reactivity, thus reproducing the only experimental measurement, while rotational excitation strongly decreased reactivity. The scattering distribution presented a forward-backward shape, associated with the presence of deep wells along the reaction path. These last two findings await experimental confirmation.

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.

On the Angular Distribution of the H+Li2 Cross Sections: a Converged Time-Independent Quantum Scattering Study

A thorough time-independent quantum scattering study is performed on a benchmark potential energy surface for the H+Li₂ reaction at the fundamental electronic state. Integral and differential cross sections are calculated along with thermal rate coefficients until convergence is reached. Our findings show that vibrational and rotational excitations of the reactant hinder reactivity, though for the latter a considerable reaction promotion was spotted as we increase the reactant rotational quantum number until the critical value of j = 4. Such a promotion then begins to retract, eventually becoming an actual inhibition for larger j. In a straightforward manner, the concept of time-independent methods implemented in this study allowed this accurate state-to-state analysis. Furthermore, a nearly isotropic behaviour of the scattering is noted to take place from the angular point of view. Remarkably, our computational protocol is ideally suited to yield converged thermal rate coefficients, revealing a non-Arrhenius pattern and showing that J-shifting approach fails to describe this particular reaction. Our results, when compared to previous and independent ones, reinforce the latest theoretical reference for future validation in the experimental field.

Time-Dependent Quantum Wave Packet Study of the Si + OH → SiO + H Reaction: Cross Sections and Rate Constants

The dynamics of the Si(³P) + OH(X²Π) → SiO(X¹Σ⁺) + H(²S) reaction is investigated by means of the time-dependent wave packet (TDWP) approach using an ab initio potential energy surface recently developed by Dayou et al. ( J. Chem. Phys. 2013 , 139 , 204305 ) for the ground X²A' electronic state. Total reaction probabilities have been calculated for the first 15 rotational states j = 0-14 of OH(v=0,j) at a total angular momentum J = 0 up to a collision energy of 1 eV. Integral cross sections and state-selected rate constants for the temperature range 10-500 K were obtained within the J-shifting approximation. The reaction probabilities display highly oscillatory structures indicating the contribution of long-lived quasibound states supported by the deep SiOH/HSiO wells. The cross sections behave with collision energies as expected for a barrierless reaction and are slightly sensitive to the initial rotational excitation of OH. The thermal rate constants show a marked temperature dependence below 200 K with a maximum value around 15 K. The TDWP results globally agree with the results of earlier quasi-classical trajectory (QCT) calculations carried out by Rivero-Santamaria et al. ( Chem. Phys. Lett. 2014 , 610-611 , 335 - 340 ) with the same potential energy surface. In particular, the thermal rate constants display a similar temperature dependence, with TDWP values smaller than the QCT ones over the whole temperature range.

Low temperature rate coefficients of the H + CH(+) → C(+) + H2 reaction: New potential energy surface and time-independent quantum scattering

The observed abundances of the methylidyne cation, CH(+), in diffuse molecular clouds can be two orders of magnitude higher than the prediction of the standard gas-phase models which, in turn, predict rather well the abundances of neutral CH. It is therefore necessary to investigate all the possible formation and destruction processes of CH(+) in the interstellar medium with the most abundant species H, H2, and e(-). In this work, we address the destruction process of CH(+) by hydrogen abstraction. We report a new calculation of the low temperature rate coefficients for the abstraction reaction, using accurate time-independent quantum scattering and a new high-level ab initio global potential energy surface including a realistic model of the long-range interaction between the reactants H and CH(+). The calculated thermal rate coefficient is in good agreement with the experimental data in the range 50 K-800 K. However, at lower temperatures, the experimental rate coefficient takes exceedingly small values which are not reproduced by the calculated rate coefficient. Instead, the latter rate coefficient is close to the one given by the Langevin capture model, as expected for a reaction involving an ion and a neutral species. Several recent theoretical works have reported a seemingly good agreement with the experiment below 50 K, but an analysis of these works show that they are based on potential energy surfaces with incorrect long-range behavior. The experimental results were explained by a loss of reactivity of the lowest rotational states of the reactant; however, the quantum scattering calculations show the opposite, namely, a reactivity enhancement with rotational excitation.

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.

Electronic nonadiabatic effects in low temperature radical-radical reactions. I. C(3P) + OH(2Π)

The formation of collision complexes, as a first step towards reaction, in collisions between two open-electronic shell radicals is treated within an adiabatic channel approach. Adiabatic channel potentials are constructed on the basis of asymptotic electrostatic, induction, dispersion, and exchange interactions, accounting for spin-orbit coupling within the multitude of electronic states arising from the separated reactants. Suitable coupling schemes (such as rotational + electronic) are designed to secure maximum adiabaticity of the channels. The reaction between C((3)P) and OH((2)Π) is treated as a representative example. The results show that the low temperature association rate coefficients in general cannot be represented by results obtained with a single (generally the lowest) potential energy surface of the adduct, asymptotically reaching the lowest fine-structure states of the reactants, and a factor accounting for the thermal population of the latter states. Instead, the influence of non-Born-Oppenheimer couplings within the multitude of electronic states arising during the encounter markedly increases the capture rates. This effect extends up to temperatures of several hundred K.

A global ab initio potential energy surface for the X2A' ground state of the Si + OH → SiO + H reaction

We report the first global potential energy surface (PES) for the X(2)A' ground electronic state of the Si((3)P) + OH(X(2)Π) → SiO(X(1)Σg(+)) + H((2)S) reaction. The PES is based on a large number of ab initio energies obtained from multireference configuration interaction calculations plus Davidson correction (MRCI+Q) using basis sets of quadruple zeta quality. Corrections were applied to the ab initio energies in the reactant channel allowing a proper description of long-range interactions between Si((3)P) and OH(X(2)Π). An analytical representation of the global PES has been developed by means of the reproducing kernel Hilbert space method. The reaction is found barrierless. Two minima, corresponding to the SiOH and HSiO isomers, and six saddle points, among which the isomerization transition state, have been characterized on the PES. The vibrational spectra of the SiOH/HSiO radicals have been computed from second-order perturbation theory and quantum dynamics methods. The structural, energetic, and spectroscopic properties of the two isomers are in good agreement with experimental data and previous high quality calculations.

Chemical reaction versus vibrational quenching in low energy collisions of vibrationally excited OH with O

Quantum scattering calculations are reported for state-to-state vibrational relaxation and reactive scattering in O + OH(v = 2 - 3, j = 0) collisions on the electronically adiabatic ground state (2)A'' potential energy surface of the HO2 molecule. The time-independent Schrödinger equation in hyperspherical coordinates is solved to determine energy dependent probabilities and cross sections over collision energies ranging from ultracold to 0.35 eV and for total angular momentum quantum number J = 0. A J-shifting approximation is then used to compute initial state selected reactive rate coefficients in the temperature range T = 1 - 400 K. Results are found to be in reasonable agreement with available quasiclassical trajectory calculations. Results indicate that rate coefficients for O2 formation increase with increasing the OH vibrational level except at low and ultralow temperatures where OH(v = 0) exhibits a slightly different trend. It is found that vibrational relaxation of OH in v = 2 and v = 3 vibrational levels is dominated by a multi-quantum process.

State-to-state quantum dynamics of the N(4S) + CH(X2Π) → CN(X2Σ+,A2Π) + H(2S) reactions

The reactions between N((4)S) and CH(X(2)Π) lead to H((2)S) plus CN in its two lowest electronic states (X(2)Σ(+) and A(2)Π), which are responsible for the interstellar CN formation. Accurate quantum dynamics of these reactions are investigated on new global potential energy surfaces of the two lowest-lying triplet states of HCN (1(3)A' and 1(3)A") fitted to more than 37,000 points at the internally contracted multi-reference configuration interaction level with the Davidson correction. The pathways for these highly exothermic and barrierless reactions feature both the HCN and HNC wells. Long-lived resonances supported by these wells manifest in reaction probabilities as numerous oscillations, particularly for low J partial waves. The 1(3)A" state is found to be more reactive than the 1(3)A' state, due apparently to its more attractive nature in the entrance channel. The CN products in both electronic states are highly excited in both vibrational and rotational degrees of freedom. The near forward-backward symmetric differential cross sections are consistent with a complex-forming mechanism.

Global potential energy surface, vibrational spectrum, and reaction dynamics of the first excited (Ã (2)A(')) state of HO(2)

The authors report extensive high-level ab initio studies of the first excited (Ã (2)A(')) state of HO(2). A global potential energy surface (PES) was developed by spline-fitting 17 000 ab initio points at the internal contracted multireference configuration interaction (icMRCI) level with the AVQZ basis set. To ascertain the spectroscopic accuracy of the PES, the near-equilibrium region of the molecule was also investigated using three interpolating moving least-squares-based PESs employing dynamically weighted icMRCI methods in the complete basis set limit. Vibrational energy levels on all four surfaces agree well with each other and a new assignment of some vibrational features is proposed. In addition, the dynamics of both the forward and reverse directions of the H+O(2)(ã (1)Δ(g))↔OH+O reaction (J=0) were studied using an exact wave packet method. The reactions are found to be dominated by sharp resonances.