Bose-Einstein condensation has been experimentally found to take place in finite trapped systems when one of the confining frequencies is increased until the gas becomes effectively two-dimensional (2D). We confirm the plausibility of this result by performing path-integral Monte Carlo (PIMC) simulations of trapped Bose gases of increasing anisotropy and comparing them to the predictions of finite-temperature many-body theory. PIMC simulations provide an essentially exact description of these systems; they yield the density profile directly and provide two different estimates for the condensate fraction. For the ideal gas, we find that the PIMC column density of the squeezed gas corresponds quite accurately to that of the exact analytic solution and, moreover, is well mimicked by the density of a 2D gas at the same temperature; the two estimates for the condensate fraction bracket the exact result. For the interacting case, we find 2D Hartree-Fock solutions whose density profiles coincide quite well with the PIMC column densities and whose predictions for the condensate fraction are again bracketed by the PIMC estimates.
Bose-Einstein condensation has been experimentally found to take place in finite trapped systems when one of the confining frequencies is increased until the gas becomes effectively two-dimensional (2D). We confirm the plausibility of this result by performing path-integral Monte Carlo (PIMC) simulations of trapped Bose gases of increasing anisotropy and comparing them to the predictions of finite-temperature many-body theory. PIMC simulations provide an essentially exact description of these systems; they yield the density profile directly and provide two different estimates for the condensate fraction. For the ideal gas, we find that the PIMC column density of the squeezed gas corresponds quite accurately to that of the exact analytic solution and, moreover, is well mimicked by the density of a 2D gas at the same temperature; the two estimates for the condensate fraction bracket the exact result. For the interacting case, we find 2D Hartree-Fock solutions whose density profiles coincide quite well with the PIMC column densities and whose predictions for the condensate fraction are again bracketed by the PIMC estimates.
The fact that two-dimensional interacting trapped systems do not undergo Bose-Einstein Condensation (BEC) in the thermodynamic limit, though rigorously proved, is somewhat mysterious because all relevant limiting cases (zero temperature, small atom numbers, noninteracting particles) suggest otherwise. We study the possibility of condensation in finite trapped two-dimensional systems. We first consider the ideal gas, which incorporates the inhomogeneity and finite size of experimental systems and can be solved exactly. A semiclassical self-consistent approximation gives us a feel for the temperature scales; diagonalization of the one-body density matrix confirms that the condensation is into a single state. We squeeze a three-dimensional system and see how it crosses over into two dimensions. Mean-field theory, our main tool for the study of interacting systems, prescribes coupled equations for the condensate and the thermal cloud: the condensate receives a full quantum-mechanical treatment, while the noncondensate is described by different schemes of varying sophistication. We analyze the T = 0 case and its approach to the thermodynamic limit, finding a criterion for the dimensionality crossover and obtaining the coupling constant of the two-dimensional system that results from squeezing a three-dimensional trap. We next apply a semiclassical Hartree-Fock approximation to purely two-dimensional finite gases and find that they can be described either with or without a condensate; this continues to be true in the thermodynamic limit. The condensed solutions have a lower free energy at all temperatures but neglect the presence of phonons within the system and cease to exist when we allow for this possibility. The uncondensed solutions, in turn, are valid under a more rigorous scheme but have consistency problems of their own. Path-integral Monte Carlo simulations provide an essentially exact description of finite interacting gases and enable us to study highly anisotropic systems at finite temperature. We find that our two-dimensional Hartree-Fock solutions accurately mimic the surface density profile and predict the condensate fraction of these systems; the equivalent interaction parameter is smaller than that dictated by the T = 0 analysis. We conclude that, in two-dimensional isotropic finite trapped systems and in highly compressed three-dimensional gases, there is a phenomenon resembling a condensation into a single state.
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.
We compare the semiclassical Hartree-Fock approximation for a trapped Bose gas to a direct path-integral quantum Monte Carlo simulation. The chosen parameters correspond to current 87Rb experiments. We observe corrections to the mean-field density profile: the path-integral calculation reveals an increase of the number of condensed particles, which is of the same order as a previously computed result for a homogeneous system. We discuss the experimental observability of the effect and propose a method to analyze data of in situ experiments.
We analyze the results of a recent experiment with bosonic rubidium atoms harmonically confined in a quasi-two-dimensional (2D) geometry. In this experiment a well-defined critical point was identified, which separates the high-temperature normal state characterized by a single component density distribution, and the low-temperature state characterized by a bimodal density distribution and the emergence of high-contrast interference between independent 2D clouds. We first show that this transition cannot be explained in terms of conventional Bose–Einstein condensation of the trapped ideal Bose gas. Using the local density approximation (LDA), we then combine the mean-field (MF) Hartree–Fock theory with the prediction for the Berezinskii–Kosterlitz–Thouless (BKT) transition in an infinite uniform system. If the gas is treated as a strictly 2D system, the MF predictions for the spatial density profiles significantly deviate from those of a recent quantum Monte Carlo (QMC) analysis. However, when the residual thermal excitation of the strongly confined degree of freedom is taken into account, excellent agreement is reached between the MF and the QMC approaches. For the interaction strength corresponding to the experiment, we predict a strong correction to the critical atom number with respect to the ideal gas theory (factor ∼2). Quantitative agreement between theory and experiment is reached concerning the critical atom number if the predicted density profiles are used for temperature calibration.
Using a finite-temperature path integral Monte Carlo simulation (PIMC) method and finite-size scaling, we have investigated the interaction-induced shift of the phase-transition temperature for Bose-Einstein condensation of homogeneous weakly interacting Bose gases in three dimensions, which is given by a proposed analytical expression ${T}_{c}={T}_{c}^{0}{1+{c}_{1}a{n}^{1∕3}+[{c}_{2}^{\ensuremath{'}}\phantom{\rule{0.2em}{0ex}}\mathrm{ln}(a{n}^{1∕3})+{c}_{2}^{\ensuremath{''}}]{a}^{2}{n}^{2∕3}+\mathrm{O}({a}^{3}n)}$, where Tc0 is the critical temperature for an ideal gas, a is the s-wave scattering length, and n is the number density. We have used smaller number densities and more time slices than in the previous PIMC simulations [Gruter et al., Phys. Rev. Lett. 79, 3549 (1997)] in order to understand the difference in the value of the coefficient c1 between their results and the (apparently) other reliable results in the literature. Our results show that ${({T}_{c}\ensuremath{-}{T}_{c}^{0})∕{T}_{c}^{0}}∕(a{n}^{1∕3})$ depends strongly on the interaction strength an1∕3 while the previous PIMC results are considerably flatter and smaller than our results. We obtain ${c}_{1}=1.32\ifmmode\pm\else\textpm\fi{}0.14$, in agreement with results from recent Monte Carlo methods of three-dimensional O(2) scalar ${\ensuremath{\phi}}^{4}$ field theory and variational perturbation theory.
We study thermal properties of a trapped Bose-Bose mixture in a dilute regime using quantum Monte Carlo methods. Our main aim is to investigate the dependence of the superfluid density and the condensate fraction on temperature for the mixed and separated phases. To this end we use the diffusion Monte Carlo method in the zero-temperature limit and the path-integral Monte Carlo method for finite temperatures. The results obtained are compared with solutions of the coupled Gross-Pitaevskii equations for the mixture at zero temperature. We notice the existence of an anisotropic superfluid density in some phase-separated mixtures. Our results also show that the temperature evolution of the superfluid density and condensate fraction is slightly different, showing noteworthy situations where the superfluid fraction is smaller than the condensate fraction.