Synergetic Approach for Superior Thermoelectric Performance in PbTe-PbSe-PbS Quaternary Alloys and Composites

: Thermoelectric power generation is an energy conversion technology from heat to electric energy, which can be applied to waste heat power conversion. Among thermoelectric materials (TE), PbTe-PbSe-PbS quaternary alloys and composites are promising candidates for thermoelectric power generation applications in the mid-temperature operating range from 500 to ~850 K. Besides, the thermoelectric performance of quaternary alloys and composites is not fully optimized regarding its composition and synthesis process. In the quaternary system, PbTe-PbSe-PbS, it was found that PbS will form nanoprecipitation in the matrix of quaternary alloy for a small content of PbS (≤0.07), which reduces the lattice thermal conductivity. The power factor of PbTe-PbSe-PbS quaternary alloys can be significantly enhanced by using a band convergence in PbTe 1−x Se x . The band structure modifications, with the result of simultaneous PbS nanoprecipitation, give rise to a high 𝑍𝑇 value of 2.3 at 800 K for (PbTe) 0.95−x (PbSe) x (PbS) 0.05 . The chemical potential tuning by effective K-doping ( 𝑥 = 0.02) and PbS substitution reveals a high power factor and low thermal conductivity, resulting in a comparatively high 𝑍𝑇 value of 1.72 at 800 K. The combination of a high Seebeck coefficient and low thermal conductivity results in a very high 𝑍𝑇 value of 1.52 at 700 K as n-type materials for low Cl-doped ( 𝑥 = 0.0005) (PbTe 0.93−x Se 0.07 Cl x ) 0.93 (PbS) 0.07 composites. Therefore, this review presents the simultaneous emergence of effective chemical potential tuning, band convergence, and nanoprecipitation, giving rise to a significant enhancement of the thermoelectric performance of both 𝑝 - and 𝑛 -type PbTe-PbSe-PbS quaternary alloy and composite TE materials. n-type (PbTe 0.93−x Se 0.07 Cl x ) 0.93 (PbS) 0.07 . The results showed that the randomly distributed interface potential, driven by the random distribution of nanoprecipitates, is an effective way to


Introduction
Recently, much effort has been devoted to searching for materials that can convert waste heat into electricity owing to the growing global demand for green energy. Thermoelectric power generation is a direct heat to electric energy conversion technology and can be applied to waste heat power conversion as well. Effective thermoelectric materials are expected to have low thermal conductivity,  , and electrical resistivity, ρ , and a high Seebeck coefficient, , which is defined by the dimensionless thermoelectric figure-of-merit, = 2 /( ). The trade-off relationship among the parameters, , ρ, and  , makes it difficult to have a high value. During the decade, the development of thermoelectric materials has progressed significantly in terms of nanostructure and new material design. Employing a nanostructure reduces the lattice thermal conductivity by scattering the acoustic phonons [1]. It has been suggested that the low-dimensional nanostructure can give rise to the quantum confinement effect, resulting in a high power factor [1], but the experimental realization in bulk compounds is elusive. In addition, the phonon-glass-electron-crystal (PGEC) concept has long been sought out for the scattering of phonons while preserving electrical transport [2].
Among thermoelectric materials, PbTe is one of the most promising and efficient thermoelectric material for mid-temperature range applications. The thermoelectric performance of PbTe-based materials has been improved significantly through electronic structure modifications [3][4][5][6]. Semiconductor PbTe of p-type can be produced by doping with Tl, Na, and K on Pb sites [5,[7][8][9][10]. On the other hand, n-type PbTe can be produced by Cl doping on the Te site [11]. In addition, its low content in the Earth's crust and the high cost of tellurium inhibit the usage of PbTe in large-scale production. A lot of efforts have been made to replace Te by Se, and S to form PbSe [12,13], PbS [14], or the ternary alloys, PbTe-PbSe [3,15,16] and PbSe-PbS [17,18]. In the PbTe-PbSe system, high thermoelectric performance is mainly attributed to band structure engineering. Band structure engineering in PbTe1−xSex is beneficial to gain a high Seebeck coefficient while maintaining high electrical conductivity by doping [3]. It is believed that a low lattice thermal conductivity, L  , in ternary PbTe1-xSex alloys is explained mainly by the formation of point defects, created by Te/Se's mixed occupation in the rock salt structure [3]. In addition, the thermoelectric performance in PbTe1−xSex ternary alloys can be increased by nanostructuring, resulting in low lattice thermal conductivity. Nanostructuring in the PbTe-PbS system can be achieved using bulk phase separation either by nucleation or spinodal decomposition depending on the relative phase fraction [19].
Recent advances report that the lattice thermal conductivity in the quaternary system of (PbTe)1−x−y(PbSe)x(PbSe)y can be reduced by point defect, which is produced by triple disorder in the rock salt structure [20][21][22]. High ≈ 2.2 at 800 K was obtained in p-type (PbTe)1−2x(PbSe)x(PbS)x quaternary alloys due to band engineering and phonon scattering from point defects [20,23,24]. Quaternary alloys PbTe-PbSe-PbS manifest themselves as effective n-type materials as well. For example, n-type (PbTe)0.75(PbSe)0.1(PbS)0.15 exhibited a high value of 1.1 at 800 K [22]. Therefore, the quaternary alloy system, PbTe-PbSe-PbS, is a promising candidate for use in thermometric power generation devices because alloy materials exhibit both n-and p-type properties with high thermoelectric performance [25]. The PbS precipitates in the matrix, resulting in reduction of the lattice thermal conductivity. The size, shape, and distribution of the nanoparticles contribute to the thermoelectric properties [26]. Here, this review briefly introduces the simultaneous emergence of band convergence, nanostructuring, and chemical potential tuning in PbTe-PbSe-PbS quaternary alloys and composites, which provides an extremely high thermoelectric figure-of-merit in p-and n-type materials.

Nanostructure in PbTe-PbSe-PbS Quaternary Alloys
In spite of high values for Pb-based quaternary alloys, the thermoelectric performances of reported alloys are partially optimized in terms of the content of binary compounds and synthesis process [27][28][29]. In the systematic investigation of (PbTe)0.95−x(PbSe)x(PbS)0.05 ( = 0.0, 0.05, 0.10, 0.15, 0.20, 0.35, and 0.95) quaternary alloys with 1 at. % Na-doping, PbS content is fixed as 5 at. %, which is enough for nanoprecipitation while maintaining nanoprecipitation and without changing the electronic band structure. The electronic band structure modification was associated with content of PbSe in (PbTe)0.95−x(PbSe)x(PbS)0.05. Figure 1 shows the X-ray diffraction (XRD) pattern of (PbTe)0.95−x(PbSe)x(PbS)0.05 having a single phase with a cubic structure. The lattice parameters are decreased with increasing Se content, which follows the Vegard's law, indicating a solid solution of PbSe and PbTe, as shown in Figure 1b. The decrease of the lattice parameter with increasing Se content is due to the fact that the atomic radius of Se is smaller compared to Te.   The lattice parameters of PbTe, PbSe, and PbS are 3.22, 3.07, and 2.965 Å , respectively, along the  (200) plane [31,32]. Numerous nanoprecipitates are also found in K-doped (Pb0.98K0.02Te)0.7(PbSe)0.25(PbS)0.05 [7] and Cl-doped (PbTe0.93−xSe0.07Clx)0.93(PbS)0.07 ( = 0.0005) [11], which are synthesized by same process with (PbTe)0. 75 Figures 3d and 4f, respectively. IFFT images along the plane show that there are many line dislocations at coherent interfaces between nanoprecipitates within the PbTe matrix. In order to identify the lattice dislocation more clearly, geometric phase analysis (GPA) of TEM images is reported. The strain field map profile around the dislocation area shows clearly the lattice strain due to lattice parameter misfit between the matrix and the precipitates. It indicates that the lattice dislocation may generate a strain field between the matrix and nanoprecipitates, as shown in Figure 4f. A positive high strain filed is observed at the edge dislocation, where there are interfaces between the matrix and peripheral precipitates.  Several key questions arise when considering the nanoprecipitates in complex PbTe-PbSe-PbS quaternary alloys. The substitution of multiple elements in the PbTe matrix with the rapid quenching process results in lower free energy for the creation of nanoprecipitates in the matrix of PbTe [8]. A similar phenomenon is also found in the ternary system, PbTe-SrTe. Nanoprecipitates are formed in PbTe-SrTe due to the difference in lattice parameters between PbTe and SrTe. The limited SrTe doping is segregated in the PbTe matrix, resulting in nanoprecipitations [32]. In PbTe-PbSe-PbS quaternary alloys, the miscibility gap between the PbTe and PbS phase is expected in the thermodynamic phase diagram, implying phase separation via metastable nucleation during grain growth or spinodal decomposition of nanoprecipitates [33] because the thermodynamic nanoscale precipitation of PbS in the PbTe-PbSe matrix is likely [7,31,34].
In order to identify the phase of the nanoprecipitates, elemental analysis is reported by energy dispersive X-ray spectroscopy (EDS) in (Pb0.98K0.02Te)0.7(PbSe)0.25(PbS)0.05 [7], as shown in Figure 5. Figures 5a and 5b show the nanoprecipitates and phase boundary due to spinodal decomposition, displayed as the dark region. The form of nanoprecipitates of (Pb0.98K0.02Te)0.7(PbSe)0.25(PbS)0.05 exhibits two different morphologies: (1) Nanodot and (2) stripe phase due to spinodal decomposition, as shown in Figure 5b. The nanoprecipitates (region 1 and region 3) are a PbS-rich phase within the PbTe matrix (region 1), as shown in Figure 5d, from the EDS spectrum.

Simultaneous Occurrence of Nanostructuring and Band Convergence in p-Type
Earlier reports show high thermoelectric performance in p-type (PbTe)0.95−x(PbSe)x(PbS)0.05 alloys by adopting simultaneous occurrence of band structure engineering and nanostructuring due to nanoprecipitation [8]. Band convergence is effective in increasing the power factor while maximizing the power factor in both high electrical conductivity and the Seebeck coefficient by highly dispersive and degenerated energy bands, respectively [35]. According to theoretical band structure calculation of PbTe, PbSe, and PbS, those compounds are known to have two valence bands, which are light hole, , and heavy hole, Σ , bands. Experimental realization of band convergence in PbTe-PbSe alloy enhances the power factor by high valley degeneracy with dispersive energy bands [36]. The energy difference between the conduction band minimum and valence, , band maximum, Δ − , and  band maximum, Δ −Σ , changes with the temperature, , and content, , as follows [37]: 0.36 0.10 The total Seebeck coefficient is defined by coefficients and electric conductivity provided by two valence bands as follows [8,36]: where , Σ , , and Σ are the Seebeck coefficients and electric conductivity due to the light hole, , and heavy hole, Σ, valence bands, respectively. The light hole, , band has a dispersive energy band, resulting in the high mobility of charge carriers, while the heavy hole, Σ, band is a flat band with a heavy effective mass of charge carriers. When we align the light, , and heavy, Σ, hole bands, highly dispersive and degenerated energy bands give rise to a high Seebeck coefficient while maintaining the high mobility of charge carriers.
Previous reports demonstrate that the effect of band alignment of and Σ bands in 1% Na-doped (PbTe)0.95−x(PbSe)x(PbS)0.05 gives rise to the increase in the value of the Seebeck coefficient and power factor. Figure 6 shows the temperature-dependences of the Seebeck coefficient and electrical conductivity of (PbTe)0.95−x(PbSe)x(PbS)0.05. S(T) of the compounds increases with increasing the temperature, indicating metallic or degenerated semiconducting behavior, as presented in Figure 6a. Figure 6a shows a broad shoulder in S(T) near 700 K for alloys with < 0.2. The plateau around 700 K is caused by band alignment of and Σ bands in (PbTe)0.95−x(PbSe)x(PbS)0.05 due to hole transfer from the to the Σ band [8]. However, the plateau disappears in the alloys (PbTe)0.95−x(PbSe)x(PbS)0.05 with ≥ 0.2 due to a increase in the band gap with growth of the Se content, implying that hole excitation from the to the Σ band may occur at temperatures higher than 800 K [8]. Electrical resistivity, ( ), shows behavior similar to the Seebeck coefficient, as shown in Figure 6b, where broad shoulders of ( ) are observed near 700 K for alloys with < 0.2 [8]. Figure 6c shows the temperature-dependent power factor for series of ternary alloys (PbTe)0.95−x(PbSe)x(PbS)0.05. The maximum power factor of 28.7 µ W m −1 K −2 at 800 K is achieved in alloy samples with = 0.15, which is higher than those of pristine PbTe and PbTe-based alloys [31,[36][37][38][39][40][41][42]. So, a high power factor is due to the contribution of band convergence of and Σ bands. The lowest power factor is observed in alloy samples with = 0.35 because an increase in the Se content results in a decrease in the mobility of charge carriers due to the scattering of charge carriers by impurities, which is not beneficial to the increase in the power factor [36]. In order to confirm the contribution of and Σ bands' convergence to thermoelectric properties, the Pisarenko plot is shown based on two valence bands' model compared with PbTe [20,43,44], Pb1−xKxTe [10], and Na-doped PbTe [43], as shown in Figure 6d. Instead of the single parabolic band model, the Pisarenko plot corresponds to the two valence bands' model with effective masses of light, , and heavy, Σ, bands equal to 0.36 and 1.6 , respectively, which means that the Fermi levels lie deep within the valence band and the two valance bands contribute significantly to the Seebeck coefficient [8].  The lowest ( ) is obtained at = 0.95. The significant decrease in ( ) due to Se incorporation is reported to arise from the phonon scattering caused by alloying. Figure 7b shows the lattice thermal conductivity, ( ), of alloys having the lowest ( )  0.40 W m −1 K −1 at 800 K that is measured in the alloy with = 0.20. Values of ( ) in the alloy with = 0.20 are lower on 60% at 300 K and 55% at 800 K compared to the pristine PbTe compound, and at 800 K, this approaches the glass limit for the bulk PbTe compound [45].
In order to clarify the mechanism of thermal conductivity in(PbTe)0.95−x(PbSe)x(PbS)0.05 alloys, the investigation is based on the Debye-Callaway model using the expression [32,46]: where , ℏ , , , θ, and are the Boltzmann constant, reduced Plank constant, absolute temperature, sound velocity, Debye temperature, and = ℏ / , respectively. and are the relaxation times due to normal phonon-phonon scattering and the combined relaxation time. From Matthiessen's rule, the combined relaxation time is obtained by accounting for the relaxation times from various scattering processes. In this case, the combined relaxation time given as: where , , , , , and are the relaxation times corresponding to scattering from the Umklapp process, normal process, boundaries, strain, dislocations, and precipitates [14,[47][48][49][50][51]. In the previous reports, the authors consider four different scattering mechanisms, such as point defect scattering, phonon-phonon scattering, grain boundary scattering, and nanoprecipitates scattering with the radius, , and volume fraction, . Therefore, the phonon relaxation time can be expressed as: 42 13 2 where Lgr is the grain's average size, and coefficients and are constant prefactors. The value of depends only on the crystal structure. Thus, the authors obtain the value of by fitting the above equation to the lattice thermal conductivity of undoped (PbTe)0.95−x(PbSe)x(PbS)0.05 alloy. The grain boundaries' scattering is not distinguished, as it is convenient to define an effective mean free path [14]: Based on TEM observation, the average size of nanoparticles is about 10 nm as shown in Figure  2 and other appropriate parameters of the Debye temperature and sound velocity are obtained from [14]. The calculated theoretical, ( ), is fitted with the experimental data.  Low thermal conductivity due to nanoprecipitates and a high power factor caused by band engineering directly affects the thermoelectric performance of (PbTe)0.95−x(PbSe)x(PbS)0.05 alloys, which is displayed clearly in the dimensionless figure-of-merit, . Figure 9a shows the temperature dependences of for (PbTe)0.95−x(PbSe)x(PbS)0.05 alloys. The alloy with = 0.20 shows the highest value 2.3 at 800 K. Biswas et al. [32] reported high value of 2.0 at 800 K in hierarchical architecture of the atomic and mesoscale in PbTe alloyed with 2% SrTe. In order to achieve a hierarchical architecture, additional processing is needed to produce fine powder and a mesoscale microstructure. On the other hand, the direct synthesis of (PbTe)0.95−x(PbSe)x(PbS)0.05 alloys realizes high without a multiple sintering process. For practical applications, we also estimated the average ZTavg value as defined by [52]: is the figure-of-merit, Tc is the cold side temperature, Th is the hot side temperature, Tavg is the average temperature (Th + Tc)/2, and ΔT is the temperature difference between hot and cold sides ( ℎ − ). The calculated ZTavg values of PbTe, (PbTe)0.75(PbS)0.25 [34], (PbTe)0.88(PbS)0.12 [36], and (PbTe)0.84(PbSe)0.07(PbS)0.07 doped with 2% Na [20] is shown in Figure 9b. It shows the highest average ZTavg of the (PbTe)0.75(PbSe)0.2(PbS)0.05 alloy compared to the other compounds quoted in the previous reports [20,34,36].

Chemical Potential Tuning and Nanostructuring in p-
Ahmad et al. suggested that the power factor of the PbTe compound can be enhanced by resonant-like density of state (DOS) distortions in p-type Pb1−xAxTe ( = K, Rb, and Cs but not Na) [53]. However, Androulakis et al. showed that K doping does not form resonance states but can control the energy difference of the maxima of two primary valence sub-bands in PbTe [54]. The role of K in PbTe has been extensively analyzed and values of 1.3 at 673 K [10] were obtained, which is comparable with that of Na-doped PbTe at the same temperature [43]. Furthermore, a high thermoelectric performance ( ≈ 1.6 at 773 K) was obtained in K-doped PbTe1−xSex ternary alloys due to an increase in the DOS around the Fermi levels, resulting in a higher Seebeck coefficient for the two valence bands of PbTe1−xSex [10]. Given the valence bands' convergence in PbTe-PbSe [10] and the presence of a nanostructure in PbTe-PbS [33,34], these are effective for increasing the power factor and decreasing thermal conductivity, respectively. Therefore, the role of K is examined in enhancing of quaternary alloys PbTe-PbSe-PbS with the combination of band converge and nanostructure.
The room temperature Hall effect measurements show the p-type conductivity of Pb1−xKxTe0.7Se0.25S0.05 alloys. Figure 10 shows the Hall concentration, nH, and Hall mobility, μH, of charge carriers at room temperature in Pb1−xKxTe0.7Se0.25S0.05 alloys depending on the K-doping level. The Hall concentration of charge carriers is reduced until = 0.02 and then increased significantly at higher doping levels > 0.02 from 4.33 × 10 19 cm −3 ( = 0.02) to 7.08 × 10 19 cm −3 ( = 0.03), which indicates that K doping effectively tunes charge carriers' concentration in Pb1−xKxTe0.7Se0.25S0.05 alloys. Hall mobility is monotonically decreased with an increasing K-doping level from 350.2 cm 2 ( = 0.01) to 120 cm 2 V −1 s −1 ( = 0.03). The reduction of mobility can be understood in terms of the defect and alloy scattering due to K doping.  Figure 11a shows the temperature-dependent Seebeck coefficient, S(T), of Pb1−xKxTe0.7Se0.25S0.05 and Pb0.98K0.02Te alloys. S(T)increased with temperature and showed a broad maximum near 700 K. The maximum Seebeck coefficient reached 313 μV/K at 750 K for doping = 0.02, which is much higher even than that previously reported for Na-doped PbTe (≈ 260 μV/K at 775 K), which may be associated with the high effective mass, m * . Here, the estimated effective masses in terms of a single parabolic band model are [55]: where F1/2(η) is the Fermi integral, η is the reduced Fermi energy,  1 is the Hall factor, and kB is the Boltzmann constant. The effective masses of charge carriers in Pb1−xKxTe0.7Se0.25S0.05 alloys are presented in Table 1. The effective masses increased with an increase in the K-doping level over ≥ 0.015. The effective masses of charge carriers in alloys are higher than in the parent compound, PbTe m * ~ (0.25-0.27)me [21]. The larger effective masses are related to two valence bands [3,10]. Because and Σ bands come closer to each other with increasing temperature, holes are transferred from the to the Σ band, which is the cause of the broad maximum near 700 K.
The electrical resistivity, ρ(T), of Pb1−xKxTe0.7Se0.25S0.05 and Pb0.98K0.02Te alloys increases with temperature and then reaches a broad maximum, as shown in Figure 11b; variation of ρ(T) with the K-doping level is attributed to the changing concentration of the charge carriers (Table 1). Another reason is impurity or alloy scattering due to K-doping, which is confirmed by the decrease in Hall mobility ( Table 1). The power factor of Pb1−xKxTe0.7Se0.25S0.05 and Pb0.98K0.02Te alloys displays the maximum value equal to 27.78 µ W m −1 K −2 at 600 K when = 0.01, which is attributed to the significant reduction in electrical resistivity. The Pisarenko plot made on two valence bands' model for Pb1−xKxTe0.7Se0.25S0.05 alloys, Na-doped PbTe, and Pb1-xKxTe is shown in Figure 11d [10,20,43,44] and the values of effective masses of the heavy hole ℎ * = 1.2 me (Σ band) and light hole * = 0.36 me ( band) with an energy band gap ΔE ≈ 0.12 eV is shown. The experimental data points of Pb1−xKxTe0.7Se0.25S0.05 alloys lie below the theoretical curve plotted using the two valence bands' model, indicating a lower effective mass of the charge carriers in the alloys. Based on the Pisaranko plot made on the two valence bands' model ( Figure 11d), it is clear that the second valence band contributes to the Seebeck coefficient and there is no effect of the resonance level. The temperature-dependent thermal conductivity, κ(T), of Pb1−xKxTe0.7Se0.25S0.05 alloys decreases with an increase in temperature as shown in Figure 12a. κ(T)decreased with an increase of the K-doping level until ≤ 0.02 and then increased again for > 0.02 as a result of decreased electric resistivity. The thermal conductivity, κ(T), for the alloy with ≤ 0.02 showed a weak increase at the high temperature region (T ≥ 700 K), implying bipolar diffusion of charge carriers. This small upturn of κ(T) was also found in 3% Na-doped (PbTe)1−x(PbS)x alloy [31]. The lowest thermal conductivity reached 1.96 W m −1 K −1 at room temperature and 0.95 W m −1 K −1 at the high temperature for (PbTe)1−x(PbS)x ( = 0.02). The extracted lattice thermal conductivity, , was found by subtracting the electronic thermal conductivity, , determined by the Wiedemann-Franz law as = 0 / , where 0 is the Lorenz number, is the absolute temperature, and ρ is the electrical resistivity. Figure 12b shows very low in Pb1−xKxTe0.7Se0.25S0.05 alloys of 0.91 W m −1 K −1 at room temperature and 0.69 W m −1 K −1 at high temperature.
In a composite of two semiconductors with different energy band gaps or metal/semiconductor composites, a band bending effect exists at the interfaces. The band bending (energy barrier) induced by different Fermi levels of two materials can selectively scatter charge carriers due to the energy-dependent scattering time, resulting in enhancement of the Seebeck coefficient. In addition, nanoinclusions in a matrix will scatter effectively on phonons. Therefore, the reports demonstrate the effects of electron and phonon scattering by nanoinclusions in n-type (PbTe0.93-xSe0.07Clx)0.93(PbS)0.07 ( = 0.0005, 0.01, 0.1, and 0.2) composites. In the matrix, nanoinclusions provide an increase in the Seebeck coefficient while reducing the lattice thermal conductivity in (PbTe0.93−xSe0.07Clx)0.93(PbS)0.07. Doping by Cl lifts the Fermi level towards the bottom of the conduction band, resulting in an increased electron concentration. The simultaneous emergence of a high Seebeck coefficient and low thermal conductivity resulted in a high value of 1.52 at 700 K for a low Cl-doping level ( = 0.0005), which is a very high value in n-type thermoelectric materials. The randomly distributed interface potential induced by Fermi level tuning with nanoinclusions is a new technique for investigating thermoelectric properties [11]. The nanoprecipitates ( Figure 4) show a strong lattice strain near precipitates in (PbTe0.93−xSe0.07Clx)0.93(PbS)0.07 alloys from TEM images, and their electron diffraction and strain field map from IFFT mapping of TEM are also shown. In view of the band gap difference between PbTe matrix and PbS nanoprecipitates the considered chemical potential is tuned to the conduction band side by Cl doping. Figure 16 shows schematic image of inclusion in matrix (left panel) and induced band bending effect on interface between matrix and nanoinclusion (right panel). The PbTe1−xSex matrix has intrinsic point defects and dislocations. Additional PbS nanoprecipitation further decreases the lattice thermal conductivity by scattering phonons. Cl doping aims to lift the chemical potential of the conduction band bottom for n-type materials showing metallic behavior in electrical transport. Because PbS has a sizable energy band gap, there should be a band bending effect near the interface between the matrix and PbS nanoprecipitates. In that case, filtering of charge carriers with energy lower than the appropriate barrier height near the interface is possible [66]. Figure 17a shows the Hall concentration, nH, and Hall mobility, μH, of charge carriers in (PbTe0.93−xSe0.07Clx)0.93(PbS)0.07 ( = 0.0005, 0.01, 0.1, 0.2).
A negative Hall coefficient is indicates n-type conductivity. The electron concentration, nH, increased systematically with an increasing Cl-doping level. The Hall mobility, μH, decreased with increased Cl doping from 954 to lower than 200 cm 2 V −1 s −1 with increasing Cl doping, which is due to enhanced charge carrier scattering. Using a single parabolic band model, Hall mobility can be expressed as [67]: where m * is the effective mass of charge carriers, and Fn(η) is the Fermi integral defined as follows: where η is the reduced electrochemical potential, is the longitudinal elastic constant, is the valley degeneracy, and Edef is the deformation potential.  [68], PbSe1−xClx [69], PbS1−xClx [14], PbTe0.9Se0.1 [63], PbSe0.84Se0.16 [29], and PbTe0.75−xClxSe0.1S0.15 [22]. Reproduced with permission from [11], copyright 2017 the Royal Society of Chemistry. Figure 17b shows the Pisarenko plot of the Hall mobility versus the Hall concentration of charge carriers in terms of the single parabolic band model for Pb1−xIxTe [68], PbSe1-xClx [69], PbS1−xClx [14], PbTe0.9Se0.1 [63], PbSe0.84Se0.16 [29], and PbTe0.75−xClxSe0.1S0.15 [22]. The Pisarenko plot of (PbTe0.93−xSe0.07Clx)0.93(PbS)0.07 lies below the value expected from the single parabolic model. It can be understood by the scattering of charge carriers by microscopic cracks, grain boundaries, disorder, and defect, leading to high residual resistivity and low Hall mobility. Figure 18 shows the temperature-dependent thermoelectric properties of composites (PbTe0.93−xSe0.07Clx)0.93(PbS)0.07 ( = 0.0005, 0.01, 0.1, 0.2). Figure 18a displays the electric conductivity, σ(T), in the temperature range from 300 to 800 K. The monotonic decrease in electric conductivity with temperature indicates typical degenerated semiconducting or metallic behavior in all the concentrations. The electric conductivity increases with an increase in the Cl-doping content due to higher concentration of charge carriers. The temperature-dependent Seebeck coefficients, S(T) (Figure 18b), show negative values, indicating ntype conductivity. Because the Seebeck coefficient follows the Mott relation [70], in which a higher concentration of charge carriers leads to a lower Seebeck coefficient, the behavior of electrical conductivity showed a trade-off relationship with the Seebeck coefficient. The most negative value of the Seebeck coefficient drops to −84.12 μV/K at 300 K for = 0.0005, which corresponds to an electron concentration of 1.62 × 10 19 cm −3 , while the smallest negative value of the Seebeck coefficient is observed as −11.14 μV/K at 300 K for = 0.1, which corresponds to an electron concentration of 1.83 × 10 20 cm −3 . The Seebeck coefficient of the composite with = 0.0005 is higher than those of n-type PbTe, PbSe, and PbS binary compounds over an entire temperature range. The peak Seebeck coefficient equals −236.63 μV/K for = 0.0005, which is higher than PbTe, PbS, and PbSe by as much as 7.57%, 15.77%, and 4.60%, respectively. There is a broad shoulder in the Seebeck coefficient at ≥ 750 K for = 0.0005, which could be mainly ascribed to the thermal excitation of the bipolar diffusion effect. The Seebeck coefficient for bipolar transport effect is given by [71]: where the subscript , ℎ, , and denotes the partial properties of the electron, hole, electron density, and hole density, respectively. The is the mobility ratio defined by = / ℎ . The Seebeck coefficient at high temperatures, ≥ 750 K for = 0.0005, come from the thermal activation of minority charge carriers. Figure 18c shows the temperature-dependent power factor, 2 . The maximum power factor is about 21.72 μW cm −1 K −2 at 500 K for = 0.0005 owing to the high absolute value of the Seebeck coefficient. Even though the room temperature power factor decreased with an increasing Cl-doping content, attributed to the decrease in the Seebeck coefficient, the power factor at high temperatures is not sensitive to the Cl-doping level. Temperature-dependent thermal conductivity decreased when increasing the temperature, which is a typical behavior of acoustic phonon scattering as presented in Figure 18d. The lowest thermal conductivity, κ, is found at = 0.0005 (κ ≈ 0.84 W m −1 K −1 at 750 K). The increase in κ over a wide temperature range with an increase in the Cl-doping level is due to the significant electron contribution to the thermal conductivity.
Here, the literature argues that the interface energy barrier between the matrix and nanoprecipitates gives rise to charge carriers' energy filtering effect, resulting in the enhancement of the Seebeck coefficient. Few authors employed the theoretical Boltzmann transport calculation with relaxation time approximation as the following relations of the electrical conductivity, σ, Seebeck coefficient, , and electron thermal conductivity, κel [72]: and (d) total thermal conductivity, ( ), versus temperature. Reproduced with permission from [11], copyright 2017 the Royal Society of Chemistry.
The measured σ(T) of composite with = 0.0005 and theoretical calculation values are lower than those of PbTe over the entire temperature range. Such a reduction in σ(T) is caused by a decrease in the electron mobility. The electron mobility of PbTe and PbSe at room temperature equals to 1730 cm 2 V −1 s −1 ( = 1.08 × 10 19 cm −3 ) and 1045 cm 2 V −1 s −1 ( = 2.38 × 10 19 cm −3 ), respectively [25]. Compared with the pristine binary compound, the scattering produced by Se and Cl substitution in the PbTe matrix scatters charge carriers additionally [22]. In addition, the interfaces and defects resulting from the distribution of PbS nanoinclusions increase the scattering of charge carriers near the interfaces [74], resulting in a decrease in the charge carriers' mobility to 954 cm 2 V −1 s −1 for = 0.0005 ( = 1.63 × 10 19 cm −3 ). As shown in Figure 20b, the theoretical calculation of S(T) is in good agreement with the experimental data ( = 0.0005) for ≥ 750 K, where the deviation for > 750 K originates from the bipolar diffusion effect. The power factor, S 2 σ(T), of the composite with nanoinclusions ( = 0.0005) is also comparable to that of n-type PbTe, as presented in Figure 20c. Figure 20d shows ( ) of the composite with = 0.0005 and the n-type PbTe compound and theoretical calculation of the Boltzmann transport. The theoretical calculation (red line) for the = 0.0005 composite fits relatively well with the experimental data. The small deviation for the high temperature ( ≥ 750 K) also comes from the bipolar effect. In the ( ) calculation, the phonon thermal conductivity as well as the electronic thermal conductivity are calculated as follows [72]: where is the speed of sound in PbTe, ℏω is the phonon energy, is the scattering time of the phonon by the Umklapp process, is the scattering time of the phonon by nanoinclusion within the geometric limit, nincl is the concentration of inclusions, is the radius of nanoinclusion, is the mole fraction of nanoinclusions, and C is determined from ℎ using the experimental value of ℎ = 2.0 W m −1 K −1 for PbTe at = 300 K. Equation (18) is only valid at the high temperature region (T > ΘD). The value of is reduced significantly in the composite with = 0.0005 compared to nanoinclusion-free n-type PbTe.
In order to understand in detail the role of nanoinclusions in enhancing the thermoelectric properties, reports shows the calculated values versus the interface potential, Vb, as shown in Figure 21a. Figure 16 shows a schematic energy band diagram of the band bending effect near the interfaces. The zero-point interface potential, Vb (0 eV), indicates that there is no phonon scattering by nanoinclusions. Theoretical values with accounting electron and electron + phonon scattering processes produced by nanoinclusions were calculated from the energy-dependent relaxation time and Born approximation. The electrostatic potential, V(r), for a single nanoinclusion can be solved by Poisson's equation, as follows [72]: where the assumed spherical nanoinclusion with a radius = 2.4 nm is = 2.5 × 10 19 cm −3 , and a 7% volume fraction of nanoinclusions. The total relaxation time of the phonon due to nanoinclusions, τ, is expressed by: where the relaxation time of the phonon for bulk is given by Equation (17) and is the relaxation time of the phonon due the to scattering by V(r) at randomly distributed nanoinclusions. The relaxation time of the phonon due the scattering can be expressed as: Figure 21. Calculated values considering electron (green dashed line) and electron + phonon scatterings by nanoinclusions (purple dashed line) and experimental data (this work for = 0.0005 and n-type PbTe) with respect to the interface potential, Vb (a), and temperature-dependent values of the experimental values and theoretical calculation (b). Reproduced with permission from [11], copyright 2017 the Royal Society of Chemistry.
where is the concentration of inclusions, is the electron transport scattering cross-section, and is the electron velocity with the p-quantum number. The electric conductivity and (electronic and lattice) thermal conductivity decreased with an increasing mole fraction of nanoinclusions, , in the composite. The Seebeck coefficient is more sensitive to the interface potential height rather than the nanoinclusion content. In order to get the maximum thermoelectric figure-of-merit, , electrical transport should be retained while the thermal conductivity should be minimized. In many cases, the decrease in electrical conductivity more significantly affects the reduction of than the decrease in thermal conductivity. Therefore, the content of nanoinclusions should not be a source of charge carriers' scattering as much as possible. In other words, the mean distance between nanoinclusions should be longer than the mean electron-free path. In this sense, the optimal mole fraction of nanoinclusions, , depends on the mean electron-free path of the material.
From the ab-initio band structure calculation, the electron-phonon interaction has a very important role in the enhancement of thermoelectric properties [75][76][77]. The electron-phonon interaction with long-range and short-range contributions, and the intravalley and intervalley transitions contribute to the electronic and thermal conductivity. The deformation-potential scattering for electron scattering by acoustic phonons makes a significant contribution to the carrier mobility due to electron-phonon interaction even at room temperature [76].
The value of n-type PbTe without nanoinclusions corresponds to the zero point of the interface potential. On the other hand, the value of the composite with = 0.0005 corresponds to the value that takes into account the electron + phonon scattering produced by nanoinclusions with a negative interface potential of Vb = −0.11 eV. In addition, theoretical values, which consider electron and phonon scatterings by nanoinclusions, are significantly higher than values, which consider electron scattering (15% or higher) only.
Therefore, the significant enhancement in the value (as much as 80.4%) for the composite with = 0.0005 is attributed to electron and phonon scattering by PbS nanoinclusions. The theoretical calculation of the value for the composite with = 0.0005 is the same as the experimental data except the high temperature range ≥ 700 K, due to the bipolar diffusion effect (Figure 21b).

Conclusion
Combining the emergence of nanoprecipitation and band engineering in bulk materials is an effective way to enhance thermoelectric performance. In p-type (PbTe)0.95−x(PbSe)x(PbS)0.5 quaternary alloys, PbS forms nanoprecipitates. The effective alloy scattering and phonon scattering by lattice dislocation and lattice distortion in the nanostructured phase produces very low lattice thermal conductivity. The observed band convergence of the two valence heavy hole, Σ, and light hole, , bands can increase the power factor (28.7 µ W cm −1 K −2 at 800 K for = 0.15). The high ZT value of 2.3 at 800 K for = 0.20 comes from the high power factor and low thermal conductivity. The alloys also have the highest average ZTavg value, which is important for many practical applications.
As the n-type counterpart of p-type (PbTe)0.95−x(PbSe)x(PbS)0.5 quaternary alloys, composites with nanoinclusion with chemical potential tuning by Cl-doping in n-type (PbTe0.93−xSe0.07Clx)0.93(PbS)0.07 composites showed high Seebeck coefficients due to the randomly distributed interface potential and energy-dependent scattering by the interface potential. The reports have shown that nanoinclusions reduce thermal conductivity due to both electron and phonon scatterings by nanoinclusions. Due to a high Seebeck coefficient and extremely low thermal conductivity, the figure-of-merit, , reaches 1.52 at 700 K for a very low Cl-doped ( = 0.0005) composite in n-type (PbTe0.93−xSe0.07Clx)0.93(PbS)0.07. The results showed that the randomly distributed interface potential, driven by the random distribution of nanoprecipitates, is an effective way to enhance thermoelectric performances. Therefore, combining strategies among band convergence, nanoprecipitation, and chemical potential tuning will be beneficiary to the enhancement of thermoelectric performance. Fine tuning and property optimization may enhance the ZT value further.
Author Contributions: Original works on PbTe based nano composites such as sample synthesis, characterization, physical properties measurements, D.G. and C.C.L; Data analysis, manuscript preparation, and supervision on the work, J.S.R. All authors have read and agreed to the published version of the manuscript.
Funding: This work was supported by the Materials and Components Technology Development Program of MOTIE/KEIT (10063286).

Conflicts of Interest:
The authors declare no conflict of interest.