Transport coefficients including electron-phonon scattering: Difference between revisions

Latest revision as of 15:09, 1 October 2024

   Warning: This page is under construction!
   
   This page contains information about the upcoming VASP 6.5.0 release.
   Content on this page is subject to change and may be missing important information.

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 k-point density required to obtain a precise value.

The most basic INCAR to compute the transport coefficients is

 PREC = Accurate
 EDIFF = 1e-8
 ISMEAR = -15; SIGMA = 0.01
 LREAL = .FALSE.
 LWAVE = .FALSE.
 LCHARG = .FALSE.
 
 #run electron-phonon calculation
 ELPH_RUN = .TRUE.
 ELPH_MODE = TRANSPORT
 
 # for the determination of the chemical potential
 EFERMI_NEDOS = 101 # number of integration points for ELPH_ISMEAR = -24
 ELPH_ISMEAR = -24 # smearing method 
 
 # for the computation of transport coefficients
 TRANSPORT_NEDOS = 101
 ELPH_SELFEN_TEMPS = 0 100 200 300 400 500

Conductivity for metals

Conductivity $\sigma _{\alpha \beta }$ relates the current to the applied electric field

J_{\alpha }=\sigma _{\alpha \beta }E_{\beta }

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 $\tau$ . To avoid this, it is common to report the transport quantities in a different set of units, where $\tau =1$ . To perform a CRTA calculation, you need to add the following variables to the INCAR

 ELPH_SCATTERING_APPROX = CRTA
 TRANSPORT_RELAXATION_TIME = 1e-14

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

Self-energy relaxation-time approximation

The k-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 k-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 q-point mesh. Currently, it is not possible to choose the two separately, so we are left with the option to increase the sampling in KPOINTS_ELPH until the transport quantities are converged.

There are a few factors determining the performance of these computations, in particular the value of 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 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.

\sigma =n_{e}\mu _{e}+n_{h}\mu _{h}

with $n_{e}$ being the density of electron carriers, $n_{h}$ the density of hole carriers and $\mu _{e}=\sigma _{e}/n_{e}$ and $\mu _{h}=\sigma _{h}/n_{h}$ with $\sigma _{e}$ being the contribution of the conduction bands to the electron conductivity and $\sigma _{h}$ 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

 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 $\propto {\sqrt {e}}$ 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 k-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 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 $q\rightarrow 0$ .

Thermoelectric coefficients and the ZT figure of merit

The ZT figure of merit is given as a function of the transport coefficients

ZT={\frac {\sigma S^{2}}{\kappa _{e}+\kappa _{l}}}

Both the conductivity $\sigma$ , Seebeck coefficient $S$ and the electronic contribution to the thermal conductivity $\kappa _{e}$ are reported in the OUTCAR file. The lattice thermal conductivity $\kappa _{l}$ can be computed with an external package such as phono3py.

Because the Seebeck coefficient $S$ and electronic contribution to the thermal conductivity $\kappa _{e}$ depend on the first and second momentum of the transport function, a larger number of points (TRANSPORT_NEDOS) is often required to obtain an accurate integration of these transport quantities than for the conductivity or mobility.

@@ Line 1: / Line 1: @@
-In the framework of the linearized Boltzman equations, we can compute a few electronic transport observables.
+{{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).
+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
+   {{TAGBL|PREC}} = Accurate
-   {{TAG|EDIFF}} = 1e-8
+   {{TAGBL|EDIFF}} = 1e-8
-   {{TAG|ISMEAR}} = -15; SIGMA = 0.01
+   {{TAGBL|ISMEAR}} = -15; SIGMA = 0.01
-   {{TAG|LREAL}} = .FALSE.
+   {{TAGBL|LREAL}} = .FALSE.
-   {{TAG|LWAVE}} = .FALSE.
+   {{TAGBL|LWAVE}} = .FALSE.
-   {{TAG|LCHARG}} = .FALSE.
+   {{TAGBL|LCHARG}} = .FALSE.
    #run electron-phonon calculation
-   {{TAG|ELPH_RUN}} = .TRUE.
+   {{TAGBL|ELPH_RUN}} = .TRUE.
-   {{TAG|ELPH_MODE}} = TRANSPORT
+   {{TAGBL|ELPH_MODE}} = TRANSPORT
    # for the determination of the chemical potential
-   {{TAG|EFERMI_NEDOS}} = 101 # number of integration points for {{TAG|ELPH_ISMEAR}} = -24
+   {{TAGBL|EFERMI_NEDOS}} = 101 # number of integration points for {{TAG|ELPH_ISMEAR}} = -24
-   {{TAG|ELPH_ISMEAR}} = -24 # smearing method
+   {{TAGBL|ELPH_ISMEAR}} = -24 # smearing method
    # for the computation of transport coefficients
-   {{TAG|TRANSPORT_NEDOS}} = 101
+   {{TAGBL|TRANSPORT_NEDOS}} = 101
-   {{TAG|ELPH_SELFEN_TEMPS}} = 0 100 200 300 400 500
+   {{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
-=== Constant relaxation time approximation ===
+::<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>.
@@ Line 40: / Line 54: @@
 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 ===
+=== 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 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.
+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}}.
@@ Line 48: / Line 62: @@
 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
+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 ==
@@ Line 55: / Line 69: @@
 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>\mu_{\alpha\beta} = \frac{\sigma_{\alpha\beta}}{n_e}</math>
+::<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
@@ Line 63: / Line 79: @@
 with the negative values indicating hole doping and the positive values, electron doping.
-=== Constant relaxation time approximation ===
+=== 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.
@@ Line 69: / Line 85: @@
 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 ===
+=== 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>.
@@ Line 81: / Line 97: @@
 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.
+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.