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.

Generation of Large Vortex-Free Superfluid Helium Nanodroplets

Superfluid helium nanodroplets are an ideal environment for the formation of metastable, self-organized dopant nanostructures. However, the presence of vortices often hinders their formation. Here, we demonstrate the generation of vortex-free helium nanodroplets and explore the size range in which they can be produced. From x-ray diffraction images of xenon-doped droplets, we identify that single compact structures, assigned to vortex-free aggregation, prevail up to 10^{8} atoms per droplet. This finding builds the basis for exploring the assembly of far-from-equilibrium nanostructures at low temperatures.

Formation of heterogeneous clusters in superfluid helium nanodroplets: phthalocyanine and water

Clusters consisting of a single phthalocyanine molecule and a single water molecule are investigated by means of electronic spectroscopy in superfluid helium droplets. A recent spectroscopic study of those clusters [J. Fischer, F. Schlaghaufer, E.-M. Lottner, A. Slenczka, L. Christiansen, H. Stapelfeldt, M. Karra, B. Friedrich, T. Mullan, M. Schütz and D. Usvyat, J. Phys. Chem. A, 2019, 123, 10057-10064] which all exhibit a water induced electronic shift to the red is now complemented by the corresponding clusters exhibiting a water induced shift to the blue. These clusters will be analyzed by means of fluorescence excitation spectra, dispersed emission spectra, and additional experimental observations made feasible by helium droplets as cryogenic reactor. Together with the blue shifted clusters a total number of at least 6 isomeric variants could be identified in helium droplets. This contrasts to a number of only three isomeric variants obtained from quantum chemical calculations [J. Fischer, F. Schlaghaufer, E.-M. Lottner, A. Slenczka, L. Christiansen, H. Stapelfeldt, M. Karra, B. Friedrich, T. Mullan, M. Schütz and D. Usvyat, J. Phys. Chem. A, 2019, 123, 10057-10064] disregarding the helium environment and to a single isomer identified in a molecular beam experiment [J. Menapace and E. Bernstein, J. Chem. Phys., 1987, 87, 6877-6889]. The discrepancy in the number of isomers provides evidence of a profound involvement of helium in clustering. Moreover, the discrepancies between the gas phase experiment and quantum chemical calculations similarly reveal the influence of the dynamics of cluster formation on the population of global and local minima that are accessible as isomeric variants.

Clustering, collision, and relaxation dynamics in pure and doped helium nanoclusters: Density- vs particle-based approaches

The clustering, collision, and relaxation dynamics of pristine and doped helium nanodroplets is theoretically investigated in cases of pickup and clustering of heliophilic argon, collision of heliophobic cesium atoms, and coalescence of two droplets brought into contact by their mutual long-range van der Waals interaction. Three approaches are used and compared with each other. The He time-dependent density functional theory method considers the droplet as a continuous medium and accounts for its superfluid character. The ring-polymer molecular dynamics method uses a path-integral description of nuclear motion and incorporates zero-point delocalization while bosonic exchange effects are ignored. Finally, the zero-point averaged dynamics approach is a mixed quantum–classical method in which quantum delocalization is described by attaching a frozen wavefunction to each He atom, equivalent to classical dynamics with effective interaction potentials. All three methods predict that the growth of argon clusters is significantly hindered by the helium host droplet due to the impeding shell structure around the dopants and kinematic effects freezing the growing cluster in metastable configurations. The effects of superfluidity are qualitatively manifested by different collision dynamics of the heliophilic atom at high velocities, as well as quadrupole oscillations that are not seen with particle-based methods, for droplets experiencing a collision with cesium atoms or merging with each other.

Aggregation of solutes in bosonic versus fermionic quantum fluids

Quantum fluid droplets made of helium-3 (³He) or helium-4 (⁴He) isotopes have long been considered as ideal cryogenic nanolabs, enabling unique ultracold chemistry and spectroscopy applications. The droplets were believed to provide a homogeneous environment in which dopant atoms and molecules could move and react almost as in free space but at temperatures close to absolute zero. Here, we report ultrafast x-ray diffraction experiments on xenon-doped ³He and ⁴He nanodroplets, demonstrating that the unavoidable rotational excitation of isolated droplets leads to highly anisotropic and inhomogeneous interactions between the host matrix and enclosed dopants. Superfluid ⁴He droplets are laced with quantum vortices that trap the embedded particles, leading to the formation of filament-shaped clusters. In comparison, dopants in ³He droplets gather in diffuse, ring-shaped structures along the equator. The shapes of droplets carrying filaments or rings are direct evidence that rotational excitation is the root cause for the inhomogeneous dopant distributions.

Vibrational energy relaxation of a diatomic molecule in a superfluid helium nanodroplet: influence of the nanodroplet size, interaction energy and energy gap

The influence of the nanodroplet size, molecule-helium interaction potential energy and ν = 1 - ν = 0 vibrational energy gap on the vibrational energy relaxation (VER) of a diatomic molecule (X₂) in a superfluid helium nanodroplet [HeND or (⁴He)_N; finite quantum solvent at T = 0.37 K] has been studied using a hybrid quantum approach recently proposed by us and taking as a reference the VER results on the I₂@(⁴He)₁₀₀ doped nanodroplet (Vilà et al., Phys. Chem. Chem. Phys., 2018, 20, 118, which corresponds to the first theoretical study on the VER of molecules embedded in a HeND). This has allowed us to obtain a deeper insight into the vibrational relaxation dynamics. The nanodroplet size has a very small effect on the VER, as this process mainly depends on the interaction between the molecule and the nanodroplet first solvation shell. Regarding the interaction potential energy and the energy gap, both factors play an important and comparable role in the VER time properties (global relaxation time, lifetime and transition time). As the former becomes stronger the relaxation time properties decrease in a significant way (their inverse follows a linear dependence with respect to the ν = 1 - ν = 0 coupling term) and they also decrease in a significant manner when the energy gap diminishes (linear dependence on the ν = 1 - ν = 0 energy difference). We expect that this study will motivate further work on the vibrational relaxation process in HeNDs.

A Path Integral Molecular Dynamics Simulation of a Harpoon-Type Redox Reaction in a Helium Nanodroplet

We present path integral molecular dynamics (PIMD) calculations of an electron transfer from a heliophobic Cs2 dimer in its (3Σu) state, located on the surface of a He droplet, to a heliophilic, fully immersed C60 molecule. Supported by electron ionization mass spectroscopy measurements (Renzler et al., J. Chem. Phys.2016, 145, 181101), this spatially quenched reaction was characterized as a harpoon-type or long-range electron transfer in a previous high-level ab initio study (de Lara-Castells et al., J. Phys. Chem. Lett.2017, 8, 4284). To go beyond the static approach, classical and quantum PIMD simulations are performed at 2 K, slightly below the critical temperature for helium superfluidity (2.172 K). Calculations are executed in the NVT ensemble as well as the NVE ensemble to provide insights into real-time dynamics. A droplet size of 2090 atoms is assumed to study the impact of spatial hindrance on reactivity. By changing the number of beads in the PIMD simulations, the impact of quantization can be studied in greater detail and without an implicit assumption of superfluidity. We find that the reaction probability increases with higher levels of quantization. Our findings confirm earlier, static predictions of a rotational motion of the Cs2 dimer upon reacting with the fullerene, involving a substantial displacement of helium. However, it also raises the new question of whether the interacting species are driven out-of-equilibrium after impurity uptake, since reactivity is strongly quenched if a full thermal equilibration is assumed. More generally, our work points towards a novel mechanism for long-range electron transfer through an interplay between nuclear quantum delocalization within the confining medium and delocalized electronic dispersion forces acting on the two reactants.

Metal clusters synthesized in helium droplets: structure and dynamics from experiment and theory

Metal clusters have drawn continuous interest because of their high potential for the assembly of matter with special properties that may significantly differ from the corresponding bulk. Controlled combination of particular elements in one nanoparticle can increase the options for the creation of new materials for photonic, catalytic, or electronic applications. Superfluid helium droplets provide confinement and ultralow temperature, i.e. an ideal environment for the atom-by-atom aggregation of a new nanoparticle. This perspective presents a review of the current research progress on the synthesis of tailored metal and metal oxide clusters including core-shell designs, their characterization within the helium droplet beam, deposition on various solid substrates, and analysis via surface diagnostics. Special attention is given to the thermal properties of mixed metal clusters and questions about alloy formation on the nanoscale. Experimental results are accompanied by theoretical approaches employing computational chemistry, molecular dynamics simulations and He density functional theory.

Fragmentation dynamics of Ar4He1000 upon electron impact ionization: Competition between ion ejection and trapping

The fragmentation upon electron impact ionization of Ar₄He₁₀₀₀ is investigated by means of mixed quantum-classical dynamics simulations. The Ar₄ ⁺ dopant dynamics is described by a surface hopping method coupled with a diatomics-in-molecules model to properly take into account the multiple Ar₄ ⁺ electronic surfaces and possible transitions between them. Helium atoms are treated individually using zero-point averaged dynamics, a method based on the building of an effective He-He potential. Fast electronic relaxation is observed from less than 2 ps to ∼30 ps, depending on initial conditions. The main fragments observed are Ar₂ ⁺He_q and Ar₃ ⁺He_q (q ≤ 1000), with a strong contribution of the bare Ar₂ ⁺ ion, and neither Ar⁺ nor Ar⁺He_q fragments are found. The smaller fragments (q ≤ 50) are found to mostly come from ion ejection, whereas larger fragments (q > 500) originate from long-term ion trapping. Although the structure of the trapped Ar₂ ⁺ ions is the same as in the gas phase, trapped Ar₃ ⁺ and Ar₄ ⁺ are rather slightly bound Ar₂ ⁺⋯Ar and Ar₂ ⁺⋯Ar⋯Ar structures (i.e., an Ar₂ ⁺ core with one or two argon atoms roaming within the droplet). These loose structures can undergo geminate recombination and release Ar₃ ⁺He_q or Ar₄ ⁺He_q (q ≤ 50) in the gas phase and/or induce strong helium droplet evaporation. Finally, the translational energy of the fragment center of mass was found to be suitable to provide a clear signature of the broad variety of processes at play in our simulations.

Alkali atoms attached to vortex-hosting helium nanodroplets

Light absorption or fluorescence excitation spectroscopy of alkali atoms attached to ⁴He droplets is investigated as a possible way for detecting the presence of vortices. To this end, we have calculated the equilibrium configuration and energetics of alkali atoms attached to a ⁴He₁₀₀₀ droplet hosting a vortex line using ⁴He density functional theory. We use them to study how the dipole absorption spectrum of the alkali atom is modified when the impurity is attached to a vortex line. Spectra are found to be blue-shifted (higher frequencies) and broadened compared to vortex-free droplets because the dimple in which the alkali atom sits at the intersection of the vortex line and the droplet surface is deeper. This effect is smaller for lighter alkali atoms and all the more so when using a quantum description since, in this case, they sit further away from the droplet surface on average due to their zero-point motion. Spectral modifications due to the presence of a vortex line are minor for np ← ns excitation and therefore insufficient for vortex detection. In the case of higher n'p ← ns or n's ← ns (n' > n) excitations, the shifts are larger as the excited state orbital is more extended and therefore more sensitive to changes in the surrounding helium density.