Numerical Simulation of a One Dimensional Two Band Model

I. Morgenstern

Hideki Matsumoto, Ferdinando Mancini

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.

Sandro Sorella, Michele Parrinello

Mihir Arjunwadkar, Padmakar V. Panat, D. G. Kanhere

The question of spin-charge separation in two-dimensional lattices has been addressed by numerical simulations of the motion of one hole in a half-filled band. The calculations have been performed on finite clusters with Hubbard and t-J models. By comparing the time evolution of spin and charge polarization currents in one and two dimensions, evidence in favor of spin-charge separation in two dimensions is presented. In contrast with this, spin-charge separation is absent in a highly doped, metallic, system.

Steven R. White, Douglas J. Scalapino, Robert L. Sugar +3 more

We report on a numerical study of the two-dimensional Hubbard model and describe two new algorithms for the simulation of many-electron systems. These algorithms allow one to carry out simulations within the grand canonical ensemble at significantly lower temperatures than had previously been obtained and to calculate ground-state properties with fixed numbers of electrons. We present results for the two-dimensional Hubbard model with half- and quarter-filled bands. Our results support the existence of long-range antiferromagnetic order in the ground state at half-filling and its absence at quarter-filling. Results for the magnetic susceptibility and the momentum occupation along with an upper bound to the spin-wave spectrum are given. We find evidence for an attractive effective d-wave pairing interaction near half-filling but have not found evidence for a phase transition to a superconducting state.

Wenjun Zhang, J.L. Shen