Accurate combined-hyperbolic-inverse-power-representation of ab initio potential energy surface for the hydroperoxyl radical and dynamics study of 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.

Accurate diabatization based on combined-hyperbolic-inverse-power-representation: 1,2 2A' states of BeH2

A diabatic potential energy matrix (DPEM) for the two lowest states of BeH2+ has been constructed using the combined-hyperbolic-inverse-power-representation (CHIPR) method. By imposing symmetry constraints on the coefficients of polynomials, the complete nuclear permutation inversion symmetry is correctly preserved in the CHIPR functional form. The symmetrized CHIPR functional form is then used in the diabatization by ansatz procedure. The ab initio energies are reproduced with satisfactory accuracy. In addition, the CHIPR-based DPEM also reproduces the local topology of a conical intersection. Future work will focus on a complete four-state diabatic representation with emphasis on the long-range interactions and spin-orbit couplings, which will enable accurate quantum scattering calculations for the Be+(2P) + H2 → BeH+(X1Σ+) + H(2S) reaction.

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.

Reconciling spectroscopy with dynamics in global potential energy surfaces: the case of the astrophysically relevant SiC_{2}

SiC_{2} is a fascinating molecule due to its unusual bonding and astrophysical importance. In this work, we report the first global potential energy surface (PES) for ground-state SiC_{2} using the combined-hyperbolic-inverse-power-representation (CHIPR) method and accurate ab initio energies. The calibration grid data is obtained via a general dual-level protocol developed afresh herein that entails both coupled-cluster and multireference configuration interaction energies jointly extrapolated to the complete basis set limit. Such an approach is specially devised to recover much of the spectroscopy from the PES, while still permitting a proper fragmentation of the system to allow for reaction dynamics studies. Besides describing accurately the valence strongly-bound region that includes both the cyclic global minimum and isomerization barriers, the final analytic PES form is shown to properly reproduce dissociation energies,diatomic potentials, and long-range interactions at all asymptotic channels, in addition to naturally reflect the correct permutational symmetry of the potential. Bound vibrational state calculations have been carried out, unveiling an excellent match of the available experimental data on c-SiC_{2}(^{1}A_{1}). To further exploit the global nature of the PES, exploratory quasi-classical trajectory calculations for the endothermic C2 +Si → SiC+C reaction are also performed, yielding thermalized rate coefficients for temperatures up to 5000 K. The results hint for the prominence of this reaction in the innermost layers of the circumstellar envelopes around carbon-rich stars, thence conceivably playing therein a key contribution to the gas-phase formation of SiC, and eventually, solid SiC dust.

From six to eight Π-electron bare rings of group-XIV elements and beyond: can planarity be deciphered from the "quasi-molecules" they embed?

Ab initio molecular orbital theory is used to study the structures of six and eight π-electron bare rings of group-XIV elements, and even larger [n]annulenes up to C₁₈H₁₈, including some of their mono-, di-, tri-, and tetra-anions. While some of the above rings are planar, others are nonplanar. A much spotlighted case is cyclo-octatetraene (C₈H₈), which is predicted to be nonplanar together with its heavier group-XIV analogues Si₈H₈ and Ge₈H₈, with the solely planar members of its family having the stoichiometric formulas C₄Si₄H₈ and C₄Ge₄H₈. A similar situation arises with the six π-electron bare rings, where benzene and substituted ones up to C₃Si₃H₆ or so are planar, while others are not. However, the explanations encountered in the literature find support in ab initio calculations for such species, often rationalized from distinct calculated features. Using second-order Møller-Plesset perturbation theory and, when affordable (particularly tetratomics, which may allow even higher levels), the coupled-cluster method including single, double, and perturbative triple excitations, a common rationale is suggested based on a novel concept of quasi-molecules or the (3+4)-atom partition scheme. Any criticism of tautology is therefore avoided. The same analysis has also been successfully applied to even larger [n]annulenes, to their mixed family members involving silicon and germanium atoms, and to the C₁₈ carbon ring. Furthermore, it has been extended to annulene anions to check the criteria of the popular Hückel rule for planarity and aromaticity. Exploratory work on cycloarenes is also reported. Besides a partial study of the involved potential energy surfaces, equilibrium geometries and harmonic vibrational frequencies have been calculated anew, for both the parent and the actual prototypes of the quasi-molecules.

An accurate NH2(X2A′′) CHIPR potential energy surface via extrapolation to the complete basis set limit and dynamics of the N(2D) + H2(X1+g) reaction

An accurate CHIPR potential energy surface for NH2(X2A′′) is structured to study the N(2D) + H2(X1Σ+g) reaction using the time-dependent wave packet and quasi-classical trajectory method.

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.

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.

A global CHIPR potential energy surface for ground-state C3H and exploratory dynamics studies of reaction C2 + CH → C3 + H

A full-dimensional global potential-energy surface (PES) is first reported for ground-state doublet C₃H using the combined-hyperbolic-inverse-power-representation (CHIPR) method and accurate ab initio energies extrapolated to the complete basis set limit. The PES is based on a many-body expansion-type development where the two-body and three-body energy terms are from our previously reported analytic potentials for C₂H(²A') and C₃(¹A',³A'), while the effective four-body one is calibrated using an extension of the CHIPR formalism for tetratomics. The final form is shown to accurately reproduce all known stationary structures of the PES, some of which are unreported thus far, and their interconversion pathways. Moreover, it warrants by built-in construction the appropriate permutational symmetry and describes in a physically reasonable manner all long-range features and the correct asymptotic behavior at dissociation. Exploratory quasi-classical trajectory calculations for the reaction C₂ + CH → C₃ + H are also performed, yielding thermalized rate coefficients for temperatures up to 4000 K.

Accurate CHIPR Potential Energy Surface for the Lowest Triplet State of C3

We report the first global ab initio-based potential energy surface (PES) for ground-state triplet C₃(³A') based on accurate energies extrapolated to the complete basis set (CBS) limit, and using the combined-hyperbolic-inverse-power-representation method for the analytical modeling. By relying on a cost-effective CBS(D,T) protocol, we ensure that the final form reproduces all topographical features of the PES, including its cyclic-linear isomerization barrier, with CBS(5,6)-quality. To partially account for the incompleteness of the N-electron basis and other minor effects, the available accurate experimental data on the relevant diatomics were used to obtain direct-fit curves that replace the theoretical ones in the many-body expansion. Besides describing properly long-range interactions at all asymptotic channels and permutational symmetry by built-in construction, the PES reported here reproduces the proper exothermicities at dissociation regions as well as the spectroscopy of the diatomic fragments. Bound vibrational state calculations in both linear and cyclic isomers have also been carried out, unveiling a good match of the available data on C₃(ã ³Π_u), while assisting with IR band positions for C₃(³A₂^') that may serve as a guide for its laboratory and astronomical detection.

Global Potential Energy Surface for HO2+ Using the CHIPR Method

An analytical potential energy function for the title ion based on the combined hyperbolic inverse power representation (CHIPR) method and its characteristics are discussed at length in the present work. The curves of two diatomic ions, O₂⁺ and OH⁺, are also obtained within the same approach. The model PES so obtained exhibits extraordinary flexibility in describing with subchemical accuracy even the weak topological features near the higher energy regions. Thus, structural properties predicted by the model may help spectroscopists who want to compare their experimental values with the ones from theory. The relaxed PESs in various coordinates have been calculated by relaxing the O₂ bond distance using the present model, thus throwing light on all the possible isomers and their interconversions. The latest estimates of IR frequencies for three vibrational modes have been compared with the computed frequencies using the present model, and the agreement seems encouraging.

Fully coupled (J > 0) time-dependent wave-packet calculations using hyperspherical coordinates for the H + O2 reaction on the CHIPR potential energy surface

ICS calculation by time dependent wavepacket approach for H + O₂ reaction using non-zero J values.

Modelling potential energy surfaces for small clusters using Shepard interpolation with Gaussian-form nodal functions

A new interpolation method based on Gaussian functions to reliably generate potential energy surfaces.

The geometric phase controls ultracold chemistry

The geometric phase is shown to control the outcome of an ultracold chemical reaction. The control is a direct consequence of the sign change on the interference term between two scattering pathways (direct and looping), which contribute to the reactive collision process in the presence of a conical intersection (point of degeneracy between two Born-Oppenheimer electronic potential energy surfaces). The unique properties of the ultracold energy regime lead to an effective quantization of the scattering phase shift enabling maximum constructive or destructive interference between the two pathways. By taking the O+OH→H+O2 reaction as an illustrative example, it is shown that inclusion of the geometric phase modifies ultracold reaction rates by nearly two orders of magnitude. Interesting experimental control possibilities include the application of external electric and magnetic fields that might be used to exploit the geometric phase effect reported here and experimentally switch on or off the reactivity.

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.

Dynamics of the O + ClO reaction: reactive and vibrational relaxation processes

Classical trajectories have been integrated to study the O + ClO reaction, both reactive and vibrational energy transfer processes, for the range of temperatures 100 ≤ T/K ≤ 500 using momentum Gaussian binning. The employed potential energy surface is the recently proposed single-sheeted double many-body expansion potential energy surface for the (2)A" ground-state of ClO2 based on multireference ab initio data. A capture-type regime with a room-temperature rate constant of (17.8 ± 0.5) × 10(-12) cm(3) s(-1) and temperature dependence of k(T/K)/cm(3) s(-1) = 22.4 × 10(-12) × T(-0.81) exp(-39.2/T) has been found. Although the value reported here is half of the experimental and recommended one, tentative explanations are given. Other dynamical attributes are also examined for the title reaction, with state-to-all and state-to-state vibrational relaxation and excitation rate constants reported for temperatures of relevance in stratospheric chemistry.

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.

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.