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

Kinetic Study of the Reactions of Ground State Atomic Carbon and Oxygen with Nitrogen Dioxide over the 50-296 K Temperature Range

The kinetics of the reactions of nitrogen dioxide, NO2, with atomic oxygen and atomic carbon in their ground triplet states (3P) have been studied at room temperature and below using a supersonic flow (Laval nozzle) reactor. O(3P) and C(3P) atoms (hereafter O and C respectively) were created in situ by the pulsed laser photolysis of the precursor molecules NO2 at 355 nm and CBr4 at 266 nm, respectively. While the progress of the O + NO2 reaction was followed by detecting O atoms by a chemiluminescent tracer method, progress of the C + NO2 reaction was followed by detecting C atoms directly by vacuum ultraviolet laser-induced fluorescence at 116 nm. The measured rate constants for the O + NO2 reaction are found to be in excellent agreement with earlier work at higher temperatures and extend the available kinetic data for this process down to 50 K. The present work represents the first kinetics study of the C + NO2 reaction. Although both reactions display rate constants that increase as the temperature falls, a more substantial rate increase is observed for the O + NO2 reaction. The effects of these reactions on the simulated abundances of interstellar NO2 and related compounds were tested using a gas-grain model of the dense interstellar medium, employing expressions for the rate constants of the form, k(T) = α(T/300)β, with α = 1 × 10-11 cm3 s-1 and β = -0.65 for the O + NO2 reaction and α = 2 × 10-10 cm3 s-1 and β = -0.11 for the C + NO2 reaction. Although these simulations predict that gas-phase NO2 abundances are low in dense interstellar clouds, NO2 abundances on interstellar dust grains are predicted to reach reasonably high levels, indicating the potential for detection of this species in warmer regions.

Quantum rates in dissipative systems with spatially varying friction

We investigate whether making the friction spatially dependent on the reaction coordinate introduces quantum effects into the thermal reaction rates for dissipative reactions. Quantum rates are calculated using the numerically exact multi-configuration time-dependent Hartree method, as well as the approximate ring-polymer molecular dynamics (RPMD), ring-polymer instanton methods, and classical molecular dynamics. By conducting simulations across a wide range of temperatures and friction strengths, we can identify the various regimes that govern the reactive dynamics. At high temperatures, in addition to the spatial-diffusion and energy-diffusion regimes predicted by Kramer's rate theory, a (coherent) tunneling-dominated regime is identified at low friction. At low temperatures, incoherent tunneling dominates most of Kramer's curve, except at very low friction, when coherent tunneling becomes dominant. Unlike in classical mechanics, the bath's influence changes the equilibrium time-independent properties of the system, leading to a complex interplay between spatially dependent friction and nuclear quantum effects even at high temperatures. More specifically, a realistic friction profile can lead to an increase (or decrease) of the quantum (classical) rates with friction within the spatial-diffusion regime, showing that classical and quantum rates display qualitatively different behaviors. Except at very low frictions, we find that RPMD captures most of the quantum effects in the thermal reaction rates.

Reactions dynamics for X + H2 insertion reactions (X = C(1D), N(2D), O(1D), S(1D)) with Cayley propagator ring-polymer molecular dynamics

In this work, rate coefficients of four prototypical insertion reactions, X + H2 → H + XH (X = C(1D), N(2D), O(1D), S(1D)), and associated isotope reactions are calculated based on ring polymer molecular dynamics (RPMD) with Cayley propagator (Cayley-RPMD). The associated kinetic isotope effects are systematically studied too. The Cayley propagator used in this work increases the stability of numerical integration in RPMD calculations and also supports a larger evolution time interval, allowing us to reach both high accuracy and efficiency. So, our results do not only provide chemical kinetic data for the title reactions in an extended temperature range but also consist of experimental results, standard RPMD, and other theoretical methods. The results in this work also reflect that Cayley-RPMD has strong consistency and high reliability in its investigations of chemical dynamics for insertion 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.

Kinetic Study of the Gas-Phase O(1D) + CH3OH and O(1D) + CH3CN Reactions: Low-Temperature Rate Constants and Atomic Hydrogen Product Yields

Atomic oxygen in its first excited singlet state, O(¹D), is an important species in the photochemistry of several planetary atmospheres and has been predicted to be a potentially important reactive species on interstellar ices. Here, we report the results of a kinetic study of the reactions of O(¹D) with methanol, CH₃OH, and acetonitrile, CH₃CN, over the 50-296 K temperature range. A continuous supersonic flow reactor is used to attain these low temperatures coupled with pulsed laser photolysis and pulsed laser-induced fluorescence to generate and monitor O(¹D) atoms, respectively. Secondary experiments examining the atomic hydrogen product channels of these reactions are also performed, through laser-induced fluorescence measurements of H(²S) atom formation. On the kinetic side, the rate constants for these reactions are seen to be large (>2 × 10^-10 cm³ s^-1) and consistent with barrierless reactions, although they display contrasting dependences as a function of temperature. On the product formation side, both reactions are seen to yield non-negligible quantities of atomic hydrogen. For the O(¹D) + CH₃OH reaction, the derived yields are in good agreement with the conclusions of previous experimental and theoretical works. For the O(¹D) + CH₃CN reaction, whose H-atom formation channels had not previously been investigated, electronic structure calculations of several new product formation channels are performed to explain the observed H-atom yields. These calculations demonstrate the barrierless and exothermic nature of the relevant exit channels, confirming that atomic hydrogen is also an important product of the O(¹D) + CH₃CN reaction.

An Experimental and Theoretical Investigation of the Gas-Phase C(3P) + N2O Reaction. Low Temperature Rate Constants and Astrochemical Implications

The reaction between atomic carbon in its ground electronic state, C(³P), and nitrous oxide, N₂O, has been studied below room temperature due to its potential importance for astrochemistry, with both species considered to be present at high abundance levels in a range of interstellar environments. On the experimental side, we measured rate constants for this reaction over the 50-296 K range using a continuous supersonic flow reactor. C(³P) atoms were generated by the pulsed photolysis of carbon tetrabromide at 266 nm and were detected by pulsed laser-induced fluorescence at 115.8 nm. Additional measurements allowing the major product channels to be elucidated were also performed. On the theoretical side, statistical rate theory was used to calculate low temperature rate constants. These calculations employed the results of new electronic structure calculations of the ³A″ potential energy surface of CNNO and provided a basis to extrapolate the measured rate constants to lower temperatures and pressures. The rate constant was found to increase monotonically as the temperature falls (k_C(³P)+N2O (296 K) = (3.4 ± 0.3) × 10^-11 cm³ s^-1), reaching a value of k_C(³P)+N2O (50 K) = (7.9 ± 0.8) × 10^-11 cm³ s^-1 at 50 K. As current astrochemical models do not include the C + N₂O reaction, we tested the influence of this process on interstellar N₂O and other related species using a gas-grain model of dense interstellar clouds. These simulations predict that N₂O abundances decrease significantly at intermediate times (10³ - 10⁵ years) when gas-phase C(³P) abundances are high.

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.

Experimental and theoretical studies of the gas-phase reactions of O(1D) with H2O and D2O at low temperature

Here we report the results of an experimental and theoretical study of the gas-phase reactions between O(¹D) and H₂O and O(¹D) and D₂O at room temperature and below. On the experimental side, the kinetics of these reactions have been investigated over the 50-127 K range using a continuous flow Laval nozzle apparatus, coupled with pulsed laser photolysis and pulsed laser induced fluorescence for the production and detection of O(¹D) atoms respectively. Experiments were also performed at 296 K in the absence of a Laval nozzle. On the theoretical side, the existing full-dimensional ground X ¹A potential energy surface for the H₂O₂ system involved in this process has been reinvestigated and enhanced to provide a better description of the barrierless H-atom abstraction pathway. Based on this enhanced potential energy surface, quasiclassical trajectory calculations and ring polymer molecular dynamics simulations have been performed to obtain low temperature rate constants. The measured and calculated rate constants display similar behaviour above 100 K, showing little or no variation as a function of temperature. Below 100 K, the experimental rate constants increase dramatically, in contrast to the essentially temperature independent theoretical values. The possible origins of the divergence between experiment and theory at low temperatures are discussed.

Quantum Dynamics of Nonadiabatic Renner-Teller Effects in Atom + Diatom Collisions

We review the quantum nonadiabatic dynamics of atom + diatom collisions due to the Renner-Teller (RT) effect, i.e., to the Hamiltonian operators that contain the total spinless electronic angular momentum L̂. As is well-known, this rovibronic effect is large near collinear geometries when at least one of the interacting states is doubly degenerate. In general, this occurs in insertion reactions and at short-range, where the potential wells exhibit deep minima and support metastable complexes. Initial-state-resolved reaction probabilities, integral cross sections, and thermal rate constants are calculated via the real wavepacket method, solving the equation of motion with an approximated or with an exact spinless RT Hamiltonian. We present the dynamics of 10 single-channel or multichannel reactions showing how RT effects depend on the product channels and comparing with the Born-Oppenheimer (BO) approximation or coexisting conical-intersection (CI) interactions. RT effects not only can significantly modify the adiabatic dynamics or correct purely CI results, but also they can be very important in opening collision channels which are closed at the BO or CI level, as in electronic-quenching reactions. In the OH(A²Σ⁺) + Kr electronic quenching, where both nonadiabatic effects (CI and RT) coexist, they are in competition because CI dominates the reactivity but RT couplings reduce the large CI cross section and open a CI-forbidden evolution toward products, so that CI + RT quantum results are in good agreement with experimental or semiclassical findings. The different roles of these couplings are due to the unlike nuclear geometries where they are large: rather far from or near to linearity for CI or RT, respectively. The OH(A²Σ⁺) + Kr electronic quenching was investigated with the exact RT Hamiltonian, validating the approximated one, which was employed for all other collisions.

Dynamics and kinetics of the Si(1D) + H2/D2 reactions on a new global ab initio potential energy surface

Recent studies on the exothermic complex-forming reactions have improved our understanding of complex-forming reactions greatly, however, so far a similar level of study on endothermic ones has been rather limited. In this work, the endothermic complex-forming reaction Si(1D) + H2 → SiH + H and its deuterated isotopic variant are investigated by quantum dynamics (QD) and ring polymer molecular dynamics (RPMD) calculations on a new global ab initio potential energy surface (PES) for the ground electronic state, which is constructed based on 8996 symmetry unique points computed at the icMRCI+Q/aug-cc-pVQZ level. The PES reproduces our ab initio data very well in the dynamically important regions, on which the ro-vibrational energy levels of SiH2 are calculated and general good agreement with experiment is obtained. The integral cross sections and product angular and state distributions are computed in a wide range of collision energies, and interesting dynamics behaviors are revealed. The rate coefficients are also investigated, which display an exponential rise from 2.09 × 10-20 to 6.00 × 10-12 cm3 s-1 for the Si(1D) + H2 reaction as the temperature increases from 300 to 1500 K, in contrast to the nearly temperature-independent behavior of exothermic complex-forming reactions. In addition, the applicability of the RPMD approach is demonstrated, and the kinetic isotope effect is investigated, the ratio of which decreases from 7.89 (300 K) to 1.70 (1500 K). The effects of tunneling and initial rotational excitation are also discussed.

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.

A kinetic study of the N(2D) + C2H4 reaction at low temperature

Electronically excited nitrogen atoms N(²D) are important species in the photochemistry of N₂ based planetary atmospheres such as Titan. Despite this, few N(²D) reactions have been studied over the appropriate low temperature range. During the present work, rate constants were measured for the N(²D) + ethene (C₂H₄) reaction using a supersonic flow reactor at temperatures between 50 K and 296 K. Here, a chemical reaction was used to generate N(²D) atoms, which were detected directly by laser induced fluorescence in the vacuum ultraviolet wavelength region. The measured rate constants displayed very little variation as a function of temperature, with substantially larger values than those obtained in previous work. Indeed, considering an average temperature of 170 K for the atmosphere of Titan leads to a rate constant that is almost seven times larger than the currently recommended value. In parallel, electronic structure calculations were performed to provide insight into the reactive process. While earlier theoretical work at a lower level predicted the presence of a barrier for the N(²D) + C₂H₄ reaction, the present calculations demonstrate that two of the five doublet potential energy surfaces correlating with reagents are likely to be attractive, presenting no barriers for the perpendicular approach of the N atom to the C[double bond, length as m-dash]C bond of ethene. The measured rate constants and new product channels taken from recent dynamical investigations of this process are included in a 1D coupled ion-neutral model of Titan's atmosphere. These simulations indicate that the modeled abundances of numerous nitrogen bearing compounds are noticeably affected by these changes.

Quantum Dynamics Study of the C(1D) + HD Reaction on the ã1A' and b̃1A″ Potential Energy Surfaces

We present an in-depth theoretical study of the C(¹D) + HD (v = 0, j = 0) → CD (CH) (v', j') + H (D) reaction using a time-dependent wave packet method with full Coriolis coupling on the Zhang-Ma-Bian potential energy surfaces (PESs) recently constructed by our group. The integral cross sections (ICS), differential cross sections, CD/CH branching ratios, and product state distributions are calculated over a wide range of collision energies. We find that the vibrational branching ratio defined as ICS(v'=1)/ICS(v'=0) obtained from the b̃¹A″ PES is much smaller than that from the ã¹A' PES for both product channels, which may be attributed to the dynamical effects of the conical intersection regulated (CI-R) intermediate on the b̃¹A″ PES. The collision energy dependence of CD/CH branching ratios displays oscillatory structures, which may be caused by the resonance states supported by the wells on the PESs. The high-temperature rate coefficients are also obtained and compared with previous results. The role of the excited-state PESs is also discussed.

Experimental and theoretical studies of the N(2D) + H2 and D2 reactions

This study reports the first kinetic measurements of the N(²D) + H₂, D₂ reactions below 200 K.

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.

Dynamics of H + HeH+(v = 0, j = 0) → H2+ + He: Insight on the Possible Complex-Forming Behavior of the Reaction

The H + HeH⁺→ He + H₂⁺ reaction has been studied by means of a combination of theoretical approaches: a statistical quantum method (SQM), ring polymer molecular dynamics (RPMD), and the quasiclassical trajectory (QCT) method. Cross sections and rate constants have been calculated in an attempt to investigate the dynamics of the process. The comparison with previous calculations and experimental results reveals that despite the fact that statistical predictions seem to reproduce some of the overall observed features, the analysis at a more detailed state-to-state level shows noticeable deviations from a complex-forming dynamics. We find some differences in cross sections and rate constants obtained in the QCT calculation with a Gaussian binning procedure with respect to previous works in which the standard histogram binning was employed.

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.

Experimental and Theoretical Study of the O(1D) + HD Reaction

This work addresses the kinetics and dynamics of the gas-phase reaction between O(¹D) and HD molecules down to low temperature. Here, measurements were performed by using a supersonic flow (Laval nozzle) reactor coupled with pulsed laser photolysis for O(¹D) production and pulsed-laser-induced fluorescence for O(¹D) detection to obtain rate constants over the 50-300 K range. Additionally, temperature-dependent branching ratios (OD + H/OH + D) were obtained experimentally by comparison of the H/D atom atom yields with those of a reference reaction. In parallel, theoretical rate constants and branching ratios were calculated by using three different techniques; mean potential phase space theory (MPPST), the statistical quantum mechanical method (SQM), and ring polymer molecular dynamics (RPMD). Although the agreement between experimental and theoretical rate constants is reasonably good, with differences not exceeding 30% over the entire temperature range, the theoretical branching ratios derived by the MPPST and SQM methods are as much as 50% larger than the experimental ones. These results are presented in the context of earlier work, while the possible origins of the discrepancies between experiment and theory are discussed.

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.

Low-Temperature Rate Constants and Product-Branching Ratios for the C(1D) + H2O Reaction

The gas-phase reaction between atomic carbon in its first electronically excited ¹D state and water has been studied over the 50-296 K temperature range using a supersonic flow apparatus. C(¹D) atoms were produced by pulsed ultraviolet multiphoton dissociation of carbon tetrabromide; a process that also generates ground-state atomic carbon C(³P). The reaction was followed by detecting product H-atoms by pulsed vacuum ultraviolet laser-induced fluorescence. Two types of experiment were performed. First, temperature-dependent rate constants were derived by recording H-atom formation curves at various gas-phase water concentrations at each temperature. Secondly, temperature-dependent H-atom yields were extracted by comparing the H-atom fluorescence intensities generated by the target C(¹D) + H₂O reaction with those of a reference reaction. The second-order rate constants are large and increase to low temperature, whereas the measured H-atom yields are close to the theoretical maximum value of 2 above 100 K. At 50 K, neither rate constants nor H-atom yields could be derived because of H-atom formation by quantum tunneling in the activated C(³P) + H₂O reaction. The present results are discussed in the context of earlier work on the C(¹D)/C(³P) + H₂O reactions.

Conical intersection-regulated intermediates in bimolecular reactions: Insights from C(1D) + HD dynamics

The importance of conical intersections (CIs) in electronically nonadiabatic processes is well known, but their influence on adiabatic dynamics has been underestimated. Here, through combined experimental and theoretical studies, we show that CIs induce a barrier and regulate conversion from a precursor metastable intermediate (CI-R) to a deep well one. This results in bond-selective activation, influencing the adiabatic dynamics markedly in the C(1D) + HD reaction. Theory is validated by experiment; quantum dynamics calculations on highly accurate ab initio potential energy surfaces yield rate coefficients and product branching ratios in excellent agreement with the experiment. Quasi-classical trajectory calculations reveal that the CI-R intermediate leads to unusual reaction mechanisms (designated as C─H activation complex conversion and cyclic complex), which are responsible for large branching ratios. We also reveal that CI-R intermediates exist in other reactive systems, and the dynamical effects uncovered here may have general significance.

A low temperature investigation of the N(2D) + CH4, C2H6 and C3H8 reactions

Low temperature rate constants for the N(²D) + C₂H₆, C₃H₈ reactions are shown to be much smaller than previously thought.

Rate constants for the N(2D) + C2H2 reaction over the 50–296 K temperature range

The reactions of metastable atomic nitrogen N(²D) are important processes in the gas-phase chemistry of several planetary atmospheres.

First-principles dynamics of collisional intersystem crossing: resonance enhanced quenching of C(1D) by N2

Intersystem crossing is a common and important nonadiabatic process in molecular systems. Revealed with a first-principles investigation, the quenching of C(¹D) by N₂ is efficient due to multiple passages via long-lived collisional resonances, despite relatively small spin–orbit couplings.

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 Ring Polymer Molecular Dynamics Approach to Study the Transition between Statistical and Direct Mechanisms in the H2 + H3+ → H3+ + H2 Reaction

Because of its fundamental importance in astrochemistry, the H₂ + H₃⁺ → H₃⁺ + H₂ reaction has been studied experimentally in a wide temperature range. Theoretical studies of the title reaction significantly lag primarily because of the challenges associated with the proper treatment of the zero-point energy (ZPE). As a result, all previous theoretical estimates for the ratio between a direct proton-hop and indirect exchange (via the H₅⁺ complex) channels deviate from the experiment, in particular, at lower temperatures where the quantum effects dominate. In this work, the ring polymer molecular dynamics (RPMD) method is applied to study this reaction, providing very good agreement with the experiment. RPMD is immune to the shortcomings associated with the ZPE leakage and is able to describe the transition from direct to indirect mechanisms below room temperature. We argue that RPMD represents a useful tool for further studies of numerous ZPE-sensitive chemical reactions that are of high interest in astrochemistry.

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.

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.

A low temperature investigation of the gas-phase N(2D) + NO reaction. Towards a viable source of N(2D) atoms for kinetic studies in astrochemistry

The gas-phase reaction of metastable atomic nitrogen N(²D) with nitric oxide has been investigated over the 296–50 K temperature range using a supersonic flow reactor.

Dynamical importance of van der Waals saddle and excited potential surface in C(1D)+D2 complex-forming reaction

Encouraged by recent advances in revealing significant effects of van der Waals wells on reaction dynamics, many people assume that van der Waals wells are inevitable in chemical reactions. Here we find that the weak long-range forces cause van der Waals saddles in the prototypical C(¹D)+D₂ complex-forming reaction that have very different dynamical effects from van der Waals wells at low collision energies. Accurate quantum dynamics calculations on our highly accurate ab initio potential energy surfaces with van der Waals saddles yield cross-sections in close agreement with crossed-beam experiments, whereas the same calculations on an earlier surface with van der Waals wells produce much smaller cross-sections at low energies. Further trajectory calculations reveal that the van der Waals saddle leads to a torsion then sideways insertion reaction mechanism, whereas the well suppresses reactivity. Quantum diffraction oscillations and sharp resonances are also predicted based on our ground- and excited-state potential energy surfaces.

It is commonly held that van der Waals wells are inevitable in chemical reactions. Here, the authors show that weak van der Waals forces in the entrance channel of a prototypical complex-forming reaction cause a van der Waals saddle instead, with different dynamical effects from a well at low collision energies.

Quasiclassical trajectory study of the C(1D) + HD reaction

Isotopic branching ratios are investigated by detailed quasiclassical trajectory calculations on our recent singlet ground and excited potential energy surfaces.

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