We investigate the thermodynamic properties of a dilute Bose gas in a correlated random potential using exact path integral Monte Carlo methods. The study is carried out in continuous space and disorder is produced in the simulations by a three-dimensional (3D) speckle pattern with tunable intensity and correlation length. We calculate the shift of the superfluid transition temperature due to disorder and we highlight the role of quantum localization by comparing the critical chemical potential with the classical percolation threshold. The equation of state of the gas is determined in the regime of strong disorder, where superfluidity is suppressed and the normal phase exists down to very low temperatures. We find a T 2 dependence of the energy in agreement with the expected behavior in the Bose glass phase. We also discuss the major role played by the disorder correlation length and we make contact with a Hartree–Fock mean-field approach that holds if the correlation length is very large. The density profiles are analyzed as a function of temperature and interaction strength. Effects of localization and the depletion of the order parameter are emphasized in the comparison between local condensate and total density. At very low temperature, we find that the energy and the particle distribution of the gas are very well described by the T =0 Gross–Pitaevskii theory, even in the regime of very strong disorder.
We investigate the thermodynamic properties of a dilute Bose gas in a correlated random potential using exact path integral Monte Carlo methods. The study is carried out in continuous space and disorder is produced in the simulations by a three-dimensional (3D) speckle pattern with tunable intensity and correlation length. We calculate the shift of the superfluid transition temperature due to disorder and we highlight the role of quantum localization by comparing the critical chemical potential with the classical percolation threshold. The equation of state of the gas is determined in the regime of strong disorder, where superfluidity is suppressed and the normal phase exists down to very low temperatures. We find a T 2 dependence of the energy in agreement with the expected behavior in the Bose glass phase. We also discuss the major role played by the disorder correlation length and we make contact with a Hartree–Fock mean-field approach that holds if the correlation length is very large. The density profiles are analyzed as a function of temperature and interaction strength. Effects of localization and the depletion of the order parameter are emphasized in the comparison between local condensate and total density. At very low temperature, we find that the energy and the particle distribution of the gas are very well described by the T =0 Gross–Pitaevskii theory, even in the regime of very strong disorder.
The superfluid transition of a three-dimensional gas of hard-sphere bosons in a disordered medium is studied using quantum Monte Carlo methods. Simulations are performed in continuous space both in the canonical and in the grand-canonical ensemble. At fixed density we calculate the shift of the transition temperature as a function of the disorder strength, while at fixed temperature we determine both the critical chemical potential and the critical density separating normal and superfluid phases. In the regime of strong disorder the normal phase extends up to large values of the degeneracy parameter, and the critical chemical potential exhibits a linear dependence in the intensity of the random potential. The role of interactions and disorder correlations is also discussed.
We present precise path-integral Monte Carlo results for the thermodynamics of a homogeneous dilute Bose gas. Pressure and energy are calculated as a function of temperature both below and above the Bose-Einstein transition. Specifically, we address interaction effects, focusing on deviations from the ideal-gas law in the thermodynamic limit. We also calculate the isothermal compressibility and the contact parameter, which provide a clear signature of the role played by interactions. In particular, we obtain indications of a discontinuity of the compressibility at the transition point. To gain physical insight, numerical results are systematically compared with the predictions of first-order Hartree-Fock and second-order Popov theories, both giving an approximate description of the gas thermodynamics. The comparison shows the extension of the critical region around the transition point, where the inaccuracies of the perturbative expansions are more pronounced.
By using exact path-integral Monte Carlo methods we calculate the equation of state of an interacting Bose gas as a function of temperature both below and above the superfluid transition. The universal character of the equation of state for dilute systems and low temperatures is investigated by modeling the interatomic interactions using different repulsive potentials corresponding to the same s-wave scattering length. The results obtained for the energy and the pressure are compared to the virial expansion for temperatures larger than the critical temperature. At very low temperatures we find agreement with the ground-state energy calculated using the diffusion Monte Carlo method.
Weakly interacting Bose gases in a disorder environment have long been a challenging topic in the field of solid-state physics due to the intriguing interplay between superfluidity and localization. In this thesis we investigate theoretically the equilibrium properties of one- and three-dimensional harmonically trapped dirty Bose-Einstein condensates (BECs). In particular we focus on the decomposition of the total particle density into three components, namely, the condensate density, the thermal density, and the density of fragmented local Bose-Einstein condensates within the respective minima of the random potential landscape. In order to determine the different density components we develop a Hartree-Fock mean-field theory for a weakly interacting trapped BEC in a random environment, then specialize it to contact interaction and delta-correlated disorder. This non-perturbative theory is worked out on the basis of the replica method, which represents a well- established tool to deal with disorder problems. With this we calculate the corresponding free energy at finite temperature and derive from it the underlying self-consistency equations for the respective density components. Then, as a first step, we investigate a quasi one-dimensional Bose-Einstein condensed gas in a harmonic trapping potential with an additional delta- correlated disorder potential at zero temperature. To this end we solve the self-consistency equations within the Thomas-Fermi approximation and find the emergence of a Bose-glass region, where the condensate vanishes. We corroborate this analysis by an elaborate numerical treatment, where the corresponding one-dimensional time-independent Gross-Pitaevskii equation is numerically solved and disorder-ensemble averages of the condensate wave function are performed. The variance of the condensate wave function quantifies the number of bosons which condense in the local minima of the random potential. For weak disorder these mini-condensates turn out to occur preferentially at the border of the condensate, while for intermediate disorder strength this happens in the trap center. Additionally we use a variational ansatz in order to describe analytically the numerically observed redistribution of the fragmented mini-condensates with increasing disorder strength. In close analogy to the one-dimensional case we then treat the three-dimensional dirty Bose gas in an isotropic harmonic trap and a delta- correlated disorder potential at zero temperature using the corresponding self-consistency equations obtained via the Hartree-Fock mean-field theory within the Thomas-Fermi approximation. Additionally we use a variational ansatz, whose results turn out to coincide qualitatively with those obtained from the Thomas-Fermi approximation. Inparticular, a first order quantum phase transition from the superfluid phase, where the condensate density contributes to the total density, to the Bose-glass phase, where all particles are in the mini-condensates, is detected at a critical disorder strength, which agrees with findings in the literature. Furthermore, in a generaltriaxial harmonic trap, we investigate the geometric effect of different trap aspect ratios on the respective properties of the dirty BEC system. Finally, we consider the three- dimensional dirty BEC at finite temperature. This allows us to study the impact of both temperature and disorder fluctuations on the respective components of the density as well as their Thomas-Fermi radii. In particular, we find that the superfluid region, the Bose-glass region, and the thermal region coexist. Furthermore, depending on the system parameters, three phase transitions are detected, namely, one from the superfluid to the Bose-glass phase, one from the Boseglass to the thermal phase, and finally one directly from the superfluid to the thermal phase. All these results could be particularly useful for a quantitative analysis of world-wide on-going experiments with dirty bosons in quasi one-dimensional and three-dimensional harmonic traps.
We revisit the problem of Anderson localization in a trapped Bose–Einstein condensate in 1D and 3D in a disordered potential, applying Quantum Monte Carlo technique because the disorder cannot be treated accurately in a perturbative way as even a small amount of disorder can produce dramatic changes in the physical properties of the system under investigation. Till date no unambiguous evidence of localization has been observed for matter waves in 3D. Matter waves made up of cold atoms are good candidates for such investigations. Simulations are performed for Rb gas in continuous space using canonical ensemble in the case of random and quasi-periodic potentials. To realize random and quasiperiodic potentials numerically we use speckle and bichromatic potentials, respectively. Owing to the high degree of control over the system parameters we specifically study the interplay of disorder and interaction in the system. A dilute Bose gas placed in a random environment falls into a fragmented localized state and the ergodicity (the repetitiveness of the wave function) is lost. An arbitrary Interaction can slowly overcome the effect of disorder and restore the ergodicity again. We observe that as the interaction strength increases, the wave functions become more and more delocalized. Since vanishing of Lyapunov exponent is only a necessary but not a sufficient condition for delocalization for probing the localization we calculate the mean square displacements as an alternative measure of localization. The path integral Monte Carlo technique in this paper numerically establishes the existing predictions of the scaling theory so far and paves a clear path for the further investigation of scaling theory to calculate more complicated properties like ‘critical exponents’ etc. in disordered quantum gases.
We study the two-dimensional weakly repulsive Bose gas at zero temperature in the presence of correlated disorder. Using large-scale simulations, we show that the low-energy Bogoliubov cumulative density of states remains quadratic up to a critical disorder strength, beyond which a power law with disorder-dependent exponent $\ensuremath{\beta}<2$ sets in. We associate this threshold behavior with the transition from superfluid to Bose glass, and compare the resulting mean-field phase diagram with scaling laws and the Thomas-Fermi percolation threshold of the mean-field density profile.