On the statistical behavior of the O+OH→H+O[sub 2] reaction: A comparison between quasiclassical trajectory, quantum scattering, and statistical calculations

Influence of Singlet Oxygen in H + O2 Collision Reaction

The collision reaction of H + O2 = OH + O is a pivotal step in combustion. To investigate the influence of singlet oxygen on this reaction, we computed potential energy surfaces (PESs) for all six lowest states using high-level ab initio methods and coupled them with embedded atom neural network (EANN) fitting. By integrating quasi-classical trajectory (QCT) with trajectory surface hopping (TSH) based on the fitted PESs, we simulated the dynamics of both ground- and excited-states to derive the reaction rate constants for the forward and reverse processes. The results reveal that the forward reaction facilitates radical generation, promoting combustion reactions. Furthermore, calculations of reverse reaction rate constants indicate that all electronic states ultimately yield ground-state oxygen, leading to radical deactivation and exerting an inhibitory effect on combustion processes.

Quantum and quasi-classical dynamics of the C(3P) + O2(3Σ-g) → CO(1Σ+) + O(1D) reaction on its electronic ground state

The dynamics of the C(³P) + O₂(³Σ⁻_g) → CO(¹Σ⁺) + O(¹D) reaction on its electronic ground state is investigated by using time-dependent wave packet propagation (TDWP) and quasi-classical trajectory (QCT) simulations. For the moderate collision energies considered (E_c = 0.001 to 0.4 eV, corresponding to a range from 10 K to 4600 K) the total reaction probabilities from the two different treatments of the nuclear dynamics agree very favourably. The undulations present in P(E) from the quantum mechanical treatment can be related to stabilization of the intermediate CO₂ complex with lifetimes on the 0.05 ps time scale. This is also confirmed from direct analysis of the TDWP simulations and QCT trajectories. Product diatom vibrational and rotational level resolved state-to-state reaction probabilities from TDWP and QCT simulations agree well except for the highest product vibrational states (v′ ≥ 15) and for the lowest product rotational states (j′ ≤ 10). Opening of the product vibrational level CO(v′ = 17) requires ∼0.2 eV from QCT and TDWP simulations with O₂(j = 0) and decreases to 0.04 eV if all initial rotational states are included in the QCT analysis, compared with E_c > 0.04 eV obtained from experiments. It is thus concluded that QCT simulations are suitable for investigating and realistically describe the C(³P) + O₂(³Σ⁻_g) → CO(¹Σ⁺) + O(¹D) reaction down to low collision energies when compared with results from a quantum mechanical treatment using TDWPs.

The dynamics of the C(³P) + O₂(³Σ⁻_g) → CO(¹Σ⁺) + O(¹D) reaction on its electronic ground state is investigated by using time-dependent wave packet propagation (TDWP) and quasi-classical trajectory (QCT) simulations.

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.

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.

Dynamics of the D(+) + H2 → HD + H(+) reaction at the low energy regime by means of a statistical quantum method

The D(+) +H2(v = 0, j = 0, 1) → HD+H(+) reaction has been investigated at the low energy regime by means of a statistical quantum mechanical (SQM) method. Reaction probabilities and integral cross sections (ICSs) between a collisional energy of 10(-4) eV and 0.1 eV have been calculated and compared with previously reported results of a time independent quantum mechanical (TIQM) approach. The TIQM results exhibit a dense profile with numerous narrow resonances down to Ec ~ 10(-2) eV and for the case of H2(v = 0, j = 0) a prominent peak is found at ~2.5 × 10(-4) eV. The analysis at the state-to-state level reveals that this feature is originated in those processes which yield the formation of rotationally excited HD(v' = 0, j' > 0). The statistical predictions reproduce reasonably well the overall behaviour of the TIQM ICSs at the larger energy range (Ec ≥ 10(-3) eV). Thermal rate constants are in qualitative agreement for the whole range of temperatures investigated in this work, 10-100 K, although the SQM values remain above the TIQM results for both initial H2 rotational states, j = 0 and 1. The enlargement of the asymptotic region for the statistical approach is crucial for a proper description at low energies. In particular, we find that the SQM method leads to rate coefficients in terms of the energy in perfect agreement with previously reported measurements if the maximum distance at which the calculation is performed increases noticeably with respect to the value employed to reproduce the TIQM results.

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.

A classical trajectory study of the intramolecular dynamics, isomerization, and unimolecular dissociation of HO2

The classical dynamics and rates of isomerization and dissociation of HO2 have been studied using two potential energy surfaces (PESs) based on interpolative fittings of ab initio data: An interpolative moving least-squares (IMLS) surface [A. Li, D. Xie, R. Dawes, A. W. Jasper, J. Ma, and H. Guo, J. Chem. Phys. 133, 144306 (2010)] and the cubic-spline-fitted PES reported by Xu, Xie, Zhang, Lin, and Guo (XXZLG) [J. Chem. Phys. 127, 024304 (2007)]. Both PESs are based on similar, though not identical, internally contracted multi-reference configuration interaction with Davidson correction (icMRCI+Q) electronic structure calculations; the IMLS PES includes complete basis set (CBS) extrapolation. The coordinate range of the IMLS PES is limited to non-reactive processes. Surfaces-of-section show similar generally regular phase space structures for the IMLS and XXZLG PESs with increasing energy. The intramolecular vibrational energy redistribution (IVR) at energies above and below the threshold of isomerization is slow, especially for O-O stretch excitations, consistent with the regularity in the surfaces-of-section. The slow IVR rates lead to mode-specific effects that are prominent for isomerization (on both the IMLS and XXZLG) and modest for unimolecular dissociation to H + O2 (accessible only on the XXZLG PES). Even with statistical distributions of initial energy, slow IVR rates result in double exponential decay for isomerization, with the slower rate correlated with slow IVR rates for O-O vibrational excitation. The IVR and isomerization rates computed for the IMLS and XXZLG PESs are quantitatively, but not qualitatively, different from one another with the largest differences ascribed to the ~2 kcal/mol difference in the isomerization barrier heights. The IMLS and XXZLG results are compared with those obtained using the global, semi-empirical double-many-body expansion DMBE-IV PES [M. R. Pastrana, L. A. M. Quintales, J. Brandão, and A. J. C. Varandas, J. Chem. Phys. 94, 8073 (1990)], for which the surfaces-of-section display more irregular phase space structure, much faster IVR rates, and significantly less mode-specific effects in isomerization and unimolecular dissociation. The calculated IVR results for all three PESs are reasonably well represented by an analytic, coupled three-mode energy transfer model.

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

We report extensive, fully quantum, time-independent (TID) calculations of cross sections at low collision energies and rate constants at low temperatures for the O+OH reaction, of key importance in the production of molecular oxygen in cold, dark, interstellar clouds and in the chemistry of the Earth's atmosphere. Our calculations are compared with TID calculations within the J-shifting approximation, with wave-packet calculations, and with quasiclassical trajectory calculations. The fully quantum TID calculations yield rate constants higher than those from the more approximate methods and are qualitatively consistent with a low-temperature extrapolation of earlier experimental values but not with the most recent experiments at the lowest temperatures.