Path Integral Monte Carlo was used to calculate the Bose-Einstein condensate fraction at the surface of a helium film at T = 0:77 K, as a function of density. Moving from the center of the slab to the surface, the condensate fraction was found to initially increase with decreasing density to a maximum value of 0.9, before decreasing. Long wavelength density correlations were observed in the static structure factor at the surface of the slab. A surface dispersion relation was calculated from imaginary-time density-density correlations. Similar calculations of the superfluid density throughout He-4 droplets doped with linear impurities (HCN)(sub n) are presented. After deriving a local estimator for the superfluid density distribution, we find a decreased superfluid response in the first solvation layer. This effective normal fluid exhibits temperature dependence similar to that of a two-dimensional helium system.
Path Integral Monte Carlo was used to calculate the Bose-Einstein condensate fraction at the surface of a helium film at T = 0:77 K, as a function of density. Moving from the center of the slab to the surface, the condensate fraction was found to initially increase with decreasing density to a maximum value of 0.9, before decreasing. Long wavelength density correlations were observed in the static structure factor at the surface of the slab. A surface dispersion relation was calculated from imaginary-time density-density correlations. Similar calculations of the superfluid density throughout He-4 droplets doped with linear impurities (HCN)(sub n) are presented. After deriving a local estimator for the superfluid density distribution, we find a decreased superfluid response in the first solvation layer. This effective normal fluid exhibits temperature dependence similar to that of a two-dimensional helium system.
A path integral Monte Carlo method was used to calculate the Bose-Einstein condensate fraction at the surface of a helium film at T = 0.77 K, as a function of density. Moving from the center of the slab to the surface, the condensate fraction was found to initially increase with decreasing density to a maximum value of 0.9 before decreasing. Long wavelength density correlations were observed in the static structure factor at the surface of the slab. Finally, a surface dispersion relation was calculated from imaginary-time density-density correlations.
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.
We report our results about the fluctuation effects on the surface static structure factor and on the local Bose-Einstein condensate (BEC) in the low density surface region of 4 He. This information has been obtained by a variational Monte Carlo calculation for a slab of 4 He modelled by a novel shadow wave function with a glue term (glue-SWF). This glue-SWF describes the self-binding of 4 He only via local interparticle correlations. BEC increases from the bulk-like value well inside the slab to a much larger value in the surface region (up to about 50%). Further out of the surface the condensate decreases and correspondingly there is an enhanced population of small momentum states. This is very different from the full BEC that was predicted using a simplified variational treatment of these 4 He systems. This different behaviour is correlated with the presence in the glue-SWF of enhanced density fluctuations in the surface region as shown by the behaviour of the static structure factor at small wave vector. This is interpreted as an effect due to the zero point motion of ripplons.
We evaluate the Bose-Einstein condensate density and the superfluid fraction of bosons in a periodic external potential using path-integral Monte Carlo methods. The periodic lattice consists of a cubic cell containing a potential well that is replicated along one dimension using periodic boundary conditions. The aim is to describe bosons in a one-dimensional optical lattice or helium confined in a periodic porous medium. The one-body density matrix is evaluated and diagonalized numerically to obtain the single boson natural orbitals and the occupation of these orbitals. The condensate fraction is obtained as the fraction of bosons in the orbital that has the highest occupation. The superfluid density is obtained from the winding number. From the condensate orbital and superfluid fraction, we investigate (1) the impact of the periodic external potential on the spatial distribution of the condensate, and (2) the correlation of localizing the condensate into separated parts and the loss of superflow along the lattice. For high-density systems, as the well depth increases, the condensate becomes depleted in the wells and confined to the plateaus between successive wells, as in pores between necks in a porous medium. For low-density systems, as the well depth increases the BEC is localized at the center of the wells (tight binding) and depleted between the wells. In both cases, the localization of the condensate suppresses superflow leading to a superfluid-insulator crossover. The impact of the external potential on the temperature dependence of the superfluidity is also investigated. The external potential suppresses the superfluid fraction at all temperatures, with a superfluid fraction significantly less than one at low temperature. The addition of an external potential does not, however, significantly reduce the transition temperature.
The experimental realization of a thin layer of spin-polarized hydrogen H↓ adsorbed on top of the surface of superfluid (4)He provides one of the best examples of a stable, nearly two-dimensional (2D) quantum Bose gas. We report a theoretical study of this system using quantum Monte Carlo methods in the limit of zero temperature. Using the full Hamiltonian of the system, composed of a superfluid (4)He slab and the adsorbed H↓ layer, we calculate the main properties of its ground state using accurate models for the pair interatomic potentials. Comparing the results for the layer with the ones obtained for a strictly 2D setup, we analyze the departure from the 2D character when the density increases. Only when the coverage is rather small the use of a purely 2D model is justified. The condensate fraction of the layer is significantly larger than in 2D at the same surface density, being as large as 60% at the largest coverage studied.
Several microscopic theories point out that Bose-Einstein condensation (BEC), i.e., a macroscopic occupation of the lowest energy single particle state in many-boson systems, may appear also in quantum fluids and solids and that it is at the origin of the phenomenon of superfluidity. Nevertheless, the connection between BEC and superfluidity is still matter of debate, since the experimental evidences indicating a non zero condensate fraction in superfluid helium are only indirect. In the theoretical study of BEC in quantum fluids and solids, perturbative approaches are useless because of the strong correlations between the atoms, arising both from the interatomic potential and from the quantum nature of the system. Microscopic Quantum Monte Carlo simulations provide a reliable description of these systems. In particular, the Path Integral Monte Carlo (PIMC) method is very suitable for this purpose. This method is able to provide exact results for the properties of the quantum system, both at zero and finite temperature, only with the definition of the Hamiltonian and of the symmetry properties of the system, giving an easy picture for superfluidity and BEC in many-boson systems. In this thesis, we apply PIMC methods to the study of several quantum fluids and solids. We describe in detail all the features of PIMC, from the sampling methods to the estimators of the physical properties. We present also the most recent techniques, such as the high-order approximations for the thermal density matrix and the worm algorithm, used in PIMC to provide reliable simulations. We study the liquid phase of condensed 4He, providing unbiased estimations of the one-body density matrix g1(r). We analyze the model for g1(r) used to fit the experimental data, highlighting its merits and its faults. In particular we see that, even if it presents some difficulties in the description of the overall behavior of g1(r), it can provide an accurate estimation of the kinetic energy K and of the condensate fraction n0 of the system. Furthermore, we show that our results for n0 as a function of the pressure are in a good agreement with the most recent experimental results. The study of the solid phase of 4He is the most significant part of this thesis. The recent observation of non classical rotational inertia (NCRI) effects in solid helium has generated big interest in the study of an eventual supersolid phase, characterized at the same time by crystalline order and superfluidity. Nevertheless, until now it has been impossible to give a theoretical model able to describe all the experimental evidences. In this work, we perform PIMC simulations of 4He at high densities, according to different microscopic configurations of the atoms. In commensurate crystals we see that BEC does not appear, our model being able to reproduce the momentum distribution obtained form neutron scattering experiments. In a crystal with vacancies, we have been able to see a transition to a superfluid phase at temperatures in agreement with experimental results if the vacancy concentration is low enough. In amorphous solids, superfluid effects are enhanced but appear at temperatures higher than the experimental estimation for the transition temperature. Finally, we study also metastable disordered configurations in molecular para-hydrogen at low temperature. The aim of this study is to investigate if a Bose liquid other than helium can display superfluidity. Choosing accurately a ¿quantum liquid¿ initial configuration and the dimensions of the simulation box, we have been able to frustrate the formation of the crystal and to calculate the temperature dependence of the superfluid density, showing a transition to a superfluid phase at temperatures close to 1 K.