Written by Michael Y. Toriyama (MichaelToriyama2024 [at] u [dot] northwestern [dot] edu)
Thermoelectric properties in a material can often be described using Boltzmann transport theory under low applied field. The electrical and thermal current densities can be described in terms of the applied electric field and a thermal gradient. The relevant charge transport properties, namely the electrical conductivity, Seebeck coefficient, and the electronic contribution to the thermal conductivity, can be derived from electrical and thermal current densities.
Thermoelectric Performance
The efficiency of a thermoelectric device, from a material engineering perspective, is parametrized by the figure-of-merit , which can be expressed as
where is the electrical conductivity, is the Seebeck coefficient, and are the electronic and lattice contributions to the thermal conductivity, and is the temperature.
Single Parabolic Band Model
The single parabolic band assumption forms the simplest, first-order description of a material's electronic structure. Often, this is a reasonable approximation in regards to charge transport properties.
We do not delve into the derivation of the transport parameters here. You can refer to the theory section of for more details. Rather, we list the final forms of the transport properties, namely the electrical conductivity (), Seebeck coefficient (), and Lorenz number ():
The electronic thermal conductivity can be expressed as
Fundamental constants are needed to describe transport; namely, is the electric charge, is the free electron mass, is the Boltzmann constant, and is Planck's constant.
The transport properties are all dependent on the reduced Fermi level () and the scattering parameter (). The reduced Fermi level is given as
where is the Fermi level and is the temperature. is referenced to the band edge; in other words, at the band edge, when the Fermi level is in the band, and when the Fermi level is in the band gap. This description holds regardless of whether the majority carriers are electrons (conduction band) or holes (valence band).
is a parameter representing the scattering mechanism of charge carriers. Some common scattering mechanisms are:
Acoustic phonon scattering:
Polar optical phonon scattering:
Ionized impurity scattering:
The Fermi-Dirac integral is given as
Using the Fermi-Dirac integral, the carrier concentration can also be expressed as
The electrical conductivity is also dependent on the weighted mobility , given as
where is the effective band mass, and is the valley degeneracy.
Given that the transport properties can be written in closed forms, we can express in a more compartmentalized way:
where is the quality factor of the majority carriers defined as
Notice that , as expressed this way, is determined by two principal quantities of the material: the intrinsic property , and an extrinsic factor .
The two-band model described here involves the transport of both electrons and holes. In other words, information about the conduction and valence bands are necessary, as well as the band gap. It is convenient to represent the electron and hole contributions separately (using the single parabolic band model of transport as described above) and derive the total electrical conductivity, Seebeck coefficient, and electronic thermal conductivity from the two contributions.
If we reference the Fermi level to the conduction band minimum (i.e. at the CBM), then the reduced Fermi levels of electrons () and holes () are given as
where is the band gap. The form of can be justified as follows: if is in the conduction band (i.e. ), then from the perspective of holes, is far below the valence band edge (i.e. ). On the other hand, if is below the conduction band edge, then is at the valence band edge, and therefore .
In general, the reference energy for is arbitrary, meaning that can be defined as the CBM, VBM, mid-gap, etc. It is more convenient to define this reference point by the majority carrier type in the material. For example, can be defined as the CBM when the material is n-type, whereas it can be defined as the VBM when the material is p-type. The mathematical formalism for the transport properties can be generalized by using the "majority" and "minority" carrier labels; in other words, we can express the reduced Fermi levels as
where at the edge of the majority carrier band. The thermoelectric transport properties can then be written as
As before, the electronic thermal conductivity is given as .
The important aspect of the two-band model to keep in mind is that the model can be built from the individual contributions of electrons and holes. Therefore, a rational way to write a script for the two-band model would be to write functions relevant for the single-band model and to subsequently use those functions to build the two-band model.
Computational Details
Necessary Modules
Scientific computing can be handled solely by numpy. The Fermi-Dirac integral can be evaluated using the module.
From fdint, we use the function fdk(k=i, phi=eta) to evaluate Fermi integrals of the form .
Fundamental Constants
We need to define fundamental constants (in SI units) to describe transport properties.
Transport Properties for a Single Band
As suggested in the , we build the transport model from two single-band transport models (one for the valence band, and the other for the conduction band). The transport properties (, , and ) within a single-band model can be codified as follows:
Here, eta is the reduced Fermi level andr is the scattering parameter. If needed, the electronic thermal conductivity () can be expressed straightforwardly as
For example, if we set the following parameters for our model system:
we obtain the following plots:
Transport Properties in a Two-Band Model
When both holes and electrons contribute to thermoelectric transport, the transport coefficients need to be evaluated using the two-band model.
The electrical conductivity is simply the sum of the conductivities of the electrons and holes:
Note that Sigma_SPB is the electrical conductivity from a single parabolic band.
The Seebeck coefficient is a sum of the Seebeck coefficients of the electrons and holes, weighted by their respective conductivities:
Note that Seebeck_SPB is the Seebeck coefficient from a single parabolic band.
The electronic thermal conductivity can be codified as
Note that Lorenz_SPB is the Lorenz number from a single parabolic band.
If we set the following parameters in our model system:
then we obtain the following plots:
