Beyond the Fröhlich Hamiltonian: Path integral treatment of large polarons in anharmonic solids

The properties of an electron in a typical solid are modified by the interaction with the crystal ions, leading to the formation of a quasiparticle: the polaron. Such polarons are often described using the Fröhlich Hamiltonian, which assumes the underlying lattice phonons to be harmonic. However, this approximation is invalid in several interesting materials, including the recently discovered high-pressure hydrides which superconduct at temperatures above 200K. In this paper, we show that Fröhlich theory can be extended to eliminate this problem. We derive four additional terms in the Fröhlich Hamiltonian to account for anharmonicity up to third order. We calculate the energy and effective mass of the new polaron, using both perturbation theory and Feynman's path integral formalism. It is shown that the anharmonic terms lead to significant additional trapping of the electron. The derived Hamiltonian is well-suited for analytical calculations, due to its simplicity and since the number of model parameters is low. Since it is a direct extension of the Fröhlich Hamiltonian, it can readily be used to investigate the effect of anharmonicity on other polaron properties, such as the optical conductivity and the formation of bipolarons.