We develop a method for calculating the electron-phonon vertex in polar semiconductors and insulators from first principles. The present formalism generalizes the Fröhlich vertex to the case of anisotropic materials and multiple phonon branches, and can be used either as a post-processing correction to standard electron-phonon calculations, or in conjunction with ab initio interpolation based on maximally localized Wannier functions. We demonstrate this formalism by investigating the electron-phonon interactions in anatase TiO2, and show that the polar vertex significantly reduces the electron lifetimes and enhances the anisotropy of the coupling. The present work enables ab initio calculations of carrier mobilities, lifetimes, mass enhancement, and pairing in polar materials.PACS numbers: 71.38.-k, 63.20.dk The electron-phonon interaction (EPI) is a cornerstone of condensed matter physics, and plays important roles in a diverse array of phenomena. Recent years have witnessed a surge of interest in ab initio calculations of EPIs, leading to new techniques and many innovative applications in the case of metals and non-polar semiconductors [1][2][3][4][5][6][7][8][9][10][11][12]. In contrast to this fast-paced progress, in the case of polar semiconductors and insulators the study of EPIs from first principles has not gone very far, owing to the prohibitive computational costs of EPI calculations for polar materials. For example, a fully ab initio calculation of the carrier mobility of a polar semiconductor has not been performed yet, while such calculations have recently been reported for non-polar semiconductors such as silicon [13] and graphene [14]. Given the fast-growing technological importance of polar semiconductors, from light-emitting devices to transparent electronics, solar cells and photocatalysts [15][16][17], developing accurate and efficient computational methods for studying EPIs in these systems is of primary importance.At variance with metals and non-polar semiconductors, in polar materials two or more atoms in the unit cell carry nonzero Born effective charge tensors [18]. As a consequence, the fluctuations of the ionic positions corresponding to longitudinal optical (LO) phonons at long wavelength generate macroscopic electric fields which can couple strongly to electrons and holes, leading to the so-called Fröhlich interaction [19]. Up to now this interaction has not been taken into account in ab initio calculations of EPIs; the two key obstacles towards a description of Fröhlich coupling from first principles are (i) the Fröhlich coupling was designed to describe simple isotropic systems with one LO phonon, and (ii) the electron-phonon vertex diverges for q → 0, where q is the phonon wavevector. The first obstacle relates to the fundamental question on how to define the Fröhlich coupling in the most general way. The second obstacle renders first-principles calculations extremely demanding, since a correct description of the singularity requires a very fine sampling of the Brillouin zone.In...