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.

Intermolecular dynamics of NH3-rare gas complexes in the ν2 umbrella region of NH3 investigated by rovibrational laser jet-cooled spectroscopy and ab initio calculations.

High resolution jet-cooled spectroscopy experiments have been realized to investigate the intermolecular dynamics of van der Waals (vdW) heterodimers between NH3 and rare gas (Rg) atoms in the ν2 umbrella mode region of NH3. With respect to a previous study dedicated to NH3-Ar [Asselin et al. Mol. Phys. 116, 3642 (2018)], the sensitivity and spectral resolution of our laser spectrometer coupled to a pulsed supersonic jet have been significantly improved to derive more accurate excited state spectroscopic parameters from rovibrational analyses. In addition, we calculated the ground and ν2 excited vibration-rotation-tunneling (VRT) states of these complexes on the four-dimensional ab initio potential energy surfaces from Loreau et al. [J. Chem. Phys. 141, 224303 (2014), ibid. 143, 184303 (2015).] Transition frequencies and intensities of the allowed ν2 = 1 ← 0 transitions obtained from the calculated energy levels and wave functions agree well with the experimental data and are helpful in their analysis. By means of a pseudodiatomic model with the assumption of weak Coriolis coupling, the rovibrational analysis of both the Πe/f(j = 1,k = 0) ←Σf(j = 0,k = 0) and Σf(j = 1,k = 0) ←Σf(j = 0,k = 0) transitions in ortho NH3-Rg (Rg = Ne, Ar, Kr, Xe) complexes enabled us to determine reliable excited state parameters and derive accurate values of the effective vdW bond length Reff, force constant ks, and vdW stretching frequency νs. Comparison between the experimental structural parameters and those from the ab initio calculated VRT levels shows good agreement for NH3-Ne, NH3-Ar and NH3-Kr, and a similar variation of Reff, ks, and νs with the polarizability of Rg in the ground and ν2 excited states. Anomalously small values of νs and ks derived for NH3-Xe in the Πe/f(j = 1,k = 0) state suggest that the applied model is not valid in this case, due to the presence of another state coupling to the perturbed Πf state. Such a state could not be found, however.

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.

Photodissociation as a probe of the H3+ avoided crossing seam.

Experiments are conducted to investigate the role of the avoided crossing seam in the photodissociation of H3+. Three-dimensional imaging of dissociation products is used to determine the kinetic energy release and branching ratio among the fragmentation channels. Vibrational distributions are measured by dissociative charge transfer of H2+ products. It is found that the photodissociation of hot H3+ in the near-ultraviolet produces cold H2+, but hot H2. Modelling the wavepacket dynamics along the repulsive potential energy surface accounts for the repopulation of the ground potential energy surface. The role of the avoided crossing seam is emphasized and its importance for the astrophysically relevant charge transfer reactions underlined. This article is part of a discussion meeting issue 'Advances in hydrogen molecular ions: H3+, H5+ and beyond'.

Scattering of CO with H2O: Statistical and classical alternatives to close-coupling calculations.

Energy transfer in inelastic atom-molecule and molecule-molecule collisions can be described theoretically using the quantum-mechanical close-coupling method. Unfortunately, for bimolecular collisions implying heavy colliders and/or for which the potential energy surface has a deep well, the resulting coupled equations become numerically intractable and approximate methods have to be employed. H2O-CO collisions provide an important example for which close-coupling calculations are not feasible. In this paper, we investigate the accuracy of three approximate methods (the coupled states method, the quasi-classical trajectory method, and the statistical adiabatic channel model) to describe inelastic collisions of H2O with CO. We perform scattering calculations on a recent 5D potential energy surface, and we compare the results of the three approximate methods to fully converged close-coupling calculations at energies below 300 cm-1 and at low values of the total angular momentum. We show that the statistical method provides an attractive alternative to fully quantum mechanical close-coupling calculations at low collision energies, while the quasi-classical method is more advantageous at high energies.

Photochemistry of forbidden oxygen lines in the inner coma of 67P/Churyumov-Gerasimenko.

Observations of the green and red-doublet emission lines have previously been realized for several comets. We present here a chemistry-emission coupled model to study the production and loss mechanisms of the O(1S) and O(1D) states, which are responsible for the emission lines of interest for comet 67P/Churyumov-Gerasimenko. The recent discovery of O2 in significant abundance relative to water 3.80 ± 0.85% within the coma of 67P has been taken into consideration for the first time in such models. We evaluate the effect of the presence of O2 on the green to red-doublet emission intensity ratio, which is traditionally used to assess the CO2 abundance within cometary atmospheres. Model simulations, solving the continuity equation with transport, show that not taking O2 into account leads to an underestimation of the CO2 abundance within 67P, with a relative error of about 25%. This strongly suggests that the green to red-doublet emission intensity ratio alone is not a proper tool for determining the CO2 abundance, as previously suggested. Indeed, there is no compelling reason why O2 would not be a common cometary volatile, making revision of earlier assessments regarding the CO2 abundance in cometary atmospheres necessary. The large uncertainties of the CO2 photodissociation cross section imply that more studies are required in order to better constrain the O(1S) and O(1D) production through this mechanism. Space weather phenomena, such as powerful solar flares, could be used as tools for doing so, providing additional information on a good estimation of the O2 abundance within cometary atmospheres.

Scattering of NH3 and ND3 with rare gas atoms at low collision energy.

We present a theoretical study of elastic and rotationally inelastic collisions of NH3 and ND3 with rare gas atoms (He, Ne, Ar, Kr, Xe) at low energy. Quantum close-coupling calculations have been performed for energies between 0.001 and 300 cm(-1). We focus on collisions in which NH3 is initially in the upper state of the inversion doublet with j = 1, k = 1, which is the most relevant in an experimental context as it can be trapped electrostatically and Stark-decelerated. We discuss the presence of resonances in the elastic and inelastic cross sections, as well as the trends in the inelastic cross sections along the rare gas series and the differences between NH3 and ND3 as a colliding partner. We also demonstrate the importance of explicitly taking into account the umbrella (inversion) motion of NH3 in order to obtain accurate scattering cross sections at low collision energy. Finally, we investigate the possibility of sympathetic cooling of ammonia using cold or ultracold rare gas atoms. We show that some systems exhibit a large ratio of elastic to inelastic cross sections in the cold regime, which is promising for sympathetic cooling experiments. The close-coupling calculations are based on previously reported ab initio potential energy surfaces for NH3-He and NH3-Ar, as well as on new, four-dimensional, potential energy surfaces for the interaction of ammonia with Ne, Kr, and Xe, which were computed using the coupled-cluster method and large basis sets. We compare the properties of the potential energy surfaces corresponding to the interaction of ammonia with the various rare gas atoms.

Photodissociation of the carbon monoxide dication in the (3)Σ(-) manifold: Quantum control simulation towards the C(2+) + O channel.

The photodissociation and laser assisted dissociation of the carbon monoxide dication X(3)Π CO(2+) into the (3)Σ(-) states are investigated. Ab initio electronic structure calculations of the adiabatic potential energy curves, radial nonadiabatic couplings, and dipole moments for the X (3)Π state are performed for 13 excited (3)Σ(-) states of CO(2+). The photodissociation cross section, calculated by time-dependent methods, shows that the C(+) + O(+) channels dominate the process in the studied energy range. The carbon monoxide dication CO(2+) is an interesting candidate for control because it can be produced in a single, long lived, v = 0 vibrational state due to the instability of all the other excited vibrational states of the ground (3)Π electronic state. In a spectral range of about 25 eV, perpendicular transition dipoles couple this (3)Π state to a manifold of (3)Σ(-) excited states leading to numerous C(+) + O(+) channels and a single C(2+) + O channel. This unique channel is used as target for control calculations using local control theory. We illustrate the efficiency of this method in order to find a tailored electric field driving the photodissociation in a manifold of strongly interacting electronic states. The selected local pulses are then concatenated in a sequence inspired by the "laser distillation" strategy. Finally, the local pulse is compared with optimal control theory.

Potential energy surface and bound states of the NH3-Ar and ND3-Ar complexes.

A new, four-dimensional potential energy surface for the interaction of NH3 and ND3 with Ar is computed using the coupled-cluster method with single, double, and perturbative triple excitations and large basis sets. The umbrella motion of the ammonia molecule is explicitly taken into account. The bound states of both NH3-Ar and ND3-Ar are calculated on this potential for total angular momentum values from J = 0 to 10, with the inclusion of Coriolis interactions. The energies and splittings of the rovibrational levels are in excellent agreement with the extensive high-resolution spectroscopic data accumulated over the years in the infrared and microwave regions for both complexes, which demonstrates the quality of the potential energy surface.

Potential energy curves for the interaction of Ag(5s) and Ag(5p) with noble gas atoms.

We investigate the interaction of ground and excited states of a silver atom with noble gases (NG), including helium. Born-Oppenheimer potential energy curves are calculated with quantum chemistry methods and spin-orbit effects in the excited states are included by assuming a spin-orbit splitting independent of the internuclear distance. We compare our results with experimentally available spectroscopic data, as well as with previous calculations. Because of strong spin-orbit interactions, excited Ag-NG potential energy curves cannot be fitted to Morse-like potentials. We find that the labeling of the observed vibrational levels has to be shifted by one unit.

Photodissociation and radiative association of HeH+ in the metastable triplet state.

We investigate the photodissociation of HeH(+) in the metastable triplet state as well as its formation through the inverse process, radiative association. In models of astrophysical plasmas, HeH(+) is assumed to be present only in the ground state, and the influence of the triplet state has not been explored. It may be formed by radiative association during collisions between a proton and metastable helium, which are present in significant concentrations in nebulae. The triplet state can also be formed by association of He(+) and H, although this process is less likely to occur. We compute the cross sections and rate coefficients corresponding to the photodissociation of the triplet state by UV photons from a central star using a wave packet method. We show that the photodissociation cross sections depend strongly on the initial vibrational state and that the effects of excited electronic states and nonadiabatic couplings cannot be neglected. We then calculate the cross section and rate coefficient for the radiative association of HeH(+) in the metastable triplet state.

Radiative lifetime of the a 3Σ+ state of HeH+ from ab initio calculations.

The first metastable triplet state of HeH(+) was found to be present in ion beam experiments, with its lifetime estimated to be between hundreds of milliseconds and thousand of seconds. In this work, we use ab initio methods to evaluate the radiative lifetimes of the six vibrational levels of the a (3)Σ(+) of HeH(+). The transition a (3)Σ(+)→X (1)Σ(+) is spin-forbidden, but acquires intensity through spin-orbit interaction with the singlet and triplet Π states. Large scale CASSCF/MRCI calculations using an adapted basis set were performed to determine the potential energy curves of the relevant states of HeH(+) as well as the matrix elements of the dipole and spin-orbit operators. The wave functions and energies of the vibrational levels of the a (3)Σ(+) and X (1)Σ(+) states are obtained using a B-spline method and compared to previous works. We find that the radiative lifetime of the vibrational levels increases strongly with v, the lifetime of the v=0 state being 150 s. We also analyze the contributions from discrete and continuum parts of the spectrum. With such a long lifetime, the a (3)Σ(+) state could have astrophysical implications.