Quantum Monte Carlo simulations of confined bosonic atoms in optical lattices

We study properties of ultracold bosonic atoms in one-, two-, and three-dimensional optical lattices by large scale quantum Monte Carlo simulations of the Bose-Hubbard model in parabolic confining potentials. Our results indicate that local properties of the atoms can be accessed by probing the system's response to local potential perturbations. Furthermore, we show how the formation of Mott insulating regions is reflected in the momentum distribution of the atoms, amenable to experimental detection. We disprove previous claims concerning the relevance of fine structure in the momentum distribution function. Furthermore, we discuss limitations of local density approximations for confined systems, and demonstrate the absence of quantum criticality due to the inhomogenous potential. Instead, we show that quantum critical behavior can be observed in flat confining potentials. Our results indicate that the experimental detection of the Mott transition in moderately sized optical lattices would be significantly eased in flat confinement potentials.