Algorithm for computing perturbation series of dynamical mean field theory

We show how to use diagrammatic techniques to compute the weak-coupling perturbation series of the self-consistent solution to a dynamical mean field theory (DMFT) problem. This approach constitutes an alternative to using diagrammatic techniques directly as an impurity solver. It allows one to bypass the need of multiple perturbative series resummations within the DMFT self-consistency loop. It can be applied at or out of equilibrium, with any diagrammatic formalism, such as real times, imaginary times, or Matsubara frequencies formalisms. As a proof of principle, we illustrate our method with the half-filled Hubbard model on the Bethe lattice in the DMFT approximation, using quantum quasi Monte Carlo (QQMC) to obtain the impurity perturbation series on the real time axis.