In the absence of a confining potential, the boson-Hubbard model exhibits a superfluid to Mott insulator quantum phase transition at commensurate fillings and strong coupling. We use quantum Monte Carlo simulations to study the ground state of the one-dimensional bosonic Hubbard model in a trap. Some, but not all, aspects of the Mott insulating phase persist. Mott behavior occurs for a continuous range of incommensurate fillings, very different from the unconfined case, and the establishment of the Mott phase does not proceed via a traditional quantum phase transition. These results have important implications for interpreting experiments on ultracold atoms on optical lattices.

We use quantum Monte Carlo simulations and exact diagonalization to explore the phase diagram of the Bose-Hubbard model with an additional superlattice potential. We first analyze the properties of superfluid and insulating phases present in the hard-core limit where an exact analytic treatment is possible via the Jordan-Wigner transformation. The extension to finite on-site interactions is achieved by means of quantum Monte Carlo simulations. We determine insulator and superfluid phase diagrams as functions of the on-site repulsive interaction, superlattice potential strength, and filling, finding that insulators with fractional occupation numbers, which are present in the hard-core case, extend deep into the soft-core region. Furthermore, at integer fillings, we find that the competition between the on-site repulsion and the superlattice potential can produce a phase transition between a Mott insulator and a charge-density-wave insulator, with an intermediate superfluid phase. Our results are relevant to the behavior of ultracold atoms in optical superlattices which are beginning to be studied experimentally.

Recent progress in experiments on trapped ultracold atoms has made it possible to study the interplay between magnetism and superfluid-insulator transitions in the boson Hubbard model. We report on quantum Monte Carlo simulations of the spin-1 boson Hubbard model in the ground state. For antiferromagnetic interactions favoring singlets, we present exact numerical evidence that the superfluid-insulator transition is first (second) order for even (odd) Mott lobes. Inside even lobes, we search for nematic-to-singlet first order transitions. In the ferromagnetic case where transitions are all continuous, we map the phase diagram and show the superfluid to be ferromagnetic. We compare the quantum Monte Carlo phase diagram with a third order perturbation calculation.

We present an exact quantum Monte Carlo study of the attractive one-dimensional Hubbard model with imbalanced fermion population. The pair-pair correlation function, which decays monotonically in the absence of polarization P, develops oscillations when P is nonzero, characteristic of Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) phase. The pair momentum distribution peaks at a momentum equal to the difference in the Fermi momenta. At strong coupling, the minority and majority momentum distributions are shown to be deformed, reflecting the presence of the other species and its Fermi surface. The FFLO oscillations survive the presence of a confining potential, and the local polarization at the trap center exhibits a marked dip, similar to that observed experimentally.

We use quantum Monte Carlo simulations to obtain zero-temperature state diagrams for strongly correlated lattice bosons in one and two dimensions under the influence of a harmonic confining potential. Since harmonic traps generate a coexistence of superfluid and Mott insulating domains, we use local quantities such as the quantum fluctuations of the density and a local compressibility to identify the phases present in the inhomogeneous density profiles. We emphasize the use of the "characteristic density" to produce a state diagram that is relevant to experimental optical lattice systems, regardless of the number of bosons or trap curvature and of the validity of the local-density approximation. We show that the critical value of U/t at which Mott insulating domains appear in the trap depends on the filling in the system, and it is in general greater than the value in the homogeneous system. Recent experimental results by Spielman et al. [Phys. Rev. Lett. 100, 120402 (2008)] are analyzed in the context of our two-dimensional state diagram, and shown to exhibit a value for the critical point in good agreement with simulations. We also study the effects of finite, but low (T