New implementation of hybridization expansion quantum impurity solver based on Newton-Leja interpolation polynomial

We introduce a new implementation of hybridization expansion continuous time quantum impurity solver which is relevant to dynamical mean-field theory. It employs Newton interpolation at a sequence of real Leja points to compute the time evolution of the local Hamiltonian efficiently. Since the new algorithm avoids not only computationally expansive matrix-matrix multiplications in conventional implementations but also huge memory consumptions required by Lanczos/Arnoldi iterations in recently developed Krylov subspace approach, it becomes advantageous over the previous algorithms for quantum impurity models with five or more bands. In order to illustrate the great superiority and usefulness of our algorithm, we present realistic dynamical mean-field results for the electronic structures of representative correlated metal SrVO$_{3}$.