Diagrammatic Monte Carlo method for many-polaron problems.

We introduce the first bold diagrammatic Monte Carlo approach to deal with polaron problems at a finite electron density nonperturbatively, i.e., by including vertex corrections to high orders. Using the Holstein model on a square lattice as a prototypical example, we demonstrate that our method is capable of providing accurate results in the thermodynamic limit in all regimes from a renormalized Fermi liquid to a single polaron, across the nonadiabatic region where Fermi and Debye energies are of the same order of magnitude. By accounting for vertex corrections, the accuracy of the theoretical description is increased by orders of magnitude relative to the lowest-order self-consistent Born approximation employed in most studies. We also find that for the electron-phonon coupling typical for real materials, the quasiparticle effective mass increases and the quasiparticle residue decreases with increasing the electron density at constant electron-phonon coupling strength.