Seebeck coefficient explained

The Seebeck coefficient (also known as thermopower,[1] thermoelectric power, and thermoelectric sensitivity) of a material is a measure of the magnitude of an induced thermoelectric voltage in response to a temperature difference across that material, as induced by the Seebeck effect.[2] The SI unit of the Seebeck coefficient is volts per kelvin (V/K),[2] although it is more often given in microvolts per kelvin (μV/K).

The use of materials with a high Seebeck coefficient[3] is one of many important factors for the efficient behaviour of thermoelectric generators and thermoelectric coolers. More information about high-performance thermoelectric materials can be found in the Thermoelectric materials article. In thermocouples the Seebeck effect is used to measure temperatures, and for accuracy it is desirable to use materials with a Seebeck coefficient that is stable over time.

Physically, the magnitude and sign of the Seebeck coefficient can be approximately understood as being given by the entropy per unit charge carried by electrical currents in the material. It may be positive or negative. In conductors that can be understood in terms of independently moving, nearly-free charge carriers, the Seebeck coefficient is negative for negatively charged carriers (such as electrons), and positive for positively charged carriers (such as electron holes).

Definition

One way to define the Seebeck coefficient is the voltage built up when a small temperature gradient is applied to a material, and when the material has come to a steady state where the current density is zero everywhere. If the temperature difference ΔT between the two ends of a material is small, then the Seebeck coefficient of a material is defined as:

S=-{\DeltaV\over\DeltaT}

where ΔV is the thermoelectric voltage seen at the terminals. (See below for more on the signs of ΔV and ΔT.)

Note that the voltage shift expressed by the Seebeck effect cannot be measured directly, since the measured voltage (by attaching a voltmeter) contains an additional voltage contribution, due to the temperature gradient and Seebeck effect in the measurement leads. The voltmeter voltage is always dependent on relative Seebeck coefficients among the various materials involved.

Most generally and technically, the Seebeck coefficient is defined in terms of the portion of electric current driven by temperature gradients, as in the vector differential equation

J=-\sigma\boldsymbol\nablaV-\sigmaS\boldsymbol\nablaT

where

\scriptstyleJ

is the current density,

\scriptstyle\sigma

is the electrical conductivity,

\scriptstyle\boldsymbol\nablaV

is the voltage gradient, and

\scriptstyle\boldsymbol\nablaT

is the temperature gradient. The zero-current, steady state special case described above has

\scriptstyleJ=0

, which implies that the two electrical conductivity terms have cancelled out and so

\boldsymbol\nablaV=-S\boldsymbol\nablaT.

Sign convention

The sign is made explicit in the following expression:

S=-

V\rm-V\rm
T\rm-T\rm
Thus, if S is positive, the end with the higher temperature has the lower voltage, and vice versa. The voltage gradient in the material will point against the temperature gradient.

The Seebeck effect is generally dominated by the contribution from charge carrier diffusion (see below) which tends to push charge carriers towards the cold side of the material until a compensating voltage has built up. As a result, in p-type semiconductors (which have only positive mobile charges, electron holes), S is positive. Likewise, in n-type semiconductors (which have only negative mobile charges, electrons), S is negative. In most conductors, however, the charge carriers exhibit both hole-like and electron-like behaviour and the sign of S usually depends on which of them predominates.

Relationship to other thermoelectric coefficients

\scriptstyle\Pi

by the exact relation

S=

\Pi
T

,

where

T

is the thermodynamic temperature.

\scriptstylelK

by

S=\int

lK
T

dT.

The constant of integration is such that

\scriptstyleS=0

at absolute zero, as required by Nernst's theorem.

Measurement

Relative Seebeck coefficient

See also: Thermocouple. In practice the absolute Seebeck coefficient is difficult to measure directly, since the voltage output of a thermoelectric circuit, as measured by a voltmeter, only depends on differences of Seebeck coefficients. This is because electrodes attached to a voltmeter must be placed onto the material in order to measure the thermoelectric voltage. The temperature gradient then also typically induces a thermoelectric voltage across one leg of the measurement electrodes. Therefore, the measured Seebeck coefficient is a contribution from the Seebeck coefficient of the material of interest and the material of the measurement electrodes. This arrangement of two materials is usually called a thermocouple.

The measured Seebeck coefficient is then a contribution from both and can be written as:

SAB=SB-SA={\DeltaVB\over\DeltaT}-{\DeltaVA\over\DeltaT}

Absolute Seebeck coefficient

Although only relative Seebeck coefficients are important for externally measured voltages, the absolute Seebeck coefficient can be important for other effects where voltage is measured indirectly. Determination of the absolute Seebeck coefficient therefore requires more complicated techniques and is more difficult, but such measurements have been performed on standard materials. These measurements only had to be performed once for all time, and for all materials; for any other material, the absolute Seebeck coefficient can be obtained by performing a relative Seebeck coefficient measurement against a standard material.

A measurement of the Thomson coefficient

l{K}

, which expresses the strength of the Thomson effect, can be used to yield the absolute Seebeck coefficient through the relation:

S(T)=

T
\int
0

{l{K}(T')\overT'}dT'

, provided that

l{K}

is measured down to absolute zero. The reason this works is that

S(T)

is expected to decrease to zero as the temperature is brought to zero—a consequence of Nernst's theorem. Such a measurement based on the integration of

l{K}/T

was published in 1932,[4] though it relied on the interpolation of the Thomson coefficient in certain regions of temperature.

Superconductors have zero Seebeck coefficient, as mentioned below. By making one of the wires in a thermocouple superconducting, it is possible to get a direct measurement of the absolute Seebeck coefficient of the other wire, since it alone determines the measured voltage from the entire thermocouple. A publication in 1958 used this technique to measure the absolute Seebeck coefficient of lead between 7.2 K and 18 K, thereby filling in an important gap in the previous 1932 experiment mentioned above.[5]

The combination of the superconductor-thermocouple technique up to 18 K, with the Thomson-coefficient-integration technique above 18 K, allowed determination of the absolute Seebeck coefficient of lead up to room temperature. By proxy, these measurements led to the determination of absolute Seebeck coefficients for all materials, even up to higher temperatures, by a combination of Thomson coefficient integrations and thermocouple circuits.[6]

The difficulty of these measurements, and the rarity of reproducing experiments, lends some degree of uncertainty to the absolute thermoelectric scale thus obtained. In particular, the 1932 measurements may have incorrectly measured the Thomson coefficient over the range 20 K to 50 K. Since nearly all subsequent publications relied on those measurements, this would mean that all of the commonly used values of absolute Seebeck coefficient (including those shown in the figures) are too low by about 0.3 μV/K, for all temperatures above 50 K.[7]

Seebeck coefficients for some common materials

In the table below are Seebeck coefficients at room temperature for some common, nonexotic materials, measured relative to platinum.[8] The Seebeck coefficient of platinum itself is approximately −5 μV/K at room temperature,[9] and so the values listed below should be compensated accordingly. For example, the Seebeck coefficients of Cu, Ag, Au are 1.5 μV/K, and of Al −1.5 μV/K. The Seebeck coefficient of semiconductors very much depends on doping, with generally positive values for p doped materials and negative values for n doping.

Material Seebeck coefficient
relative to platinum (μV/K)
900
500
440
330
47
25
19
10
7.5
6.5
6.0
4.5
4.0
3.5
3.0
0.6
0 (definition)
-2.0
-9.0
-15
-35
-72

Physical factors that determine the Seebeck coefficient

A material's temperature, crystal structure, and impurities influence the value of thermoelectric coefficients. The Seebeck effect can be attributed to two things:[10] charge-carrier diffusion and phonon drag.

Charge carrier diffusion

On a fundamental level, an applied voltage difference refers to a difference in the thermodynamic chemical potential of charge carriers, and the direction of the current under a voltage difference is determined by the universal thermodynamic process in which (given equal temperatures) particles flow from high chemical potential to low chemical potential. In other words, the direction of the current in Ohm's law is determined via the thermodynamic arrow of time (the difference in chemical potential could be exploited to produce work, but is instead dissipated as heat which increases entropy). On the other hand, for the Seebeck effect not even the sign of the current can be predicted from thermodynamics, and so to understand the origin of the Seebeck coefficient it is necessary to understand the microscopic physics.

Charge carriers (such as thermally excited electrons) constantly diffuse around inside a conductive material. Due to thermal fluctuations, some of these charge carriers travel with a higher energy than average, and some with a lower energy. When no voltage differences or temperature differences are applied, the carrier diffusion perfectly balances out and so on average one sees no current:

\scriptstyleJ=0

. A net current can be generated by applying a voltage difference (Ohm's law), or by applying a temperature difference (Seebeck effect). To understand the microscopic origin of the thermoelectric effect, it is useful to first describe the microscopic mechanism of the normal Ohm's law electrical conductance—to describe what determines the

\scriptstyle\sigma

in

\scriptstyleJ=-\sigma\boldsymbol\nablaV

. Microscopically, what is happening in Ohm's law is that higher energy levels have a higher concentration of carriers per state, on the side with higher chemical potential. For each interval of energy, the carriers tend to diffuse and spread into the area of device where there are fewer carriers per state of that energy. As they move, however, they occasionally scatter dissipatively, which re-randomizes their energy according to the local temperature and chemical potential. This dissipation empties out the carriers from these higher energy states, allowing more to diffuse in. The combination of diffusion and dissipation favours an overall drift of the charge carriers towards the side of the material where they have a lower chemical potential.[11]

For the thermoelectric effect, now, consider the case of uniform voltage (uniform chemical potential) with a temperature gradient. In this case, at the hotter side of the material there is more variation in the energies of the charge carriers, compared to the colder side. This means that high energy levels have a higher carrier occupation per state on the hotter side, but also the hotter side has a lower occupation per state at lower energy levels. As before, the high-energy carriers diffuse away from the hot end, and produce entropy by drifting towards the cold end of the device. However, there is a competing process: at the same time low-energy carriers are drawn back towards the hot end of the device. Though these processes both generate entropy, they work against each other in terms of charge current, and so a net current only occurs if one of these drifts is stronger than the other. The net current is given by

\scriptstyleJ=-\sigmaS\boldsymbol\nablaT

, where (as shown below) the thermoelectric coefficient

\scriptstyle\sigmaS

depends literally on how conductive high-energy carriers are, compared to low-energy carriers. The distinction may be due to a difference in rate of scattering, a difference in speeds, a difference in density of states, or a combination of these effects.

Mott formula

The processes described above apply in materials where each charge carrier sees an essentially static environment so that its motion can be described independently from other carriers, and independent of other dynamics (such as phonons). In particular, in electronic materials with weak electron-electron interactions, weak electron-phonon interactions, etc. it can be shown in general that the linear response conductance is

\sigma=\intc(E)(-

df(E)
dE

)dE,

and the linear response thermoelectric coefficient is

\sigmaS=

k\rm
-e

\int

E-\mu
k\rmT

c(E)(-

df(E)
dE

)dE

where

\scriptstylec(E)

is the energy-dependent conductivity, and

\scriptstylef(E)

is the Fermi–Dirac distribution function. These equations are known as the Mott relations, of Sir Nevill Francis Mott.[12] The derivative
-df(E)
dE

=

1
4k\rmT

\operatorname{sech}2\left(

E-\mu
2k\rmT

\right)

is a function peaked around the chemical potential (Fermi level)

\mu

with a width of approximately

3.5k\rmT

. The energy-dependent conductivity (a quantity that cannot actually be directly measured — one only measures

\sigma

) is calculated as

c(E)=e2D(E)\nu(E)

where

D(E)

is the electron diffusion constant and

\nu(E)

is the electronic density of states (in general, both are functions of energy).

In materials with strong interactions, none of the above equations can be used since it is not possible to consider each charge carrier as a separate entity. The Wiedemann–Franz law can also be exactly derived using the non-interacting electron picture, and so in materials where the Wiedemann–Franz law fails (such as superconductors), the Mott relations also generally tend to fail.[13]

The formulae above can be simplified in a couple of important limiting cases:

Mott formula in metals

In semimetals and metals, where transport only occurs near the Fermi level and

\scriptstylec(E)

changes slowly in the range

E\mu\pmk\rmT

, one can perform a Sommerfeld expansion

\scriptstylec(E)=c(\mu)+c'(\mu)(E-\mu)+O[(E-\mu)2]

, which leads to

S\rm=

2
\pi
2
k
\rmB
T
-3e
c'(\mu)
c(\mu)

+O[(k\rmT)3],\sigma\rm=c(\mu)+O[(k\rmT)2].

This expression is sometimes called "the Mott formula", however it is much less general than Mott's original formula expressed above.

In the free electron model with scattering, the value of

\scriptstylec'(\mu)/c(\mu)

is of order

\scriptstyle1/(k\rmT\rm)

, where

T\rm

is the Fermi temperature, and so a typical value of the Seebeck coefficient in the Fermi gas is

\scriptstyleS\rm\tfrac{\pi2k\rm

} T/T_ (the prefactor varies somewhat depending on details such as dimensionality and scattering). In highly conductive metals the Fermi temperatures are typically around 104 – 105 K, and so it is understandable why their absolute Seebeck coefficients are only of order 1 – 10 μV/K at room temperature. Note that whereas the free electron model predicts a negative Seebeck coefficient, real metals actually have complicated band structures and may exhibit positive Seebeck coefficients (examples: Cu, Ag, Au).

The fraction

\scriptstylec'(\mu)/c(\mu)

in semimetals is sometimes calculated from the measured derivative of

\scriptstyle\sigma\rm

with respect to some energy shift induced by field effect. This is not necessarily correct and the estimate of

\scriptstylec'(\mu)/c(\mu)

can be incorrect (by a factor of two or more), since the disorder potential depends on screening which also changes with field effect.[14]
Mott formula in semiconductors

In semiconductors at low levels of doping, transport only occurs far away from the Fermi level. At low doping in the conduction band (where

\scriptstyleE\rm-\mu\ggk\rmT

, where

\scriptstyleE\rm

is the minimum energy of the conduction band edge), one has

\scriptstyle-

df(E)
dE

\tfrac{1}{k\rmT}

-(E-\mu)/(k\rmT)
e
. Approximating the conduction band levels' conductivity function as

\scriptstylec(E)=A\rm(E-E\rm

a\rm
)

for some constants

\scriptstyleA\rm

and

\scriptstylea\rm

,

S\rm=

k\rm
-e

[

E\rm-\mu
k\rmT

+a\rm+1],\sigma\rm=A\rm(k\rm

a\rm
T)
-E\rm-\mu
k\rmT
e

\Gamma(a\rm+1).

whereas in the valence band when

\scriptstyle\mu-E\rm\ggkT

and

\scriptstylec(E)=A\rm(E\rm-

a\rm
E)

,

S\rm=

k
e

[

\mu-E\rm
k\rmT

+a\rm+1],\sigma\rm=A\rm(k\rm

a\rm
T)
-\mu-E\rm
k\rmT
e

\Gamma(a\rm+1).

The values of

\scriptstylea\rm

and

\scriptstylea\rm

depend on material details; in bulk semiconductor these constants range between 1 and 3, the extremes corresponding to acoustic-mode lattice scattering and ionized-impurity scattering.[15]

In extrinsic (doped) semiconductors either the conduction or valence band will dominate transport, and so one of the numbers above will give the measured values. In general however the semiconductor may also be intrinsic in which case the bands conduct in parallel, and so the measured values will be

S\rm=

\sigma\rmS\rm+\sigma\rmS\rm
\sigma\rm+\sigma\rm

,\sigma\rm=\sigma\rm+\sigma\rm

\scriptstyle\sigmaS2

,[16] or the thermoelectric figure of merit, and the optimum generally occurs at high doping levels.[17]

Phonon drag

See main article: Phonon drag. Phonons are not always in local thermal equilibrium; they move against the thermal gradient. They lose momentum by interacting with electrons (or other carriers) and imperfections in the crystal. If the phonon-electron interaction is predominant, the phonons will tend to push the electrons to one end of the material, hence losing momentum and contributing to the thermoelectric field. This contribution is most important in the temperature region where phonon-electron scattering is predominant. This happens for

T{1\over5}\thetaD

where

\scriptstyle\theta\rm

is the Debye temperature. At lower temperatures there are fewer phonons available for drag, and at higher temperatures they tend to lose momentum in phonon-phonon scattering instead of phonon-electron scattering. At lower temperatures, material boundaries also play an increasing role as the phonons can travel significant distances.[18] Practically speaking, phonon drag is an important effect in semiconductors near room temperature (even though well above

\scriptstyle\theta\rm/5

), that is comparable in magnitude to the carrier-diffusion effect described in the previous section.[18]

This region of the thermopower-versus-temperature function is highly variable under a magnetic field.

Relationship with entropy

The Seebeck coefficient of a material corresponds thermodynamically to the amount of entropy "dragged along" by the flow of charge inside a material; it is in some sense the entropy per unit charge in the material.[19]

Notes and References

  1. Thermopower is a misnomer as this quantity does not actually express a power quantity: Note that the unit of thermopower (V/K) is different from the unit of power (watts).
  2. Book: Blundell . Stephen . Blundell . Stephen J. . Blundell . Katherine M. . 2010 . Concepts in Thermal Physics . Oxford University Press . 978-0-19-956210-7 . 415 .
  3. 10.1002/anie.200900598 . 19866458 . New and Old Concepts in Thermoelectric Materials . Joseph R. Sootsman . Duck Young Chung . Mercouri G. Kanatzidis . 2009 . Angewandte Chemie. 48 . 46 . 8616–8639 .
  4. Borelius. G.. Keesom. W. H.. Johannson. C. H.. Linde. J. O.. 1932. Establishment of an Absolute Scale for the Thermo-electric Force. Proceedings of the Royal Academy of Sciences at Amsterdam. 35. 1. 10.
  5. Christian . J. W. . Jan . J.-P. . Pearson . W. B. . Templeton . I. M. . Thermo-Electricity at Low Temperatures. VI. A Redetermination of the Absolute Scale of Thermo-Electric Power of Lead . 10.1098/rspa.1958.0078 . Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences . 245 . 1241 . 213 . 1958 . 1958RSPSA.245..213C . 96708128 .
  6. Cusack . N. . Kendall . P. . 10.1088/0370-1328/72/5/429 . The Absolute Scale of Thermoelectric Power at High Temperature . Proceedings of the Physical Society . 72 . 5 . 898 . 1958 . 1958PPS....72..898C .
  7. Roberts . R. B. . Absolute scales for thermoelectricity . 10.1016/0263-2241(86)90016-3 . Measurement . 4 . 3 . 101–103 . 1986 . 1986Meas....4..101R .
  8. http://www.electronics-cooling.com/2006/11/the-seebeck-coefficient/ The Seebeck Coefficient, Electronics Cooling.com (accessed 2013-Feb-01)
  9. Moore . J. P. . Absolute Seebeck coefficient of platinum from 80 to 340 K and the thermal and electrical conductivities of lead from 80 to 400 K . 10.1063/1.1662324 . Journal of Applied Physics . 44 . 3 . 1174–1178 . 1973 . 1973JAP....44.1174M .
  10. Book: Kong. Ling Bing. Waste Energy Harvesting. 24. Springer. 978-3-642-54634-1. 263–403. 10.1007/978-3-642-54634-1. Lecture Notes in Energy. 2014.
  11. Book: 978-0-521-63145-7 . Quantum Transport: Atom to Transistor . Datta . Supriyo . Supriyo Datta . 2005 . Cambridge University Press .
  12. Cutler . M. . Mott . N. . 10.1103/PhysRev.181.1336 . Observation of Anderson Localization in an Electron Gas . Physical Review . 181 . 3 . 1336 . 1969 . 1969PhRv..181.1336C .
  13. Jonson . M. . Mahan . G. . 10.1103/PhysRevB.21.4223 . Mott's formula for the thermopower and the Wiedemann-Franz law . Physical Review B . 21 . 10 . 4223 . 1980 . 1980PhRvB..21.4223J .
  14. Hwang . E. H. . Rossi . E. . Das Sarma . S. . Theory of thermopower in two-dimensional graphene . 10.1103/PhysRevB.80.235415 . Physical Review B . 80 . 23 . 235415 . 2009 . 0902.1749 . 2009PhRvB..80w5415H . 8125966 .
  15. Semiconductor Physics: An Introduction, Karlheinz Seeger
  16. Imai . H. . Shimakawa . Y. . Kubo . Y. . Large thermoelectric power factor in TiS2 crystal with nearly stoichiometric composition . Physical Review B . 10 December 2001 . 64 . 24 . 241104 . 10.1103/PhysRevB.64.241104 . cond-mat/0111063 . 2001PhRvB..64x1104I . 119389373 .
  17. G. Jeffrey Snyder, "Thermoelectrics". http://www.its.caltech.edu/~jsnyder/thermoelectrics/
  18. Mahan . G. D. . Lindsay . L. . Broido . D. A. . The Seebeck coefficient and phonon drag in silicon . Journal of Applied Physics . 28 December 2014 . 116 . 24 . 245102 . 10.1063/1.4904925 . 2014JAP...116x5102M . 1185754 .
  19. Bulusu . A. . Walker . D. G. . 10.1016/j.spmi.2008.02.008 . Review of electronic transport models for thermoelectric materials . Superlattices and Microstructures . 44 . 1 . 1 . 2008 . 2008SuMi...44....1B .