We present in this work the study of small (4)He(N)-Cs(2)((3)Σ(u)) aggregates (2 ≤ N ≤ 30) through combined variational, diffusion Monte Carlo (DMC), and path integral Monte Carlo (PIMC) calculations. The full surface is modeled as an addition of He-Cs(2) interactions and He-He potentials. Given the negligible strength and large range of the He-Cs(2) interaction as compared with the one for He-He, a propensity of the helium atoms to pack themselves together, leaving outside the molecular dopant is to be expected. DMC calculations determine the onset of helium gathering at N = 3. To analyze energetic and structural properties as a function of N, PIMC calculations with no bosonic exchange, i.e., Boltzmann statistics, at low temperatures are carried out. At T = 0.1 K, although acceptable one-particle He-Cs(2) distributions are obtained, two-particle He-He distributions are not well described, indicating that the proper symmetry should be taken into account. PIMC distributions at T = 1 K already compare well with DMC ones and show minor exchange effects, although binding energies are still far from having converged in terms of the number of quantum beads. As N increases, the He-He PIMC pair correlation function shows a clear tendency to coincide with the experimental boson-liquid helium one at that temperature. It supports the picture of a helium droplet which carries the molecular impurity on its surface, as found earlier for other triplet dimers.
We present in this work the study of small (4)He(N)-Cs(2)((3)Σ(u)) aggregates (2 ≤ N ≤ 30) through combined variational, diffusion Monte Carlo (DMC), and path integral Monte Carlo (PIMC) calculations. The full surface is modeled as an addition of He-Cs(2) interactions and He-He potentials. Given the negligible strength and large range of the He-Cs(2) interaction as compared with the one for He-He, a propensity of the helium atoms to pack themselves together, leaving outside the molecular dopant is to be expected. DMC calculations determine the onset of helium gathering at N = 3. To analyze energetic and structural properties as a function of N, PIMC calculations with no bosonic exchange, i.e., Boltzmann statistics, at low temperatures are carried out. At T = 0.1 K, although acceptable one-particle He-Cs(2) distributions are obtained, two-particle He-He distributions are not well described, indicating that the proper symmetry should be taken into account. PIMC distributions at T = 1 K already compare well with DMC ones and show minor exchange effects, although binding energies are still far from having converged in terms of the number of quantum beads. As N increases, the He-He PIMC pair correlation function shows a clear tendency to coincide with the experimental boson-liquid helium one at that temperature. It supports the picture of a helium droplet which carries the molecular impurity on its surface, as found earlier for other triplet dimers.
In this paper, we study the energy and geometric properties of Cs 2 ( 3 Σ u )–( 4 He) N clusters, 2⩽ N ⩽20, N even, through a diffusion Monte Carlo methodology. Considering the results for doped clusters in which the He–impurity interaction dominates over the He–He one, our aim is to investigate the case when this assumption is not fulfilled anymore and the helium–helium potential becomes dominant. We find a scenario where a pure helium subcomplex is formed, leaving out the alkaline dimer, with the largest species gathering the helium atoms in a two-shell-like structure, the first of which is filled with ten particles. Our results are in agreement with previous theoretical and experimental findings.
Helium systems at zero temperature are strongly correlated neutral systems which can only be described by means of quantum mechanics. The existence of the two isotopes 4He and 3He allows the comparison between the quantum effects due to different statistics, bosonic and fermionic. In this Thesis, an analysis of Helium droplets has been done. Throughout this work, both energetic and structural properties of systems with different number of particles N are obtained by means of the Quantum Monte Carlo techniques (Variational and Diffusion Monte Carlo - VMC and DMC-).
Several studies on Helium droplets have been published in the last years, most of them using the mean field theory (as the Density Functional) and variational approaches (as VMC). In the case of 4He, also exact DMC results have already been reported, actually, a great number of works have described with high accuracy the ground state of this bosonic system. The appearance of magic numbers from experimental point of view has motivated the study of 4He droplets, and one chapter of this Thesis is focused on this bosonic system. We have carried out a one-by-one energetic study in the range from 20 to 30 particles using the DMC technique to discuss the existence of magic numbers. On the other hand, a brief study about vortices is also done; using out code, one can calculate the energy of a droplet with one vortex inside (EV) and, therefore, the energy of a vortex (EV -E0).
About 3He droplets, the sign problem makes their study quite more difficult, for this reason, the number of works is considerably smaller than in the bosonic case. One aim of this thesis was to write a DMC code using both Fixed and Released Node approximations (FN and RN) able to perform high accurate calculations of energetic and structural properties of these finite fermionic systems.
A new full-dimension potential energy surface of the three-body He-Rb₂(³Σ(u)(+)) complex and a quantum study of small (⁴He)(N)-Rb₂(³Σ(u)(+)) clusters, 1 ≤ N ≤ 4, are presented. We have accurately fitted the ab initio points of the interaction to an analytical form and addressed the dopant's vibration, which is found to be negligible. A Variational approach and a Diffusion Monte Carlo technique have been applied to yield energy and geometric properties of the selected species. Our quantum structure calculations show a transition in the arrangements of the helium atoms from N = 2, where they tend to be separated across the diatomic bond, to N = 4, in which a closer packing of the rare gas particles is reached, guided by the dominance of the He-He potential over the weaker interaction of the latter adatoms with the doping dimer. The deepest well of the He-Rb₂ interaction is placed at the T-shape configuration, a feature which causes the dopant to be located as parallel to the helium "minidroplet". Our results are shown to agree with previous findings on this and on similar systems.
Spectroscopy of doped helium nanodroplets has become, in recent years, a popular tool for the study of bothfinite size samples of the helium quantum fluid and molecules isolated in this sub Kelvin environment. Whilemany dopants migrate to the center of the droplet, alkali (Ak) atoms in their ground electronic state, becauseof their so weak interaction with He, remain at the surface, the latter being in this way slightly distortedexhibiting a dimple. However, when in their first excited states, Ak atoms are forming exciplexes. These areanisotropic with the central atom surrounded by a ring of He atoms located in the nodal plane of the valenceelectron pz orbital. Structural properties, obtained by DMC and PIMC, are presented here for some Ak -Hen exciplexes.Ak2 dimers in helium environment have also gained interest in this last decade, in connexion with the generationof ultracold molecules. When attached to helium droplets, Ak atoms can float on the surface andform molecules in cold collisions. DMC and PIMC methods are used to compute energetics and structuralinformations for Rb2 in helium. In particular, the orientation of the dimer relative to the droplet surface isexamined.
Below Tλ = 2.17 K, bulk 4He is a superfluid and has a non-zero Bose-Einstein condensate fraction. This work will focus primarily on how phenomena such as superfluidity, Bose condensation and superfluid vortices are manifested in microscopic, inhomogeneous helium systems. Path Integral Monte Carlo is a powerful method for calculating the equilibrium properties of quantum systems at finite temperature. We have achieved linear scaling of computer time with number of particles through the use of neighbor lists, allowing us to simulate systems of several thousand atoms. We have derived a local superfluid estimator and used it to examine the microscopic superfluid response around a molecule rotating in a helium nanodroplet. We found that the first solvation layer is well-described by a two dimensional superfluid, with thermal excitations occuring at a lower temperature than in bulk helium. The effective moment of inertia of a linear impurity in a helium droplet is calculated, and compared with experimental scattering results. In addition, we calculated the vortex formation energy for both pure droplets and droplets doped with linear impurities, and found that the linear impurities had a negligible impact on the formation energy. A possible spectroscopic signature of vortices in doped helium droplets was suggested. After deriving a local estimator, we calculated the condensate fraction throughout the free helium surface of a semi-infinite slab. These results, along with densitydensity correlation functions, were used to characterize the surface excitations and calculate the extent to which ripplons are present. In addition, the ripplon dispersion relation was calculated using imaginary-time correlation functions, and found to be lll in good agreement with experimental results. Finally, we have calculated the Dey be-Waller factor in solid helium for a range of temperatures and densities, and compared the scaling behavior with the predictions of harmonic theory. The first non-Gaussian contribution to the density distribution was calculated.
Small helium (4He) clusters containing the lighter isotope He3 are studied by means of quantum Monte Carlo methods. Accurate ground state energies and structural properties are obtained using accurate trial wave functions and the Tang–Tonnies–Yiu (TTY) helium–helium pair potential. The dimer He4–3He is not bound; as well as the trimer He4He23. The smallest cluster containing He3 is He24He3 with a nonrigid structure having a marked linear contribution. Interestingly, this weakly bound system, with an energy one order of magnitude less than the He34 trimer, is able to bind another He3 atom, forming the tetramer He24He23, which shows the odd feature of having five out of six unbound pairs. In general, the substitution of a single He4 atom in a pure cluster with a He3 atom leads to an energetic destabilization, as the pair He4–3He is not bound. The isotopic impurity is found to perturb only weakly the distributions of the remaining He4 atoms, which retain the high floppiness already found in the pure clusters. As the number of atoms increases the isotopic impurity has the marked tendency to stay on the surface of the cluster. This behavior is consistent with the formation of the so-called “Andreev states” of a single He3 in liquid He4 helium and droplets, where the impurity tends to form single-particle states on the surface of the pure He4.