Temperature dependence of the electronic structure of semiconductors and insulators

The renormalization of electronic eigenenergies due to electron-phonon coupling is sizable in many materials with light atoms. This effect, often neglected in ab-initio calculations, can be computed using the perturbation-based Allen-Heine-Cardona theory in the adiabatic or non-adiabatic harmonic approximation. After a short description of the numerous recent progresses in this field, and a brief overview of the theory, we focus on the issue of phonon wavevector sampling convergence, until now poorly understood. Indeed, the renormalization is obtained numerically through a q-point sampling inside the BZ. For q-points close to G, we show that a divergence due to non-zero Born effective charge appears in the electron-phonon matrix elements, leading to a divergence of the integral over the BZ for band extrema. Although it should vanish for non-polar materials, unphysical residual Born effective charges are usually present in ab-initio calculations. Here, we propose a solution that improves the coupled q-point convergence dramatically. For polar materials, the problem is more severe: the divergence of the integral does not disappear in the adiabatic harmonic approximation, but only in the non-adiabatic harmonic approximation. In all cases, we study in detail the convergence behavior of the renormalization as the q-point sampling goes to infinity and the imaginary broadening parameter goes to zero. This allows extrapolation, thus enabling a systematic way to converge the renormalization for both polar and non-polar materials. Finally, the adiabatic and non-adiabatic theory, with corrections for the divergence problem, are applied to the study of five semiconductors and insulators: a-AlN, b-AlN, BN, diamond and silicon. For these five materials, we present the zero-point renormalization, temperature dependence, phonon-induced lifetime broadening and the renormalized electronic bandstructure.