Microsolvation of an Ionic Dopant in Small4He Clusters: OH+(3Σ)(4He)N via Genetic Algorithm Optimizations

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.

Quantum effects on the stability of the He5 I2 van der Waals conformers

We present 15-dimensional quantum multiconfiguration time-dependent Hartree calculations of the vibrational levels of the He₅ I₂ van der Waals (vdW) complex employing an ab initio-based potential energy surface (PES). The energies and spatial features of such bound structures are analyzed, providing predictions on the structures and relative stabilities of its three lowest isomers. We found that the most stable isomer corresponds to all five He atoms encircling the I₂ molecule, indicating that in this case the anharmonic quantum effects do not stabilize the isomers involving a He atom in a linear configuration as reported previously for the smaller He_N I₂ systems. Such finding provides information on the overall structuring of the finite-size-solvent systems, highlighting the intriguing interplay between weak intermolecular interactions and quantum effects. © 2019 Wiley Periodicals, Inc.

Microsolvation of CH+ in helium: An ab initio study

In the present study, microsolvation and interaction of the CH[Formula: see text] cation with He[Formula: see text] clusters are investigated by means of ab initio calculations at MP2/aug-cc-pVTZ and QCISD/aug-cc-pVTZ levels. Stabilization energies of the studied complexes including basis set superposition error (BSSE) and zero point energy (ZPE) corrections are in the range of [Formula: see text][Formula: see text]kJ/mol and [Formula: see text][Formula: see text]kJ/mol. A good linear correlation is found between the stabilization energy and stretching frequency shift ([Formula: see text]) in the studied complexes. According to energy decomposition analysis, it is found that polarization effects are the major source of the attraction in these complexes.

IR spectroscopy of protonated leu-enkephalin and its 18-crown-6 complex embedded in helium droplets

Ultracold IR spectra of protonated leu-enkephalin and its 18-crown-6 complex embedded in superfluid helium droplets have been recorded using a free-electron laser as radiation source.

Energetics and solvation structure of a dihalogen dopant (I2) in (4)He clusters

The energetics and structure of small HeNI2 clusters are analyzed as the size of the system changes, with N up to 38. The full interaction between the I2 molecule and the He atoms is based on analytical ab initio He-I2 potentials plus the He-He interaction, obtained from first-principle calculations. The most stable structures, as a function of the number of solvent He atoms, are obtained by employing an evolutionary algorithm and compared with CCSD(T) and MP2 ab initio computations. Further, the classical description is completed by explicitly including thermal corrections and quantum features, such as zero-point-energy values and spatial delocalization. From quantum PIMC calculations, the binding energies and radial/angular probability density distributions of the thermal equilibrium state for selected-size clusters are computed at a low temperature. The sequential formation of regular shell structures is analyzed and discussed for both classical and quantum treatments.

Weakly bound finite systems: (⁴He)N-Rb₂(³Σu), clustering structures from a quantum Monte Carlo approach

We report here ((4)He)(N)-Rb(2)((3)Σ(u)) complexes, 2 ≤ N ≤ 20, analysed through a quantum diffusion Monte Carlo stochastic approach. The calculations show that the spin stretched dimer molecule is bound outside the pure He sub-complex, due to the stronger He-He potential as compared with the He-Rb(2) interaction, while the rare gas atom moiety presents, in turn, a shell-like structure with ten He adatoms completing the first shell. Our results agree with previous findings on this and similarly weakly interacting systems.

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.

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.

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.

Energetics and structures of charged helium clusters: comparing stabilities of dimer and trimer cationic cores

We present accurate ab initio calculations of the most stable structures of He(n)(+) clusters in order to determine the more likely ionic core arrangements existing after reaching structural equilibrium of the clusters. Two potential energy surfaces are presented: one for the He(2)(+) and the other with the He(3)(+) linear ion, both interacting with one He atom. The two computed potentials are in turn employed within a classical structure optimization where the overall interaction forces are obtained within the sum-of-potentials approximation described in the main text. Because of the presence of many-body effects within the ionic core, we find that the arrangements with He(3)(+) as a core turn out to be energetically preferred, leading to the formation of He(3)(+)(He)(n-3) stable aggregates. Nanoscopic considerations about the relative stability of clusters with the two different cores are shown to give us new information on the dynamical processes observed in the impact ionization experiments of pure helium clusters and the importance of pre-equilibrium evaporation of the ionic dimers in the ionized clusters.

Microsolvation of Cationic Dimers in 4He Droplets: Geometries of (He)N (A = Li, Na, K) from Optimized Energies

Ab initio computed interaction forces are employed to describe the microsolvation of the A+2(2Sigma) (A=Li, Na, K) molecular ion in 4He clusters of small variable size. The minimum energy structures are obtained by performing energy minimization based on a genetic algorithm approach. The symmetry features of the collocation of solvent adatoms around the dimeric cation are analyzed in detail, showing that the selective growth of small clusters around the two sides of the ion during the solvation process is a feature common to all three dopants.