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.

• Received 20 February 1992

DOI:https://doi.org/10.1103/PhysRevLett.69.1240