During the last decade, many exciting phenomena have been experimentally observed and theoretically predicted for ultracold atoms in optical lattices. This paper reviews these rapid developments concentrating mainly on the theory. Different types of
the bosonic systems in homogeneous lattices of different dimensions as well as in the presence of harmonic traps are considered. An overview of the theoretical methods used for these investigations as well as of the obtained results is given. Available experimental techniques are presented and discussed in connection with theoretical considerations. Eigenstates of the interacting bosons in homogeneous lattices and in the presence of harmonic confinement are analysed. Their knowledge is essential for understanding of quantum phase transitions at zero and finite temperature.
We study a quantum quench in the Bose-Hubbard model where the tunneling rate $J$ is suddenly switched from zero to a finite value in the Mott regime. In order to solve the many-body quantum dynamics far from equlibrium, we consider the reduced densit
y matrices for a finite number (one, two, three, etc.) of lattice sites and split them up into on-site density operators, i.e., the mean field, plus two-point and three-point correlations etc. Neglecting three-point and higher correlations, we are able to numerically simulate the time-evolution of the few-site density matrices and the two-point quantum correlations (e.g., their effective light-cone structure) for a comparably large number ${cal O}(10^3)$ of lattice sites.
We study the Bose and Fermi Hubbard model in the (formal) limit of large coordination numbers $Zgg1$. Via an expansion into powers of $1/Z$, we establish a hierarchy of correlations which facilitates an approximate analytical derivation of the time-e
volution of the reduced density matrices for one and two sites etc. With this method, we study the quantum dynamics (starting in the ground state) after a quantum quench, i.e., after suddenly switching the tunneling rate $J$ from zero to a finite value, which is still in the Mott regime. We find that the reduced density matrices approach a (quasi) equilibrium state after some time. For one lattice site, this state can be described by a thermal state (within the accuracy of our approximation). However, the (quasi) equilibrium state of the reduced density matrices for two sites including the correlations cannot be described by a thermal state. Thus, real thermalization (if it occurs) should take much longer time. This behavior has already been observed in other scenarios and is sometimes called ``pre-thermalization. Finally, we compare our results to numerical simulations for finite lattices in one and two dimensions and find qualitative agreement.
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.
