Nonlinear electron-phonon interactions from first principles

Abstract

Electron-phonon interactions underpin a variety of phenomena, ranging from transport and superconductivity to polarons and ultrafast carrier dynamics. Despite being one of the most intensely studied subjects in condensed matter physics, research on electron-phonon physics mostly focused on linear, first-order couplings. Second- and higher-order nonlinear couplings are commonly ignored because their calculations are too demanding and we lack computational frameworks that can provide both diagonal and off-diagonal coupling matrix elements. In this work, we report a theory and computational method for computing nonlinear electron-phonon interactions of any order in real materials. Our approach combines the advantages of unit-cell calculations of electron wavefunctions and supercell calculations of phonon perturbations, is systematically improvable, can be used with either semilocal or nonlocal exchange-correlation functionals, and is amenable to Wannier-Fourier interpolation. As a first proof of concept, we illustrate this method by computing second-order electron-phonon coupling matrix elements in diamond, lithium fluoride, and graphite as representative nonpolar semiconductors, polar semiconductors, and metals, respectively. Furthermore, we generalize the ab initio polaron equations to second-order electron-phonon couplings, and we show that second-order couplings are essential to achieve quantitative accuracy in polaron formation energies and hopping barriers. The present methodology will find application in the study of all properties and phenomena that are currently being investigated within the linear electron-phonon coupling approximation, from phonon-mediated superconductivity to excited-states dynamics, both in harmonic and anharmonic systems.