We present real-time inchworm quantum Monte Carlo results for single-site dynamical mean field theory on an infinite coordination number Bethe lattice. Our numerically exact results are obtained on the L-shaped Keldysh contour and, being evaluated in real time, avoid the analytic continuation issues typically encountered in Monte Carlo calculations. Our results show that inchworm Monte Carlo methods have now reached a state where they can be used as dynamical mean field impurity solvers and the dynamical sign problem can be overcome. As nonequilibrium problems can be simulated at the same cost, we envisage the main use of these methods as dynamical mean field solvers for time-dependent problems far from equilibrium.
We present clear numerical evidence for the coexistence of metallic and insulating dynamical mean-field theory (DMFT) solutions in a half-filled single-band Hubbard model at finite temperature. The quantum Monte Carlo (QMC) method is used to solve the DMFT equation. We discuss important technical aspects of the DMFT-QMC method, which should be properly dealt with in order to obtain reliable DMFT solutions at low temperatures. Among these technical aspects are the critical slowing down near the boundaries of the coexistence region, the convergence criteria for the DMFT iterations, and the correct interpolation of the discretized Green's functions.
Dynamical mean-field theory (DMFT) is one of the most widely-used methods to treat accurately electron correlation effects in ab-initio real material calculations. Many modern large-scale implementations of DMFT in electronic structure codes involve solving a quantum impurity model with a continuous-time quantum Monte Carlo (CT-QMC) solver [Rubtsov et al., Phys. Rev. B 72, 035122 (2005); Werner et al., Phys. Rev. Lett. 97, 076405 (2006); Werner and Millis, Phys. Rev. B 74, 155107 (2006); Gull et al., Rev. Mod. Phys. 83, 349 (2011)]. The main advantage of CT-QMC is that, unlike standard quantum Monte Carlo approaches, it is able to generate the local Green's functions $G(\ensuremath{\tau})$ of the correlated system on an arbitrarily fine imaginary time $\ensuremath{\tau}$ grid, and is free of any systematic errors. In this work, we extend a hybrid QMC solver proposed by Khatami et al. [Phys. Rev. E 81, 056703 (2010)] and Rost et al. [Phys. Rev. E 87, 053305 (2013)] to a multiorbital context. This has the advantage of enabling impurity solver QMC calculations to scale linearly with inverse temperature $\ensuremath{\beta}$, and permit its application to
d- and
f-band materials. In addition, we present a Green's-function processing scheme which generates accurate quasicontinuous imaginary time solutions of the impurity problem which overcome errors inherent to standard QMC approaches. This solver and processing scheme are incorporated into a full DFT+DMFT calculation using the castep DFT code [Clark et al., Z. Kristallogr. 220, 567 (2005)]. Benchmark calculations for
SrVO3 properties are presented. The computational efficiency of this method is also demonstrated.