Quantum and quasiclassical dynamical simulations for the Ar2H+ on a new global analytical potential energy surface

A kernel-based machine learning potential and quantum vibrational state analysis of the cationic Ar hydride (Ar2H+)

One of the most fascinating discoveries in recent years, in the cold and low pressure regions of the universe, was the detection of ArH+ and HeH+ species. The identification of such noble gas-containing molecules in space is the key to understanding noble gas chemistry. In the present work, we discuss the possibility of [Ar2H]+ existence as a potentially detectable molecule in the interstellar medium, providing new data on possible astronomical pathways and energetics of this compound. As a first step, a data-driven approach is proposed to construct a full 3D machine-learning potential energy surface (ML-PES) via the reproducing kernel Hilbert space (RKHS) method. The training and testing data sets are generated from CCSD(T)/CBS[56] computations, while a validation protocol is introduced to ensure the quality of the potential. In turn, the resulting ML-PES is employed to compute vibrational levels and molecular spectroscopic constants for the cation. In this way, the most common isotopologue in ISM, [36Ar2H]+, was characterized for the first time, while simultaneously, comparisons with previously reported values available for [40Ar2H]+ are discussed. Our present data could serve as a benchmark for future studies on this system, as well as on higher-order cationic Ar-hydrides of astrophysical interest.

Ar+ ArH+ Reactive Collisions of Astrophysical Interest: The Case of 36 Ar

The reactive collision between 36 Ar and the 36 ArH+ species has been investigated by means of quantum mechanical (QM), quasiclassical trajectories (QCT) and statistical quantum mechanical (SQM) approaches. Reaction probabilities, cross sections as a function of the energy and rate constants in terms of the temperature have been obtained. Cumulative distributions as a function of the collision time and the inspection of selected QCT corresponding to specific dynamical mechanisms have been analysed. Predictions by means of the SQM method are in good agreement with the QM results, thus supporting the complex-forming nature of the process.

Atom-Diatom Reactive Scattering Collisions in Protonated Rare Gas Systems

The study of the dynamics of atom-diatom reactions involving two rare gas (Rg) atoms and protons is of crucial importance given the astrophysical relevance of these processes. In a series of previous studies, we have been investigating a number of such Rg(1)+ Rg(2)H+→ Rg(2)+ Rg(1)H+ reactions by means of different numerical approaches. These investigations comprised the construction of accurate potential energy surfaces by means of ab initio calculations. In this work, we review the state-of-art of the study of these protonated Rg systems making special emphasis on the most relevant features regarding the dynamical mechanisms which govern these reactive collisions. The aim of this work therefore is to provide an as complete as possible description of the existing information regarding these processes.

Fermi resonance switching in KrH+Rg and XeH+Rg (Rg = Ne, Ar, Kr, and Xe)

Matrix isolation experiments have been successfully employed to extensively study the infrared spectrum of several proton-bound rare gas complexes. Most of these studies have focused on the spectral signature for the H⁺ stretch (ν₃) and its combination bands with the intermolecular stretch coordinate (ν₁). However, little attention has been paid to the Fermi resonance interaction between the H⁺ stretch (ν₃) and H⁺ bend overtone (2ν₂) in the asymmetric proton-bound rare gas dimers, RgH⁺Rg'. In this work, we have investigated this interaction on KrH⁺Rg and XeH⁺Rg with Rg = (Ne, Ar, Kr, and Xe). A multilevel potential energy surface (PES) was used to simulate the vibrational structure of these complexes. This PES is a dual-level comprising of second-order Møller-Plesset perturbation theory and coupled-cluster singles doubles with perturbative triples [CCSD(T)] levels of ab initio theories. We found that when both the combination bands (nν₁ + ν₃) and bend overtone 2ν₂ compete to borrow intensity from the ν₃ band, the latter wins over the former, which then results in the suppression of the nν₁ + ν₃ bands. The current simulations offer new assignments for the ArH⁺Xe and KrH⁺Xe spectra. Complete basis set (CBS) binding energies for these complexes were also calculated at the CCSD(T)/CBS level.