Cationic Noble-Gas Hydrides: From Ion Sources to Outer Space

Microsolvation of a Proton by Ar Atoms: Structures and Energetics of ArnH+ Clusters

We present a computational investigation on the structural arrangements and energetic stabilities of small-size protonated argon clusters, Ar nH +. Using high-level ab initio electronic structure computations, we determined that the linear symmetric triatomic ArH +Ar ion serves as the molecular core for all larger clusters studied. Through harmonic normal-mode analysis for clusters containing up to seven argon atoms, we observed that the proton-shared vibration shifts to lower frequencies, consistent with measurements in gas-phase IRPD and solid Ar-matrix isolation experiments. We explored the sum-of-potentials approach by employing kernel-based machine-learning potential models trained on CCSD(T)-F12 data. These models included expansions of up to two-body, three-body, and four-body terms to represent the underlying interactions as the number of Ar atoms increases. Our results indicate that the four-body contributions are crucial for accurately describing the potential surfaces in clusters with n> 3. Using these potential models and an evolutionary programming method, we analyzed the structural stability of clusters with up to 24 Ar atoms. The most energetically favored Ar nH + structures were identified for magic size clusters at n = 7, 13, and 19, corresponding to the formation of Ar-pentagon rings perpendicular to the ArH +Ar core ion axis. The sequential formation of such regular shell structures is compared to ion yield data from high-resolution mass spectrometry measurements. Our results demonstrate the effectiveness of the developed sum-of-potentials model in describing trends in the nature of bonding during the single proton microsolvation by Ar atoms, encouraging further quantum nuclear studies.

Noble Gas Anions: An Overview of Strategies and Bonding Motifs

This review article aims to provide an overview of the strategies employed to prepare noble gas anions under different environments and experimental conditions, and of the bonding motifs typically occurring in these species. Observed systems include anions fixed into synthesized salts, detected in the gas phase or in high-pressure devices. The major role of the theoretical calculations is also highlighted, not only in support of the experiments, but also as effective in predicting still unreported species. The chemistry of noble gas anions overall appears as a varied and rich paint, offering fascinating opportunities for both experimentalists and theoreticians.

Chemical Aristocracy: He3 Dication and Analogous Noble-Gas-Exclusive Covalent Compounds

Herein, we predict the first set of covalently bonded triatomic molecular compounds composed exclusively of noble gases. Using a combination of double-hybrid DFT, CCSD(T), and MRCI+Q calculations and a range of bonding analyses, we explored a set of 270 doubly charged triatomics, which included various combinations of noble gases and main group elements. This extensive exploration uncovered nine noble-gas-exclusive covalent compounds incorporating helium, neon, argon, or combinations thereof, exemplified by cases such as He32+ and related systems. This work brings to light a previously uncharted domain of noble gas chemistry, demonstrating the potential of noble gases in forming covalent molecular clusters.

Metastable Liquid Properties and Surface Flow Patterns of Ultrahigh Temperature Alloys Explored in Outer Space

The metastable liquid properties and chemical bonds beyond 2000 K remain a huge challenge for ground-based research on liquid materials chemistry. We show the strong undercooling capability, metastable liquid properties and surface wave patterns of refractory Nb-Si and Zr-V binary alloys explored in space environment. The floating droplet of Nb82.7Si17.3 eutectic alloy superheated up to 2338 K exhibited an extreme undercooling of 437 K, approaching the 0.2TE threshold for homogeneous nucleation of liquid-solid reaction. The microgravity state endowed alloy droplets with nearly perfect sphericity and thus ensured the high accuracy to determine metastable undercooled liquid properties. A special kind of swirling flow was induced for liquid alloy owing to Marangoni convection, which resulted in the spiral microstructures on Zr64V36 alloy surface during liquid-solid phase transition. The coupled impacts of surface nucleation and surface flow brought in a novel olivary shape for these binary alloys. Furthermore, the chemical bonds and atomic structures of high temperature liquids were revealed to understand the liquid properties in outer space circumstances.

Noble gas hydrides: theoretical prediction of the first group of anionic species

The first group of anionic noble-gas hydrides with the general formula HNgBeO- (Ng = Ar, Kr, Xe, Rn) is predicted through MP2, Coupled-Cluster, and Density Functional Theory computations employing correlation-consistent atomic basis sets. We derive that these species are stable with respect to the loss of H, H-, BeO, and BeO-, but unstable with respect to Ng + HBeO-. The energy barriers of the latter process are, however, high enough to suggest the conceivable existence of the heaviest HNgBeO- species as metastable in nature. Their stability arises from the interaction of the H- moiety with the positively-charged Ng atoms, particularly with the σ-hole ensuing from their ligation to BeO. This actually promotes relatively tight Ng-H bonds featuring a partially-covalent character, whose degree progressively increases when going from HArBeO- to HRnBeO-. The HNgBeO- compounds are also briefly compared with other noble-gas anions observed in the gas phase or isolated in crystal lattices.

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.

Trihydrogen Cation Helium Clusters: A New Potential Energy Surface

We present a new analytical potential energy surface (PES) for the interaction between the trihydrogen cation and a He atom,H 3 + - H e ${{H}_{3}^{+}-He}$ , in its electronic ground state. The proposed PES has been built as a sum of two contributions: a polarization energy term due to the electric field generated by the molecular cation at the position of the polarizable He atom, and an exchange-repulsion and dispersion interactions represented by a sum of "atom-bond" potentials between the three bonds ofH 3 + ${{H}_{3}^{+}}$ and the He atom. All parameters of this new PES have been chosen and fitted from data obtained from high-level ab-initio calculations. Using this new PES plus the Aziz-Slaman potential for the interaction between Helium atoms and assuming pair-wise interactions, we carry out classical Basin-Hopping (BH) global optimization, semiclassical BH with Zero Point Energy corrections, and quantum Diffusion Monte Carlo simulations. We have found the minimum energy configurations of small He clusters doped withH 3 + ${{H}_{3}^{+}}$ ,H 3 + H e N ${{H}_{3}^{+}{\left(He\right)}_{N}}$ , with N=1-16. The study of the energies of these clusters allows us to find a pronounced anomaly for N=12, in perfect agreement with previous experimental findings, which we relate to a greater relative stability of this aggregate.

Quantum computations in heavy noble-gas hydride cations: Reference energies and new spectroscopic data

Computational quantum chemistry has become a powerful tool with a wide range of possibilities to solve chemical-physical problems. As a result of this, the interest in the applications of computational quantum chemistry has expanded considerably, and has opened up novel research opportunities. In particular, those related to the characterization of heavy-atoms complexes, as most electronic structure calculations for such systems struggle with the problem posed by the large number of electrons present in them, and consequently, the introduction of relativistic effects. The present study performed an exhaustive assess to characterized the uncommon NgH+ (Ng = Kr, Xe, and Rn) hydride cations in order to provide accurate rovibrational data of their isotopes to assist in the laboratory characterization or even their astronomical detection. Scalar relativistic effects were included, and the ground and first electronically exited states potential curves were obtained from benchmark ab initio CCSD(T)/CBS and MRCI+Q electronic structure calculations. Next, such interaction potentials, correctly extended to long-range asymptotic regions, were employed in quantum bound state calculations and molecular spectroscopic constants were determined for the most abundance 84Kr, 132Xe, and 222Rn isotopes. Our results were discussed in comparison with available experimental and previous theoretical estimates, aiming to treat accuracy issues. The new sets provide reference data that could serve for spectroscopic characterization of such low abundance and high radioactive species.

Hitting the Bull's Eye: Stable HeBeOH+ Complex

It is now known that the heavier noble gases (Ng=Ar-Rn) show some varying degrees of reactivity with a gradual increase in reactivity along Ar-Rn. However, because of their very small size and very high ionization potential, helium and neon are the hardest targets to crack. Although few neon complexes are isolated at very low temperatures, helium needs very extreme situations like very high pressure. Here, we find that protonated BeO, BeOH⁺ can bind helium and neon spontaneously at room temperature. Therefore, extreme conditions like very low temperature and/or high pressure will not be required for their experimental isolation. The Ng-Be bond strength is very high for their heavier homologs and the bond strength shows a gradual increase from He to Rn. Moreover, the Ng-Be attractive energy is almost exclusively originated from the orbital interaction which is composed of one Ng(s/p_σ )→BeOH⁺ σ-donation and two weaker Ng(p_π )→BeOH⁺ π-donations, except for helium. Helium uses its low-lying vacant 2p orbitals to accept π-electron density from BeOH⁺ . Previously, such electron-accepting ability of helium was used to explain a somewhat stronger helium bond than neon for neutral complexes. However, the present results indicate that such π-back donations are too weak in nature to decide any energetic trend between helium and neon.

Superstrong Chemical Bonding of Noble Gases with Oxidoboron (BO+) and Sulfidoboron (BS+)

Inspired by the overwhelming exploration of noble gas-boron (Ng-B) bond containing chemical compounds, the stability of the Ng bound BY⁺ and AlY⁺ (Y = O and S) has been investigated by using various ab initio based quantum chemical methods. Ng atoms are found to form exceptionally strong bonds with BO⁺ species in the predicted NgBO⁺ (Ng = He-Rn) complexes with remarkably high Ng-B dissociation energies ranging from 138.0 to 462.2 kJ mol^-1 for the He-Rn series. It is the highest ever Ng-B binding energy in conjunction with the smallest Ng-B bond length for any of the cationic species involving a Ng-B bond as reported until today. More importantly, the calculated Ng-B bond lengths have been found to be much lower than the respective covalent limits in both NgBO⁺ and NgBS⁺ ions. The electronegativity difference between O and S atoms has been reflected nicely in the Ng-B and Ng-Al binding energies, which are found to be 91.9-346.5, 9.6-169.2, and 6.8-142.1 kJ mol^-1 in NgBS⁺, NgAlO⁺, and NgAlS⁺, respectively. The strong covalent bonding between Ng and B/Al atoms in the predicted chemical systems has also been supported by the natural bonding orbital (NBO) and electron density based atoms-in-molecule (AIM) analysis. In addition, the energy decomposition analysis (EDA) in combination with the natural bond orbital for chemical valence (NOCV) indicates that the orbital interaction term is the prime contributor to the total attraction energy in the Ng-B and Ng-Al bonds. Furthermore, Ng-B and Ng-Al bonding can be assessed using the donor-acceptor model where the σ-electron donation that takes place from Ng (HOMO) → XY⁺ (LUMO) (X = B and Al; Y = O and S) is the major contributor to the orbital interaction energy. All the computational results along with the very recent experimental observation of ArOH⁺ and NgMX (Ng = Ar-Xe; M = Cu, Ag, Au; X = F, Cl) clearly indicate that it might be possible to synthesize and characterize these superstrong complexes, NgXY⁺ (Ng = He-Rn; X = B and Al; Y = O and S), under suitable experimental technique(s).

Noble Gas-Silicon Cations: Theoretical Insights into the Nature of the Bond

The structure, stability, and bonding situation of some exemplary noble gas-silicon cations were investigated at the MP2/aVTZ level of theory. The explored species include the mono-coordinated NgSiX3+ (Ng = He-Rn; X = H, F, Cl) and NgSiF22+ (Ng = He-Rn), the di-coordinated Ar2SiX3+ (X = H, F, Cl), and the “inserted” FNgSiF2+ (Ng = Kr, Xe, Rn). The bonding analysis was accomplished by the method that we recently proposed to assay the bonding situation of noblegas compounds. The Ng-Si bonds are generally tight and feature a partial contribution of covalency. In the NgSiX3+, the degree of the Ng-Si interaction mirrors the trends of two factors, namely the polarizability of Ng that increases when going from Ng = He to Ng = Rn, and the Lewis acidity of SiX3+ that decreases in the order SiF3+ > SiH3+ > SiCl3+. For the HeSiX3+, it was also possible to catch peculiar effects referable to the small size of He. When going from the NgSiF3+ to the NgSiF22+, the increased charge on Si promotes an appreciable increase inthe Ng-Si interaction, which becomes truly covalent for the heaviest Ng. The strength of the bond also increases when going from the NgSiF3+ to the “inserted” FNgSiF2+, likely due to the cooperative effect of the adjacent F atom. On the other hand, the ligation of a second Ar atom to ArSiX3+ (X = H, F, Cl), as to form Ar2(SiX3+), produces a weakening of the bond. Our obtained data were compared with previous findings already available in the literature.

Stability-Order Reversal in FSiY and FYSi (Y = N and P) Molecules after the Insertion of a Noble Gas Atom

Recent theoretical prediction and experimental identification of fluorinated noble gas cyanides and isocyanides motivate us to explore a unique novel series of neutral noble gas-inserted heavier cyanofluoride isomers, FNgYSi and FNgSiY (Ng = Kr, Xe, and Rn; Y = N and P), theoretically using quantum chemical calculations. The concerned minima and saddle point geometries have been optimized using DFT, MP2, and CCSD(T) methods. The precursor molecule FSiY is more stable than its isomer FYSi, and the stability order is found to be reversed after the insertion of a noble gas (Ng) atom into them which is in contrast to the previously reported FCN/FNC systems where the stability order in the precursors remains intact after the insertion of a Ng atom into them. The predicted FNgYSi molecules are metastable in nature as they are kinetically stable but thermodynamically unstable with respect to the global minima products (FYSi and Ng). All the calculations for the corresponding FNgSiY molecules clearly indicate that the less stable FNgSiY behaves similarly to the FNgYSi in all respects. The energetics, force constant, and spectroscopic data strongly reinforce the possibility of occurrence of these predicted FNgYSi and FNgSiY molecules which might be experimentally realized under suitable cryogenic condition(s).

Noble-gas compounds: A general procedure of bonding analysis

This paper accounts for a general procedure of bonding analysis that is, expectedly, adequate to describe any type of interaction involving the noble-gas (Ng) atoms. Building on our recently proposed classification of the Ng-X bonds (X = binding partner) [New J. Chem. 44, 15536 (2020)], these contacts are first distinguished into three types, namely, A, B, or C, based on the topology of the electron energy density H(r) and on the shape of its plotted form. Bonds of type B or C are, then, further assigned as B-loose (B^l) or B-tight (B^t) and C-loose (C^l) or C-tight (C^t) depending on the sign that H(r) takes along the Ng-X bond path located from the topological analysis of ρ(r), particularly at around the bond critical point (BCP). Any bond of type A, B^l/B^t, or C^l/C^t is, finally, assayed in terms of contribution of covalency. This is accomplished by studying the maximum, minimum, and average value of H(r) over the volume enclosed by the low-density reduced density gradient (RDG) isosurface associated with the bond (typically, the RDG isosurface including the BCP) and the average ρ(r) over the same volume. The bond assignment is also corroborated by calculating the values of quantitative indices specifically defined for the various types of interactions (A, B, or C). The generality of our taken approach should encourage its wide application to the study of Ng compounds.

A new krypton complex – Experimental and computational investigation of the krypton sulphur pentafluoride cation, [KrSF5]+, in the gas phase

The gas-phase reactions of noble gas (Ng) cations, namely Kr+ and Xe+, with SF6 were investigated experimentally by Fourier transform ion cyclotron resonance mass spectrometry and computationally using MP2 and...

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.

Computational Characterization of Astrophysical Species: The Case of Noble Gas Hydride Cations

Theoretical-computational studies together with recent astronomical observations have shown that under extreme conditions in the interstellar medium (ISM), complexes of noble gases may be formed. Such observations have generated a wide range of possibilities. In order to identify new species containing such atoms, the present study gathers spectroscopic data for noble gas hydride cations, NgH+ (Ng = He, Ne, Ar) from high-level ab initio quantum chemistry computations, aiming to contribute in understanding the chemical bonding and electron sharing in these systems. The interaction potentials are obtained from CCSD(T)/CBS and MRCI+Q calculations using large basis sets, and then employed to compute vibrational levels and molecular spectroscopic constants for all known stable isotopologues of ground state NgH+ cations. Comparisons with previously reported values available are discussed, indicating that the present data could serve as a benchmark for future studies on these systems and on higher-order cationic noble gas hydrides of astrophysical interest.

Dynamics of HeHHe+ Rotational State Changes Induced by Collision with He: A Possible New Path in Early Universe Chemistry

Ab initio calculations are employed to generate the rigid rotor (RR) potential energy surface (PES) describing the interaction of the linear molecular cation HeHHe+, at its equilibrium geometry, with the neutral He atom. The resulting interaction is employed to investigate the efficiency of rotational state-changing collisions at the temperatures relevant to the early universe conditions, where the latter molecule has been postulated to exist, albeit not yet observed. The inelastic rate coefficients are found to be fairly large and are compared with those found for another important cation just recently observed in the interstellar medium: the HeH+ polar molecule. The possibility for this cation to provide new options to energy dissipation routes under early universe conditions after the recombination era is briefly discussed.

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.

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

A new analytical potential energy surface (PES) has been constructed for the Ar₂H⁺ system from a dataset consisting of a large number of ab initio energies computed using the coupled-cluster singles, doubles and perturbative triples method and aug-cc-pVQZ basis set. The long-range interaction is added to the diatomic potentials using a standard long range expansion form to better describe the asymptotic regions. The vibrational states for the most stable structures of the Ar₂H⁺ system have been calculated, and few low lying states are assigned to quantum numbers. Reactive scattering studies have been performed for the Ar + Ar'H⁺ → Ar' + ArH⁺ proton exchange reaction on the newly generated PES. Reaction probability, cross sections, and rate constants are calculated for the Ar + Ar'H⁺(v = 0, j = 0) collisions within 0.01 eV-0.6 eV of relative translational energy using exact quantum dynamical simulations as well as quasiclassical trajectory (QCT) calculations. The effect of vibrational excitation of the reactants is also explored for the reaction. State averaged rate constants are calculated for the proton exchange reaction at different temperatures using the QCT method. The mechanistic pathways for the reaction are understood by analyzing the quasiclassical trajectories.

Spectroscopic signatures of HHe2+ and HHe3. Phys Chem Chem Phys 2020;22:22885-22888. [PMID: 33034329 DOI: 10.1039/d0cp04649c]

Using two different action spectroscopic techniques, a high-resolution quantum cascade laser operating around 1300 cm-1 and a cryogenic ion trap machine, the proton shuttle motion of the cations HHe2+ and HHe3+ has been probed at a nominal temperature of 4 K. For HHe3+, the loosely bound character of this complex allowed predissociation spectroscopy to be used, and the observed broad features point to a lifetime of a few ps in the vibrationally excited state. For He-H+-He, a fundamental linear molecule consisting of only three nuclei and four electrons, the method of laser-induced inhibition of complex growth (LIICG) enabled the measurement of three accurate rovibrational transitions, pinning down its molecular parameters for the first time.