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 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 effect of Anderson localization on a Bose-Einstein condesate in 3d in a disordered potential by Feynman-Kac path integral technique. Simulations are performed in continuous space using canonical ensemble. Owing to the high degree of control over the system parameters we also study the interplay of disorder and interaction in the system. We numerically compute the localization length, mobility edge and the density profile of the condensate. We observe that as the interaction strength increases, the wave functions become more and more delocalized.
Anderson localization of waves in disordered media was originally predicted[1] fifty years ago, in the context of transport of electrons in crystals[2]. The phenomenon is much more general[3] and has been observed in a variety of systems, including light waves[4, 5]. However, Anderson localization has not been observed directly for matter waves. Owing to the high degree of control over most of the system parameters (in particular the interaction strength), ultracold atoms offer opportunities for the study of disorder-induced localization[6]. Here we use a non-interacting Bose-Einstein condensate to study Anderson localization. The experiment is performed with a one-dimensional quasi-periodic lattice - a system that features a crossover between extended and exponentially localized states, as in the case of purely random disorder in higher dimensions. Localization is clearly demonstrated through investigations of the transport properties and spatial and momentum distributions. We characterize the crossover, finding that the critical disorder strength scales with the tunnelling energy of the atoms in the lattice. This controllable system may be used to investigate the interplay of disorder and interaction (ref. 7 and references therein), and to explore exotic quantum phases[8, 9].
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.
By numerical and variational solutions of the Gross-Pitaevskii equation, we studied the localization of a noninteracting and weakly interacting Bose-Einstein condensate (BEC) in a disordered cold atom lattice and a speckle potential. In the case of a single BEC fragment, the variational analysis produced good results. For a weakly disordered potential, the localized BEC's are found to have an exponential tail as in the weak Anderson localization. We also investigated the expansion of a noninteracting BEC in these potentials. We find that the BEC will be locked in an appropriate localized state after an initial expansion and will execute breathing oscillation around a mean shape when a BEC at equilibrium in a harmonic trap is suddenly released into a disorder potential.
We study the interplay of disorder localization and strong local interactions within the Anderson-Hubbard model. Taking into account local Mott-Hubbard physics and static screening of the disorder potential, the system is mapped onto an effective single-particle Anderson model, which is studied within the self-consistent theory of electron localization. For fermions, we find rich nonmonotonic behavior of the localization length ξ, particularly in two-dimensional systems, including an interaction-induced exponential enhancement of ξ for small and intermediate disorders and a strong reduction of ξ due to hopping suppression by strong interactions. In three dimensions, we identify for half filling a Mott-Hubbard-assisted Anderson localized phase existing between the metallic and the Mott-Hubbardgapped phases. For small U there is re-entrant behavior from the Anderson localized phase to the metallic phase. For bosons, the unrestricted particle occupation number per lattice site yields a monotonic enhancement of ξ as a function of decreasing interaction, which we assume to persist until the superfluid Bose- Einstein condensate phase is entered. Besides, we study cold atomic gases expanding, by a diffusion process, in a weak random potential. We show that the density-density correlation function of the expanding gas is strongly affected by disorder and we estimate the typical size of a speckle spot, i.e., a region of enhanced or depleted density. Both a Fermi gas and a Bose-Einstein condensate (in a mean-field approach) are considered.
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.