Contrasting Thermoelectric Transport Behaviors of <i>p</i>-Type PbS Caused by Doping Alkali Metals (Li and Na)

Contrasting Thermoelectric Transport Behaviors of p-Type PbS Caused by Doping Alkali Metals (Li and Na)

PDF

Zhenghao Hou¹, Dongyang Wang¹, Jinfeng Wang², Guangtao Wang², Zhiwei Huang¹, Li-Dong Zhao¹

Research. Vol 2020 Article ID 4084532

Less

Research. Vol 2020 Article ID 4084532

• Research Article •

Contrasting Thermoelectric Transport Behaviors of p-Type PbS Caused by Doping Alkali Metals (Li and Na)

Full

Zhenghao Hou¹, Dongyang Wang¹, Jinfeng Wang², Guangtao Wang², Zhiwei Huang¹, Li-Dong Zhao¹

Affiliations

¹School of Materials Science and Engineering, Beihang University, Beijing 100191, China

²School of Physics, Henan Normal University, Xinxiang 453007, China

Published: 2020-12-20 doi: 10.34133/2020/4084532

Outline

Abstract

Less

PbS is a latent substitute of PbTe thermoelectric materials, which is on account of its superiority in low cost and earth abundance. Here, the thermoelectric transport properties of p-type PbS by doping alkali metals (Na and Li) are investigated and it is verified that Li is a more effective dopant than Na. By introducing Li, the electrical and thermal transport properties were optimized collectively. The electrical transport properties were boosted remarkably via adjusting carrier concentration, and the maximum power factor (PF_max) of ~11.5 μW/cmK² and average power factor (PF_ave) ~9.9 μW/cmK² between 423 and 730 K in Pb_0.99Li_0.01S were achieved, which are much higher than those (~9.5 and ~7.7 μW/cmK²) of Pb_0.99Na_0.01S. Doping Li and Na can weaken the lattice thermal conductivity effectively. Combining the enlarged PF with suppressed total thermal conductivity, a maximum ZT ~0.5 at 730 K and a large average ZT ~0.4 at 423-730 K were obtained in p-type Pb_0.99Li_0.01S, which are higher than ~0.4 and ~0.3 in p-type Pb_0.99Na_0.01S, respectively.

Cite this Article

Zhenghao Hou, Dongyang Wang, Jinfeng Wang, Guangtao Wang, Zhiwei Huang, Li-Dong Zhao. Contrasting Thermoelectric Transport Behaviors of p-Type PbS Caused by Doping Alkali Metals (Li and Na)[J]. Research, 2020 , 2020 (12) : 4084532 - . DOI: 10.34133/2020/4084532

Full Text

Less

1. 　Introduction

Less

The search for reliable and environmentally friendly new energy has attracted worldwide attention because of the shortage of fossil energy. A thermoelectric device is capable of transforming heat into electric energy immediately, which has shown great prospect in clean energy field [1–5]. The thermoelectric device efficiency is positively associated with the dimensionless figure of merit [6–9],

ZT = S 2 σ T / κ tot

, where

T

represents absolute temperature,

S

expresses the Seebeck coefficient,

σ

denotes electrical conductivity, and

κ tot

represents total thermal conductivity comprising electronic (

κ ele

) and lattice (

κ lat

) contributions [2, 10, 11].

Lead telluride- (PbTe-) based materials, as a kind of medium-temperature operating materials, have attracted extensive research interests on account of outstanding thermoelectric performance [12–14]. However, considering the high costs and low earth abundance of the Te element, the thermoelectric materials with rich resources should be developed. To date, one research hotpot in this field is to find an alternative material to substitute PbTe which possesses prominent thermoelectric properties [14–16]. As a similar alternative of PbTe, PbS possesses a NaCl structure and an alike band structure. Nevertheless, the poor electrical properties and large lattice thermal conductivity historically make PbS become an inferior thermoelectric material [11]. Aimed at solving the shortcomings of PbS, the approaches including carrier concentration optimization [17], band manipulation [18], and microstructure engineering [19–22] have been proved as effective strategies to manipulate electrical properties and thermal conductivity; the achievements realized through the above strategies well elucidate the potential performance of PbS.

Usually, doping is a powerful tactic to enhance ZT. Doping is essential in adjusting carrier concentration, and it is the prerequisite to gain a high ZT as all of those thermoelectric properties are interlinked by carrier concentration. On the assumption that the scattering or band structure is not modified obviously by a dopant, the Seebeck coefficient and electrical conductivity of degenerated semiconductor with a single parabolic band can be given using Equations (1) and (2) [23, 24].

(1) S = 8 π 2 k B 2 3 e h 2 m ∗ T π 3 n 2 / 3, (2) σ = n e μ n,

where

k B

displays the Boltzmann constant,

e

shows the electron charge,

h

denotes the Planck constant,

m ∗

expresses effective mass,

n

is the carrier concentration, and

μ n

represents the carrier mobility. Apparently,

S

and

σ

are in an inversely proportional relationship. Therefore, adjusting carrier concentration through balancing the relationship between

S

and

σ

is an important method to boost power factors. Doping with different elements may induce diverse impact on carrier concentration optimization and band structure manipulation. For example, p-type Na-doped PbS with CdS as second phases attains a large ZT ~1.3 at 923 K owing to extensive phonon scattering by nanophase precipitates and better electrical transport [25]; the ZT of p-type Tl-doped PbTe reaches ~1.5 at 773 K due to greatly enhanced Seebeck coefficients by deformation of electronic density of states [17]. It is meaningful to speculate the impacts by doping other alkali metals on thermoelectric performance in PbS.

In this article, we focused on PbS, which contains highly earth-abundant elements and owns higher melting point compared to PbTe. The thermoelectric properties in PbS doped by Li and Na were investigated systematically. The consequences reveal that the electrical and thermal properties were optimized synchronously through alkali metal doping. The electrical properties were improved through adjusting carrier concentration, and the PF_max of Pb_0.99Li_0.01S reached ~11.5 μW/cmK², which is far greater than ~9.5 μW/cmK² in Pb_0.99Na_0.01S. Both Li and Na can bring down the

κ lat

because of point defects in PbS matrix. Li was more effective than Na in reducing

κ lat

on account of larger mass and strain field fluctuations. Integrating enhanced PF and reduced

κ tot

, a higher ZT value ~0.5 at 730 K and average ZT ~0.4 at 423-730 K can be reached in Pb_0.99Li_0.01S, which are higher than ~0.4 and~0.3 in p-type Pb_0.99Na_0.01S, respectively, indicating that Li doping can improve thermoelectric performance of PbS more effectively than Na doping.

2. 　Experimental Section

Less

2.1. 　Preparation Method

High-purity chemicals of Pb particle (99.99%), S (99.99%), Na (99.99%), and Li (99.99%) were weighed and loaded into carbon-coated quartz ampules under a N₂-filled glove box. The ampules of the chemicals were evacuated under vacuum and flame-sealed. The pure chemicals were gradually warmed up to 723 K in 12 h, elevated to 1423 K in 7 h before keeping stable at 1423 K for 6 h, and finally naturally cooled to indoor temperature. The prepared specimens were pulverized and filtered with 400-mesh sieves for sintering through spark plasma sintering (SPS-211LX) using a pressure of 50 MPa at 923 K for 10 min.

2.2. 　Thermoelectric Properties

The acquired cylinder bulk materials were incised for measuring relevant thermoelectric properties. The CTA system was applied to measure electrical parameters (

S

and

σ

) at 300-730 K under He gas, and the samples were polished in a rectangular shape of

10 mm × 3 mm × 3 mm

. The surfaces of measured samples were sprayed with thin-layer BN, which can inhibit volatilization and protect instrument [26]. The cylindrical disks with thickness of 1 mm and diameter of 6 mm were used to measure the thermal diffusivity (

D

). The thermal conductivity is computed through

κ tot = D C p ρ

, and the thermal diffusivity was characterized using a Netzsch LFA457 instrument with a laser flash method [27]. A thin graphite film on the surface of samples was utilized to cut down errors of emissivity for testing

D

. The density (

ρ

) was obtained based on mass and volume. All the densities of samples are around 7.2 g/cm³. The heat capacity (

C p

) was computed using the Debye model [28].

2.3. 　X-Ray Diffraction

The phase structure was investigated using an X-ray diffraction technique with D/MAX2200pc system with CuKα at

2 θ = 20 – 80 °

(Rigaku, Japan, 40 kV, 40 mA). The scanning speed and step size of the XRD measurement are 6° min^-1 and 0.02°, respectively.

2.4. 　Theoretical Calculations

The density functional theory (DFT) calculations were acquired through a projector-augmented wave (PAW) strategy [29, 30] with the Vienna Ab initio Simulation Package (VASP) [31]. The Perdew-Burke-Ernzerhof (PBE) exchange-correlation functional was used to model crystal and electronic structure. The used kinetic cutoff energy of plane waves is 500 eV. A

3 × 3 × 3

supercell (Pb₂₇S₂₇) was constructed to evaluate the defect formation energy of Li- (Pb₂₆LiS₂₇) and Na- (Pb₂₆NaS₂₇) doped systems. The internal coordinates of all atoms are entirely relaxed while the maximum residual ionic force is lower than 0.01 eV Å^-1, and the total energy difference approaches 10^-6 eV.

The formation energy of the defect

Q

(

Q = Na

, Li) in charge

q

is calculated by [32, 33]

(3) ∆ H f PbS, Q q = E PbS, Q q − E PbS + ∑ i n i E i + μ i + q E F + E V + ∆ V,

where

E PbS, Q q

and

E PbS

refer to the total energy of defect

α

contained and undoped PbS in supercells with the same dimension, separately.

μ i

E i

, and

n i

are the chemical potential, average energy of element in its most stable crystal structure, and the number of the atom

i

added to (

n i < 0

) or taken from (

n i > 0

) the host, respectively.

E F

presents the Fermi level relative to energy location of valence band maximum (

E v

), which changes between 0 and band gap in PbS. The correction term

Δ V

is adopted to arrange the reference potential between the defect-containing and pure supercells with the same size [34]. All characteristic values are recombined to the 1s core level of the atom farthest from the defect [32, 35].

The formation energy depends on the chemical potential of each element, which is related to the off-stoichiometric degree (Pb- or S-rich condition). The different off-stoichiometric degrees will result in different chemical potential and formation energy. The upper and lower boundary chemical potential (

μ i

) is determined by the off-stoichiometric degree and the stability against precipitation of elemental Pb, S, Li, and Na:

(4) μ Pb, μ S, μ Li, μ Na ≤ 0 .

The host compounds are obtained from the sum of the chemical potentials of Pb and S:

(5) μ Pb + μ S = ∆ H f PbS,

where

Δ H f PbS

is the formation energy of PbS in a rock-salt structure.

The second phase of Na₂S, Li₂S, and PbS₂ should be avoided, and the corresponding formation energy should be larger than the sum of elemental chemical potential:

(6) 2 μ Na + μ S < ∆ H f Na 2 S, 2 μ Li + μ S < ∆ H f Li 2 S, μ Pb + 2 μ S < ∆ H f PbS 2 .

3. 　Results and Discussion

Less

Figure 1 demonstrates the detailed information of XRD results. All specimens possess a single phase of cubic PbS. The data peak transfers to low angle range as the Li and Na content was added, which indicates that Li and Na are doped into PbS lattice.

Figure 2 depicts electrical properties in PbS with Li and Na doping. It can be clearly observed from Figures 2(a) and 2(b) that the

σ

falls off when temperature rises, except for low doping samples of Pb_1-xNa_xS (

x = 0.005

and 0.0075).

For Li-doped samples, the

σ

possesses a tendency to increase first and then decrease with the stoichiometry of Li increasing and reaches to its maximum 777 S/cm in Pb_0.995Li_0.005S, as shown in Figure 2(a). However, the Na doping presents different results. As presented in Figure 2(b), the

σ

has an increased trend with increasing of Na content, and the maximum

σ

of 1274 S/cm can be realized in Pb_0.98Na_0.02S. The

σ

is positively correlated with

n H

and

μ n

from Equation (2), which are determined by the solid solubility and the different scattering mechanisms, respectively. The continued increase in

σ

of Pb_1-xNa_xS is mainly due to the fact that the higher solubility of Na than Li was caused by better ion radius matching (

r Li + = 0.76 Å

r Na + = 1.02 Å

r Pb 2 + = 1.26 Å

As displayed in Figures 2(c) and 2(d), different from the undoped PbS, the Seebeck coefficients (

S

) for all doped samples are positive, indicating that Li and Na are effective p-type dopants in PbS. For Li-doped samples, the

S

present the trend of first decreasing and then increasing with an increasing Li content. For Na-doped samples, the trend is reversed that the

S

increase first and then decrease with an increasing Na content. These diametrically opposite trends reflect the contrary change of carrier concentration (

n H

) in those materials since the

S

are negatively correlated with

n H

To evaluate the doping efficiency of Li and Na in PbS, the formation energy of potential defect was calculated and shown in Figures 3(a) and 3(b). The lower formation energy of Na_Pb indicates the spontaneous formation of Na_Pb in any conditions, which is even lower than that in V_Pb. However, the Li_Pb has higher formation energy under Pb- and S- rich situations. In an equilibrium theory, the defect concentration can be evaluated by the formation energy

Δ H

, expressed as

n i = N i × e − ∆ H / k T

[36]. Thus, the larger formation energy of Li leads to a lower

n H

and a larger

S

(Figure 2(c)).

As presented in Figure 2(e), for Li-doped samples, a higher PF can be obtained in Pb_0.99Li_0.01S in a broad temperature range, and the peak value can reach 11.5 μW/cmK² at 450 K. The peak PF for Pb_0.99Li_0.01S is much higher than Pb_0.99Na_0.01S which is ascribed to the lower

n H

, namely, adjusting carrier concentration to an optimized scope. Figures 3(c) and 3(d) show the carrier mobility and carrier concentration at room temperature which are calculated by the carrier effective mass of PbS (

m ∗ = 0.38 m 0

) [25]. According to Rowe and Bhandari’s study [37], the

S

decreases and the

σ

increases as the

n H

increases and the PF maximizes at a suitable

n H

for a semiconductor. Therefore, adjusting the

n H

to a reasonable range is the key factor to obtain higher PF. Compared with Pb_0.99Na_0.01S, Li doping leads to a relative lower

n H

and higher PF. Figures 2(g) and 2(h) show the maximum power factor (PF_max) at 300-730 K and average power factor (PF_ave) within 423-730 K of p-type PbS samples. The PF_ave is calculated by Equation (7) in which the

T h

and

T c

are the temperatures of hot and cold ends. The PF_ave represents the overall capacity and level of electrical transports over a specified wide temperature range. The PF_max and PF_ave of the Pb_0.99Li_0.01S sample are 11.5 and 9.9 μW/cmK², respectively, higher than those of Pb_0.99Na_0.01S which are 9.5 and 7.7 μW/cmK². The present results reveal that the different dopants can reach the PF_max under their proper

n H

, which is strictly determined by the solid solubility of dopants in PbS.

(7) PF ave = 1 T h − T c ∫ T c T h PF d T .

Figures 4(a) and 4(b) depict the

κ tot

which decreases monotonically with the increase of temperature. The

κ tot

of Pb_1-xLi_xS is lower than that in the undoped PbS, which is different from the larger content of Pb_1-xNa_xS since

σ

is higher. The

C p

of Pb_0.99Li_0.01S and Pb_0.99Na_0.01S is 0.2094 J/g·K and 0.2092 J/g·K at 730 K, respectively. The

C p

of Li-doped samples is similar to the C_p of Na-doped samples at the same content and temperature. The

κ tot

includes lattice thermal conductivity and electronic thermal conductivity (

κ tot = κ lat + κ ele

) [22, 25], where the relationship between

κ ele

σ

, and Lorenz number (

L

) described in Equation (8) indicates that the

κ ele

is proportional to

σ

[38, 39].

(8) κ ele = L σ T .

The Lorenz number was obtained through calculating the Seebeck coefficient and integral chemical potentials [40]. Figures 4(c) and 4(d) show the Lorenz number in all samples as function of temperature. Higher

L

and

σ

lead to the larger

κ ele

than those in undoped PbS, as revealed through Figures 4(e) and 4(f). Thus, the reduction in

κ tot

is primarily caused by the decrease of

κ lat

Figures 5(a) and 5(b) show that the

κ lat

of all doped samples is lower than that of the undoped sample. The point defect scattering presumably reduces the

κ lat

by Li and Na doping. Obviously, the

κ lat

decreases with the increasing dopant content. More importantly, Li and Na are both effective in reducing the

κ lat

. To understand the phonon transports in Li- (Na) doped PbS, we adopted the Callaway model to evaluate point defect scattering caused by Li and Na doping [28, 41, 42].

When the temperature is higher than the Debye temperature, the point defect is an intensive scattering center to reduce the

κ lat

. According to the Callaway model [28, 42, 43], the ratio of the

κ lat

between the defect-containing material and host material can be written as

(9) κ lat κ lat, p = tan − 1 u u,

in which

κ lat

and

κ lat, p

represent the lattice thermal conductivities in doped and parent materials, separately. The parameter

u

is described using

(10) u = π 2 θ D Ω h v a 2 K L, p Γ 1 / 2,

in which

h

Ω

v a

, and

θ D

express the Planck constant, average atom volume, average sound velocity, and Debye temperature, separately. The imperfection scaling parameter (

Γ

) indicates that the phonon scattering intensity by atomic scale defects contains mass fluctuation

Γ M

and strain field fluctuation

Γ S

. The phenomenological adjustable parameter (

ε

) regulates the uncertainty of

Γ S

. The imperfection scaling parameter

Γ

and the phenomenological adjustable parameter

ε

are expressed by the following equations [42]:

(11) Γ = Γ M + ε Γ S, ε = 29 6.4 × γ 1 + ν p 1 − ν p 2,

where

ν p

displays the Poisson ratio, which is calculated using the longitudinal (

ν l

) and transverse (

ν s

) acoustic velocities. The acoustic velocity of PbS was adopted in Poisson ratio and Grüneisen parameter (

γ

) calculation by the following equations:

(12) ν p = 1 − 2 ν s / ν l 2 2 − 2 ν s / ν l 2, γ = 23 1 + ν p 2 − 3 ν p .

When Pb sites are replaced by Li (Na), no change happens on the position of S,

Γ S = 0

, which is defined by [44, 45]

(13) Γ Pb 1 − x Q x S = 12 M Pb, Q M ¯ 2 Γ Pb, Q, Γ Pb, Q = Γ M, Pb, Q + ε Γ S, Pb, Q, Γ M, Pb, Q = x 1 − x ∆ M M Pb, Q 2,

where

∆ M = M Pb − M Q

, and

M Pb, Q = 1 − x M Pb + x M Q

(14) Γ S, Pb, Q = x 1 − x ∆ r r Pb, Q 2,

where

∆ r = r P b − r Q

, and

r Pb, Q = 1 − x r Pb + x r Q

Then,

(15) Γ Pb 1 − x Q x S = 12 M Pb, Q M ¯ 2 x 1 − x ∆ M M Pb, Q 2 + ε ∆ r r Pb, Q 2 .

The calculated mass fluctuations

Γ M, Pb, Q

and strain field fluctuations

Γ S, Pb, Q

have been given by Table 1. The higher deviations in atomic radius and mass between Pb and Li lead to larger Γ_(Pb,Li)

than

Γ Pb, Na

, indicating more effective decreasing of

κ lat

by Li doping. In Figures 5(c) and 5(d), the calculated results based on the Callaway model exhibit the same trend with the experimental data. The huge deviation may result from the formation of nanostructure even though in moderate doping concentration. This phenomenon confirms that Li and Na could both play effective roles in suppressing

κ lat

The temperature-dependent ZT of PbS doped by Li and Na are presented in Figures 6(a) and 6(b). Pb_1-xLi_xS samples exhibit larger ZT than Pb_1-xNa_xS samples. The maximum ZT (ZT_max) in Pb_0.99Li_0.01S attained ~0.5 when

T = 730 K

, which is higher than that in Pb_0.99Na_0.01S. The better thermoelectric performance of Li-doped samples is mainly due to the higher PF which results from the obtained proper

n H

range and the slightly lower

κ tot

from a more effective point defect scattering.

The variation trends of maximum ZT (ZT_max) and average ZT (ZT_ave), calculated by Equation (16), are consistent with PF_max and PF_ave, as displayed in Figures 6(c) and 6(d). The ZT_max from room temperature to 730 K and ZT_ave within 423-730 K are 0.5 and 0.4 in Pb_0.99Li_0.01S, which is much higher than Pb_0.99Na_0.01S (0.4 and 0.3). The quality factor

B

is a parameter for estimating the optimal thermoelectric properties of a specific material according to the effective mass model, and the quality factor

B

is obtained by Equation (16). The weighted mobility

μ w

is calculated by the electrical conductivity and Seebeck coefficient according to Equation (17) [46, 47].

(16) B = 9 μ W κ lat T 300 5 / 2, (17) μ W = 3 σ 8 π e F 0 η h 2 2 m e k B T 3 / 2,

in which

m e

and

e

are unit mass of free electron and the electron charge, respectively.

F n η

represents the Fermi integral with

n = 0

and is calculated by the following equations.

(18) F n η = ∫ 0 ∞ x n 1 + e x − η d x, S = ± k B e r + 5 / 2 F r + 3 / 2 η r + 3 / 2 F r + 1 / 2 η − η,

in which

r

shows the scattering factor and equals -1/2 here and

η

is the reduced chemical potential [46].

The calculated quality factors of Pb_0.99Li_0.01S and Pb_0.99Na_0.01S at 730 K are 0.4 and 0.2, respectively. The quality factor of Pb_0.99Li_0.01S is about twice higher than that of Pb_0.99Na_0.01S, so the ZT of Pb_0.99Li_0.01S is higher, which is caused by the enhanced PF by adjusting

n H

in a reasonable range. The thermoelectric conversion efficiencies are calculated by Equation (20) [28]:

(19) ZT ave = 1 T h − T c ∫ T C T h ZT d T, (20) η = T h − T c T h 1 + ZT ave − 1 1 + ZT ave + T c / T h,

in which

T h

and

T c

represent the temperature in hot and cold end, respectively. The maximum calculated thermoelectric conversion efficiency based on single leg is ~4.8% in Pb_0.99Li_0.01S which is higher than Pb_0.99Na_0.01S (~3.4%) when

T h = 730

K and

T c = 423

K, indicating Li is a valid dopant to regulate the thermoelectric performance through tuning

n H

4. 　Conclusion

Less

This work indicates that Li doping is more effective than Na doping in thermoelectric performance optimization in PbS. The boosted thermoelectric performance of Li-doped PbS is completed by enhancing the PF through regulating

n H

in a reasonable range. The PF_max and PF_ave between 423 and 730 K of Pb_0.99Li_0.01S reached ~11.5 and~9.9 μW/cmK², which are much better compared with ~9.5 and ~7.7 μW/cmK² of Pb_0.99Na_0.01S. Pb_1-xLi_xS samples possess slightly smaller

κ lat

than that of Pb_1-xNa_xS because of larger mass and strain field fluctuations. At last, higher ZT_max ~0.5 at 730 K and ZT_ave ~0.4 at 423 K-730 K can be obtained in Pb_0.99Li_0.01S. The calculated thermoelectric conversion efficiency ~4.8% is achieved in Pb_0.99Li_0.01S with

T h = 730

K and

T c = 423

K. In the future, the ZT for Li-doped PbS can also be raised through nanostructuring, manipulating band structures, and other approaches.

Funding

Less

China Postdoctoral Science Foundation(2019M650429)
Henan Normal University
Beihang University
National Postdoctoral Program for Innovative Talents(BX20200028)
Academic Excellence Foundation of BUAA for PhD Students
National Science Fund for Distinguished Young Scholars(51925101)
111 Project(B17002)
Beijing Municipal Natural Science Foundation(JQ18004)
National Key Research and Development Program of China(2018YFB0703600; 2018YFA0702100)
National Natural Science Foundation of China(51671015; 51772012)

References

Less

Bell

L. E.

, “Cooling, heating, generating power, and recovering waste heat with thermoelectric systems,” Science, vol. 321, no. 5895, pp. 1457–1461, 2008

Tan

, Zhao

L.-D.

, and Kanatzidis

M. G.

, “Rationally designing high-performance bulk thermoelectric materials,” Chemical Reviews, vol. 116, no. 19, pp. 12123–12149, 2016

Chen

Z. G.

, Shi

, Zhao

L. D.

, and Zou

, “High-performance SnSe thermoelectric materials: progress and future challenge,” Progress in Materials Science, vol. 97, pp. 283–346, 2018

Xiao

, Chang

, Pei

, Wu

, Peng

, Zhou

, Gong

, He

, Zhang

, Zeng

, and Zhao

L. D.

, “Origin of low thermal conductivity in SnSe,” Physical Review B, vol. 94, no. 12, article 125203, 2016

Tang

, Wei

, Zhang

, Li

, Wang

, Xu

, Chang

, Wang

, du

, and Zhao

L. D.

, “Realizing high figure of merit in phase-separated polycrystalline Sn_1-xPb_xSe,” Journal of the American Chemical Society, vol. 138, no. 41, pp. 13647–13654, 2016

Liu

, Zhao

L.-D.

, Zhu

, Liu

, Li

, Yu

, Liu

D. B.

, Xu

, Lin

Y. H.

, and Nan

C. W.

, “Synergistically optimizing electrical and thermal transport properties of BiCuSeO via a dual-doping approach,” Advanced Energy Materials, vol. 6, no. 9, article 1502423, 2016

Zhao

L.-D.

, Zhang

, Wu

, Tan

, Pei

, Xiao

, Chang

, Wu

, Chi

, Zheng

, Gong

, Uher

, He

, and Kanatzidis

M. G.

, “Enhanced thermoelectric properties in the counter-doped SnTe system with strained endotaxial SrTe,” Journal of the American Chemical Society, vol. 138, no. 7, pp. 2366–2373, 2016

Xiao

, Wu

, Li

, Yin

, Pei

, Zhang

, Fu

, Chen

, Pennycook

S. J.

, Huang

, He

, and Zhao

L. D.

, “Remarkable roles of Cu to synergistically optimize phonon and carrier transport in n-type PbTe-Cu₂Te,” Journal of the American Chemical Society, vol. 139, no. 51, pp. 18732–18738, 2017

Qin

B. C.

, Xiao

, Zhou

Y. M.

, and Zhao

L. D.

, “Thermoelectric transport properties of Pb-Sn-Te-Se system,” Rare Metals, vol. 37, no. 4, pp. 343–350, 2018

Chang

, Wu

, He

, Pei

, Wu

C. F.

, Wu

, Yu

, Zhu

, Wang

, Chen

, Huang

, Li

J. F.

, He

, and Zhao

L. D.

, “3D charge and 2D phonon transports leading to high out-of-planeZTin n-type SnSe crystals,” Science, vol. 360, no. 6390, pp. 778–783, 2018

Pei

Y. L.

, and Liu

, “Electrical and thermal transport properties of Pb-based chalcogenides: PbTe, PbSe, and PbS,” Journal of Alloys and Compounds, vol. 514, pp. 40–44, 2012

Xiao

, Wu

, Cui

, Wang

, Fu

, Zhang

, Chen

, He

, Pennycook

S. J.

, and Zhao

L. D.

, “Realizing high performance n-type PbTe by synergistically optimizing effective mass and carrier mobility and suppressing bipolar thermal conductivity,” Energy & Environmental Science, vol. 11, no. 9, pp. 2486–2495, 2018

Xiao

, Wu

, Wang

, Niu

, Pei

, Zhang

, Spanopoulos

, Witting

I. T.

, Li

, Pennycook

S. J.

, Snyder

G. J.

, Kanatzidis

M. G.

, and Zhao

L. D.

, “Amphoteric indium enables carrier engineering to enhance the power factor and thermoelectric performance inn‐Type AgnPb100InnTe100+2n(LIST),” Advanced Energy Materials, vol. 9, no. 17, article 1900414, 2019

Zhao

L.-D.

, Wu

H. J.

, Hao

S. Q.

, Wu

C. I.

, Zhou

X. Y.

, Biswas

, He

J. Q.

, Hogan

T. P.

, Uher

, Wolverton

, Dravid

V. P.

, and Kanatzidis

M. G.

, “All-scale hierarchical thermoelectrics: MgTe in PbTe facilitates valence band convergence and suppresses bipolar thermal transport for high performance,” Energy & Environmental Science, vol. 6, no. 11, pp. 3346–3355, 2013

Hsu

, Loo

, Guo

, Chen

, Dyck

J. S.

, Uher

, Hogan

, Polychroniadis

E. K.

, and Kanatzidis

M. G.

, “Cubic AgPb_mSbTe_2+m: bulk thermoelectric materials with high figure of merit,” Science, vol. 303, no. 5659, pp. 818–821, 2004

Biswas

, He

, Blum

I. D.

, Wu

C. I.

, Hogan

T. P.

, Seidman

D. N.

, Dravid

V. P.

, and Kanatzidis

M. G.

, “High-performance bulk thermoelectrics with all-scale hierarchical architectures,” Nature, vol. 489, no. 7416, pp. 414–418, 2012

Heremans

J. P.

, Jovovic

, Toberer

E. S.

, Saramat

, Kurosaki

, Charoenphakdee

, Yamanaka

, and Snyder

G. J.

, “Enhancement of thermoelectric efficiency in PbTe by distortion of the electronic density of states,” Science, vol. 321, no. 5888, pp. 554–557, 2008

Tan

, Shi

, Hao

, Chi

, Zhao

L. D.

, Uher

, Wolverton

, Dravid

V. P.

, and Kanatzidis

M. G.

, “Codoping in SnTe: enhancement of thermoelectric performance through synergy of resonance levels and band convergence,” Journal of the American Chemical Society, vol. 137, no. 15, pp. 5100–5112, 2015

Zhao

, Chang

, Xiao

, Gu

, He

, and Zhao

L. D.

, “Investigations on distinct thermoelectric transport behaviors of Cu in n-type PbS,” Journal of Alloys and Compounds, vol. 781, pp. 820–830, 2019

Wang

, Wen

, Fan

, Bao

, Huang

, Tu

, and Wang

, “Energy-filtering-induced high power factor in PbS-nanoparticles-embedded TiS₂,” AIP Advances, vol. 5, no. 4, article 047126, 2015

Chang

, Xiao

, Zhang

, Pei

, Li

, Ma

, Yuan

, Liu

, Gong

, and Zhao

L. D.

, “High performance thermoelectrics from earth-abundant materials: enhanced figure of merit in PbS through nanostructuring grain size,” Journal of Alloys and Compounds, vol. 664, pp. 411–416, 2016

Qian

, Zheng

, Xiao

, Chang

, and Zhao

L. D.

, “Enhancing thermoelectric performance of n-type PbSe via additional meso-scale phonon scattering,” Inorganic Chemistry Frontiers, vol. 4, no. 4, pp. 719–726, 2017

Cutler

, Leavy

J. F.

, and Fitzpatrick

R. L.

, “Electronic transport in semimetallic cerium sulfide,” Physical Review, vol. 133, no. 4A, pp. A1143–A1152, 1964

Poudel

, Hao

, Ma

, Lan

, Minnich

, Yu

, Yan

, Wang

, Muto

, Vashaee

, Chen

, Liu

, Dresselhaus

M. S.

, Chen

, and Ren

, “High-thermoelectric performance of nanostructured bismuth antimony telluride bulk alloys,” Science, vol. 320, no. 5876, pp. 634–638, 2008

Zhao

L.-D.

, He

, Hao

, Wu

C. I.

, Hogan

T. P.

, Wolverton

, Dravid

V. P.

, and Kanatzidis

M. G.

, “Raising the thermoelectric performance of p-type PbS with endotaxial nanostructuring and valence-band offset engineering using CdS and ZnS,” Journal of the American Chemical Society, vol. 134, no. 39, pp. 16327–16336, 2012

Zhao

L.-D.

, He

, Wu

C. I.

, Hogan

T. P.

, Zhou

, Uher

, Dravid

V. P.

, and Kanatzidis

M. G.

, “Thermoelectrics with earth abundant elements: high performance p-type PbS nanostructured with SrS and CaS,” Journal of the American Chemical Society, vol. 134, no. 18, pp. 7902–7912, 2012

Johnsen

, He

, Androulakis

, Dravid

V. P.

, Todorov

, Chung

D. Y.

, and Kanatzidis

M. G.

, “Nanostructures boost the thermoelectric performance of PbS,” Journal of the American Chemical Society, vol. 133, no. 10, pp. 3460–3470, 2011

Qin

, Wang

, He

, Zhang

, Wu

, Pennycook

S. J.

, and Zhao

L. D.

, “Realizing high thermoelectric performance in p-type SnSe through crystal structure modification,” Journal of the American Chemical Society, vol. 141, no. 2, pp. 1141–1149, 2019

Blöchl

P. E.

, “Projector augmented-wave method,” Physical Review B, vol. 50, no. 24, pp. 17953–17979, 1994

Kresse

, and Joubert

, “From ultrasoft pseudopotentials to the projector augmented-wave method,” Physical Review B, vol. 59, no. 3, pp. 1758–1775, 1999

Kresse

, and Furthmüller

, “Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set,” Physical Review B, vol. 54, no. 16, pp. 11169–11186, 1996

Zhou

, Li

, Wu

, Zhao

L. D.

, He

, Wei

S. H.

, and Huang

, “Influence of defects on the thermoelectricity in SnSe: a comprehensive theoretical study,” Physical Review B, vol. 97, no. 24, article 245202, 2018

Wei

S. H.

, “Overcoming the doping bottleneck in semiconductors,” Computational Materials Science, vol. 30, no. 3-4, pp. 337–348, 2004

Yuan

Z.-K.

, Chen

, Xie

, Park

J. S.

, Xiang

, Gong

X. G.

, and Wei

S. H.

, “Na-diffusion enhanced p-type conductivity in Cu(In, Ga) Se₂: a new mechanism for efficient doping in semiconductors,” Advanced Energy Materials, vol. 6, no. 24, article 1601191, 2016

Xia

, Jia

, and Zhang

, “First-principles electronic structure and formation energies of group V and VII impurities in the α-Fe₂O₃ alloys,” Journal of Applied Physics, vol. 116, no. 11, article 113706, 2014

Zhang

S. B.

, Wei

S.-H.

, and Zunger

, “Overcoming doping bottlenecks in semiconductors and wide-gap materials,” Physica B-Condensed Matter, vol. 273-274, pp. 976–980, 1999

Moore

G. A.

, “Modern thermoelectrics,” Electronics and Power, vol. 30, no. 9, p. 733, 1984

Mori

, Usui

, Ochi

, and Kuroki

, “Temperature- and doping-dependent roles of valleys in the thermoelectric performance of SnSe: a first-principles study,” Physical Review B, vol. 96, no. 8, p. 10, 2017

Pei

, LaLonde

A. D.

, Heinz

N. A.

, Shi

, Iwanaga

, Wang

, Chen

, and Snyder

G. J.

, “Stabilizing the optimal carrier concentration for high thermoelectric efficiency,” Advanced Materials, vol. 23, no. 47, pp. 5674–5678, 2011

Zhao

L.-D.

, Dravid

V. P.

, and Kanatzidis

M. G.

, “The panoscopic approach to high performance thermoelectrics,” Energy & Environmental Science, vol. 7, no. 1, pp. 251–268, 2014

Ahn

, Biswas

, He

, Chung

, Dravid

, and Kanatzidis

M. G.

, “Enhanced thermoelectric properties of p-type nanostructured PbTe-MTe (M = Cd, Hg) materials,” Energy & Environmental Science, vol. 6, no. 5, pp. 1529–1537, 2013

Zhao

, Chang

, Xiao

, and Zhao

L. D.

, “High performance of n-type (PbS)_1-x-y(PbSe)_x(PbTe)_y thermoelectric materials,” Journal of Alloys and Compounds, vol. 744, pp. 769–777, 2018

Wan

C. L.

, Pan

, Xu

, Qin

Y. X.

, Wang

J. D.

, Qu

Z. X.

, and Fang

M. H.

, “Effect of point defects on the thermal transport properties of(La_xGd_1−x)₂Zr₂O₇: experiment and theoretical model,” Physical Review B, vol. 74, no. 14, article 144109, 2006

Pei

Y.-L.

, He

, Li

J. F.

, Li

, Liu

, Pan

, Barreteau

, Berardan

, Dragoe

, and Zhao

L. D.

, “High thermoelectric performance of oxyselenides: intrinsically low thermal conductivity of Ca-doped BiCuSeO,” NPG Asia Materials, vol. 5, no. 5, article e47, 2013

Tan

, Shi

, Sun

, Zhao

L. D.

, Uher

, Dravid

V. P.

, and Kanatzidis

M. G.

, “SnTe–AgBiTe₂ as an efficient thermoelectric material with low thermal conductivity,” Journal of Materials Chemistry A, vol. 2, no. 48, pp. 20849–20854, 2014

Xiao

, Wang

, Zhang

, Chen

, Zhang

, Wang

, Pennycook

S. J.

, Snyder

G. J.

, Wu

, and Zhao

L. D.

, “Band sharpening and band alignment enable high quality factor to enhance thermoelectric performance in n-type PbS,” Journal of the American Chemical Society, vol. 142, no. 8, pp. 4051–4060, 2020

Qin

, He

, and Zhao

L.-D.

, “Estimation of the potential performance in p-type SnSe crystals through evaluating weighted mobility and effective mass,” Journal of Materiomics, vol. 6, no. 4, pp. 671–676, 2020

Appendix

Less

Year 2020 volume 2020 Issue 12

PDF

Cite this Article

BibTeX

Article Info

doi: 10.34133/2020/4084532

Receive Date：2020-09-26
Online Date：2026-10-27
Published：2020-12-20

Article Data

Affiliations

History

Received：2020-09-26
Accepted：2020-10-25

Funding

China Postdoctoral Science Foundation(2019M650429)

Henan Normal University

Beihang University

National Postdoctoral Program for Innovative Talents(BX20200028)

Academic Excellence Foundation of BUAA for PhD Students

National Science Fund for Distinguished Young Scholars(51925101)

111 Project(B17002)

Beijing Municipal Natural Science Foundation(JQ18004)

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

National Natural Science Foundation of China(51671015; 51772012)

Affiliations

¹School of Materials Science and Engineering, Beihang University, Beijing 100191, China

²School of Physics, Henan Normal University, Xinxiang 453007, China

References

Share

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

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

Table 1: Calculated imperfection scaling parameters and $κ lat$ (W/mK) of Li- and Na-doped PbS based on the Callaway model.

Samples	$Γ M, Pb, Q$	$Γ S, Pb, Q$	$Γ Pb, Q$	$κ lat$
Pb_0.995Li_0.0050S	0.004693	0.000072	0.010254	2.03
Pb_0.995Li_0.0075S	0.007055	0.000107	0.015382	1.87
Pb_0.99Li_0.0100S	0.009429	0.000143	0.020510	1.75
Pb_0.98Li_0.0200S	0.019038	0.000284	0.041028	1.45
Pb_0.995Na_0.0050S	0.003967	0.000042	0.007189	2.19
Pb_0.995Na_0.0075S	0.005962	0.000062	0.010780	2.06
Pb_0.99Na_0.0100S	0.007965	0.000083	0.014370	1.96
Pb_0.98Na_0.0200S	0.016056	0.000163	0.028714	1.68

Figure 1: Powder XRD patterns of (a) Pb_1-xLi_xS and (b) Pb_1-xNa_xS.

Figure 2: Temperature dependence of thermoelectric transport properties of Pb_1-xLi_xS and Pb_1-xNa_xS: (a, b) electrical conductivity, (c, d) Seebeck coefficient, (e, f) power factor (PF), (g) comparisons of maximum PF, and (h) averaged PF at 423-730 K.

Figure 3: The calculated formation energy of potential defect in PbS under (a) Pb-rich and (b) S-rich conditions. The Fermi level with respect to the valence band maximum (VBM) of PbS. (c, d) The carrier mobility and carrier concentration of Pb_1-xLi_xS and Pb_1-xNa_xS at room temperature.

Figure 4: Temperature dependence of thermal transport properties of Pb_1-xLi_xS and Pb_1-xNa_xS: (a, b) total thermal conductivity, (c, d) Lorenz number, and (e, f) electronic thermal conductivity.

Figure 5: Temperature dependence of thermal transport properties of Pb_1-xLi_xS and Pb_1-xNa_xS: (a, b) lattice thermal conductivity and (c, d) comparisons of $κ lat$ between experimental and calculated values.

Figure 6: Temperature dependence of ZT of (a) Pb_1-xLi_xS and (b) Pb_1-xNa_xS and comparisons of (c) maximum ZT and (d) averaged ZT at 423-730 K.

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