We report on a non-perturbative approach to the 1D and 2D Hubbard models that is capable of recovering both strong and weak-coupling limits. We first show that even when the on-site Coulomb repulsion, U, is much smaller than the bandwith, the Mott-Hubbard gap never closes at half-filling in both 1D and 2D. Consequently, the Hubbard model at half-filling is always in the strong-coupling non-perturbative regime. For both large and small U, we find that the population of nearest-neighbour singlet states approaches a value of order unity as $T\to 0$ as would be expected for antiferromagnetic order. We also find that the double occupancy is a smooth monotonic function of U and approaches the anticipated non-interacting limit and large U limits. Finally, in our results for the heat capacity in 1D differ by no more than 1% from the Bethe ansatz predictions. In addition, we find that in 2D, the heat capacity vs T for different values of U exhibits a universal crossing point at two characteristic temperatures as is seen experimentally in a wide range of strongly-correlated systems such as $^3He$, $UBe_3$, and $CeCu_{6-x}Al_x$. The success of this method in recovering well-established results that stem fundamentally from the Coulomb interaction suggests that local dynamics are at the heart of the physics of strongly correlated systems.
Read full abstract