Rate Coefficients of the HCl + OH → Cl + H2O Reaction from Ring Polymer Molecular Dynamics

First-principles surface reaction rates by ring polymer molecular dynamics and neural network potential: role of anharmonicity and lattice motion

Elementary gas-surface processes are essential steps in heterogeneous catalysis. A predictive understanding of catalytic mechanisms remains challenging due largely to difficulties in accurately characterizing the kinetics of such steps. Experimentally, thermal rates for elementary surface reactions can now be measured using a novel velocity imaging technique, providing a stringent testing ground for ab initio rate theories. Here, we propose to combine ring polymer molecular dynamics (RPMD) rate theory with state-of-the-art first-principles-determined neural network potential to calculate surface reaction rates. Taking NO desorption from Pd(111) as an example, we show that the harmonic approximation and the neglect of lattice motion in the commonly-used transition state theory overestimates and underestimates the entropy change during the desorption process, respectively, leading to opposite errors in rate coefficient predictions and artificial error cancellations. Including anharmonicity and lattice motion, our results reveal a generally neglected surface entropy change due to significant local structural change during desorption and obtain the right answer for the right reasons. Although quantum effects are found to be less important in this system, the proposed approach establishes a more reliable theoretical benchmark for accurately predicting the kinetics of elementary gas-surface processes.

New Stable and Fast Ring-Polymer Molecular Dynamics for Calculating Bimolecular Rate Coefficients with an Example of OH + CH4

The accurate and efficient calculation of the rate coefficients of chemical reactions is a key issue in the research of chemical dynamics. In this work, by applying the dimension-free ultrastable Cayley propagator, the thermal rate coefficients of a prototypic high dimensional chemical reaction OH + CH₄ → H₂O + CH₃ in the temperature range of 200 to 1500 K are investigated with ring polymer molecular dynamics (RPMD) on a highly accurate full-dimensional potential energy surface. Kinetic isotope effects (KIEs) for three isotopologues of the title reaction are also studied. The results demonstrate excellent agreement with experimental data, even in the deep tunneling region. Especially, the Cayley propagator shows a high calculation efficiency with little loss of accuracy. The present results confirmed the applicability of the RPMD method, particularly the speed-up using a Cayley propagator, in theoretical calculations of bimolecular reaction rates.

Theory of vibrational polariton chemistry in the collective coupling regime

We theoretically demonstrate that the chemical reaction rate constant can be significantly suppressed by coupling molecular vibrations with an optical cavity, exhibiting both the collective coupling effect and the cavity frequency modification of the rate constant. When a reaction coordinate is strongly coupled to the solvent molecules, the reaction rate constant is reduced due to the dynamical caging effect. We demonstrate that collectively coupling the solvent to the cavity can further enhance this dynamical caging effect, leading to additional suppression of the chemical kinetics. This effect is further amplified when cavity loss is considered.

Strong non-Arrhenius behavior at low temperatures in the OH + HCl → H2O + Cl reaction due to resonance induced quantum tunneling

The OH + HCl reaction possesses many Feshbach resonances trapped in the hydrogen bond well in the entrance channel, which substantially enhance the reaction rates at low temperatures.

Theoretical H + O3 rate coefficients from ring polymer molecular dynamics on an accurate global potential energy surface: assessing experimental uncertainties

Thermal rate coefficients and kinetic isotope effects have been calculated for an important atmospheric reaction H/D + O3 → OH/OD + O2 based on an accurate permutation invariant polynomial-neural network potential energy surface, using ring polymer molecular dynamics (RPMD), quasi-classical trajectory (QCT) and variational transition-state theory (VTST) with multidimensional tunneling. The RPMD approach yielded results that are generally in better agreement with experimental rate coefficients than the VTST and QCT ones, especially at low temperatures, attributable to its capacity to capture quantum effects such as tunneling and zero-point energy. The theoretical results support one group of existing experiments over the other. In addition, rate coefficients for the D + O3 → OD + O2 reaction are also reported using the same methods, which will allow a stringent assessment of future experimental measurements, thus helping to reduce the uncertainty in the recommended rate coefficients of this reaction.

Kinetic study of the OH + HO2 → H2O + O2 reaction using ring polymer molecular dynamics and quantum dynamics

The reaction OH + HO2 → H2O + O2 is a prototype of radical-radical reactions. It plays an important role in interstellar/atmospheric chemistry and combustion, and considerable attention has thus been dedicated to its kinetics. In our previous work, we reported an accurate full-dimensional potential energy surface for the title reaction on the ground triplet electronic state. The quasi-classical trajectory (QCT) approach was employed to investigate its kinetics. Although the QCT rate coefficients were in good agreement with some experimental and theoretical results, QCT cannot account for the quantum mechanical effects, such as zero-point vibrational energy, recrossing, and tunneling, which may significantly affect the rate coefficients, particularly at low temperatures. In this work, the reduced-dimensional quantum dynamics and ring polymer molecular dynamics calculations were carried out to examine these effects and their impact on rate coefficients over the temperature range of 300-1300 K.

Probing the Partial Activation of Water by Open-Shell Interactions, Cl(H2O)1-4

The partial chemical activation of water by reactive radicals was examined computationally for small clusters of chlorine and water, Cl^•(H₂O)_n=1-4. Using an automated isomer-search procedure, dozens of unique, stable structures were computed. Among the resulting structural classes were intact, hydrated-chlorine isomers, as well as hydrogen-abstracted (HCl)(OH)(H₂O)_n-1 configurations. The latter showed increased stability as the degree of hydration increased, until n = 4, where a new class of structures was discovered with a chloride ion bound to an oxidized water network. The electronic structure of these three structural classes was investigated, and spectral signatures of this hydration-based evolution were connected to these electronic properties. An ancillary outcome of this detailed computational analysis, including coupled-cluster benchmarks, was the calibration of cost-effective quantum chemistry methods for future studies of these radical-water complexes.

New Stress Test for Ring Polymer Molecular Dynamics: Rate Coefficients of the O(3P) + HCl Reaction and Comparison with Quantum Mechanical and Quasiclassical Trajectory Results

In the past decade, ring polymer molecular dynamics (RPMD) has emerged as a very efficient method to determine thermal rate coefficients for a great variety of chemical reactions. This work presents the application of this methodology to study the O(³P) + HCl reaction, which constitutes a stringent test for any dynamical calculation due to rich resonant structure and other dynamical features. The rate coefficients, calculated on the ³A' and ³A″ potential energy surfaces (PESs) by Ramachandran and Peterson [ J. Chem. Phys. 2003 , 119 , 9590 ], using RPMD and quasiclassical trajectories (QCT) are compared with the existing experimental and the quantum mechanical (QM) results by Xie et al. [ J. Chem. Phys. 2005 122 , 014301 ]. The agreement is very good at T > 600 K, although RPMD underestimates rate coefficients by a factor between 4 and 2 in the 200-500 K interval. The origin of these discrepancies lies in the large contribution from tunneling on the ³A″ PES, which is enhanced by resonances due to quasibound states in the van der Waals wells. Although tunneling is fairly well accounted for by RPMD even below the crossover temperature, the effect of resonances, a long-time effect, is not included in the methodology. At the highest temperatures studied in this work, 2000-3300 K, the RPMD rate coefficients are somewhat larger than the QM ones, but this is shown to be due to limitations in the QM calculations and the RPMD are believed to be more reliable.

Effects of water, ammonia and formic acid on HO2 + Cl reactions under atmospheric conditions: competition between a stepwise route and one elementary step

Quantum chemical calculations at M06-2X and CCSD(T) levels of theory have been performed to investigate the effects of H₂O, NH₃, and HCOOH on the HO₂ + Cl → HCl + O₂ reaction. The results show that catalyzed reactions with three catalysts could proceed through two different mechanisms, namely a stepwise route and one elementary step, where the former reaction is more favorable than the latter. Meanwhile, for the stepwise route, a single hydrogen atom transfer pathway in the presence of all catalysts has more advantages than the respective double hydrogen atom transfer pathway. Then, the relative impacts of catalysts under tropospheric conditions were investigated by considering the temperature dependence of the rate constants and the altitude dependence of catalyst concentrations. The calculated results show that at 0 km altitude, the HO₂ + Cl → HCl + O₂ reaction with catalysts, such as H₂O, NH₃, or HCOOH, cannot compete with the reaction without a catalyst, as the effective rate constant with a catalyst is smaller by 2–6 orders of magnitude than the naked reaction within the temperature range 280–320 K. The calculated results also show that at altitudes of 5, 10 and 15 km, the effective rate constant of the HCOOH-catalyzed reaction increases obviously with an increase in altitude. At 15 km altitude, its value is up to 9.63 × 10⁻¹¹ cm³ per molecule per s, which is close to the corresponding value of the reaction without a catalyst, showing that the contribution of HCOOH to the HO₂ + Cl → HCl + O₂ reaction cannot be neglected at high altitudes. The new findings in this investigation are not only of great necessity and importance for elucidating the gas-phase reaction of HO₂ with Cl in the presence of acidic, neutral and basic catalysts, but are also of great interest for understanding the importance of other types of hydrogen abstraction in the atmosphere.

The effects of acidic (FA), neutral (WM) and basic (AM) catalysts on the energetic and kinetic aspects of the HO₂ + Cl reaction have been studied. At 298 K, the catalytic order of FA, WM and AM is WM > FA > AM.

Kinetics of the OH+HCl→H2 O+Cl reaction: Rate determining roles of stereodynamics and roaming and of quantum tunneling

The OH + HCl → H₂ O + Cl reaction is one of the most studied four-body systems, extensively investigated by both experimental and theoretical approaches. Here, as a continuation of our previous work on the OH + HBr and OH + HI reactions, which manifest an anti-Arrhenius behavior that was explained by stereodynamic and roaming effects, we extend the strategy to understand the transition to the sub-Arrhenius behavior occurring for the HCl case. As previously, we perform first-principles on-the-fly Born-Oppenheimer molecular dynamics calculations, thermalized at four temperatures (50, 200, 350, and 500 K), but this time we also apply a high-level transition-state-theory, modified to account for tunneling conditions. We find that the theoretical rate constants calculated with Bell tunneling corrections are in good agreement with extensive experimental data available for this reaction in the ample temperature range: (i) simulations show that the roles of molecular orientation in promoting this reaction and of roaming in finding the favorable path are minor than in the HBr and HI cases, and (ii) dominating is the effect of quantum mechanical penetration through the energy barrier along the reaction path on the potential energy surface. The discussion of these results provides clarification of the origin on different non-Arrhenius mechanisms observed along this series of reactions. © 2018 Wiley Periodicals, Inc.

Quantum dynamics of ClH2O- photodetachment: Isotope effect and impact of anion vibrational excitation

Photodetachment of the ClH₂O^- anion is investigated using full-dimensional quantum mechanics on accurate potential energy surfaces of both the anion and neutral species. Detailed analysis of the photoelectron spectrum and the corresponding wavefunctions reveals that the photodetachment leads to, in the product channel of the exothermic HCl + OH → Cl + H₂O reaction, the formation of numerous Feshbach resonances due apparently to slow energy transfer from H₂O vibrational modes to the dissociation coordinate. These long-lived resonances can be grouped into two broad peaks in the low-resolution photoelectron spectrum, which is in good agreement with available experiments, and they are assigned to the ground and first excited OH stretching vibrational manifolds of H₂O complexed with Cl. In addition, effects of isotope substitution on the photoelectron spectrum were small. Finally, photodetachment of the vibrationally excited ClH₂O^- in the ionic hydrogen bond mode is found to lead to Feshbach resonances with higher stretching vibrational excitations in H₂O.

Automated calculation of thermal rate coefficients using ring polymer molecular dynamics and machine-learning interatomic potentials with active learning

We propose a methodology for the fully automated calculation of thermal rate coefficients of gas phase chemical reactions, which is based on combining ring polymer molecular dynamics (RPMD) and machine-learning interatomic potentials actively learning on-the-fly.

Accurate Determination of Tunneling-Affected Rate Coefficients: Theory Assessing Experiment

The thermal rate coefficients of a prototypical bimolecular reaction are determined on an accurate ab initio potential energy surface (PES) using ring polymer molecular dynamics (RPMD). It is shown that quantum effects such as tunneling and zero-point energy (ZPE) are of critical importance for the HCl + OH reaction at low temperatures, while the heavier deuterium substitution renders tunneling less facile in the DCl + OH reaction. The calculated RPMD rate coefficients are in excellent agreement with experimental data for the HCl + OH reaction in the entire temperature range of 200-1000 K, confirming the accuracy of the PES. On the other hand, the RPMD rate coefficients for the DCl + OH reaction agree with some, but not all, experimental values. The self-consistency of the theoretical results thus allows a quality assessment of the experimental data.

A global coupled cluster potential energy surface for HCl + OH ↔ Cl + H2O

A new and more accurate full-dimensional global potential energy surface (PES) for the ground electronic state of the ClH₂O system is developed by using the permutation invariant polynomial-neural network (PIP-NN) method to fit 15 777 points obtained using an explicitly correlated unrestricted coupled-cluster method with single, double, and perturbative triple excitations (UCCSD(T)-F12b).

An experimental and theoretical investigation of the C(1D) + D2 reaction

Rate constants derived from ring polymer molecular dynamics calculations confirm the validity of this method for studying low-temperature complex-forming reactions

A ring polymer molecular dynamics study of the OH + H2(D2) reaction

Using ring polymer molecular dynamics we have calculated the rate coefficients for the OH + H₂ reaction.