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

Theoretical Study of the Thermal Rate Coefficients of the H3+ + C2H4 Reaction: Dynamics Study on a Full-Dimensional Potential Energy Surface

Ion-molecular reactions play a significant role in molecular evolution within the interstellar medium. In this study, the entrance channel reaction, H3+ + C2H4 → H2 + C2H5+, was investigated using classical molecular dynamic (classical MD) and ring polymer molecular dynamic (RPMD) simulation techniques. We developed an analytical potential energy surface function with a permutationally invariant polynomial basis, specifically employing the monomial symmetrized approach. Our dynamic simulations reproduced the rate coefficient of 300 K for H3+ + C2H4 → H2 + C2H5+, aligning reasonably well with the values in the kinetic database commonly utilized in astrochemistry. The thermal rate coefficients obtained using both the classical MD and RPMD techniques exhibited an increase from 100 K to 300 K as the temperature rose. Additionally, we analyzed the excess energy distribution of the C2H5+ fragment with respect to temperature to investigate the indirect reaction pathway of C2H5+ → H2 + C2H3+. This result suggests that the indirect reaction pathway of C2H5+ → H2 + C2H3+ holds minor significance, although the distribution highly depends on the collisional temperature.

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.

An efficient algorithm for capturing quantum effects in classical reactive scattering: application to D + H+3 → H2D+ + H

Motivated by a recent semiclassical analysis of chemical reaction thresholds [Bonnet et al., J. Chem. Phys., 2022, 157, 094114], we present an efficient algorithm for including zero-point energy (ZPE) effects in classical reactive scattering. The algorithm is an extension of the quasi-classical trajectory (QCT) Gaussian binning method. We apply it to the astrophysically important D + H+3 reaction, where there are significant quantum effects and where application of other methods is problematic [Braunstein et al., Phys. Chem. Chem. Phys., 2022, 24, 5489]. The rate constants computed with the new, general algorithm closely match recent Ring Polymer Molecular Dynamics (RPMD) [Bulut et al., J. Phys. Chem. A, 2019, 123, 8766] and experimentally derived [Bowen et al., J. Chem. Phys., 2021, 154, 084307] ones spanning ∼4 orders of magnitude from 70 to 1500 K.

Low-temperature reaction dynamics of paramagnetic species in the gas phase

Radicals are abundant in a range of important gas-phase environments. They are prevalent in the atmosphere, in interstellar space, and in combustion processes. As such, understanding how radicals react is essential for the development of accurate models of the complex chemistry occurring in these gas-phase environments. By controlling the properties of the colliding reactants, we can also gain insights into how radical reactions occur on a fundamental level. Recent years have seen remarkable advances in the breadth of experimental methods successfully applied to the study of reaction dynamics involving paramagnetic species-from improvements to the well-known crossed molecular beams approach to newer techniques involving magnetically guided and decelerated beams. Coupled with ever-improving theoretical methods, quantum features are being observed and interesting insights into reaction dynamics are being uncovered in an increasingly diverse range of systems. In this highlight article, we explore some of the exciting recent developments in the study of chemical dynamics involving paramagnetic species. We focus on low-energy reactive collisions involving neutral radical species, where the reaction parameters are controlled. We conclude by identifying some of the limitations of current methods and exploring possible new directions for the field.

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.

Kinetic and dynamic studies of the NH2+ + H2 reaction on a high-level ab initio potential energy surface

The dynamics and kinetics of the NH2+ + H2 reaction are investigated on a newly developed ab initio potential energy surface using the quasi-classical trajectory method.

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.

Study on the kinetics and dynamics of the H2 + NH2- reaction on a high-level ab initio potential energy surface

Gas-phase ion-molecule reactions play major roles in many fields of chemistry and physics. The reaction of an amino radical anion with a hydrogen molecule is one of the simplest proton transfer reactions involving anions. A globally accurate full-dimensional potential energy surface (PES) for the NH₂^- + H₂ reaction is developed by the fundamental invariant-neural network method, resulting in a root mean square error of 0.116 kcal mol^-1. Quasi-classical trajectory calculations are then carried out on the newly developed PES to give integral cross sections, differential cross sections and thermal rate coefficients. This reaction has two reaction channels, proton transfer and hydrogen exchange. The reactivity of the proton transfer channel is about one or two orders of magnitude stronger than that of the hydrogen exchange channel in the energy range studied. Vibrational excitation of H₂ promotes the proton transfer reaction, while fundamental excitation of each vibrational mode of NH₂^- has a negligible effect. In addition, the theoretical rate coefficients of the proton transfer reaction on the PES show inverse temperature dependence from 150 to 750 K, in accordance with the available experimental results.

The symmetric C-D stretching spectator mode in the H + CHD3 → H2 + CD3 reaction and its effect on dynamical modeling

The symmetric C-D stretching mode is a spectator mode in the H + CHD3 → H2 + CD3 reaction. Effects of multiple vibrational excitations of the CHD3 reactant are studied with the quantum transition-state (QTS) framework and an eight-dimensional (8D) model Hamiltonian developed by Palma and Clary. By including many thermal flux eigenstates, results have been obtained up to high energies, allowing the study of the symmetric C-D stretching spectator mode. A new concept of a state-specific thermal flux operator is proposed to analyze the C-D stretching spectator mode in detail, providing a new and insightful venue for studying transition-state control of chemical reactions. Furthermore, as a spectator mode, whether the C-D stretching motion can be excluded in a seven-dimensional (7D) model has not been fully interrogated, although the 7D model is a reasonable approximation and has provided accurate theoretical predictions. By comparing with available results of full-dimensional calculations, both the 7D and 8D models predict reasonably accurate results. However, the 7D model underestimates the mixing of two vibrational states that are in Fermi resonance. Despite its spectator nature, the C-D stretch is important in the dynamical modeling of chemical reaction systems affected by state mixing.

Large Anharmonic Effects on Tunneling and Kinetics: Reaction of Propane with Muonium

Calculations of kinetic isotope effects (KIEs) provide challenging tests of quantal mass effects on reaction rates, and muonium KIEs are the most challenging. Here, we show that it can be very important to include reaction-coordinate-dependent vibrational anharmonicity along the whole reaction path to calculate tunneling probabilities and KIEs. For the reaction of propane with Mu, this decreases both the height and width of the vibrationally adiabatic potential barrier, with both effects increasing the rate constants. Our results agree well with the experimental observations.

Neural network potential energy surface for the low temperature ring polymer molecular dynamics of the H2CO + OH reaction

A new potential energy surface (PES) and dynamical study of the reactive process of H₂CO + OH toward the formation of HCO + H₂O and HCOOH + H are presented. In this work, a source of spurious long range interactions in symmetry adapted neural network (NN) schemes is identified, which prevents their direct application for low temperature dynamical studies. For this reason, a partition of the PES into a diabatic matrix plus a NN many-body term has been used, fitted with a novel artificial neural network scheme that prevents spurious asymptotic interactions. Quasi-classical trajectory (QCT) and ring polymer molecular dynamics (RPMD) studies have been carried on this PES to evaluate the rate constant temperature dependence for the different reactive processes, showing good agreement with the available experimental data. Of special interest is the analysis of the previously identified trapping mechanism in the RPMD study, which can be attributed to spurious resonances associated with excitations of the normal modes of the ring polymer.

Dynamics of the isotope exchange reaction of D with H3 +, H2D+, and D2H. J Chem Phys 2021;154:084307. [PMID: 33639774 DOI: 10.1063/5.0038434]

We have measured the merged-beams rate coefficient for the titular isotope exchange reactions as a function of the relative collision energy in the range of ∼3 meV-10 eV. The results appear to scale with the number of available sites for deuteration. We have performed extensive theoretical calculations to characterize the zero-point energy corrected reaction path. Vibrationally adiabatic minimum energy paths were obtained using a combination of unrestricted quadratic configuration interaction of single and double excitations and internally contracted multireference configuration interaction calculations. The resulting barrier height, ranging from 68 meV to 89 meV, together with the various asymptotes that may be reached in the collision, was used in a classical over-the-barrier model. All competing endoergic reaction channels were taken into account using a flux reduction factor. This model reproduces all three experimental sets quite satisfactorily. In order to generate thermal rate coefficients down to 10 K, the internal excitation energy distribution of each H₃ ⁺ isotopologue is evaluated level by level using available line lists and accurate spectroscopic parameters. Tunneling is accounted for by a direct inclusion of the exact quantum tunneling probability in the evaluation of the cross section. We derive a thermal rate coefficient of <1×10^-12 cm³ s^-1 for temperatures below 44 K, 86 K, and 139 K for the reaction of D with H₃ ⁺, H₂D⁺, and D₂H⁺, respectively, with tunneling effects included. The derived thermal rate coefficients exceed the ring polymer molecular dynamics prediction of Bulut et al. [J. Phys. Chem. A 123, 8766 (2019)] at all temperatures.