Category:Electron-phonon interactions

In most of the features implemented in VASP, the electron and vibrational degrees of freedom (phonons) are treated separately. Electron-phonon coupling refers to the coupling of these two systems, and can be used to include the effect of each other when treating the two systems separately.

The inclusion of the effects of the ionic degrees of freedom in the electronic structure is important in the determination of many physical observables such as the bandgap, spectral functions, electronic conductivity, Seebeck coefficient or electronic thermal conductivity to name a few.

With VASP, one can use a few different approaches to include these effects: by statistically averaging observables during a molecular dynamics calculation, by stochastically sampling the ionic degrees of freedom or using many-body perturbation theory. In this page we outline the two later approaches, their advantages and disadvantages and their range of application.

Theoretical background

For a very extensive overview of the theory of electron-phonon coupling, we recommend reading the Rev. Mod. Phys. Here we outline a few equations that are important to describe the features implemented in VASP.

The stochastic approach allows obtaining the bandgap renormalization and an approximation of the electronic spectral function due to the ionic degrees of freedom under the static approximation using a supercell approach. This has the advantage that it can be easily implemented and used with different levels of theory to describe the electronic states, such as different exchange correlation functional or even the GW approximation.

Another approach to include the coupling between the electronic and ionic degrees of freedom is by using the methods and language of many-body perturbation theory, where the coupling is included as a perturbation of the electronic or phononic states. In the case of the perturbation of the electronic states, we can access the bandgap renormalization as well as electronic transport coefficients with the inclusion of phonon scattering which is dominant in a wide range of applications such as mobility of semiconductors or conductivity of metals at room temperature.

How to

The stochastic approach implies the knowledge of the phonons on a supercell, these are then used to displace the atoms along a set of random or a single specially chosen direction [cite ZG] with the displacement length determined by the temperature of the ionic system. The bandgap can be directly obtained by averaging over the set of randomly displaced supercells, or directly from a supercell with the ZG displacements.

• Electron-phonon interactions from Monte-Carlo sampling

When treating the electron-phonon interaction as a perturbation the central quantity that we need to compute becomes the electron-phonon matrix element which within the framework of DFT requires the knowledge of the change of Kohn-Sham potential with an ionic perturbation as well as the initial and final electronic Kohn-Sham states. This requires the calculation to be done in two steps: first, the derivatives of the KS potential are computed on a representative supercell and then the electron-phonon matrix elements are computed on a dense `k`-point grid.

• Zero-point renormalization from a perturbation theory approach

The transport coefficients are computed within the framework of the linearized Boltzmann transport equation, which in turn requires the knowledge of the electronic lifetimes due to every possible scattering event. There are a few approximations for these scattering rates: in the constant-relaxation time approximation (CRTA), an empirical parameter is used to approximate the scattering rate throughout the Brillouin zone, in the self-energy relaxation time approximation (SERTA) or energy relaxation time approximation (ERTA), the lifetimes are approximated from the electron self-energy and finally in the momentum relaxation time approximation (MRTA) an additional geometric factor is included to account for the correct balance between inwards and outwards scattering.

• Phonon-limited electron transport coefficients