Dynamical mean field theory of the bilayer Hubbard model with inchworm Monte Carlo

Dynamical mean field theory allows access to the physics of strongly correlated materials with nontrivial orbital structure, but relies on the ability to solve auxiliary multiorbital impurity problems. The most successful approaches to date for solving these impurity problems are the various continuous time quantum Monte Carlo algorithms. Here, we consider perhaps the simplest realization of multiorbital physics: the bilayer Hubbard model on an infinite-coordination Bethe lattice. Despite its simplicity, the majority of this model's phase diagram cannot be predicted by using traditional Monte Carlo methods. We show that these limitations can be largely circumvented by recently introduced inchworm Monte Carlo techniques. We then explore the model's phase diagram at a variety of interaction strengths, temperatures, and filling ratios.