Two-site correlation in analysis of the Hubbard model

Abstract

Electronic states of a two-dimensional Hubbard model near half filling is investigated by use of the composite-operator method, which is improved by including two-site correlations in combined use with the resolvent method. The densities of states obtained for the upper and lower Hubbard bands have three-peek structures at half filling and give shapes above a certain temperature similar to the ones obtained in the numerical simulation by Dagotto, Ortolani, and Scalapino [Phys. Rev. B 46, 3183 (1992)], but show narrowing of the bands at low temperature due to the formation of spin-singlet states among the nearest-neighbor electrons. A doped case shows a coherent peak at the Fermi level and the feature of the density of states is very similar to the one from numerical simulation. The one-particle spectra in the (k,k) direction is also presented. The results show that, at half filling, the spectral weight is not evenly distributed between two bands; most of the weight of the lower Hubbard band appears between (0,0) and (\ensuremath{\pi}/2a,\ensuremath{\pi}/2a), while that of the upper Hubbard band lies in the region between (\ensuremath{\pi}/2a,\ensuremath{\pi}/2a) and (\ensuremath{\pi}/a,\ensuremath{\pi}/a), keeping roughly the original band distribution. With hole doping, a flat band develops at the Fermi level with a large Fermi surface and the lower Hubbard band extends through the whole Brillouin zone.