We examine the range of validity of the second-order cumulant expansion (CE) for the calculation of spectral functions, quasiparticle properties, and mobility of the Holstein polaron. We devise an efficient numerical implementation that allows us to make comparisons in a broad interval of temperature, electron-phonon coupling, and phonon frequency. For a benchmark, we use the dynamical mean-field theory which gives, as we have recently shown, rather accurate spectral functions in the whole parameter space, even in low dimensions. We find that in one dimension, the CE resolves well both the quasiparticle and the first satellite peak in a regime of intermediate coupling. At high temperatures, the charge mobility assumes a power law $\ensuremath{\mu}\ensuremath{\propto}{T}^{\ensuremath{-}2}$ in the limit of weak coupling and $\ensuremath{\mu}\ensuremath{\propto}{T}^{\ensuremath{-}3/2}$ for stronger coupling. We find that, for stronger coupling, the CE gives slightly better results than the self-consistent Migdal approximation (SCMA), while the one-shot Migdal approximation is appropriate only for a very weak electron-phonon interaction. We also analyze the atomic limit and the spectral sum rules. We derive an analytical expression for the moments in CE and find that they are exact up to the fourth order, as opposed to the SCMA where they are exact to the third order. Finally, we analyze the results in higher dimensions.
Read full abstract