Fast multi-orbital equation of motion impurity solver for dynamical mean field theory

We propose a fast multi-orbital impurity solver for the dynamical mean field theory (DMFT). Our DMFT solver is based on the equations of motion (EOM) for local Green's functions and constructed by generalizing from the single-orbital case to the multi-orbital case with inclusion of the inter-orbital hybridizations and applying a mean field approximation to the inter-orbital Coulomb interactions. The two-orbital Hubbard model is studied using this impurity solver within a large range of parameters. The Mott metal-insulator transition and the quasiparticle peak are well described. A comparison of the EOM method with the QMC method is made for the two-orbital Hubbard model and a good agreement is obtained. The developed method hence holds promise as a fast DMFT impurity solver in studies of strongly correlated systems.