Solvation of ions in helium

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.

Electronic and vibrational spectroscopies of aromatic clusters with He in a supersonic jet: The case of neutral and cationic phenol-Hen (n = 1 and 2)

Van der Waals clusters composed of He and aromatic molecules provide fundamental information about intermolecular interactions in weakly bound systems. In this study, phenol-helium clusters (PhOH-Hen with n ≤ 2) are characterized for the first time by UV and IR spectroscopies. The S1 ← S0 origin and ionization energy both show small but additive shifts, suggesting π-bound structures of these clusters, a conclusion supported by rotational contour analyses of the S1 origin bands. The OH stretching vibrations of the PhOH moiety in the clusters match with those of bare PhOH in both the S0 and D0 states, illustrating the negligible perturbation of the He atoms on the molecular vibration. Matrix shifts induced by He attachment are discussed based on the observed band positions with the help of complementary quantum chemical calculations. For comparison, the UV and ionization spectra of PhOH-Ne are reported as well.

Probing the binding and activation of small molecules by gas-phase transition metal clusters via IR spectroscopy

Isolated transition metal clusters have been established as useful models for extended metal surfaces or deposited metal particles, to improve the understanding of their surface chemistry and of catalytic reactions. For this objective, an important milestone has been the development of experimental methods for the size-specific structural characterization of clusters and cluster complexes in the gas phase. This review focusses on the characterization of molecular ligands, their binding and activation by small transition metal clusters, using cluster-size specific infrared action spectroscopy. A comprehensive overview and a critical discussion of the experimental data available to date is provided, reaching from the initial results obtained using line-tuneable CO₂ lasers to present-day studies applying infrared free electron lasers as well as other intense and broadly tuneable IR laser sources.

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.

Structure and formation of copper cluster ions in multiply charged He nanodroplets

The structure of cationic and anionic Cu clusters grown in multiply charged superfluid He nanodroplets was investigated using He tagging as a chemical probe. Further, the structure assignment was done based on the magic-numbered ions, representing the most energetically favorable structures. The exact geometry of the cluster and positions of He is verified by calculations. It was found that the structure of the clusters grown in the He droplets is similar to that produced with a laser ablation source and the lowest energy structures predicted by theoretical investigations. The only difference is the structure of the Cu₅⁺, which in our experiments has a twisted-X geometry, rather than a bipyramid or planar half-wheel geometry suggested by previous studies. This might be attributed to the different cluster formation mechanisms, the absence of the Ar-tag and the ultracold environment. It was also found that He tends to bind to partially more electro-negative or positive areas of the anionic or cationic clusters, respectively.

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.

Solvation of Large Polycyclic Aromatic Hydrocarbons in Helium: Cationic and Anionic Hexabenzocoronene

The adsorption of helium on charged hexabenzocoronene (Hbc, C₄₂H₁₈), a planar polycyclic aromatic hydrocarbon (PAH) molecule of D_6h symmetry, was investigated by a combination of high-resolution mass spectrometry and classical and quantum computational methods. The ion abundance of He_nHbc⁺ complexes versus size n features prominent local anomalies at n = 14, 38, 68, 82, and a weak one at 26, indicating that for these “magic” sizes, the helium evaporation energies are relatively large. Surprisingly, the mass spectra of anionic He_nHbc⁻ complexes feature a different set of anomalies, namely at n = 14, 26, 60, and 62, suggesting that the preferred arrangement of the adsorbate atoms depends on the charge of the substrate. The results of our quantum calculations show that the adsorbate layer grows by successive filling of concentric rings that surround the central benzene ring, which is occupied by one helium atom each on either side of the substrate. The helium atoms are fairly localized in filled rings and they approximately preserve the D_6h symmetry of the substrate, but helium atoms in partially filled rings are rather delocalized. The first three rings contain six atoms each; they account for magic numbers at n = 14, 26, and 38. The size of the first ring shrinks as atoms are filled into the second ring, and the position of atoms in the second ring changes from hollow sites to bridge sites as atoms are filled into the third ring. Beyond n = 38, however, the arrangement of helium atoms in the first three rings remains essentially frozen. Presumably, another ring is filled at n = 68 for cations and n = 62 for anions. The calculated structures and energies do not account for the difference between charge states, although they agree with the measurements for the cations and show that the first solvation shell of Hbc^± is complete at n = 68. Beyond that size, the adsorbate layer becomes three-dimensional, and the circular arrangement of helium changes to hexagonal.

Adsorption of Helium and Hydrogen on Triphenylene and 1,3,5-Triphenylbenzene

The adsorption of helium or hydrogen on cationic triphenylene (TPL, C18H12), a planar polycyclic aromatic hydrocarbon (PAH) molecule, and of helium on cationic 1,3,5-triphenylbenzene (TPB, C24H18), a propeller-shaped PAH, is studied by a combination of high-resolution mass spectrometry and classical and quantum computational methods. Mass spectra indicate that HenTPL+ complexes are particularly stable if n = 2 or 6, in good agreement with the quantum calculations that show that for these sizes, the helium atoms are strongly localized on either side of the central carbon ring for n = 2 and on either side of the three outer rings for n = 6. Theory suggests that He14TPL+ is also particularly stable, with the helium atoms strongly localized on either side of the central and outer rings plus the vacancies between the outer rings. For HenTPB+, the mass spectra hint at enhanced stability for n = 2, 4 and, possibly, 11. Here, the agreement with theory is less satisfactory, probably because TPB+ is a highly fluxional molecule. In the global energy minimum, the phenyl groups are rotated in the same direction, but when the zero-point harmonic correction is included, a structure with one phenyl group being rotated opposite to the other two becomes lower in energy. The energy barrier between the two isomers is very small, and TPB+ could be in a mixture of symmetric and antisymmetric states, or possibly even vibrationally delocalized.

Chemistry and physics of dopants embedded in helium droplets

Helium droplets represent a cold inert matrix, free of walls with outstanding properties to grow complexes and clusters at conditions that are perfect to simulate cold and dense regions of the interstellar medium. At sub-Kelvin temperatures, barrierless reactions triggered by radicals or ions have been observed and studied by optical spectroscopy and mass spectrometry. The present review summarizes developments of experimental techniques and methods and recent results they enabled.

Mass Spectra Resulting from Collision Processes

A new database and viewer for mass spectra resulting from collision processes is presented that follows the standards of the Virtual Atomic and Molecular Data Centre (VAMDC). A focus was placed on machine read and write access, as well as ease of use. In a browser-based viewer, mass spectra and all parameters related to a given measurement can be shown. The program additionally enables a direct comparison between two mass spectra, either by plotting them on top of each other or their difference to identify subtle variations in the data.

Microsolvation of H2O+, H3O+, and CH3OH2+ by He in a cryogenic ion trap: structure of solvation shells

Due to the weak interactions of He atoms with neutral molecules and ions, the preparation of size-selected clusters for the spectroscopic characterization of their structures, energies, and large amplitude motions is a challenging task. Herein, we generate H₂O⁺He_n (n ≤ 9) and H₃O⁺He_n (n ≤ 5) clusters by stepwise addition of He atoms to mass-selected ions stored in a cryogenic 22-pole ion trap held at 5 K. The population of the clusters as a function of n provides insight into the structure of the first He solvation shell around these ions given by the anisotropy of the cation-He interaction potential. To rationalize the observed cluster size distributions, the structural, energetic, and vibrational properties of the clusters are characterized by ab initio calculations up to the CCSD(T)/aug-cc-pVTZ level. The cluster growth around both the open-shell H₂O⁺ and closed-shell H₃O⁺ ions begins by forming nearly linear and equivalent OH⋯He hydrogen bonds (H-bonds) leading to symmetric structures. The strength of these H-bonds decreases slightly with n due to noncooperative three-body induction forces and is weaker for H₃O⁺ than for H₂O⁺ due to both enhanced charge delocalization and reduced acidity of the OH protons. After filling all available H-bonded sites, addition of further He ligands around H₂O⁺ (n = 3-4) occurs at the electrophilic singly occupied 2p_z orbital of O leading to O⋯He p-bonds stabilized by induction and small charge transfer from H₂O⁺ to He. As this orbital is filled for H₃O⁺, He atoms occupy in the n = 4-6 clusters positions between the H-bonded He atoms, leading to a slightly distorted regular hexagon ring for n = 6. Comparison between H₃O⁺He_n and CH₃OH₂⁺He_n illustrates that CH₃ substitution substantially reduces the acidity of the OH protons, so that only clusters up to n = 2 can be observed. The structure of the solvation sub-shells is visible in both the binding energies and the predicted vibrational OH stretch and bend frequencies.

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.

Helium Nanodroplet Infrared Action Spectroscopy of the Proton-Bound Dimer of Hydrogen Sulfate and Formate: Examining Nuclear Quantum Effects

The proton-bound dimer of hydrogen sulfate and formate is an archetypal structure for ionic hydrogen-bonding complexes that contribute to biogenic aerosol nucleation. Of central importance for the structure and properties of this complex is the location of the bridging proton connecting the two conjugate base moieties. The potential energy surface for bridging proton translocation features two local minima, with the proton localized at either the formate or hydrogen sulfate moiety. However, electronic structure methods reveal a shallow potential energy surface governing proton translocation, with a barrier on the order of the zero-point energy. This shallow potential complicates structural assignment and necessitates a consideration of nuclear quantum effects. In this work, we probe the structure of this complex and its isotopologues, utilizing infrared (IR) action spectroscopy of ions captured in helium nanodroplets. The IR spectra indicate a structure in which a proton is shared between the hydrogen sulfate and formate moieties, HSO₄^-···H⁺···^-OOCH. However, because of the nuclear quantum effects and vibrational anharmonicities associated with the shallow potential for proton translocation, the extent of proton displacement from the formate moiety remains unclear, requiring further experiments or more advanced theoretical treatments for additional insight.

Adsorption of Helium on Small Cationic PAHs: Influence of Hydrocarbon Structure on the Microsolvation Pattern

The adsorption of up to ∼100 helium atoms on cations of the planar polycyclic aromatic hydrocarbons (PAHs) anthracene, phenanthrene, fluoranthene, and pyrene was studied by combining helium nanodroplet mass spectrometry with classical and quantum computational methods. Recorded time-of-flight mass spectra reveal a unique set of structural features in the ion abundance as a function of the number of attached helium atoms for each of the investigated PAHs. Path-integral molecular dynamics simulations were used with a polarizable potential to determine the underlying adsorption patterns of helium around the studied PAH cations and in good general agreement with the experimental data. The calculated structures of the helium-PAH complexes indicate that the arrangement of adsorbed helium atoms is highly sensitive toward the structure of the solvated PAH cation. Closures of the first solvation shell around the studied PAH cations are suggested to lie between 29 and 37 adsorbed helium atoms depending on the specific PAH cation. Helium atoms are found to preferentially adsorb on these PAHs following the 3 × 3 commensurate pattern common for graphitic surfaces, in contrast to larger carbonaceous molecules like corannulene, coronene, and fullerenes that exhibit a 1 × 1 commensurate phase.

Size and Velocity Distribution of Negatively Charged Helium Nanodroplets

Precharged helium nanodroplets can be used in doping experiments with the advantage that they are amenable to size selection with electrostatic fields, therefore adding a useful tuning parameter for dopant growth. For all these applications, the knowledge of the size distribution of charged droplets is an essential parameter, which we have so far assumed would be equivalent to that of their neutral precursors. Here, this assumption is experimentally investigated for negatively charged clusters for temperatures between 4 and 9 K at a stagnation pressure of 2 MPa. We observe a dependency of the velocity of the droplets on mass per charge, especially at the lowest temperatures of the investigated range, and values 20% lower than those known from the literature. Below 6 K, a large deviation from the literature is also found for the average droplet sizes. This information has to be taken into consideration in future experiments where large, charged droplets are sought to produce large dopant clusters. Possible origins for this deviation are discussed in the text.

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.

Complexes with Atomic Gold Ions: Efficient Bis-Ligand Formation

Complexes of atomic gold with a variety of ligands have been formed by passing helium nanodroplets (HNDs) through two pickup cells containing gold vapor and the vapor of another dopant, namely a rare gas, a diatomic molecule (H₂, N₂, O₂, I₂, P₂), or various polyatomic molecules (H₂O, CO₂, SF₆, C₆H₆, adamantane, imidazole, dicyclopentadiene, and fullerene). The doped HNDs were irradiated by electrons; ensuing cations were identified in a high-resolution mass spectrometer. Anions were detected for benzene, dicyclopentadiene, and fullerene. For most ligands L, the abundance distribution of AuL_n⁺ versus size n displays a remarkable enhancement at n = 2. The propensity towards bis-ligand formation is attributed to the formation of covalent bonds in Au⁺L₂ which adopt a dumbbell structure, L-Au⁺-L, as previously found for L = Xe and C₆₀. Another interesting observation is the effect of gold on the degree of ionization-induced intramolecular fragmentation. For most systems gold enhances the fragmentation, i.e., intramolecular fragmentation in AuL_n⁺ is larger than in pure L_n⁺. Hydrogen, on the other hand, behaves differently, as intramolecular fragmentation in Au(H₂)_n⁺ is weaker than in pure (H₂)_n⁺ by an order of magnitude.

Rotation of CH3+ Cations in Helium Droplets

Observation of the free rotation of molecules in helium droplets enabled microscopic study of interaction of quantum rotors with a superfluid environment at T = 0.4 K. This work extends studies of rotation in helium to molecular cations, such as methenium, CH₃⁺. The spectrum of the v₃ band of CH₃⁺ around 3130 cm^-1 has three prominent peaks assigned to the rotational structure of the band. While the free CH₃⁺ is an oblate top, in helium it behaves as a prolate top. This effect is ascribed to the strong binding of two He atoms along the figure axis of the ion. Our results indicate that the other He atoms within the first solvation shell remain fluxional and in disparity with the widely accepted model of a rigid He "snowball" surrounding ions.

SF6+: Stabilizing Transient Ions in Helium Nanodroplets

There are myriad ions that are deemed too short-lived to be experimentally accessible. One of them is SF₆⁺. It has never been observed, although not for lack of trying. We demonstrate that long-lived SF₆⁺ can be formed by doping charged helium nanodroplets (HNDs) with sulfur hexafluoride; excess helium is then gently stripped from the doped HNDs by collisions with helium gas. The ion is identified by high-resolution mass spectrometry (resolution m/Δm = 15000), the close agreement between the expected and observed yield of ions that contain minor sulfur isotopes, and collision-induced dissociation in which mass-selected He_nSF₆⁺ ions collide with helium gas. Under optimized conditions, the yield of SF₆⁺ exceeds that of SF₅⁺. The procedure is versatile and suitable for stabilizing many other transient molecular ions.

Bimodal velocity and size distributions of pulsed superfluid helium droplet beams

We report detailed measurements of velocities and sizes of superfluid helium droplets produced from an Even-Lavie pulse valve at stagnation pressures of 20-60 atm and temperatures between 5.7 and 18.0 K. By doping neutral droplets with Rhodamine 6G cations produced from an electrospray ionization source and detecting the positively charged droplets at two different locations along the beam path, we determine the velocities of the different groups of droplets. By subjecting the doped droplet beam to a retardation field, size distributions can then be analyzed. We discover that at stagnation temperatures above 8.0 K, a single group of droplets is observed at both locations, but at 8.0 K and below, two different groups of droplets with different velocities are detectable. The slower group, considered from fragmentation of liquid helium, cannot be deterred by the retardation voltage at 9 kV, implying an exceedingly large size. The faster group, considered from condensation of gaseous helium, has a bimodal distribution when the stagnation temperatures are below 12.3 K at 20 and 40 atm, or 16.1 K at 60 atm. We also report similar size measurements using low energy electrons for impact ionization, and this latter method can be used for facile in situ characterization of pulsed droplet beams. The mechanism of the bimodal size distribution of the condensation group and the reason for the coexistence of both the condensation and fragmentation groups remain elusive.

Electronic transitions in Rb2+ dimers solvated in helium

We have measured depletion spectra of the heteronuclear (85Rb87Rb+) dimer cation complexed with up to 10 He atoms. Two absorption bands are observed between 920 and 250 nm. The transition into the repulsive 12Σu+state of HeRb2+gives rise to a broad feature at 790 nm (12,650 cm−1); it exhibits a blueshift of 98 cm−1per added He atom. The transition into the bound 12Πustate of HeRb2+reveals vibrational structure with a band head at ≤ 15,522 cm−1, a harmonic constant of 26 cm−1, and a spin–orbit splitting of ≤ 183 cm−1. The band experiences an average redshift of − 38 cm−1per added He atom. Ab initio calculations rationalize the shape of the spectra and spectral shifts with respect to the number of helium atoms attached. For a higher number of solvating helium atoms, symmetric solvation on both ends of the Rb2+ion is predicted.