A combined experimental and theoretical investigation of Cs+ ions solvated in HeN clusters

Computational Modeling: Up-to-Date Approaches and Cutting-Edge Applications from Clusters, Nanostructures to Bulk Systems

The contributions in this special theme collection, in honor to Prof. P. Villarreal, cover a broad variety of computational methodologies and experimental techniques, containing studies on gas phase, clusters and condensed phase systems.

Time-resolved solvation of alkali ions in superfluid helium nanodroplets

The sinking of alkali cations in superfluid 4He nanodroplets is investigated theoretically using liquid 4He time-dependent density functional theory at zero temperature. The simulations illustrate the dynamics of the buildup of the first solvation shell around the ions. The number of helium atoms in this shell is found to linearly increase with time during the first stages of the dynamics. This points to a Poissonian capture process, as concluded in the work of Albrechtsen et al. on the primary steps of Na+ solvation in helium droplets [Albrechtsen et al., Nature 623, 319 (2023)]. The energy dissipation rate by helium atom ejection is found to be quite similar between all alkalis, the main difference being a larger energy dissipated per atom for the lighter alkalis at the beginning of the dynamics. In addition, the number of helium atoms in the first solvation shell is found to be lower at the end of the dynamics than at equilibrium for both Li+ and Na+, pointing to a kinetic rather than thermodynamical control of the snowball size for small and strongly attractive ions.

Structures and stability of K+ cation solvated in Arn clusters

The solvation of K+ cation plays an important role in various phenomena such as biological procedures, geological time, and archaeological properties. Monte Carlo (MC) simulation and DFT method are employed to study the structural and energetic characteristics of the K + Arn (n = 1-14) clusters. The potential model (PM) and the Basin-Hopping (BH) method are the foundation of the MC simulation. The pairwise PM (PW-PM) is improved by introducing the N-body interactions via the polarizable potential model (PPM). On the other side, the DFT functional M05-2X, combined with the 6-311++G(3d2f,2p) basis set, and the Grimme dispersion correction GD3 was used to deeply investigate the geometrical properties and the relative stability of the K + Arn clusters. Starting from n = 12, a structural transition from square antiprism (SA) to icosahedron (ICOS) form is detected. Additionally, the PPM allows us to examine the largest sizes (n = 15-54). Herein, the first ICOS layers are found for n = 12 and 54 cluster sizes, respectively. The binding energy and the second energy difference as a function of cluster size are used to evaluate the relative stability of K + Arn clusters. The obtained data are in concordance with the available results in the literature.

Observation of Multiple Ordered Solvation Shells in Doped Helium Droplets: The Case of HeNCa2

In this Letter, we report the experimental detection of likely the largest ordered structure of helium atoms surrounding a monatomic impurity observed to date using a recently developed technique. The mass spectrometry investigation of He_NCa²⁺ clusters, formed in multiply charged helium nanodroplets, reveals magic numbers at N = 12, 32, 44, and 74. Classical optimization and path integral Monte Carlo calculations suggest the existence of up to four shells surrounding the calcium dication which are closed with well-ordered Mozartkugel-like structures: He₁₂Ca²⁺ with an icosahedron, the second at He₃₂Ca²⁺ with a dodecahedron, the third at He₄₄Ca²⁺ with a larger icosahedron, and finally for He₇₄Ca²⁺, we find that the outermost He atoms form an icosidodecahedron which contains the other inner shells. We analyze the reasons for the formation of such ordered shells in order to guide the selection of possible candidates to exhibit a similar behavior.

Helium nanodroplets as an efficient tool to investigate hydrogen attachment to alkali cations

We report a novel method to reversibly attach and detach hydrogen molecules to positively charged sodium clusters formed inside a helium nanodroplet host matrix. It is based on the controlled production of multiply charged helium droplets which, after picking up sodium atoms and exposure to H2 vapor, lead to the formation of Nam+(H2)n clusters, whose population was accurately measured using a time-of-flight mass spectrometer. The mass spectra reveal particularly favorable Na+(H2)n and Na2+(H2)n clusters for specific "magic" numbers of attached hydrogen molecules. The energies and structures of these clusters have been investigated by means of quantum-mechanical calculations employing analytical interaction potentials based on ab initio electronic structure calculations. A good agreement is found between the experimental and the theoretical magic numbers.

Microsolvation of lithium cation in xenon clusters: An octahedral growth pattern

The structural and energetic proprieties for the Li ⁺ Xe_n (n = 1-18) clusters are investigated using both Basin-Hopping combined with Potential Model description (BH-PM) and DFT methods. A structural transition from tetrahedral (4 coordination) form to octahedral (6 coordination) one is observed for n = 6. Above this size, all structures have an octahedral core. The cubic-face-centered arrangement for xenon atoms is detected for Li ⁺ Xe₁₄. To the best of our knowledge, the Li ⁺ Xe_n (n = 1-18) clusters are studied in the present work for the first time using the DFT theoretical approach. The M062X functional combined with aug-cc-pVDZ (for Li) and def2-TZVP (for Xe) basis sets reproduces accurately the CCSD(T) potential energy curve of Li + Xe system. Atom-Centered Density Matrix Propagation (ADMP) molecular dynamic calculations have been carried. Moreover, we investigate the larger sizes n = 31-35, 44, and 55 for the first time using the BH-PM theoretical approach. The closing of the first and second octahedron shells are proved for the n = 6 and 34 sizes, respectively. The relative stabilities of the Li ⁺ Xe_n molecules are also studied by computing the total energy, the binding energy per atoms for each size n. Then, the second energy difference between the size n and its two near neighbors allows identifying the magic number series. Our present data are analyzed, discussed and compared with the available theoretical and experimental data.

Helium structures around SF5+ and SF6+: novel intermolecular potential and mass spectrometry experiments

Helium clusters around the recently experimentally observed sulphur hexafluoride SF₆⁺ and sulphur pentafluoride SF₅⁺ ions are investigated using a combined experimental and theoretical effort. Mass spectrometry ion yields are obtained and the energetics and structure of the corresponding He_N–SF₆⁺ and He_N–SF₅⁺ clusters are analyzed using path integral molecular dynamics calculations as a function of N, the number of He atoms, employing a new intermolecular potential describing the interaction between the dopant and the surrounding helium. The new force field is optimized on benchmark potential energy ab initio calculations and represented by improved Lennard-Jonnes analytical expressions. This procedure improves the previous potentials employed in similar simulations for neutral SF₆ attached to helium nanodroplets. The theoretical analysis explains the characteristic features observed in the experimental ion yields which suggest the existence of stable configurations at specific sizes.

The structure of the He atoms around SF₅⁺ and SF₆⁺ is investigated both experimentally and theoretically.

Solvation of Isoelectronic Halide and Alkali Metal Ions by Argon Atoms

This work systematically examines the interactions of alkali metal cations and their isoelectronic halide counterparts with up to six solvating Ar atoms (M⁺Ar_n and X^-Ar_n, where M = Li, Na, K, and Rb; X = H, F, Cl, and Br; and n = 1-6) via full geometry optimizations with the MP2 method and robust, correlation-consistent quadruple-ζ (QZ) basis sets. 116 unique M⁺Ar_n and X^-Ar_n stationary points have been characterized on the MP2/QZ potential energy surface. To the best of our knowledge, approximately two dozen of these stationary points have been reported here for the first time. Some of these new structures are either the lowest-energy stationary point for a particular cluster or energetically competitive with it. The CCSD(T) method was employed to perform additional single-point energy computations upon all MP2/QZ-optimized structures using the same basis set. CCSD(T)/QZ results indicate that internally solvated structures with the ion at/near the geometric center of the cluster have appreciably higher energies than those placing the ion on the periphery. While this study extends the prior investigations of M⁺Ar_n clusters found within the literature, it notably provides one of the first thorough characterizations of and comparisons to the corresponding negatively charged X^-Ar_n clusters.

Ca+ Ions Solvated in Helium Clusters

We present a combined experimental and theoretical investigation on Ca+ ions in helium droplets, HeNCa+. The clusters have been formed in the laboratory by means of electron-impact ionization of Ca-doped helium nanodroplets. Energies and structures of such complexes have been computed using various approaches such as path integral Monte Carlo, diffusion Monte Carlo and basin-hopping methods. The potential energy functions employed in these calculations consist of analytical expressions following an improved Lennard-Jones formula whose parameters are fine-tuned by exploiting ab initio estimations. Ion yields of HeNCa+ -obtained via high-resolution mass spectrometry- generally decrease with N with a more pronounced drop between N=17 and N=25, the computed quantum HeNCa+ evaporation energies resembling this behavior. The analysis of the energies and structures reveals that covering Ca+ with 17 He atoms leads to a cluster with one of the smallest energies per atom. As new atoms are added, they continue to fill the first shell at the expense of reducing its stability, until N=25, which corresponds to the maximum number of atoms in that shell. Behavior of the evaporation energies and radial densities suggests liquid-like cluster structures.

Leading Interaction Components in the Structure and Reactivity of Noble Gases Compounds

The nature, strength, range and role of the bonds in adducts of noble gas atoms with both neutral and ionic partners have been investigated by exploiting a fine-tuned integrated phenomenological–theoretical approach. The identification of the leading interaction components in the noble gases adducts and their modeling allows the encompassing of the transitions from pure noncovalent to covalent bound aggregates and to rationalize the anomalous behavior (deviations from noncovalent type interaction) pointed out in peculiar cases. Selected adducts affected by a weak chemical bond, as those promoting the formation of the intermolecular halogen bond, are also properly rationalized. The behavior of noble gas atoms excited in their long-life metastable states, showing a strongly enhanced reactivity, has been also enclosed in the present investigation.

de Almeida M, Marques JMC, Prudente FV. A thermodynamic view on the microsolvation of ions by rare gas: application to Li+ with argon

Parallel tempering Monte Carlo calculations on the Li⁺Ar_n microsolvation clusters have shown that the two peaks appearing in the heat capacity curve as a function of temperature correspond to the melting of the second and first solvation shells.

Considerable matrix shift in the electronic transitions of helium-solvated cesium dimer cation Cs2He

We investigate the photodissociation of helium-solvated cesium dimer cations using action spectroscopy and quantum chemical calculations. The spectrum of Cs2He+ shows three distinct absorption bands into both bound and dissociative states. Upon solvation with further helium atoms, considerable shifts of the absorption bands are observed, exceeding 0.1 eV (850 cm-1) already for Cs2He10+, along with significant broadening. The shifts are highly sensitive to the character of the excited state. Our calculations show that helium atoms adsorb on the ends of Cs2+. The shifts are particularly pronounced if the excited state orbitals extend to the area occupied by the helium atoms. In this case, Pauli repulsion leads to a deformation of the excited state orbitals, resulting in the observed blue shift of the transition. Since the position of the weakly bound helium atoms is ill defined, Pauli repulsion also explains the broadening.

Atomic Gold Ions Clustered with Noble Gases: Helium, Neon, Argon, Krypton, and Xenon

High-resolution mass spectra of helium droplets doped with gold and ionized by electrons reveal HenAu+ cluster ions. Additional doping with heavy noble gases results in NenAu+, ArnAu+, KrnAu+, and XenAu+ cluster ions. The high stability predicted for covalently bonded Ar2Au+, Kr2Au+, and Xe2Au+ is reflected in their relatively high abundance. Surprisingly, the abundance of Ne2Au+, which is predicted to have zero covalent bonding character and no enhanced stability, features a local maximum, too. The predicted size and structure of complete solvation shells surrounding ions with essentially nondirectional bonding depends primarily on the ratio σ* of the ion-ligand versus the ligand-ligand distance. For Au+ solvated in helium and neon, the ratio σ* is slightly below 1, favoring icosahedral packing in agreement with a maximum observed in the corresponding abundance distributions at n = 12. HenAu+ appears to adopt two additional solvation shells of Ih symmetry, containing 20 and 12 atoms, respectively. For ArnAu+, with σ* ≈ 0.67, one would expect a solvation shell of octahedral symmetry, in agreement with an enhanced ion abundance at n = 6. Another anomaly in the ion abundance at Ar9Au+ matches a local maximum in its computed dissociation energy.

Complexes of Alkali Metal Cations and Molecular Hydrogen: Potential Energy Surfaces and Bound States

Complexes between metal cations and molecular hydrogen are systems quite amenable for precise spectroscopic and theoretical studies, and at the same time, they are relevant for applications in hydrogen storage and astrochemistry. In this work, we report new intermolecular potential energy surfaces and rovibrational states calculations for complexes involving molecular hydrogen and alkaline metal cations, M⁺-H₂ (M⁺ = Na⁺, K⁺, Rb⁺, Cs⁺). The intermolecular potentials, formulated in an internally consistent way to emphasize differences in the properties of the systems, are represented by simple analytical expressions whose parameters have been optimized from comparison with accurate ab initio calculations. Properties of the low-lying bound states-binding energies, frequencies, and rotational constants-are compared with previous measurements or computations and an overall good agreement is achieved, supporting the reliability of the present formulation. Variations of these properties as a function of the cation size and isotopic substitution, with a proper sequence of ortho and para rotational levels, are also discussed.