Optimizing thermocouple’s ZT through design innovation | Scientific Reports

Scientific Reports volume 11, Article number: 19338 (2021) Cite this article

1820 Accesses

3 Citations

Metrics details

This work demonstrates that in parallel with the one existed at high doping concentration, there also exists an optimal combination of the transport properties of a thermoelectric material at low doping concentration as the curve of the relation between electrical conductivity and doping concentration is rigidly shifted toward that direction without disturbing the Seebeck coefficient and the thermal conductivity. Based on this finding, a new thermocouple design that uses low doping legs and high doping semiconductors as the external carrier injectors surrounding the legs is developed. The analytical model developed for the new thermocouple indicated that its efficiency and power output could be more than tripled as compared to those of the original design. A single thermocouple made of Silicon semiconductors was simulated numerically using different sets of input parameters. The results showed that the density of the externally injected carriers played a significant role in enhancing the thermocouple’s efficiency and power output.

Enhancing energy-use efficiency through thermoelectric generator, a solid-state device that can reliably convert heat released from vehicle exhaust and various manufacturing processes directly into electricity without the need for any moving mechanical components, constitutes one of the important measures toward building a sustainable environment. Such a technique can not only reduce fuel consumption that is depleting the limited oil reserves but also lower greenhouse gas emissions that have a tremendous detrimental impact on our environment. More than half of the energy supplied to the current thermal-process-based power generation systems or combustion engines is lost mostly as heat1,2. The potential to convert the released heat into electricity through thermoelectic generator is significant. Figure 1 shows the architecture of a flat planar thermoelectric power generator used in most of the contemporary designs. Evidently, the power output and the efficiency of a generator is determined by the atomic structure of the semiconductor legs, the design of the leg, thermocouple or module, the heat exchangers, and the circuit of the electricity harnessing unit. In the current generator design, the atomic structure of a semiconductor leg plays a limiting role in determining the generator’s efficiency. The quality of a semiconductor as a thermoelectric leg is measured by a dimensionless parameter called “figure of merit” defined as \({zT}=\alpha ^2\,\sigma \,T/\kappa\), where \(\alpha\) is the Seebeck coefficient, \(\sigma\) the electrical conductivity, and \(\kappa\) the thermal conductivity. For a thermocouple or a module, it is transformed into \({{ZT}}=\alpha ^2\,T/(R\,K)\), where R is the electrical resistance and K the thermal conductance, to take into account of the effects of its geometry and size.

A typical flat planar thermoelectric generator. It is usually assembled from a numbers of thermocouples made of a p- and an n-type semiconductor legs. These legs are connected through interconnectors (metal strips in brown) and sandwiched between a hot and a cold electrical insulators (ceramics in gray) so that the electric current flows in series and the heat current flows in parallel through these legs. These modules are sandwiched between heat exchangers (not shown) that supply heat to the hot side and remove the exit heat from the cold side to maintain a temperature difference between the two sides. The external terminals of these modules are connected to an electrical load (bulb) to form a closed system for electricity generation.

For a thermoelectric leg with temperature-independent transport properties, its efficiency under a given temperature difference can be expressed as a function of its zT by:

where \({T_h}\) and \({T_c}\) are the absolute temperatures at the hot- and the cold-ends of the leg, respectively, and \(T_m\) is the mean temperature measured between the two ends. Clearly, a large zT of a thermoelectric leg is crucial for the efficiency of a generator and is the ultimate goal in the search for a good thermoelectric material besides its abundance and cost. “Untuned” bulk thermoelectric materials have the zTs much less than 1.0 and can hardly be useful for high-efficient power generation. These materials have to be “tuned” to have (1) large power factor \(\alpha ^2\,\sigma\) and (2) lower thermal conductivity \(\kappa\) to attain large zTs. However, it has proved to be a challenging task to meet both of these two requirements in a material since its transport properties are intrinsically related and are functions of the band structures, doping concentrations, and scattering processes among others. As it will be showed in later section that \(\alpha\) and \(\sigma\) generally vary in a reciprocal manner, which makes any improvement on zT difficult.

Clearly, an ideal thermoelectric material should have a large Seebeck coefficient and behave like a crystalline solid for its charge carrier transport and like an amorphous solid for its phonon transport, so called “phonon glass/electron crystal, PGEC3,4. Through many decades’ search, a deep insight of the intrinsic nature of a good thermoelectric material has been gained, theories, various synthesis technologies, and treatment techniques have been developed toward realizing the ideal thermoelectric material. These include low-dimensional theory5,6, band engineering, structural modification, phase manipulation, etc. (for a comprehensive review in this regard, refer to Ref7). These efforts have led to the developments of many new materials and the achievement of zT as high as 3.1. These advances in the quality of thermoelectric materials represented by the improvement made on zTs have been compiled into historical charts for typical materials7,8 . Based on these charts, the upper envelopes of the zTs of several classes of thermoelectric materials can be drawn as shown in Fig. 2.

Upper zT-envelopes constructed from the existing data for: SnSe9–11, Cu2Se12,13,14,15,16, GeTe17,18,2.3 in Ge $$_{1-x-y}$$ 1 - x - y Sb $$_x$$ x In $$_y$$ y Te via reducing the phase-transition temperature and introducing resonant energy doping. Adv. Mater. 30, 1705942 (2018)." href="#ref-CR19" id="ref-link-section-d50013821e1283_2">19,20,21,22, PbTe23,24,25,26,27, PbS28,29,30,31, SnTe32,33,34,35,36,37,38,39, Bi\(_2\)Te\(_3\)40,41,42,43,44, Skutterudites45,46,47,48,49, BiCuSeO50,51,52, SiGe53,54,55,56, Half-Heuslers700 $$^0\,$$ 0 C) p-type thermoelectric materials. Appl. Phys. Lett. 93, 022105 (2008)." href="#ref-CR57" id="ref-link-section-d50013821e1384">57,58,59,60,61. The lines are only for differentiating one group of materials from the other.

Translating high zT thermoelectric materials into a high-efficient generator has proven to be challenging not only because the transport properties of both the n- and the p-type semiconductor legs are temperature-dependent but also because their optimum zTs are different both in the magnitude and at the temperature level and are only attainable within a small fraction of the operating temperature range of the device. Furthermore, construction of a generator needs other components such as the interconnectors, thermal insulating ceramics, heat exchangers, and electricity extracting unit. These components introduce additional electrical and thermal resistance that have to be minimized to make the device as efficient as its legs. Consequently, module design plays an important role in realizing the efficiency of thermoelectric legs at the generator level and currently lags far behind the advances of thermoelectric materials.

Design optimizations so far are limited to the parameters including leg length, cross-sectional area, and shape; module fill factor and the electrical and thermal resistance induced by the interfaces between the leg and the interconnector and between the interconnector and the thermal insulator. It was found that an optimal ZT could be obtained when the relation \({A_p}{L_n}/({A_n} {L_p}) = \sqrt{{{\kappa }_n}{{\rho }_p}/({{\kappa }_p}{{\rho }_n})}\) for the two legs was satisfied62, where A and L are, respectively, the cross-sectional area and length of the leg and the subscripts n and p are leg types. Numerical simulation showed that this analytical relation was approximately true even when the temperature-dependent transport properties of the materials and the electrical and thermal contact resistances were taken into account63. The numerical simulation further showed that a similar power output per unit module area obtained with a greater number of longer legs could be realized with a smaller number of shorter legs when the thermal contact resistance was neglected.

Early investigation64 indicated that different shapes of a constant cross-sectional area along a leg had little effect on power generation and efficiency of the module. However, a leg with variable cross-sectional area showed improved power generation and efficiency65,66. Among different leg shapes investigated, pyramidal leg66,67 showed a higher output power density than cuboid, exponential68,69, and quadratic shaped legs70. Subsequently, a laboratory-scale thermoelectric module with pyramid-shaped leg was recently built for proof-of-concept71. The experimental data confirmed the analytical prediction. A pyramid-shaped leg can increase the temperature gradient along the leg, thus the Seebeck voltage, by lowering its thermal conductance. Furthermore, it can take into account of the Thomson effect which is usually ignored in the leg with a constant cross-sectional area. The laboratory-scale module showed a 67% increase in power generation compared with the module having cuboid legs. More leg shapes , including cuboid, trapezoid, hourglass, Y-shaped, and hollow legs and their different arrangements were modeled using the finite element method72. The numerical results showed that the hourglass shaped leg had more than doubled electrical potential and the maximum power output compared to the cuboid leg when a constant temperature is maintained at the hot end and a natural convection is applied at the cold end.

Three approaches, i.e., functionally graded thermoelectric material (FGTEM), segmented thermoelectric legs, and cascaded thermoelectric modules, are commonly used to improve the low average zT within a large operating temperature range. The fundamental concept of FGTEM is to synthesize a leg material to have an optimum zT all the way from the hot-end to the cold-end through locally selecting a particular material, microstructure, or composition. Several techniques, including tuning the amount of dopant at different spatial positions of73 and creating a carrier concentration or a composition gradient in a material73,74,75,76,77,78,79,80,81,82, have been proposed to synthesize FGTEMs. FGTEM has the advantages to have a gradual phase transition without inducing any interface resistances. However, the synthesizing technology so far is challenging. A segmented FGTEM that joins the existing thermoelectric materials to fit the zT profile of a thermoelectric leg with their optimum zTs was proposed by Ioffe in his patent application83 and gains increased interest since83,84,85,86,87. Several promising designs88,89 have been studied and showed improved efficiencies. A shortcoming of the segmented legs is the presence of the large interface resistances at the junctions between consecutive sections. Furthermore, the materials of different segments have to be compatible to make a high-efficient leg90,91. Cascade modules can avoid the electrical contact resistance problem by stacking multiple stages of thermocouples with different materials and using separated electric circuit for each stage to improve its performance within large temperature difference. Cascade module is complicated to implement.

Clearly, current device design optimizations can only produce an optimal efficiency within the range limited by its leg’s zT which is yet below the anticipated level, \(zT>3\) which is needed to make the device potentially competitive with the traditional thermal-process-based energy conversion technology. In this work, a new design approach is introduced to develop a thermoelectric generator with conversion efficiency and power output that can potentially exceed those of the legs’ materials. The existence of an optimal combination of the transport properties of a thermoelectric material at low doping concentration under an ideal condition is firstly illustrated. A new thermocouple design that uses low doping legs and high doping semiconductors as external carrier injectors surrounding the legs is then introduced. Finally, an example is provided to demonstrate the potential of the new approach in enhancing the device efficiency.

In addition to depicting the complex nature of the intrinsically dependent transport properties as stated in the Introduction, Fig. 3a also shows clearly that lower doping concentration yields two best transport properties, i.e., Seebeck coefficient and thermal conductivity, while higher doping concentration results in only one best transport property, electrical conductivity, for thermoelectric applications. Such a pattern is also observed in other materials92. In this regard, using low doping concentration as a starting material to optimize its thermal-to-electrical conversion efficiency could have certain advantages over the high-doping concentration approach. However, low doping material has limited charge carriers to work with to improve its electrical conductivity. This could be a major barrier in searching for an ideal thermoelectric material with low doping concentration. In this section, it will be shown that an optimal combination of the transport properties also exists at lower doping concentration for certain materials with their transport properties have the similar trends as those shown in Fig. 3a.

Optimizing zT through tuning carrier concentration. (a). Thermoelectric transport properties as functions of carrier concentration. Reproduced from Ref.93. (b). Shifting the \(\sigma\)-curve by 10\(^{18}\) cm\(^{-3}\) toward lower doping concentration without disturbing the \(\alpha\)- and \(\kappa\)-curves leading to a 83% increase in the electrical conductivity, electric current, power factor, and zT, and a 68% increase in the conversion efficiency at \(N_D=10^{18}\). (c). Shifting the \(\sigma\)-curve by \(3\cdot 10^{18}\) cm\(^{-3}\) leading to a more than twofold increase in the electrical conductivity, electric current, power factor, and zT, and a more than one and half fold increase in the conversion efficiency. (d). Shifting the \(\sigma\)-curve by \(5\cdot 10^{18}\) cm\(^{-3}\) leading to a more than three and half fold increase in the electrical conductivity, electric current, power factor, and zT and nearly two and half increase in the conversion efficiency. It can be noticed that the optimum zT also shifted toward low doping concentration.

Assuming the electrical conductivity curve shown in Fig. 3a can be rigidly shifted toward lower doping concentration without disturbing the Seebeck coefficient and the thermal conductivity curves, a new set of transport properties can be obtained as shown in Fig. 3. Figure 3b to d are the results obtained by shifting the \(\sigma\)-curve by \(10^{18}\), \(3\cdot 10^{18}\), and \(5\cdot 10^{18}\) cm\(^{-3}\), respectively, toward low doping concentration. For a thermoelectric leg with \(A_n\)=0.267 cm\(^2\), \(L_n\)=2 mm, and initial doping concentration \(N_D=10^{18}\) cm\(^{-3}\), the improvemed electrical conductivity, power factor, and zT are extracted from those shifted curves shown in Fig. 3 and are presented in Table 1. The intrinsic electrical resistance is estimated by \({R_I}=L_n/\sigma /A_n\). The load resistance is estimated by \({R_L}=10^2/\sigma\) to approximately match the intrinsic electrical resistance. The electric current is then calculated by \({I}=\alpha \,{\Delta T}/(R_L+R_I)\). The conversion efficiency is determined by Eq. (1). The electrical conductivity increased by about 83% for the case of \({\Delta n}=10^{18}\), by about 238% for \(\Delta n=3\cdot 10^{18}\), and by about 378% for \(\Delta n=5\cdot 10^{18}\), respectively. Similar increase scales are also found for electric current, power factor, and figure of merit. The conversion efficiency increased by about 68% from 0.026 to 0.044 in the first case, by about 169% from 0.026 to 0.07 in the second case, and by about 242% from 0.026 to 0.089 in the third case. It can be noticed from Fig. 3d that the optimum zT is also shifted toward low doping concentration, which indicates that the optimum zT can also be attainable in low doping materials through tuning the electrical conductivity.

Shifting \(\sigma\)-curve without disturbing \(\alpha\)- and \(\kappa\)-curves implies that the charge carriers are injected into the leg without altering the lattice structure of the material. Consequently, the Seebeck coefficient and the thermal conductivity remain the same as those before shifting the \(\sigma\)-curve. Applications of such a concept can be found in most of the semiconductor devices, typically in bipolar junction transistors (BJTs) and metal-oxide-semiconductor field effect transistors (MOSFETs). In BJT, nearly 100% of the charge carriers of emitter is injected into the collector through the two p-n junctions to increase the collector’s electric current without altering its lattice structure. In MOSFET, a significant amount of carriers of the n- or p-channel is injected into the drain to increase its electric current without altering its lattice structure. Injecting charge carriers into a thermoelectric leg without disturbing its lattice structure is one of the aspects on which the new design concept is based. The another aspect is to use lower doping concentration as the starting material that has both the best Seebeck coefficient and the best thermal conductivity, leaving the electrical conductivity to be improved by injecting carriers from the leg’s surrounding boundary using higher doping concentration materials. The theory behind this concept can be reasoned by the following three processes.

First, when a semiconductor with high carrier concentration is joined with a leg that has lower carrier concentration, the concentration gradient at the interface causes the carriers to diffuse from the higher concentration material to the lower concentration leg. As a result, the carrier concentration near the boundary of the leg increases, which leads to higher electric current density. From Drude’s model, \(\sigma ={q}{n_e}{\mu }\), where q is electric charge, \(n_e\) carrier concentration, and \(\mu\) carrier mobility, it is known that the electrical conductivity increases as the carrier concentration increases. Since the higher carrier concentration only exists near the leg boundary, the interaction between these injected carriers and the carriers and phonons of the leg is supposed to be minimal. It can be reasonably anticipated that the injected carriers cause little disturbance to the leg’s lattice structure, thus the Seebeck coefficient and the thermal conductivity remain unchanged. Consequently, the figure of merit of the leg increases with the increase of injected carrier concentration.

Second, as carriers continue to move across the interface into the leg, some uncompensated positive ions remain near the interface on the side of the material with high carrier concentration. An electric field with a direction from the positive charge toward the negative charge is build-up at the interface. This electric field creates a barrier potential preventing further carrier diffusion until an equilibrium state is reached. This electric field pulls the injected carriers toward leg’s exterior boundary. As a result, the carrier concentration at the boundary is significant higher than that of the leg, which can prevents electrons of a n-type leg move toward the boundary where, otherwise, an infinitive carrier recombination rate exists. Consequently, this carrier concentration gradient reduces the carrier recombination rate at the boundary such that electron mobility near the boundary is improved. From Drude’s model, it is also known that the electrical conductivity increases with electron mobility.

Third, as the higher carrier concentration material is joined with the leg at the leg’s exterior boundary, the originally free boundary of the leg becomes constrained. Thus, its phonon waves also become constrained at the boundary. Consequently, the higher carrier concentration material can be an additional freedom to modulate the lattice vibration mode of the leg to further lower its lattice thermal conductivity.

To investigate the effects of the injected carriers on the transport properties and on the efficiency of a thermocouple, a common-emitter bipolar junction transistor is used in this work to substitute the high-carrier-concentration semiconductor. The BJTs are joined or embedded on the exterior boundaries of both the n- and the p-type legs. Figure 4 shows a schematics of the design of the new thermocouple. By adjusting the base current of the transistor, the amount of the injected carriers into the thermoelectric legs can be modulated to obtain an optimal injected carrier concentration that maximizes the efficiency and power output of the thermocouple.

(a) Schematics of the new thermocouple design including a p-type and an n-type thermoelectric legs, conductor strips at the hot end and the cold ends and a electrical load resistance, \(R_L\). An n-p-n bipolar junction resistor (not to scale) is joined with the n-type leg at the exterior surface and a p-n-p BJT is joined with the p-type leg. These two BJTs are common-emitter configured. \({V_{be}}\) is the bias between the base and the emitter. \({V_{ce}}\) is the bias between the collector and the emitter. \({R_{nc}}\) and \({R_{pc}}\) are the load resistance for the n-p-n and p-n-p transistors, respectively. (b) Symbols of an n-p-n transistor. (c) Symbols of a p-n-p transistor.

Models to determine the current-voltage characteristics of an n-p-n bipolar transistor is briefly introduced in this section (see Methods for details and variable definitions). A bipolar junction transistor is made of three regions and two p-n junctions as shown in Fig. 4b for an n-p-n transistor and 4c for a p-n-p transistor. These three regions are designated as: emitter, base, and collector with the base sandwiched between the emitter and the collector. The emitter and the base regions form the first p-n junction. The second p-n junction is formed between the base and the collector regions. The emitter has the highest doping, the collector has the lowest doping. The doping concentration in the base is smaller than that in the emitter but higher than that in the collector. Three terminals connect, respectively, to each region. Based on the bias applied at these terminals, a BJT can be operated in four different modes: forward active, saturation, cutoff, and inverted. For each operating mode, the bias sets the level of barrier potentials at the two metallurgical junctions to modulate the carrier flow from one region to another.

In the new thermocouple design, the normal operating condition, i.e., the forward active mode is used, in which the emitter-base junction is forward biased, whereas the collector-based junction is reverse biased. Based on ideal p-n junction theory, the emitter, base, and collector currents, \({I_e}\), \({I_b}\), and \({I_c}\), respectively, can be expressed in the following compact form93,94,96

where

When \(I_b\) and \(V_{bc}\) are given, the common-emitter collector current can be obtained from Eq. (2) and expressed as:

The current gain of a common-base bipolar junction transistor, denoted as \({\alpha _F}\), is defined as the ratio of collector current to emitter current, i.e.:

where \({\alpha _T}\) is defined as the base transport factor and \({\gamma }\) is the emitter injection efficiency. These current gain factors can be evaluated by:

The common-emitter current gain, \({\beta }\), is defined as the ratio of collector current to base current, i.e., \(\beta = {\partial I_c}/{\partial I_b}\). Using Eq. (2c), it can be expressed in the common-base current gain as:

In most modern bipolar junction transistors, \(\alpha _F \approx 1\), hence \(\beta\) is significantly larger than 1.

A schematics of a simple amplifier. It is composed of a common-emitter n-p-n bipolar transistor and two resistors. \({V_{cc}}\) is the power supply hold at a constant positive voltage. \({V_{in}}\) is the input signal. \(V_{out}=V_{ce}\) is the output signal. The \(V_{cc}\), \({R_c}\), collector, and emitter form an output loop. The \({V_{in}}\), \({R_b}\), base, and emitter form an input loop.

Voltages and currents can be amplified by a bipolar transistor in combination with other components. This amplification functionality of a bipolar transistor is used here to modulate the amount of charge carriers injected into the thermoelectric legs in the new thermocouple design. Figure 5 shows a schematics of a simple amplifier isolated from Fig. 4. The amplifier consists of an n-p-n bipolar transistor in a common-emitter configuration and two resistors to modulate the voltages at the base and the collector. \(V_{cc}\) is the supplied positive voltage hold at a constant value. The input signal is delivered by the source voltage, \(V_{in}\). The output signal, \(V_{out}\) which is equal to \(V_{ce}\), is measured between collector and emitter terminals. Using Ohm’s and Kirchhoff’s laws, the relationship between the output voltage and the input voltage can be determined, respectively, from the collector-emitter loop and the base-emitter loop as:

For a common-emitter transistor in forward active mode, \(I_c=\beta \,I_b\) and \(V_{be}\) is a constant. The amplified voltage, i.e., \(V_{out}\) can be expressed as97:

Since \(V_{be}\) is a constant, any change in the input signal \(V_{in}\) leads to a change in the output signal \({V_{out}}\). The variation of the output voltage is proportional to the change of the input voltage by an amplification factor expressed as:

Multiplying the second equation of Eq. (6) by the collector current, the supplied power distribution between the resistor \(R_c\) and the transistor is obtained as:

where \(I_c\,V_{cc}\) is the supplied power, \(R_c\,I_c^2\) is the power dissipated in the resistor, and \(I_c\,V_{ce}\) is the power used by the transistor to amplify the input signal. The later term is the power used in the new thermocouple design to improve the leg’s electrical conductivity and the terminal voltage.

From heat and electric current coupled equations, heat balance equations can be reduced for a thermoelectric leg. The heat balance equations at the hot and the cold ends of the leg are, respectively, expressed as98,99:

The voltage of the leg between the hot and the cold terminals is given by \(V_{TEG}=\alpha \,\Delta T\), where \(\Delta T=T_h-T_c\). The electric current flowing through the leg is calculated by Ohm’s law as \(I_{TEG}=V_{TEG}/(R_I+R_L)\), where \(R_I\) and \(R_L\) are, respectively, the intrinsic and load resistances. The output power of the leg is:

For the leg of the new thermocouple, there exist two electric currents flowing through it. One is from the thermocouple loop and the other is from the bipolar transistor loop as shown in Figs 4 and 5. As it is pointed out earlier that the power \(I_c\,V_{ce}\) used by the transistor to amplify the input signal can improve both the electrical conductivity and the terminal voltage in the new thermocouple design. As a first-order approximation, the current and the voltage of the leg can be expressed as:

where \({f_{I}}\) and \({f_{V}}\) are, respectively, the fractions of the collector current and the voltage between collector and emitter of the transistor, contributing to the output power of the leg. The output power of the new leg is then computed by:

Introducing Eq. (7) into Eq. (13a), the output power can be expressed as:

Since both the second and the third terms of Eq. (13b) are larger than 0, these two terms can be represented by \(\delta P (> 0)\). Equation (13b) can then be approximated by:

where \({\alpha _P}\) is denoted as the power gain of the leg.

The improvement in the output power can be measured by the power gain rate defined as:

For a single thermocouple, \(I_{TEG}\) is usually smaller than \(I_c\). If the term \(f_I\,I_c\) in Eq. (14) can be designed to be greater than \(I_{TEG}\), the power gain rate can be greater than 3. If the heat rate used to generate \(P_0\) and \(P_e\) are the same, the efficient gain rate would be the same as \(\eta _P\).

The ratio between the output power and the supplied power, i.e.,

can be considered as the effectiveness of the leg for power generation. The effectiveness depends on the fractions of the collector current and the collector-emitter voltage contributed to the output power of the leg. It is worth to note that \(I_c\) is the total collector current injected by the emitter. If the transistor is designed such that large portion of the injected carriers enters the thermoelectric loop, the effectiveness of the leg would be large. In addition, the effectiveness is only relevant for the design using a bipolar transistor. When charge carriers are injected through high-doping semiconductors joined with the leg at the leg’s exterior boundary, the power supply is not needed.

A simple example is given here to illustrate the potential of using external carrier injection in improving the efficiency of a thermoelectric leg. The carriers are injected using a bipolar junction transistor as shown in Fig. 4. The uni-leg thermocouple is made of n-type Silicon semiconductor and has a cross-sectional area (\(A_n\)) of 0.267 cm\(^2\), a length (\(L_n\)) of 2 mm, and an initial doping concentration (\(N_D\)) of 10\(^{15}\)cm\(^{-3}\). The temperature is 300 K at the cold end and 500 K at the hot end. The bipolar transistor is embedded on the exterior surface of the leg. A schematics of the transistor is shown in Fig. 5. The parameters used for each region of the n-p-n transistor are given in Table 2. These parameters are estimated based on those used in Ref94,95. The cross-sectional area (\(A_t\)) of the transistor is 5 \(\mu\)m\(^2\). The supplied voltage \(V_{cc}\) is fixed at 3.5 V. Three load resistance (\(R_c\)), 3 \(\Omega\), 1.95 \(\Omega\), and 1.35 \(\Omega\) are used for the collector-emitter loop. The load resistance (\(R_b\)) for the emitter-base loop is 1 k\(\Omega\). The electric current flowing through the thermoelectric loop, including the leg and a load resistance \(R_L\), is 0.09 A calculated using a thermoelectric model. The calculated mean electrical conductivity of the leg is 32 S/m. The mean intrinsic resistance of the leg is 2.34 \(\Omega\) estimated by \(R_I=L_n/(\sigma \,A_n)\). The load resistance is \(R_L=3.12\,\Omega\) calculated by \(10^2/\sigma\) to approximately match the leg intrinsic resistance.

The three electric current equations given by Eqs. (2a), (2b), (2c) and the voltage equations given by Eqs. (6) plus the relation of these voltages given by:

are used to determine the collector current along the load line defined by the second equation of (6). Substituting \(V_{ce}\) of Eq. (16) into the second equation of (6), the supplied voltage can be expressed as:

Introducing \(V_{be}\) obtained from Eq. (2b) and \(I_c\) of Eq. (3) into Eq. (17), the interception point between current-voltage characteristic curve and the load line can be determined for a given base current, \(I_b\).

Figure 6 shows the current-voltage curve, the load line, and their interception points for different base currents. The interception point with maximum \(I_c\) and \(I_b\) is the saturation point below which both the base-emitter and the base-collector junctions are forward biased and the transistor is in saturation mode. The interception point with \(I_c=0\) and \(V_{ce}=V_{cc}\) is the cutoff point. All the points intercepted by different \(I_b\) currents between the saturation and the cutoff points are the quiescent points (Q-point). A transistor always operates at the quiescent points for the given \(V_{cc}\) and \(R_c\).

Current-voltage characteristics, load line and their interception points for several specific base currents \(I_b\)s. The collector doping concentration used to obtain these results is \(N_D=2.5\cdot \,10^{15}\)cm\(^{-3}\). The rest parameters of the transistor are given in the text.

Substituting \(V_{be}\) of the first equation of (6) into Eq. (17) and using Eq. (3) for \(I_b\), the common-emitter collector current along the load line can be determined by:

Setting \(V_{be}\) of Eq. (2b) equals to that of the first equation of (6) and in combination with Eq. (18), \(V_{bc}\) along the load line can be obtained by solving the following equation:

Figure 7 shows the load lines, the output power, and the power gain rate of the thermoelectric leg obtained using Eqs. (18), (19), (13), (14), and (15) for 3 different \(R_c\)s: 3.2 \(\Omega\). 1.95 \(\Omega\), and 1.35 \(\Omega\). Figure 7a shows the common-emitter collector currents along the three load lines. As it is known from Fig. 6, the base current, \(I_b\), varies along load line and \(I_b\) decreases with increase of the collector-emitter voltage, \(V_{ce}\). Based on this knowledge, it can be recognized from Fig. 7a that \(I_c\) level depends on both the base current \(I_b\) and the slope (\(1/R_c\)) of the load line. Larger \(I_b\) and smaller \(R_c\) produce a larger \(I_c\). Figure 7b shows the power generated by this new thermoelectric leg. The amount of the generated power has a similar trend as the load line but nonlinear, that is, increases nonlinearly with \(I_b\) and \(1/R_c\). Figure 7b also shows the power generated solely through Seebeck effect, \(P_0\), in gray line. Figure 7c shows the power gain rate defined by Eq. (14). As it is showed in Fig. 7a that the magnitude of \(I_c\) is significantly larger than the thermoelectric current \(I_{TEG}\) due to Seebeck effect, therefore, the power gain rate is significant. The power gain rate also depends on the base current (\(I_b\)) and the slope of the load line (\(1/R_c\)). Figure 7d shows the ratio between the generated power and supplied power, i.e., \(P_e/[V_{cc}\,(1-f_I)\,I_c]\). As noted earlier, the collector current is significantly larger than the thermoelectric current \(I_{TEG}\) of the leg and the fraction that \(I_c\) contributed to \(P_e\) is small (\(f_I=0.1\) is assumed). Consequently, the effectiveness is smaller for most of the \(I_c\) on the load line. For smaller \(I_c\) near the cutoff point of the load line, the effectiveness can be larger than 1 but the power gain rate is low. Evidently, there exists a trade off between the power gain ratio and the effectiveness when designing such a thermocouple. It worth to emphases that when high doping semiconductor is used as the carrier injector rather than the BJTs, the effectiveness parameter is irrelevant.

Performance of the new thermoelectric leg. (a) Load lines of the amplifier obtained for three \(R_c\)s. (b) Power generated by the thermoelectric leg. \(P_0\) power produced solely through Seebeck effect. (c) Power gain rate. (d) The effectiveness of the new thermoelectric leg. The parameters of the thermoelectric leg and the amplifier used to obtain these results are given in the text.

From the relations between thermoelectric properties (Seebeck coefficient, electrical conductivity, thermal conductivity, power factor, and figure of merit) and doping concentration, an optimal figure of merit can be attainable by rigidly shifting the electrical conductivity curve toward low doping concentration without disturbing the Seebeck coefficient and the thermal conductivity curves. This is in contrast to the current high-doping concentration approach in finding an optimal figure of merit among various class of materials.

The proposed new thermocouple design that uses low doping Silicon as the thermoelectric legs and bipolar junction transistors as the external carrier injectors can triple the power output and efficiency of the original design when the electric current in the leg contributed by the injected carriers is equal to or greater than the current induced by the Seebeck effect. The improvement on the power output and the efficiency of the new thermocouple significantly depends on the base current and the load line slope of the transistor.

This work presents a typical example to improve the figure of merit of a thermocouple through its design such that its efficiency and power output can potentially exceed those of the materials of the thermoelectric legs. Using low-doping material as thermoelectric leg and high-doping semiconductors as external carrier injectors is a new concept for thermocouple design and can potentially lead to the development of new thermoelectric materials and high-efficient thermoelectric modules.

Under normal operating condition, i.e., the forward active mode, the emitter-base junction is forward biased, whereas the collector-based junction is reverse biased. In this mode, the electron concentration in the base at the emitter-base junction is much higher than the equilibrium concentration of minority carriers in the base. Based on ideal p-n junction theory, the electron concentration in the base at the emitter-base junction can be expressed as:

where \({n_{b}}\) is the electron concentration in the base, \({x_{e}}\) the boundary of the depletion region of the emitter-base junction at the emitter side, \({V_{th}}=k_B\,T/q\) is the thermal voltage, \(V_{be}\) is the voltage drop across the junction (\(V_{be} \gg V_{th}\)). The equilibrium concentration of minority carriers in the base, \({n_{b0}}\), is calculated by \(n_{b0}=n_{ib}^2/N_{Ab}\), where \({n_{ib}}\) and \(N_{Ab}\) are, respectively, the intrinsic carrier and the acceptor concentrations of the base.

To the contrary, the electron concentration in the base at the collector-base junction is much smaller than the equilibrium concentration of minority carriers in the base and is given by:

where \(V_{bc}\) is the collector-base voltage (\(V_{bc} < 0\) for forward active mode), and W is the base width. Here, it is assumed that the widths of the depletion regions at the emitter-base and collector-base interface are sufficiently small and can be neglected comparing to the base width.

The distribution of minority carrier in the base region can be determined by the one-dimensional steady-state continuity equation:

in which \({\Delta n_{b}}(x)=n_b(x)-n_{b0}\) and \({L_{nb}}=\sqrt{D_{nb}\,\tau _{nb}}\) is the diffusion length, where \({D_{nb}}\) and \({\tau _{nb}}\) are, respectively, the diffusivity and the lifetime of electron in the base.

The solution of the continuity equation, Eq. (22), is straightforward. For the two boundaries given by Eqs. (20) and (21), the solution can be written as:

Using a one-dimensional continuity equation identical to Eq. (22) and using the minority carrier concentrations at the boundaries of the emitter: \(p_e(0)={p_{e0}}\) and \(p_e(x_e)=p_{e0}\,e^{V_{be}/V_{th}}\), the distribution of hole concentration in the emitter, \({p_{e}}(x)\), can be determined by:

where \({\Delta p_{e}}(x)=p_e(x)-p_{e0}\), \({L_{pe}}=\sqrt{D_{pe}\,\tau _{pe}}\) is the hole diffusion length, in which \({D_{pe}}\) and \(\tau _{pe}\) are, respectively, the diffusivity and lifetime of hole in the emitter, \(p_{e0}={p_{ie}}^2/N_{De}\), \(p_{ie}\) the intrinsic hole concentration, and \(N_{De}\) the donor concentration in the emitter.

Similarly, using the minority carrier concentrations at the boundaries of the collector: \(p_c(x_e+W)=p_{c0}\,e^{V_{bc}/V_{th}}\) and \({p_{c}}(x_e+W+{x_c})={p_{c0}}\), the hole concentration in the collector, \(p_c(x)\), can be found and is given by:

where \({\Delta p_{c}}(x)=p_c(x)-p_{c0}\), \({L_{pc}}=\sqrt{D_{pc}\,\tau _{pc}}\) is the hole diffusion length, in which \({D_{pc}}\) and \({\tau _{pc}}\) are, respectively, the diffusivity and lifetime of hole in the collector, \(p_{c0}={p_{ic}}^2/N_{Dc}\) the equilibrium hole concentration, and \(N_{Dc}\) the donor concentration of the collector.

Electric current induced by carrier diffusion is proportional to the spatial derivative of the carrier concentration. The diffusion current in the base at the emitter edge is thus given by:

where S is the cross-sectional area of the base region.

The diffusion current in the base at the collector edge is determined by:

The electric current caused by hole diffusion in the emitter at the emitter-base junction is determined by:

Similarly, the electric current caused by hole diffusion in the collector at the collector-base junction is calculated by:

These emitter, base, and collector currents can be obtained and expressed in the following compact form:

where

When \(I_b\) and \(V_{bc}\) are given, the common-emitter collector current can be obtained from Eqs. (30a), (30b), and (30c) as:

Current gain of a common-base transistor

Seebeck coefficient

Power gain of the new thermoelectric module

Base transport factor

Common-emitter current gain

Thermal-to-electrical energy conversion efficiency

Effectiveness of the new thermoelectric module

Voltage ratio: qVbc/(kBT) with kB being Boltzmann constant

Voltage ratio: qVbc/(kBT)

Power gain rate of the new thermoelectric module

Emitter injection efficiency

Thermal conductivity

Thermal conductivity of an n-type leg

Thermal conductivity of a p-type leg

Carrier mobility

Electrical resistivity of an n-type leg

Electrical resistivity of a p-type leg

Electrical conductivity

Lifetime of electron in the base

Lifetime of hole in the collector

Lifetime of hole in the emitter

EURELECTRIC. Efficiency in Electricity Generation (EURELECTRIC, Brussels, 2003).

Google Scholar

Energetics-Incorporated. Manufacturing energy and carbon footprints. http://energy.gov/eere/amo/manufacturing-energy-and-carbon-footprints-2010-mecs [Accessed 15 April 2014] (2010).

Slack, G. A. New materials and performance limits for thermoelectric cooling. In CRC Handbook of Thermoelectrics (ed. Rowe, D. M.) 407 (CRC Press, Boca Raton, 1995).

Google Scholar

Nolas, G. S., Morelli, D. T. & Tritt, T. M. Skutterudites: A phonon-glass-electron-crystal approach to advanced thermoelectric energy conversion applications. Annu. Rev. Mater. Sci. 29, 89–116 (1999).

Article ADS CAS Google Scholar

Hicks, L. D. & Dresselhaus, M. Figure of merit of a one-dimensional conductor. Phys. Rev. B 47, 16631–34 (1993).

Article ADS CAS Google Scholar

Hicks, L. D. & Dresselhaus, M. Effect of quantum well structures on the thermoelectric figure of merit. Phys. Rev. B 47, 12727–31 (1993).

Article ADS CAS Google Scholar

Shi, X. L., Zou, J. & Chen, Z.-G. Advanced thermoelectric design: From materials and structures to devices. Chem. Rev. 120, 7399–515 (2020).

Article CAS PubMed Google Scholar

Zhang, X. & Chen, L. D. Thermoelectric materials: Energy conversion between heat and electricity. J Materiomics 1, 92–105 (2015).

Article Google Scholar

Zhao, L. D. et al. Ultralow thermal conductivity and high thermoelectric figure of merit in SnSe crystals. Nature 508, 373 (2014).

Article ADS CAS PubMed Google Scholar

Chang, C. et al. 3D charge and 2D phonon transports leading to high out-of-plane ZT in n-type SnSe crystals. Science 360, 778–783 (2018).

Article CAS PubMed Google Scholar

Zhou, C. et al. Polycrystalline SnSe with a thermoelectric figure of merit greater than the single crystal. Nat. Mater. DOI:https://doi.org/10.1038/s41563-021-01064-6 (2021).

Article PubMed PubMed Central Google Scholar

Liu, H. L. et al. Copper ion liquid-like thermoelectrics. Nat. Mater. 11, 422–425 (2012).

Article ADS PubMed CAS Google Scholar

Yu, B. et al. Thermoelectric properties of copper selenide with ordered selenium layer and disordered copper layer. Nano Energy 1, 472–478 (2012).

Article CAS Google Scholar

Liu, H. L. et al. Ultrahigh thermoelectric performance by electron and phonon critical scattering in Cu\(_{2}\)Se\(_{1-x}\)I\(_x\). Adv. Mater. 25, 6607–12 (2013).

Article CAS PubMed Google Scholar

Zhong, B. et al. High superionic conduction arising from aligned large lamellae and large figure of merit in bulk Cu\(_{1.94}\)Al\(_{0.02}\)Se. Appl Phys Lett 105 (2014).

Yang, D. W. et al. Blocking ion migration stabilizes the high thermoelectric performance in Cu\(_2\)Se composites. Adv. Mater. 32, e2003730 (2020).

Article PubMed CAS Google Scholar

Gelbstein, Y., Davidow, J., Girard, S. N., Chumg, D. Y. & Kanatzidis, M. Controlling metallurgical phase separation reactions of the Ge\(_{0.87}\)Pb\(_{0.13}\)Te alloy for high thermoelectric performance. Adv. Energy Mater. 3, 815–820 (2013).

Article CAS Google Scholar

Li, J. et al. Simultaneous optimization of carrier concentration and alloy scattering for ultrahigh performance GeTe thermoelectrics. Adv. Sci. 4, 1700341 (2017).

Article CAS Google Scholar

Hong, M. et al. Realizing zT>2.3 in Ge\(_{1-x-y}\)Sb\(_x\)In\(_y\)Te via reducing the phase-transition temperature and introducing resonant energy doping. Adv. Mater. 30, 1705942 (2018).

Article CAS Google Scholar

Hong, M. et al. Arrays of planar vacancies in superior thermoelectric Ge\(_{1-x-y}\)Cd\(_x\)Bi\(_y\)Te with band convergence. Adv. Energy Mater. 8, 1801837 (2018).

Article CAS Google Scholar

Li, J. et al. High-performance GeTe thermoelectrics in both rhombohedral and cubic phases. J. Am. Chem. Soc. 140, 16190–16197 (2018).

Article CAS PubMed Google Scholar

Bu, Z. L. et al. Realizing a 14% single-leg thermoelectric efficiency in GeTe alloys. Sci. Adv. 7, eabf2738 (2021).

Article ADS CAS PubMed PubMed Central Google Scholar

Hsu, K. F. et al. Cubic AgPb\(_m\)SbTe\(_{2+m}\): Bulk thermoelectric materials with high figure of merit. Science 303, 818–21 (2004).

Article ADS CAS PubMed Google Scholar

Biswas, K. et al. High-performance bulk thermoelectrics with all-scale hierarchical architectures. Nature 489, 414–8 (2012).

Article ADS CAS PubMed Google Scholar

Wu, D. et al. Superior thermoelectric performance in PbTe–PbS pseudobinary: Extremely low thermal conductivity and modulated carrier concentration. Energy Environ. Sci. 8, 2056–68 (2015).

Article CAS Google Scholar

Tan, G. et al. Non-equilibrium processing leads to record high thermoelectric figure of merit in PbTe–SrTe. Nat. Commun. 7, 12167 (2016).

Article ADS CAS PubMed PubMed Central Google Scholar

Wu, Y. et al. Lattice strain advances thermoelectrics. Joule 3, 1276–1288 (2019).

Article CAS Google Scholar

Zhao, L. D. et al. High performance thermoelectrics from earth-abundant materials: Enhanced figure of merit in PbS by second phase nanostructures. J. Am. Chem. Soc. 133, 20476–87 (2011).

Article CAS PubMed Google Scholar

Zhao, L. D. et al. Raising the thermoelectric performance of p-type PbS with endotaxial nanostructuring and valence-band offset engineering using cds and zns. J. Am. Chem. Soc. 134, 16327–36 (2012).

Article CAS PubMed Google Scholar

Ibanez, M. et al. High-performance thermoelectric nanocomposites from nanocrystal building blocks. Nat. Commun. 7, 10766 (2016).

Article ADS CAS PubMed PubMed Central Google Scholar

Jiang, B. et al. Realizing high-efficiency power generation in low-cost PbS-based thermoelectric materials. Energy Environ. Sci. 13, 579–591 (2020).

Article CAS Google Scholar

Zhang, Q. et al. High thermoelectric performance by resonant dopant indium in nanostructured SnTe. PNAS 110, 1326–6 (2013).

Google Scholar

Tan, G. et al. High thermoelectric performance of p-type SnTe via a synergistic band engineering and nanostructuring approach. J. Am. Chem. Soc. 136, 7006–17 (2014).

Article CAS PubMed Google Scholar

Tan, G. et al. Codoping in snte: Enhancement of thermoelectric performance through synergy of resonance levels and band convergence. J. Am. Chem. Soc. 137, 5100–12 (2015).

Article CAS PubMed Google Scholar

Li, W. et al. Promoting snte as an eco-friendly solution for p-PbTe thermoelectric via band convergence and interstitial defects. Adv. Mater. 29, 1605887 (2017).

Article CAS Google Scholar

Bhat, D. K. & Shenoy, U. S. Enhanced thermoelectric performance of bulk tin telluride: Synergistic effect of calcium and indium co-doping. Mater. Today Phys. 4, 12–18 (2018).

Article Google Scholar

Tang, J. et al. Manipulation of band structure and interstitial defects for improving thermoelectric SnTe. Adv. Funct. Mater. 28, 1803586 (2018).

Article CAS Google Scholar

Hong, M. et al. Nanoscale pores plus precipitates rendering high-performance thermoelectric SnTe Se with refined band structures. Nano Energy 60, 1–7 (2019).

Article CAS Google Scholar

Hwang, J. et al. Gigantic phonon-scattering cross section to enhance thermoelectric performance in bulk crystals. ACS Nano 13, 8347–8355 (2019).

Article CAS PubMed Google Scholar

Poudel, B. et al. High-thermoelectric performance of nanostructured bismuth antimony telluride bulk alloys. Science 320, 634–8 (2008).

Article ADS CAS PubMed Google Scholar

Xie, W. J., Tang, X. F., Yan, Y. G., Zhang, Q. J. & Tritt, T. M. Unique nanostructures and enhanced thermoelectric performance of melt-spun BiSbTe alloys. Appl. Phys. Lett. 94, 102111 (2009).

Article ADS CAS Google Scholar

Fan, S. et al. p-type Bi0.4Sb1.6Te3 nanocomposites with enhanced figure of merit. Appl. Phys. Lett. 96, 182104 (2010).

Article ADS CAS Google Scholar

Liu, C.-J., Lai, H.-C., Liu, Y.-L. & Chen, L.-R. High thermoelectric figure-of-merit in p-type nanostructured (Bi, Sb) Te fabricated via hydrothermal synthesis and evacuated-and-encapsulated sintering. J. Mater. Chem. 22, 4825–4831 (2012).

Article CAS Google Scholar

Kim, S. I. et al. Dense dislocation arrays embedded in grain boundaries for high-performance bulk thermoelectrics. Science 348, 109–114 (2015).

Article ADS CAS PubMed Google Scholar

Katsuyama, S., Okada, H. & Miyajima, K. Thermoelectric properties of CeFe\(_3\)CoSb\(_{12}\)-MoO\(_2\) composit. Mater. Trans. 49, 1731–6 (2008).

Article CAS Google Scholar

Shi, X. et al. Multiple-filled skutterudites: High thermoelectric figure of merit through separately optimizing electrical and thermal transports. J. Am. Chem. Soc. 133, 7837–7846 (2011).

Article CAS PubMed Google Scholar

Rogl, G. et al. In-doped multifilled n-type skutterudites with ZT = 1.8. Acta Mater. 95, 201–211 (2015).

Article ADS CAS Google Scholar

Zhao, W. et al. Magnetoelectric interaction and transport behaviours in magnetic nanocomposite thermoelectric materials. Nat. Nanotechnol 12, 55–60 (2017).

Article ADS CAS PubMed Google Scholar

Zhao, W. et al. Superparamagnetic enhancement of thermoelectric performance. Nature 549, 247–251 (2017).

Article ADS CAS PubMed Google Scholar

Zhao, L. D. et al. Bi\(_{1-x}\)Sr\(_x\)CuSeO oxyselenides as promising thermoelectric materials. Appl. Phys. Lett. 97, 092118 (2010).

Article ADS CAS Google Scholar

Sui, J. et al. Texturation boosts the thermoelectric performance of BiCuSeO oxyselenides. Energy Environ. Sci. 6, 2916–20 (2013).

Article CAS Google Scholar

Gu, Y. et al. Rational electronic and structural designs advance BiCuSeO thermoelectrics. Adv. Funct. Mater. 31, 2101289 (2021).

Article CAS Google Scholar

Joshi, G. et al. Enhanced thermoelectric figure-of-merit in nanostructured p-type silicon germanium bulk alloys. Nano Lett 8, 4670–4 (2008).

Article ADS CAS PubMed Google Scholar

Zebarjadi, M. et al. Power factor enhancement by modulation doping in bulk nanocomposites. Nano Lett 11, 2225–30 (2011).

Article ADS CAS PubMed Google Scholar

Yu, B. et al. Enhancement of thermoelectric properties by modulation-doping in silicon germanium alloy nanocomposites. Nano Lett 12, 2077–82 (2012).

Article ADS CAS PubMed Google Scholar

Bathula, S. et al. Enhanced thermoelectric figure-of-merit in spark plasma sintered nanostructured n-type SiGe alloys. Appl. Phys. Lett. 101, 213902 (2012).

Article ADS CAS Google Scholar

Culp, S. R. et al. (Zr, Hf)Co(Sb, Sn) half-heusler phases as high-temperature (>700\(^0\,\)C) p-type thermoelectric materials. Appl. Phys. Lett. 93, 022105 (2008).

Article ADS CAS Google Scholar

Yan, X. A. et al. Enhanced thermoelectric figure of merit of p-type half heuslers. Nano Lett. 11, 556–60 (2011).

Article ADS CAS PubMed Google Scholar

Joshi, G. et al. NbFeSb-based p-type half-heuslers for power generation applications. Energy Environ. Sci. 7, 4070–6 (2014).

Article CAS Google Scholar

Fu, C. G. et al. Realizing high figure of merit in heavy-band p-type half-heusler thermoelectric materials. Nat. Commun. 6, 8144 (2015).

Article ADS PubMed Google Scholar

Zhu, H. T. et al. Discovery of TaFeSb-based half-heuslers with high thermoelectric performance. Nat. Commun. 10, 270 (2019).

Article ADS PubMed PubMed Central CAS Google Scholar

Angrist, W. Direct Energy Conversion (Allyn and Bacon, Boston, 1977).

Google Scholar

Fateh, H., Baker, C. A., Hall, M. J. & Shi, L. High fidelity finite difference model for exploring multi-parameter thermoelectric generator design space. Appl. Energy 129, 373–83 (2014).

Article CAS Google Scholar

Gehlhoff, P., Justi, E. & Kohler, M. Verfeinerte theorie der elektrothermischen k\(\ddot{a}\)lteerzeugung. Abhandlungen, Braunschweigischen Wissenschaftlischen Gesellschaft 2, 149–64 (1950).

Google Scholar

Sahin, A. & Yilbas, B. The thermoelement as thermoelectric power generator: Effect of leg geometry. Energy Convers. Manag. 65, 26–32 (2013).

Article Google Scholar

Oki, S., Ito, K. & Natsui, S. Numerical optimization of trapezoidal thermoelectric elements. J. Electro Mater. 45, 1358–64 (2016).

Article ADS CAS Google Scholar

Xuan, X. C., Ng, K. C., Yap, C. & Chua, H. T. Optimization of two-stage thermoelectric coolers with two design configurations. Energy Convers. Manag. 43, 2041–52 (2002).

Article Google Scholar

Ali, H., Sahin, A. Z. & Yilbas, B. S. Thermodynamic analysis of a thermoelectric power generator in relation to geometric configuration device pins. Energy Convers. Manag. 78, 634 (2014).

Article Google Scholar

Yilbas, B. S. & Ali, H. Thermoelectric generator performance analysis: Influence of pin tapering on the first and second law efficiencies. Energy Convers. Manag. 100, 138–46 (2015).

Article Google Scholar

Shi, Y. et al. Nominal power density analysis of thermoelectric pins with non-constant cross sections. Energy Convers. Manag. 97, 1–6 (2015).

Article Google Scholar

Fabi\(\acute{a}\)n-Mijangos, A., Gao, M. & Alvarez-Quintana, J. Enhanced performance thermoelectric module having asymmetrical legs. Energy Convers Manag 148, 1372–81 (2017).

Sisik, B. & LeBlanc, S. The influence of leg shape on thermoelectric performance under constant temperature and heat flux boundary conditions. Front Mater. 7, 595955 (2020).

Article Google Scholar

Jin, Z. H., Wallace, T. T., Lad, R. J. & Su, J. Energy conversion efficiency of an exponentially graded thermoelectric material. J. Electron. Mater. 43, 308–13 (2014).

Article ADS CAS Google Scholar

Koizumi, M. FGM activities in Japan. Compos. Part B Eng. 28, 1–4 (1997).

Article Google Scholar

Kaliazin, A. E., Kuznetsov, V. L. & Rowe, D. M. Rigorous calculations related to functionally graded and segmented thermoelements. In Proceedings of the Twentieth International Conference on Thermoelectrics 286–292 (2001).

Dashevsky, Z., Gelbstein, Y., Edry, I., Drabkin, I. & Dariel, M. P. Optimization of thermoelectric efficiency in graded materials. In Proceedings ICT’03. 22nd International Conference on Thermoelectrics 421–424 (2003).

Kuznetsov, V. L., Kuznetsova, L. A., Kaliazin, A. E. & Rowe, D. M. High performance functionally graded and segmented Bi\(_2\)Te\(_3\) based materials for thermoelectric power generation. J. Mater. Sci. 37, 2893–7 (2002).

Article ADS CAS Google Scholar

Hedegaard, E. M. J., Mamakhel, A. A. H., Reardon, H. & Iversen, B. Functionally graded (PbTe)\(_{1-x}\)(SnTe)\(_x\) thermoelectrics. Chem. Mater. 30, 280–7 (2018).

Article CAS Google Scholar

Hedegaard, E. M. J., Johnsen, S., Bjerg, L., Borup, K. A. & Iversen, B. Functionally graded Ge\(_{1-x}\)Si\(_x\) thermoelectrics by simultaneous band gap and carrier density engineering. Chem. Mater. 26, 4992–7 (2014).

Article CAS Google Scholar

Dariel, M. P., Dashevsky, Z., Jarashnely, A., Shusterman, S. & Horowitz, A. Carrier concentration gradient generated in \(p\) type PbTe. J. Cryst. Growth 234, 164 (2002).

Article ADS CAS Google Scholar

Dashevsky, Z., Shusterman, S., Dariel, M. P. & Drabkin, I. Thermoelectric efficiency in graded indium-doped PbTe crystals. J. Appl. Phys. 92, 1425–30 (2002).

Article ADS CAS Google Scholar

Gelbstein, Y., Dashevsky, Z. & Dariel, M. P. High performance \(n\)-type PbTe-based materials for thermoelectric applications. Phys. B Condens. Matter. 363, 196–205 (2005).

Article ADS CAS Google Scholar

Ioffe, A. F. Thermoelectric battery, SU Author’s Certificate No. 126158, claimed 19/08/1949. Invention Bulletin No. 4 (1960).

Schilz, J., Helmers, L., Kang, Y. S., Noda, Y. & Niino, M. Bismuth-telluride/iron disilicide segmented thermoelectric elements: Pattering, preparation and properties. In Proceedings of the 16th International Conference on Thermoelectrics 375–378 (1997).

Noda, Y., Kang, Y. S. & Niino, M. Preparation and characterization of segmented-type thermoelectric branches of Bi\(_2\)Te\(_3\)/PbTe. In Proceedings of the 16th International Conference on Thermoelectrics 371–374 (1997).

Shiota, I. & Nishida, I. Development of FGM thermoelectric materials in Japan - the state of the art. In Proceedings of the 16th International Conference on Thermoelectrics 364–70 (1997).

Helmers, L., M\(\ddot{u}\)ller, E., Schilz, J. & Kaysser, W. A. Graded and stacked thermoelectric generators—numerical description and maximisation of output power. Mater. Sci .Eng. B 56(1), 60–8 (1998).

Miyamoto, Y., Niino, M. & Koizumi, M. FGM research programs in Japan—from structural to functional uses. In Shiota I, Miyamoto Y, editors. Proceedings of the Fourth International Symposium on FGM 1–8 (1996).

El-Genk, M. S. & Saber, H. Efficient segmented thermoelectric unicouples for space power applications. Energy Convers. Manag. 44, 1069 (2003).

Article CAS Google Scholar

Snyder, G. J. & Ursell, T. S. Thermoelectric efficiency and compatibility. Phys. Rev. Lett. 91, 148301 (2003).

Article ADS PubMed CAS Google Scholar

Seifert, W., M\(\ddot{u}\)ller, E., Snyder, G. J. & Walczak, S. Compatibility factor for the power output of a thermogenerator. Phys Status Solidi (RRL) Rapid Res Lett 1, 250–2 (2007).

Neophytou, N. et al. Simultaneous increase in electrical conductivity and seebeck coefficient in highly boron-doped nanocrystalline Si. Nanotechnology 24, 205402 (2013).

Article ADS PubMed CAS Google Scholar

Snyder, G. J. & Toberer, E. S. Complex thermoelectric materials. Nat. Mater. 7, 105–14 (2008).

Article ADS CAS PubMed Google Scholar

Shur, M. Semiconductor Devices, Physics and Technology (Prentice Hall, Upper Saddle River, 1990).

Google Scholar

Sze, S. Semiconductor Devices, Physics and Technology 2nd edn. (Wiley, New York, 2002).

Google Scholar

Sze, S. & Ng, K. Physics of Semiconductor Devices 3rd edn. (Wiley, Hoboken, 2007).

Google Scholar

Colinge, J. P. & Colinge, C. A. Physics of Semiconductor Devices (Kluwer Academic Publishers, New York, 1997).

Google Scholar

Ioffe, A. Semiconductor Thermoelements and Thermoelectric Cooling (Infosearch, ltd (London UK, 1957).

Goupil, C. Thermodynamics of thermoelectricity. In Thermodynamics (ed. Tadashi, M.) 275–92 (Intech, Shanghai, 2011).

Download references

Formerly, University of Alaska Fairbanks, 306 Tanana Drive, Duckering Building, Fairbanks, AK, 99775, USA

Tinggang Zhang

You can also search for this author in PubMed Google Scholar

T.G. conceptualized the current work, developed the methodology, and conducted the analysis, drafted, reviewed and edited the manuscript.

Correspondence to Tinggang Zhang.

The author declares no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Cross-sectional area of an n-type leg.

Cross-sectional area of a p-type leg.

Diffusivity of electron in the base.

Diffusivity of hole in the collector.

Diffusivity of hole in the emitter.

Base electric current.

Collector electric current.

Emitter electric current.

Electric current of a thermoelectric module.

Electron diffusion current in the base at the base-collector junction.

Electron diffusion current in the base at the base-emitter junction.

Hole diffusion current in the collector at the base-collector junction.

Hole diffusion current in the emitter at the base-emitter junction.

Total electric current of a thermoelectric module.

Electric current.

Thermal conductance of a thermoelectric module.

Thermal conductance.

Length of an n-type leg.

Length of a p-type leg.

Diffusion length of electron in the base.

Diffusion length of holes in the collector.

Diffusion length of holes in the emitter.

Thermoelectric power factor.

Output power of the original thermoelectric module.

Output power of the new thermoelectric module.

Heat flux at cold side of a thermoelectric module.

Heat flux at hot side of a thermoelectric module.

Intrinsic electrical resistance.

Electrical load resistance.

Base circuit load resistance.

Supply circuit load resistance.

Circuit load resistance for n-p-n transistor.

Circuit load resistance for p-n-p transistor.

Electrical resistance.

Cross-sectional area of the base region.

Absolute temperature at cold side.

Absolute temperature at hot side.

Base-emitter bias.

Supplied positive voltage.

Collector-emitter bias.

Input signal voltage.

Output signal voltage equals Vce.

kBT/q—thermal voltage.

Total terminal voltage of a thermoelectric module.

Width of the base.

Dimensionless figure of merit of thermocouple.

Absolute temperature difference.

Variation of input signal voltage.

Variation of output signal voltage.

Excess electron concentration in the base.

Incremental carrier concentration.

Excess hole concentration in the collector.

Excess hole concentration in the emitter.

Fraction of the collector current.

Fraction of the collector-emitter bias.

Carrier concentration.

Equilibrium concentration of minority carrier in the base.

Electron concentration in the base at the base-collector junction.

Electron concentration in the base at base-emitter junction.

Electron concentration of the base.

Intrinsic carrier concentration of the base.

Hole concentration at the far-end of the collector (relative to the base-collector junction).

Hole concentration in the collector.

Hole concentration at the far-end of the emitter (relative to the base-emitter junction).

Hole concentration in the emitter.

Intrinsic hole concentration in the collector.

Intrinsic hole concentration in the emitter.

Electric charge.

Length of the collector region.

Length of the emitter region.

Dimensionless figure of merit of thermoelectric material.

Bipolar junction transistor.

Functional graded thermoelectric material.

Metal-oide-semiconductor field efffect transistor.

Phonon glass electron crystal.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

Zhang, T. Optimizing thermocouple’s ZT through design innovation. Sci Rep 11, 19338 (2021). https://doi.org/10.1038/s41598-021-98562-z

Download citation

Received: 16 June 2021

Accepted: 08 September 2021

Published: 29 September 2021

DOI: https://doi.org/10.1038/s41598-021-98562-z

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

Scientific Reports (2022)