Abstract
Non-local, inhomogeneous and retarded response similar to that observed in experiments is studied theoretically by introducing the Inhomogeneous Momentum Average (IMA) approximation for single-polaron problems with disorder in the on-site potential and/or spatial variations of the electron-phonon couplings and/or phonon frequencies. We show that the electron-phonon coupling gives rise to an additional inhomogeneous, strongly retarded potential. This potential describes essential physics ignored by "instanteneous" approximations. The accuracy of IMA is demonstrated by comparison with single-impurity results from the approximation-free Diagrammatic Monte Carlo (DMC) method. Its simplicity allows for easy study of many problems that were previously unaccessible.