Exact calculation of the condensate fraction in multi-dimensional inhomogeneous interacting Bose systems which do not possess continuous symmetries is a difficult computational problem. We have developed an iterative procedure which allows to calcula
te the condensate fraction as well as the corresponding eigenfunction of the one-body density matrix. We successfully validate this procedure in diffusion Monte Carlo simulations of a Bose gas in an optical lattice at zero temperature. We also discuss relation between different criteria used for testing coherence in cold Bose systems, such as fraction of particles that are superfluid, condensed or are in the zero-momentum state.
Quantum phases of ultracold bosons with repulsive interactions in lattices in the presence of quenched disorder are investigated. The disorder is assumed to be caused by the interaction of the bosons with impurity atoms having a large effective mass.
The system is described by the Bose-Hubbard Hamiltonian with random on-site energies which have a discrete binary probability distribution. The phase diagram at zero temperature is calculated using several methods like a strong-coupling expansion, an exact numerical diagonalization, and a Bose-Fermi mapping valid in the hard-core limit. It is shown that the Mott-insulator phase exists for any strength of disorder in contrast to the case of continuous probability distribution. We find that the compressibility of the Bose glass phase varies in a wide range and can be extremely low. Furthermore, we evaluate experimentally accessible quantities like the momentum distribution, the static and dynamic structure factors, and the density of excited states. The influence of finite temperature is discussed as well.
Bose-Einstein-condensed gases in external spatially random potentials are considered in the frame of a stochastic self-consistent mean-field approach. This method permits the treatment of the system properties for the whole range of the interaction s
trength, from zero to infinity, as well as for arbitrarily strong disorder. Besides a condensate and superfluid density, a glassy number density due to a spatially inhomogeneous component of the condensate occurs. For very weak interactions and sufficiently strong disorder, the superfluid fraction can become smaller than the condensate fraction, while at relatively strong interactions, the superfluid fraction is larger than the condensate fraction for any strength of disorder. The condensate and superfluid fractions, and the glassy fraction always coexist, being together either nonzero or zero. In the presence of disorder, the condensate fraction becomes a nonmonotonic function of the interaction strength, displaying an antidepletion effect caused by the competition between the stabilizing role of the atomic interaction and the destabilizing role of the disorder. With increasing disorder, the condensate and superfluid fractions jump to zero at a critical value of the disorder parameter by a first-order phase transition.
