Coexistence of solutions in dynamical mean-field theory of the Mott transition

In this paper, I discuss the finite-temperature metal-insulator transition of the paramagnetic Hubbard model within dynamical mean-field theory. I show that coexisting solutions, the hallmark of such a transition, can be obtained in a consistent way both from Quantum Monte Carlo (QMC) simulations and from the Exact Diagonalization method. I pay special attention to discretization errors within QMC. These errors explain why it is difficult to obtain the solutions by QMC close to the boundaries of the coexistence region.