High Thermoelectric Performance of Cu-Doped PbSe-PbS System Enabled by High-Throughput Experimental Screening

High Thermoelectric Performance of Cu-Doped PbSe-PbS System Enabled by High-Throughput Experimental Screening

PDF

Li You¹, Zhili Li¹, Quanying Ma¹, Shiyang He¹, Qidong Zhang¹, Feng Wang¹, Guoqiang Wu¹, Qingyi Li¹, Pengfei Luo¹, Jiye Zhang¹, Jun Luo¹^,²

Research. Vol 2020 Article ID 1736798

Less

Research. Vol 2020 Article ID 1736798

• Research Article •

High Thermoelectric Performance of Cu-Doped PbSe-PbS System Enabled by High-Throughput Experimental Screening

Full

Li You¹, Zhili Li¹, Quanying Ma¹, Shiyang He¹, Qidong Zhang¹, Feng Wang¹, Guoqiang Wu¹, Qingyi Li¹, Pengfei Luo¹, Jiye Zhang¹, Jun Luo¹^,²

Affiliations

¹School of Materials Science and Engineering, Shanghai University, 99 Shangda Road, Shanghai 200444, China

²Materials Genome Institute, Shanghai University, 99 Shangda Road, Shanghai 200444, China

About Author:

Authors’ Contributions

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.

Published: 2020-03-20 doi: 10.34133/2020/1736798

Outline

Abstract

Less

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 $Se / S$ ratio showing high thermoelectric performance has been efficiently screened out. Subsequently, based on the screened $Se / S$ ratio, the doping content of Cu has been subtly adjusted to reach the optimum carrier concentration. As a result, an outstanding peak $zT ~ 1.6$ is achieved at 873 K for a 1.8 at% Cu-doped PbSe_0.6S_0.4 sample, which is the superior value among the $n$ -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.

Cite this Article

Li You, Zhili Li, Quanying Ma, Shiyang He, Qidong Zhang, Feng Wang, Guoqiang Wu, Qingyi Li, Pengfei Luo, Jiye Zhang, Jun Luo. High Thermoelectric Performance of Cu-Doped PbSe-PbS System Enabled by High-Throughput Experimental Screening[J]. Research, 2020 , 2020 (3) : 1736798 - . DOI: 10.34133/2020/1736798

Full Text

Less

1. 　Introduction

Less

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 [1]. By virtue of the combined high-throughput (HTP) computational and experimental techniques, the discovery, optimization, and deployment of new materials will be greatly accelerated [1]. 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 [2], lithium battery [3], photovoltaics [4], 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 [10]. 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 [11]. 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,

zT

, which is defined as

zT = σ S 2 T / κ e + κ L

, where

σ

S

κ e

κ L

, and

T

are the electrical conductivity, Seebeck coefficient, electrical thermal conductivity, lattice thermal conductivity, and absolute temperature, respectively [16]. 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

zT

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 [24]. Januszko et al. prepared gradient samples using sedimentation of atoms under a strong gravitational field and obtained a large diffusion area in a Bi

1 ‐ x

x

alloy [25]. Hedegaard et al. prepared the (PbTe)

1 ‐ x

(SnTe)

x

[26] and Ge

1 ‐ x

x

[27] 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)

1 ‐ x

(GeTe)

x

sample and investigated the influence of the Ge content on the thermoelectric transport properties [28]. In addition to the Bridgman method, Gelbstein et al. prepared functionally graded (PbTe)

1 ‐ x

(SnTe)

x

samples using cold pressing followed by an annealing process and investigated the distribution of the Seebeck coefficients under different heat treatment conditions [29]. Recently, the same group paid attention to the influence of the PbI₂ content on the thermoelectric properties of (Pb_0.95Sn_0.05Te)_0.92(PbS)_0.08 by a similar approach [30].

In our previous study, we found that dynamic doping effect of interstitial Cu was highly effective in enhancing the thermoelectric performance of

n

-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 [33]. 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

Se / S

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

Se / S

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

zT ~ 1.6

at 873 K for a 1.8 at% Cu-doped PbSe_0.6S_0.4 sample is achieved, which is a superior value among the

n

-type Te-free lead chalcogenides.

2. 　Results and Discussion

Less

To screen out a satisfactory

Se / S

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

Se / S

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 PbCu_0.02S and PbCu_0.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 (

1 ‐ x

)(PbCu_0.02Se):

x

(PbCu_0.02S) (

x = 0.1, 0.2, ⋯, 0.7

). The detail of the operating principle for the system can be found in our previous report [24]. 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

Fm 3 ¯ m

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 [34]. 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

Se / S

ratio, the absolute value of the average Seebeck coefficient of each layer increases from

80 × 10 − 6

V·K^-1 to

128 × 10 − 6

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

t = 1370

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 [31]. It is also found that spot 4 with the

Se / S

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 (

m d ∗

, 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 [32]. 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

zT

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 [35]. 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 (

x ~ 0.2

) and then becomes approximately unchanged in the region of

x = 0.2 ~ 0.8

. 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

m d ∗

. 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,

Se / S

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,

m d ∗

, and lattice thermal conductivity. Thus, the

Se / S ratio ~ 1.5

of spot 4 is screened out for further thermoelectric performance modulation.

Subsequently, precise carrier concentration optimization of the PbSe_0.6S_0.4Cu_y system was carefully carried out by adjusting the Cu content, and the results are presented in Figure 3. Temperature-dependent electrical-transport properties for PbSe_0.6S_0.4Cu

y

(

y = 0.005, 0.014, 0.016, 0.018, 0.02

) 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 [31]. This reveals clearly the dynamic doping effect in the PbSe_0.6S_0.4Cu

y

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

30 × 10 − 6

Ω·m to

22 × 10 − 6

Ω·m, whereas the absolute values of the Seebeck coefficients decrease from

213 × 10 − 6

V·K^-1 to

198 × 10 − 6

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

17 × 10 − 4

W·m^-1·K^-2 at 873 K for a 1.8 at% Cu-doped PbSe_0.6S_0.4 sample.

The temperature-dependent thermal-transport properties of Cu-doped PbSe_0.6S_0.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

κ e

from

κ tot

κ e

is determined by the Wiedemann-Franz law, ,

κ e = L σ T

where

L

is the Lorenz number and

σ

is the electrical conductivity. It should be noted that

L

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 [36]. As shown in Figure 3(e), a 2 at% Cu-doped PbSe_0.6S_0.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 PbSe_0.6S_0.4 sample shows outstanding

zT

values at high temperature. As shown in Figure 3(f), a peak

zT ~ 1.6

at 873 K is achieved for the 1.8 at% Cu-doped PbSe_0.6S_0.4 sample, superior to most of the state-of-the-art

n

-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 PbSe

1 ‐ x

x

(

x = 0.1, 0.2

) samples show inferior peak

zT values ~ 1.2

and 1.3 at 873 K, which further confirms the validity and effectiveness of our HTP screening method.

3. 　Conclusions

Less

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

Se / S

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

m d ∗

, 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

Se / S

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

zT

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

n

-type Te-free lead chalcogenides. It is worth noting that, for the HTP experimental screening of

Se / S

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

Less

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 PbSe

1 ‐ x

x

Cu_0.02 (

x = 0, 1

) and PbSe_0.6S_0.4Cu

y

(

y = 0.005, 0.014, 0.106, 0.018

). The silica tube was then sealed under vacuum. To ensure the

n

-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 (

40 kV × 30 mA

). The SEM images and elemental mapping of the HTP thin slab and 2 at% Cu-doped PbSe_0.6S_0.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

κ tot = C p × D × α

, where

C p

is the heat capacity,

D

is the actual density of the sample, and

α

is the thermal diffusivity. The heat capacity was estimated via

C p k B / atom = 3.07 + 0.00047 × T ‐ 300

, 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.

Funding

Less

National Natural Science Foundation of China(51371194; 51632005; 51772186)
National Key Research and Development Program of China(2018YFA0702100; 2018YFB0703600)

References

Less

de Pablo

J. J.

, Jackson

N. E.

, Webb

M. A.

, Chen

L.-Q.

, Moore

J. E.

, Morgan

, Jacobs

, Pollock

, Schlom

D. G.

, Toberer

E. S.

, Analytis

, Dabo

, DeLongchamp

D. M.

, Fiete

G. A.

, Grason

G. M.

, Hautier

, Mo

, Rajan

, Reed

E. J.

, Rodriguez

, Stevanovic

, Suntivich

, Thornton

, and Zhao

J.-C.

, “New frontiers for the materials genome initiative,” npj Computational Materials, vol. 5, no. 1, p. 41, 2019

Norskov

J. K.

, Bligaard

, Rossmeisl

, and Christensen

C. H.

, “Towards the computational design of solid catalysts,” Nature Chemistry, vol. 1, no. 1, pp. 37–46, 2009

, Jain

, Rajput

N. N.

, Cheng

, Zhang

, Ong

S. P.

, Brafman

, Maginn

, Curtiss

L. A.

, and Persson

K. A.

, “The electrolyte genome project: a big data approach in battery materials discovery,” Computational Materials Science, vol. 103, pp. 56–67, 2015

Deng

, Wei

, Sun

, Kieslich

, Cheetham

A. K.

, and Bristowe

P. D.

, “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

Wang

, Zhang

, Snoussi

, and Zhang

, “Machine Learning Approaches for Thermoelectric Materials Research,” Advanced Functional Materials, vol. 30, no. 5, article 1906041, 2019

Yang

, Li

, Wu

, Zhang

, Chen

, and 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

, Pan

, Li

, Xu

, Ni

, Sun

, Yang

, Luo

, Xi

, Zhu

, Li

, Di Jiang

R. D.

, Shi

, Snyder

G. J.

, and Zhang

, “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

Bang

, Kim

, Wee

, Samsonidze

, and 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

Shuai

, Mao

, Song

, Zhang

, Chen

, and Ren

, “Recent progress and future challenges on thermoelectric Zintl materials,” Materials Today Physics, vol. 1, pp. 74–95, 2017

Arute

, Arya

, Babbush

, Bacon

, Bardin

J. C.

, Barends

, Biswas

, Boixo

, Brandao

F. G. S. L.

, Buell

D. A.

, Burkett

, Chen

, Chiaro

, Collins

, Courtney

, Dunsworth

, Farhi

, Foxen

, Fowler

, Gidney

, Giustina

, Graff

, Guerin

, Habegger

, Harrigan

M. P.

, Hartmann

M. J.

, Ho

, Hoffmann

, Huang

, Humble

T. S.

, Isakov

S. V.

, Jeffrey

, Jiang

, Kafri

, Kechedzhi

, Kelly

, Klimov

P. V.

, Knysh

, Korotkov

, Kostritsa

, Landhuis

, Lindmark

, Lucero

, Lyakh

, Mandra

, McClean

J. R.

, McEwen

, Megrant

, Mi

, Michielsen

, Mohseni

, Mutus

, Naaman

, Neeley

, Neill

, Niu

M. Y.

, Ostby

, Petukhov

, Platt

J. C.

, Quintana

, Rieffel

E. G.

, Roushan

, Rubin

N. C.

, Sank

, Satzinger

K. J.

, Smelyanskiy

, Sung

K. J.

, Trevithick

M. D.

, Vainsencher

, Villalonga

, White

, Yao

Z. J.

, Yeh

, Zalcman

, Neven

, and Martinis

J. M.

, “Quantum supremacy using a programmable superconducting processor,” Nature, vol. 574, no. 7779, pp. 505–510, 2019

Zhu

, Mao

, Li

, Sun

, Wang

, Zhu

, Li

, Song

, Zhou

, Fu

, He

, Tong

, Liu

, Ren

, You

, Wang

, Luo

, Sotnikov

, Bao

, Nielsch

, Chen

, Singh

D. J.

, and Ren

, “Discovery of TaFeSb-based half-Heuslers with high thermoelectric performance,” Nature Communications, vol. 10, no. 1, p. 270, 2019

Liu

, Mao

, Liu

T.-H.

, Chen

, and Ren

, “Nano-microstructural control of phonon engineering for thermoelectric energy harvesting,” MRS Bulletin, vol. 43, no. 3, pp. 181–186, 2018

Song

S. W.

, Mao

, Bordelon

, He

, Wang

Y. M.

, Shuai

, Sun

J. Y.

, Lei

X. B.

, Ren

Z. S.

, Chen

, Wilson

, Nielsch

, Zhang

Q. Y.

, and Ren

Z. F.

, “Joint effect of magnesium and yttrium on enhancing thermoelectric properties of n-type Zintl Mg_{3+ _δ_} Y_0.02Sb_1.5Bi_0.5,” Materials Today Physics, vol. 8, pp. 25–33, 2019

Ren

, Zhu

, Mao

, You

, Song

, Tong

, Bao

, Luo

, Wang

, and Ren

, “Manipulation of Ni interstitials for realizing large power factor in TiNiSn-based materials,” Advanced Electronic Materials, vol. 5, no. 7, article 1900166, 2019

Shang

P.-P.

, Dong

, Pei

, Sun

F.-H.

, Pan

, Tang

, Zhang

B.-P.

, Zhao

L.-D.

, and Li

J.-F.

, “Highly textured N-type SnSe polycrystals with enhanced thermoelectric performance,” Research, vol. 2019, article 9253132, pp. 1–10, 2019

Ren

, Zhu

, Saparamadu

, He

, Liu

, Mao

, Wang

, Nielsch

, Wang

, and Ren

, “Ultrahigh power factor in thermoelectric system Nb_0.95M_0.05FeSb (M = Hf, Zr, and Ti),” Advanced Science, vol. 5, no. 7, article 1800278, 2018

Liu

, Sun

, Mao

, Zhu

, Ren

, Zhou

, Wang

, Singh

D. J.

, Sui

, Chu

C.-W.

, and Ren

, “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

, Zhang

, Chen

, Lin

, Li

, Shen

, Witting

I. T.

, Faghaninia

, Chen

, Jain

, Chen

, Snyder

G. J.

, and Pei

, “Low-symmetry rhombohedral GeTe thermoelectrics,” Joule, vol. 2, no. 5, pp. 976–987, 2018

, Zheng

, Ge

, Lin

, Zhang

, Chen

, Chang

, and Pei

, “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

Tang

, Gao

, Lin

, Li

, Chen

, Xiong

, Li

, Chen

, and Pei

, “Manipulation of band structure and interstitial defects for improving thermoelectric SnTe,” Advanced Functional Materials, vol. 28, no. 34, article 1803586, 2018

Luo

, You

, Zhang

, Guo

, Zhu

, Gu

, Yang

, Li

, Yang

, and Zhang

, “Enhanced average thermoelectric figure of merit of the PbTe-SrTe-MnTe alloy,” ACS Applied Materials & Interfaces, vol. 9, no. 10, pp. 8729–8736, 2017

Zhang

, Song

, Wang

, Sun

, Zhu

, Dahal

, Lin

, Cao

, Zhou

, Chen

, Mao

, and Ren

, “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

, Nan

, Chen

, Zeng

, Lin

, Zhang

, Dong

, Chen

, Gu

, Li

, Chen

, Ge

, and Pei

, “Manipulation of band degeneracy and lattice strain for extraordinary PbTe thermoelectrics,” Research, vol. 2020, article 8151059, pp. 1–12, 2020

H. Y.

, Xie

R. Q.

, Peng

, Yang

, He

S. Y.

, Luo

, Sun

, Xie

S. R.

, and Luo

, “Accelerating sample preparation of graded thermoelectric materials using an automatic powder feeding system,” Advanced Manufacturing, vol. 7, no. 3, pp. 278–287, 2019

Januszko

, Stabrawa

, Ogata

, Tokuda

, Khandaker

J. I.

, Wojciechowski

, and Mashimo

, “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

Hedegaard

E. M. J.

, Mamakhel

A. A. H.

, Reardon

, and Iversen

B. B.

, “Functionally graded (PbTe)1–x(SnTe)xThermoelectrics,” Chemistry of Materials, vol. 30, no. 1, pp. 280–287, 2018

Hedegaard

E. M. J.

, Johnsen

, Bjerg

, Borup

K. A.

, and Iversen

B. B.

, “Functionally graded Ge1–xSixThermoelectrics by simultaneous band gap and carrier density engineering,” Chemistry of Materials, vol. 26, no. 17, pp. 4992–4997, 2014

Kohri

, Nishida

I. A.

, and Shiota

, “Improvement of thermoelectric properties for n-type PbTe by adding Ge,” Materials Science Forum, vol. 423-425, pp. 381–384, 2003

Gelbstein

, Dashevsky

, and Dariel

M. P.

, “Powder metallurgical processing of functionally graded p-Pb₁₋ __x_ Sn __x_ Te materials for thermoelectric applications,” Physica B, vol. 391, no. 2, pp. 256–265, 2007

Hazan

, Ben-Yehuda

, Madar

, and Gelbstein

, “Functional graded germanium-lead chalcogenide-based thermoelectric module for renewable energy applications,” Advanced Energy Materials, vol. 5, no. 11, article 1500272, 2015

You

, Liu

, Li

, Nan

, Ge

, Jiang

, Luo

, Pan

, Pei

, Zhang

, Snyder

G. J.

, Yang

, Zhang

, and Luo

, “Boosting the thermoelectric performance of PbSe through dynamic doping and hierarchical phonon scattering,” Energy & Environmental Science, vol. 11, no. 7, pp. 1848–1858, 2018

You

, Zhang

, Pan

, Jiang

, Wang

, Yang

, Pei

, Zhu

, Agne

M. T.

, Snyder

G. J.

, Ren

, Zhang

, and Luo

, “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

Volykhov

A. A.

, Yashina

L. V.

, and Shtanov

V. I.

, “Phase relations in pseudobinary systems of germanium, tin, and lead chalcogenides,” Inorganic Materials, vol. 42, no. 6, pp. 596–604, 2006

Wang

, Wang

, Cao

, and Snyder

G. J.

, “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

Wang

, LaLonde

A. D.

, Pei

, and Snyder

G. J.

, “The criteria for beneficial disorder in thermoelectric solid solutions,” Advanced Functional Materials, vol. 23, no. 12, pp. 1586–1596, 2013

Kim

H.-S.

, Gibbs

Z. M.

, Tang

, Wang

, and Snyder

G. J.

, “Characterization of Lorenz number with Seebeck coefficient measurement,” APL Materials, vol. 3, no. 4, article 041506, 2015

Appendix

Less

Year 2020 volume 2020 Issue 3

PDF

Cite this Article

BibTeX

Article Info

doi: 10.34133/2020/1736798

Receive Date：2019-12-30
Online Date：2026-10-27
Published：2020-03-20

Article Data

Affiliations

History

Received：2019-12-30
Accepted：2020-02-07

Funding

National Natural Science Foundation of China(51371194; 51632005; 51772186)

National Key Research and Development Program of China(2018YFA0702100; 2018YFB0703600)

Affiliations

¹School of Materials Science and Engineering, Shanghai University, 99 Shangda Road, Shanghai 200444, China

²Materials Genome Institute, Shanghai University, 99 Shangda Road, Shanghai 200444, China

References

Share

https://castjournals.cast.org.cn/joweb/research/EN/10.34133/2020/1736798

Share to

Scan QR to access full text

Cite this article

BibTeX

Citations

表12种不同金属材料的力学参数

科 Family	属数 Number of genus	种数 Number of species	占总种数比例 Percentage of total species (%)	属 Genus	种数 Number of species	占总种数比例 Percentage of total species (%)
鹅膏菌科Amanitaceae	2	11	5.26	鹅膏菌属 Amanita	10	4.78
小菇科 Mycenaceae	2	12	5.74	丝盖伞属 Inocybe	5	2.39
多孔菌科 Polyporaceae	8	14	6.70	蜡蘑属 Laccaria	5	2.39
红菇科 Russulaceae	3	23	11.00	小皮伞属 Marasmius	6	2.87
				小菇属 Mycena	11	5.26
				光柄菇属 Pluteus	5	2.39
				红菇属 Russula	17	8.13
				栓菌属 Trametes	5	2.39

关闭全屏

BibTeX
EndNote
RefWorks
TxT

Figure 1: Experimental high-throughput screening strategy for the Cu-doped PbSe-PbS system. (a) HTP sample preparation process. (b) HTP structure and composition characterization. (c) Seebeck coefficient mapping of the HTP thin slab. (d) Thermal-transport property evaluation of the HTP thin slab. The red opened symbol presented in (b) was taken from Ref. [34]. The figures in (b) are the lattice parameters derived from micro XRD measurements at the labeled spots (A) and the SEM-EDS mapping results (B) for the HTP thin slab along the gradient direction. The data presented in (c) are the Seebeck coefficient distribution and the average value for each section (separated by red dotted lines). (d) Represents the upper surface temperature distribution and the time-dependent temperature variation of the marked spots for the HTP thin slab. The image of the HTP thin slab was captured by an infrared camera at $t = 1373$ s (the corresponding bottom temperature of the thin slab is ~250°C).

Figure 2: Influence of S alloying on physical properties of PbSe $1 ‐ x$ S $x$ predicted by theoretical considerations. S content dependence of (a) density-of-state effective mass of the conduction band, (b) Hall mobility, (c) lattice thermal conductivity, and (d) quality factor β. The opened symbol presented in (a) was taken from Ref. [34]. The solid lines in (b) are theoretical results of the Hall mobility based on the SKB model with the assumption that acoustic phonon and alloying scattering dominate electron transport. The Pisarenko relations for PbSe-PbS at 300 K and 850 K can be found in Figure S3. The solid lines in (c) are obtained from Klemens’ model. The solid lines in (d) are the calculated results based on (a–c). The calculation detail and the physical parameters for modeling (Table S3) can be found in the Supplementary material.

Figure 3: Temperature-dependent thermoelectric transport properties for PbSe_0.6S_0.4Cu $y$ ( $y = 0.005, 0.014, 0.016, 0.018, 0.02$ ) samples: (a) Electrical resistivities; (b) Seebeck coefficients; (c) power factors; (d) total thermal conductivities; (e) lattice thermal conductivity; (f) Figure-of-merit $zT$ . The solid lines in (f) are the temperature-dependent $zT$ values for Cu-doped PbTe and PbSe with optimized Cu contents from Ref. [31, 32].

Articles: Latest Articles; Most Read; Collections

Updates: Events; News; Multimedia

About: About Us

Contact

No. 86 Xueyuan South Road, Haidian District, Beijing

100081

010-62662699(Operation and Technical Support)

qkjq@cast.org.cn

Copyright © 2025 China Association for Science and Technology. All rights reserved. For all open access content, the relevant licensing terms apply.
Sponsored by the Office of the Leading Group for Cybersecurity and Informatization of CAST, and supported by Science and Technology Review Publishing House