Theoretical Kinetics Study of the F(2P) + NH3 Hydrogen Abstraction Reaction

Accurate ab initio based potential energy surface and kinetics of the Cl + NH3 → HCl + NH2 reaction

The gas-phase reaction Cl + NH3 → HCl + NH2 is a prototypical hydrogen abstraction reaction, whose minimum energy path involves several intermediate complexes. In this work, a full-dimensional, spin-orbit corrected potential energy surface (SOC PES) is constructed for the ground electronic state of the Cl + NH3 reaction. About 52 000 energy points are sampled and calculated at the UCCSD(T)-F12a/aug-cc-pVTZ level, in which the data points located in the entrance channel are spin-orbit corrected. The spin-orbit corrections are predicted by a fitted three-dimensional energy surface from about 7520 energy points in the entrance channel at the level of CASSCF (15e, 11o)/aug-cc-pVTZ. The fundamental-invariant neural network method is utilized to fit the SOC PES, resulting in a total root mean square error of 0.12 kcal mol-1. The calculated thermal rate constants of the Cl + NH3 → HCl + NH2 reaction on the SOC PES with the soft-zero-point energy constraint agree reasonably well with the available experimental values.

Kinetic study of the CN + C2H6 hydrogen abstraction reaction based on an analytical potential energy surface

The temperature dependence of the thermal rate constants and kinetic isotope effects (KIE) of the CN + C2H6 gas-phase hydrogen abstraction reaction was theoretically determined within the 25-1000 K temperature range, i.e., from very low- to high-temperature regimes. Based on a recently developed full-dimensional analytical potential energy surface fitted to highly accurate explicitly correlated ab initio calculations, three different kinetic theories were used: canonical variational transition state theory (CVT), quasiclassical trajectory theory (QCT), and ring polymer molecular dynamics (RPMD) method for the computation of rate constants. We found that the thermal rate constants obtained with the three theories show a V-shaped temperature dependence, with a pronounced minimum near 200 K, qualitatively reproducing the experimental measurements. Among the three methods used in this work, the QCT and RPMD methods have the best agreement with the experiment at low and high temperatures, respectively, while the CVT model shows the largest discrepancies. The significant increase in the rate constant at very low temperatures in this very exothermic and practically barrierless reaction could be attributed to the large value of the impact parameter, possibly ruling out the role of the tunneling effect and the intermediate complexes in the entrance channel. The theoretical H/D KIE depicted a "normal" behaviour, i.e., values greater than unity, emulating the experimental measurements and improving previous theoretical results. Finally, the discrepancies between theory and experiments were analysed as a function of several factors, such as limitations of the kinetics theories and the potential energy surface, as well as the uncertainties in the experimental measurements.

Theoretical Insights into the Dynamics of Gas-Phase Bimolecular Reactions with Submerged Barriers

Much attention has been paid to the dynamics of both activated gas-phase bimolecular reactions, which feature monotonically increasing integral cross sections and Arrhenius kinetics, and their barrierless capture counterparts, which manifest monotonically decreasing integral cross sections and negative temperature dependence of the rate coefficients. In this Perspective, we focus on the dynamics of gas-phase bimolecular reactions with submerged barriers, which often involve radicals or ions and are prevalent in combustion, atmospheric chemistry, astrochemistry, and plasma chemistry. The temperature dependence of the rate coefficients for such reactions is often non-Arrhenius and complex, and the corresponding dynamics may also be quite different from those with significant barriers or those completely dominated by capture. Recent experimental and theoretical studies of such reactions, particularly at relatively low temperatures or collision energies, have revealed interesting dynamical behaviors, which are discussed here. The new knowledge enriches our understanding of the dynamics of these unusual reactions.

Dynamics study of the post-transition-state-bifurcation process of the (HCOOH)H+ → CO + H3O+/HCO+ + H2O dissociation: application of machine-learning techniques

The process of protonated formic acid dissociating from the transition state was studied using ring-polymer molecular dynamics (RPMD), classical MD, and quasi-classical trajectory (QCT) simulations. Temperature had a strong influence on the branching fractions for the HCO⁺ + H₂O and CO + H₃O⁺ dissociation channels. The RPMD and classical MD simulations showed similar behavior, but the QCT dynamics were significantly different owing to the excess energies in the quasi-classical trajectories. Machine-learning analysis identified several key features in the phase information of the vibrational motions at the transition state. We found that the initial configuration and momentum of a hydrogen atom connected to a carbon atom and the shrinking coordinate of the CO bond at the transition state play a role in the dynamics of HCO⁺ + H₂O production.

Observation of resonances in the transition state region of the F + NH3 reaction using anion photoelectron spectroscopy

The transition state of a chemical reaction is a dividing surface on the reaction potential energy surface (PES) between reactants and products and is thus of fundamental interest in understanding chemical reactivity. The transient nature of the transition state presents challenges to its experimental characterization. Transition-state spectroscopy experiments based on negative-ion photodetachment can provide a direct probe of this region of the PES, revealing the detailed vibrational structure associated with the transition state. Here we study the F + NH₃ → HF + NH₂ reaction using slow photoelectron velocity-map imaging spectroscopy of cryogenically cooled FNH₃^- anions. Reduced-dimensionality quantum dynamical simulations performed on a global PES show excellent agreement with the experimental results, enabling the assignment of spectral structure. Our combined experimental-theoretical study reveals a manifold of vibrational Feshbach resonances in the product well of the F + NH₃ PES. At higher energies, the spectra identify features attributed to resonances localized across the transition state and into the reactant complex that may impact the bimolecular reaction dynamics.

Current Status of the X + C2H6 [X ≡ H, F(2P), Cl(2P), O(3P), OH] Hydrogen Abstraction Reactions: A Theoretical Review

This paper is a detailed review of the chemistry of medium-size reactive systems using the following hydrogen abstraction reactions with ethane, X + C₂H₆ → HX + C₂H₅; X ≡ H, F(²P), Cl(²P), O(³P) and OH, and focusing attention mainly on the theoretical developments. These bimolecular reactions range from exothermic to endothermic systems and from barrierless to high classical barriers of activation. Thus, the topography of the reactive systems changes from reaction to reaction with the presence or not of stabilized intermediate complexes in the entrance and exit channels. The review begins with some reflections on the inherent problems in the theory/experiment comparison. When one compares kinetics or dynamics theoretical results with experimental measures, one is testing both the potential energy surface describing the nuclei motion and the kinetics or dynamics method used. Discrepancies in the comparison may be due to inaccuracies of the surface, limitations of the kinetics or dynamics methods, and experimental uncertainties that also cannot be ruled out. The paper continues with a detailed review of some bimolecular reactions with ethane, beginning with the reactions with hydrogen atoms. The reactions with halogens present a challenge owing to the presence of stabilized intermediate complexes in the entrance and exit channels and the influence of the spin-orbit states on reactivity. Reactions with O(³P) atoms lead to three surfaces, which is an additional difficulty in the theoretical study. Finally, the reactions with the hydroxyl radical correspond to a reactive system with ten atoms and twenty-four degrees of freedom. Throughout this review, different strategies in the development of analytical potential energy surfaces describing these bimolecular reactions have been critically analyzed, showing their advantages and limitations. These surfaces are fitted to a large number of ab initio calculations, and we found that a huge number of calculations leads to accurate surfaces, but this information does not guarantee that the kinetics and dynamics results match the experimental measurements.

Capturing quantum effects with quasi-classical trajectories in the D + H+3 → H2D+ + H reaction

We present quasi-classical trajectory (QCT) cross sections, rate constants, and product state distributions for the D + H+3 → H₂D⁺ + H reaction. Using the same H+4 potential surface, the rate constants obtained from several QCT-based methods correcting for zero-point effects by Gaussian binning the product H₂D⁺ are compared to ring polymer molecular dynamics (RPMD) rate constants [Bulut et al., J. Phys. Chem. A, 2019, 123, 8766] which include quantum effects and to recent experimentally derived rate constants [Bowen et al., J. Chem. Phys., 2021, 154, 084307]. QCT with standard binning predicts rate constants that increase slowly as the temperature decreases from 1500 to 100 K. In contrast, the RPMD rate constants decrease rapidly with decreasing temperature. By 100 K, the QCT standard binning rate constant is more than 3 orders of magnitude larger than the RPMD rate constant. We show that QCT with Gaussian binning and proper normalization captures the zero-point effects and reproduces the RPMD rate constants over a large temperature range. Furthermore, the simple technique of counting only reactive trajectories with vibrational energy above the product zero-point energy matches the RPMD results well down to ∼300 K. The present Gaussian binned rate constants are in fair agreement with new experimentally derived rate constants from 100 to 1500 K. However, because the Gaussian binned rate constants do not include tunneling, important at lower temperatures, and the RPMD and experimentally derived rate constants have significant differences, the roles of the competing effects of zero-point energy, internal excitation of the H+3, and quantum tunneling are not simple and require further study for a consistent picture of the dynamics. Since rate constants for complex forming reactions, such as the title reaction, are difficult to converge with RPMD, alternative QCT-based methods, which include quantum effects and in addition provide product state distributions as described here, are highly desirable.

Ring-Polymer Molecular Dynamics Calculations of Thermal Rate Coefficients and Branching Ratios for the Interstellar H3+ + CO → H2 + HCO+/HOC+ Reaction and Its Deuterated Analogue

The reaction between H₃⁺ and CO is important in understanding the H₃⁺ destruction mechanism in the interstellar medium. In this work, thermal rate coefficients for the H₃⁺ + CO and D₃⁺ + CO reactions are calculated using ring-polymer molecular dynamics (RPMD) on a high-level machine-learning potential energy surface. The RPMD results agree well with the classical molecular dynamics results, where nuclear quantum effects are completely ignored, whereas the agreement between the RPMD results and the previous quasi-classical trajectory is good only at low temperatures. The calculated [HCO⁺]/[HOC⁺] product branching ratios decrease as the temperature increases, and the product branching is exclusively determined by the initial collisional orientation, which governs the formation of an ion-dipole complex, H₃⁺···CO or H₃⁺···OC, that dissociates into H₂ + HCO⁺ or H₂ + HOC⁺, respectively, via a direct mechanism. However, the contribution of the indirect mechanism via the rearrangement between H₃⁺···CO and H₃⁺···OC increases as the temperature increases, although its absolute fraction is small.

Effects of bending excitation on the reaction dynamics of fluorine atoms with ammonia

Vibrational excitation has been established as an efficient way to control the chemical reaction outcome. Stretching vibration of polyatomic molecules is believed to be efficient to promote abstraction reactions since energy is placed directly into the breaking bond. In this work, we report on a counterexample showing that exciting the low-frequency umbrella bending mode of ammonia enhances its reaction with fluorine atoms much more than exciting the high-frequency symmetric or asymmetric stretching mode over a wide range of collision energy, validated using both quasiclassical trajectory simulations and full-dimensional quantum dynamics calculations under the centrifugal-sudden approximation. This interesting mode-specific reaction dynamic originates from the increased chance of capturing the fluorine atom by ammonia due to the enlarged attractive interaction between them and the enhancement of the direct stripping reaction mediated by two submerged barriers.

Final-State-Resolved Dynamics of the H3+ + CO → H2 +HCO+/HOC+ Reaction: A Quasi-Classical Trajectory Study

The ion-molecule reaction H₃⁺ + CO → H₂ + HCO⁺/HOC⁺, which initiates the formation of crucial organic molecules, plays a key role in interstellar and circumstellar environments. In this work, the quasi-classical trajectory method is employed to study the reaction dynamics on a recently developed full-dimensional global potential energy surface (PES). The calculated product internal energy distributions and relative internal excited fractions agree reasonably well with the experimental measurements. For the two reaction channels, most of the available energy flows into the vibrational modes of HCO⁺ or HOC⁺ at low collision energies, followed by the translational mode and the rotational modes of HCO⁺ or HOC⁺. As the collision energy increases, the proportion of the product translational energy increases while the proportion of the product vibrational energy decreases. Furthermore, the CH and CO stretching modes and their combination bands are effectively excited for the product HCO⁺ while the bending mode is remarkably excited for the product HOC⁺.

VTST and RPMD kinetics study of the nine-body X + C2H6 (X ≡ H, Cl, F) reactions based on analytical potential energy surfaces

Thermal rate constants of nine-atom hydrogen abstraction reactions, X + C2H6 → HX + C2H5 (X ≡ H, Cl, F) with qualitatively different reaction paths, have been investigated using two kinetics approaches - variational transition state theory with multidimensional tunnelling (VTST/MT) and ring polymer molecular dynamics (RPMD) - and full dimensional analytical potential energy surfaces. For the H + C2H6 reaction, which proceeds through a noticeable barrier height of 11.62 kcal mol-1, kinetics approaches showed excellent agreement between them (with differences less than 30%) and with the experiment (with differences less than 60%) in the wide temperature range of 200-2000 K. For X = Cl and F, however, the situation is very different. The barrier height is either low or very low, 2.44 and 0.23 kcal mol-1, respectively, and the presence of van der Waals complexes in the entrance channel leads to a very flat topography and, consequently, imposes theoretical challenges. For the Cl(2P) reaction, VTST/MT underestimates the experimental rate constants (with differences less than 86%), and RPMD demonstrates better agreement (with differences less than 47%), although the temperature dependence is opposite to the experiment at low temperatures. Finally, for the F(2P) reaction, available experimental information shows discrepancies, both in the absolute values of the rate constants and also in the temperature dependence. Unfortunately, kinetics theories did not resolve this discrepancy. Different possible causes of these theory/experiment discrepancies were analyzed.

An accurate potential energy surface and ring polymer molecular dynamics study of the Cl + CH4→ HCl + CH3reaction

Thermal rate coefficients for the Cl + CH₄/CD₄reactions were studied on a new full-dimensional accurate potential energy surface with the spin–orbit corrections considered in the entrance channel.

Quantum Effects on the D + H3+ → H2D+ + H Deuteration Reaction and Isotopic Variants

The title reaction and its isotopic variants are studied using quasi-classical trajectory (QCT) (without taking into account corrections to account for the possible zero point energy breakdown) and ring polymer molecular dynamics (RPMD) methods with a full dimensional and accurate potential energy surface which presents an exchange barrier of approximately 0.144 eV. The QCT rate constant increases when the temperature decreases from 1500 to 10 K. On the contrary, the RPMD rate constant decreases with decreasing temperature, in semiquantitative agreement with recent experimental results. The present RPMD results are in between the thermal and translational experimental rate constants, extracted from the measured data to eliminate the initial vibrational excitation of H₃⁺, obtained in an arc discharge. The difference between the present RPMD results and experimental values is attributed to the possible existence of non thermal vibrational excitation of H₃⁺, not completely removed by the semiempirical model used for the analysis of the experimental results. Also, it is found that, below 200 K, the RPMD trajectories are trapped, forming long-lived collision complexes, with lifetimes longer than 1 ns. These collision complexes can fragment by either redissociating back to reactants or react to products, in the two cases tunneling through the centrifugal and reaction barriers, respectively. The contribution of the formation of the complex to the total deuteration rate should be calculated with more accurate quantum methods, as has been found recently for reactions of larger systems, and the present four atoms system is a good candidate to benchmark the adequacy of RPMD method at temperatures below 100 K.

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.

Theoretical study of the F(2P) + NH3→ HF + NH2 reaction on an accurate potential energy surface: dynamics and kinetics

The highly exothermic hydrogen abstraction reaction of the F atom with NH₃ is investigated using the quasi-classical trajectory method on a newly developed potential energy surface (PES) for the ground electronic state. The full-dimensional PES is constructed by fitting 41 282 ab initio energy points at the level of UCCSD(T)-F12/aug-cc-pVTZ. The flexible fundamental invariant-neural network method is applied in the fitting, resulting in a total root mean square error of 0.13 kcal mol^-1. On one hand, the calculated differential cross sections agree reasonably well with the experimental results and indicate that the reaction is dominated by the direct abstraction and stripping mechanisms while a considerable amount of reaction takes place by the indirect "yo-yo" mechanism. The product energy partition also reproduces well the experimental result, which can be understood according to the geometry change along the minimum energy path. On the other hand, the obtained vibrational state distribution of the product HF follows P_νHF=2≈P_νHF=1 > P_νHF=0 > P_νHF=3, less consistent with the scattered experimental results. In addition, the calculated thermal rate coefficients have practically no temperature dependence within the statistical errors.

The hydrogen abstraction reaction H + C2H6 → H2(v,j) + C2H5. Part II. Theoretical kinetics and dynamics study

Two important issues motivated the present study: the role of the tunnelling contribution at low temperatures and the role of the alkyl fragment in the dynamics.

The low temperature D+ + H2→ HD + H+ reaction rate coefficient: a ring polymer molecular dynamics and quasi-classical trajectory study

The reaction between D⁺ and H₂ plays an important role in astrochemistry at low temperatures and also serves as a prototype for a simple ion-molecule reaction. Its ground X[combining tilde]¹A' state has a very small thermodynamic barrier (up to 1.8 × 10^-2 eV) and the reaction proceeds through the formation of an intermediate complex lying within the potential well with a depth of at least 0.2 eV, thus representing a challenge for dynamical studies. In the present work, we analyze the title reaction within the temperature range of 20-100 K by means of ring polymer molecular dynamics (RPMD) and quasi-classical trajectory (QCT) methods over the full-dimensional global potential energy surface developed by Aguado et al. [A. Aguado, O. Roncero, C. Tablero, C. Sanz and M. Paniagua, J. Chem. Phys., 2000, 112, 1240]. The computed thermal RPMD and QCT rate coefficients are found to be almost independent of temperature and fall within the range of 1.34-2.01 × 10^-9 cm³ s^-1. They are also in very good agreement with previous time-independent quantum mechanical and statistical quantum method calculations. Furthermore, we observe that the choice of asymptotic separation distance between the reactants can markedly alter the rate coefficient in the low temperature regime (20-50 K). Therefore it is of utmost importance to correctly assign the value of this parameter for dynamical studies, particularly at very low temperatures of astrochemical importance. We finally conclude that the experimental rate measurements for the title reaction are highly desirable in future.

A combined theoretical and experimental investigation of the kinetics and dynamics of the O(1D) + D2 reaction at low temperature

The O(¹D) + H₂ reaction is a prototype for simple atom-diatom insertion type mechanisms considered to involve deep potential wells. While exact quantum mechanical methods can be applied to describe the dynamics, such calculations are challenging given the numerous bound quantum states involved. Consequently, efforts have been made to develop alternative theoretical strategies to portray accurately the reactive process. Here we report an experimental and theoretical investigation of the O(¹D) + D₂ reaction over the 50-296 K range. The calculations employ three conceptually different approaches - mean potential phase space theory, the statistical quantum mechanical method and ring polymer molecular dynamics. The calculated rate constants are in excellent agreement over the entire temperature range, exhibiting only weak temperature dependence. The agreement between experiment and theory is also very good, with discrepancies smaller than 26%. Taken together, the present and previous theoretical results validate the hypothesis that long-lived complex formation dominates the reaction dynamics at low temperature.

F(2P) + C2H6 → HF + C2H5 kinetics study based on a new analytical potential energy surface

An exhaustive kinetics study was performed for the title reaction using two theoretical approaches: variational transition-state theory and quasi-classical trajectory calculations, based on an original new analytical full-dimensional potential energy surface, named PES-2018, which has been fitted to high-level ab initio calculations.

Kinetics study of the CN + CH4 hydrogen abstraction reaction based on a new ab initio analytical full-dimensional potential energy surface

We have developed an analytical full-dimensional potential energy surface, named PES-2017, for the gas-phase hydrogen abstraction reaction between the cyano radical and methane.

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.

Ring-Polymer Molecular Dynamics for the Prediction of Low-Temperature Rates: An Investigation of the C((1)D) + H2 Reaction

Quantum mechanical calculations are important tools for predicting the rates of elementary reactions, particularly for those involving hydrogen and at low temperatures where quantum effects become increasingly important. These approaches are computationally expensive, however, particularly when applied to complex polyatomic systems or processes characterized by deep potential wells. While several approximate techniques exist, many of these have issues with reliability. The ring-polymer molecular dynamics method was recently proposed as an accurate and efficient alternative. Here, we test this technique at low temperatures (300-50 K) by analyzing the behavior of the barrierless C((1)D) + H2 reaction over the two lowest singlet potential energy surfaces. To validate the theory, rate coefficients were measured using a supersonic flow reactor down to 50 K. The experimental and theoretical rates are in excellent agreement, supporting the future application of this method for determining the kinetics and dynamics of a wide range of low-temperature reactions.

Stress Test for Quantum Dynamics Approximations: Deep Tunneling in the Muonium Exchange Reaction D + HMu → DMu + H

Quantum effects play a crucial role in chemical reactions involving light atoms at low temperatures, especially when a light particle is exchanged between two heavier partners. Different theoretical methodologies have been developed in the last decades attempting to describe zero-point energy and tunneling effects without abandoning a classical or semiclassical framework. In this work, we have chosen the D + HMu → DMu + H reaction as a stress test system for three well-established methods: two representative versions of transition state theory (TST), canonical variational theory and semiclassical instanton, and ring polymer molecular dynamics (RPMD). These calculations will be compared with accurate quantum mechanical results. Despite its apparent simplicity, the exchange of the extremely light muonium atom (0.114 u) becomes a most challenging reaction for conventional methods. The main result of this work is that RPMD provides an overall better performance than TST-based methods for such a demanding reaction. RPMD might well turn out to be a useful tool beyond TST applicability.

Quantum rate coefficients and kinetic isotope effect for the reaction Cl + CH4 → HCl + CH3 from ring polymer molecular dynamics

Thermal rate coefficients and kinetic isotope effect have been calculated for prototypical heavy-light-heavy polyatomic bimolecular reactions Cl + CH4/CD4 → HCl/DCl + CH3/CD3, using a recently proposed quantum dynamics approach: ring polymer molecular dynamics (RPMD). Agreement with experimental rate coefficients, which are quite scattered, is satisfactory. However, differences up to 50% have been found between the RPMD results and those obtained from the harmonic variational transition-state theory on one of the two full-dimensional potential energy surfaces used in the calculations. Possible reasons for such discrepancy are discussed. The present work is an important step in a series of benchmark studies aimed at assessing accuracy for RPMD for chemical reaction rates, which demonstrates that this novel method is a quite reliable alternative to previously developed techniques based on transition-state theory.

Ring-Polymer Molecular Dynamics Rate Coefficient Calculations for Insertion Reactions: X + H2 → HX + H (X = N, O)

The thermal rate constants of two prototypical insertion-type reactions, namely, N/O + H2 → NH/OH + H, are investigated with ring polymer molecular dynamics (RPMD) on full-dimensional potential energy surfaces using recently developed RPMDrate code. It is shown that the unique ability of the RPMD approach among the existing theoretical methods to capture the quantum effects, e.g., tunneling and zero-point energy, as well as recrossing dynamics quantum mechanically with ring-polymer trajectories leads to excellent agreement with rigorous quantum dynamics calculations. The present result is encouraging for future applications of the RPMD method and the RPMDrate code to complex-forming chemical reactions involving polyatomic reactants.