Thermoelectric Properties of Carbon Nanotubes

: Thermoelectric (TE) material is a class of materials that can convert heat to electrical energy directly in a solid-state-device without any moving parts and that is environmentally friendly. The study and development of TE materials have grown quickly in the past decade. However, their development goes slowly by the lack of cheap TE materials with high Seebeck coefﬁcient and good electrical conductivity. Carbon nanotubes (CNTs) are particularly attractive as TE materials because of at least three reasons: (1) CNTs possess various band gaps depending on their structure, (2) CNTs represent unique one-dimensional carbon materials which naturally satisﬁes the conditions of quantum conﬁnement effect to enhance the TE efﬁciency and (3) CNTs provide us with a platform for developing lightweight and ﬂexible TE devices due to their mechanical properties. The TE power factor is reported to reach 700–1000 µW/mK 2 for both p-type and n-type CNTs when puriﬁed to contain only doped semiconducting CNT species. Therefore, CNTs are promising for a variety of TE applications in which the heat source is unlimited, such as waste heat or solar heat although their ﬁgure of merit ZT is still modest (0.05 at 300 K). In this paper, we review in detail from the basic concept of TE ﬁeld to the fundamental TE properties of CNTs, as well as their applications. Furthermore, the strategies are discussed to improve the TE properties of CNTs. Finally, we give our perspectives on the tremendous potential of CNTs-based TE materials and composites.


Introduction
According to a report by Lawrence Livermore National Laboratory [1], it is said that only about 1/3 of the world primary energy consumption (oil, coal, gas, etc.) is effectively utilized, while most of the remainder is lost as waste heat.A solid-state thermoelectric (TE) device, which converts waste heat directly to electrical energy by utilizing the Seebeck effect, is expected as a solution for this problem [2].However, TE device has still not achieved widespread-practical applications [3] due to two main reasons: (1) TE performance of existing TE materials is not good enough and (2) commercial TE materials are based on the scarce and/or toxic elements.In order to solve these issues, many strategies have been proposed in recent years.
To understand the possible strategies for improving the TE performance, we need to control the basic parameters in thermoelectrics.The TE performance is determined by how much the electricity can be obtained for a given input heat.Power factor (PF) and figure-of-merit ZT, which depend on the Seebeck coefficient S, electrical conductivity σ and thermal conductivity κ for ZT of the materials, are usually considered for assessing the performance of the TE materials [2,4].If the source of heat is "free", such as waste-heat recovery from cars and solar heat, a maximum PF corresponds to the minimum cost of generating power.On the other hand, in the case that heat source is costly such as fossil fuel combustion, the cost of generating power is reduced by a large ZT value [5][6][7].
The fundamental problem of physics that has hindered the search for TE materials is a conventional trade-off relationship between S and σ [8].In addition, the situation in which electricity (heat) is well (poorly) conducted is paradoxical due to the Wiedemann-Franz law [8] for metal.In 1993, Hicks and Dresselhaus published two theoretical papers [9,10] to overcome the trade-off relationship between S and σ by using low-dimensional materials.In the one-dimensional (1D) or two-dimensional (2D) materials, the density of states (DOS) per unit volume is enhanced at bottom (or top) of the valence (or conduction) band, which would lead to increase the power factor, PF = Sσ 2 , over the bulk (3D) values.However, some previous experiments showed that the PF is not always enhanced with reducing the dimensional size of materials [11,12].Hung et al. [13,14] have proposed a theory suggesting that the discrepancy between theory and exceptions in experiments comes from a competition between the confinement length L (i.e., diameter or thickness of the 1D or 2D materials) and the thermal de Broglie wavelength Λ.They showed that the PF is only enhanced when L < Λ.This condition gives a challenge for experimentalists to obtain the enhancement PF of the 1D materials, while the L values of many 2D materials such as InSe [15], GaN [16] and SrTiO 3 [17,18] are early reduced to satisfy L < Λ.For 1D nanowires, their diameters can be reduced to several tens of nanometers, but most of them are still much larger than their thermal de Broglie wavelengths, which are about few nanometers [13].Carbon nanotubes (CNTs), on the other hand, with a wide range of diameters (0.4∼40 nm) [19], are reasonable candidates for 1D TE materials that can satisfy the condition L < Λ.
Other challenge in TE materials is to find a good TE material with low cost, nontoxic and natural abundance.The commercial TE materials are PbTe or Bi 2 Te 3 and their alloys which has a property of ZT ∼ 0.5-1 and PF ∼ 1500-4500 µW/mK 2 at room temperature [20,21].The scarcity of Te and the toxicity of Pb make newly emerging organic materials such as polymers as promising TE materials that do not contain these elements, despite their lower efficiencies.Due to their intrinsic low thermal conductivity (∼1.5 W/mK), a semiconducting polymer showed decent ZT values about 0.1-0.42 at room temperature, while the PF (∼78-500 µW/mK 2 ) is much smaller than that of Bi 2 Te 3 [22][23][24].Individual CNTs are not so promising as a TE material with high ZT values because of their intrinsic high thermal conductivity (∼2000-3000 W/mK) [25].However, by controlling of the doping level of the semiconducting CNT, their PF can reach 700-1000 µW/mK 2 [26,27].The higher PF of CNTs compared with the polymer makes that CNTs can apply for a variety of TE applications, in which the heat source is unlimited such as waste heat or solar heat.Furthermore, like the polymers, the CNTs can allow a low cost manufacturing method [28], which can reduce the total cost of the TE device.
In this review, we focus to explain fundamental TE properties of CNTs, in particular the semiconducting single-wall carbon nanotubes (s-SWNTs).Because many experimentalists are working on CNTs but not all are working with thermoelectrics, the backgrounds about the basic parameters of CNTs and the performance metric of CNTs-based TE materials and devices might be useful.Thus, we will discuss how diameter, electronic structure, doping level and temperature influence the TE properties of CNTs.After discussing the TE properties from a theoretical point of view, a series of recent experimental results is reviewed, where we show the potential to enable flexible, lightweight and low-cost TE devices.Finally, we give our perspectives on the the strategies to improve the TE properties of CNTs and the niche applications for CNT-based TEs as sensors, solar TE conversion, or bio-devices.

Basic Parameters for Thermoelectric Materials
In this section, we review important factors in thermoelectrics that are useful to evaluate the TE performance of CNTs.First, we examine a fundamental effect in thermoelectrics, the so-called Seebeck effect.Then, we show two general parameters for TE performance assessment, namely the power factor PF and the figure-of-merit ZT.The reason why we choose the 1D CNT as a TE material will also be discussed in the final part of this section.

Seebeck's Effect
In 1821, Thomas J. Seebeck, a German physicist made an important observation related to thermoelectrics.Connecting a bismuth wire to an antimony wire, he found a voltage generated between two free ends of a pair of wires by heating the junction region between the two wires.This was the first TE device and the phenomenon is known as Seebeck's effect, in which the Seebeck coefficient (or thermopower), S, is defined by how much the voltage, ∆V, can be obtained from an applied temperature difference, ∆T.The expression for the Seebeck coefficient is with the units V/K (volts per kelvin) in the SI unit.
One can understand the origin of the Seebeck effect from the illustration shown in Figure 1.Let us consider a material, in which one side is touched to a cold source while the other side to a hot source.There is a large distribution of electrons above the Fermi energy E F at the hot region according to Fermi-Dirac's distribution function (Figure 1 left), while there is small distribution of electrons above E F at the cold region (Figure 1 right).In addition, the carriers have the higher kinetic energy (KE) at the hot region because the average KE value is about 3  2 k B T for an ideal gas, where k B is the Boltzmann constant and T is the absolute temperature.From this fact, carriers in the hot region are more agile than those at the cold region.As a result, the carriers (either electrons above E F or holes below E F ) will diffuse from the hot region to the cold region, which may generate ∆V in the material.However, if the electrons and holes move in the same direction, we do not get a net electronic current because they will cancel charges each other.Thus a useful TE device should consists of p-type and n-type semiconductors.Since CNTs allow us to have both p-type and n-type semiconductors [27,29,30], a TE device with a single chirality CNT [31] used for both legs is feasible.

Power Factor PF and Figure of Merit ZT
Performance metrics of a TE device should be evaluated at maximum power and/or maximum efficiency.For maximizing the output power, it is required that one has to optimize the electrical power density Q in W/m 3 , which is defined as [6] where PF is the power factor that is defined by PF = S 2 σ in W/mK 2 , h in m is the leg length of TE devices, T c (in K) is the temperature of the cold region and T h (in K) is the temperature of the hot region.A high Q value requires a high PF value because the term (T h − T c ) 2 /4h is given by the boundary condition.
The maximum efficiency η of the TE device, on the other hand, is determined by the Carnot efficiency as follows [2]: where ZT is the dimensionless figure-of-merit of the material and is given by [2] where κ is the thermal conductivity and T = (T h + T c )/2 is the average temperature.Equation (3) shows that a larger ZT leads to a higher efficiency.To achieve the large ZT values, both low κ and high PF are required for a given T.An efficient way to reduce κ is using the low-dimensional structures, in which the κ is dominated by phonon heat transport.For example, two independent studies [11,12] reported that Si nanowires have the κ values below the theoretical limit of bulk Si (0.99 W/mK).Therefore, ZT > 1 at room temperature is obtained for a 1D silicon nanowire, which is much higher than bulk silicon (ZT = 0.01) [11,12].Another way to reduce κ is using the phonon glass system [32,33].Glassy materials have a lower κ compared with crystalline materials because phonons are scattered in the amorphous structure with mean-free paths in atomic dimensions.The phonon-glass characteristic also can be found in a perfectly crystalline solid such as semiconducting clathrates [34] or filled skutterudites [32].For example, the filled skutterudites with AgSbTe 2 showed a glassy thermal conductivity κ = 0.7 W/mK at 300 K with the mean-free path about 3.2-6.5 nm [35][36][37].Recently, Mukhopadhyay et al.
[38] showed that Ti 3 VSe 4 has ultralow κ = 0.3 W/mK at 300 K, in which the mean-free path is approaching the Ioffe-Regel limit (i.e., the average atomic distance in Ti 3 VSe 4 ).Because of the success in reducing κ, the next challenge of the TE is how to enhance the PF.

Low-Dimensional Thermoelectric Energy Conversion
A fundamental aspect in the PF is a trade-off relationship between S and σ [8].In early 1990s, Hicks and Dresselhaus [9,10] pointed out that the low-dimensional (1D & 2D) materials provide a way to solve the trade-off relationship between S and σ based on the quantum mechanics.In low-dimensional materials, the motion of electrons is limited in certain direction, in other words, electrons are confined in a potential in infinitely tall walls.In this case, it is possible to change the DOS of systems as shown in Figure 2a when the dimension of the systems, L, reduces to nanometre scale and this is the so-called "the confinement effect".The origin for the enhancement of the PF comes from the improved energy-dependent electrical conductivity σ(E) due to the enhanced DOS at the Fermi energy E F at the bottom (or top) of conduction (or valence) band (see Figure 2a).According to the Drude model, σ(E) can be expressed by σ(E) = n(E)q 2 τ/m, where q, m and τ are the unit charge, the effective mass and the relaxation time of carrier, respectively, and n(E) is the carrier density, which is given as n(E) = g(E) f (E), where f (E) is the Fermi-Dirac distribution function and g(E) is the DOS per unit volume and per unit energy.For 1D and 2D materials, g(E) is inversely proportional to L 2 and L, respectively.Therefore, the PF = S 2 σ is inversely proportional to L 2 and L for the 1D and 2D materials, respectively, as shown in Figure 2b.This prediction has made a lot of progresses on low-dimensional nanostructure and their applications in thermoelectrics [39][40][41][42][43].However, they also noticed that some common materials such as Si nanowires did not show enhancement of the PF although their confinement length L is sufficiently small (∼50 nm) [11,12].By re-evaluating the Hick-Dresselhaus theory, we show that the discrepancy between theory and experiments comes from an interplay between L and the thermal de Broglie wavelength Λ [13,14], in which Λ is defined by where h is the Planck constant.

3D 2D 1D
L In Figure 3a, we show the optimized power factor, PF opt , as a function of the ratio L/Λ for the 1D, 2D and 3D materials.The PF opt values are scaled by the PF opt of a 3D material (PF 3D opt ).The enhancement factor (PF opt /PF 3D opt ) is inversely proportional to (L/Λ) 2 and L/Λ for the 1D and 2D materials, respectively [13,14].By considering Λ, theoretical prediction is consistent with most of the experimental results as show by symbols in Figure 3a.Furthermore, recent experimental data of the 2D materials such as InSe [15] and GaN [18] showed that the PF is enhanced about several times compared with bulk case when the condition L < Λ is satisfied.However, this condition make a challenge for experimentalists to obtain the enhancement PF of the 1D materials since reducing the diameter of the 1D materials to sub-10 nanometer is still challenging even for the current technology.In this sense, carbon nanotube (CNTs) with a wide range of natural diameters (0.4∼40 nm) [19] show a potential to achieve the enhancement PF.The diameter of CNTs can be well controlled by using the metallic nanoparticles to grow CNTs [44,45].In Figure 3b, Λ of semiconducting single wall carbon nanotubes (s-SWNTs) is plotted as a function of diameter L for 0.5 < L < 1.5 nm [46].It shows that Λ is much larger than L. Thus, if s-SWNTs is well separated, single s-SWNTs with a small diameter should behave as a good TE material.

Thermoelectric Properties of SWNTs
In Tables 1 and 2, we list the observed TE properties of p-type and n-type SWNTs at room temperature including the Seebeck coefficient S, the electrical conductivity σ, the power factor PF, the thermal conductivity κ and the figure of merit ZT.In general, the TE properties of CNTs are considered for SWNTs.Before 2014, most of the SWNTs samples contain only metallic SWNTs (m-SWNTs) or the SWNTs with both m-and s-SWNTs, in the measured Seebeeck coefficients reported in the range of |S| = 14 − 60 µV/K for both p-type and n-type, respectively.These magnitude of the Seebeeck coefficients depend on the form of samples (bundles, individual or films), processing method, and dopants.In the case of m-SWNTs, the DOS suggests a low S since m-SWNTs have a nearly constant DOS around the Fermi energy, while semiconducting SWNTs (s-SWNTs) with a van Hove singularity near the Fermi energy can show a high S.These properties of m-SWNTs and s-SWNTs are observed by experiment for individual SWNTs [47], in which S values of m-SWNTs and s-SWNTs are 35 µV/K and 260 µV/K, respectively.Therefore, after 2014, most of the studies have focused on the TE properties of the purified s-SWNTs films.The TE power factors of s-SWNTs are reported over 1000 µW/mK 2 for both n-type and p-type at room temperature, which can compete with traditional telluride-based TE materials as BiTe 2 and PbTe.In this section, we will focus to explain the TE properties of s-SWNT.First, we will show a brief electronic properties of SWNTs, which depends on a specific (n, m) structure that determines the diameter of SWNTs.Then, the TE properties such as S, σ, PF and ZT are given as a function of the diameter of SWNTs.

Electronic Properties of SWNTs
One can imagine a SWNT as a single sheet of graphene which is rolled up into a cylinder [19].Graphene itself is a an atomic layer of carbon atoms arranged in a 2D honeycomb lattice.Since graphene is a building block of other graphitic materials in different dimensions, the electronic properties of SWNTs can be derived from those of graphene.First we start with the geometrical structure of the SWNT according to the graphene geometry.The SWNT circumference can be calculated from the chiral vector C h , which is written in terms of graphene's lattice vectors a 1 and a 2 [19], where (n, m) is an integer pair requiring n ≥ m, n > 0 and m ≥ 0. The formulas for the circumferential length |C h | and diameter L are given by: where  From top to bottom, the chiral index is (5, 0), (4, 2), (3,3).Thick solid lines in (a) and (c) are to stress the definition of "zigzag" and "armchair" structures, respectively.Redrawn with different style by the authors adapting from Reference [62].
The 1D unit cell of a SWNT is expressed by C h and T (a translational vector in the direction of tube axis [19]), while the corresponding vectors in reciprocal space are K 1 along the tube circumference and K 2 along the tube axis.Note that only K 2 is a 1D reciprocal lattice vector, while K 1 gives discrete wave vector values in the direction of C h .The formulas for K 1 and K 2 can be obtained from their relations with C h and T: It follows, where b 1 and b 2 are graphene's reciprocal lattice vectors.The vectors K 1 and K 2 are shown in Figure 5 for the (4, 2) SWNT.The 1D Brillouin zone (BZ) of the SWNT is then constructed by the N hex line segments with length of K 2 .These lines are often referred to as the "cutting lines" [63].Due to the presence of cutting lines, the wave vectors k allowed in a SWNT are given by where µ = 0, 1, . . ., N hex − 1 is the index of each cutting line and the 1D k value is within Before showing the SWNT dispersion relation, we first show the electronic structure of graphene in Figure 6a.We can see that graphene exhibits linear energy dispersions around the K and K ′ points in the 2D BZ of graphene [Figure 6b].In Figure 6c, we show the contour of graphene's energy dispersion near the K point.We can see that the energy surface changes from circle to triangle by increasing distance from the K point.This behavior is known as the trigonal warping effect [65], an important feature in the electronic structure of SWNTs.In SWNTs, the allowed wave vectors k around the SWNT circumference is quantized.Therefore, the SWNT energy dispersions can be obtained from the corresponding graphene's energy dispersions along the cutting lines.When we superimpose the 1D cutting lines of a SWNT, that is, in Equation (11), on the 2D energy dispersion surface of graphene, we can obtain N hex pairs of energy dispersions the SWNT: where E b 2D denotes 2D energy dispersion of the π bonds of graphene.For a specific (n, m) SWNT, if a cutting line passes through the K or K ′ point of the 2D BZ of graphene, where the valence and conduction bands touch each other, the 1D energy bands of the SWNT will have a zero gap, therefore, it is classified as a metallic SWNT.On the other hand, if the cutting line does not pass through the K or K ′ point, the SWNT is semiconducting with a finite energy gap. Figure 7 shows two examples of the SWNT dispersion relations.In Figure 7a the armchair (6,6) SWNT is metallic, while in Figure 7b, the zigzag (10, 0) SWNT is metallic, Beside their structure, we can further classify SWNTs based on their metallicity.To obtain this classification, let us first see Figure 8a.If we project the ΓK vector pointing toward the K point onto the K 1 direction perpendicular to the cutting lines, that can be denoted by Note that, if (2n + m)/3 is an integer, ΓK has an integer number of K 1 components, thus one of the cutting lines will pass through the K point, giving rise to a metallic SWNT.On the other hand, if (2n + m)/3 is not an integer, i.e, the remainder of the division is either 1 or 2, the K point will lie at one third or two third of the spacing between two adjacent cutting lines near the K point, thus giving rise to a semiconducting SWNT, as shown in Figure 8b.Therefore, three types of SWNTs are obtained based on their metallicity.They are often referred to as M-, S1-and S2-SWNTs, which are given by the following conditions: respectively.The S1-and S2-SWNTs are also often written as type-I and type-II semiconducting SWNTs.

Seebeck Coefficient of SWNTs
The Seebeck coefficient S of SWNTs are obtained by solving the linearized Boltzmann equations as [66] where T is the average absolute temperature and L i (i = 0, 1) is the transport integral that is defined by where E is the energy of electron.In Equation ( 18), v x (E) is the group velocity of electron in the x direction, τ(E) is the relaxation time of electron and g(E) is the DOS.
is the Fermi-Dirac distribution function, where µ is the chemical potential, which is referred to the Fermi energy.From Equations ( 17) and ( 18), S depend on three parameters including the energy band structure E(k), the chemical potential µ and the temperature T, in which E(k) depend on a specific (n, m) structure (or diameter) of SWNTs as shown in Section 3.1.
In Figure 9a, we show the calculated results of S for (11, 0) s-SWNT, as a function of µ and T. The Seebeck coefficient increases when temperature decreases because S ∝ 1/T as shown in Equation (17).The maximum S of (11, 0) s-SWNT is about 1420 µV/K at |µ| = ±0.05eV and T = 300 K, which is much larger than that of bundled SWNTs (S∼100-200 µV/K [48,67]).In Figure 9b, we plot S at T = 300 K as a function of µ for the (11, 0), (12,4) and (15,5) SWNTs.For all chiralities, the maximum |S| values at the negative and positive sign of µ arise due to the p-type and n-type s-SWNTs, respectively, which is consistent with an experimental observation by Yanagi et al. [67].Since the Seebeck coefficient depends on the chemical potential, it is possible to control the Seebeck coefficient of s-SWNTs by applying a gate voltage [47].In order to analyze the numerical results of the Seebeck coefficient (see Equation ( 17) and Figure 9), we have derived an analytical formula for the Seebeck coefficient by a two-band model [4,68].It is noted that the contribution of two bands is required to get the optimum of S for s-SWNTs since S is a linear function with a large µ value within the one-band model [13],.By using the two-band model, S of s-SWNTs is given as [31] where E g is the band gap of s-SWNT and r is a characteristic exponent of the scattering mechanism [8].By using the constant-relaxation-time approximation (CRTA) [69] with r = 0, the Seebeck coefficient of s-SWNT is given by In Figure 9b, the solid and dashed lines represent the numerical and analytical results of the Seebeck coefficient using Equation (17) and Equation (20), respectively, for the (11, 0), (12,4) and (15, 5) s-SWNTs.As shown in Figure 9b, the analytical results from Equation (20) reproduce the numerical results near µ = 0, in which S for the p-type and n-type s-SWNTs can be well-reproduced at near µ = 0 (i.e., the Fermi energy is located at the middle of band gap).At larger |µ| values, the two-band model is no longer valid for S at a higher doping level.Therefore, the analytical results deviate from the numerical results at larger |µ| values, as shown in Figure 9b.However, if we focus on the optimum values of S, the two-band model is sufficient to calculate S of the s-SWNTs.
In Figure 10a,b, we plot, respectively, the optimum S values and the band gaps of many s-SWNTs as the functions of diameter L for 0.5 < L < 1.5 nm.In Figure 10a, S increases when the diameter L decreases.For the s-SWNTs with L < 0.6 nm, such as s-SWNTs with 2n + m = 13 (or the (5, 3) and (6, 1) s-SWNTs), S can reach the values more than 2000 µV/K [31].These values are about 6 − 10 times larger than the Seebeck coefficient of the common thermoelectric materials [11,12,20,70].By using the density function theory, Avery et al. [26] showed that S = 1285 µV/K for (7, 5) s-SWNT and S = 705 µV/K for (16, 0) s-SWNT, in which the diameter of (7, 5) s-SWNT (∼ 0.82 nm) is smaller than that of (16, 0) s-SWNT (∼ 1.26 nm).The larger S of s-SWNTs with the smaller diameter can be explained by the relation of S CNT to E g as shown in Equation (20) and by the fact that E g ∝ 1/L [65].The 1D character of DOS of a SWNT enhances the Seebeck coefficient [9,10].Here, it is pointed out that S of s-SWNTs depends on not only the diameter but also the family pattern, in which mod(2n + m, 3) = 1 and mod(2n + m, 3) = 2 represent the nanotube SI and SII families, respectively [65].The family pattern dependence of S is very similar to that found in Figure 10b, which is often referred to as the Kataura plot [65,71,72].Since the band gap is directly connected the Seebeck coefficient as can be seen in the S CNT formula of Equation (20), it suggests that the measurement of the Seebeck coefficient of a single s-SWNT sample might be able to predict their band gap value.(20)] are denoted by circles and plus symbols, respectively.(b) Band gap of many s-SWNTs are plotted as a function of SWNT diameter that shows the family pattern of s-SWNTs, in which the solid lines connect the s-SWNTs with the same family number 2n + m.The SI and SII tubes correspond to the s-SWNTs with mod(2n + m, 3) = 1 and 2, respectively.Adapted with permission from Reference [31].Copyright 2015, American Physical Society.

Thermoelectric Power Factor of SWNTs
The power factor is defined by PF = S 2 σ (see Section 2.2), in which S is given by Equation ( 17) and the electrical conductivity is expressed by σ = L 0 .The problem to optimize the PF is a trade-off between maximizing S and maximizing σ.This because a maximum value of S in a semiconductor usually corresponds to a minimum values of σ.The relationship between S and σ of a semiconductor could be qualitatively understood using the modified Mott-Keike formula: [73] where β and σ max are fitting parameters.Note that σ max can be derived from Mott-Keike formula [74] as follows: where c ≪ 1 is the fraction of carriers available in the energy bands.σ max might be interpreted as an possible maximum conductivity at S = 0.However, σ max = σ is inherently not achievable from Equation (22).
In Figure 11a,b, we show the S − σ curves for several β and σ max values, respectively.By increasing β and σ max values, S − σ curve decreases and increases, respectively.Therefore, the PF decreases by increasing β as shown in Figure 11c, while the PF increases by increasing σ max as shown in Figure 11d.Equation ( 21) thus tells us that a smaller β value and higher σ max value will give a larger PF value.For the SWNTs, the values of β and σ max can be controlled by the preparation of samples or the dopants.In Figure 12a, the S − σ curves are plotted for thin films of high-pressure disproportionation of carbon monoxide (HiPCO) s-SWNTs extracted with supramolecular polymer (SMP) or poly[(9,9-dihexylfluorenyl-2,7-diyl)-co-(9,10-anthracene)] (PFH-A) [27], in which s-SWNTs are treated with increasing levels of triethyloxonium hexachloroantimonate (OA) as a p-type semiconductor [75].The solid and dashed trend lines are plotted using Equation (21).The maximum achievable of σ of the s-SWNTs films with the SMP is about 4 × 10 5 1/Ωm, which is much larger than that of the s-SWNTs films with the PFH-A (1.5 × 10 5 1/Ωm).Since the S − σ curve of s-SWNT films from the SMP shifts to the right compared with that from the PFH-A, the σ max value of s-SWNT films from the SMP is higher than that from the PFH-A.Therefore, the maximum PF for the s-SWNTs film from SMP (560 µW/mK 2 ) is improved relative to the film from PFH-A (340 µW/mK 2 ), as shown in Figure 12b.In Figure 12c,d, we show a similar OA doping series to Figure 12a,b, but comparing two plasma-torch (PT) s-SWNTs films dispersed with the SMP and poly[(9,9-di-n-dodecyl-2,7-fluorendiyl-dimethine)-(1,4-phenylene-dinitrilomethine)] (PFPD) [27].In Figure 12c, the maximum σ for the PFPD-based film reaches 6 × 10 5 1/Ωm, doubling that of the SMP-based film (∼ 3 × 10 5 1/Ωm).Therefore, the PF increases from 300 µW/mK 2 in the case of SMP film to 705 µW/mK 2 in the case of PFPD film.In Figure 12e, we show the relationship between the maximum PF and the maximum σ that are obtained at the optimal and maximum doping, respectively.Since, the maximum PF a linear function of the maximum σ, it suggests that the maximum σ can be used as a quickly metric for screening the TE potential of s-SWNT samples [27].It notes that in order to obtain the maximum PF, both σ and S for a wide range of chemical potential are needed that requires a significant time and effort.In contrast, the maximum σ can be easily obtained by a heavy doping, that can be achieved by long exposure (> 10 min) of the s-SWNT to a concentrated dopant solution [27].
In Figure 13a, the maximum PF value is plotted as a function of the diameter of the s-SWNTs films, in which these films are extracted by the cleavable SMP and non-removable polymers (PFO-BPy or PFH-A) [27].Both cleavable and non-cleavable s-SWNTs films show a similar diameter dependence of the PF, and it suggests that the maximum PF increases s-SWNT diameter decreases.Since s-SWNT diameter is an important parameter for optimizing PF, the next question is which specific chiral-vector (n, m) will give a highest PF.In order to analyze the dependence of s-SWNTs diameter and chiral vector (n, m) on the maximum PF, we have derived an analytical formula for the PF by a one-band model as [13,14] where B(x, y) is the Beta function, d = 1, 2 or 3 denotes the dimension of the material and µ 0 is the carrier mobility is defined by [13] where Γ(x) is the Gamma function.For 1D s-SWNTs (d = 1), r and τ 0 are set to 0 and 1 ps, respectively, in the CRTA case for the sake of simplicity.In Figure 13b, we show the optimum PF values of many s-SWNTs as a function of diameter L for 0.5 < L < 1.5 nm.The PF is found to increase as the diameter L decreases.Diameter-dependent trends in the PF is consistent with the observed experiment results in Figure 13a.More importantly, it is pointed out that the PF of the s-SWNTs shows the family pattern, which is very similar to that found in S as a function of diameter in Figure 10a.For the s-SWNTs with L < 0.6 nm and 2n + m = 14, that is, the (6, 2) and (7, 0) s-SWNTs, the PF can reach a value more than 1100 µW/mK 2 for both n-type and p-type s-SWNTS.We note that the s-SWNTs with 2n + m = 13 show the highest S values, while the s-SWNTs with 2n + m = 14 show the highest PF values because of the correlation between S and σ of the s-SWNTs.Therefore, not only the nanotube diameter but also the selection of the chiralty (n, m) may be an important rule for optimizing the power factor of the s-SWNTs.

Thermal Conductivity and Figure of Merits of SWNTs
The figure of merit ZT is determined by both the PF and the total thermal conductivity as ZT = PF × T/(κ e + κ ph ), where κ e and κ ph are the electronic and phonon contribution to total thermal conductivity, respectively.The Wiedemann-Franz law is commonly used to calculate the magnitude of κ e at temperature T, which is defined by [76] where known as the Lorenz number.In Equation (25), κ e is proportional to σ.Since σ ∼ 30 − 3000 1/Ωcm of the s-SWNTs (see Tables 1 and 2) is much smaller than that of copper (5.96 × 10 5 1/Ωcm [77]), κ e of the SWNTs might be a small value.Hone et al. [78] showed that the total thermal conductivity of the SWNTs are orders of magnitude higher than expected from Equation (25), implying that the thermal transport in the SWNT is dominated by phonons.Zhao et al. [79] also showed that κ e is measured to be less than 1.1 W/mK for the SWNT films.
The thermal conductivity of a phonon mode in the SWNTs is given by κ ph ∝ Cvl [80], where C is the heat capacity, v is the group velocity of the phonon and l is the mean free path (MFP) of the phonon.Pop et al. [81] showed that the thermal conductivity is nearly 3500 W/mK at room temperature (T = 300 K) for the SWNT with length about 2.6 µm and diameter about 1.7 nm.This value is close to the thermal conductivity of graphene about 4000 W/mK at T = 300 K [82].κ ph is large in the SWNTs and graphene for the following two reasons: (1) the carbon atoms in the SWNTs have a small mass and the sp 2 covalent bonds have large force constants of their vibration, which leads to large value of v; and (2) the structure has high value of MFP since many ahharmonic terms are suppressed or give small values.
In Figure 14a, we show the thermal conductivity of graphene as a function of temperature T. Since κ ph is proportional to 1/T at near the room temperature, κ ph may be limited by the three-phonon processes [83].In graphene, the three-phonon processes must follow the selection rule.Therefore, all three-phonon processes having an odd number (1 or 3) of out-of-plane phonon modes (ZA and ZO) are excluded [84].In particular, the three-phonon processes such as ZA + LA ↔ LA, ZA + TA ↔ LA, ZA + TA ↔ TA, ZA + ZA ↔ ZA, ZA + ZO ↔ ZO and ZA + ZA ↔ ZO do not occur.The selection rule show that about 60% of the three-phonon scattering of the ZA phonon is forbidden.Therefore, the dominant contribution ∼ 75% to κ ph comes from the ZA phonon, which is higher than other acoustic phonon (TA and LA) [84,85].In Figure 14b, the thermal conductivity of SWNTs is plotted as a function of nanotube diameter at T = 300 K.For very small diameter of the SWNTs, the optic phonons has sufficiently high frequencies where the optic phonon modes contribute little to the scattering processes with acoustic phonons.By increasing the diameter, the optic phonons shift to lower frequencies and they significantly scatter acoustic phonons, which leads to reduce κ ph of the SWNTs.The thermal conductivity is found to reach a minimum for the (20, 0) SWNT with the diameter about 1.6 nm, as shown in Figure 14b.As the diameters continue to increase, the optic phonons not only scatter the acoustic phonons, but also begin to contribute to κ ph due to higher phonon populations.In addition, the graphene selection rule begins to play a role for ZA phonon contribution in large diameter.In the limit that diameter → ∞, the thermal transport properties of the SWNTs should approach that of graphene.
The high thermal conductivity leads to a small ZT values of the SWNTs in the range of 10 −3 to 10 −4 [87].In order to enhance the ZT value, we need to reduce the thermal conductivity of the SWNTs, while the electronic properties are still maintained to obtain the good power factor.Since the commercial thermoelectric materials have a thermal conductivity about 1-10 W/mK, a significant reduction in the thermal conductivity of the SWNTs is to be done and it is still a challenge in the TE field.In the next section, we will show some strategies to improve the TE properties of the SWNTs.

Optimizing Power Factor by Controlling Doping
As discussed in Section 3.3, the thermoelectric power factor depends on S and σ, in which both S and σ are functions of the chemical potential, temperature and energy dispersion.Since the energy dispersion depends on the chirality (n, m) (or diameter) of the SWNTs and the temperature depends on the specific application, changing the chemical potential is an easy way to optimize the PF.For the SWNTs, the molecular dopants can be used to tune the chemical potential by ground state charge transfer interactions that inject holes (p-type) or electrons (n-type) into the SWNTs [75,88].Furthermore, the SWNTs have some advantages for using charge transfer dopants.First, the SWNTs with their 1D structure have very high surface areas, where the adsorption molecules can easy access on the tube surface.Second, the SWNTs have the strongly of the covalent bonds between C atoms in-plane and the sensitivity of the π bonds out-of-plane, so the SWNTs can interact with the molecular dopants without modification of structure or electronic structure.Therefore, the SWNT can be doped simply be soaking in solutions of the molecular dopants, which is in marked contrast to most conventional TE materials that often create defects by atom substitution [21].
In Figure 15a, we show many molecular dopants, which are available for p-type or n-type doping for the s-SWNT films.In general, Lewis acids are used for the p-type SWNTs, while Lewis bases are used for the n-type SWNT [53,75,88].In Figure 15b, the power factor are obtained for many p-type and n-type molecular dopants in Figure 15a.The power factor can reach about 26 µW/mK 2 for both p-type and n-type SWNTs with TCNQ and dpmp dopants, respectively.We note that the pristine SWNT film has the positive value of S = 49 µV/K, indicating that pristine SWNT exhibit a p-type characteristic [89].
The power factor can be improved by selecting a specific molecular dopant, but it might not be the optimized value.In order to optimize the power factor, a wide range of chemical doping is needed to find the optimum chemical potential.Baker et al. [26] used a doping strategy with the triethyloxonium hexachloroantimonate (OA) for the s-SWNT films.In this case, the chemical potential (or the Fermi energy) is controlled by the varying concentrations of the OA, in which the Fermi-energy shift is measured by the bleaching of the first S 11 and second S 22 exciton in the normalized absorption spectrum, as shown in Figure 16a,b.This strategy shows a large value of the Fermi-energy shift, resulting in optimization value of the power factor up to 705 µW/mK 2 in the s-SWNT films, as shown in Figure 12d.
Another strategy toward to control the Fermi energy is using an electric field.In particular, the ionic gating has been used to control of the Fermi energy over a wide range of carrier densities [67].As shown in Figure 17a, the gate voltage of a electric-double-layer transistors is used to tune the carrier density within the s-SWNT films [67] for both n-type and p-type cases.This device has also been applied for a single s-SWNT sample [47].In Figure 17b,c, the Seebeck coefficient and power factor are showed as a function of the applied gate voltage.Shimizu et al. [90] showed that the optimum power factor can be obtained when the Fermi energy is located within the second sub-band of the s-SWNTs.The chemical doping also show that the optimum power factor is observed at heavy doping (or high carrier concentration), as shown in Figure 16.Therefore, both strategies suggest that control of the Fermi energy around the second sub-band of the s-SWNTs is important to optimize the power factor.

Reducing Thermal Conductivity by Using SWNT Networks
As discussed in Section 3.3, an individual s-SWNT has high thermal conductivity (3500 W/mK) due to the very long MFP (∼0.5-1.5 µm [78]) of the phonons.In contrast, the s-SWNT film with a random network, as shown in Figure 18a, has very low thermal conductivity about 1-5 W/mK [27] compared with the individual s-SWNT.We note that the thermal conductivity of the s-SWNT film with aligned sample (∼ 250 W/mK [91]) is higher than that with a random network.Therefore, the crossed SWNT junctions of the network play a important role for low thermal conductivity.Chalopin et al. [92] derived the thermal conductivity of a s-SWNT network as where h is the length of the SWNT, d is the average distance between junctions and G is tube-tube contact conductance that is estimated about 5 pW/K [93].Since d can not be smaller than the nanotube diameter, Equation (26) shows an upper bound by a given h.By considering a typical nanotube length of 1 µm, Chalopin et al. [92] showed that all κ of many SWNT networks for different nanotube chiralities are of the same order and smaller than ∼ 5 W/mK.This is also consistent with the experimental data for various s-SWNT networks, as shown in Figure 18b.The very low upper bound for κ of the SWNT networks is abour 700 times smaller than that of the individual SWNT, which leads to ZT ∼ 0.12 in the SWNT networks [27].
An important fact from Equation ( 26) is the relationship between κ of the SWNT network and the diameter of the individual SWNTs making up the network.To see this, it is enough to realize that G is roughly proportional to the contact area.Therefore, G is proportional to the square of the diameter [92] of the SWNTs, which leads to reduce κ by decreasing the diameter.MacLeod et al. [27] found that κ = 1.3 W/mK for a s-SWNT network with average diameter of 1.1 nm, while κ = 2.9 W/mK for a s-SWNT network with average diameter of 1.3 nm, which is consistent with the above story.This behavior may look the opposite to the behavior of the thermal conductivity as discussed in Figure 14b for the individual s-SWNTs with the diameter < 1.6 nm.However, in the case of Figure 14b, they theoretically obtain the thermal conductivity as a function of diameter, while the thermal conductivity is discussed as a function of the contact area between two SWNTs in the case of the SWNTs networks.Therefore, for the s-SWNTs networks, the phonon scattering at the contact area might be more important than the intrinsic phonon scattering in the each individual s-SWNT.

Figure of Merit of Bi 2 Te 3 -SWNT Hybrid
Although the SWNT networks with low thermal conductivity (1 − 5 W/mK) can show a promising value of ZT about 0.12 [27], it is still far for industrial application.On the other hand, the traditional thermoelectric materials such as Bi 2 Te 3 show a best performance at near room temperature with ZT ∼ 0.8 [94].However, intrinsic brittleness and hardness of Bi 2 Te 3 limit their use as flexible TE device.Recently, Jin et al. [95] show that sputtering bismuth telluride onto the SWNTs makes a high-performing and flexible TE composite.By using the transmission electron microscopy (TEM) analyses, they showed that the Bi 2 Te 3 grows on the SWNT surfaces with the [1210] direction of Bi 2 Te 3 aligning with the length of the SWNTs.Interestingly, the Bi 2 Te 3 -SWNT hybrid material shows no significant changes in the Seebeck coefficient and the electrical conductivity due to the presence of the SWNT framework, while the thermal conductivity is about half that of bulk Bi 2 Te 3 due to more phonon scattering effect.Therefore, this hybrid material has the PF of ∼ 1600 µW/mK 2 and ZT of 0.9 at room temperature [95], which is comparable to that of a commercial bulk Bi 2 Te 3 , and is much higher than other flexible TE materials [96].Furthermore, the Bi 2 Te 3 -SWNT hybrid material can maintain the TE performance after 100 bending cycles with a bending radius about 4 mm and no obvious fractures or microcracks [95].

Summary and Perspective
In summary, many experiments and theories have demonstrated the high TE power factor for both p-type and n-type carbon nanotubes.When the sample contains only doped semiconducting SWNT species, the power factor is reported to reach 700-1000 µW/mK 2 .The higher PF of the SWNTs compared with the polymer allows the SWNTs to be applied for a variety of TE applications, in which the heat source is unlimited such as solar heat.The power factor can also be further improved by using a small diameter and selected chirality of the s-SWNT.At present, most of the TE devices is using the s-SWNT films with a 2D or 3D structure.Hopefully, the TE devices based on a single 1D s-SWNT can be made in the near future.
Furthermore, like the polymers, the SWNT shows a low cost manufacturing method that can replicate hundreds of thermocouples in a TE device.A simple comparison with solar cells shows the importance of material costs.Currently, a TE device has the efficiency five times lower than solar cells, but it requires fifty times more active material.Therefore, enhancement the power factor is an important step towards the widespread use of the TE.
Another point to overcome is that the SWNTs show a small ZT value about 0.12 at T = 300 K compared with the traditional TE materials such as Bi 2 Te 3 (ZT = 0.8 at T = 300 K).At present, this problem has not yet been solved, mainly due to the high thermal conductivity of the SWNTs.However, one can foresee many applications to take advantages of both the SWNTs and the traditional TE materials.For example, Bi 2 Te 3 -SWNT hybrid material showed a high-performing and highly flexible TE composite with the power factor and ZT of about 1600 µW/mK 2 and 0.9, respectively [95].A flexible TE device can be used for many applications such as clothes capable of powering various devices from body heat.

Figure 1 .
Figure 1.Illustration of the Seebeck effect.By moving electron from left to right, voltage appears.

Figure 2 .
Figure 2. (a) Electronic density of states for a bulk (3D), two-dimensional (2D) and one-dimensional (1D) materials.(b) Enhancement power factor PF is potted as a function of confinement length L for low-dimensional materials.

Figure 3 .
Figure 3. (a) Enhancement of power factor PF scaled by that for 3D materials is plotted as a function of L/Λ, where L and Λ are the confinement length and the thermal de Broglie length, respectively.The lines and dots, respectively, denote theory and experiments for several materials.It is clear that enhancement PF occurs when L < Λ (or L/Λ < 1) for the 1D and 2D materials.Adapted with permission from Reference [13].Copyright 2016, American Physical Society.(b) Thermal de Broglie wavelength Λ of many (n, m) s-SWNTs is plotted as a function of s-SWNT diameter L. Adapted with permission from Reference [46].Copyright 2019, Springer Publishing.
a ≡ |a 1 | = |a 2 |.The chiral angle θ is defined as the angle between C h and a 1 with 0 ≤ |θ| ≤ 30 • .One can classify the SWNT according to its (n, m) value that follows a certain symmetry as shown in Figure 4, from which three types of SWNTs emerge: (a) zigzag, (b) chiral and (c) armchair SWNTs.Various SWNT geometries will give rich physical properties of SWNTs.n , m; n≠m Chiral: =30 °;  n , n ; n=m Armchair:

Figure 5 .
Figure 5. Vectors K 1 and K 2 in reciprocal space and the BZ of a (4, 2) SWNT, which is represented by N hex = 28 parallel cutting lines.Reproduced from the Ph.D. thesis of one of the authors [64].

Figure 6 .
Figure 6.Graphene's energy dispersion in its first BZ.The energy dispersion is shown for the whole region of the 2D BZ in panel (a), while in panel (b) the energy dispersion is shown near the K point.Panel (c) gives the contour plot of the dispersion around the K point.Reproduced from the Ph.D. thesis of one of the authors [64].

Figure 7 .
Figure 7. Examples of 1D dispersion relations for SWNTs: (a) armchair SWNT with chirality(6,6), and (b) zigzag with chirality (10, 0).There are no bandgaps can in the case (a), thus the SWNT is metallic, while the SWNT in (b) is semiconducting due to the appearance of a gap.Redrawn with different style by the authors adapting from Reference[62].

Figure 8 .
Figure 8.(a) Condition for metallic energy bands is related to the ratio of the length of vector YK to that of K 1 .If the ratio is an integer, metallic energy bands are obtained [19].(b) Three possible configurations of the cutting lines in the vicinity of the K point depending on the value of mod(2n + m, 3).From left to right, the nanotube type is M-(metallic), S1-(type-I semiconducting) and S2-(type-II semiconducting) SWNT, respectively.The solid lines represent the cutting lines and the dashed lines indicate the KM directions, which are the boundaries of the first BZ of the SWNT.Redrawn with different style by the authors adapting from Reference [62].

Figure 10 .
Figure 10.(a) Optimum thermopower (or Seebeck coefficient) values of many s-SWNTs are plotted as a function of SWNT diameter at T = 300 K. Numerical and analytical results [Equation(20)] are denoted by circles and plus symbols, respectively.(b) Band gap of many s-SWNTs are plotted as a function of SWNT diameter that shows the family pattern of s-SWNTs, in which the solid lines connect the s-SWNTs with the same family number 2n + m.The SI and SII tubes correspond to the s-SWNTs with mod(2n + m, 3) = 1 and 2, respectively.Adapted with permission from Reference[31].Copyright 2015, American Physical Society.

5 × 10 5 max = 5 × 10 5 Figure 11 .
Figure 11.Effect of the fitting parameters (a) β and (b) σ max on the Seebeck coefficient S -electrical conductivity σ curves described by the modified Mott-Heike relationship in Equation (21).Effect of the fitting parameters (c) β and (d) σ max on the power factor PF -electrical conductivity σ curves.

Figure 12 .
Figure 12.(a) Seebeck coefficient (or thermopower) and (b) power factor are plotted as functions of electrical conductivity for high-pressure disproportionation of carbon monoxide (HiPCO) s-SWNT networks.(c) Thermopower and (d) power factor are plotted as functions of electrical conductivity for plasma-torch (PT) s-SWNT networks.(e) Relationship between the maximum power factor and the maximum conductivity for several s-SWNT networks.Adapted with permission from Reference [27], Copyright 2017, The Royal Society of Chemistry.

Figure 13 .
Figure 13.(a) Power factor is plotted as a function of diameter for several s-SWNTs films.Adapted with permission from Reference [27].Copyright 2017, The Royal Society of Chemistry.(b) Optimum of the power factor is plotted as a function of diameter for many s-SWNTs at room temperature (T = 300 K), in which the solid lines connect the s-SWNTs with the same family number 2n + m.

Figure 14 .
Figure14.(a) Thermal conductivity is plotted as a function of temperature for graphene with several concentrations of13 C. The lines and dots, respectively, denote theory and experiments.Adapted with permission from Reference[85].Copyright 2018, American Physical Society.(b) Thermal conductivity of the SWNTs is plotted as a function of the nanotube diameter.Blue squares, yellow circles and green triangles denote the armchair, zigzag and chiral SWNTs, respectively.Length of the SWNTs is set to 3 µm and temperature T is set to 300 K for all cases.All values are scaled by thermal conductivity of graphene κ graphene .Adapted with permission from Reference[86].Copyright 2010, American Physical Society.

Figure 15 .
Figure 15.(a) Chemical structures of p-type and n-type molecular dopants.(b) Power factors of pristine and doped SWNT films at 310 K. Adapted with permission from Reference [53].Copyright 2013, Nature Publishing Group.

Figure 16 .
Figure 16.(a) Normalized absorbance spectra is plotted as a function of wavelength for the laser vaporization (LV) SWNT films.(b) Electrical conductivity is plotted as a function of the doping-induced bleach of the S 11 and S 22 peaks enveloped for the s-SWNT films.Adapted with permission from Reference [57].Copyright 2016, American Chemical Society.

Figure 17 .Figure 18 .
Figure 17.(a) Schematic structure of ion-gated transistor devices for the n-type and p-type SWNTs films.(b) Seebeck coefficient and (c) power factor are plotted as functions of channel voltage, in which the dotted lines are inserted to identify the regions of the p-type and n-type SWNTs films.Adapted with permission from Reference [67].Copyright 2014, American Chemical Society.(a)(b)

Table 1 .
Thermoelectric properties of p-type single wall carbon nanotubes (SWNTs) at room temperature including the Seebeck coefficient S, the electrical conductivity σ, the power factor PF, the thermal conductivity κ and the figure of merit ZT.

Table 2 .
Thermoelectric properties of n-type single wall carbon nanotubes (SWNTs) at room temperature including the Seebeck coefficient S, the electrical conductivity σ, the power factor PF, the thermal conductivity κ and the figure of merit ZT.