Ion-molecule chemistry at very low temperatures: cold chemical reactions between Coulomb-crystallized ions and velocity-selected neutral molecules

Cold Ion-Molecule Reactions in the Extreme Environment of a Coulomb Crystal

Coulomb crystals provide a unique environment in which to study ion-neutral gas-phase reactions. In these cold, trapped ensembles, we are able to study the kinetics and dynamics of small molecular systems. These measurements have connections to chemistry in the Interstellar Medium (ISM) and planetary atmospheres. This Feature Article will describe recent work in our laboratory that uses Coulomb crystals to study translationally cold, ion-neutral reactions. We provide a description of how the various affordances of our experimental system allow for detailed studies of the reaction mechanisms and the corresponding products. In particular, we will describe quantum-state resolved reactions, isomer-dependent reactions, and reactions with a rarely studied, astrophysically relevant ion, CCl+.

Reactions of Acetonitrile with Trapped, Translationally Cold Acetylene Cations

The reaction of the acetylene cation (C₂H₂⁺) with acetonitrile (CH₃CN) is measured in a linear Paul ion trap coupled to a time-of-flight mass spectrometer. C₂H₂⁺ and CH₃CN are both noted for their astrochemical abundance and predicted relevance for understanding prebiotic chemistry. The observed primary products are c-C₃H₃⁺, C₃H₄⁺, and C₂NH₃⁺. The latter two products react with excess CH₃CN to form the secondary product C₂NH₄⁺, protonated acetonitrile. The molecular formula of these ionic products can be verified with the aid of isotope substitution via deuteration of the reactants. Primary product reaction pathways and thermodynamics are investigated with quantum chemical calculations and demonstrate exothermic pathways to two isomers of C₂NH₃⁺, two isomers of C₃H₄⁺, and the cyclopropenyl cation c-C₃H₃⁺. This study deepens our understanding of the dynamics and products of a pertinent ion-molecule reaction between two astrochemically abundant molecules in conditions that mimic those of the interstellar medium.

Reaction of an Ion and a Free Radical near 0 K: He+ + NO → He + N+ + O

The reactions between ions and free radicals are among the fastest chemical reactions. They are predicted to proceed with large rates, even near 0 K, but so far, this prediction has not been verified experimentally. We report on measurements of the rate coefficient of the reaction between the ion He⁺ and the free radical NO at collision energies in the range between 0 and ∼ k_B·10 K. To avoid heating of the ions by stray electric fields, the reaction is observed within the large orbit of a Rydberg electron of principal quantum number n ≥ 30, which shields the ion from external electric fields without affecting the reaction. Low collision energies are reached by merging a supersonic beam of He Rydberg atoms with a supersonic beam of NO molecules and adjusting their relative velocity using a chip-based Rydberg-Stark decelerator and deflector. We observe a strong enhancement of the reaction rate at collision energies below ∼k_B·2 K. This enhancement is interpreted on the basis of adiabatic-channel capture-rate calculations as arising from the near-degenerate rotational levels of opposite parity resulting from the Λ-doubling in the X ²Π_1/2 ground state of NO. With these new results, we examine the reliability of broadly used approximate analytic expressions for the thermal rate constants of ion-molecule reactions at low temperatures.

Ion-molecule reactions in the HBr+ + HCl (DCl) system: a combined experimental and theoretical study

Reactions in the system HBr⁺ + HCl (DCl) were investigated inside a guided ion-beam apparatus under single-collision conditions. In the HBr⁺ + HCl system, the proton transfer (PT_HCl) and charge transfer (CT) are observable. In the HBr⁺ + DCl system, proton transfer (PT_DCl) and deuterium abstraction (DA) are accessible. The cross sections for all reaction channels were measured as a function of the collision energy E_cm and of the rotational energy E_rot of the ion. The rotationally state-selective formation of the ionic species was realized by resonance-enhanced multiphoton ionization (REMPI). As expected, the PT-channels are exothermic, and the cross section decreases with increasing collision energy for both PT_HCl and PT_DCl. The cross section for DA also decreases with an increasing E_c.m.. In the case of a considerably endothermic CT-channel, the reaction efficiency increases with increasing collision energy but has an overall much smaller cross sections compared to PT and DA reactions. Both PT-reactions are hindered by ion rotation, whereas DA is independent of E_rot. The CT-channel shows a rotational enhancement near the thermochemical threshold. The experiment is complemented by theory, using ab initio molecular dynamics (AIMD, also known as direct dynamics) simulations and taking the rotational enhancement of HBr⁺ into account. The simulations show good agreement with the experimental results. The cross section of PT_HCl decreases with an increase of the rotational energy. Furthermore, the absolute cross sections are in the same order of magnitude. The CT channel shows no reactions in the simulation.

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.

Capture Theory Models: An overview of their development, experimental verification, and applications to ion-molecule reactions

Since Arrhenius first proposed an equation to account for the behaviour of thermally activated reactions in 1889, significant progress has been made in our understanding of chemical reactivity. A number of capture theory models have been developed over the past several decades to predict the rate coefficients for reactions between ions and molecules-ranging from the Langevin equation (for reactions between ions and non-polar molecules) to more recent fully quantum theories (for reactions at ultra-cold temperatures). A number of different capture theory methods are discussed, with the key assumptions underpinning each approach clearly set out. The strengths and limitations of these capture theory methods are examined through detailed comparisons between low-temperature experimental measurements and capture theory predictions. Guidance is provided on the selection of an appropriate capture theory method for a given class of ion-molecule reaction and set of experimental conditions-identifying when a capture-based model is likely to provide an accurate prediction. Finally, the impact of capture theories on fields such as astrochemical modelling is noted, with some potential future directions of capture-based approaches outlined.

Multipole-moment effects in ion-molecule reactions at low temperatures: part II - charge-quadrupole-interaction-induced suppression of the He+ + N2 reaction at collision energies below kB·10 K

We report on an experimental and theoretical investigation of the He+ + N2 reaction at collision energies in the range between 0 and kB·10 K. The reaction is studied within the orbit of a highly excited Rydberg electron after merging a beam of He Rydberg atoms (He(n), n is the principal quantum number), with a supersonic beam of ground-state N2 molecules using a surface-electrode Rydberg-Stark decelerator and deflector. The collision energy Ecoll is varied by changing the velocity of the He(n) atoms for a fixed velocity of the N2 beam and the relative yields of the ionic reaction products N+ and N2+ are monitored in a time-of-flight mass spectrometer. We observe a reduction of the total reaction-product yield of ∼30% as Ecoll is reduced from ≈kB·10 K to zero. An adiabatic capture model is used to calculate the rotational-state-dependent interaction potentials experienced by the N2 molecules in the electric field of the He+ ion and the corresponding collision-energy-dependent capture rate coefficients. The total collision-energy-dependent capture rate coefficient is then determined by summing over the contributions of the N2 rotational states populated at the 7.0 K rotational temperature of the supersonic beam. The measured and calculated rate coefficients are in good agreement, which enables us to attribute the observed reduction of the reaction rate at low collision energies to the negative quadrupole moment, Qzz, of the N2 molecules. The effect of the sign of the quadrupole moment is illustrated by calculations of the rotational-state-dependent capture rate coefficients for ion-molecule reactions involving N2 (negative Qzz value) and H2 (positive Qzz value) for |J, M〉 rotational states with J ≤ 5 (M is the quantum number associated with the projection of the rotational angular momentum vector J⃑ on the collision axis). With decreasing value of |M|, J⃑ gradually aligns perpendicularly to the collision axis, leading to increasingly repulsive (attractive) interaction potentials for diatomic molecules with positive (negative) Qzz values and to reaction rate coefficients that decrease (increase) with decreasing collision energies.

Velocity-tunable beam of continuously decelerated polar molecules for cold ion-molecule reaction studies

Producing high densities of molecules is a fundamental challenge for low-temperature, ion-molecule reaction studies. Traveling-wave Stark decelerators promise to deliver high density beams of cold, polar molecules but require non-trivial control of high-voltage potentials. We have overcome this experimental challenge and demonstrate continuous deceleration of ND₃ from 385 to 10 m/s, while driving the decelerator electrodes with a 10 kV amplitude sinewave. In addition, we test an alternative slowing scheme, which increases the time delay between decelerated packets of ND₃ and non-decelerated molecules, allowing for better energy resolution of subsequent reaction studies. We characterize this source of neutral, polar molecules suitable for energy-resolved reaction studies with trapped ions at cold translational temperatures. We also propose a combined apparatus consisting of the traveling-wave decelerator and a linear ion trap with a time-of-flight mass spectrometer and discuss to what extent it may achieve cold, energy-resolved, ion-neutral reactions.

Inverse kinetic isotope effects in the charge transfer reactions of ammonia with rare gas ions

In the absence of experimental data, models of complex chemical environments rely on predicted reaction properties. Astrochemistry models, for example, typically adopt variants of capture theory to estimate the reactivity of ionic species present in interstellar environments. In this work, we examine astrochemically-relevant charge transfer reactions between two isotopologues of ammonia, NH3 and ND3, and two rare gas ions, Kr+ and Ar+. An inverse kinetic isotope effect is observed; ND3 reacts faster than NH3. Combining these results with findings from an earlier study on Xe+ (Petralia et al., Nat. Commun., 2020, 11, 1), we note that the magnitude of the kinetic isotope effect shows a dependence on the identity of the rare gas ion. Capture theory models consistently overestimate the reaction rate coefficients and cannot account for the observed inverse kinetic isotope effects. In all three cases, the reactant and product potential energy surfaces, constructed from high-level ab initio calculations, do not exhibit any energetically-accessible crossing points. Aided by a one-dimensional quantum-mechanical model, we propose a possible explanation for the presence of inverse kinetic isotope effects in these charge transfer reaction systems.

Towards chemistry at absolute zero

The prospect of cooling matter down to temperatures that are close to absolute zero raises intriguing questions about how chemical reactivity changes under these extreme conditions. Although some types of chemical reaction still occur at 1 μK, they can no longer adhere to the conventional picture of reactants passing over an activation energy barrier to become products. Indeed, at ultracold temperatures, the system enters a fully quantum regime, and quantum mechanics replaces the classical picture of colliding particles. In this Review, we discuss recent experimental and theoretical developments that allow us to explore chemical reactions at temperatures that range from 100 K to 500 nK. Although the field is still in its infancy, exceptional control has already been demonstrated over reactivity at low temperatures.

Self-Reactions in the HBr+ (DBr+) + HBr System: A State-Selective Investigation of the Role of Rotation

Self-reactions observed in the HBr⁺ (DBr⁺) + HBr system have been investigated using a guided ion-beam experiment under single-collision conditions. The reaction channels observed are proton transfer/hydrogen abstraction (PT/HA) in the case of HBr⁺ and deuteron transfer/hydrogen abstraction (DT/HA) and charge transfer (CT) in the case of DBr⁺. HBr⁺/DBr⁺ ions have been formed with rotational energies selected using the resonance-enhanced multiphoton ionization (REMPI) formation process. Cross sections have been measured as a function of the rotational energy of the ion, E_rot, and of the center-of-mass collision energy, E_cm. In the region of low rotational energies, the cross section for both PT/HA and DT/HA decreases with increasing ion rotation. In this region, the cross section for CT increases with increasing ion rotation. For higher rotational energies, the cross section increases with increasing ion rotation for PT/HA and less pronounced for DT/HA. The cross section for CT becomes independent of ion rotation for high rotational energies. Since all reaction channels are exothermic, all cross sections decrease with increasing E_cm.

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.

Strong inverse kinetic isotope effect observed in ammonia charge exchange reactions

Isotopic substitution has long been used to understand the detailed mechanisms of chemical reactions; normally the substitution of hydrogen by deuterium leads to a slower reaction. Here, we report our findings on the charge transfer collisions of cold [Formula: see text] ions and two isotopologues of ammonia, [Formula: see text] and [Formula: see text]. Deuterated ammonia is found to react more than three times faster than hydrogenated ammonia. Classical capture models are unable to account for this pronounced inverse kinetic isotope effect. Moreover, detailed ab initio calculations cannot identify any (energetically accessible) crossing points between the reactant and product potential energy surfaces, indicating that electron transfer is likely to be slow. The higher reactivity of [Formula: see text] is attributed to the greater density of states (and therefore lifetime) of the deuterated reaction complex compared to the hydrogenated system. Our observations could provide valuable insight into possible mechanisms contributing to deuterium fractionation in the interstellar medium.

Conformer Selection by Matter-Wave Interference

We establish that matter-wave diffraction at near-resonant ultraviolet optical gratings can be used to spatially separate individual conformers of complex molecules. Our calculations show that the conformational purity of the prepared beam can be close to 100% and that all molecules remain in their electronic ground state. The proposed technique is independent of the dipole moment and the spin of the molecule and thus paves the way for structure-sensitive experiments with hydrocarbons and biomolecules, such as neurotransmitters and hormones, which have evaded conformer-pure isolation so far.

An ion trap time-of-flight mass spectrometer with high mass resolution for cold trapped ion experiments

Trapping molecular ions that have been sympathetically cooled with laser-cooled atomic ions is a useful platform for exploring cold ion chemistry. We designed and characterized a new experimental apparatus for probing chemical reaction dynamics between molecular cations and neutral radicals at temperatures below 1 K. The ions are trapped in a linear quadrupole radio-frequency trap and sympathetically cooled by co-trapped, laser-cooled, atomic ions. The ion trap is coupled to a time-of-flight mass spectrometer to readily identify product ion species and to accurately determine trapped ion numbers. We discuss, and present in detail, the design of this ion trap time-of-flight mass spectrometer and the electronics required for driving the trap and mass spectrometer. Furthermore, we measure the performance of this system, which yields mass resolutions of m/Δm ≥ 1100 over a wide mass range, and discuss its relevance for future measurements in chemical reaction kinetics and dynamics.

Development of a wavy Stark velocity filter for studying interstellar chemistry

Cold polar molecules are key to both the understanding of fundamental physics and the characterization of the chemical evolution of interstellar clouds. To facilitate such studies over a wide range of temperatures, we developed a new type of Stark velocity filter for changing the translational and rotational temperatures of velocity-selected polar molecules without changing the output beam position. The translational temperature of guided polar molecules can be significantly varied by exchanging the wavy deflection section with one having a different radius of the curvature and a different deflection angle. Combining in addition a temperature variable gas cell with the wavy Stark velocity filter enables to observe the translational and rotational temperature dependence of the reaction-rate constants of cold ion-polar molecule reactions over the interesting temperature range of 10-100 K.

Low-temperature kinetics and dynamics with Coulomb crystals

Coulomb crystals-as a source of translationally cold, highly localized ions-are being increasingly utilized in the investigation of ion-molecule reaction dynamics in the cold regime. To develop a fundamental understanding of ion-molecule reactions, and to challenge existing models that describe the rates, product branching ratios, and temperature dependence of such processes, investigators need to exercise full control over the experimental reaction parameters. This requires not only state selection of the reactants, but also control over the collision process (e.g., the collisional energy and angular momentum) and state-selective product detection. The combination of Coulomb crystals in ion traps with cold neutral-molecule sources is enabling the measurement of state-selective reaction rates in a diverse range of systems. With the development of appropriate product detection techniques, we are moving toward the ultimate goal of examining low-energy, state-to-state ion-molecule reaction dynamics.

Getting a grip on the transverse motion in a Zeeman decelerator

Zeeman deceleration is an experimental technique in which inhomogeneous, time-dependent magnetic fields generated inside an array of solenoid coils are used to manipulate the velocity of a supersonic beam. A 12-stage Zeeman decelerator has been built and characterized using hydrogen atoms as a test system. The instrument has several original features including the possibility to replace each deceleration coil individually. In this article, we give a detailed description of the experimental setup, and illustrate its performance. We demonstrate that the overall acceptance in a Zeeman decelerator can be significantly increased with only minor changes to the setup itself. This is achieved by applying a rather low, anti-parallel magnetic field in one of the solenoid coils that forms a temporally varying quadrupole field, and improves particle confinement in the transverse direction. The results are reproduced by three-dimensional numerical particle trajectory simulations thus allowing for a rigorous analysis of the experimental data. The findings suggest the use of a modified coil configuration to improve transverse focusing during the deceleration process.

Cold state-selected molecular collisions and reactions

Over the past decade, and particularly the past five years, a quiet revolution has been building at the border between atomic physics and experimental quantum chemistry. The rapid development of techniques for producing cold and even ultracold molecules without a perturbing rare-gas cluster shell is now enabling the study of chemical reactions and scattering at the quantum scattering limit with only a few partial waves contributing to the incident channel. Moreover, the ability to perform these experiments with nonthermal distributions comprising one or a few specific states enables the observation and even full control of state-to-state collision rates in this computation-friendly regime: This is perhaps the most elementary study possible of scattering and reaction dynamics.

Taming molecular collisions using electric and magnetic fields

In molecular collision experiments, studying the collision process in high detail requires controlling molecular degrees of freedom before the collision.

Millikelvin reactive collisions between sympathetically cooled molecular ions and laser-cooled atoms in an ion-atom hybrid trap

We report on a study of cold reactive collisions between sympathetically cooled molecular ions and laser-cooled atoms in an ion-atom hybrid trap. Chemical reactions were studied at average collision energies <E(coll)>/k(B)>/~20 mK, about 2 orders of magnitude lower than has been achieved in previous experiments with molecular ions. Choosing N(2)(+)+Rb as a prototypical system, we find that the reaction rate is independent of the collision energy within the range studied, but strongly dependent on the internal state of Rb. Highly efficient charge exchange four times faster than the Langevin rate was observed with Rb in the excited (5p) (2)P(3/2) state. This observation is rationalized by a capture process dominated by the charge-quadrupole interaction and a near resonance between the entrance and exit channels of the system. Our results provide a test of classical models for reactions of molecular ions at the lowest energies reached thus far.

Light-assisted ion-neutral reactive processes in the cold regime: radiative molecule formation versus charge exchange

We present a combined experimental and theoretical study of cold reactive collisions between laser-cooled Ca+ ions and Rb atoms in an ion-atom hybrid trap. We observe rich chemical dynamics which are interpreted in terms of nonadiabatic and radiative charge exchange as well as radiative molecule formation using high-level electronic structure calculations. We study the role of light-assisted processes and show that the efficiency of the dominant chemical pathways is considerably enhanced in excited reaction channels. Our results illustrate the importance of radiative and nonradiative processes for the cold chemistry occurring in ion-atom hybrid traps.

Direct and inverse reactions of LiH+ with He(1S) from quantum calculations: mechanisms and rates

The gas-phase reaction of LiH(+) (X(2)Σ) with He((1)S) atoms, yielding Li(+)He with a small endothermicity for the rotovibrational ground state of the reagents, is analysed using the quantum reactive approach that employs the Negative Imaginary Potential (NIP) scheme discussed earlier in the literature. The dependence of low-T rates on the initial vibrational state of LiH(+) is analysed and the role of low-energy Feshbach resonances is also discussed. The inverse destruction reaction of LiHe(+), a markedly exothermic process, is also investigated and the rates are computed in the same range of temperatures. The possible roles of these reactions in early universe astrophysical networks, in He droplets environments or in cold traps are briefly discussed.

Theoretical characterization of laser- and sympathetically-cooled ions in surface-electrode ion traps

Using molecular dynamics simulations we characterize theoretically Coulomb clusters of laser- and sympathetically-cooled ions in a five-wire surface-electrode ion trap. We show that the asymmetry of the trapping potential leads to significantly different cluster structures and ion energy distributions in comparison to conventionally used linear Paul traps and to an asymmetric segregation of the ions in bi-component Coulomb clusters. We explore the impact of our results on the implementation of sympathetic cooling of molecular ions in surface-electrode traps and discuss possible challenges for the realization of such experiments.

Dynamics of individual rotational states in an electrostatic guide for neutral molecules

The guiding properties of individual rotational states of deuterated ammonia inside an electrostatic hexapole guide are presented. The guide is combined with resonance enhanced multiphoton ionization detection to assess the guiding probabilities and velocity distributions as a function of the rotational quantum numbers J and K. Due to the differences in the effective dipole moment these states are prepared at significantly different translational temperatures. A model is presented that describes the velocity-distribution for individual M-sublevels, and this model is also used to determine a rotational-state dependent translational temperature. Furthermore, the hexapole field has been replaced by a dipole field in order to obtain a band-pass velocity filter. However, the resulting change in the final velocity distribution is similar to that obtained from a hexapole guide but with increased backing pressure, leading to collisional acceleration of the slow molecules.

Cold chemistry with electronically excited Ca+ Coulomb crystals

Rate constants for chemical reactions of laser-cooled Ca(+) ions and neutral polar molecules (CH(3)F, CH(2)F(2), or CH(3)Cl) have been measured at low collision energies (<E(coll)>/k(B)=5-243 K). Low kinetic energy ensembles of (40)Ca(+) ions are prepared through Doppler laser cooling to form "Coulomb crystals" in which the ions form a latticelike arrangement in the trapping potential. The trapped ions react with translationally cold beams of polar molecules produced by a quadrupole guide velocity selector or with room-temperature gas admitted into the vacuum chamber. Imaging of the Ca(+) ion fluorescence allows the progress of the reaction to be monitored. Product ions are sympathetically cooled into the crystal structure and are unambiguously identified through resonance-excitation mass spectrometry using just two trapped ions. Variations of the laser-cooling parameters are shown to result in different steady-state populations of the electronic states of (40)Ca(+) involved in the laser-cooling cycle, and these are modeled by solving the optical Bloch equations for the eight-level system. Systematic variation of the steady-state populations over a series of reaction experiments allows the extraction of bimolecular rate constants for reactions of the ground state ((2)S(1/2)) and the combined excited states ((2)D(3/2) and (2)P(1/2)) of (40)Ca(+). These results are analyzed in the context of capture theories and ab initio electronic structure calculations of the reaction profiles. In each case, suppression of the ground state rate constant is explained by the presence of a submerged or real barrier on the ground state potential surface. Rate constants for the excited states are generally found to be in line with capture theories.