Numerical solution of thed=∞ Hubbard model: Evidence for a Mott transition

We present a numerical solution of the infinite-dimensional Hubbard model at finite temperature in the paramagnetic phase. The problem reduces to a single-impurity Anderson model supplemented by a self-consistency condition. Using Monte Carlo methods and complete enumeration we determine the imaginary-time Green's function, the density of doubly occupied sites, and the compressibility close to half filling. All three quantities present direct evidence for a Mott insulating phase above a critical value of U.