Multi-site mean-field theory for cold bosonic atoms in optical lattices

We present a detailed derivation of a multi-site mean-field theory (MSMFT) used to describe the Mott-insulator to superfluid transition of bosonic atoms in optical lattices. The approach is based on partitioning the lattice into small clusters which are decoupled by means of a mean field approximation. This approximation invokes local superfluid order parameters defined for each of the boundary sites of the cluster. The resulting MSMFT grand potential has a non-trivial topology as a function of the various order parameters. An understanding of this topology provides two different criteria for the determination of the Mott insulator superfluid phase boundaries. We apply this formalism to $d$-dimensional hypercubic lattices in one, two and three dimensions, and demonstrate the improvement in the estimation of the phase boundaries when MSMFT is utilized for increasingly larger clusters, with the best quantitative agreement found for $d=3$. The MSMFT is then used to examine a linear dimer chain in which the on-site energies within the dimer have an energy separation of $Delta$. This system has a complicated phase diagram within the parameter space of the model, with many distinct Mott phases separated by superfluid regions.

Exact diagonalization analysis of the Anderson-Hubbard model and comparison to real-space self-consistent Hartree-Fock solutions

We have obtained the exact ground state wave functions of the Anderson-Hubbard model for different electron fillings on a 4x4 lattice with periodic boundary conditions - for 1/2 filling such ground states have roughly 166 million states. When compared to the uncorrelated ground states (Hubbard interaction set to zero) we have found strong evidence of the very effective screening of the charge homogeneities due to the Hubbard interaction. We have successfully modelled these local charge densities using a non-interacting model with a static screening of the impurity potentials. In addition, we have compared such wave functions to self-consistent real-space unrestricted Hartree-Fock solutions and have found that these approximate ground state wave functions are remarkably successful at reproducing the local charge densities, and may indicate the role of dipolar backflow in producing a novel metallic state in two dimensions.