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)).
A simple and commonly employed approximate technique with which one can examine spatially disordered systems when strong electronic correlations are present is based on the use of real-space unrestricted self-consistent Hartree-Fock wave functions. I
n such an approach the disorder is treated exactly while the correlations are treated approximately. In this report we critique the success of this approximation by making comparisons between such solutions and the exact wave functions for the Anderson-Hubbard model. Due to the sizes of the complete Hilbert spaces for these problems, the comparisons are restricted to small one-dimensional chains, up to ten sites, and a 4x4 two-dimensional cluster, and at 1/2 filling these Hilbert spaces contain about 63,500 and 166 million states, respectively. We have completed these calculations both at and away from 1/2 filling. This approximation is based on a variational approach which minimizes the Hartree-Fock energy, and we have completed comparisons of the exact and Hartree-Fock energies. However, in order to assess the success of this approximation in reproducing ground-state correlations we have completed comparisons of the local charge and spin correlations, including the calculation of the overlap of the Hartree-Fock wave functions with those of the exact solutions. We find that this approximation reproduces the local charge densities to quite a high accuracy, but that the local spin correlations, as represented by < S_i . S_j >, are not as well represented. In addition to these comparisons, we discuss the properties of the spin degrees of freedom in the HF approximation, and where in the disorder-interaction phase diagram such physics may be important.
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 compare
d 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.
