Bosonic helium droplets with cationic impurities: onset of electrostriction and snowball effects from quantum calculations

Transitional Supersolidity in Ion Doped Helium Droplets

^{4}He nanodroplets doped with an alkali ion feature a snowball of crystallized layers surrounded by superfluid helium. For large droplets, we predict that a transitional supersolid layer can form, bridging between the solid core and the liquid bulk, where the ^{4}He density displays modulations of icosahedral group symmetry. To identify the different phases, we combine density functional theory with the semiclassical Gaussian time-dependent Hartree method for localized many-body systems. This hybrid approach can handle large particle numbers and provides insight into the physical origin of the supersolid layer. For small droplets, we verify that the predictions of our approach are in excellent agreement with Path-Integral Monte Carlo calculations.

Femtosecond-and-atom-resolved solvation dynamics of a Na+ ion in a helium nanodroplet

Recently, it was shown how the primary steps of solvation of a single Na+ ion, instantly created at the surface of a nanometer-sized droplet of liquid helium, can be followed at the atomic level [Albrechtsen et al., Nature 623, 319 (2023)]. This involved measuring, with femtosecond time resolution, the gradual attachment of individual He atoms to the Na+ ion as well as the energy dissipated from the local region of the ion. In this current work, we provide a more comprehensive and detailed description of the experimental findings of the solvation dynamics and present an improved Poisson-statistical analysis of the time-resolved yields of the Na+Hen ions recorded. For droplets containing an average of 5200 He atoms, this analysis gives a binding rate of 1.84 ± 0.09 atoms/ps for the binding of the first five He atoms to the Na+ ion. In addition, thanks to accurate theoretical values for the evaporation energies of the Na+Hen ions, obtained by path integral Monte Carlo methods using a new potential energy surface presented here for the first time, we improve the determination of the time-dependent removal of the solvation energy from the region around the sodium ion. We find that it follows Newton's law of cooling for the first 5 ps. Measurements were carried out for three different average droplet sizes, ⟨ND⟩ = 9000, 5200, and 3600 helium atoms, and differences between these results are discussed.

Microsolvation of cationic alkali dimers in helium: quantum delocalization and solid-like/liquid-like behaviors of He shells

We performed path-integral molecular dynamics (PIMD) simulations in the NVT ensemble to investigate the quantum solvation of Li2+ in He nanoclusters at a low temperature of 2 K. The interaction potentials were modeled using a sum-of-potentials approach, incorporating automated learning ab initio-based models up to three-body terms. Additionally, the semiclassical quadratic Feynman-Hibbs approach was applied to incorporate quantum effects into classical computations effectively, enabling the study of HeNLi2+ complexes with up to 50 He atoms. The quantum simulations revealed strong evidence of local solid-like behavior in the He atoms within the first solvation shell surrounding the Li2+ dimer cation. In contrast, the second and third solvation shells displayed delocalized He densities, allowing for the interchange of He atoms between these layers, indicative of a liquid-like structure. Our findings align with earlier studies of He-doped clusters, particularly in systems where the charged impurity interacts strongly with the solvent medium, significantly impacting the helium environment at the microscopic level.

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.

Observing the primary steps of ion solvation in helium droplets

Solvation is a ubiquitous phenomenon in the natural sciences. At the macroscopic level, it is well understood through thermodynamics and chemical reaction kinetics1,2. At the atomic level, the primary steps of solvation are the attraction and binding of individual molecules or atoms of a solvent to molecules or ions of a solute1. These steps have, however, never been observed in real time. Here we instantly create a single sodium ion at the surface of a liquid helium nanodroplet3,4, and measure the number of solvent atoms that successively attach to the ion as a function of time. We found that the binding dynamics of the first five helium atoms is well described by a Poissonian process with a binding rate of 2.0 atoms per picosecond. This rate is consistent with time-dependent density-functional-theory simulations of the solvation process. Furthermore, our measurements enable an estimate of the energy removed from the region around the sodium ion as a function of time, revealing that half of the total solvation energy is dissipated after four picoseconds. Our experimental method opens possibilities for benchmarking theoretical models of ion solvation and for time-resolved measurements of cation-molecule complex formation.

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.

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.

Modelling interactions of cationic dimers in He droplets: microsolvation trends in HenK2+ clusters

We report the results of a detailed theoretical investigation of small K₂⁺-doped He clusters. The structural characteristics and stabilities of such cations are determined from ab initio electronic structure calculations at the MRCI+Q level of theory. The underlying interactions show a multireference character and such effects are analyzed. The interaction potentials are constructed employing an interpolation technique within the inverse problem theory method, while the nuclear quantum effects are computed for the trimers, their spatial arrangements are discussed, and information was extracted on the orientational anisotropy of the forces. We found that energetically the most stable conformer corresponds to linear arrangements that are taking place under large amplitude vibrations, with high zero-point energy. We have further looked into the behavior of higher-order species with various He atoms surrounding the cationic dopant. By using a sum of potentials approach and an evolutionary programming method, we analyzed the structural stability of clusters with up to six He atoms in comparison with interactions energies obtained from MRCI+Q quantum chemistry computations. Structures containing He_n motifs that characterize pure rare gas clusters, appear for the larger K₂⁺-doped He clusters, showing selective growth during the microsolvation process of the alkali-dimer cation surrounded by He atoms. Such results indicate the existence of local solvation microstructures in these aggregates, where the cationic impurity could get trapped for a short time, contributing to the slow ionic mobility observed experimentally in ultra-cold He-droplets.

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

Solvation of Cs⁺ ions inside helium droplets has been investigated both experimentally and theoretically. On the one hand, mass spectra of doped helium clusters ionized with a crossed electron beam, He_NCs⁺, have been recorded for sizes up to N = 60. The analysis of the ratio between the observed peaks for each size N reveals evidences of the closure of the first solvation shell when 17 He atoms surround the alkali ion. On the other hand, we have obtained energies and geometrical structures of the title clusters by means of basin-hopping, diffusion Monte Carlo (DMC), and path integral Monte Carlo (PIMC) methods. The analytical He-Cs⁺ interaction potential employed in our calculations is represented by the improved Lennard-Jones expression optimized on high level ab initio energies. The weakness of the existing interaction between helium and Cs⁺ in comparison with some other alkali ions such as Li⁺ is found to play a crucial role. Our theoretical findings confirm that the first solvation layer is completed at N = 17 and both evaporation and second difference energies obtained with the PIMC calculation seem to reproduce a feature observed at N = 12 for the experimental ion abundance. The analysis of the DMC probability distributions reveals the important contribution from the icosahedral structure to the overall configuration for He₁₂Cs⁺.

Snowball formation for Cs+ solvation in molecular hydrogen and deuterium

Interactions of atomic cations with molecular hydrogen are of interest for a wide range of applications in hydrogen technologies. These interactions are fairly strong despite being non-covalent, hence one can ask whether hydrogen molecules would form dense, solid-like, solvation shells around the ion (snowballs) or rather a more weakly bound compound. In this work, the interactions between Cs+ and H2 are studied both experimentally and computationally. Isotopic substitution of H2 by D2 is also investigated. On the one hand, helium nanodroplets doped with cesium and hydrogen or deuterium are ionized by electron impact and the (H2/D2)nCs+ (up to n = 30) clusters formed are identified via mass spectrometry. On the other hand, a new analytical potential energy surface, based on ab initio calculations, is developed and used to study cluster energies and structures by means of classical and quantum-mechanical Monte Carlo methods. The most salient features of the measured ion abundances are remarkably mimicked by the computed evaporation energies, particularly for the clusters composed of deuterium. This result supports the reliability of the present potential energy surface and allows us to recommend its use in related systems. Clusters with either twelve H2 or D2 molecules stand out for their stability and quasi-rigid icosahedral structures. However, the first solvation shell involves thirteen or fourteen molecules for hydrogenated or deuterated clusters, respectively. This shell retains its internal structure when extra molecules are added to the second shell and is nearly solid-like, especially for the deuterated clusters. The role played by three-body induction interactions as well as the rotational degrees of freedom is analyzed and they are found to be significant (up to 15% and 18%, respectively) for the molecules belonging to the first solvation shell.

Lithium ions solvated in helium

We report on a combined experimental and theoretical study of Li+ ions solvated by up to 50 He atoms. The experiments show clear enhanced abundances associated with HenLi+ clusters where n = 2, 6, 8, and 14. We find that classical methods, e.g. basin-hopping (BH), give results that qualitatively agree with quantum mechanical methods such as path integral Monte Carlo, diffusion Monte Carlo and quantum free energy, regarding both energies and the solvation structures that are formed. The theory identifies particularly stable structures for n = 4, 6 and 8 which line up with some of the most abundant features in the experiments.

Hyperspherical coupled channel calculations of energy and structure of (4)He-(4)He-Li(+) and its isotopic combinations

The ground state vibrational energy and spatial features of (4)He-(4)He-Li(+) and its triatomic isotopic complexes are studied using the slow variable discretization (SVD) method in the hyperspherical coordinates for the zero total angular momentum. Our results show that the dominant structure of the system is an isosceles triangle with the shorter side associated with the two Li(+)-He distances using the sum-of-potential approximation. Corrections caused by the induced dipole-induced dipole interactions on the He atoms are also investigated. The effects are seen to be small and have a minor influence on the binding energy and the structure of present system. The results are also compared with the full ab initio calculations including all the three-body interactions and information of three-body corrections is obtained.

Effect of kinetic energy on the doping efficiency of cesium cations into superfluid helium droplets

We present an experimental investigation of the effect of kinetic energy on the ion doping efficiency of superfluid helium droplets using cesium cations from a thermionic emission source. The kinetic energy of Cs(+) is controlled by the bias voltage of a collection grid collinearly arranged with the droplet beam. Efficient doping from ions with kinetic energies from 20 eV up to 480 V has been observed in different sized helium droplets. The relative ion doping efficiency is determined by both the kinetic energy of the ions and the average size of the droplet beam. At a fixed source temperature, the number of doped droplets increases with increasing grid voltage, while the relative ion doping efficiency decreases. This result implies that not all ions are captured upon encountering with a sufficiently large droplet, a deviation from the near unity doping efficiency for closed shell neutral molecules. We propose that this drop in ion doping efficiency with kinetic energy is related to the limited deceleration rate inside a helium droplet. When the source temperature changes from 14 K to 17 K, the relative ion doping efficiency decreases rapidly, perhaps due to the lack of viable sized droplets. The size distribution of the Cs(+)-doped droplet beam can be measured by deflection and by energy filtering. The observed doped droplet size is about 5 × 10(6) helium atoms when the source temperature is between 14 K and 17 K.

On the size and structure of helium snowballs formed around charged atoms and clusters of noble gases

Helium nanodroplets doped with argon, krypton, or xenon are ionized by electrons and analyzed in a mass spectrometer. HenNgx(+) ions containing up to seven noble gas (Ng) atoms and dozens of helium atoms are identified; the high resolution of the mass spectrometer combined with advanced data analysis make it possible to unscramble contributions from isotopologues that have the same nominal mass but different numbers of helium or Ng atoms, such as the magic He20(84)Kr2(+) and the isobaric, nonmagic He41(84)Kr(+). Anomalies in these ion abundances reveal particularly stable ions; several intriguing patterns emerge. Perhaps most astounding are the results for HenAr(+), which show evidence for three distinct, solid-like solvation shells containing 12, 20, and 12 helium atoms. This observation runs counter to the common notion that only the first solvation shell is solid-like but agrees with calculations by Galli et al. for HenNa(+) [J. Phys. Chem. A 2011, 115, 7300] that reveal three shells of icosahedral symmetry. HenArx(+) (2 ≤ x ≤ 7) ions appear to be especially stable if they contain a total of n + x = 19 atoms. A sequence of anomalies in the abundance distribution of HenKrx(+) suggests that rings of six helium atoms are inserted into the solvation shell each time a krypton atom is added to the ionic core, from Kr(+) to Kr3(+). Previously reported strong anomalies at He12Kr2(+) and He12Kr3(+) [Kim , J. H.; et al. J. Chem. Phys. 2006, 124, 214301] are attributed to a contamination. Only minor local anomalies appear in the distributions of HenXex(+) (x ≤ 3). The distributions of HenKr(+) and HenXe(+) show strikingly similar, broad features that are absent from the distribution of HenAr(+); differences are tentatively ascribed to the very different fragmentation dynamics of these ions.

On the stability of cationic complexes of neon with helium--solving an experimental discrepancy

Helium nanodroplets are doped with neon and ionized by electrons. The size-dependence of the ion abundance of HenNex(+), identified in high-resolution mass spectra, is deduced for complexes containing up to seven neon atoms and dozens of helium atoms. Particularly stable ions are inferred from anomalies in the abundance distributions. Two pronounced anomalies at n = 11 and 13 in the HenNe(+) series confirm drift-tube data reported by Kojima et al. [T. M. Kojima et al., Z. Phys. D, 1992, 22, 645]. The discrepancy with previously published spectra of neon-doped helium droplets, which did not reveal any abundance anomalies [T. Ruchti et al., J. Chem. Phys., 1998, 109, 10679-10687; C. A. Brindle et al., J. Chem. Phys., 2005, 123, 064312], is most likely due to limited mass resolution, which precluded unambiguous analysis of contributions from different ions with identical nominal mass. However, calculated dissociation energies of HenNe(+) reported so far do not correlate with the present data, possibly because of challenges in correctly treating the linear, asymmetric [He-Ne-He](+) ionic core in HenNe(+). Anomalies identified in the distributions of HenNex(+) for x > 1, including prominent ones at He12Ne2(+) and He14Ne2(+), may help to better understand solvation of Ne(+) and Nex(+) in helium.

Solvation of Na+, K+, and their dimers in helium

Helium atoms bind strongly to alkali cations which, when embedded in liquid helium, form so-called snowballs. Calculations suggest that helium atoms in the first solvation layer of these snowballs form rigid structures and that their number (n) is well defined, especially for the lighter alkalis. However, experiments have so far failed to accurately determine values of n. We present high-resolution mass spectra of Na⁺He_n, K⁺He_n, Na₂⁺He_n and K₂⁺He_n, formed by electron ionization of doped helium droplets; the data allow for a critical comparison with several theoretical studies. For sodium and potassium monomers the spectra indicate that the value of n is slightly smaller than calculated. Na₂⁺He_n displays two distinct anomalies at n=2 and n=6, in agreement with theory; dissociation energies derived from experiment closely track theoretical values. K₂⁺He_n distributions are fairly featureless, which also agrees with predictions.

Higher order diffusion Monte Carlo propagators for linear rotors as diffusion on a sphere: development and application to O2@He(n)

Exploiting the theoretical treatment of particles diffusing on corrugated surfaces and the isomorphism between the "particle on a sphere" and a linear molecule rotation, a new diffusion kernel is introduced to increase the order of diffusion Monte Carlo (DMC) simulations involving linear rotors. Tests carried out on model systems indicate the superior performances of the new rotational diffusion kernel with respect to the simpler alternatives previously employed. In particular, it is evidenced a second order convergence toward exact results with respect to the time step of dynamical correlation functions, a fact that guarantees an identical order for the diffusion part of the DMC projector. The algorithmic advantages afforded by the latter are discussed, especially with respect to the "a posteriori" and "on the fly" extrapolation schemes. As a first application to the new algorithm, the structure and energetics of O(2)@He(n) (n = 1-40) clusters have been studied. This was done to investigate the possible cause of the quenching of the reaction between O(2) and Mg witnessed upon increasing the size of superfluid He droplets used as a solvent. With the simulations on O(2) indicating a strong localization in the cluster core, the behaviour as a function of n is ascribed to the extremely fluxional comportment of Mg@He(n), which dwells far from the droplet center, albeit being solvated, when n is large.

Path integral Monte Carlo study of 4He clusters doped with alkali and alkali-earth ions

Path integral Monte Carlo calculations of (4)He nanodroplets doped with alkali (Na(+), K(+) and Cs(+)) and alkali-earth (Be(+) and Mg(+)) ions are presented. We study the system at T = 1 K and between 14 and 128 (4)He atoms. For all studied systems, we find that the ion is well localized at the center of the droplet with the formation of a "snowball" of well-defined shells of localized (4)He atoms forming solid-like order in at least the first surrounding shell. The number of surrounding helium shells (two or three) and the number of atoms per shell and the degree of localization of the helium atoms are sensitive to the type of ion. The number of (4)He atoms in the first shell varies from 12 for Na(+) to 18 for Mg(+) and depends weakly on the size of the droplet. The study of the density profile and of the angular correlations shows that the local solid-like order is more pronounced for the alkali ions with Na(+) giving a very stable icosahedral order extending up to three shells.

Water as a solute: competitive shell formation in (He,Ne) mixed microdroplets

Quantum Monte Carlo (QMC) stochastic calculations are carried out for a series of mixed rare gas clusters containing He and Ne which further include one H(2)O molecule as a single dopant. The ab initio potentials employed in the calculations, and the structural details provided by the QMC results, clearly reveal that the differences in the interaction forces which exist between the two solvent adatoms and the molecular solute are causing strongly competing environments that generate preferential shell structuring when surrounding the water molecule. The different behaviour of the two solvents, and the energetics of mixing, are analyzed in detail for small aggregates and for larger mixtures, revealing structural effects which originate from the different networking between solvent adatoms.

Review: gas-phase ion chemistry of the noble gases: recent advances and future perspectives

This review article surveys recent experimental and theoretical advances in the gas-phase ion chemistry of the noble gases. Covered issues include the interaction of the noble gases with metal and non-metal cations, the conceivable existence of covalent noble-gas anions, the occurrence of ion-molecule reactions involving singly-charged xenon cations, and the occurrence of bond-forming reactions involving doubly-charged cations. Research themes are also highlighted, that are expected to attract further interest in the future.

Forming Rb(+) snowballs in the center of He nanodroplets

Helium nanodroplets doped with rubidium atoms are ionized by applying a resonant two-step ionization scheme. Subsequent immersion of rubidium ions is observed in time-of-flight mass spectra. While alkali-metal atoms usually desorb from the surface of a helium nanodroplet upon electronic excitation, rubidium in its excited 5(2)P(1/2) state provides an exception from this rule (Auböck et al., Phys. Rev. Lett., 2008, 101, 35301). In our new experiment, Rb atoms are selectively excited either to the 5(2)P(1/2) or to the 5(2)P(3/2) state. From there they are ionized by a laser pulse. Time-of-flight mass spectra of the ionization products reveal that the intermediate population of the 5(2)P(1/2) state does not only make the ionization process Rb-monomer selective, but also gives rise to a very high yield of Rb(+)-He(N) complexes. Ions with masses of up to several thousand amu have been monitored, which can be explained by an immersion of the single Rb ion into the He nanodroplet, where most likely a snowball is formed in the center of the He nanodroplet. As the most stable position for an ion is in the center of a He nanodroplet, our results agree well with theory.

Snowballs, quantum solvation and coordination: lead ions inside small helium droplets

Ab initio calculations are used to construct an analytical many-body potential for Pb(2+)He(n) and Pb(+)He(n) clusters which accounts for non pairwise additive interactions. The potential surface reproduces the global minima for cluster sizes ranging from n = 1 to n = 16 obtained from explicit ab initio calculations and found in a previous search for ultrahigh coordination numbers. Ground state energies and structures obtained by accurate diffusion quantum Monte Carlo calculations are used to investigate if quantum effects qualitatively affect the formation of coordination shells. For Pb(2+) doped clusters a first solvation shell is closed at n = 12 and gradually softened by additional helium atoms which start to form a distinct second shell only at n = 16. Spin-orbit coupling profoundly influences the structure of Pb(+)He(n) clusters and causes a gradual structural evolution without pronounced solvation shells.

Alkali-helium snowball complexes formed on helium nanodroplets

We systematically investigate the formation and stability of snowballs formed by femtosecond photoionization of small alkali clusters bound to helium nanodroplets. For all studied alkali species Ak = (Na,K,Rb,Cs) we observe the formation of snowballs Ak(+)He(N) when multiply doping the droplets. Fragmentation of clusters Ak(N) upon ionization appears to enhance snowball formation. In the case of Na and Cs we also detect snowballs Ak(2) (+)He(N) formed around Ak dimer ions. While the snowball progression for Na and K is limited to less than 11 helium atoms, the heavier atoms Rb and Cs feature wide distributions at least up to Ak(+)He(41). Characteristic steps in the mass spectra of Cs-doped helium droplets are found at positions consistent with predictions on the closure of the first shell of helium atoms around the Ak(+) ion based on variational Monte Carlo simulations.

On the size of ions solvated in helium clusters

Helium nanodroplets are doped with SF(6), C(4)F(8), CCl(4), C(6)H(5)Br, CH(3)I, and I(2). Upon interaction with free electrons a variety of positively and negatively charged cluster ions X(+/-)He(n) are observed where X(+/-) = F(+/-), Cl(+/-), Br(+/-), I(+), I(2) (+), or CH(3)I(+). The yield of these ions versus cluster size n drops at characteristic sizes n(s) that range from n(s) = 10.2+/-0.6 for F(+) to n(s) = 22.2+/-0.2 for Br(-). n(s) values for halide anions are about 70% larger than for the corresponding cations. The steps in the ion yield suggest closure of the first solvation shell. We propose a simple classical model to estimate ionic radii from n(s). Assuming the helium density in the first solvation shell equals the helium bulk density one finds that radii of halide anions in helium are nearly twice as large as in alkali halide crystals, indicating the formation of an anion bubble due to the repulsive forces that derive from the exchange interaction. In spite of the simplicity of our model, anion radii derived from it agree within approximately 10% with values derived from the mobility of halide anions in superfluid bulk helium, and with values computed by quantum Monte Carlo methods for X(-)He(n) cluster anions.

Spin-driven structural effects in alkali doped (4)He clusters from quantum calculations

In this paper, we carry out variational Monte Carlo and diffusion Monte Carlo (DMC) calculations for Li(2)((1)Sigma(g) (+))((4)He)(N) and Li(2)((3)Sigma(u) (+))((4)He)(N) with N up to 30 and discuss in detail the results of our computations. After a comparison between our DMC energies with the "exact" discrete variable representation values for the species with one (4)He, in order to test the quality of our computations at 0 K, we analyze the structural features of the whole range of doped clusters. We find that both species reside on the droplet surface, but that their orientation is spin driven, i.e., the singlet molecule is perpendicular and the triplet one is parallel to the droplet's surface. We have also computed quantum vibrational relaxation rates for both dimers in collision with a single (4)He and we find them to differ by orders of magnitude at the estimated surface temperature. Our results therefore confirm the findings from a great number of experimental data present in the current literature and provide one of the first attempts at giving an accurate, fully quantum picture for the nanoscopic properties of alkali dimers in (4)He clusters.

Chemical solutions in a quantum solvent: anionic electrolytes in 4He nanodroplets

Variational and diffusion Monte Carlo (VMC and DMC) calculations are presented for anionic electrolytes solvated in (4)He. The electrolytes have the general structure X(-)(He)(N), with X=F, Cl, Br and I, and N varying up to 40 (41 for I(-)). The overall interaction potential is obtained from accurate ab initio data for the two-body components and then using the sum-of-potentials approximation. Our computational scheme is a robust procedure, giving us accurate trial wavefunctions that can be used to perform high-quality DMC calculations. The results indicate very marked delocalization and permanence of the liquid-like quantum features of the solvent adatoms surrounding the anionic impurities. This finding stands in contrast to the more structured, solid-like behavior of the quantum solutions with alkali metal cations embedded in He nanodroplets. While other negatively charged species such as H(-) have shown an overall repulsive interaction with He, the present calculations clearly indicate that the halogen anions remain solvated within liquid-like solvent "bubbles" of species-dependent size.