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.

The temperature dependence of electronic eigenenergies in the adiabatic harmonic approximation

The renormalization of electronic eigenenergies due to electron-phonon interactions (temperature dependence and zero-point motion effect) is important in many materials. We address it in the adiabatic harmonic approximation, based on first principles (e.g. Density-Functional Theory), from different points of view: directly from atomic position fluctuations or, alternatively, from Janaks theorem generalized to the case where the Helmholtz free energy, including the vibrational entropy, is used. We prove their equivalence, based on the usual form of Janaks theorem and on the dynamical equation. We then also place the Allen-Heine-Cardona (AHC) theory of the renormalization in a first-principle context. The AHC theory relies on the rigid-ion approximation, and naturally leads to a self-energy (Fan) contribution and a Debye-Waller contribution. Such a splitting can also be done for the complete harmonic adiabatic expression, in which the rigid-ion approximation is not required. A numerical study within the Density-Functional Perturbation theory framework allows us to compare the AHC theory with frozen-phonon calculations, with or without the rigid-ion terms. For the two different numerical approaches without rigid-ion terms, the agreement is better than 7 $mu$eV in the case of diamond, which represent an agreement to 5 significant digits. The magnitude of the non rigid-ion terms in this case is also presented, distinguishing specific phonon modes contributions to different electronic eigenenergies.