Rotational and Isotopic Effects in the H2 + OH+ → H + H2O+ Reaction

Higher-Order Split Operator Schemes for Solving Tetratomic Reactions Using the Time-Dependent Wave Packet Method

In this work, using the time-dependent quantum wave packet method, quite a few typical higher-order split operators (HOSOs) were for the first time applied to calculate the tetratomic reactive scattering processes in the hyperspherical coordinate. It was found that the HOSOs were hardly efficient for a tetratomic reaction calculation, unlike those for a triatomic reactive scattering calculation. We proposed an efficient HOSO with a force gradient (denoted as 2G1 in the main text) for efficiently and accurately calculating a tetratomic reaction using the quantum wave packet method. Several typical tetratomic reactions, such as H2 + OH, HF + OH, and H2 + OH+, are calculated for demonstrating the effectiveness of the proposed 2G1 in terms of (product state-resolved) reaction probability and inelastic probability, by comparing with the performance of the previously reported various HOSOs. We suggest that the 2G1 propagator could be applied to efficiently calculate a general tetratomic reaction.

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.

Theoretical studies on the kinetics and dynamics of the BeH+ + H2O reaction: comparison with the experiment

The reaction of BeH+ with background gaseous H2O may play a role in qubit loss for quantum information processing with Be+ as trapped ions, and yet its reaction mechanism has not been well understood until now. In this work, a globally accurate, full-dimensional ground-state potential energy surface (PES) for the BeH+ + H2O reaction was constructed by fitting a total of 170 438 ab initio energy points at the level of RCCSD(T)-F12/aug-cc-pVTZ using the fundamental invariant-neural network method. The total root-mean-square error of the final PES was 0.178 kcal mol-1. For comparison, quasi-classical trajectory calculations were carried out on the PES at an experimental temperature of 150 K. The obtained thermal rate constant and product branching ratio of the BeD+ + H2O reaction agreed quite well with experimental results. In addition, the vibrational state distributions and energy disposals of the products were calculated and rationalized using the sudden vector projection model.

Validating experiments for the reaction H2 + NH2- by dynamical calculations on an accurate full-dimensional potential energy surface

Ion-molecule reactions play key roles in the field of ion related chemistry. As a prototypical multi-channel ion-molecule reaction, the reaction H₂ + NH₂^- → NH₃ + H^- has been studied for decades. In this work, we develop a new globally accurate potential energy surface (PES) for the title system based on hundreds of thousands of sampled points over a wide dynamically relevant region that covers long-range interacting configuration space. The permutational invariant polynomial-neural network (PIP-NN) method is used for fitting and the resulting total root mean squared error (RMSE) is extremely small, 0.026 kcal mol^-1. Extensive dynamical and kinetic calculations are carried out on this PIP-NN PES. Impressively, a unique phenomenon of significant reactivity suppression by exciting the rotational mode of H₂ is reported, supported by both the quasi-classical trajectory (QCT) and quantum dynamics (QD) calculations. Further analysis uncovers that exciting the H₂ rotational mode would prevent the formation of the reactant complex and thus suppress the reactivity. The calculated rate coefficients for H₂/D₂ + NH₂^- agree well with the experimental results, which show an inverse temperature dependence from 50 to 300 K, consistent with the capture nature of this barrierless reaction. The significant kinetic isotope effect observed by experiments is well reproduced by the QCT computations as well.

Mode-specific dynamics in multichannel reaction NH+ + H2

The vibrational- and rotational-mode specificity of the multichannel NH⁺ + H₂ reaction is studied on a recently constructed ab initio-based global potential energy surface using an initial state selected quasi-classical trajectory method, and the trajectories are analyzed using an isometric feature mapping and k-means approach. All excitation modes promote two reactions (R1: NH'⁺ + H₂ → NH⁺ + HH' and R4: NH'⁺ + H₂ → NH₂⁺ + H') where both NH and HH bonds are broken, but reduce the reactivity of the proton-transfer reaction R2 (NH'⁺ + H₂ → N + H'H₂⁺) at low collision energies. For the hydrogen-transfer reaction R3 (NH'⁺ + H₂ → HNH'⁺ + H), the rotational excitation of NH⁺ enhances the reactivity remarkably, while its vibrational excitation has an inhibiting effect on the reaction. The trajectory analyses show that the vibrational and rotational excitations of NH⁺ make R3 tend to go over a submerged saddle point instead of extracting hydrogen atoms directly. On the other hand, the motions of the H₂ reactant facilitate the enhancement of the reactivity but they do not affect the mechanism of R3. In addition, the results suggest that the coupling of the isometric feature mapping and the k-means approach in the trajectory analysis is an appropriate tool for reaction-dynamics studies.

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.

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.

High-resolution collision energy control through ion position modulation in atom-ion hybrid systems

We demonstrate an ion shuttling technique for high-resolution control of atom-ion collision energy by translating an ion held within a radio-frequency trap through a magneto-optical atom trap. The technique is demonstrated both experimentally and through numerical simulations, with the experimental results indicating control of ion kinetic energies from 0.05 to 1 K with a fractional resolution of ∼10 and the simulations demonstrating that kinetic energy control up to 120 K with a maximum predicted resolution of ∼100 is possible, offering order-of-magnitude improvements over most alternative techniques. Finally, we perform a proof-of-principle chemistry experiment using this technique and outline how the method may be refined in the future and applied to the study of molecular ion chemistry.

Low temperature rates for key steps of interstellar gas-phase water formation

The gas-phase formation of water molecules in the diffuse interstellar medium (ISM) proceeds mainly via a series of reactions involving the molecular ions OH⁺, H₂O⁺, and H₃O⁺ and molecular hydrogen. These reactions form the backbone for the chemistry leading to the formation of several complex molecular species in space. A comprehensive understanding of the mechanisms involved in these reactions in the ISM necessitates an accurate knowledge of the rate coefficients at the relevant temperatures (10 to 100 K). We present measurements of the rate coefficients for two key reactions below 100 K, which, in both cases, are significantly higher than the values used in astronomical models thus far. The experimental rate coefficients show excellent agreement with dedicated theoretical calculations using a novel ring-polymer molecular dynamics approach that offers a first-principles treatment of low-temperature barrierless gas-phase reactions, which are prevalent in interstellar chemical networks.

Comparison of experimental and theoretical quantum-state-selected integral cross-sections for the H2O+ + H2 (D2) reactions in the collision energy range of 0.04–10.00 eV

A combined experimental–theoretical study of the rovibrationally state-selected ion–molecule reactions H₂O⁺(X²B₁; v₁⁺v₂⁺v₃⁺; N_Ka+Kc+⁺) + H₂ (D₂) → H₃O⁺ (H₂DO⁺) + H (D), where (v₁⁺v₂⁺v₃⁺) = (000), (020), and (100) and N_Ka+Kc+⁺ = 0₀₀, 1₁₁, and 2₁₁.

Mode specific dynamics in the H2 + SH → H + H2S reaction

Full-dimensional quantum dynamics and quasi-classical trajectory studies indicate strong mode selectivity in the H₂ + SH reaction.