Semiclassical approximation solved by Monte Carlo integration as an efficient impurity solver for dynamical mean field theory and its cluster extensions

Abstract

We propose that a combination of the semiclassical approximation with Monte Carlo simulations can be an efficient and reliable impurity solver for dynamical mean field theory equations and their cluster extensions with large cluster sizes. In order to show the reliability of the method, we consider two test cases: (i) the single-band Hubbard model within the dynamical cluster approximation with 4- and 8-site clusters and (ii) the anisotropic two-orbital Hubbard model with orbitals of different band width within the single-site dynamical mean field theory. We compare our results with those obtained from solving the dynamical mean field equations with continuous time and determinant quantum Monte Carlo. In both test cases we observe reasonable values of the metal-insulator critical interaction strength $U_c/t$ and, while some details of the spectral functions cannot be captured by the semiclassical approximation due to the freezing of dynamical fluctuations, the main features are reproduced by the approach.