Open access peer-reviewed chapter

Thermoelectric Effect and Application of Organic Semiconductors

By Nianduan Lu, Ling Li and Ming Liu

Submitted: April 6th 2016Reviewed: September 20th 2016Published: December 21st 2016

DOI: 10.5772/65872

Downloaded: 1098


Human development and society progress require solving many pressing issues, including sustainable energy production and environmental conservation. Thermoelectric power generation looks like promising opportunity converting huge heat from the sun and waste heat from industrial sector, housing appliances and infrastructure and automobile and other fuel combustion exhaust directly to electrical energy. Thermoelectric power generation will be of high demand, when technology will be affordable, providing low price, high conversion efficiency, reliability, easy applicability and advanced ecological properties of end products. In this context, organic thermoelectric materials attract great interest caused by non-scarcity of raw materials, non-toxicity, potentially low costs in high-scale production, low thermal conductivity and wide capabilities to control thermoelectric properties. In this chapter, we focus mainly on thermoelectric effect in several organic semiconductors, both crystalline and disordered. We present theory of some transport phenomena determining thermoelectric properties of organic semiconductors, including general expression of thermoelectric effect, percolation theory of Seebeck coefficient, hybrid model of Seebeck coefficient, Monte Carlo simulation and first-principle theory. Finally, a future outlook of this field is briefly discussed.


  • organic semiconductor
  • thermoelectric effect
  • theoretical model

1. Introduction

1.1. Organic semiconductors

1.1.1. History

Organic semiconductors have revealed promising performance and received considerable attentions due to large area, low-end, lightweight and flexible electronics applications [1]. Currently, organic semiconductors have appealed for a broad range of devices including sensors, solar cells, light-emitting devices and thermoelectric application [2, 3]. Historically, organic materials were viewed as insulators with applications commonly seen in inactive packaging, coating, containers, moldings and so on. The earliest research on electrical behavior of organic materials dates back to the 1960s [4]. In the 1970s, photoconductive organic materials were recognized and were used in solar cells and xerographic sensors, etc. [5]. In about the same age, organic thermoelectric materials, such as conducting polymers, had been investigated [6], although First European Conference on Thermoelectric in 1988 had no mention of organic materials in their proceeding. Since proof of concept for organic semiconductors occurred in the 1980s, remarkable development of organic semiconductors has promoted improvement of performance that maybe competitive with amorphous silicon (a-Si), increasing their suitability for commercial applications [7]. Appearance of conductive polymers in the late 1970s, and of conjugated semiconductors and photoemission polymers in the 1980s, greatly accelerated development in the field of organic electronics [8]. Polyacetylene was one of the earliest polymer materials known to be potential as conducting electricity [9], and one could find that oxidative dopant with iodine could greatly increase conductivity by 12 orders of magnitude [10]. This discovery and development of highly conductive organic materials were attributed to three scientists: Alan J. Heeger, Alan G. MacDiarmid and Hideki Shirakawa, who were jointly awarded the Nobel Prize in Chemistry in 2000 for their discovery in 1977 and development of oxidized, iodine-doped polyacetylene. After that, plenty of organic semiconductor materials were synthesized, and research field of organic electronics matured over the years from proof-of-principle phase into major interdisciplinary research area, involving physics, chemistry and other disciplines. As important branch of organic semiconductors, in the past decade, organic thermoelectric effect has received much attention. Figure 1 shows Thomson Reuters Web of Science publication report for the topic “organic thermoelectric Seebeck effect” for the last 16 years [3]. Research interest in organic thermoeletrics Seebeck effect has been growing remarkably over the last 5 years.

Figure 1.

Thomson Reuters web of science publication report for the topic “organic thermoelectric Seebeck effect” from 2000 to 2015.

1.1.2. Structure

Organic semiconductors can be usually classified as two types: crystalline and amorphous materials, in terms of crystalline fraction (also static disorder) [3, 11].

All organic semiconductor materials are generally characterized by weak van der Waals bonding, which leads to weak intermolecular interactions. This weak coupling of molecules would induce weak interaction energy to give narrow electronic bandwidths. Otherwise, existing narrow electronic bands will be eliminated by statistical variation of width in energy level distribution of molecules, which, hence, creates Anderson charge localization. For crystalline organic semiconductor materials, localization of charge carriers is attributed to intermolecular thermal fluctuations (dynamic disorder), at which size of localized wave function is expected to be on the order of molecular spacing and charge carrier transport in this weakly localized field is treated as “intermediate hopping transport regime” [12]. Because of high density of crystal imperfections in disordered organic semiconductors, such as impurities, grain boundaries, dangling bonds and periodicity loss of crystal, localization of charge carriers is attributed to spatial and energetic disorder due to weak intermolecular interactions [13]. Disorder in organic semiconductors results in basic charge transport mechanism, common for very rich variety of such materials [14], incoherent tunneling (hopping) of charge carriers between localized states [15].

1.2. Charge transport mechanism in organic semiconductors

1.2.1. Dispersive transport

Charge carriers are always slowing down during conduction processes in dispersive transport regime. It happens only when charge carriers’ distribution is thermally nonequilibrium. In dispersive transport regime, energy fluctuations allow release (emission) of charge carrier captured by localized state (trap), when state (trap) becomes temporarily shallow during fluctuations [16]. Otherwise, charge carrier release time is strongly dependent on both temperature and energy. Charge carriers localized on shallow states are usually released before the states can change their energies noticeably. In contrast, charge carriers, always being trapped in energetically deeper and deeper states, have to perform longer and longer tunneling transitions to hop to the next destination site. As a result, temperature-dependent distribution of effective activation energies just follows the density of states (DOS) function within domain of shallow states, while these distributions appear to be very different for deep traps.

1.2.2. Nearest-neighbor hopping

Transport in disordered organic semiconductors is generally characterized by charge carrier’s localization and hopping transport mechanism [17]. Localized states, randomly distributed in energy and space, form discrete array of sites in hopping space. The most probable hop for charge carrier on site with particular energy is to the closest empty site, that is, to its nearest-neighbor site in hopping space. In conjunction with hopping probability rate, it gives mobility for carriers at this energy. In a word, nearest-neighbor hopping describes hopping regime, in which tunneling part of hopping rates in Eq. (1) referred as hopping rates of Miller-Abrahams is so much slower than energy contribution, that only the nearest neighbors are addressed in hops [18]:


where ν0is attempt-to-jump frequency, Ris hopping range, αis inverse localized length describing extension of wave function of localized state, Rijis distance between site iand site j, Eiand Ejare energies of sites iand j, respectively. As long as charge carrier can find shallow and empty sites with energies below its current state, it will perform nearest-neighbor hopping to energetically lower sites, since in this case rates are limited by spatial tunneling distances only.

1.2.3. Variable-range hopping

Variable-range hopping (VRH) theory was first proposed by Neville Mott in 1971 [19] and hence was called Mott VRH, which is model describing low temperature conduction in strongly disordered systems with localized states. VRH transport has characteristic temperature dependence of 3D electrical conductance:


here kBT0=βg(Ef)×α3, σ0is prefactor, α1is localization length, kBis Boltzmann constant, g(Ef)is DOS function at Fermi energy Ef, and βis constant coefficient with value in interval 10.0–37.2 according to different theories.

For 3D electrical conductance, and in general for d-dimensions, VRH transport is expressed as:


here d is dimensionality.

1.2.4. Multiple trapping and release theory

Multiple trapping and release theory assume that charge carrier transport occurs in extended states, and that most of charge carriers are trapped in localized states. Energy of localized state is separated from mobility edge energy. When localized state energy is slightly lower mobility edge, then localized state acts as shallow trap, from which charge carrier can be released (emitted) by thermal excitations. But, if that energy is far below mobility edge energy, then charge carrier cannot be thermally excited (emitted). In multiple trapping and release theory, total charge carriers’ concentration ntotalis equal to sum of concentrations nein extended states and in localized states as in Ref. [20]:


where f(E)is Fermi-Dirac distribution function.

2. Organic thermoelectric materials

2.1. Polymer-based thermoelectric materials

Polymers as thermoelectric materials recently have attracted much attention due to easy fabrication processes and low material cost [21, 22], as well as, low thermal conductivity, which is highly desirable for thermoelectric applications. Different types of polymers have been used in thermoelectric devices [2325], such as polyaniline (PANI), poly(p-phenylenevinylene) (PPV), poly(3,4-ethylenedioxythiophene) (PEDOT), tosylate(tos), poly(styrenesulfonate) (PSS), and poly(2,5-dimethoxy phenylenevinylene) (PMeOPV), poly[2-methoxy-5-(2-ethylhexyloxy)-1,4-phenylenevinylene] (MEHPPV) and poly(3-hexylthiophene-2,5-diyl) (P3HT). Figure 2 shows chemical structure of some polymers. These polymers are chosen as thermoelectric materials for their conductive nature. Takao Ishida’s group has reported high power factor (PF) values of over 100 μW/(m K2) on PEDOT films through chemical reduction with different chemicals in their review [25], as shown in Table 1.

MaterialsS (μV/K)PF (μW/m K2)ZT
PEDOT:tos (dedoped)2003240.25
PEDOT:PSS (dedoped)~501120.093
PEDOT:PSS (dedoped)431160.2

Table 1.

Thermoelectric properties of polymer-based materials (tos: tosylate, and PSS: poly(stryrenesulfonate) [25].

Figure 2.

Molecular structures of some typical polymers [23].

Thermoelectric properties of 1,1,2,2-ethenetetrathiolate(ett)–metal coordination polymers poly[Ax(M–ett)] (A = Na, K; M = Ni, Cu) also have been studied, as shown in Figure 3a [26]. P-type poly[Cux(Cu–ett)] exhibited best ZT of 0.014 at 380 K with electrical conductivity of ~15 S/cm, Seebeck coefficient of 80 μV/K and thermal conductivity of 0.45 W/(m K); n-type poly[Kx(Ni–ett)] showed best ZT of 0.2 at 440 K with electrical conductivity of ~60 S/cm, Seebeck coefficient of 150 μV/K and thermal conductivity of 0.25 W/(m K). Otherwise, thermoelectric module based on p-type poly[Cux(Cu–ett)] and n-type poly[Nax(Ni–ett)] (i.e., ZT of 0.1 at 440 K) was built (Figure 3b and c ). More recently, Pipe et al. reported thermoelectric measurements of PEDOT: PSS with 5% of dimethylsulfoxide (DMSO) and ethylene glycol (EG) after submerging films in EG several times (from 0 to 450 min) [27]. Insulating polymer (PSS) is removed, and consequently, electrical conductivity and Seebeck coefficient increase simultaneously. ZT value of 0.42 reported in Pipe et al.’s work is the highest ever obtained for polymer until date. Although thermoelectric performance of organic semiconductors has been increased by using different fabricated methods or dopant, as compared with inorganic thermoelectric materials, organic thermoelectric materials still exhibit lower ZT so far. Nevertheless, researchers are making their great efforts in organic semiconductors instead of inorganic materials, due to several more advantages in organic thermoelectric materials, for example, the non-scarcity of raw materials, non-toxicity and large area applications.

Figure 3.

(a) Synthetic route of poly[Ax(M-ett)]s, (b) module structure and (c) photograph of module and measurement system with a hot plane and cooling fan.

2.2. Small molecule-based thermoelectric materials

Currently, investigation of small molecule-based organic thermoelectric materials is lagging behind that of polymers-based organic materials. However, small molecule-based organic thermoelectric materials exhibit plenty of attractions. For example, small molecules are easier to be purified and crystallized and may be more feasible to achieve n-type conduction. Figure 4 shows chemical structures of some small molecules studied as organic thermoelectric materials [24].

Figure 4.

Chemical structures of some small molecules studied as organic thermoelectric materials.

As the benchmark of p-type organic semiconductors, pentacene is the best-known small molecule studied for organic thermoelectric applications. A remarkable attraction for pentacene is attributed to its high mobility up to 3 cm2/V s in thin film transistors (TFT). Due to low charge carriers’ concentration of pentacene in neutral state, appropriate doping treatment is critical to optimization of its thermoelectric properties. So far, F4TCNQ and iodine are efficient dopant for pentacene [24]. F4TCNQ is strong electron acceptor that is frequently used as p-type dopant in organic electronics [28]. Harada et al. have achieved maximum electrical conductivity of 4.1 × 10−2 S/cm and maximum PF of 0.16 μW/(m K2) by using dopant of 2.0 mol.%, as shown in Figure 5a. Furthermore, when thickness of pentacene layer was varied, electrical conductivity could be optimized, while Seebeck coefficient was unaffected (around 200 μV/K). Finally, maximum PF of 2.0 μW/(m K2) was obtained in 6-nm-thick pentacene sample, as shown in Figure 5b. For using iodine dopant, Minakata et al. have achieved highest electrical conductivity of 60 S/cm with Seebeck coefficient in the range of 40–60 μV/K [29]. As a result, the highest PF is 13 μW/(m K2), which is more than six times that of pentacene/F4TCNQ bilayer sample. In a word, as compared with polymer-based thermoelectric materials, small molecule-based materials are relatively less explored and need to put in more effort in the future.

Figure 5.

Thermoelectric properties of pentacene samples doped by (a) co-evaporation with F4TCNQ and (b) forming pentacene/F4TCNQ bilayer structures.

3. Thermoelectric transport theory of organic semiconductors

3.1. General expression of thermoelectric effect

Since organic semiconductors consist of amorphous and crystal structure, theoretical model of charge carrier thermoelectric transport should be more general. Derivation of present expression of thermoelectric effect was inspired by work of Cutler and Mott [30]. Basic expression was in terms of electrical conductivity σ. Based on definition of Cutler and Mott, for hopping system in disordered lattice at zero and finite temperature, σwas expressed as:


Otherwise, one can start out by writing electrical conductivity as integral over single states neglecting electron correlation effects [31]:


Then, energy dependence of electrical conductivity is written as:


here g(E)is density of states, μ(E)is mobility, and f(E)is Fermi-Dirac distribution function. Seebeck coefficient Sis related to Peltier coefficient Πas follows:


Peltier coefficient Πis energy carried by electrons per unit charge. Carried energy is characterized connected with Fermi energy Ef. Contribution to Πof each electron is in proportion to its correlative contribution to total conduction. Weighting factor for electrons in interval dEat energy Eis thus σ(E)σdE, with energy dependence of conductivity σ(E)as Eq. (7). Therefore, one can obtain general expression of Seebeck coefficient as [32]:


To distinguish crystalline solids, general Seebeck coefficient also can be defined as shape of transport energy with mean energy of conducting charge carriers as in Ref. [33]:


where transport energy Etransis defined as averaged energy weighted by electrical conductivity distribution:


Usual expression of Seebeck coefficient was used early in doped organic materials by Roland Schmechel in 2003 [33]. In Schmechel’s article, a detailed method to complex hopping transport in doped system (p-doped zinc-phthalocyanine) was proposed and used to discuss experimental data on effect of doping on conductivity, mobility and Seebeck coefficient.

3.2. Percolation theory of Seebeck coefficient

Percolation theory is considered as the best way known to analytically investigate charge carriers hopping transport characteristics. Percolation problem for charge carriers transport properties in disordered semiconductors has been argued early by Ziman and a number of researchers [34, 35], at which charge transport should be in proportion to percolation probability P(p). A simple definition is that approximates firstly electrical conductivity as [36]:


where P(p)percolation probability is fraction of the volume allowed, but not isolated, and σ0denotes a large allowed value of the material. P(p)is known to vanish for pless than critical value Bc, but drops sharply to zero as pBc. To make percolation question simple, researchers have put forward two kinds of standpoints for critical Bc, that is, Bc=1, and Bc=2.8or Bc=2.7. Although researchers have not achieved unified agreement, percolation theory in hopping system was generally established to explain charge carrier transport characteristics.

Generally speaking, charge carrier transport is described by a four-dimensional (4D) hopping space, including three spatial coordinates and one energy coordinate, at which probability of charge carrier hopping between localized sites associated with these four coordinates. Therefore, charge carrier transport would be more complex based on percolation approach addressing, if all of positions, that is, spatial positions of sites and their energies and occupation of sites, must be included. In Eq. (12), to simulate electrical conductivity σ(E), the key is to seek out percolation path in hopping space. Thus, a random resistor network linking all of molecular sites under percolation model is essential. Figure 6 shows schematic diagram of charge carrier transport in hopping system and corresponding percolation current through polymer matrix for charge carrier to travel through [37].

Figure 6.

(a) Schematic diagram of charge carrier transport in hopping space with density of states and (b) corresponding percolation current in disordered organic semiconductor.

Based on the following general definition through Kelvin-Onsager relation to Peltier coefficient Π, percolation theory to calculate Seebeck coefficient Sin hopping transport is expressed as in Eq. (8), where Πis generally identified with average site energy on percolation cluster and can be written as:


where P(Ei)is probability that site of energy Eiis on current-carrying percolation cluster and was further given by:


where g(Ei)is density of states per unit volume, Emis maximum site energy, and P1(Zm|Ei)is probability that second smallest resistance emanating from site with energy Eiis not larger than maximum resistance on percolation cluster Zm. Expression of probability P1(Zm|Ei)is written as:


where P(Zm|Ei)is bonds’ density, which means average number of resistance of Zmor less connected to site energy Ei. To calculate Peltier coefficient (or Seebeck coefficient), an expression for P(Zm|Ei)is necessary.

Based on percolation theory, disordered organic material is regarded as a random resistor network (see Figure 6b). To address total electrical conductivity in disordered system, the initial is to obtain reference conductance Zand eliminate all conductive pathways between sites iand jwith Zij<Z. Conductance between sites iand jis given by Zijexp(Sij)with [38]:


The density of bonds P(Zm|Ei)then can be written as:


If the density of participating sites is Ps, then critical parameter Scis found by solving equation:


Based on numerical studies for three-dimensional amorphous system, the formation of infinite cluster corresponds to Bc=2.7[38]. By connecting Eqs. (16)–(18), bonds’ density can be formulated as:


here is normalized energy as =EkBT. This expression has been split into two regimes of i>fand i<f, which are corresponding to contributions of iabove or below Fermi energy to P(Zm|Ei)and, therefore, Seebeck coefficient S, respectively. Above Fermi energy, charge carriers in shallow states will move by hopping to other shallow states. While below Fermi energy, charge carriers in deep states will move by thermal excitation to shallower states. By substituting Eqs. (13)–(15) and Eq. (19) into Eq. (8), one can obtain the final result of Seebeck coefficient.

Figure 7 shows simulated and experimental dependences of Seebeck coefficient on charge carriers’ density; simulation is based on percolation theory, experimental data measured by using field-effect transistor (FET) from three kinds of conjugated polymers, that is, IDTBT, PBTTT and PSeDPPBT [39]. Model of percolation theory can reasonably reproduce experimental data under the whole range of charge carriers’ density for different conjugated polymers.

Figure 7.

Charge carrier’s density dependence of Seebeck coefficient for different materials at room temperature [37]. Symbols and solid lines are experimental and simulated results, respectively.

3.3. Hybrid model of Seebeck coefficient

Usual behavior of Seebeck coefficient is to decrease with increasing charge carriers’ density. However, Germs et al. have observed unusual thermoelectric behavior for pentacene in TFT [40], indeed, at room temperature, increasing charge carriers’ density results in expected decrease in S, while with decreasing temperature to values below room temperature, Sdemonstrates growth with increasing charge carriers’ density at T = 250 K and even more pronounced at T = 200 K, as shown in Figure 8.

Figure 8.

Measurements (symbols) and calculation (lines) of Seebeck coefficient as function of gate bias in pentacene TFT. Gate voltage Vg is corrected for threshold voltage Vth of TFT.

To explain this unusual thermoelectric behavior, Germs et al. developed simplified hybrid model that incorporates both variable-range hopping (VRH) and mobility edge (ME) transport [40]. Charge carrier and energy transport can be described independently by two processes: VRH-type process that occurs within exponential tail of localized states and band-like type transport that occurs within band-like states above mobility edge. Then, Seebeck coefficient of hybrid model is expressed as conductivity-weighted average of two contributing transport channels:


where SMEand σMEare Seebeck coefficient and electrical conductivity in ME part, and SVRHand σVRHare Seebeck coefficient and electrical conductivity in VRH part, respectively.

Then, general expression of Seebeck coefficient in Eq. (9) reduces:




where ε=EEc, Ecis energy value at mobility edge.

In Eq. (21), A accounts for excitations beyond the band edge with 1–20% of SME. Similarly, within VRH model, it is assumed that transport is determined by hopping event from equilibrium energy state to relatively narrow transport energy E[41], and Eq. (9) becomes:


Electrical conductivity in ME part is calculated as:


with power law dependence on temperature, μfree(T)=μ0Tm.

VRH part is described by Mott-Martens model that assumes transport to be dominated by hops from Fermi energy to transport level E. Electrical conductivity in VRH part is subsequently calculated by optimizing Miller-Abrahams expression as:


where position of transport level Eand typical hopping distance Ris connected via percolation argument:


where Bc=2.8is critical number of bonds, g(E)represents DOS, which here is simplified to single exponential trap tail below mobility edge and constant density of extended states above Ec:


where n0is divided by kBT0for dimensional reasons. The number of charge carriers above Ec, nfreefollows from Fermi-Dirac distribution.

Figure 9 shows measured and calculated dependences of Seebeck coefficient on difference between gate voltage Vgand threshold voltage Vthon TFT at T = 200 K. One can see that at 200 K heat transported at mobility edge (EcEf)and heat transported at transport level (EEf)both decrease with increasing charge carriers’ density, accounting for downward trends in SMEand SVRH. Consequently, weight-averaged Seebeck coefficient SHybshifts from SVRHvalues at small gate bias up toward SMEfor large gate bias. Relatively large value of SMEat lower temperatures follows from Eq. (21) and temperature independence of EC.

Figure 9.

Measured (symbols) and calculated (lines) dependences of Seebeck coefficient on gate bias at T = 200 K [40]. SHyb is conductivity-weighted average of SVRH and SME.

3.4. Monte Carlo simulation

As compared with numeric model, analytical thermoelectric transport models exhibit more context and physical property, but they also have inevitable shortcoming due to the use of plenty of free parameters during simulation and calculation. In order to eliminate these hindrances, universal method has been used based on Monte Carlo (MC) simulation for describing hopping transport and insuring validity and accuracy of thermoelectric transport.

Kinetic Monte Carlo simulation generally includes six steps as follows [42]:

  1. Initializing site energy Ei. Random energy Eiat site iderives from Gaussian or exponential DOS.

  2. Initial placement of charge carriers. Fermi-Dirac occupation probabilities will be used in randomly placing Nchcharge carriers.

  3. Choosing hopping events. If neglecting the polaron effect, hopping rates νшоfrom site ito site jare based on Miller-Abrahams expression in Eq. (1). Corresponding setting is: (a) hopping rates equal to 0 to prevent hopping into already occupied sites; (b) introducing cut-off distance and set νшо=0for jumps longer than this distance.

    Otherwise, renormalizing hopping rates Γijas pij=Γiji,jΓij. Sum of rates includes only jumps from occupied to unoccupied sites, that is, Γij=0for occupied site jor unoccupied site i. For every pair ij, index k(i.e., ijkand pij=pk), where k{1,...,kmax}with kmaxbeing total number of all possible hopping events. Then, partial sum Skis defined for every index k:


    Apparently, for every k, extent from interval [Sk1,Sk]equals to probability pkfor kthjump, and total extent of all intervals equals to 1, for example Skmax=1. Then, determining a random real number rfrom interval [0, 1] and finding index k, here Sk1<<r<<Sk, and one can find hopping event. After determining hoping event, one would move charge carrier between corresponding sites iand j.

  4. Calculation of waiting time. After finding each hopping transport, total simulation time tand waiting time τwould be added that has passed until the event took place. This time is determined by describing a random number from exponential waiting time distribution P(τ)=νiexp(νiτ)with νi=jνijbeing the total rate for charge carrier hopping from site i. It is, therefore, written as:


    where random number xis drawn from interval [0, 1].

  5. Calculating current density. Every time, when predefined numbers of jumps have occurred, current density J(t)can be expressed as:


    where N+and Nare the total number of jumps in and opposite direction of electric field for cross-sectional slice in yzplane, and a is lattice constant.

  6. Calculating Seebeck coefficient. Seebeck coefficient is given by expression as in Eq. (10), where transport energy is defined as averaged energy weighted by electrical conductivity distribution:


    with σ(T)=σ(E,T)(fE)dE.

Although kinetic MC technique gives a direct simulation of thermoelectric transport in organic semiconductor materials and, therefore, it is accurate for the most description of electronic conductivity, its negative factor is to require extensive computational resources, which leads to difficultly analyze and fit experimental data. Figure 10 shows comparison of analytical model with MC simulation for Seebeck coefficient [42]. It is seen in Figure 10, that results exhibit qualitative agreement for all values of parameter α. For large αin Figure 10a, Ssaand SMCshow not only qualitative, but even relatively good quantitative agreement in energy interval E<0(corresponding to relative charge carriers’ concentration n/N0<0.5). For higher energies (and thus for higher concentrations), difference between Ssaand SMCincreases. As parameter αdecreases, functional dependencies Ssaand SMCremain very close to each other, but Ssagets shifted with respect to SMC(Figure 10b).

Figure 10.

Monte Carlo and semi-analytical calculations of Seebeck coefficient for different values of localization length (a) α−1 = 1 nm and (b) α−1 = 0.2 nm. Ef is in units of Δ, Δ=4kBT and T = 300 K.

3.5. First-principle calculation theory

First-principle (ab initio) theory would be deemed to the best type of theory for hopping charge transport in organic semiconductors, since it starts from particular chemical and geometrical structure of the system, and it starts directly at the level of established science and does not make assumptions, such as empirical model and fitting parameters. So far, a few researches on charge carriers transport properties based on first-principle theory are hardly beyond the scope of crystalline. Otherwise, direct calculation of Seebeck coefficient is hardly realistic. Current method combines generally first-principle calculations with transport theory. For example, Gao et al. have investigated theoretically Seebeck coefficient of narrow bandgap crystalline polymers, including crystalline solids β-Zn4Sb3 and AuIn2 and these polymers, based on muffin-tin orbital and full-potential linearized augmented plane-wave (FLAPW) electronic structure code [43]. In essence, Gao et al.’s method for calculation of transport properties of crystalline solid is firstly based on semiclassical Boltzmann theory, following as:


where e, τ, fand νrepresent free electron charge, electronic relaxation time, Fermi-Dirac distribution function and Fermi velocity, respectively. If relaxation time for electron scattering processes is assumed to be constant, that is, τ(E,T)=const, which may yield reasonable simulated results in a wide range of materials, then temperature dependence of σ0(T)can be simulated based on constant relaxation time τ:


Then, Seebeck coefficient is calculated from ratio of the zeroth and first moments of electrical conductivity with respect to energy:




Product of the density of states N(E)and arbitrary quantity g, which is relative to energy and k-vector as in Eqs. (33) and (35), can be calculated by using integration on constant energy surface Sin kspace. Electronic band structure can be calculated by using WIEN97 and WIEN2K FLAPW [44] or pseudopotential plane-wave code in Vienna ab initio simulation package (VASP) [45].

Figure 11a and b shows simulated energy band structure of polythiophene polymer. Here, internal structural parameters of the polymer are fully optimized, and electronic band structure is calculated by pseudopotential plane-wave calculations employing ultrasoft Vanderbilt pseudopotential and generalized gradient scheme. Simulated results display that band structure of polythiophene is very simple, which exhibits semiconductor performance with band gap of 0.9 eV. Neutral polythiophene is electrical insulator. Otherwise, by inspecting band structure, one can find that, except very close to zone center, where the density of states is high, band dispersions are considerable. The special band structure, thus, induces very low value of Seebeck coefficient (~20 μV/K), as shown in Figure 11.

Figure 11.

Calculated energy band structure (a), (b), and Seebeck coefficient (c), (d), for polythiophene (left) and polyaminosquarine (right), respectively [43]. Isolated planar polymer chains were used for this calculation.

Similar to calculated method from Gao et al., Shuai et al. have combined first-principles band structure calculations and Boltzmann transport theory to study thermoelectric in pentacene and rubrene crystals [46]. Electronic contribution to Seebeck coefficient is obtained in approximations of constant relaxation time and rigid band, as shown in Figure 12. Calculation results exhibited also the similar trend compared with experimental Seebeck coefficient.

Figure 12.

Seebeck coefficient calculated as function of charge carrier’s concentration (a) for pentacene and rubrene at room temperature and compared to FET measurements. Calculated Seebeck coefficients have been averaged over three crystal directions (b) for rubrene at temperature in the range between 200 and 300 K [46].

Afterward, Shuai et al. applied also combining method to calculate thermoelectric properties of organic materials, which is used to calculating α-form phthalocyanine crystals H2Pc, CuPc, NiPc and TiOPc [47]. This combining method includes first-principles band structure calculations, Boltzmann transport theory and deformation potential theory. They used first-principles calculations in VASP to calculate Seebeck coefficient. After obtaining band structure, Boltzmann transport theory was performed to calculate properties related to charge carrier transport, as in Eqs. (32) and (33).

Being different from Gao et al.’s and Shuai et al.’s previous works, it is assumed that relaxation time is a constant, which can be estimated by deformation potential theory for treating electron-phonon scattering. In terms of corresponding articles [46, 48], acoustic phonon scattering in both pristine and doped system was simulated by this theory including scattering matrix element for electrons from kstate to kstate expressing as:


where E1is deformation potential constant that represents energy band shift caused by crystal lattice deformation, and Ciiis elastic constant in the direction of lattice wave’s propagation. Then, relaxation time can be expressed by scattering probability:


where δ[ε(i,k)ε(i,k)]is Dirac delta function, and θis angle between kand k. In Bardeen and Shockley’s method, scattering is assumed to be isotropic, and matrix element of interactions M(k,k)is not relative with kand k.

Calculated dependences of Seebeck coefficient on charge carriers’ concentration are shown in Figure 13. Calculated dependences of Seebeck coefficient display different properties such as positive Svalues for holes and negative Svalues for electrons. Seebeck coefficient value is isotropic at first glance, and it decreases rapidly as charge carriers’ concentration increases.

Figure 13.

Seebeck coefficient for (a) H2Pc, (b) CuPc, (c) NiPc, (d) TiOPc calculated as a function of the charge carriers concentration at 298 K [47]. a, b and c denote a, b and c crystal axes, respectively.

4. Conclusion and outlook

In the past decades, the research on organic thermoelectric materials has made great progress. Rich variety of novel organic materials has been synthesized and applied in thermoelectric devices, and thermoelectric performances of organic semiconductors have been promoted greatly. However, as compared to inorganic thermoelectric materials, organic thermoelectric materials still exhibit lower ZT so far. However, situation looks like that progress in theoretical study of organic thermoelectric effect lags far behind experimental investigation in the last 30 years, but has been changed remarkably until the recent five years. Here, we have tried to describe organic thermoelectric materials and theoretical approaches, which allow to calculate characteristics of charge carrier transport processes responsible for thermoelectric effect in organic semiconductors. We hope that these contexts can be helpful to improve thermoelectric effect in organic materials and provide motivation for growth of thermoelectric applications of organic semiconductors. We believe that quest for green energy sources will stimulate intensive research and development works in the field of novel organic thermoelectric materials and devices that will result in serious improvement in thermoelectric efficiency of organic thermoelectric materials and enable development of novel high-performance and affordable organic thermoelectric devices.



density of states


variable-range hopping,










poly(2,5-dimethoxy phenylenevinylene),








ethylene glycol,


Seebeck coefficient,


Peltier coefficient,


thin film transistors,


power factor,


field-effect transistor,


mobility edge,


gate voltage,


threshold voltage,


Monte Carlo,


full-potential linearized augmented plane-wave,


Vienna ab initio simulation package.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Nianduan Lu, Ling Li and Ming Liu (December 21st 2016). Thermoelectric Effect and Application of Organic Semiconductors, Thermoelectrics for Power Generation, Sergey Skipidarov and Mikhail Nikitin, IntechOpen, DOI: 10.5772/65872. Available from:

Embed this chapter on your site Copy to clipboard

<iframe src="" />

Embed this code snippet in the HTML of your website to show this chapter

chapter statistics

1098total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Review of Research on the Thermoelectric Material ZnSb

By Xin Song and Terje G. Finstad

Related Book

First chapter

The Phenomenon of Wireless Energy Transfer: Experiments and Philosophy

By Héctor Vázquez-Leal, Agustín Gallardo-Del-Angel, Roberto Castañeda-Sheissa and Francisco Javier González-Martínez

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More about us