We give a variational Monte Carlo description at T = 0 K of 4He filling a porous glass within the shadow wave function technique. The confining media is modeled by a smooth cylindrical pore, and we have considered two different pore radii: R = 13 A and R = 22.5 A. The radial density profiles show a strong layering of the 4He atoms. In the R = 13 A pore, as the density is increased, solidification takes place layer by layer, starting from the pore wall, and the radius is too small to allow a bulk‐like solid to nucleate in the small liquid region at the center of the pore. For the R = 22.5 A pore we find an hcp lattice accommodated in the central region surrounded by four layers which turn out to be solid at the density of the computation. Computing the one‐body density matrix (in the R = 13 A case) we are able to estimate the Bose‐Einstein condensate fraction, which is still non‐zero even when the whole system is in the solid phase.
We have computed the one-body density matrix ${\ensuremath{\rho}}_{1}$ in solid
4He at
T=0K using the shadow wave-function (SWF) variational technique. The accuracy of the SWF has been tested with an exact projector method. We find that off-diagonal long-range order is present in ${\ensuremath{\rho}}_{1}$ for a perfect hcp and bcc solid
4He for a range of densities above the melting one, at least up to 54 bars. This is first microscopic indication that Bose Einstein condensation (BEC) is present in perfect solid
4He. At melting, the condensate fraction in the hcp solid is $5\ifmmode\times\else\texttimes\fi{}{10}^{\ensuremath{-}6}$ and it decreases by increasing the density. The key process giving rise to BEC is the formation of vacancy interstitial pairs. We also present values for Leggett's upper bound on the superfluid fraction deduced from the exact local density.
We present path integral Monte Carlo (PIMC) calculations of the superfluid fraction, ${\ensuremath{\rho}}_{S}/\ensuremath{\rho}$, and the one-body density matrix (OBDM) [Bose-Einstein condensation (BEC)] of liquid
4He confined in nanopores. Liquid
4He in nanopores represents a dense Bose liquid at reduced dimension and in disorder. The goal is to determine the effective dimensions of the liquid in the pores. It is to test whether observed properties, such as a very low onset temperature for superflow,
Tc, can be predicted by a standard, static PIMC ${\ensuremath{\rho}}_{S}/\ensuremath{\rho}$. We simulate a cylinder of liquid of diameter
dL surrounded by 5 \AA{} of inert solid
4He in a nanopore of diameter
d;
d=dL+10A˚. We find a PIMC ${\ensuremath{\rho}}_{S}(T)/\ensuremath{\rho}$ and OBDM that scales as a 1D fluid Luttinger liquid at extremely small liquid pore diameters only,
dL=6A˚. At this
dL, the liquid fills the pore in a 1D line at the center of the pore and there is no PIMC superflow. In the range $8\ensuremath{\le}{d}_{L}\ensuremath{\le}22\phantom{\rule{0.28em}{0ex}}\AA{}$ the PIMC ${\ensuremath{\rho}}_{S}(T)/\ensuremath{\rho}$ scales as a 2D liquid. In this
dL range the liquid fills the pores in 2D-like cylindrical layers. The crossover from no superflow at
d=16A˚ to superflow at $d\ensuremath{\ge}18\phantom{\rule{0.28em}{0ex}}\AA{}$ agrees with experiment. There is a crossover to 3D scaling at larger ${d}_{L}\ensuremath{\simeq}22\phantom{\rule{0.28em}{0ex}}\AA{}$. In the range $8\ensuremath{\le}{d}_{L}\ensuremath{\le}22\phantom{\rule{0.28em}{0ex}}\AA{}$, the
Tc predicted using the Kosterlitz-Thouless 2D scaling criterion of the OBDM agrees well with that obtained from ${\ensuremath{\rho}}_{S}(T)/\ensuremath{\rho}$. These results suggest that the superflow observed in small pore media is standard static superflow with the low
Tc arising from its 2D character. An operational onset temperature,
TBEC, for BEC can be defined as the temperature at which there is a crossover from exponential to algebraic decay in the OBDM. This definition leads to a ${T}_{\mathrm{BEC}}\ensuremath{\ge}{T}_{c}$ as observed in larger pore media.
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.
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.
There is a growing interest in the relation between Bose-Einstein condensation (BEC) and the superfluidity. A Bose system confined in random media such as porous glass is suitable for studying this relation because BEC and superfluidity can be suppressed and controlled in such a disordered environment. However, it is not clear how this relation is affected by disorder and there are few theoretical studies that can be quantitatively tested by experiment. In this work we develop a dilute Bose gas model with a random potential that takes into account the pore-size dependence of porous glass. Then we compare our model with the measured low-temperature specific heat, condensate density, and superfluid density of
4He in Vycor glass. This comparison uses no free parameters. We predict phenomena at low temperatures: First, the random potential causes a T-linear specific heat instead of the
T3 dependence that is usually caused by phonons. Second, the BEC can remain even when the superfluidity disappears at low densities. Third, the system makes a reentrant transition at low densities; that is, the superfluid phase changes to the normal phase again as the temperature is reduced. This reentrant transition is more likely to be observed when the strength of the random potential is increased.