View source for Transport coefficients including electron-phonon scattering

{{elph release}}

In the framework of the linearized Boltzmann 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

  {{TAGBL|PREC}} = Accurate
  {{TAGBL|EDIFF}} = 1e-8
  {{TAGBL|ISMEAR}} = -15; SIGMA = 0.01
  {{TAGBL|LREAL}} = .FALSE.
  {{TAGBL|LWAVE}} = .FALSE.
  {{TAGBL|LCHARG}} = .FALSE.
  
  #run electron-phonon calculation
  {{TAGBL|ELPH_RUN}} = .TRUE.
  {{TAGBL|ELPH_MODE}} = TRANSPORT
  
  # for the determination of the chemical potential
  {{TAGBL|EFERMI_NEDOS}} = 101 # number of integration points for {{TAG|ELPH_ISMEAR}} = -24
  {{TAGBL|ELPH_ISMEAR}} = -24 # smearing method 
  
  # for the computation of transport coefficients
  {{TAGBL|TRANSPORT_NEDOS}} = 101
  {{TAGBL|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>
<!--The conductivity is computed from the Onsager coefficient

::<math> \quad \sigma = \mathcal{L}_{11} </math>

which in turn is defined as
::<math>
    \mathcal{L}_{ij}= \int_{-\infty}^{+\infty} d \epsilon \, \frac{\partial  f }{\partial \mu} \Big( \frac{\epsilon-\mu}{-e}\Big)^{i+j-2} \mathcal{T}(\epsilon) </math>

with <math>f</math> the fermi dirac ocupation and <math>\mathcal{T}(\epsilon)</math> the transport function defined as
::<math>
 \mathcal{T}(\epsilon)=\frac{e^2}{V N} \sum_{\mathbf{k}n} \tau_{\mathbf{k}n} \mathbf{v}_{\mathbf{k}n} \otimes \mathbf{v}_{\mathbf{k}n} \delta(\epsilon-\epsilon_{\mathbf{k}n}) </math>

with <math>V</math> the volume of the lattice, and N the number of k-point in the full Brillouin zone.-->
=== 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.

=== Self-energy relaxation-time approximation === 

The <b>k</b>-mesh used in the previous calculation is a good starting point for a computation in the SERTA 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_ELPH}} 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.
::<math>\sigma = n_e \mu_e + n_h \mu_h</math>
with <math>n_e</math> being the density of electron carriers, <math>n_h</math> the density of hole carriers and <math>\mu_e = \sigma_e/n_e</math> and <math>\mu_h = \sigma_h/n_h</math> with <math>\sigma_e</math> being the contribution of the conduction bands to the electron conductivity and <math>\sigma_h</math> of the valence bands.

For the computation of a realistic mobility, it is usual to compute the conductivity for different doping levels and extrapolate to the limit of the lowest possible carrier concentration. The value of the mobility should be relatively constant in the limit of low carrier concentrations.
Additionally, one can specify electron or hole doping. A convenient setting is

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

with the negative values indicating hole doping and the positive values, electron doping.

=== 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.

=== Self-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>

Both the conductivity <math>\sigma</math>, Seebeck coefficient <math>S</math> and the electronic contribution to the thermal conductivity <math>\kappa_e</math> are reported in the {{FILE|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 these transport quantities than for the conductivity or mobility.