Superfluid and insulating phases in an interacting-boson model: mean-field theory and the RPA

The bosonic Hubbard model is studied via a simple mean-field theory. At zero temperature, in addition to yielding a phase diagram that is qualitatively correct, namely a superfluid phase for non-integer fillings and a Mott transition from a superfluid to an insulating phase for integer fillings, this theory gives results that are in good agreement with Monte Carlo simulations. In particular, the superfluid fraction obtained as a function of the interaction strength U for both integer and non-integer fillings is close to the simulation results. In all phases the excitation spectra are obtained by using the random phase approximation (RPA): the spectrum has a gap in the insulating phase and is gapless (and linear at small wave vectors) in the superfluid phase. Analytic results are presented in the limits of large U and small superfluid density. Finite-temperature phase diagrams and the Mott-insulator-normal-phase crossover are also described