The ground state of the half-filled Hubbard model is studied using a real-space renormalization-group technique. We first study in detail the one-dimensional case. We find the ground state to be insulating for all finite coupling, in agreement with the exact solution. We compute the ground-state energy, localization length, energy gap, magnitude of the local moment, and spin-density autocorrelation function. For those quantities that are exactly known we find good agreement with the exact results. Using a simple extension of our one-dimensional calculation, we are able to study approximately two- and three-dimensional lattices. We find a Mott transition at finite interaction for these cases. The critical exponents for these transitions are found to satisfy an approximate interdimensional scaling law.
We study the ground state properties of the doped Hubbard model with strong interactions on the honeycomb lattice by the density matrix renormalization group (DMRG) method. At half-filling, due to the absence of minus sign problem, it is now well established by large-scale quantum Monte Carlo calculations that a Dirac semimetal to antiferromagnetic Mott insulator transition occurs with the increase of the interaction strength
U for the Hubbard model on the honeycomb lattice. However, an understanding of the fate of the antiferromagnetic Mott insulator when holes are doped into the system is still lacking. In this work, by calculating the local spin and charge density for width-4 cylinders with DMRG, we discover a half-filled stripe order in the doped Hubbard model on the honeycomb lattice. We also perform complementary large-scale mean-field calculations with renormalized interaction strength. We observe half-filled stripe order and find stripe states with filling close to one half are nearly degenerate in energy.
The one-band two-dimensional Hubbard Hamiltonian is analyzed using a real-space renormalization-group block method. The renormalization-group method, previously applied to the one-dimensional case has been extended to two dimensions. It is also shown how to avoid the proliferation of new terms during the process if the shape of the blocks and the symmetries of the kept states are chosen conveniently. We characterize the ground state by its energy per site, the gap of charge excitations, the double occupancy, and the effective hopping. The system shows an insulating behavior in the case of half filling for all nonzero interaction parameter, U, if the hopping is isotropic (
ty=
tx); that is to say, the gap opens for all positive U, as in the one-dimensional case. If the hopping anisotropy, \ensuremath{\alpha}=
ty/
tx, is different from the unity and zero the ground state is a conductor up to a critical value (U/t
)c1, which depends on \ensuremath{\alpha}. We have analyzed in the same way another filling, n=7/9. The system behaves as a conductor, paramagnetic up to a value (U/t
)c2\ensuremath{\sim}2 and weakly ferromagnetic for higher values of U.