  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.
