Calculations on the rate of the ion-molecule reaction between NH3+ and H2

Tunnelling measured in a very slow ion-molecule reaction

Quantum tunnelling reactions play an important role in chemistry when classical pathways are energetically forbidden¹, be it in gas-phase reactions, surface diffusion or liquid-phase chemistry. In general, such tunnelling reactions are challenging to calculate theoretically, given the high dimensionality of the quantum dynamics, and also very difficult to identify experimentally^2-4. Hydrogenic systems, however, allow for accurate first-principles calculations. In this way the rate of the gas-phase proton-transfer tunnelling reaction of hydrogen molecules with deuterium anions, H₂ + D^- → H^- + HD, has been calculated⁵, but has so far lacked experimental verification. Here we present high-sensitivity measurements of the reaction rate carried out in a cryogenic 22-pole ion trap. We observe an extremely low rate constant of (5.2 ± 1.6) × 10^-20 cm³ s^-1. This measured value agrees with quantum tunnelling calculations, serving as a benchmark for molecular theory and advancing the understanding of fundamental collision processes. A deviation of the reaction rate from linear scaling, which is observed at high H₂ densities, can be traced back to previously unobserved heating dynamics in radiofrequency ion traps.

Rotational Cooling Effect on the Rate Constant in the CH3F + Ca+ Reaction at Low Collision Energies

The rotational cooling effect on the reaction rate constant of the gas-phase ion-polar-molecule reaction CH₃F + Ca⁺ → CH₃ + CaF⁺ was experimentally studied at low collision energies. Fluoromethane molecules showed higher reactivity as the rotational temperature decreased. The experimental rate constants were compared with the capture rate constants which were obtained by the Perturbed Rotational State (PRS) theory assuming the rotational level distribution corresponding to the experimental conditions. The PRS result shows a strong dependence of the capture rate constants on the rotational level distribution in accordance with the experimental findings. However, the PRS capture rate constants deviate from the measurement values as the average collision energy increases especially when the fluoromethane molecules are rotationally cooled far below room temperature. The present paper suggests that the rotational state distribution significantly affects the rate constants of ion-polar-molecule reactions and is one of the important issues to be considered in the study of molecular synthesis in the interstellar medium, where the thermal equilibrium is not necessarily established.

A study of the translational temperature dependence of the reaction rate constant between CH3CN and Ne+ at low temperatures

We have measured the translational temperature dependence of the reaction rate constant for CH₃CN + Ne⁺ → products at low temperatures. A cold Ne⁺ ensemble was embedded in Ca⁺ Coulomb crystals by a sympathetic laser cooling technique, while cold acetonitrile (CH₃CN) molecules were produced by two types of Stark velocity filters to widely change the translational temperatures. The measured reaction rate constant gradually increases with the decrease in the translational temperature of the velocity-selected CH₃CN molecules from 60 K down to 2 K, and thereby, a steep increase was observed at temperatures lower than 5 K. A comparison between experimental rate constants and the ion-dipole capture rate constants by the Perturbed Rotational State (PRS) theory was performed. The PRS capture rate constant reproduces well the reaction rate constant at a few kelvin but not for temperatures higher than 5 K. The result indicates that the reaction probability is small compared to typical ion-polar molecule reactions at temperatures above 5 K.

Atom tunnelling in the reaction NH3+ + H2 → NH4+ + H and its astrochemical relevance

The title reaction is involved in the formation of ammonia in the interstellar medium. We have calculated thermal rates including atom tunnelling using different rate theories. Canonical variational theory with microcanonically optimised multidimensional tunnelling was used for bimolecular rates, modelling the gas-phase reaction and also a surface-catalysed reaction of the Eley-Rideal type. Instanton theory provided unimolecular rates, which model the Langmuir-Hinshelwood type surface reaction. The potential energy was calculated on the CCSD(T)-F12 level of theory on the fly. We report thermal rates and H/D kinetic isotope effects. The latter have implications for observed H/D fractionation in molecular clouds. Tunnelling causes rate constants to be sufficient for the reaction to play a role in interstellar chemistry even at cryogenic temperature. We also discuss intricacies and limitations of the different tunnelling approximations to treat this reaction, including its pre-reactive minimum.

Pressure dependent low temperature kinetics for CN + CH3CN: competition between chemical reaction and van der Waals complex formation

The gas phase reaction between the CN radical and acetonitrile CH₃CN was investigated experimentally with a CRESU apparatus and a slow flow reactor as well as theoretically to explore the temperature and pressure dependence of its rate coefficient from 354 K down to 23 K.

Concluding remarks: astrochemistry of dust, ice and gas

In this closing article, we first introduce the topics of dust and ice chemistry and their role in astrochemistry. We then discuss the invited contributions and discussions concerning these topics, dividing the papers into groupings by subject: (i) astronomical sources, (ii) basic properties of dust, (iii) processes on bare grains, (iv) processes on and in ice mantles, and (v) complex organic molecules. A sample of poster contributions is included in the text, when they complement the discussion. The article ends with some suggestions for future research.

Why HOC+ is detectable in interstellar clouds: the rate of the reaction between HOC+ and H2

The recent confirmation by Ziurys and Apponi of the detection of HOC+ toward Sgr B2 (OH), and their identification of the ion in Orion-KL and several other sources show that HOC+ is far more abundant than predicted by previous ion-molecule models. In these models, the reaction HOC(+) + H2 --> HCO(+) + H2 is assumed to rapidly destroy HOC+, based on the results of a prior calculation. We have recalculated the rate of this reaction as a function of temperature using a new ab initio potential surface and a phase space approach to the dynamics which includes tunneling. The newly calculated rate is small (< or = 1 x 10(-10) cm3 s-1) at temperatures under 100 K.

Calculations on the competition between association and reaction for C3H(+) + H2

A potential energy surface has been calculated for the competing associative and reactive ion-molecule processes involving the reactants C3H(+) + H2. Our ab initio results show that the linear ion C3H+ and H2 can directly access the deep potential well of the propargyl ion H2CCCH+, which is calculated to lie 390 kJ mol-1 below the zero-point energy of the reactants. Isomerization between the propargyl ion and the lower energy, cyclic C3H3+ ion, calculated to lie 501 kJ mol-1 below the zero-point energy of reactants, can subsequently occur via two pathways. One of these pathways involves a transition state lying 22 kJ mol-1 below the energy of the reactants while the other, which occurs at much lower energies, involves two transition states and an intermediate. The dissociation of c-C3H3+ into c-C3H2(+) + H is calculated to occur directly, without any intermediate potential energy maximum, but the energy of the products lies 7.3 kJ mol-1 above the energy of the reactants. Using the minimum energy potential pathway and properties of the stationary point structures determined via ab initio methods, we have calculated both the association rate coefficient to produce C3H3+ as a function of density and the branching ratio between the propargyl and cyclic structures of the ion. Our results are in good agreement with some experimental results and in conflict with others. Specifically, we agree with the 1:1 branching ratio measured for the propargyl and cyclic isomers of C3H3+ at 80 and 300 K and we agree with the rate coefficient for radiative association measured at 80 K. We cannot reproduce reported measurements that the reactive channel (C3H2(+) + H) is the dominant channel at 80 K and at low gas densities, or that the association channel at high densities saturates at an effective rate coefficient well below the Langevin value -2x10(-11) cm3 s-1 at 300K and 1x10(-10) cm3 s-1 at 80K.