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

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.

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.

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.

Dynamical interference in the vibronic bond breaking reaction of HCO

First-principles treatments of quantum molecular reaction dynamics have reached the level of quantitative accuracy even in cases with strong non-Born-Oppenheimer effects. This achievement permits the interpretation of puzzling experimental phenomena related to dynamics governed by multiple coupled potential energy surfaces. We present a combined experimental and theoretical study of the photodissociation of formyl radical (HCO). Oscillations observed in the distribution of product states are found to arise from the interference of matter waves-a manifestation analogous to Young's double-slit experiment.

Dynamics of Complex-Forming Bimolecular Reactions: A Comparative Theoretical Study of the Reactions of H Atoms with O2((3)Σg(-)) and O2((1)Δg)

The atomic-level mechanism of the reaction of H atoms with triplet and singlet molecular oxygen, H((2)S) + O2((3)Σg(-)) → O((3)P) + OH((2)Πg) ( R1 ) and H((2)S) + O2((1)Δg) → O((3)P) + OH((2)Πg) ( R2 ) is analyzed in terms of the topology of the potential energy surfaces (PES) of the two reactions. Both PES exhibit a deep potential well corresponding to the ground and first excited electronic state of HO2. The ground-state reaction is endothermic with no barrier on either side of the well; the excited-state reaction is exothermic with a barrier in the entrance valley of the PES. The differences of the PES are manifested in properties such as the excitation functions, which show reaction R1 to be much slower and the effect of rotational excitation on reactivity, which speeds up reaction R1 and has little effect on R2 . Numerous common dynamics features arise from the presence of the deep potential well on the PES. Such are the significant role of isomerization (for example, 90% of reactive collisions in R2 involve at least one H atom transfer from one of the O atoms to the other in reaction R2 ), which is shown to give rise to a significant rotational excitation of the product OH radicals. Common is the significant sideways scattering of the products that originates from collisions in propeller-type arrangements induced by the presence of two bands of acceptance around the O2 molecule. The HO2 complex in both reactions proves to behave nonstatistically, with signatures of the dynamics in lifetime distributions, angular distributions, opacity functions, and product quantum-state distributions.

The cis- and trans-formylperoxy radical: fundamental vibrational frequencies and relative energies of the X̃ 2A′′ and à 2A′ states

Acylperoxy radicals [RC(O)OO˙] play an important catalytic role in many atmospheric and combustion 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 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.

Theoretical analysis of reaction kinetics with singlet oxygen molecules

A comparative analysis of predictive ability of three approaches to estimate the rate constants of reactions of H(2), H, H(2)O and CH(4) with electronically excited O(2)(a(1)Δ(g)) and O(2)(b(1)Σ(g)(+)) molecules is conducted. The first approach is based on a detailed ab initio study of potential energy surfaces. The second one is known as the "bond energy-bond order" method, and the third approach is a modification of the updated method of vibronic terms that makes it possible to evaluate the activation energy of reactions involving electronically excited species. The comparison showed that the estimates of the energy barrier by the updated method of vibronic terms for some reactions can be in good agreement with ab initio calculations and available experimental data. It was revealed that reactions of O(2)(b(1)Σ(g)(+)) molecules with H(2), H(2)O and CH(4) molecules and with the H atom result in the formation of electronically excited species. The reactivity of O(2)(b(1)Σ(g)(+)) molecules is smaller than that of O(2)(a(1)Δ(g)) ones, but much higher as compared to the reactivity of ground state O(2) molecules. For each reaction under study involving oxygen molecules in the excited electronic states O(2)(a(1)Δ(g)) and O(2)(b(1)Σ(g)(+)) the recommended temperature-dependent rate constants are presented.