View source for Transport coefficients including electron-phonon scattering

In the framework of the linearized Boltzman equations, we can compute a few electronic transport observables.
The transport coefficients can be evaluated rather straightforwardly under the approximation of the constant relaxation time.
The most computationally demanding part of the calculation is the electronic linewidths due to the electron-phonon scattering.
As such, it is instructive to start by computing the transport coefficients in the constant-relaxation time approximation (CRTA).
This gives us an idea of the <b>k</b>-point density required to obtain a precise value.

The most basic {{FILE|INCAR}} to compute the transport coefficients is

  {{TAG|PREC}} = Accurate
  {{TAG|EDIFF}} = 1e-8
  {{TAG|ISMEAR}} = -15; SIGMA = 0.01
  {{TAG|LREAL}} = .FALSE.
  {{TAG|LWAVE}} = .FALSE.
  {{TAG|LCHARG}} = .FALSE.
  
  #run electron-phonon calculation
  {{TAG|ELPH_RUN}} = .TRUE.
  {{TAG|ELPH_MODE}} = TRANSPORT
  
  # for the determination of the chemical potential
  {{TAG|EFERMI_NEDOS}} = 101 # number of integration points for {{TAG|ELPH_ISMEAR}} = -24
  {{TAG|ELPH_ISMEAR}} = -24 # smearing method 

  # for the computation of transport coefficients
  {{TAG|TRANSPORT_NEDOS}} = 101
  {{TAG|ELPH_SELFEN_TEMPS}} = 0 100 200 300 400 500

== Conductivity for metals ==
Conductivity <math>\sigma_{\alpha\beta}</math> relates the current to the applied electric field
::<math> J_\alpha = \sigma_{\alpha\beta} E_\beta</math>

=== Constant relaxation time approximation ===
The constant relaxation time is often a good approximation for metals because the density of states often changes little around the Fermi level.
What remains challenging is to determine a value for the constant relaxation time <math>\tau</math>. To avoid this, it is common to report the transport quantities in a different set of units, where <math>\tau=1</math>.
To perform a CRTA calculation, you need to add the following variables to the {{FILE|INCAR}}

  {{TAG|ELPH_SCATTERING_APPROX}} = CRTA
  {{TAG|TRANSPORT_RELAXATION_TIME}} = 1e-14

Note that this calculation can be done without the presence of the {{FILE|phelel_params.hdf5}} that contains the derivatives of the Kohn-Sham potential.

=== Energy relaxation time approximation === 

The <b>k</b>-mesh used in the previous calculation is a good starting point for a computation in the ERTA approximations (or MRTA) since the integral of the transport function is a requirement in both cases. This <b>k</b>-point mesh is, however, not a guarantee of a precise result: the calculation of the linewidth for each Kohn-Sham state requires an integration in a <b>q</b>-point mesh. Currently, it is not possible to choose the two separately, so we are left with the option to increase the sampling in {{FILE|KPOINTS_DENSE}} until the transport quantities are converged.

There are a few factors determining the performance of these computations, in particular the value of {{TAG|TRANSPORT_NEDOS}}.
This tag chooses the number of points to be used in a Gauss-Legendre integral used to integrate the spectral transport function and obtain the Onsager coefficients, which in turn are used to determine the final transport coefficients.
The number of points in turn determines an energy window for which the linewidths will contribute and need to be computed.
This selection of the states for which the linewidths are computed is crucial to attain a nearly linear scaling of the calculation with the number of k-points.
The size of the energy window also depends on the temperature: for larger temperatures the required energy windows are larger, which means that vasp will compute the linewidths for more Kohn-Sham states. The states selected for computation are reported in the {{TAG|OUTCAR}} file

== Mobility for semiconductors ==

In semiconductors, it is more usual to compute the mobility instead of the conductivity.
This is because the conductivity is proportional to the number of free carriers, which means that by adding more free carriers, one obtains a larger conductivity.
By defining the mobility as the conductivity divided by the number of carriers, obtain a more intrinsic property of the material.
For the computation of a realistic mobility, it is usual to compute the conductivity for different doping levels and extrapolate to the limit of low carrier concentration. Additionally, one can specify electron or hole doping. A convenient settings is

  {{TAG|ELPH_SELFEN_CARRIER_DEN}} = -1e17 -1e16 -1e15 -1e14 1e14 1e15 1e16 1e17

=== Constant relaxation time approximation ===
The constant relaxation time is a very crude approximation in the case of semiconductors because it is proportional to the electronic density of states which have a <math>\propto \sqrt{e}</math> behavior around the valence band maximum or conduction band minimum.
It is however still instructive to perform these calculations, such as to get an idea of the <b>k</b>-sampling required to converge the mobility.

Because of the sharp jump in the density of states in the regions of interest for the transport computation, it is crucial to ensure that the integration mesh {{TAG|TRANSPORT_NEDOS}} has enough points. More points are often required for semiconductors than for metals.

=== Energy relaxation time approximation ===

The same observations apply here than in the case of metals, with the particularity that it is more difficult to converge the linewidths for semiconductors. This is specially true in the case of polar materials, where the electron-phonon matrix elements have an integrable divergence when <math>q\rightarrow0</math>.

== Thermoelectric coefficients and the ZT figure of merit ==

The ZT figure of merit is given as a function of the transport coefficients
::<math>ZT=\frac{\sigma S^2}{\kappa_e+\kappa_l}</math>

These different quantities (except for the lattice thermal conductivity <math>\kappa_l</math>) are reported in the OUTCAR file.
The lattice thermal conductivity <math>\kappa_l</math> can be computed with an external package such as [https://phonopy.github.io/phono3py/ phono3py].

Because the Seebeck coefficient <math>S</math> and electronic contribution to the thermal conductivity <math>\kappa_e</math> depend on the first and second momentum of the transport function, a larger number of points ({{TAG|TRANSPORT_NEDOS}}) is often required to obtain an accurate integration of thes transport quantities than for the conductivity or mobility.