Recent advances in high-throughput (HTP) computational power and machine learning have led to great achievements in exploration of new thermoelectric materials. However, experimental discovery and optimization of thermoelectric materials have long relied on the traditional Edisonian trial and error approach. Herein, we demonstrate that ultrahigh thermoelectric performance in a Cu-doped PbSe-PbS system can be realized by HTP experimental screening and precise property modulation. Combining the HTP experimental technique with transport model analysis, an optimal ratio showing high thermoelectric performance has been efficiently screened out. Subsequently, based on the screened ratio, the doping content of Cu has been subtly adjusted to reach the optimum carrier concentration. As a result, an outstanding peak is achieved at 873 K for a 1.8 at% Cu-doped PbSe0.6S0.4 sample, which is the superior value among the -type Te-free lead chalcogenides. We anticipate that current work will stimulate large-scale unitization of the HTP experimental technique in the thermoelectric field, which can greatly accelerate the research and development of new high-performance thermoelectric materials.
In 2011, the Obama administration launched the “Materials Genome Initiative” (MGI) project, which was aimed at reducing the cost and shortening the research and development cycle for exploring new materials . By virtue of the combined high-throughput (HTP) computational and experimental techniques, the discovery, optimization, and deployment of new materials will be greatly accelerated . Recent advances in computational power and machine learning have led to great achievements in theoretical prediction of new functional materials in the fields of catalysis , lithium battery , photovoltaics , and thermoelectrics [5–9]. In particular, Google recently announced the realization of quantum supremacy by using a programmable superconducting processor, indicating great advancement in developing next-generation computers . It is reasonable to believe that in the foreseeable future, the ability of HTP computational power will be boosted by using the next-generation supercomputer. However, the development of the HTP experimental technique is relatively slow, which seriously lags behind that of the theoretical prediction. Historically, discovery and optimization of functional materials have long depended on the traditional Edisonian trial and error approach, which leads to costly and time-consuming procedures in verifying the massive theoretical HTP results . Therefore, it is of great theoretical and practical significance to exploit experimental HTP techniques, which may bring about revolutionary breakthrough in material research.
Optimizing the performance of a thermoelectric material is relatively difficult due to its adversely interdependent transport parameters [12–15]. Specifically, the performance of a thermoelectric material depends on the dimensionless figure of merit, , which is defined as , where , , , , and are the electrical conductivity, Seebeck coefficient, electrical thermal conductivity, lattice thermal conductivity, and absolute temperature, respectively . Thus, in order to improve thermoelectric performance, these transport parameters should be synergistically optimized. Conventionally, doping or alloying at specific lattice sites may impose great influence on both the electronic and phonon transport properties of a thermoelectric material. Up to now, the multidoping (alloying) strategy has been demonstrated to be a highly effective approach in boosting the performance of the state-of-art bulk thermoelectric materials. For example, in order to obtain ultrahigh values, doping or alloying with at least two elements is necessary for the GeTe- [17, 18], SnTe- [19, 20], and PbTe- [21–23] based compounds. Apparently, it becomes increasingly difficult to achieve satisfactory compositions for the above-mentioned system via serial experimental synthesis and characterization processes. Therefore, it should be more favorable to employ an experimental HTP technique to screen out an appropriate target composition instead of the trial and error approach before concentrating on subtle experimental optimization.
Generally, the HTP sample fabrication technique can be divided into two categories, i.e., material gene chip technology and continuously gradient composition preparation technique. The latter is usually adopted to fabricate bulk functional gradient materials (FGM) for transport property distribution characterization in thermoelectric communities due to its favorable composition features . Januszko et al. prepared gradient samples using sedimentation of atoms under a strong gravitational field and obtained a large diffusion area in a BiSb alloy . Hedegaard et al. prepared the (PbTe)(SnTe)  and GeSi  gradient samples by the Bridgman and Czochralski crystal growth method and used the Potential-Seebeck Microprobe (PSM) to investigate the electrical-transport properties. By the Bridgman method, Kohri et al. grew a Ge-graded (PbTe)(GeTe) sample and investigated the influence of the Ge content on the thermoelectric transport properties . In addition to the Bridgman method, Gelbstein et al. prepared functionally graded (PbTe)(SnTe) samples using cold pressing followed by an annealing process and investigated the distribution of the Seebeck coefficients under different heat treatment conditions . Recently, the same group paid attention to the influence of the PbI2 content on the thermoelectric properties of (Pb0.95Sn0.05Te)0.92(PbS)0.08 by a similar approach .
In our previous study, we found that dynamic doping effect of interstitial Cu was highly effective in enhancing the thermoelectric performance of -type PbSe and PbTe [31, 32]. Herein, we move forward to the Cu-doped PbSe-PbS solid solution system and intend to realize perfect doping effect through manipulating the interstitial space. According to the phase diagram, PbSe and PbS form a complete solid solution in a wide temperature range . S alloying is expected to enlarge the band gap and thus extends the working temperature range. Besides, by using more affordable elemental S, substantial reduction of the cost can be also achieved. Clearly, this system has two adjustable composition parameters, i.e., Cu content and the ratio. In order to obtain the optimized composition efficiently, we have therefore divided the experimental procedure into two steps. First, with a fixed Cu content of 2 at%, we have synthesized a bulk HTP sample with graded compositions through hot-pressing and subsequent annealing procedure. Then, by combining the HTP characterization technique with the theoretical analysis on the transport model, the best ratio was determined. Subsequently, based on the screened result, a more subtle experimental study was devoted to optimizing the carrier concentration through precise turning of Cu contents. As a result, an outstanding peak at 873 K for a 1.8 at% Cu-doped PbSe0.6S0.4 sample is achieved, which is a superior value among the -type Te-free lead chalcogenides.
2. Results and Discussion
To screen out a satisfactory ratio efficiently, we have divided the HTP screening strategy into three steps: (1) HTP sample fabrication process, (2) HTP sample characterization, and (3) transport model analysis. The probably optimized ratio is determined through combining the results of step (2) and step (3) for further thermoelectric performance optimization. As shown in Figure 1(a), the HTP sample fabrication process consists of three procedures. First of all, two samples with nominal composition PbCu0.02S and PbCu0.02Se were synthesized by a traditional vacuum melting method followed by a very short period of ball-milling (5 min) to obtain precursor powders. Then, the precursor powders were weighed and loaded into a graphite die through our home-made automatic batching and prepressing system (Movie S1) according to a composition gradient of ()(PbCu0.02Se):(PbCu0.02S) (). The detail of the operating principle for the system can be found in our previous report . After that, a high-density cylinder-like sample with composition gradient was obtained by induction hot-pressing and subsequent annealing process under vacuum. Finally, the HTP thin slab was carefully cut and polished for further characterization.
The crystal structure and composition distribution of the HTP thin slab was characterized and analyzed by microarea X-ray diffractometry and SEM-EDS. The results are presented in Figure 1(b). As shown in the upper left panel in Figure 1(b) and Figure S1(a), the XRD patterns were collected from 13 micro regions labeled from 1 to 13. The diffraction patterns for all regions can be well indexed to the space group and no obvious impurities are observed (Figure S1(a)). Furthermore, it reveals that the diffraction peaks shift to higher angles from regions 1 to 13, indicating that the lattice parameter gradually decreases along the increasing direction of S contents (Figure S1(b) and Table S1), and thus, composition-graded PbSe-PbS solid solutions have been formed. As illustrated in the upper right panel in Figure 1(b), the actual composition for each region can be roughly obtained by Vegard’s law, which is determined by the linear fitting results from literature . The composition distribution of the HTP thin slab was further determined by the SEM-EDS analysis. It is worth noting that the area for SEM-EDS characterization is approximately close to 7 mm in length, which is too large that may result in a considerable error to determine the actual composition in a one-time measurement. Thus, we have divided this area into 14 small regions and then used EDS mapping to determine elemental distributions and compositions separately. The results are shown in the bottom panel of Figure 1(b). It can be seen that Pb and Cu elements are uniformly distributed, whereas Se and S exhibit a clear composition gradient along the entire length of the sample.
As shown in Figure 1(c), the Seebeck coefficient distribution of all the samples shows a clear layered feature with gradual change along the composition gradient. More precisely, with the decreasing ratio, the absolute value of the average Seebeck coefficient of each layer increases from V·K-1 to V·K-1 (see lower region in Figure 1(c)). According to our previous study on Cu-doped lead chalcogenides, the room-temperature solubilities of the interstitial Cu in PbSe and PbS were very small [31, 32], and both Cu-doped PbSe and PbS had similar carrier concentration at room temperature. This implies that the carrier concentration of each region of the HTP sample should be also similar at room temperature (this is further confirmed by our Hall effect measurement; see the carrier concentration and mobility in Table S2). Therefore, the enhancement of the absolute value of the Seebeck coefficient with increasing S contents should be ascribed to the increased density-of-state effective mass near the conduction band edge (see the discussion below).
To screen out the favorable composition for high thermoelectric performance, both electrical and thermal properties should be measured. In the thermoelectric community, transport property measurement for the FGM was primarily focused on electrical properties, more precisely, the distribution of the Seebeck coefficient. To the best of our knowledge, HTP characterization of the thermal properties for the FGM sample remains scarce. Therefore, a home-made apparatus equipped with an in situ dynamic vacuum rapid-heating unit and an infrared camera with a 50 μm micro lens was first developed and applied to characterize thermal properties of a sample with multiple composition areas. The images of the apparatus are presented in Figure S2. Apparently, at a given baseplate temperature, the surface temperature distribution of the sample should directly reflect the thermal conductivity of the corresponding area, enabling us to qualitatively estimate the thermal conductivity of each position.
Based on the apparatus, the thermal properties of the HTP thin slab were qualitatively evaluated under the dynamic vacuum atmosphere. The bottom surface of the HTP thin slab was quickly heated up to 250°C with a heating rate of 50°C/min, kept at this temperature for 400 s, and then cooled to room temperature by forced water cooling. The time-dependent temperature variation of several composition regions for the HTP thin slab and a photograph taken at s are presented in Figure 1(d). As shown in the upper panel in Figure 1(d), the surface temperature shows obvious difference along the gradient component, indicating the obvious discrepancy of the thermal conductivity among the corresponding areas. Specifically, spots 1, 2, and spot 4 (Figure 1(d)) show relatively low surface temperatures, suggesting lower thermal conductivities for such regions. This is further confirmed by the time-dependent surface temperature variation recorded by the infrared camera. As illustrated in Figure 1(d), spots 1, 2, and 4 show obvious lower surface temperatures and more gentle heating/cooling slopes, revealing lower thermal conductivities of these areas. The low thermal conductivities of spot 1 and spot 2 could be ascribed to the hierarchical phonon scattering for the compositions with low S alloying fraction, resulting in a dramatically decreased lattice thermal conductivity over a wide temperature range . It is also found that spot 4 with the ratio approximately close to 1.5 also exhibits a lower thermal conductivity, which can be attributed to the alloying effect.
For the PbSe-PbS system, both the band gap (see Figure S4) and the density-of-state effective mass of conduction band (, see Figure 2(a)) increase upon S alloying, leading to increasing demand for optimal carrier concentration at higher temperature. According to our previous work, the space of the interstitial site is critical for a system to exhibit the dynamic doping effect. An “interstitial engineering” strategy has been thereby proposed, which is devoted to manipulating the dynamic doping effect . In the S-alloyed PbSe system, due to the smaller atomic radius of the S atom, the tetrahedral interstitial space is noticeably increased upon S alloying, which may lead to a prominent dynamic doping effect at the high temperature range and provide an ideal platform for better utilization of the dynamic doping effect in a wide temperature range. In addition, the increased solubility of interstitial Cu owing to the enlarged interstitial space will also result in enhanced point-defect scattering and reduced lattice thermal conductivity. It can be expected that the peak might be greatly boosted in a Cu-doped PbSe-PbS solid solution through the prominent dynamic doping.
However, S alloying not only leads to prominent dynamic doping effect but also results in an alloying effect, which will simultaneously enhance the scattering of electrons and phonons in crystalline solids . As shown in Figures 2(b) and 2(c), both the calculated Hall mobility and the lattice thermal conductivity show an asymmetrical variation with S content. As shown in Figure 2(b), the Hall mobility of the PbSe-PbS solid solution at 300 K drops dramatically as a small fraction of the S is alloyed () and then becomes approximately unchanged in the region of . A similar trend is also observed in the calculated lattice thermal conductivity. To theoretically evaluate the alloying effect on thermoelectric performance, the quality factor (β) has been hence calculated based on the results in thermal conductivity, Hall mobility, and . As shown in Figure 2(d), both at 300 K and 850 K, the quality factor decreases with increasing S content, indicating that the reduction of the lattice thermal conductivity cannot compensate for the deterioration in the Hall mobility. Therefore, it can be concluded that alloying effect is harmful to the overall thermoelectric performance of the PbSe-PbS system. All the calculation details can be found in the Supplementary material.
In order to balance the contributions of the dynamic doping effect and alloying effect, ratios must be carefully chosen for further performance optimization. Base on the HTP experimental results, the region near spot 4 has a moderate Seebeck coefficient and lower thermal conductivity, indicating that its related composition should have a suitable band gap, , and lattice thermal conductivity. Thus, the of spot 4 is screened out for further thermoelectric performance modulation.
Subsequently, precise carrier concentration optimization of the PbSe0.6S0.4Cuy system was carefully carried out by adjusting the Cu content, and the results are presented in Figure 3. Temperature-dependent electrical-transport properties for PbSe0.6S0.4Cu () samples are illustrated in Figures 3(a)–3(c). As shown in Figures 3(a) and 3(b), both the temperature-dependent electrical resistivities and Seebeck coefficients show a “quasi” degenerate semiconductor behavior with a visible platform around 650 K, which is similar to the Cu-doped PbSe system . This reveals clearly the dynamic doping effect in the PbSe0.6S0.4Cu system. Furthermore, a Cu-rich secondary phase is found to embed in the matrix (see SEM images in Figure S5), which is a prerequisite for Cu-doped lead chalcogenides to exhibit dynamic doping effect. Specifically, at 873 K, by varying the Cu content from 0.5 at% to 1.8 at%, the electrical resistivities substantially decrease from Ω·m to Ω·m, whereas the absolute values of the Seebeck coefficients decrease from V·K-1 to V·K-1. Further increased Cu content will lead to deterioration of electrical-transport properties, presumably due to a redundant Cu-rich secondary phase at high temperature, indicating that the solubility of the interstitial Cu at 873 K is close to 1.8 at%. Benefiting from the reinforced dynamic doping effect upon S alloying, the power factor reaches an outstanding value of W·m-1·K-2 at 873 K for a 1.8 at% Cu-doped PbSe0.6S0.4 sample.
The temperature-dependent thermal-transport properties of Cu-doped PbSe0.6S0.4 samples are illustrated in Figures 3(d) and 3(e). As shown in Figure 3(d), the total thermal conductivities decrease with increasing temperature, which is a common phenomenon in heavily doped degenerate semiconductors. The lattice thermal conductivities have been calculated by subtracting from . is determined by the Wiedemann-Franz law, , where is the Lorenz number and is the electrical conductivity. It should be noted that used in this work is calculated based on the single Kane band (SKB) model, which is lower than that calculated by the single parabolic band (SPB) model . As shown in Figure 3(e), a 2 at% Cu-doped PbSe0.6S0.4 sample possesses the lowest thermal conductivity at room temperature, which is roughly consistent with our HTP screening result. It can also be observed from Figure 3(e) that the lattice thermal conductivities decrease continuously with the increasing Cu content, which is in agreement with our previous work [31, 32].
Benefiting from the modulated interstitial space, band gap, and density-of-state effective mass, the Cu-doped PbSe0.6S0.4 sample shows outstanding values at high temperature. As shown in Figure 3(f), a peak at 873 K is achieved for the 1.8 at% Cu-doped PbSe0.6S0.4 sample, superior to most of the state-of-the-art -type Te-free lead chalcogenides. Furthermore, to demonstrate the reliability of the HTP screening result, we have also evaluated the thermoelectric performance for low S alloying fraction samples. The results are presented in Figure S6. It is demonstrated that 2 at% Cu-doped PbSeS () samples show inferior peak and 1.3 at 873 K, which further confirms the validity and effectiveness of our HTP screening method.
In summary, we have achieved an ultrahigh thermoelectric performance in a Cu-doped PbSe-PbS system through a combined strategy including experimental HTP screening, theoretical transport property analysis, and further subtle property modulation. First, a satisfactory ratio for subsequently precise property modulation has been screened out by experimental HTP technique and transport model analysis. The theoretical consideration reveals that the demand for optimal carrier concentration of the PbSe-PbS system is noticeably increased due to the increase of , interstitial space, and band gap with increasing S alloying. Thus, the PbSe-PbS system may provide an ideal platform to achieve better dynamic doping effect, especially at the high temperature range. However, according to the transport model calculation, the alloying effect results in deteriorated overall thermoelectric performance. Therefore, in order to balance the two inversely correlated effects, the ratio is finally chosen to be 1.5 based on the experimental HTP screening results. Then, the thermoelectric performance has been further optimized by subtly adjusting the Cu content. As a result, an outstanding peak 1.6 at 873 K for a 1.8 at% Cu-doped (PbSe)0.6-(PbS)0.4 sample is achieved, which is a superior value among the -type Te-free lead chalcogenides. It is worth noting that, for the HTP experimental screening of ratio, only two samples need to be presynthesized by the conventional method. This work demonstrates that the HTP experimental technique is very efficient and saves time to screen out the target compositions with high thermoelectric properties.
4. Materials and Methods
4.1. Sample Synthesis
High-purity Pb shots (>99.99%, Aladdin, China), Se chunks (>99.99% Aladdin, China), S powders (>99.9%), and Cu powder (>99.99%, Aladdin, China) were weighted and loaded into a silica tube according to the stoichiometric ratio of PbSeSCu0.02 () and PbSe0.6S0.4Cu (). The silica tube was then sealed under vacuum. To ensure the -type conduction for all samples, a small amount of excess Pb was intentionally added during the synthesis process. The silica tubes were then put in a computer-controlled furnace and slowly heated up to 1373 K at the heating rate of 1.5 K/min and then kept at this temperature for 11 h and followed by furnace cooling. The obtained ingots were ground into fine powders through an agate mortar and pestle. For the cylinder-like HTP sample with gradient compositions, it was densified by induction hot pressing in the dynamic vacuum at 823 K for 5 h under the uniaxial pressure of 50 MPa and then annealed at 673 K for 4 days. Samples for thermoelectric transport property measurement were also obtained by vacuum hot pressing under the uniaxial pressure of 70 MPa for 20 min; a high-density (>96% of theoretical density) pellet-like sample with a diameter of 10 mm was obtained for transport property measurements.
4.2. Sample Characterization
Micro area X-ray diffraction data for the HTP thin slab were collected by a Rigaku SmartLab-II diffractometer with Cu Kα radiation (). The SEM images and elemental mapping of the HTP thin slab and 2 at% Cu-doped PbSe0.6S0.4 sample were characterized by a scanning electron microscope (SEM, Gemini 300, Zeiss) equipped with an EDS detector (Oxford Instrument). The distribution of the Seebeck coefficient was scanned by a Potential-Seebeck Microprobe (PSM II, Panco Ltd., Germany). The thermal properties for the HTP thin slab was obtained by our home-made apparatus (see Figure S2). Room-temperature optical band gap was obtained by Fourier-transform infrared spectroscopy with the wavelength range from 2500 to 25000 nm (FTIR, VERTEX 70, Germany). The temperature-dependent electrical resistivity and Seebeck coefficient were characterized by the SEM-3 system (ULVAC-RIKO, Japan) under the protection of helium gas. The thermal conductivity was determined via , where is the heat capacity, is the actual density of the sample, and is the thermal diffusivity. The heat capacity was estimated via, which is commonly adopted to calculate heat capacity for lead chalcogenides. The density of the sample was determined by the Archimedes method, and the thermal diffusivity of the sample was measured by a laser flash apparatus (Netzsch LFA 457, Germany) with the Cowan model plus pulse correction.
Conflicts of Interest
The authors declare that there is no conflict of interest regarding the publication of this article.
J. Luo contributed in the conceptualization of this study; J. Luo, J.Y. Zhang, and L. You formulated the methodology of this study; L. You, Q.Y. Ma, Q.D. Zhang, and Q.Y. Li prepared the samples and conducted transport property measurements; S.Y. He, F. Wang, and L. You measured electrical and thermal properties for the HTP samples; Z.L. Li, G.Q. Wu, and L. You contributed to the SEM characterization; J.Y. Zhang, Z.L. Li, and P.F. Luo contributed to the micro XRD measurement; and L. You and J.Y. Zhang contributed to the transport model calculation. All the authors analyzed and discussed the data; J. Luo, J.Y. Zhang, and L. You wrote the manuscript with input from all authors.
We thank Prof. Y.Z. Pei and Mr. Z.L. Bu at Tongji University for their great help in the experimental measurements and discussion. This work was supported by the National Key Research and Development Program of China (Nos. 2018YFB0703600 and 2018YFA0702100) and the National Natural Science Foundation of China (Grant Nos. 51772186, 51632005, and 51371194).
Supplementary 1. Figure S1: micro area XRD patterns and derived lattice parameters for the HTP thin slab. Figure S2: images of home-made apparatus to characterize the thermal properties for an HTP thin slab. Figure S3: Pisarenko relation for PbSeS () at 300 K and 850 K. Figure S4: room-temperature XRD patterns and FTIR spectra for PbSeS () samples. Figure S5: SEM images for the 2 at% Cu-doped PbSe0.6S0.4 sample. Figure S6: temperature dependency of thermoelectric transport properties for 2 at% Cu-doped PbSeS () samples. Table S1: lattice parameters of the corresponding micro regions of the HTP sample. Table S2: room-temperature carrier concentration and mobility for ()(PbCu0.02Se):(PbCu0.02S) () samples. Table S3: physical parameters of PbSe and PbS for modeling and the details of the calculation by the SKB model and the Klemens model and the expression of the thermoelectric quality factor.
Supplementary 2. Movie S1: introduction to the automatic batching and prepressing system.
- J. J. de Pablo, N. E. Jackson, M. A. Webb et al., “New frontiers for the materials genome initiative,” npj Computational Materials, vol. 5, no. 1, p. 41, 2019.
- J. K. Norskov, T. Bligaard, J. Rossmeisl, and C. H. Christensen, “Towards the computational design of solid catalysts,” Nature Chemistry, vol. 1, no. 1, pp. 37–46, 2009.
- X. Qu, A. Jain, N. N. Rajput et al., “The electrolyte genome project: a big data approach in battery materials discovery,” Computational Materials Science, vol. 103, pp. 56–67, 2015.
- Z. Deng, F. Wei, S. Sun, G. Kieslich, A. K. Cheetham, and P. D. Bristowe, “Exploring the properties of lead-free hybrid double perovskites using a combined computational-experimental approach,” Journal of Materials Chemistry A, vol. 4, no. 31, pp. 12025–12029, 2016.
- T. Wang, C. Zhang, H. Snoussi, and G. Zhang, “Machine Learning Approaches for Thermoelectric Materials Research,” Advanced Functional Materials, vol. 30, no. 5, article 1906041, 2019.
- J. Yang, H. Li, T. Wu, W. Zhang, L. Chen, and J. Yang, “Evaluation of half-Heusler compounds as thermoelectric materials based on the calculated electrical transport properties,” Advanced Functional Materials, vol. 18, no. 19, pp. 2880–2888, 2008.
- L. Xi, S. Pan, X. Li et al., “Discovery of high-performance thermoelectric chalcogenides through reliable high-throughput material screening,” Journal of the American Chemical Society, vol. 140, no. 34, pp. 10785–10793, 2018.
- S. Bang, J. Kim, D. Wee, G. Samsonidze, and B. Kozinsky, “Estimation of electron-phonon coupling via moving least squares averaging: a method for fast-screening potential thermoelectric materials,” Materials Today Physics, vol. 6, pp. 22–30, 2018.
- J. Shuai, J. Mao, S. Song, Q. Zhang, G. Chen, and Z. Ren, “Recent progress and future challenges on thermoelectric Zintl materials,” Materials Today Physics, vol. 1, pp. 74–95, 2017.
- F. Arute, K. Arya, R. Babbush et al., “Quantum supremacy using a programmable superconducting processor,” Nature, vol. 574, no. 7779, pp. 505–510, 2019.
- H. Zhu, J. Mao, Y. Li et al., “Discovery of TaFeSb-based half-Heuslers with high thermoelectric performance,” Nature Communications, vol. 10, no. 1, p. 270, 2019.
- Z. Liu, J. Mao, T.-H. Liu, G. Chen, and Z. Ren, “Nano-microstructural control of phonon engineering for thermoelectric energy harvesting,” MRS Bulletin, vol. 43, no. 3, pp. 181–186, 2018.
- S. W. Song, J. Mao, M. Bordelon et al., “Joint effect of magnesium and yttrium on enhancing thermoelectric properties of n-type Zintl Mg3+ _δ_ Y0.02Sb1.5Bi0.5,” Materials Today Physics, vol. 8, pp. 25–33, 2019.
- W. Ren, H. Zhu, J. Mao et al., “Manipulation of Ni interstitials for realizing large power factor in TiNiSn-based materials,” Advanced Electronic Materials, vol. 5, no. 7, article 1900166, 2019.
- P.-P. Shang, J. Dong, J. Pei et al., “Highly textured N-type SnSe polycrystals with enhanced thermoelectric performance,” Research, vol. 2019, article 9253132, pp. 1–10, 2019.
- W. Ren, H. Zhu, Q. Zhu et al., “Ultrahigh power factor in thermoelectric system Nb0.95M0.05FeSb (M = Hf, Zr, and Ti),” Advanced Science, vol. 5, no. 7, article 1800278, 2018.
- Z. Liu, J. Sun, J. Mao et al., “Phase-transition temperature suppression to achieve cubic GeTe and high thermoelectric performance by Bi and Mn codoping,” Proceedings of the National Academy of Sciences of the United States of America, vol. 115, no. 21, pp. 5332–5337, 2018.
- J. Li, X. Zhang, Z. Chen et al., “Low-symmetry rhombohedral GeTe thermoelectrics,” Joule, vol. 2, no. 5, pp. 976–987, 2018.
- W. Li, L. Zheng, B. Ge et al., “Promoting SnTe as an eco-friendly solution for p-PbTe thermoelectric via band convergence and interstitial defects,” Advanced Materials, vol. 29, no. 17, article 1605887, 2017.
- J. Tang, B. Gao, S. Lin et al., “Manipulation of band structure and interstitial defects for improving thermoelectric SnTe,” Advanced Functional Materials, vol. 28, no. 34, article 1803586, 2018.
- J. Luo, L. You, J. Zhang et al., “Enhanced average thermoelectric figure of merit of the PbTe-SrTe-MnTe alloy,” ACS Applied Materials & Interfaces, vol. 9, no. 10, pp. 8729–8736, 2017.
- Q. Zhang, Q. Song, X. Wang et al., “Deep defect level engineering: a strategy of optimizing the carrier concentration for high thermoelectric performance,” Energy & Environmental Science, vol. 11, no. 4, pp. 933–940, 2018.
- Y. Wu, P. Nan, Z. Chen et al., “Manipulation of band degeneracy and lattice strain for extraordinary PbTe thermoelectrics,” Research, vol. 2020, article 8151059, pp. 1–12, 2020.
- H. Y. Pu, R. Q. Xie, Y. Peng et al., “Accelerating sample preparation of graded thermoelectric materials using an automatic powder feeding system,” Advanced Manufacturing, vol. 7, no. 3, pp. 278–287, 2019.
- K. Januszko, A. Stabrawa, Y. Ogata et al., “Influence of sedimentation of atoms on structural and thermoelectric properties of Bi-Sb alloys,” Journal of Electronic Materials, vol. 45, no. 3, pp. 1947–1955, 2016.
- E. M. J. Hedegaard, A. A. H. Mamakhel, H. Reardon, and B. B. Iversen, “Functionally graded (PbTe)1–x(SnTe)xThermoelectrics,” Chemistry of Materials, vol. 30, no. 1, pp. 280–287, 2018.
- E. M. J. Hedegaard, S. Johnsen, L. Bjerg, K. A. Borup, and B. B. Iversen, “Functionally graded Ge1–xSixThermoelectrics by simultaneous band gap and carrier density engineering,” Chemistry of Materials, vol. 26, no. 17, pp. 4992–4997, 2014.
- H. Kohri, I. A. Nishida, and I. Shiota, “Improvement of thermoelectric properties for n-type PbTe by adding Ge,” Materials Science Forum, vol. 423-425, pp. 381–384, 2003.
- Y. Gelbstein, Z. Dashevsky, and M. P. Dariel, “Powder metallurgical processing of functionally graded p-Pb1− _x_ Sn _x_ Te materials for thermoelectric applications,” Physica B, vol. 391, no. 2, pp. 256–265, 2007.
- E. Hazan, O. Ben-Yehuda, N. Madar, and Y. Gelbstein, “Functional graded germanium-lead chalcogenide-based thermoelectric module for renewable energy applications,” Advanced Energy Materials, vol. 5, no. 11, article 1500272, 2015.
- L. You, Y. Liu, X. Li et al., “Boosting the thermoelectric performance of PbSe through dynamic doping and hierarchical phonon scattering,” Energy & Environmental Science, vol. 11, no. 7, pp. 1848–1858, 2018.
- L. You, J. Zhang, S. Pan et al., “Realization of higher thermoelectric performance by dynamic doping of copper in n-type PbTe,” Energy & Environmental Science, vol. 12, no. 10, pp. 3089–3098, 2019.
- A. A. Volykhov, L. V. Yashina, and V. I. Shtanov, “Phase relations in pseudobinary systems of germanium, tin, and lead chalcogenides,” Inorganic Materials, vol. 42, no. 6, pp. 596–604, 2006.
- H. Wang, J. Wang, X. Cao, and G. J. Snyder, “Thermoelectric alloys between PbSe and PbS with effective thermal conductivity reduction and high figure of merit,” Journal of Materials Chemistry A, vol. 2, no. 9, pp. 3169–3174, 2014.
- H. Wang, A. D. LaLonde, Y. Pei, and G. J. Snyder, “The criteria for beneficial disorder in thermoelectric solid solutions,” Advanced Functional Materials, vol. 23, no. 12, pp. 1586–1596, 2013.
- H.-S. Kim, Z. M. Gibbs, Y. Tang, H. Wang, and G. J. Snyder, “Characterization of Lorenz number with Seebeck coefficient measurement,” APL Materials, vol. 3, no. 4, article 041506, 2015.