We present a multi-site formulation of mean-field theory applied to the disordered Bose-Hubbard model. In this approach the lattice is partitioned into clusters, each isolated cluster being treated exactly, with inter-cluster hopping being treated ap
proximately. The theory allows for the possibility of a different superfluid order parameter at every site in the lattice, such as what has been used in previously published site-decoupled mean-field theories, but a multi-site formulation also allows for the inclusion of spatial correlations allowing us, e.g., to calculate the correlation length (over the length scale of each cluster). We present our numerical results for a two-dimensional system. This theory is shown to produce a phase diagram in which the stability of the Mott insulator phase is larger than that predicted by site-decoupled single-site mean-field theory. Two different methods are given for the identification of the Bose glass-to-superfluid transition, one an approximation based on the behaviour of the condensate fraction, and one of which relies on obtaining the spatial variation of the order parameter correlation. The relation of our results to a recent proposal that both transitions are non self-averaging is discussed.
We present a systematic study of various forms of renormalization that can be applied in the calculation of the self-energy of the Hubbard model within the T-matrix approximation. We compare the exact solutions of the attractive and repulsive Hubbard
models, for linear chains of lengths up to eight sites, with all possible taxonomies of the T-matrix approximation. For the attractive Hubbard model, the success of a minimally self-consistent theory found earlier in the atomic limit (Phys. Rev. B 71, 155111 (2005)) is not maintained for finite clusters unless one is in the very strong correlation limit. For the repulsive model, in the weak correlation limit at low electronic densities -- that is, where one would expect a self-consistent T-matrix theory to be adequate -- we find the fully renormalized theory to be most successful. In our studies we employ a modified Hubbard interaction that eliminates all Hartree diagrams, an idea which was proposed earlier (Phys. Rev. B 63, 035104 (2000)).
Partially-projected Gutzwiller variational wavefunctions are used to describe the ground state of disordered interacting systems of fermions. We compare several different variational ground states with the exact ground state for disordered one-dimens
ional chains, with the goal of determining a minimal set of variational parameters required to accurately describe the spatially-inhomogeneous charge densities and spin correlations. We find that, for weak and intermediate disorder, it is sufficient to include spatial variations of the charge densities in the product state alone, provided that screening of the disorder potential is accounted for. For strong disorder, this prescription is insufficient and it is necessary to include spatially inhomogeneous variational parameters as well.
